1 #ifndef ALIANAELECTRON_H
\r
2 #define ALIANAELECTRON_H
\r
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
\r
4 * See cxx source for full Copyright notice */
\r
7 //_________________________________________________________________________
\r
9 // Class for the electron identification.
\r
10 // Clusters from EMCAL matched to tracks are selected
\r
11 // and kept in the AOD. Few histograms produced.
\r
14 //-- Author: J.L. Klay (Cal Poly)
\r
16 // --- ROOT system ---
\r
22 // --- ANALYSIS system ---
\r
23 #include "AliAnaPartCorrBaseClass.h"
\r
25 class AliAODMCParticle;
\r
26 class AliCaloTrackReader;
\r
30 class AliAnaElectron : public AliAnaPartCorrBaseClass {
\r
34 AliAnaElectron() ; // default ctor
\r
35 AliAnaElectron(const AliAnaElectron & g) ; // cpy ctor
\r
36 AliAnaElectron & operator = (const AliAnaElectron & g) ;//cpy assignment
\r
37 virtual ~AliAnaElectron() ; //virtual dtor
\r
39 TList * GetCreateOutputObjects();
\r
43 void MakeAnalysisFillAOD() ;
\r
45 void MakeAnalysisFillHistograms() ;
\r
48 Int_t GetDVMBtag(AliAODTrack * tr); //returns # tracks from secvtx
\r
50 //Temporary local method to get DCA because AliAODTrack is stupid
\r
51 Bool_t GetDCA(const AliAODTrack* tr,Double_t imp[2], Double_t cov[3]);
\r
53 Bool_t PhotonicPrim(const AliAODPWG4Particle* part); //check with track list
\r
54 Bool_t PhotonicV0(Int_t trackId); //check with V0 list
\r
56 //check if track has been flagged as a non-photonic or DVM electron
\r
57 //used with the jet tracks to tag bjets
\r
58 Bool_t CheckTrack(const AliAODTrack* track,const char* type);
\r
59 Bool_t IsMcBJet(Double_t x, Double_t y);
\r
60 Bool_t IsMcDJet(Double_t x, Double_t y);
\r
62 void Print(const Option_t * opt)const;
\r
64 TString GetCalorimeter() const {return fCalorimeter ; }
\r
65 Double_t GetpOverEmin() const {return fpOverEmin ; }
\r
66 Double_t GetpOverEmax() const {return fpOverEmax ; }
\r
67 Bool_t GetWriteNtuple() const {return fWriteNtuple ; }
\r
69 Double_t GetDrCut() const { return fDrCut; }
\r
70 Double_t GetPairDcaCut() const { return fPairDcaCut; }
\r
71 Double_t GetDecayLenCut() const { return fDecayLenCut; }
\r
72 Double_t GetImpactCut() const { return fImpactCut; }
\r
73 Double_t GetAssocPtCut() const { return fAssocPtCut; }
\r
74 Double_t GetMassCut() const { return fMassCut; }
\r
75 Double_t GetSdcaCut() const { return fSdcaCut; }
\r
76 Int_t GetITSCut() const { return fITSCut; }
\r
77 Int_t GetNTagTrackCut() const { return fNTagTrkCut; }
\r
78 Double_t GetIPSigCut() const { return fIPSigCut; }
\r
79 Double_t GetMinClusEne() const { return fMinClusEne; }
\r
81 void SetCalorimeter(TString det) {fCalorimeter = det ; }
\r
82 void SetpOverEmin(Double_t min) {fpOverEmin = min ; }
\r
83 void SetpOverEmax(Double_t max) {fpOverEmax = max ; }
\r
84 void SetResidualCut(Double_t cut) {fResidualCut = cut ; }
\r
85 void SetWriteNtuple(Bool_t val) {fWriteNtuple = val ; }
\r
87 void SetDrCut(Double_t dr) { fDrCut = dr; }
\r
88 void SetPairDcaCut(Double_t pdca) { fPairDcaCut = pdca; }
\r
89 void SetDecayLenCut(Double_t dlen) { fDecayLenCut = dlen; }
\r
90 void SetImpactCut(Double_t imp) { fImpactCut = imp; }
\r
91 void SetAssocPtCut(Double_t pt) { fAssocPtCut = pt; }
\r
92 void SetMassCut(Double_t mass) { fMassCut = mass; }
\r
93 void SetSdcaCut(Double_t sdca) { fSdcaCut = sdca; }
\r
94 void SetITSCut(Int_t its) { fITSCut = its; }
\r
95 void SetNTagTrackCut(Int_t ntr) { fNTagTrkCut = ntr; }
\r
96 void SetIPSigCut(Double_t ips) { fIPSigCut = ips; }
\r
97 void SetMinClusEne(Double_t ene) { fMinClusEne = ene; }
\r
99 void InitParameters();
\r
101 void Terminate(TList * outputList);
\r
102 void ReadHistograms(TList * outputList); //Fill histograms with
\r
103 //histograms in ouput list,
\r
104 //needed in Terminate.
\r
106 //For DVM B-tag method
\r
107 Double_t ComputeSignDca(AliAODTrack *track, AliAODTrack *track2 , float cut1);
\r
108 //the 2 following functions are internal methods of the b-tagging
\r
109 //based on transverse impact parameter
\r
110 Double_t GetIPSignificance(AliAODTrack *tr, Double_t jetPhi);
\r
111 void GetImpactParamVect(Double_t Pxy[2], Double_t t[2], Double_t Vxy[2], Double_t ip[2]);
\r
112 //For determining origin of electron
\r
113 Int_t GetMCSource(Int_t mctag);
\r
115 //Need a clean way to get the MC info. An AliAODMCParticle object
\r
116 //is returned from whichever source we are operating on
\r
117 AliAODMCParticle* GetMCParticle(Int_t part);
\r
118 //Get MC B Parent pt
\r
119 Double_t GetBParentPt(Int_t label);
\r
120 //Get Number of particles in AliAODMCParticle array, if it exists
\r
121 Int_t GetNumAODMCParticles();
\r
124 TString fCalorimeter; //! Which detector? EMCAL or PHOS
\r
125 Double_t fpOverEmin; //! Minimum p/E value for Electrons
\r
126 Double_t fpOverEmax; //! Maximum p/E value for Electrons
\r
127 Double_t fResidualCut; //! Track-cluster matching distance
\r
128 Double_t fMinClusEne; //! Min clus energy for matching
\r
131 Double_t fDrCut; //max dR
\r
132 Double_t fPairDcaCut; //max pair-DCA
\r
133 Double_t fDecayLenCut; //max 3d-decaylength
\r
134 Double_t fImpactCut; //max track impact param
\r
135 Double_t fAssocPtCut; //min associated pt
\r
136 Double_t fMassCut; //min Minv cut
\r
137 Double_t fSdcaCut; //min signDca
\r
138 Int_t fITSCut; //min ITS hits (both)
\r
140 Int_t fNTagTrkCut; //min number of tracks required for IP sig tag
\r
141 Double_t fIPSigCut; //min IP significance cut
\r
143 Double_t fJetEtaCut; //max eta for jets
\r
144 Double_t fJetPhiMin; //min phi for jets
\r
145 Double_t fJetPhiMax; //max phi for jets
\r
147 Bool_t fWriteNtuple; //flag for filling ntuple or not
\r
149 ///////////////////////////////////////
\r
150 //Output histograms and Ntuples
\r
152 ///////////////////////////////////////
\r
153 //RC = RECO only - these histos will be filled using only reco
\r
157 TH1F * fhImpactXY; //! XY impact parameter of all tracks to primary vertex
\r
158 TH1F * fhRefMult; //! refmult (tracks with |eta| < 0.5)
\r
159 TH1F * fhRefMult2; //! refmult2 (tracks with |eta| < 0.5 & impXY,impZ < 1.0)
\r
162 TH3F *fh3pOverE; //! p/E for track-cluster matches vs pt vs mult
\r
163 TH3F *fh3EOverp; //! E/p for track-cluster matches vs pt vs mult
\r
164 TH3F *fh3pOverE2; //! p/E for track-cluster matches vs pt vs mult
\r
165 TH3F *fh3EOverp2; //! E/p for track-cluster matches vs pt vs mult
\r
166 TH3F *fh3pOverE3; //! p/E for track-cluster matches vs pt vs mult
\r
167 TH3F *fh3EOverp3; //! E/p for track-cluster matches vs pt vs mult
\r
170 TH2F *fh2pOverE; //! p/E for track-cluster matches vs pt vs mult
\r
171 TH2F *fh2EOverp; //! E/p for track-cluster matches vs pt vs mult
\r
172 TH2F *fh2pOverE2; //! p/E for track-cluster matches vs pt vs mult
\r
173 TH2F *fh2EOverp2; //! E/p for track-cluster matches vs pt vs mult
\r
176 TH1F *fh1dR; //! distance between projected track and cluster
\r
177 TH2F *fh2EledEdx; //! dE/dx vs. momentum for electron candidates
\r
178 TH2F *fh2MatchdEdx; //! dE/dx vs. momentum for all matches
\r
179 TH2F *fh2dEtadPhi; //! DeltaEta vs. DeltaPhi of all track/cluster pairs
\r
180 TH2F *fh2dEtadPhiMatched; //! DeltaEta vs. DeltaPhi of matched track/cluster pairs
\r
181 TH2F *fh2dEtadPhiUnmatched; //! DeltaEta vs. DeltaPhi of unmatched track/cluster pairs
\r
183 TH2F* fh2TrackPVsClusterE; //!track momentum vs. cluster energy
\r
184 TH2F* fh2TrackPtVsClusterE; //!track pt vs. cluster energy
\r
185 TH2F* fh2TrackPhiVsClusterPhi; //!track phi vs. cluster phi
\r
186 TH2F* fh2TrackEtaVsClusterEta; //!track eta vs. cluster eta
\r
188 //Photonic Electron checks
\r
189 TH1F* fh1OpeningAngle; //!opening angle between pairs of photon candidates
\r
190 TH1F* fh1MinvPhoton; //!invariant mass distribution of electron pairs
\r
192 //Reconstructed electrons
\r
193 TH1F * fhPtElectron; //! Number of identified electron vs transverse momentum
\r
194 TH2F * fhPhiElectron; //! Azimuthal angle of identified electron vs transverse momentum
\r
195 TH2F * fhEtaElectron; //! Pseudorapidity of identified electron vs tranvserse momentum
\r
197 TH1F * fhPtNPE; //! Number of non-photonic electron vs transverse momentum
\r
198 TH2F * fhPhiNPE; //! Azimuthal angle of non-photonic electron vs transverse momentum
\r
199 TH2F * fhEtaNPE; //! Pseudorapidity of non-photonic electron vs tranvserse momentum
\r
201 TH1F * fhPtPE; //! Number of photonic electron vs transverse momentum
\r
202 TH2F * fhPhiPE; //! Azimuthal angle of photonic electron vs transverse momentum
\r
203 TH2F * fhEtaPE; //! Pseudorapidity of photonic electron vs tranvserse momentum
\r
206 TH2F * fhDVMBtagCut1; //! DVM B-tagging result for cut1 (minv>1.0)
\r
207 TH2F * fhDVMBtagCut2; //! DVM B-tagging result for cut2 (minv>1.5)
\r
208 TH2F * fhDVMBtagCut3; //! DVM B-tagging result for cut3 (minv>1.8)
\r
209 TH2F * fhDVMBtagQA1; //! DVM B-tagging : QA of pairDca vs decaylength
\r
210 TH2F * fhDVMBtagQA2; //! DVM B-tagging : QA of signDca vs mass
\r
211 TH1F * fhDVMBtagQA3; //! DVM B-tagging : QA number of ITS clusters
\r
212 TH1F * fhDVMBtagQA4; //! DVM B-tagging : QA prim vtx impXY
\r
213 TH1F * fhDVMBtagQA5; //! DVM B-tagging : QA prim vtx impZ
\r
215 TH1F * fhIPSigBtagQA1; //! IPSig B-tagging : QA of # tag tracks
\r
216 TH1F * fhIPSigBtagQA2; //! IPSig B-tagging : QA of IP sig
\r
217 TH1F * fhTagJetPt1x4; //! IPSig B-tagging : result for (1 track, ipSignif>4)
\r
218 TH1F * fhTagJetPt2x3; //! IPSig B-tagging : result for (2 track, ipSignif>3)
\r
219 TH1F * fhTagJetPt3x2; //! IPSig B-tagging : result for (3 track, ipSignif>2)
\r
220 TH1F * fhePlusTagJetPt1x4; //! IPSig B-tagging : eJet + result for (1 track, ipSignif>4)
\r
221 TH1F * fhePlusTagJetPt2x3; //! IPSig B-tagging : eJet + result for (2 track, ipSignif>3)
\r
222 TH1F * fhePlusTagJetPt3x2; //! IPSig B-tagging : eJet + result for (3 track, ipSignif>2)
\r
225 TH2F* fhJetType; //! How many of each tag were found vs jet pt
\r
226 TH2F* fhLeadJetType; //! How many leading of each tag were found vs jet pt
\r
227 TH2F* fhBJetXsiFF; //! B-tagged jet FF with xsi = log(pt_Jet/pt_Track)
\r
228 TH2F* fhBJetPtFF; //! B-tagged jet FF with pt_Track
\r
229 TH2F* fhBJetEtaPhi; //! B-tagged jet eta-phi distribution
\r
230 TH2F* fhNonBJetXsiFF; //! Non b-tagged jet FF with xsi = log(pt_Jet/pt_Track)
\r
231 TH2F* fhNonBJetPtFF; //! Non b-tagged jet FF with pt_Track
\r
232 TH2F* fhNonBJetEtaPhi; //! Non b-tagged jet eta-phi distribution
\r
234 ///////////////////////////////////////////////////////////////////
\r
235 //MC = From here down, the histograms use MC information, so they will
\r
236 //only be filled in simulations
\r
237 TNtuple* fEleNtuple; //! testing ntuple
\r
239 TH2F * fhPhiConversion; //! Azimuthal angle of conversion electron vs transverse momentum
\r
240 TH2F * fhEtaConversion; //! Pseudorapidity of conversion electron vs tranvserse momentum
\r
242 //Histograms for comparison to tracking detectors
\r
243 TH2F* fhPtTrack; //!Pt distribution of reco tracks with MC-ID
\r
244 TH2F* fhPtHadron; //!Pt distribution of reco charged hadrons
\r
245 //!(pi,k,p) in EMCAL acceptance
\r
246 TH2F* fhPtNPEleTPC; //!Pt distribution of non-photonic reco electrons using
\r
247 //!just TPC dEdx info in EMCAL acceptance
\r
248 TH2F* fhPtNPEleTPCTRD; //!Pt distribution of non-photonic reco electrons using
\r
249 //!pid info from tracking detectors only in EMCAL acceptance
\r
250 TH2F* fhPtNPEleTTE; //!Pt distribution of non-photonic reco
\r
251 //!electrons using pid info from TPC+TRD+EMCAL
\r
252 //!in EMCAL acceptance
\r
253 TH2F* fhPtNPEleEMCAL; //!Pt distribution of non-photonic reco
\r
254 //!electrons using EMCAL only
\r
255 //!in EMCAL acceptance
\r
257 TH2F* fhPtNPEBHadron; //!correlate our best reconstructed
\r
258 //b-electrons with the b-hadron momentum
\r
260 //For computing efficiency of IPSIG tag
\r
261 //these require that an MC b-Ancestor is present in the jet
\r
262 TH1F * fhBJetPt1x4; //! IPSig B-tagging : result for (1 track, ipSignif>4)
\r
263 TH1F * fhBJetPt2x3; //! IPSig B-tagging : result for (2 track, ipSignif>3)
\r
264 TH1F * fhBJetPt3x2; //! IPSig B-tagging : result for (3 track, ipSignif>2)
\r
266 TH1F * fhFakeJetPt1x4; //! IPSig B-tagging : fake result for (1 track, ipSignif>4)
\r
267 TH1F * fhFakeJetPt2x3; //! IPSig B-tagging : fake result for (2 track, ipSignif>3)
\r
268 TH1F * fhFakeJetPt3x2; //! IPSig B-tagging : fake result for (3 track, ipSignif>2)
\r
270 TH2F* fhDVMJet; //! DVM jet algo check
\r
272 ////////////////////////////
\r
273 //MC Only Rate histograms
\r
275 TNtuple *fMCEleNtuple; //! Ntuple of MC electrons
\r
277 TH2F* fhMCBJetElePt; //! Pt of B-Jet vs pt of electron
\r
278 TH2F* fhMCBHadronElePt; //! Pt of B-hadrons vs pt of electron
\r
279 TH1F* fhPtMCHadron; //! Pt distribution of MC charged hadrons (pi,k,p) in EMCAL acceptance
\r
280 TH2F* fhPtMCElectron; //! Pt distribution of MC electrons from various sources in EMCAL
\r
281 TH2F* fhMCXYConversion; //! XY distribution of conversion electrons
\r
282 TH2F* fhMCRadPtConversion; //! Radius vs. pT distribution of conversion electrons
\r
284 ClassDef(AliAnaElectron,12)
\r
289 #endif//ALIANAELECTRON_H
\r