]>
Commit | Line | Data |
---|---|---|
1c5acb87 | 1 | #ifndef ALICALOPID_H |
2 | #define ALICALOPID_H | |
3 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
4 | * See cxx source for full Copyright notice */ | |
5 | /* $Id: $ */ | |
6 | ||
7 | //_________________________________________________________________________ | |
bdd2a262 | 8 | // Class for PID selection with calorimeters |
9 | // The Output of the 2 main methods GetPdg is a PDG number identifying the cluster, | |
10 | // being kPhoton, kElectron, kPi0 ... as defined in the header file | |
11 | // - GetPdg(const TString calo, const Double_t * pid, const Float_t energy) | |
12 | // Reads the PID weights array of the ESDs and depending on its magnitude identifies the particle | |
0ae57829 | 13 | // - GetPdg(const TString calo,const TLorentzVector mom, const AliVCluster * cluster) |
bdd2a262 | 14 | // Recalcultes PID, the bayesian or any new one to be implemented in the future |
15 | // Right now only the possibility to recalculate EMCAL with bayesian and simple PID. | |
16 | // In order to recalculate Bayesian, it is necessary to load the EMCALUtils library | |
17 | // and do SwitchOnBayesianRecalculation(). | |
18 | // To change the PID parameters from Low to High like the ones by default, use the constructor | |
19 | // AliCaloPID(flux) | |
20 | // where flux is AliCaloPID::kLow or AliCaloPID::kHigh | |
21 | // If it is necessary to change the parameters use the constructor | |
22 | // AliCaloPID(AliEMCALPIDUtils *utils) and set the parameters before. | |
23 | // - SetPIDBits: Simple PID, depending on the thresholds fDispCut fTOFCut and even the | |
24 | // result of the PID bayesian a different PID bit is set. | |
25 | // | |
26 | // All these methods can be called in the analysis you are interested. | |
1c5acb87 | 27 | // |
28 | //*-- Author: Gustavo Conesa (INFN-LNF) | |
29 | ||
30 | // --- ROOT system --- | |
31 | #include <TObject.h> | |
32 | class TString ; | |
33 | class TLorentzVector ; | |
a33b76ed | 34 | //class TFormula ; |
bdd2a262 | 35 | class TTask; |
1c5acb87 | 36 | |
37 | //--- AliRoot system --- | |
0ae57829 | 38 | class AliVCluster; |
1c5acb87 | 39 | class AliAODPWG4Particle; |
bdd2a262 | 40 | class AliEMCALPIDUtils; |
f2ccb5b8 | 41 | class AliCalorimeterUtils; |
1c5acb87 | 42 | |
43 | class AliCaloPID : public TObject { | |
44 | ||
477d6cee | 45 | public: |
46 | ||
47 | AliCaloPID() ; // ctor | |
bdd2a262 | 48 | AliCaloPID(const Int_t particleFlux) ; // ctor, to be used when recalculating bayesian PID |
49 | AliCaloPID(const TTask * emcalpid) ; // ctor, to be used when recalculating bayesian PID and need different parameters | |
477d6cee | 50 | AliCaloPID(const AliCaloPID & g) ; // cpy ctor |
477d6cee | 51 | virtual ~AliCaloPID() ;//virtual dtor |
614701c6 | 52 | |
53 | private: | |
54 | AliCaloPID & operator = (const AliCaloPID & g) ;//cpy assignment | |
477d6cee | 55 | |
614701c6 | 56 | public: |
57 | ||
477d6cee | 58 | enum PidType { |
59 | kPhoton = 22, | |
60 | kPi0 = 111, | |
61 | kEta = 221, | |
62 | kElectron = 11, | |
63 | kEleCon = -11, | |
64 | kNeutralHadron = 2112, | |
65 | kChargedHadron = 211, | |
66 | kNeutralUnknown = 130, | |
67 | kChargedUnknown=321 | |
68 | }; | |
69 | ||
70 | enum TagType {kPi0Decay, kEtaDecay, kOtherDecay, kConversion, kNoTag = -1}; | |
71 | ||
72 | void InitParameters(); | |
73 | ||
74 | Int_t GetPdg(const TString calo, const Double_t * pid, const Float_t energy) const ; | |
75 | ||
0ae57829 | 76 | Int_t GetPdg(const TString calo,const TLorentzVector mom, const AliVCluster * cluster) const ; |
477d6cee | 77 | |
78 | TString GetPIDParametersList(); | |
79 | ||
f2ccb5b8 | 80 | void SetPIDBits(const TString calo, const AliVCluster * cluster, AliAODPWG4Particle *aodph, const AliCalorimeterUtils* cu); |
477d6cee | 81 | |
82 | void Print(const Option_t * opt)const; | |
83 | ||
84 | //Weight getters | |
85 | Float_t GetEMCALPhotonWeight() const { return fEMCALPhotonWeight ; } | |
86 | Float_t GetEMCALPi0Weight() const { return fEMCALPi0Weight ; } | |
87 | Float_t GetEMCALElectronWeight() const { return fEMCALElectronWeight ; } | |
88 | Float_t GetEMCALChargeWeight() const { return fEMCALChargeWeight ; } | |
89 | Float_t GetEMCALNeutralWeight() const { return fEMCALNeutralWeight ; } | |
90 | Float_t GetPHOSPhotonWeight() const { return fPHOSPhotonWeight ; } | |
91 | Float_t GetPHOSPi0Weight() const { return fPHOSPi0Weight ; } | |
92 | Float_t GetPHOSElectronWeight() const { return fPHOSElectronWeight ; } | |
93 | Float_t GetPHOSChargeWeight() const { return fPHOSChargeWeight ; } | |
94 | Float_t GetPHOSNeutralWeight() const { return fPHOSNeutralWeight ; } | |
95 | ||
a33b76ed | 96 | // Bool_t IsPHOSPIDWeightFormulaOn() const { return fPHOSWeightFormula ; } |
97 | // TFormula * GetPHOSPhotonWeightFormula() const { return fPHOSPhotonWeightFormula ; } | |
98 | // TFormula * GetPHOSPi0WeightFormula() const { return fPHOSPi0WeightFormula ; } | |
99 | // | |
477d6cee | 100 | //Weight setters |
101 | void SetEMCALPhotonWeight(Float_t w){ fEMCALPhotonWeight = w ; } | |
102 | void SetEMCALPi0Weight(Float_t w){ fEMCALPi0Weight = w ; } | |
103 | void SetEMCALElectronWeight(Float_t w){ fEMCALElectronWeight = w ; } | |
104 | void SetEMCALChargeWeight(Float_t w){ fEMCALChargeWeight = w ; } | |
105 | void SetEMCALNeutralWeight(Float_t w){ fEMCALNeutralWeight = w ; } | |
106 | void SetPHOSPhotonWeight(Float_t w){ fPHOSPhotonWeight = w ; } | |
107 | void SetPHOSPi0Weight(Float_t w){ fPHOSPi0Weight = w ; } | |
108 | void SetPHOSElectronWeight(Float_t w){ fPHOSElectronWeight = w ; } | |
109 | void SetPHOSChargeWeight(Float_t w){ fPHOSChargeWeight = w ; } | |
110 | void SetPHOSNeutralWeight(Float_t w){ fPHOSNeutralWeight = w ; } | |
111 | ||
a33b76ed | 112 | // void UsePHOSPIDWeightFormula(Bool_t par) { fPHOSWeightFormula = par; } |
113 | // void SetPHOSPhotonWeightFormula(TFormula * const photon) { fPHOSPhotonWeightFormula = photon; } | |
114 | // void SetPHOSPi0WeightFormula(TFormula * const pi0) { fPHOSPi0WeightFormula = pi0; } | |
477d6cee | 115 | |
116 | //PID bits setters and getters | |
5ae09196 | 117 | |
f2ccb5b8 | 118 | Bool_t IsTrackMatched(const AliVCluster * cluster, const AliCalorimeterUtils* cu) const ; |
5ae09196 | 119 | |
477d6cee | 120 | void SetDispersionCut(Float_t dcut ) {fDispCut = dcut; } |
121 | Float_t GetDispersionCut() const {return fDispCut ;} | |
122 | ||
123 | void SetTOFCut(Float_t tcut ) {fTOFCut = tcut; } | |
124 | Float_t GetTOFCut() const {return fTOFCut ;} | |
125 | ||
126 | void SetDebug(Int_t deb) {fDebug=deb;} | |
127 | Int_t GetDebug() const {return fDebug;} | |
128 | ||
bdd2a262 | 129 | void SwitchOnBayesianRecalculation() {fRecalculateBayesian = kTRUE ; } |
130 | void SwitchOffBayesianRecalculation() {fRecalculateBayesian = kFALSE; } | |
131 | enum eventType{kLow,kHigh}; | |
132 | void SetLowParticleFlux() {fParticleFlux = kLow;} | |
133 | void SetHighParticleFlux() {fParticleFlux = kHigh;} | |
134 | ||
f2ccb5b8 | 135 | private: |
477d6cee | 136 | |
137 | Float_t fEMCALPhotonWeight; //Bayesian PID weight for photons in EMCAL | |
138 | Float_t fEMCALPi0Weight; //Bayesian PID weight for pi0 in EMCAL | |
139 | Float_t fEMCALElectronWeight; //Bayesian PID weight for electrons in EMCAL | |
140 | Float_t fEMCALChargeWeight; //Bayesian PID weight for charged hadrons in EMCAL | |
141 | Float_t fEMCALNeutralWeight; //Bayesian PID weight for neutral hadrons in EMCAL | |
142 | Float_t fPHOSPhotonWeight; //Bayesian PID weight for photons in PHOS | |
143 | Float_t fPHOSPi0Weight; //Bayesian PID weight for pi0 in PHOS | |
144 | Float_t fPHOSElectronWeight; //Bayesian PID weight for electrons in PHOS | |
145 | Float_t fPHOSChargeWeight; //Bayesian PID weight for charged hadrons in PHOS | |
146 | Float_t fPHOSNeutralWeight; //Bayesian PID weight for neutral hadrons in PHOS | |
147 | ||
a33b76ed | 148 | // Bool_t fPHOSWeightFormula ; //Use parametrized weight threshold, function of energy |
149 | // TFormula * fPHOSPhotonWeightFormula ; //Formula for photon weight | |
150 | // TFormula * fPHOSPi0WeightFormula ; //Formula for pi0 weight | |
477d6cee | 151 | |
152 | Float_t fDispCut; //Cut on dispersion, used in PID evaluation | |
153 | Float_t fTOFCut; //Cut on TOF, used in PID evaluation | |
154 | ||
155 | Int_t fDebug; //Debug level | |
bdd2a262 | 156 | |
f2ccb5b8 | 157 | Bool_t fRecalculateBayesian; // Recalculate PID bayesian or use simple PID? |
158 | Int_t fParticleFlux; // Particle flux for setting PID parameters | |
159 | AliEMCALPIDUtils * fEMCALPIDUtils; // Pointer to EMCALPID to redo the PID Bayesian calculation | |
bdd2a262 | 160 | |
161 | ClassDef(AliCaloPID,4) | |
e5dbdaf0 | 162 | } ; |
1c5acb87 | 163 | |
164 | ||
165 | #endif //ALICALOPID_H | |
166 | ||
167 | ||
168 |