PWG2/FEMTOSCOPY/AliFemto/AliFemtoBasicTrackCut.cxx

   1 ////////////////////////////////////////////////////////////////////////////////
   2 //                                                                            //
   3 // AliFemtoBasicTrackCut - the basic cut for tracks.                          //
   4 // Cuts on particle identification, transverse momentum, rapidity, distance   //
   5 // of closest approach to primary vertex and charge                           //
   6 //                                                                            //
   7 ////////////////////////////////////////////////////////////////////////////////
   8
   9 #include "AliFemtoBasicTrackCut.h"
  10 #include <cstdio>
  11
  12 #ifdef __ROOT__
  13 ClassImp(AliFemtoBasicTrackCut)
  14 #endif
  15
  16 AliFemtoBasicTrackCut::AliFemtoBasicTrackCut():
  17   fCharge(0),
  18   fNTracksPassed(0),
  19   fNTracksFailed(0)
  20 {
  21   // Default constructor
  22   fNTracksPassed = fNTracksFailed = 0;
  23   fCharge = 1;  // takes both charges 0
  24   fNSigmaPion[0] = -100.0;   fNSigmaPion[1] = 100.0;
  25   fNSigmaKaon[0] = -100.0;   fNSigmaKaon[1] = 100.0;
  26   fNSigmaProton[0] = -100.0; fNSigmaProton[1] = 100.0;
  27   fNHits[0] = 10;          fNHits[1] = 180;
  28   fPt[0]=0.0;              fPt[1] = 100.0;//100
  29   fRapidity[0]=-2;       fRapidity[1]=2;//-2 2
  30   fDCA[0] = -1.0;           fDCA[1] = 20.0;
  31
  32 }
  33 //------------------------------
  34 //AliFemtoBasicTrackCut::~AliFemtoBasicTrackCut(){
  35 //  /* noop */
  36 //}
  37 //------------------------------
  38 bool AliFemtoBasicTrackCut::Pass(const AliFemtoTrack* track){
  39   // test the particle and return
  40   // true if it meets all the criteria
  41   // false if it doesn't meet at least one of the criteria
  42
  43   //  return true ;  // THIS CUT IS A STHBTDUMMY!!
  44
  45   /*
  46     cout << endl;
  47     cout << "#track " << trackCount++;
  48     cout << " * pion " << (track->NSigmaPion() > fNSigmaPion[0]) && (track->NSigmaPion() < fNSigmaPion[1]);
  49     cout << " * kaon " << (track->NSigmaKaon() > fNSigmaKaon[0]) && (track->NSigmaKaon() < fNSigmaKaon[1]);
  50     cout << " * proton " << (track->NSigmaProton() > fNSigmaProton[0]) && (track->NSigmaProton() < fNSigmaProton[1]);
  51     cout << " * charge " << (track->Charge() == fCharge);
  52   */
  53   bool goodPID = 1;
  54   /* ----- NOT DOING PID CUTS !!!! ------
  55   bool goodPID = ((track->NSigmaPion()   > fNSigmaPion[0]) &&
  56                   (track->NSigmaPion()   < fNSigmaPion[1]) &&
  57                   (track->NSigmaKaon()   > fNSigmaKaon[0]) &&
  58                   (track->NSigmaKaon()   < fNSigmaKaon[1]) &&
  59                   (track->NSigmaProton() > fNSigmaProton[0]) &&
  60                   (track->NSigmaProton() < fNSigmaProton[1]));
  61   ----- NOT DOING PID CUTS !!!! ------ */
  62   if (fCharge !=0){               // if user requests "charge=0" then that means ignore charge
  63     goodPID = (goodPID&&(track->Charge() == fCharge));
  64   }
  65   if (goodPID){
  66     float tEnergy = ::sqrt(track->P().mag2()+fMass*fMass);
  67     float tRapidity = 0.5*::log((tEnergy+track->P().z())/
  68                             (tEnergy-track->P().z()));
  69
  70     float tPt = ::sqrt((track->P().x())*(track->P().x())+
  71                     (track->P().y())*(track->P().y()));
  72
  73
  74     /*
  75       cout << " * DCAxy " << (track->DCAxy()  > fDCA[0]) && (track->DCAxy()  < fDCA[1]);
  76       cout << " * fDCA[0] " << fDCA[0];
  77       cout << " * fDCA[1] " << fDCA[1];
  78       cout << " * track->DCAxy " << track->DCAxy();
  79       cout << " * NHits " <<  (track->NHits() > fNHits[0]) && (track->NHits() < fNHits[1]);
  80       cout << " * tPt " << (tPt > fPt[0]) && (tPt < fPt[1]);
  81       cout << " * y " << (tRapidity > fRapidity[0]) && (tRapidity < fRapidity[1]);
  82       cout << endl;
  83     */
  84
  85     bool goodTrack=
  86       (//(track->DCAxy()  > fDCA[0]) &&
  87      //  (track->DCAxy()  < fDCA[1]) &&
  88   //     (track->NHits() > fNHits[0]) &&
  89     //   (track->NHits() < fNHits[1]) &&
  90        (tPt             > fPt[0]) &&
  91        (tPt             < fPt[1]) &&
  92        (tRapidity      > fRapidity[0]) &&
  93        (tRapidity      < fRapidity[1]))&&
  94        (track->PidProbPion()>0.5)&&//moje
  95        (track->PidProbMuon()<0.47)&&//moje
  96        (track->Label()>0);//moje
  97
  98     //    cout << track->DCAxy() << " " << track->NHits() << " " << Pt << " " << tRapidity << " " << tEnergy << endl;
  99
 100     goodTrack ? fNTracksPassed++ : fNTracksFailed++;
 101     return (goodTrack);
 102   }
 103   else{
 104     fNTracksFailed++;
 105     return (goodPID);
 106   }
 107 }
 108 //------------------------------
 109 AliFemtoString AliFemtoBasicTrackCut::Report(){
 110   // construct report
 111   string tStemp;
 112   char tCtemp[100];
 113   sprintf(tCtemp,"Particle mass:\t%E\n",this->Mass());
 114   tStemp=tCtemp;
 115   sprintf(tCtemp,"Particle charge:\t%d\n",fCharge);
 116   tStemp=tCtemp;
 117   sprintf(tCtemp,"Particle Nsigma from pion:\t%E - %E\n",fNSigmaPion[0],fNSigmaPion[1]);
 118   tStemp+=tCtemp;
 119   sprintf(tCtemp,"Particle Nsigma from kaon:\t%E - %E\n",fNSigmaKaon[0],fNSigmaKaon[1]);
 120   tStemp+=tCtemp;
 121   sprintf(tCtemp,"Particle Nsigma from proton:\t%E - %E\n",fNSigmaProton[0],fNSigmaProton[1]);
 122   tStemp+=tCtemp;
 123   sprintf(tCtemp,"Particle #hits:\t%d - %d\n",fNHits[0],fNHits[1]);
 124   tStemp+=tCtemp;
 125   sprintf(tCtemp,"Particle pT:\t%E - %E\n",fPt[0],fPt[1]);
 126   tStemp+=tCtemp;
 127   sprintf(tCtemp,"Particle rapidity:\t%E - %E\n",fRapidity[0],fRapidity[1]);
 128   tStemp+=tCtemp;
 129   sprintf(tCtemp,"Particle DCA:\t%E - %E\n",fDCA[0],fDCA[1]);
 130   tStemp+=tCtemp;
 131   sprintf(tCtemp,"Number of tracks which passed:\t%ld  Number which failed:\t%ld\n",fNTracksPassed,fNTracksFailed);
 132   tStemp += tCtemp;
 133   AliFemtoString returnThis = tStemp;
 134   return returnThis;
 135 }