]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGHF/hfe/AliAnalysisTaskEMCalHFEpA.h
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Why the h*ll do we make a remote commit when pulling?
[u/mrichter/AliRoot.git] / PWGHF / hfe / AliAnalysisTaskEMCalHFEpA.h
CommitLineData
c852fdae 1#ifndef AliAnalysisTaskEMCalHFEpA_cxx
2#define AliAnalysisTaskEMCalHFEpA_cxx
3
4/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
6
9a1f999b 7 ////////////////////////////////////////////////////////////////////////
8 // //
9 // Task for Heavy-flavour electron analysis in pPb collisions //
10 // (+ Electron-Hadron Jetlike Azimuthal Correlation) //
11 // //
12 // version: September 12th, 2013. //
13 // //
14 // Authors //
15 // Elienos Pereira de Oliveira Filho (epereira@cern.ch) //
16 // Cristiane Jahnke (cristiane.jahnke@cern.ch) //
17 // //
18 ////////////////////////////////////////////////////////////////////////
c852fdae 19
20class TH1F;
21class TH2F;
22class AliESDEvent;
23class AliESDtrackCuts;
24class AliESDtrack;
25class AliHFEcontainer;
26class AliHFEcuts;
27class AliHFEpid;
28class AliHFEpidQAmanager;
29class AliCFManager;
30class AliPIDResponse;
31class AliCentrality;
32class AliAODEvent;
33class AliVEvent;
34class AliAODMCHeader;
35class AliSelectNonHFE;
36class AliEventPoolManager;
37class AliEventPool;
38class TObjArray;
39
9a1f999b 40 //______________________________________________________________________
41 //Library
c852fdae 42#include "AliAnalysisTaskSE.h"
43#include "AliHFEpid.h"
44#include "AliLog.h"
9a1f999b 45 //______________________________________________________________________
c852fdae 46
9a1f999b 47 //______________________________________________________________________
c852fdae 48class AliAnalysisTaskEMCalHFEpA : public AliAnalysisTaskSE
49{
9a1f999b 50 //______________________________________________________________________
51public:
c852fdae 52 AliAnalysisTaskEMCalHFEpA();
53 AliAnalysisTaskEMCalHFEpA(const char *name);
54 virtual ~AliAnalysisTaskEMCalHFEpA();
9a1f999b 55
c852fdae 56 virtual void UserCreateOutputObjects();
57 virtual void UserExec(Option_t *option);
58 virtual void Terminate(Option_t *);
9a1f999b 59
60 //Setters
c852fdae 61 void SetHFECuts(AliHFEcuts * const cuts) {fCuts = cuts;};
62 void SetRejectKinkMother(Bool_t rejectKinkMother = kFALSE) {fRejectKinkMother = rejectKinkMother;};
a89620c8 63 void SetCorrelationAnalysis(Bool_t CorrelationFlag=kTRUE) {fCorrelationFlag = CorrelationFlag;};
c852fdae 64 void SetMCanalysis() {fIsMC = kTRUE;};
65 void SetCentrality(Double_t CentralityMin, Double_t CentralityMax) { fCentralityMin = CentralityMin; fCentralityMax = CentralityMax; fHasCentralitySelection = kTRUE; };
66 void SetAODanalysis(Bool_t IsAOD) {fIsAOD = IsAOD;};
67 void SetEventMixing(Bool_t EventMixingFlag) { fEventMixingFlag = EventMixingFlag;};
68 void SetNonHFEmassCut(Double_t MassCut) { fMassCut = MassCut; fMassCutFlag = kTRUE;};
69 void SetEtaCut(Double_t EtaCutMin,Double_t EtaCutMax ) { fEtaCutMin = EtaCutMin; fEtaCutMax = EtaCutMax; };
a6f21076 70
71 void SetdPhidEtaCut(Double_t dPhiCut, Double_t dEtaCut ) { fdPhiCut = dPhiCut;fdEtaCut = dEtaCut ;};
72
c852fdae 73 void SetEoverPCut(Double_t EoverPCutMin,Double_t EoverPCutMax ) { fEoverPCutMin = EoverPCutMin; fEoverPCutMax = EoverPCutMax; };
a6f21076 74
75 void SetM02Cut(Double_t M02CutMin,Double_t M02CutMax ) { fM02CutMin = M02CutMin; fM02CutMax = M02CutMax; };
76 void SetM20Cut(Double_t M20CutMin,Double_t M20CutMax ) { fM20CutMin = M20CutMin; fM20CutMax = M20CutMax; };
77
78
c852fdae 79 void SetNonHFEangleCut(Double_t AngleCut) { fAngleCut = AngleCut; fAngleCutFlag = kTRUE;};
80 void SetNonHFEchi2Cut(Double_t Chi2Cut) { fChi2Cut = Chi2Cut; fChi2CutFlag = kTRUE;};
81 void SetNonHFEdcaCut(Double_t DCAcut) { fDCAcut = DCAcut; fDCAcutFlag = kTRUE;};
82 void SetUseEMCal() { fUseEMCal=kTRUE;};
a89620c8 83 void SetUseShowerShapeCut(Bool_t UseShowerShapeCut=kFALSE) { fUseShowerShapeCut=UseShowerShapeCut;};
84 void SetBackground(Bool_t FillBackground=kFALSE) { fFillBackground=FillBackground;};
c852fdae 85 void SetEMCalTriggerEG1() { fEMCEG1=kTRUE; };
86 void SetEMCalTriggerEG2() { fEMCEG2=kTRUE; };
87 void SetCentralityEstimator(Int_t Estimator) { fEstimator=Estimator; }; //0 = V0A, 1 = Other
9a1f999b 88 void SetAdditionalCuts(Double_t PtMinAsso, Int_t TpcNclsAsso) {fPtMinAsso = PtMinAsso; fTpcNclsAsso = TpcNclsAsso;};
2c8e1f6c 89 void SetSPDCutForHadrons() {fAssocWithSPD = kTRUE;};
c852fdae 90
9a1f999b 91 //Getters
c852fdae 92 AliHFEpid *GetPID() const {return fPID;};
9a1f999b 93 //______________________________________________________________________
c852fdae 94
9a1f999b 95 //______________________________________________________________________
96private:
97
98 //Function to process track cuts
c852fdae 99 Bool_t ProcessCutStep(Int_t cutStep, AliVParticle *track);
9a1f999b 100 //Function to process eh analysis
c852fdae 101 void ElectronHadronCorrelation(AliVTrack *track, Int_t trackIndex, AliVParticle *vtrack);
9a1f999b 102 //Function to find non-HFE and fill histos
9a4682da 103 void Background(AliVTrack *track, Int_t trackIndex, AliVParticle *vtrack, Bool_t IsTPConly);
9a1f999b 104 //Selected Hadrons, for mixed event analysis
c852fdae 105 TObjArray* SelectedHadrons();
9a1f999b 106 //DiHadron Correlation Background
c852fdae 107 void DiHadronCorrelation(AliVTrack *track, Int_t trackIndex);
9a1f999b 108 //Find Mothers (Finde HFE and NonHFE from MC information)
c852fdae 109 Bool_t FindMother(Int_t mcIndex);
9a1f999b 110
111 //Flags for specifics analysis
c852fdae 112 Bool_t fCorrelationFlag;
113 Bool_t fIsMC;
114 Bool_t fUseEMCal;
a89620c8 115 Bool_t fUseShowerShapeCut;
116 Bool_t fFillBackground;
2c8e1f6c 117 Bool_t fAssocWithSPD;
9a4682da 118
9a1f999b 119
c852fdae 120 Bool_t fEMCEG1;
121 Bool_t fEMCEG2;
9a1f999b 122
123 //Used in the function FindMother
c852fdae 124 Bool_t fIsHFE1;
125 Bool_t fIsHFE2;
126 Bool_t fIsNonHFE;
127 Bool_t fIsFromD;
128 Bool_t fIsFromB;
129 Bool_t fIsFromPi0;
130 Bool_t fIsFromEta;
131 Bool_t fIsFromGamma;
132
9a1f999b 133 //General variables
c852fdae 134 AliESDEvent *fESD;
135 AliAODEvent *fAOD; /// new
136 AliVEvent *fVevent; /// new
137 AliESDtrackCuts *fPartnerCuts;
138 TList *fOutputList;
139 AliPIDResponse *fPidResponse;
140 AliSelectNonHFE *fNonHFE;
141
9a1f999b 142 //For the case of AOD analysis
c852fdae 143 Bool_t fIsAOD; //flag for AOD analysis
144
9a1f999b 145 //For Centrality Selection
c852fdae 146 AliCentrality *fCentrality;
147 Double_t fCentralityMin;
148 Double_t fCentralityMax;
149 Bool_t fHasCentralitySelection;
150 TH1F *fCentralityHist;
151 TH1F *fCentralityHistPass;
152 Float_t fZvtx;
153 Int_t fEstimator;
154
9a1f999b 155 //EMCal
a89620c8 156
c852fdae 157 AliVCluster *fClus;
9a1f999b 158 //AliESDCaloCluster *fClusESD;
c852fdae 159
9a1f999b 160 //Histograms
c852fdae 161 TH1F *fNevent;
162 TH1F *fPtElec_Inc;
9a1f999b 163
a6f21076 164
165 TH1F *fCharge_n;
166 TH1F *fCharge_p;
167
a89620c8 168 TH2D *fTime;
169 TH2D *fTime2;
170 TH2D *ftimingEle;
171 TH2D *ftimingEle2;
172
c852fdae 173 TH1F *fPtElec_ULS;
174 TH1F *fPtElec_LS;
9a4682da 175 TH1F *fPtElec_ULS2;
176 TH1F *fPtElec_LS2;
c852fdae 177
2c8e1f6c 178 TH1F *fPtElec_ULS_weight;
179 TH1F *fPtElec_LS_weight;
180 TH1F *fPtElec_ULS2_weight;
181 TH1F *fPtElec_LS2_weight;
182
9a1f999b 183 //PID Histograms
2c8e1f6c 184
185 TH2F *fTOF01;
186 TH2F *fTOF02;
187 TH2F *fTOF03;
188
c852fdae 189 TH1F *fpid;
190
191 TH2F **fEoverP_pt;
192 TH2F **fEoverP_tpc;
193
194 TH1F **fTPC_pt;
195 TH2F **fTPC_p;
196
197 TH1F **fTPCnsigma_pt;
198 TH2F **fTPCnsigma_p;
199 TH2F *fTPCnsigma_pt_2D;
a89620c8 200 TH2F *fShowerShapeCut;
9a4682da 201 TH2F *fShowerShapeM02_EoverP;
202 TH2F *fShowerShapeM20_EoverP;
203 TH2F *fShowerShape_ha;
204 TH2F *fShowerShape_ele;
c852fdae 205
206 TH2F *fTPCnsigma_eta;
207 TH2F *fTPCnsigma_phi;
208
9a1f999b 209
c852fdae 210 TH1F **fECluster;
211 TH2F **fEtaPhi;
212 TH1F **fVtxZ;
213 TH1F **fNTracks;
214 TH1F **fNClusters;
215 TH2F **fTPCNcls_EoverP;
aaec953d 216 TH2F **fTPCNcls_pid;
c852fdae 217
218 TH1F **fEta;
219 TH1F **fPhi;
220 TH1F **fR;
221 TH2F **fR_EoverP;
222 TH1F **fNcells;
223 TH2F **fNcells_EoverP;
224 TH1F **fNcells_electrons;
225 TH1F **fNcells_hadrons;
226 TH1F **fECluster_ptbins;
227 TH1F **fEoverP_ptbins;
228 TH1F **fEoverP_wSSCut;
229 TH2F **fM02_EoverP;
230 TH2F **fM20_EoverP;
231 TH2F **fTPCnsigma_eta_electrons;
232 TH2F **fTPCnsigma_eta_hadrons;
233
234 TH2F *fEoverP_pt_pions;
235
236 TH2F *ftpc_p_EoverPcut;
237 TH2F *fnsigma_p_EoverPcut;
238
239 TH2F *fEoverP_pt_pions2;
240 TH2F *fNcells_pt;
241 TH2F *fEoverP_pt_hadrons;
242
9a1f999b 243 //Electron-Hadron Correlation Histograms
c852fdae 244 TH2F **fCEtaPhi_Inc;
245
246 TH2F **fCEtaPhi_ULS;
247 TH2F **fCEtaPhi_LS;
248 TH2F **fCEtaPhi_ULS_NoP;
249 TH2F **fCEtaPhi_LS_NoP;
250
251 TH2F **fCEtaPhi_ULS_Weight;
252 TH2F **fCEtaPhi_LS_Weight;
253 TH2F **fCEtaPhi_ULS_NoP_Weight;
254 TH2F **fCEtaPhi_LS_NoP_Weight;
255
256 TH1F *fInvMass;
257 TH1F *fInvMassBack;
258 TH1F *fDCA;
259 TH1F *fDCABack;
260 TH1F *fOpAngle;
261 TH1F *fOpAngleBack;
262
9a4682da 263 TH1F *fInvMass2;
264 TH1F *fInvMassBack2;
265 TH1F *fDCA2;
266 TH1F *fDCABack2;
267 TH1F *fOpAngle2;
268 TH1F *fOpAngleBack2;
269
c852fdae 270 Double_t fMassCut;
271 Double_t fEtaCutMin;
272 Double_t fEtaCutMax;
a6f21076 273
274 Double_t fdPhiCut;
275 Double_t fdEtaCut;
276
c852fdae 277 Double_t fEoverPCutMin;
278 Double_t fEoverPCutMax;
a6f21076 279 Double_t fM20CutMin;
280 Double_t fM20CutMax;
281 Double_t fM02CutMin;
282 Double_t fM02CutMax;
283
c852fdae 284 Double_t fAngleCut;
285 Double_t fChi2Cut;
286 Double_t fDCAcut;
287 Bool_t fMassCutFlag;
288 Bool_t fAngleCutFlag;
289 Bool_t fChi2CutFlag;
290 Bool_t fDCAcutFlag;
291
9a1f999b 292 //Non-HFE reconstruction efficiency
c852fdae 293 TH1F *fPtBackgroundBeforeReco;
9a4682da 294 TH1F *fPtBackgroundBeforeReco2;
2c8e1f6c 295 TH1F *fPtBackgroundBeforeReco_weight;
296 TH1F *fPtBackgroundBeforeReco2_weight;
297
2235208e 298 TH2F *fpT_m_electron;
299 TH2F *fpT_gm_electron;
300
9a4682da 301 TH1F *fPtBackgroundAfterReco;
302
9a1f999b 303 Double_t fPtMinAsso;
304 Int_t fTpcNclsAsso;
9a4682da 305
9a1f999b 306 //Tracking Efficiency
c852fdae 307 TH1F *fPtMCparticleAll;
9a1f999b 308 TH1F *fPtMCparticleAll_nonPrimary;
309 TH1F *fPtMCparticleAlle_nonPrimary;
2c8e1f6c 310 TH1F *fPtMCparticleAlle_Primary;
9a1f999b 311
c852fdae 312 TH1F *fPtMCparticleReco;
9a1f999b 313 TH1F *fPtMCparticleReco_nonPrimary;
c852fdae 314 TH1F *fPtMCparticleAllHfe1;
315 TH1F *fPtMCparticleRecoHfe1;
316 TH1F *fPtMCparticleAllHfe2;
317 TH1F *fPtMCparticleRecoHfe2;
318 TH1F *fPtMCelectronAfterAll;
9a1f999b 319 TH1F *fPtMCelectronAfterAll_nonPrimary;
2c8e1f6c 320 TH1F *fPtMCelectronAfterAll_Primary;
c852fdae 321
322 TH1F *fPtMCpi0;
9a1f999b 323 TH1F *fPtMCeta;
c852fdae 324
325 TH1F *fPtMC_EMCal_All;
326 TH1F *fPtMC_EMCal_Selected;
327 TH1F *fPtMC_TPC_All;
328 TH1F *fPtMC_TPC_Selected;
329
330 TH1F *fPtMCWithLabel;
331 TH1F *fPtMCWithoutLabel;
332 TH1F *fPtIsPhysicaPrimary;
333
9a1f999b 334 //For the HFE package
c852fdae 335 AliHFEcuts *fCuts; // Cut Collection for HFE
336 AliCFManager *fCFM; // Correction Framework Manager
337 AliHFEpid *fPID; // PID
338 AliHFEpidQAmanager *fPIDqa; // PID QA manager
339
9a1f999b 340 //Others
c852fdae 341 AliStack *fMCstack; //
342 Bool_t fRejectKinkMother; //
343 TParticle *fMCtrack;
344 TParticle *fMCtrackMother;
345 TParticle *fMCtrackGMother;
346 TParticle *fMCtrackGGMother;
347 TParticle *fMCtrackGGGMother;
348 TClonesArray *fMCarray;
349 AliAODMCHeader *fMCheader;
350 AliAODMCParticle *fMCparticle;
351 AliAODMCParticle *fMCparticleMother;
352 AliAODMCParticle *fMCparticleGMother;
353 AliAODMCParticle *fMCparticleGGMother;
354 AliAODMCParticle *fMCparticleGGGMother;
355 AliMCEventHandler *fEventHandler;
356 AliMCEvent *fMCevent;
9a1f999b 357
358 //______________________________________________________________________
359 //Mixed event analysis
c852fdae 360 AliEventPoolManager *fPoolMgr;
361 AliEventPool *fPool;
362 TObjArray *fTracksClone;
363 TObjArray *fTracks;
364
365 TH2F **fCEtaPhi_Inc_EM;
366
367 TH2F **fCEtaPhi_ULS_EM;
368 TH2F **fCEtaPhi_LS_EM;
369
370 TH2F **fCEtaPhi_ULS_Weight_EM;
371 TH2F **fCEtaPhi_LS_Weight_EM;
372
373 TH1F *fPoolNevents;
374
375 Bool_t fEventMixingFlag;
9a1f999b 376 //______________________________________________________________________
377
378 //______________________________________________________________________
379 //Di-hadron correlation
c852fdae 380 TH2F **fCEtaPhi_Inc_DiHadron;
381 TH1F *fPtTrigger_Inc;
9a1f999b 382 //______________________________________________________________________
383
c852fdae 384 AliAnalysisTaskEMCalHFEpA(const AliAnalysisTaskEMCalHFEpA&); // not implemented
385 AliAnalysisTaskEMCalHFEpA& operator=(const AliAnalysisTaskEMCalHFEpA&); // not implemented
9a1f999b 386
c852fdae 387 ClassDef(AliAnalysisTaskEMCalHFEpA, 1); // example of analysis
9a1f999b 388 //______________________________________________________________________
c852fdae 389};
390
9a1f999b 391 ///_________________________________________________________________________________________________
392 ///Class copied from : $ALICE_ROOT/PWGCF/Correlations/DPhi/AliAnalysisTaskLongRangeCorrelations.h
393 ///Author: Christoph Mayer
c852fdae 394class AliEHCParticle : public TObject {
395public:
9a1f999b 396 AliEHCParticle(Double_t eta=0, Double_t phi=0, Double_t pt=0)
c852fdae 397 : fEta(eta), fPhi(phi), fPt(pt) {}
9a1f999b 398 virtual ~AliEHCParticle() {}
399
400 Double_t Eta() const { return fEta; }
401 Double_t Phi() const { return fPhi; }
402 Double_t Pt() const { return fPt; }
403
c852fdae 404protected:
405private:
9a1f999b 406 AliEHCParticle(const AliEHCParticle&);
407 AliEHCParticle& operator=(const AliEHCParticle&);
408
409 Double_t fEta;
410 Double_t fPhi;
411 Double_t fPt;
412
413 ClassDef(AliEHCParticle, 1);
c852fdae 414} ;
9a1f999b 415 ///_________________________________________________________________________________________________
c852fdae 416
417#endif