]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGGA/EMCALTasks/AliAnalysisTaskEMCALPhoton.h
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
update from Constantin
[u/mrichter/AliRoot.git] / PWGGA / EMCALTasks / AliAnalysisTaskEMCALPhoton.h
CommitLineData
f2f73202 1#ifndef AliAnalysisTaskEMCALPhoton_h\r
2#define AliAnalysisTaskEMCALPhoton_h\r
33d8da67 3\r
cd231d42 4// $Id$\r
5\r
33d8da67 6class TH1;\r
7class TH2;\r
8class TObjArray;\r
9class AliESDEvent;\r
10class AliMCEvent;\r
11class AliStack;\r
12class AliESDtrack;\r
13class AliESDtrackCuts;\r
14class AliESDCaloCells;\r
15class AliEMCALGeometry;\r
16class AliVCluster;\r
17class AliAnalysisTaskEMCALClusterizeFast;\r
18class TParticle;\r
b2d49404 19class AliPhotonHeaderObj;\r
20class AliPhotonConvObj;\r
21class AliPhotonClusterObj;\r
22class AliPhotonCellObj;\r
33d8da67 23\r
24#include "AliAnalysisTaskSE.h"\r
25\r
26class AliAnalysisTaskEMCALPhoton : public AliAnalysisTaskSE {\r
27 public:\r
f2f73202 28 AliAnalysisTaskEMCALPhoton();\r
33d8da67 29 AliAnalysisTaskEMCALPhoton(const char *name);\r
30 virtual ~AliAnalysisTaskEMCALPhoton() {}\r
31\r
32 void UserCreateOutputObjects();\r
33 void UserExec(Option_t *option);\r
34 void Terminate(Option_t *);\r
35\r
36 void SetTrackCuts(AliESDtrackCuts *c) { fTrCuts = c; }\r
37 void SetPrimTrackCuts(AliESDtrackCuts *c) { fPrTrCuts = c; }\r
38 void SetTimeResTOF(Float_t tr = 130.) { fTimeResTOF = tr; }\r
39 void SetMipResponseTPC(Float_t mr = 47.9) { fMipResponseTPC = mr; }\r
40 void SetGeoName(const char *n) { fGeoName = n; }\r
41 void SetPeriod(const char *n) { fPeriod = n; }\r
42 void SetTrainMode(Bool_t t) { fIsTrain = t; }\r
43 void SetGridMode(Bool_t g) { fIsGrid = g; }\r
44 void SetClusThreshold(Double_t et) { fClusThresh = et; }\r
45 void SetClusterizer(AliAnalysisTaskEMCALClusterizeFast *c) { fClusterizer = c; }\r
46 void SetMcMode(Bool_t mc) { fIsMC = mc; }\r
47 void FindConversions();\r
48 void FillMyCells();\r
49 void FillMyClusters();\r
50 void FillMyAltClusters();\r
3b37c011 51 void FillIsoTracks();\r
33d8da67 52 void FillMcPart(TParticle *mcP, Int_t ipart, Int_t itrack);\r
53 void GetMcParts();\r
54 Double_t GetTrackIsolation(Double_t cEta, Double_t cPhi, Double_t radius=0.2, Double_t pt=0.) const;\r
55 Double_t GetPhiBandEt(Double_t cEta, Double_t cPhi, Double_t radius=0.2, Double_t pt=0.) const;\r
56 // Double_t GetPhiBandEt(Double_t cEta, Double_t cPhi, Double_t radius=0.2, Double_t pt=0.) const;\r
57 Double_t GetCrossEnergy(const AliVCluster *cluster, Short_t &idmax);\r
58 Double_t GetMaxCellEnergy(const AliVCluster *cluster, Short_t &id) const; \r
59 \r
60 protected:\r
61 AliESDtrackCuts *fTrCuts; // track cuts\r
62 AliESDtrackCuts *fPrTrCuts; // primary track cuts\r
63 TObjArray *fSelTracks; //!pointer to selected inclusive tracks\r
64 TObjArray *fSelPrimTracks; //!pointer to selected primary tracks\r
b2d49404 65 TClonesArray *fPhotConvArray; //!array of AliPhotonConvObj\r
66 TClonesArray *fMyClusts; //!array of AliPhotonClusterObj\r
67 TClonesArray *fMyAltClusts; //!array of AliPhotonClusterObj from the alternative clusterizer\r
68 TClonesArray *fMyCells; //!array of AliPhotonCellObj\r
69 TClonesArray *fMyTracks; //!array of AliPhotonTrackObj\r
70 TClonesArray *fMyMcParts; //!array of AliPhotonMcPartObj\r
71 AliPhotonHeaderObj *fHeader; //!\r
33d8da67 72 TRefArray *fCaloClusters; //!pointer to EMCal clusters\r
73 TClonesArray *fCaloClustersNew; //!pointer to EMCal clusters v2\r
74 AliESDCaloCells *fEMCalCells; //!pointer to EMCal cells\r
75 AliEMCALGeometry *fGeom; // geometry utils\r
76 Float_t fTimeResTOF; //TOF time resolution for track PID\r
77 Float_t fMipResponseTPC; //TPC mip response for track pid\r
78 TString fGeoName; // geometry name (def = EMCAL_FIRSTYEARV1)\r
79 TString fPeriod; // string to the LHC period\r
80 Bool_t fIsTrain; //variable to set train mode\r
81 Bool_t fIsMC; //variable to switch mcparts branch on/off\r
82 Bool_t fIsGrid; //variable to set grid mode\r
83 Double_t fClusThresh; //!energy threshold for cluster be saved\r
84 AliAnalysisTaskEMCALClusterizeFast *fClusterizer; //!pointer for alternative clusterizer\r
85 TString fCaloClustersName; //alternative clusterizer name\r
86 \r
87 \r
88 private:\r
89 AliESDEvent *fESD; //! ESD object\r
90 AliMCEvent *fMCEvent; //! MC event object\r
91 AliStack *fStack; //!MC particles stack object\r
92 \r
93 TList *fOutputList; //! Output list\r
94 TTree *fTree; //!output tree\r
95 \r
96 //conversion histograms\r
97 TH2F *fNV0sBefAndAftRerun; //!check the number of V0s before and after rerun\r
98 TH2F *fConversionVtxXY; //! X x Y of found conversion vertices\r
99 TH1F *fInvMassV0; //!invariant mass from v0->GetEffMass()\r
100 TH1F *fInvMassV0KF; //!invariant mass from the v0 tracks\r
101 TH1F *fInvMassV0SS; //!invariant mass from the tracks in the "dirty" finder\r
102 TH2F *fDedxPAll; //!dE/dx vs p of all selected tracks\r
103 \r
104\r
105 \r
106 AliAnalysisTaskEMCALPhoton(const AliAnalysisTaskEMCALPhoton&); // not implemented\r
107 AliAnalysisTaskEMCALPhoton& operator=(const AliAnalysisTaskEMCALPhoton&); // not implemented\r
108 \r
109 ClassDef(AliAnalysisTaskEMCALPhoton, 1); // example of analysis\r
110};\r
111\r
112#endif\r
113\r
b2d49404 114#ifndef AliPhotonObjs_h\r
115#define AliPhotonObjs_h\r
33d8da67 116\r
b2d49404 117class AliPhotonHeaderObj : public TObject\r
33d8da67 118{\r
b2d49404 119 public: AliPhotonHeaderObj() :\r
33d8da67 120 TObject(), fTrClassMask(0), fTrCluster(0), fV0Cent(0), fV0(0), fCl1Cent(0), \r
b8b781cb 121 fCl1(0), fTrCent(0), fTr(0), fNClus(0), fNCells(0), fTrackMult(0) {;}\r
33d8da67 122 public:\r
123 ULong64_t fTrClassMask; // trigger class mask\r
124 UChar_t fTrCluster; // trigger cluster mask\r
125 Double32_t fV0Cent; //[0,0,16] v0 cent\r
126 Double32_t fV0; //[0,0,16] v0 result used for cent \r
127 Double32_t fCl1Cent; //[0,0,16] cl1 cent\r
128 Double32_t fCl1; //[0,0,16] cl1 result used for cent \r
129 Double32_t fTrCent; //[0,0,16] tr cent\r
130 Double32_t fTr; //[0,0,16] tr result used for cent \r
131 Int_t fNClus;\r
132 Int_t fNCells;\r
b8b781cb 133 Int_t fTrackMult;\r
33d8da67 134\r
b8b781cb 135 ClassDef(AliPhotonHeaderObj,3)\r
33d8da67 136};\r
137\r
b2d49404 138class AliPhotonConvObj : public TObject\r
33d8da67 139{\r
b2d49404 140 public: AliPhotonConvObj() : \r
33d8da67 141 TObject(), fPt(0), fEta(0), fPhi(0), fVR(0), fVEta(0), fVPhi(0), fMass(0), fMcLabel(-1),\r
142 fNegPt(0), fNegEta(0), fNegPhi(0), fNegDedx(0), fNegMcLabel(-1),\r
143 fPosPt(0), fPosEta(0), fPosPhi(0), fPosDedx(0), fPosMcLabel(-1) {;}\r
144 public:\r
145 Double32_t fPt; //[0,0,16] pt\r
146 Double32_t fEta; //[0,0,16] eta\r
147 Double32_t fPhi; //[0,0,16] phi\r
148 Double32_t fVR; //[0,0,16] prod r (cylinder)\r
149 Double32_t fVEta; //[0,0,16] prod eta\r
150 Double32_t fVPhi; //[0,0,16] prod phi\r
151 Double32_t fMass; //[0,0,16] if correctly filled, should be <50 MeV\r
152 Short_t fMcLabel; //corresponding MC label\r
153\r
154 //negative daughter\r
155 Double32_t fNegPt; //[0,0,16] pt\r
156 Double32_t fNegEta; //[0,0,16] eta\r
157 Double32_t fNegPhi; //[0,0,16] phi\r
158 Double32_t fNegDedx; //[0,0,16] if correctly filled, should be <50 MeV\r
159 Short_t fNegMcLabel; //corresponding MC label\r
160\r
161 //positive daughter\r
162 Double32_t fPosPt; //[0,0,16] pt\r
163 Double32_t fPosEta; //[0,0,16] eta\r
164 Double32_t fPosPhi; //[0,0,16] phi\r
165 Double32_t fPosDedx; //[0,0,16] if correctly filled, should be <50 MeV\r
166 Short_t fPosMcLabel; //corresponding MC label\r
167\r
b2d49404 168 ClassDef(AliPhotonConvObj,1) // conversion class\r
33d8da67 169\r
170};\r
b2d49404 171class AliPhotonClusterObj : public TObject\r
33d8da67 172{\r
b2d49404 173 public: AliPhotonClusterObj() : \r
33d8da67 174 TObject(), fE(0), fEt(0), fR(0), fEta(0), fPhi(0), fN(0),fEmax(0),fTmax(0), fIdmax(0), fEcross(0),fDisp(-1), \r
175 fM20(-1), fM02(-1),fTrDEta(0), fTrDPhi(0), fTrEp(0), fTrDedx(0), fTrIso01(0), fTrIso02(0), fTrIso03(0), fTrIso04(0), \r
176 fTrPhiBand01(0), fTrPhiBand02(0), fTrPhiBand03(0), fTrPhiBand04(0), fCellsAbsId(""),fMcLabel(-1)\r
177 {;}\r
178 public:\r
179 Double32_t fE;\r
180 Double32_t fEt;\r
181 Double32_t fR;\r
182 Double32_t fEta;\r
183 Double32_t fPhi;\r
184 UShort_t fN;\r
185 Double_t fEmax;\r
186 Double_t fTmax;\r
187 Short_t fIdmax;\r
188 Double_t fEcross;\r
189 Double32_t fDisp;\r
190 Double32_t fM20;\r
191 Double32_t fM02;\r
192 Double32_t fTrDEta;\r
193 Double32_t fTrDPhi;\r
194 Double32_t fTrEp;\r
195 Double32_t fTrDedx;\r
196 Double32_t fTrIso01;\r
197 Double32_t fTrIso02;\r
198 Double32_t fTrIso03;\r
199 Double32_t fTrIso04;\r
200 Double32_t fTrPhiBand01;\r
201 Double32_t fTrPhiBand02;\r
202 Double32_t fTrPhiBand03;\r
203 Double32_t fTrPhiBand04;\r
204 TString fCellsAbsId; //cluster cells absid\r
205 Short_t fMcLabel;\r
206 \r
207 \r
208 \r
b2d49404 209 ClassDef(AliPhotonClusterObj,6) // cluster class\r
33d8da67 210\r
211};\r
212\r
b2d49404 213class AliPhotonCellObj : public TObject\r
33d8da67 214{\r
b2d49404 215 public: AliPhotonCellObj() : \r
33d8da67 216 TObject(), fAbsID(-1), fE(0), fEt(0), fEta(0), fPhi(0), fTime(0)\r
217 {;}\r
218 public:\r
219 Short_t fAbsID;\r
220 Double32_t fE;\r
221 Double32_t fEt;\r
222 Double32_t fEta;\r
223 Double32_t fPhi;\r
224 Double32_t fTime;\r
225 \r
226 \r
227 \r
b2d49404 228 ClassDef(AliPhotonCellObj,1) // cell class\r
33d8da67 229\r
230};\r
231\r
b2d49404 232class AliPhotonTrackObj : public TObject\r
33d8da67 233{\r
b2d49404 234 public: AliPhotonTrackObj() :\r
33d8da67 235 TObject(), fPt(0), fEta(0), fPhi(0), fDedx(0), fCharge(0), fMcLabel(-1) {;}\r
236 public:\r
237 Double32_t fPt;\r
238 Double32_t fEta;\r
239 Double32_t fPhi;\r
240 Double32_t fDedx;\r
241 Short_t fCharge;\r
242 Short_t fMcLabel;\r
243\r
b2d49404 244 ClassDef(AliPhotonTrackObj,3)\r
33d8da67 245};\r
246\r
b2d49404 247class AliPhotonMcPartObj : public TObject\r
33d8da67 248{\r
b2d49404 249 public: AliPhotonMcPartObj() :\r
33d8da67 250 TObject(), fLabel(-1), fPdg(0), fPt(0), fEta(0), fPhi(0), \r
251 fVR(0), fVEta(0), fVPhi(0), fMother(-1) {;}\r
252 public:\r
253 Short_t fLabel;\r
254 Short_t fPdg;\r
255 Double32_t fPt;\r
256 Double32_t fEta;\r
257 Double32_t fPhi;\r
258 Double32_t fVR;\r
259 Double32_t fVEta;\r
260 Double32_t fVPhi;\r
261 Short_t fMother;\r
262\r
b2d49404 263 ClassDef(AliPhotonMcPartObj,1)\r
33d8da67 264};\r
265\r
266#endif\r