PWG2/UNICOR/AliDEventAliceESD.cxx

   1 // Author: Dariusz Miskowiec <mailto:d.miskowiec@gsi.de> 2008
   2
   3 //=============================================================================
   4 // AliDEvent-AliESD interface                                                   //
   5 //=============================================================================
   6
   7 #include <cmath>
   8 #include "AliESDEvent.h"
   9 #include "AliDEventAliceESD.h"
  10
  11 ClassImp(AliDEventAliceESD)
  12
  13 //=============================================================================
  14 AliDEventAliceESD::AliDEventAliceESD() : AliDEvent(), AliESDEvent()
  15 {
  16   // constructor
  17
  18   printf("%s object created\n",AliDEvent::ClassName());
  19 }
  20 //=============================================================================
  21 AliDEventAliceESD::~AliDEventAliceESD()
  22 {
  23   // destructor
  24
  25 }
  26 //=============================================================================
  27 Bool_t AliDEventAliceESD::Good() const
  28 {
  29   // event cuts
  30
  31   if (fabs(Zver())>1) return kFALSE;
  32   if (AliESDEvent::GetPrimaryVertex()->GetZRes()>0.1) return kFALSE;
  33   return kTRUE;
  34 }
  35 //=============================================================================
  36 Bool_t AliDEventAliceESD::ParticleGood(Int_t i, Int_t pidi) const
  37 {
  38   // track cuts and particle id cuts; pidi=0 means take all tracks
  39   // consider using the standard ESDcut
  40
  41   // track quality cuts
  42
  43   AliESDtrack *track = AliESDEvent::GetTrack(i);
  44   if (!track->IsOn(AliESDtrack::kTPCrefit)) return 0;        // TPC refit
  45   if (track->GetTPCNcls() < 120) return 0;                   // number of TPC clusters
  46   const AliExternalTrackParam *tp = track->GetTPCInnerParam();
  47   if (!tp) return 0;
  48
  49   Float_t r,z;
  50   track->GetImpactParameters(r,z);
  51   //  if (fabs(z)>0.2) return 0;                          // impact parameter in z
  52   //  if (fabs(r)>0.1) return 0;                          // impact parameter in xy
  53
  54   // pid
  55
  56   if (pidi==0) return 1;
  57   if (!track->IsOn(AliESDtrack::kTPCpid)) return 0;
  58   Double_t p[AliPID::kSPECIES];
  59   track->GetESDpid(p);
  60   Int_t q = tp->Charge();
  61   if (pidi == -211) return p[AliPID::kPion]>0.5 && q==-1;
  62   else if (pidi == 211) return p[AliPID::kPion]>0.5 && q==1;
  63   else return 0;
  64 }
  65 //=============================================================================
  66 Bool_t AliDEventAliceESD::PairGood(Double_t /*p0*/, Double_t the0, Double_t phi0,
  67                                    Double_t /*p1*/, Double_t the1, Double_t phi1) const {
  68
  69   // two-track separation cut
  70
  71   double r = 85; // TPC entrance radius in cm
  72   double x0 = r*sin(the0)*cos(phi0);
  73   double x1 = r*sin(the1)*cos(phi0);
  74   double y0 = r*sin(the0)*sin(phi0);
  75   double y1 = r*sin(the1)*sin(phi1);
  76   double z0 = r*cos(the0);
  77   double z1 = r*cos(the1);
  78   double dx = x1-x0;
  79   double dy = y1-y0;
  80   double dz = z1-z0;
  81   double dist2 = dx*dx+dy*dy+dz*dz;
  82   return dist2>2*2;
  83 }
  84 //=============================================================================