flyopunta/src/src_ana/anakpip0ee.cpp

/***************************************************************************
                          anak2pi.cpp  -  description
   GMP 27-11-2010
  
  ***************************************************************************/
  // #define MIN(x,y,z) ((x)<(y)?(x)<(z)?(x):(z):(y)<(z)?(y):(z))
  
#include "flyoh.h"
#include "anakpip0ee.h"
using namespace std;

//-------------------------------------------------------------
//-------------------------------------------------------------

 AnaKpip0ee::AnaKpip0ee()
   {tipo=31;nfit=2;titol="K--> pipi0ee";nome="Kpip0ee";
    ntrk=good=0;
     mas_pip=ptr_Pip->Get_Massa();
     mas_Kp=ptr_Kapp->Get_Massa();
     mas_pi0=ptr_Pi0->Get_Massa();
     mas_ele=ptr_Elec->Get_Massa();
    Gout<<"\n Attivato il fit di "<<titol;
    Pk_ric=new GigaBeam();
   }                  
                                                
  AnaKpip0ee::~AnaKpip0ee(){};

 //=============== C h a r g e  2 pions    F i t ========

int AnaKpip0ee::Fisica()
  {
//     int good;
//     float chiq_vert,xvtx,yvtx,zvtx, mqinv,mqgg,mqee;
//     float Pk_mom,pxkp,pykp,pzkp,pxpip,pypip,pzpip,pxelec,pyelec,pzelec,pxelep,pyelep,pzelep;
  count_call++;
//  
  good=0;;
  chiq_vert=xvtx=yvtx=zvtx=mqinv=mqgg=mqee=0.0;
   momk=pxkp=pykp=pzkp=pxpip=pypip=pzpip=pxelec=pyelec=pzelec=pxelep=pyelep=pzelep=0.0;
 
    //taglio sul vertice generato (taglia solo gli eventi che finiscono sul collimatore)da togliere!!!
    if (evento_.Gen.zv < 10000.) {
        getta_cog++;
        return -1;
    }

    count_decay++;
 
    
    if (Trkstraw==0) { cout<<"\n Tracce non attivato!!! exit "<<endl;exit(0);}
    if( Trkstraw->ntrk_done<3)
    {
        getta_trk++;
        return -1;
    }

//*****************************************
    
    Track  Trk0=Trkstraw->Traccia[0];
    Track  Trk1=Trkstraw->Traccia[1];    
    Track  Trk2=Trkstraw->Traccia[2];     

    if((Trk0.charge+Trk1.charge+Trk2.charge) != 1)
    {// wrong charge!!!!
         getta_char++;
        return -1;     
    }


    // ricavo l'impulso del K dal gigatrak
    Pk_mom=Pk_ric->Beam_Rec();   // impulso del beam dal gigatrack
    Pkin.setvn ( Pk_mom,mas_Kp ); // definizione del quadrimomento
    Pkin.putvn( pxkp,pykp,pzkp,momk );  // salvo per root

//    double mas_Kp, mas_pip, mas_pi0, mas_ele;
//    qvet  Pkin, Ppi, Ppi0, Pelec, Pelep,Pgam1,Pgam2;

// calcolo del vertice  
    Vtxstruct Xdecay(4);   
    Xdecay.P[3].setvn(0.,0.,10240.); // beam at G3 (0.,0.,10240.); Fare meglio
    Xdecay.V[3]=Pkin.NVerso();       //  versore del beam...    
    Xdecay.P[0]=Trkstraw->Traccia[0].Xmg; // coordinate della traccia 0 al centro del magnete 
    Xdecay.V[0]=Trkstraw->Traccia[0].P_in.NVerso(); // versore della traccia 0     
    Xdecay.P[1]=Trkstraw->Traccia[1].Xmg; // coordinate della traccia 1 al centro del magnete 
    Xdecay.V[1]=Trkstraw->Traccia[1].P_in.NVerso(); // versore della traccia 1     
    Xdecay.P[2]=Trkstraw->Traccia[2].Xmg; // coordinate della traccia 2 al centro del magnete 
    Xdecay.V[2]=Trkstraw->Traccia[2].P_in.NVerso(); // versore della traccia 2
    Xdecay.Vertice();
    chiq_vert=Xdecay.chiq;
    Xdecay.Vtx.putv( xvtx,yvtx,zvtx );   
    
    if(chiq_vert>100.)
    {// wrong fit!!!!
        getta_fit++;
        return -1;     
    }
    
// cinematica     
// seleziono l'indice dell'elettrone
    if(lkry->mhit<5)
      {// not enough hits in lkr
        getta_ch++;
        return -1;     
    }
    
// matching  lkr-tracks
int Copp[3]={-1,-1,-1};
int Copg[2]={-1,-1};
int ig=0,ipip=-1;

for(int i=0;i<5;i++)
{
  gvet Xl=lkry->M_Hits[i].Xlab;  
  for (int j=0;j<3;j++)
  {
    gvet Xdl= Xl-Trkstraw->Traccia[j].Xlkr;
    if( Xdl.XYNorma()<2.5)
    {
      Copp[j]=i;   
      if(lkry->me[i]/Trkstraw->Traccia[j].P_in.mom<.85) ipip=i; // è il pip
      break;
    }   
  }
  if(ig<2)Copg[ig]=i;
  ig++;
 } 
  
  if(Copp[0]==-1||Copp[1]==-1||Copp[2]==-1||
     Copg[0]==-1||Copg[1]==-1||ipip==-1)
  {
    getta_gg++;
    return -1;
  }





    int elect=0;
    if(Trk1.charge<0)elect=1;
    if(Trk2.charge<0)elect=2;
    
    
    
     
//     
    
    
    count_wnt++;
    good=1;
   return 1;

}

//======================================================================================
 

 //----------------------------------------------------
   void AnaKpip0ee::print_scale()
    {

     Gout<<"\n -----> "<<titol<<"   Analysis Summary. <-------\n";
     Gout<<"\n Fit entries          "<<count_call;
     Gout<<"\n Z cut                "<<getta_cog;
     Gout<<"\n no tracks            "<<getta_trk;
     Gout<<"\n Wrong charge         "<<getta_char;
     Gout<<"\n No good vertex       "<<getta_fit;
     Gout<<"\n Not enough lkr hits  "<<getta_ch;
     Gout<<"\n Tracks-Lkr no match  "<<getta_gg;
     Gout<<"\n Good events          "<<count_wnt;
     Gout<<"\n\n ======>           d o n e     <============ ";
    }
//-------------------------------------------