2 // All cuts for single pions in Sigma* analysis 2010,
3 // based on track quality and particle identification
5 // Author: Massimo Venaruzzo.
12 #include "AliPIDResponse.h"
13 #include "AliRsnCutDaughterSigmaStar2010PP.h"
15 ClassImp(AliRsnCutDaughterSigmaStar2010PP)
17 //__________________________________________________________________________________________________
18 AliRsnCutDaughterSigmaStar2010PP::AliRsnCutDaughterSigmaStar2010PP(const char *name, AliPID::EParticleType pid) :
19 AliRsnCut(name, AliRsnTarget::kDaughter),
21 fCutQuality(Form("%sQuality", name)),
24 //fDCARptFormula("0.0182+0.0350/pt^1.01")
29 // Initialize track quality cuts to 2010 defaults
32 fCutQuality.SetPtRange(0.15, 1E+20);
33 fCutQuality.SetEtaRange(-0.8, 0.8);
34 fCutQuality.SetDCARmax(fDCARmax);
35 //fCutQuality.SetDCARPtFormula(fDCARptFormula);
36 fCutQuality.SetDCAZmax(2.0);
37 fCutQuality.SetSPDminNClusters(1);
38 fCutQuality.SetITSminNClusters(0);
39 fCutQuality.SetITSmaxChi2(1E+20);
40 fCutQuality.SetTPCminNClusters(fMinTPCcluster);
41 fCutQuality.SetTPCmaxChi2(4.0);
42 fCutQuality.SetRejectKinkDaughters();
43 fCutQuality.SetAODTestFilterBit(5);
46 //__________________________________________________________________________________________________
47 Bool_t AliRsnCutDaughterSigmaStar2010PP::IsSelected(TObject *obj)
54 if (!TargetOK(obj)) return kFALSE;
57 AliVTrack *track = fDaughter->Ref2Vtrack();
59 if (!fDaughter->GetRef()) AliWarning("NULL ref");
63 AliDebugClass(2, "Checking status");
67 if ((track->GetStatus() & AliESDtrack::kTPCin ) == 0) return kFALSE;
68 if ((track->GetStatus() & AliESDtrack::kTPCrefit) == 0) return kFALSE;
69 if ((track->GetStatus() & AliESDtrack::kITSrefit) == 0) return kFALSE;
71 AliDebugClass(2, "Checking quality cuts");
74 if (!fCutQuality.IsSelected(obj)) return kFALSE;
76 AliDebugClass(2, "...passed");
79 // check initialization of PID object
80 AliPIDResponse *pid = fEvent->GetPIDResponse();
82 AliFatal("NULL PID response");
86 // check if TOF is matched
87 // and computes all values used in the PID cut
88 //Bool_t isTOF = MatchTOF(track);
89 //Double_t pTPC = track->GetTPCmomentum();
90 //Double_t p = track->P();
91 Double_t nsTPC = TMath::Abs(pid->NumberOfSigmasTPC(track, fPID));
92 //Double_t nsTOF = TMath::Abs(pid->NumberOfSigmasTOF(track, fPID));
93 Double_t maxTPC = 1E20;
94 //Double_t maxTOF = 1E20;
96 // applies the cut differently depending on the PID and the momentum
98 // TPC: 5sigma cut for all
99 //if (nsTPC > 5.0) return kFALSE;
100 // TOF: 3sigma below 1.5 GeV, 2sigma above
101 //if (p < 1.5) maxTOF = 3.0; else maxTOF = 2.0;
102 //return (nsTOF <= maxTOF);
105 //For the moment TOF is not used - PID ONLY WITH TPC - 3 sigmas in the whole range
110 return (nsTPC <= maxTPC);
112 AliDebugClass(2, "Good Pion (hallelujah)");