5 /* History of cvs commits:
8 * Revision 1.12 2007/02/20 20:17:43 hristov
9 * Corrected array size, removed warnings (icc)
11 * Revision 1.11 2006/12/19 08:49:35 gustavo
12 * New PID class for EMCAL, bayesian analysis done with ESD data, PID information filled when calling AliEMCALPID in AliEMCALReconstructor::FillESD()
17 ///////////////////////////////////////////////////////////////////////////////
19 ///////////////////////////////////////////////////////////////////////////////
23 #include "AliESDEvent.h"
26 class AliEMCALPID : public TTask {
31 virtual ~AliEMCALPID() { }
33 void RunPID(AliESDEvent *esd);
34 void ComputePID(Double_t energy, Double_t lambda0); // give the PID of a cluster
35 TArrayD DistLambda0(Double_t energy, Int_t nature); // compute lambda0 distributions
37 Double_t GetPID(Int_t idx) const {if (idx>=0&&idx<3) return fPID[idx]; else return 0.;}
38 Double_t GetPIDFinal(Int_t idx) const {if (idx>=0&&idx<AliPID::kSPECIESN) return fPIDFinal[idx]; else return 0.;}
39 Double_t GetPIDWeight(Int_t idx) const {if (idx>=0&&idx<3) return fPIDWeight[idx]; else return 0.;}
41 void SetPID(Double_t val, Int_t idx) {if (idx>=0&&idx<3) fPID[idx] = val;}
42 void SetPIDFinal(Double_t val, Int_t idx) {if (idx>=0&&idx<AliPID::kSPECIESN) fPIDFinal[idx] = val;}
43 void SetPIDWeight(Double_t val, Int_t idx) {if (idx>=0&&idx<3) fPIDWeight[idx] = val;}
44 void SetPrintInfo(Bool_t yesno) {fPrintInfo = yesno;}
45 void SetReconstructor(Bool_t yesno) {fReconstructor = yesno;}
48 Double_t Polynomial(Double_t x, Double_t *params);
50 Bool_t fPrintInfo; // flag to decide if details about PID must be printed
52 Double_t fGamma[6][6]; // Parameter to Compute PID
53 Double_t fHadron[6][6]; // Parameter to Compute PID
54 Double_t fPiZero5to10[6][6]; // Parameter to Compute PID
55 Double_t fPiZero10to60[6][6]; // Parameter to Compute PID
59 Float_t fPIDFinal[AliPID::kSPECIESN+1]; // final PID format
60 Float_t fPIDWeight[3]; // order: gamma, pi0, hadrons,
61 Double_t fProbGamma; // probility to be a Gamma
62 Double_t fProbPiZero; // probility to be a PiO
63 Double_t fProbHadron; // probility to be a Hadron
64 Bool_t fReconstructor; //Fill esdcalocluster when called from EMCALReconstructor
66 ClassDef(AliEMCALPID, 0)
69 #endif // ALIEMCALPID_H