2 -----------------------------------------------------------------------------------
3 AliFlowBayesianPID class implemented by F. Noferini (noferini@bo.infn.it)
4 needed to use Bayesian probability in the flow package (used in AliFlowTrackCuts)
5 -----------------------------------------------------------------------------------
7 #ifndef ALIFLOWBAYESIANPID_H
8 #define ALIFLOWBAYESIANPID_H
10 #include "AliESDpid.h"
11 #include "AliPIDResponse.h"
27 1) in your main program:
29 AliFlowBayesianPID *mypid = new AliFlowBayesianPID();
31 or (if you want a global AliESDpid object)
33 AliESDpid *PIDesd = new AliESDpid();
35 AliFlowBayesianPID *mypid = new AliFlowBayesianPID(PIDesd);
37 for data starting from the PbPb pass2 reconstruction
39 mypid->SetNewTrackParam();
42 before to loop on the on ESD tracks
43 mypid->SetDetResponse(esdEvent, centrality,AliESDpid::kTOF_T0,kFALSE); // centrality > 0 = PbPb or -1 for pp collisions
45 before to loop on the on AOD tracks
46 mypid->SetDetResponse(aodEvent, centrality); // centrality > 0 = PbPb or -1 for pp collisions
48 if you want to use a dE/dx depdence on DeltaPhi set the EP with its resolution otherwise it will be skipped
49 mypid->SetPsiCorrectionDeDx(psiEP,resEP);
52 for(...){ // track loop
53 mypid->ComputeProb(track,centrality); // both for ESD and AOD tracks
54 Float_t *prob = mypid->GetProb(); // Bayesian Probability (from 0 to 4) (Combined TPC || TOF) including a tuning of priors and TOF mismatch parameterization
59 // for the single detector weights (no priors)
60 Float_t *tpcWeight = mypid->GetWeights(0); // TPC weights (equal weights in case of no kTPCpid)
61 Float_t *tofWeight = mypid->GetWeights(1); // TOF weights (equal weights in case of no kTOFpid)
64 TSpline3 *mismatchShape = mypid->GetMismatch(); // time distribution of mismatched tracks - L/c (L is the distance of the pad from the PV)
65 TF1 *tpcprob = mypid->GetTPCprob(); // normalized TPC response function (in Nsigmas)
66 TF1 *tofprob = mypid->GetTOFprob(); // normalized TOF response function (in Nsigmas)
68 TH2D *hPr = mypid->GetHistoPriors(isp); // 2D (centrality - pT) histo for the priors of specie-isp (centrality < 0 means pp collisions)
69 // all the priors are normalized to the pion ones
73 class AliFlowBayesianPID : public AliPIDResponse{
75 AliFlowBayesianPID(AliESDpid *esdpid=NULL);
76 virtual ~AliFlowBayesianPID();
78 // virtual method of AliPIDResponse
82 void SetDetResponse(AliESDEvent *esd,Float_t centrality=-1.0,EStartTimeType_t flagStart=AliESDpid::kTOF_T0,Bool_t /*recomputeT0TOF*/=kFALSE);
83 void SetDetResponse(AliAODEvent *aod,Float_t centrality=-1.0,EStartTimeType_t flagStart=AliESDpid::kTOF_T0);
84 void SetNewTrackParam(Bool_t flag=kTRUE){fNewTrackParam=flag;};
85 void SetDetAND(Int_t idet){if(idet < fgkNdetectors && idet >= 0) fMaskAND[idet] = kTRUE;};
86 void SetDetOR(Int_t idet){if(idet < fgkNdetectors && idet >= 0) fMaskOR[idet] = kTRUE;};
87 void ResetDetAND(Int_t idet){if(idet < fgkNdetectors && idet >= 0) fMaskAND[idet] = kFALSE;};
88 void ResetDetOR(Int_t idet){if(idet < fgkNdetectors && idet >= 0) fMaskOR[idet] = kFALSE;};
89 void SetPsiCorrectionDeDx(Float_t psi,Float_t res);
90 void SetMC(Bool_t flag){fIsMC=flag;};
93 AliESDpid* GetESDpid(){return fPIDesd;};
94 const TH2D *GetHistoPriors(Int_t specie) const {if(specie >=0 && specie < fgkNspecies) return fghPriors[specie]; else return NULL;};
95 TSpline3 *GetMismatch();
96 const TF1 *GetTOFprob() const {return fTOFResponseF;};
97 const TF1 *GetTPCprob() const {return fTPCResponseF;};
98 const Float_t *GetWeights(Int_t det) const {if(det < fgkNdetectors && det >= 0){return fWeights[det];} else{return NULL;}};
99 Float_t *GetProb() {return fProb;};
100 Float_t GetTOFMismWeight() const {return fWTofMism;};
101 Float_t GetTOFMismProb() const {return fProbTofMism;};
102 Float_t GetMassOverZ() const {return fMassTOF;};
103 Float_t GetZ() const {return fZ;};
104 Bool_t GetDetANDstatus(Int_t idet) const {if(idet < fgkNdetectors && idet >= 0){return fMaskAND[idet];} else{return kFALSE;} };
105 Bool_t GetDetORstatus(Int_t idet) const {if(idet < fgkNdetectors && idet >= 0){return fMaskOR[idet];} else{return kFALSE;} };
106 Bool_t GetCurrentMask(Int_t idet) const {if(idet < fgkNdetectors && idet >= 0){return fMaskCurrent[idet];} else{return kFALSE;} };
108 Float_t GetExpDeDx(const AliVTrack *t,Int_t iS) const;
109 Float_t GetExpDeDx(const AliVTrack *t,Float_t m) const;
111 // methods for Bayesina Combined PID
112 void ComputeWeights(const AliESDtrack *t);
113 void ComputeProb(const AliESDtrack *t,Float_t); // obsolete method
114 void ComputeProb(const AliESDtrack *t){ComputeProb(t,0.0);};
115 void ComputeWeights(const AliAODTrack *t,const AliAODEvent *aod=NULL);
116 void ComputeProb(const AliAODTrack *t,const AliAODEvent *aod=NULL); // obsolete method
118 void SetTOFres(Float_t res){fTOFresolution=res;};
120 Float_t GetDeDx() const {return fDedx;};
122 void ForceOldDedx(Bool_t status=kTRUE) {fForceOldDedx=status;};
127 static const Int_t fgkNdetectors = 2; // Number of detector used for PID
128 static const Int_t fgkNspecies = 9;// 0=el, 1=mu, 2=pi, 3=ka, 4=pr, 5=deuteron, 6=triton, 7=He3
129 static TH2D* fghPriors[fgkNspecies]; // histo with priors (hardcoded)
130 static TSpline3 *fgMism; // function for mismatch
132 AliESDpid *fPIDesd;//ESDpid object
133 TDatabasePDG *fDB; // Database pdg
134 Double_t fMass[fgkNspecies]; // mass for el(0),mu(1),pi(2),K(3),p(4)
136 Bool_t fNewTrackParam; // switch for new tracking resolution TOF parameterization
137 Float_t fTOFresolution; // TOF res needed only if T0-TOF should be recomputed
139 AliFlowBayesianPID(const AliFlowBayesianPID&); // not implemented
140 AliFlowBayesianPID& operator=(const AliFlowBayesianPID&); // not implemented
142 TF1 *fTOFResponseF; // TOF Gaussian+tail response function (tail at 1.1 sigma)
143 TF1 *fTPCResponseF; // TPC Gaussian+tail response function (tail at 1.8 sigma)
145 Float_t fWeights[fgkNdetectors][fgkNspecies]; // weights: 0=tpc,1=tof
146 Float_t fProb[fgkNspecies],fWTofMism,fProbTofMism; // Bayesian Combined PID + mismatch weights and probability
148 Float_t fZ,fMassTOF; //measured charge(Z) and mass/Z
149 TF1 *fBBdata; // Bethe Bloch function (needed to compute the charge of the particle)
151 Bool_t fMaskAND[fgkNdetectors],fMaskOR[fgkNdetectors],fMaskCurrent[fgkNdetectors]; // mask detector should be used
153 Float_t fCurrCentrality; // Centrality in current event
155 Float_t fPsi,fPsiRes; // RP and its resolution for the event (999 if not available) to correct dEdx for p > 1
157 Bool_t fIsMC; // switch for MC analysis
159 Bool_t fForceOldDedx; // switch to force to use old 2010 dEdx paramterization (even if PIDResponse is available)
161 Float_t fDedx; // dE/dx tuned for MC
163 Bool_t fIsTOFheaderAOD; // check the TOF header in AOD
165 static TH1D *fgHtofChannelDist; // channel distance from IP
167 ClassDef(AliFlowBayesianPID, 10); // example of analysis