AliPhysics  608b256 (608b256)
AliAnalysisTaskSOH.h
Go to the documentation of this file.
1 #ifndef ALIANALYSISTASKSOH_H
2 #define ALIANALYSISTASKSOH_H
3 
4 // $Id$
5 
6 class TList;
7 class TH1F;
8 class TH2F;
9 class TH3F;
10 class THnSparse;
11 class TArrayI;
12 class AliESDEvent;
13 class AliMCEvent;
14 class AliESDtrack;
15 class AliESDCaloCluster;
16 class AliESDtrackCuts;
17 
18 #include "AliAnalysisTaskSE.h"
19 
21  public:
23  AliAnalysisTaskSOH(const char *name);
24  virtual ~AliAnalysisTaskSOH();
25 
27  void UserExec(Option_t *option);
28  void Terminate(Option_t *);
29 
34  void SetZvtx(Double_t zvtx) { fZVtxMax = zvtx; }
35 
36  void SetEsdTrackCuts(AliESDtrackCuts *cuts) { fEsdTrackCuts = cuts ; }
37  void SetHybridTrackCuts1(AliESDtrackCuts *cuts) { fHybridTrackCuts1 = cuts ; }
38  void SetHybridTrackCuts2(AliESDtrackCuts *cuts) { fHybridTrackCuts2 = cuts ; }
39 
40  private:
41 
42  AliESDtrack *GetAcceptTrack(AliESDtrack *esdtrack) ;
43  Bool_t IsGoodCluster(AliESDCaloCluster *cluster) ;
44  Bool_t IsGoodMcParticle(AliVParticle* vParticle, Int_t ipart);
45  Bool_t EsdVertexOk() const ;
46  void ProcessTrack() ;
47  void ProcessCluster() ;
48  void ProcessMc() ;
49  void ProcessScaleFactor() ;
50 
52  AliMCEvent *fMC;
53  Double_t fZVtxMax; // Max vertex z cut
54  AliESDtrackCuts *fEsdTrackCuts; // esd track cuts
55  AliESDtrackCuts *fHybridTrackCuts1; // hybrid track cuts
56  AliESDtrackCuts *fHybridTrackCuts2; // hybrid track cuts
64 
66 
67  TH1F *fHEventStat;
68  TH1F *fHScaleFactor;
91  THnSparse *fHClsEoverMcE_All;
92  THnSparse *fHClsEoverMcE_Photon;
93  THnSparse *fHClsEoverMcE_Elec;
94  THnSparse *fHClsEoverMcE_Pion;
111 
112  AliAnalysisTaskSOH(const AliAnalysisTaskSOH&); // not implemented
113  AliAnalysisTaskSOH& operator=(const AliAnalysisTaskSOH&); // not implemented
114 
115  ClassDef(AliAnalysisTaskSOH, 13); // Analysis task Saehanseul Oh
116 };
117 #endif
TArrayI * fClusterIndices
selected track index
Bool_t EsdVertexOk() const
double Double_t
Definition: External.C:58
TH2F * fHistEsub1Pch
cluster E vs. truth photon E
Definition: External.C:260
void Terminate(Option_t *)
TH2F * fHEMCalResponseElec
same as above for pions
TH2F * fHistEsub1PchRat
(subtracted E in 100% HC) vs. total track P, clusters with 2 matching tracks
Definition: External.C:236
TH2F * fHEMCalRecdPhidEtaposEta
same as above with negative truth charge matching
AliESDtrackCuts * fEsdTrackCuts
TObjArray * fClusterArray
cluster with two matched MC track index
AliESDtrackCuts * fHybridTrackCuts1
TH2F * fHPhotonEdiff70HC
(truth E - calculated E in 100% HC)/truth E vs. truth E with photon
TH3F * fHDetGenPion_p[3]
secondary pion tracks pt, phi, eta spectrum
TH3F * fHDetRecPion_rmInj_m[3]
minus charged generated detector level particle(pion) without injected signal pt, phi...
void SetMcProcess(Bool_t p)
TH2F * fHPhotonEdiff30HC
(truth E - calculated E in 70% HC)/truth E vs. truth E with photon
TH2F * fHPhotonEdiff0HC
(truth E - calculated E in 30% HC)/truth E vs. truth E with photon
void SetClusterProcess(Bool_t p)
TH1F * fHScaleFactor
statistics histo
TH2F * fHEMCalResponseProton
same as above for electrons
Bool_t IsGoodMcParticle(AliVParticle *vParticle, Int_t ipart)
void SetEsdTrackCuts(AliESDtrackCuts *cuts)
TH2F * fHPhotonEVsClsE
(truth E - cluster E)/truth E vs. truth E with photon
void SetSFProcess(Bool_t p)
AliAnalysisTaskSOH & operator=(const AliAnalysisTaskSOH &)
int Int_t
Definition: External.C:63
TH2F * fHPhotonEdiff100HC
same as above for negative eta
TH3F * fHDetRecPion_m[3]
minus pion mc detector level pt, phi, eta spectrum
TH3F * fHDetRecSecPion
secondary pion tracks pt, phi, eta spectrum
TH3F * fHDetRecFakePion
fake pion tracks pt, phi, eta spectrum
Bool_t fMcProcess
selected cluster array
Double_t fZVtxMax
mv event
THnSparse * fHClsEoverMcE_Pion
above for electron
THnSparse * fHClsEoverMcE_Elec
above for photon
TH2F * fHEMCalRecdPhidEtanegEta
same as above for positive eta
AliESDtrack * GetAcceptTrack(AliESDtrack *esdtrack)
TH2F * fHEMCalResponsePion
(cluster energy over reconstructed track p) vs. track pt
TH2F * fHEMCalRecdPhidEtaM_Truth
same as above with positive truth charge matching
TH3F * fHParGenPion_rmInj_p
minus pion mc truth pt, phi, eta spectrum
TH3F * fHDetGenSecPion
fake pion tracks pt, phi, eta spectrum
TH2F * fHEMCalRecdPhidEta_Truth
same as above for negative charge tracks
TH3F * fHParGenPion_rmInj_m
plus charged mc truth(pion) without injected signal pt, phi, eta spectrum
TH2F * fHEOverPVsPt
scale factor with 100% HC spectrum
TH1F * fHEventStat
output list
TH2F * fHEMCalRecdPhidEtaP
(EMCal cluster phi - track phi) vs. (EMCal cluster eta - track eta)
TH3F * fHParGenPion_m
plus pion mc truth pt, phi, eta spectrum
TH3F * fHDetRecPion_rmInj_p[3]
plus charged generated detector level particle(pion) without injected signal pt, phi, eta spectrum
void SetHybridTrackCuts2(AliESDtrackCuts *cuts)
void UserExec(Option_t *option)
TH3F * fHDetGenPion_rmInj_m[3]
plus charged reconstructed detector level pion+ track without injected signal pt, phi...
void SetHybridTrackCuts1(AliESDtrackCuts *cuts)
TH3F * fHDetGenPion_m[3]
plus pion reconstructed detector level pt, phi, eta spectrum
void SetZvtx(Double_t zvtx)
TH3F * fHDetGenPion_rmInj_p[3]
minus pion reconstructed detector level pt, phi, eta spectrum
const char Option_t
Definition: External.C:48
TH2F * fHEMCalRecdPhidEtaM
same as above for positive charge tracks
Bool_t IsGoodCluster(AliESDCaloCluster *cluster)
TH3F * fHDetGenFakePion
minus charged mc truth(pion) without injected signal pt, phi, eta spectrum
TH1F * fHScaleFactor100HC
scale factor spectrum
TH2F * fHEMCalRecdPhidEta
same as above for protons
bool Bool_t
Definition: External.C:53
TH3F * fHDetRecPion_p[3]
plus pion mc detector level pt, phi, eta spectrum
TH2F * fHEMCalRecdPhidEtaP_Truth
same as above with mc truth matching
AliESDtrackCuts * fHybridTrackCuts2
THnSparse * fHClsEoverMcE_All
(subtracted E in 100% HC)/total track P vs. total track P, clusters with 2 matching tracks ...
THnSparse * fHClsEoverMcE_Photon
cluster E/MC particle E, cluster with only one matching particle
TH3F * fHParGenPion_p
above for pion
TH2F * fHistEsub2PchRat
(subtracted E in 100% HC)/total track P vs. total track P, clusters with 1 matching track ...
TH2F * fHistEsub2Pch
(subtracted E in 100% HC) vs. total track P, clusters with 1 matching track
AliMCEvent * fMC
esd event
void SetTrackProcess(Bool_t p)