flyopunta/src/src_ana/anakpip0.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 "anakpip0.h"
using namespace std;

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

 AnaKpip0::AnaKpip0()
   {tipo=29;nfit=2;titol="K--> pi+ pi0";nome="Kpip0";
    ntrk=good=0;
     mas_pip=ptr_Pip->Get_Massa();
     mas_Kp=ptr_Kapp->Get_Massa();
     mas_pi0=ptr_Pi0->Get_Massa();
    Gout<<"\n Attivato il fit di "<<titol;
   }                  
                                                
  AnaKpip0::~AnaKpip0(){};

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

int AnaKpip0::Fisica()
  {

  count_call++;
//  
  
     if(Trkstraw!=0)  ntrk=Trkstraw->ntrk_done;
  if(ntrk<1) {good=0;return(-1);}
 
   Pkin.setvn(-0.09,0.,75.,mas_Kp);
 Pkin=SelRea->avo->P;
     
  Track  Trk=Trkstraw->Traccia[0];
  charge=Trk.charge;
  chisq=Trk.chiq_fit; 
  Ppi.setvn(Trkstraw->Traccia[0].P_in,mas_pip);
  
 //float charge,chisq,pmis,masmis, ppip; 
  ppip=Ppi.mom;
 
  Ppi0=Pkin-Ppi;
  Ppi0.Invar();
 // Ppi0.print("P_mis");
  
  masmis=Ppi0.m;
  pmis=Ppi0.mom;
  
  if(ppip<15.||ppip>35.)return -1;
  
    count_wnt++;
    good=1;
   return 1;

}

//======================================================================================
 int AnaKpip0::Vertice()
   {

     return 1;
   }

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

     Gout<<"\n -----> "<<titol<<"   Analysis Summary. <-------\n";
     Gout<<"\n Fit entries          "<<count_call;
     Gout<<"\n Track energ to low   "<<getta_ch;
     Gout<<"\n Killed by energy     "<<getta_en;
     Gout<<"\n Tau cut              "<<getta_tau;
     Gout<<"\n No good vertex       "<<getta_fit;
     Gout<<"\n Good events          "<<count_wnt;
     Gout<<"\n\n ======>           d o n e     <============ ";
    }
//-------------------------------------------