AliPhysics  2853087 (2853087)
AliAnalysisTaskMaterialHistos.h
Go to the documentation of this file.
1 #ifndef AliAnalysisTaskMaterialHistos_cxx
2 #define AliAnalysisTaskMaterialHistos_cxx
3 
4 #include "AliAnalysisTaskSE.h"
5 #include "AliConversionPhotonBase.h"
6 #include "TH1.h"
7 #include "TH2.h"
8 #include "TH3.h"
9 #include "TTreeStream.h"
10 #include "AliLog.h"
11 #include <vector>
12 #include "AliV0ReaderV1.h"
13 #include "AliConversionPhotonCuts.h"
14 #include "AliConvEventCuts.h"
15 #include "TList.h"
16 #include "TClonesArray.h"
17 
18 
19 using namespace std;
20 
21 
22 class AliAnalysisTaskMaterialHistos : public AliAnalysisTaskSE{
23 
24  public:
25 
26  AliAnalysisTaskMaterialHistos();
27  AliAnalysisTaskMaterialHistos(const char *name);
28  virtual ~AliAnalysisTaskMaterialHistos();
29 
30  virtual void UserCreateOutputObjects();
31  virtual Bool_t Notify();
32  virtual void UserExec(Option_t *option);
33  virtual void SetLogBinningXTH2(TH2* histoRebin);
34  virtual void Terminate(Option_t *);
35  void SetIsHeavyIon(Int_t flag) {fIsHeavyIon = flag;}
36  void SetIsMC(Int_t isMC) {fIsMC=isMC;}
37  void SetV0Reader(AliV0ReaderV1 *v0Reader) {fV0Reader=v0Reader;}
38  void SetV0ReaderName(TString name) {fV0ReaderName=name; return;}
39 
40  void SetEventCutList(Int_t nCuts, TList *CutArray) {fnCuts = nCuts;
41  fEventCutArray = CutArray;}
42  void SetConversionCutList(Int_t nCuts, TList *CutArray) {fnCuts = nCuts;
43  fConversionCutArray = CutArray;}
44 
45 
46 
47  /* void SetDoTreesForMaterial(Bool_t flag) { fDoTreesForMaterial = flag ;} */
48  /* void SetDoHistosForMaterial(Bool_t flag) { fDoHistosForMaterial = flag ;} */
49  private:
50 
51  void ProcessPhotons();
52  void ProcessMCPhotons();
53  void FillMCHistograms(Int_t eventPos);
54  Int_t CountTracks09();
55  Int_t CountTracks0914();
56 
57  AliV0ReaderV1 *fV0Reader; //
58  TString fV0ReaderName;
59  TClonesArray *fConversionGammas; // Reconstructed Photons;
60  TList* fGammaCandidates; //
61  TList* fConversionCutArray;
62  TList* fEventCutArray; //
63  TList** fCutFolder; //
64  TList** fESDList; //
65  TList** fTrueList; //
66  TList** fMCList; //
67  TList *fOutputList; //
68  TList *fAllMCGammaList; //
69  TList *fAllMCConvGammaList; //
70  Float_t fPrimVtxZ; //
71  Int_t fNContrVtx; //
72  Int_t fNESDtracksEta09; //
73  Int_t fNESDtracksEta0914; //
74  Int_t fNESDtracksEta14; //
75  Float_t fGammaMCPt; //
76  Float_t fGammaMCTheta; //
77  Float_t fGammaMCConvPt; //
78  Float_t fGammaMCConvTheta; //
79  Float_t fGammaPt; //
80  Float_t fGammaTheta; //
81  Float_t fGammaChi2NDF; //
82  UChar_t fKind; //
83  Int_t fIsHeavyIon; //
84  Int_t fIsMC; //
85  AliVEvent* fInputEvent; //
86  AliMCEvent* fMCEvent; //
87  Int_t fnCuts; //
88  Int_t fiCut; //
89  TH1F** hNEvents;
90  TH1F** hNGoodESDTracksEta09;
91  TH1F** hNGoodESDTracksEta14;
92  TH1F** hNGoodESDTracksEta09_14;
93  TH2F** hESDConversionRPhi;
94  TH2F** hESDConversionRZ;
95  TH2F** hESDConversionRPt;
96  TH2F** hESDConversionREta;
97  TH1F** hESDConversionDCA;
98  TH1F** hESDConversionPsiPair;
99  TH1F** hESDConversionChi2;
100  TH1F** hESDConversionMass;
101  TH1F** hESDConversionRRejSmall;
102  TH1F** hESDConversionRRejLarge;
103  TH2F** hESDConversionAsymP;
104  TH2F** hElectronRdEdx;
105  TH2F** hElectronRNSigmadEdx;
106  TH2F** hPositronRdEdx;
107  TH2F** hPositronRNSigmadEdx;
108 
109  TH2F** hMCConversionRPhi;
110  TH2F** hMCConversionRPt;
111  TH2F** hMCConversionREta;
112  TH1F** hMCConversionRRejSmall;
113  TH1F** hMCConversionRRejLarge;
114  TH1F** hMCAllGammaPt;
115 
116 
117  TH2F** hMCTrueConversionRPhi;
118  TH2F** hMCTrueConversionRZ;
119  TH2F** hMCTrueConversionRPt;
120  TH2F** hMCTrueConversionRPtMCRPt;
121  TH2F** hMCTrueConversionREta;
122  TH1F** hMCTrueConversionDCA;
123  TH1F** hMCTrueConversionPsiPair;
124  TH1F** hMCTrueConversionChi2;
125  TH1F** hMCTrueConversionMass;
126  TH2F** hMCTrueConversionAsymP;
127  TH1F** hMCTrueConversionRRejSmall;
128  TH1F** hMCTrueConversionRRejLarge;
129 
130  TH2F** hMCTruePi0DalConversionRPt;
131  TH1F** hMCTruePi0DalConversionEta;
132  TH2F** hMCTrueEtaDalConversionRPt;
133  TH1F** hMCTrueEtaDalConversionEta;
134  TH2F** hMCTrueCombinatorialConversionRPt;
135  TH1F** hMCTrueCombinatorialConversionEta;
136  //Bool_t fDoHistosForMaterial; // flag for using Trees for Material Budget evaluation
137 
138  AliAnalysisTaskMaterialHistos(const AliAnalysisTaskMaterialHistos&); // not implemented
139  AliAnalysisTaskMaterialHistos& operator=(const AliAnalysisTaskMaterialHistos&); // not implemented
140 
141 
142  ClassDef(AliAnalysisTaskMaterialHistos, 14);
143 };
144 
145 #endif
AliConversionPhotonCuts.h
AliAnalysisTaskMaterialHistos::SetV0Reader
void SetV0Reader(AliV0ReaderV1 *v0Reader)
Definition: AliAnalysisTaskMaterialHistos.h:37
TH2F
Definition: External.C:236
AliAnalysisTaskMaterialHistos::SetEventCutList
void SetEventCutList(Int_t nCuts, TList *CutArray)
Definition: AliAnalysisTaskMaterialHistos.h:40
Int_t
int Int_t
Definition: External.C:63
Float_t
float Float_t
Definition: External.C:68
AliConversionPhotonBase.h
isMC
Bool_t isMC
Definition: ReadDvsMultiplicity.C:60
AliConvEventCuts.h
TH2
Definition: External.C:220
AliAnalysisTaskMaterialHistos::SetConversionCutList
void SetConversionCutList(Int_t nCuts, TList *CutArray)
Definition: AliAnalysisTaskMaterialHistos.h:42
AliV0ReaderV1.h
Option_t
const char Option_t
Definition: External.C:48
Bool_t
bool Bool_t
Definition: External.C:53