]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWG4/PartCorrDep/AliAnaOmegaToPi0Gamma.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
coverity fix
[u/mrichter/AliRoot.git] / PWG4 / PartCorrDep / AliAnaOmegaToPi0Gamma.cxx
CommitLineData
0852acd0 1/**************************************************************************\r
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *\r
3 * *\r
4 * Author: The ALICE Off-line Project. *\r
5 * Contributors are mentioned in the code where appropriate. *\r
6 * *\r
7 * Permission to use, copy, modify and distribute this software and its *\r
8 * documentation strictly for non-commercial purposes is hereby granted *\r
9 * without fee, provided that the above copyright notice appears in all *\r
10 * copies and that both the copyright notice and this permission notice *\r
11 * appear in the supporting documentation. The authors make no claims *\r
12 * about the suitability of this software for any purpose. It is *\r
13 * provided "as is" without express or implied warranty. *\r
14 **************************************************************************/\r
15/* $Id: $ */\r
16//_________________________________________________________________________\r
17// class to extract omega(782)->pi0+gamma->3gamma\r
18//\r
19//-- Author: Renzhuo Wan (IOPP-Wuhan, China)\r
20//_________________________________________________________________________\r
21\r
22// --- ROOT system\r
23class TROOT;\r
24\r
25// --- AliRoot system\r
26//class AliVEvent;\r
27// --- ROOT system ---\r
28#include "TH2F.h"\r
29#include "TLorentzVector.h"\r
30#include "TParticle.h"\r
31#include "TCanvas.h"\r
32#include "TFile.h"\r
33//---- AliRoot system ----\r
34#include "AliAnaOmegaToPi0Gamma.h"\r
35#include "AliCaloTrackReader.h"\r
36#include "AliCaloPID.h"\r
37#include "AliStack.h"\r
38#include "AliVEvent.h"\r
39#include "AliAODEvent.h"\r
40#include "AliAODMCParticle.h"\r
41ClassImp(AliAnaOmegaToPi0Gamma)\r
42\r
43//______________________________________________________________________________\r
44AliAnaOmegaToPi0Gamma::AliAnaOmegaToPi0Gamma() : AliAnaPartCorrBaseClass(),\r
edff0964 45fInputAODPi0(0), fInputAODGammaName(""),\r
1e3689fa 46fEventsList(0x0),fNVtxZBin(0), fNCentBin(0), fNRpBin(0), fNBadChDistBin(0), fNpid(0),\r
edff0964 47fNmaxMixEv(0), fVtxZCut(0), fCent(0), fRp(0), \r
48fPi0Mass(0),fPi0MassWindow(0),fPi0OverOmegaPtCut(0),\r
49fGammaOverOmegaPtCut(0),\r
0852acd0 50fhEtalon(0),\r
d5b2a3f1 51fRealOmega0(0), fMixAOmega0(0),\r
52fMixBOmega0(0), fMixCOmega0(0),\r
0852acd0 53fRealOmega1(0), fMixAOmega1(0),\r
54fMixBOmega1(0), fMixCOmega1(0),\r
55fRealOmega2(0), fMixAOmega2(0),\r
56fMixBOmega2(0), fMixCOmega2(0),\r
57fhOmegaPriPt(0)\r
58{\r
59 //Default Ctor\r
edff0964 60 InitParameters();\r
0852acd0 61}\r
78219bac 62/*\r
0852acd0 63//______________________________________________________________________________\r
64AliAnaOmegaToPi0Gamma::AliAnaOmegaToPi0Gamma(const AliAnaOmegaToPi0Gamma & ex) : AliAnaPartCorrBaseClass(ex),\r
0852acd0 65fInputAODPi0(new TClonesArray (*ex.fInputAODPi0)),\r
66fInputAODGammaName(ex.fInputAODGammaName),\r
1e3689fa 67fEventsList(0x0), \r
0852acd0 68fNVtxZBin(ex.fNVtxZBin), fNCentBin(ex.fNCentBin), fNRpBin(ex.fNRpBin),\r
d5b2a3f1 69fNBadChDistBin(ex.fNBadChDistBin),fNpid(ex.fNpid),\r
70fNmaxMixEv(ex.fNmaxMixEv),\r
f6862ceb 71fVtxZCut(ex.fVtxZCut), fCent(ex.fCent), fRp(ex.fRp), \r
d5b2a3f1 72fPi0Mass(ex.fPi0Mass),\r
73fPi0MassWindow(ex.fPi0MassWindow),\r
74fPi0OverOmegaPtCut(ex.fPi0OverOmegaPtCut),\r
75fGammaOverOmegaPtCut(ex.fGammaOverOmegaPtCut),\r
0852acd0 76fhEtalon(ex.fhEtalon),\r
d5b2a3f1 77fRealOmega0(ex.fRealOmega0), fMixAOmega0(ex.fMixAOmega0),\r
78fMixBOmega0(ex.fMixBOmega0), fMixCOmega0(ex.fMixCOmega0),\r
0852acd0 79fRealOmega1(ex.fRealOmega1), fMixAOmega1(ex.fMixAOmega1),\r
80fMixBOmega1(ex.fMixBOmega1), fMixCOmega1(ex.fMixCOmega1),\r
81fRealOmega2(ex.fRealOmega2), fMixAOmega2(ex.fMixAOmega2),\r
82fMixBOmega2(ex.fMixBOmega2), fMixCOmega2(ex.fMixCOmega2),\r
83fhOmegaPriPt(ex.fhOmegaPriPt)\r
84{\r
85 // cpy ctor\r
86 //Do not need it\r
87}\r
78219bac 88*/\r
1e3689fa 89/*\r
0852acd0 90//______________________________________________________________________________\r
91AliAnaOmegaToPi0Gamma & AliAnaOmegaToPi0Gamma::operator = (const AliAnaOmegaToPi0Gamma & ex)\r
92{\r
93 // assignment operator\r
c5b989a8 94 \r
0852acd0 95 if(this == &ex)return *this;\r
96 ((AliAnaPartCorrBaseClass *)this)->operator=(ex);\r
0852acd0 97 fInputAODPi0 = new TClonesArray(*ex.fInputAODPi0);\r
0852acd0 98 fInputAODGammaName = ex.fInputAODGammaName;\r
99 fEventsList = ex.fEventsList;\r
100\r
d5b2a3f1 101 fNVtxZBin=ex.fNVtxZBin;\r
0852acd0 102 fNCentBin=ex.fNCentBin;\r
103 fNRpBin=ex.fNRpBin;\r
104 fNBadChDistBin=ex.fNBadChDistBin;\r
105 fNpid=ex.fNpid;\r
0852acd0 106 fNmaxMixEv =ex.fNmaxMixEv;\r
d5b2a3f1 107\r
108 fVtxZCut=ex.fVtxZCut;\r
109 fCent=ex.fCent;\r
110 fRp=ex.fRp;\r
d5b2a3f1 111\r
0852acd0 112 fPi0Mass=ex.fPi0Mass;\r
113 fPi0MassWindow=ex.fPi0MassWindow;\r
4007d0a2 114 fPi0OverOmegaPtCut=ex.fPi0OverOmegaPtCut;\r
115 fGammaOverOmegaPtCut=ex.fGammaOverOmegaPtCut;\r
0852acd0 116\r
0852acd0 117 fhEtalon=ex.fhEtalon;\r
d5b2a3f1 118 fRealOmega0=ex.fRealOmega0;\r
119 fMixAOmega0=ex.fMixAOmega0;\r
120 fMixBOmega0=ex.fMixBOmega0;\r
121 fMixCOmega0=ex.fMixCOmega0;\r
0852acd0 122 fRealOmega1=ex.fRealOmega1;\r
123 fMixAOmega1=ex.fMixAOmega1;\r
124 fMixBOmega1=ex.fMixBOmega1;\r
125 fMixCOmega1=ex.fMixCOmega1;\r
126 fRealOmega2=ex.fRealOmega2;\r
127 fMixAOmega2=ex.fMixAOmega2;\r
128 fMixBOmega2=ex.fMixBOmega2;\r
129 fMixCOmega2=ex.fMixCOmega2;\r
130 fhOmegaPriPt=ex.fhOmegaPriPt;\r
131 return *this;\r
132 \r
133}\r
1e3689fa 134*/\r
edff0964 135\r
0852acd0 136//______________________________________________________________________________\r
137AliAnaOmegaToPi0Gamma::~AliAnaOmegaToPi0Gamma() {\r
138\r
139 //dtor\r
7787a778 140// Done by the maker \r
141// if(fInputAODPi0){\r
142// fInputAODPi0->Clear();\r
143// delete fInputAODPi0;\r
144// } \r
0852acd0 145\r
146 if(fEventsList){\r
147 for(Int_t i=0;i<fNVtxZBin;i++){\r
148 for(Int_t j=0;j<fNCentBin;j++){\r
149 for(Int_t k=0;k<fNRpBin;k++){\r
150 fEventsList[i*fNCentBin*fNRpBin+j*fNRpBin+k]->Clear();\r
151 delete fEventsList[i*fNCentBin*fNRpBin+j*fNRpBin+k];\r
152 }\r
153 }\r
154 }\r
155 }\r
156 delete [] fEventsList;\r
157 fEventsList=0;\r
d5b2a3f1 158\r
159 delete [] fVtxZCut;\r
160 delete [] fCent;\r
161 delete [] fRp;\r
d5b2a3f1 162\r
0852acd0 163}\r
164\r
165//______________________________________________________________________________\r
166void AliAnaOmegaToPi0Gamma::InitParameters()\r
167{\r
168//Init parameters when first called the analysis\r
169//Set default parameters\r
edff0964 170 fInputAODGammaName="PhotonsDetector"; \r
171 fNVtxZBin=1; \r
172 fNCentBin=1; \r
173 fNRpBin=1; \r
174 fNBadChDistBin=3; \r
175 fNpid=9; \r
176 fNmaxMixEv=8; \r
177 \r
178 fPi0Mass=0.1348; \r
179 fPi0MassWindow=0.015; \r
180 fPi0OverOmegaPtCut=0.8; \r
181 fGammaOverOmegaPtCut=0.2; \r
182}\r
183\r
184\r
185//______________________________________________________________________________\r
186TList * AliAnaOmegaToPi0Gamma::GetCreateOutputObjects()\r
187{\r
188 //\r
0852acd0 189 fVtxZCut = new Double_t [fNVtxZBin];\r
4007d0a2 190 for(Int_t i=0;i<fNVtxZBin;i++) fVtxZCut[i]=10*(i+1);\r
f8006433 191 \r
0852acd0 192 fCent=new Double_t[fNCentBin];\r
c5b989a8 193 for(int i = 0;i<fNCentBin;i++)fCent[i]=0;\r
f8006433 194 \r
0852acd0 195 fRp=new Double_t[fNRpBin];\r
c5b989a8 196 for(int i = 0;i<fNRpBin;i++)fRp[i]=0;\r
edff0964 197 //\r
4007d0a2 198 Int_t nptbins = GetHistoPtBins();\r
199 Float_t ptmax = GetHistoPtMax();\r
200 Float_t ptmin = GetHistoPtMin();\r
f8006433 201 \r
4007d0a2 202 Int_t nmassbins = GetHistoMassBins();\r
203 Float_t massmin = GetHistoMassMin();\r
204 Float_t massmax = GetHistoMassMax();\r
f8006433 205 \r
4007d0a2 206 fhEtalon = new TH2F("hEtalon","Histo with binning parameters", nptbins,ptmin,ptmax,nmassbins,massmin,massmax) ;\r
207 fhEtalon->SetXTitle("P_{T} (GeV)") ;\r
208 fhEtalon->SetYTitle("m_{inv} (GeV)") ;\r
f8006433 209 \r
210 // store them in fOutputContainer\r
211 fEventsList = new TList*[fNVtxZBin*fNCentBin*fNRpBin];\r
212 for(Int_t i=0;i<fNVtxZBin;i++){\r
213 for(Int_t j=0;j<fNCentBin;j++){\r
214 for(Int_t k=0;k<fNRpBin;k++){\r
af7b3903 215 fEventsList[i*fNCentBin*fNRpBin+j*fNRpBin+k] = new TList();\r
216 fEventsList[i*fNCentBin*fNRpBin+j*fNRpBin+k]->SetOwner(kFALSE); \r
f8006433 217 }\r
218 }\r
219 }\r
0852acd0 220 \r
f8006433 221 TList * outputContainer = new TList() ; \r
222 outputContainer->SetName(GetName());\r
223 const Int_t buffersize = 255;\r
224 char key[buffersize] ;\r
225 char title[buffersize] ;\r
226 const char * detector= fInputAODGammaName.Data();\r
227 Int_t ndim=fNVtxZBin*fNCentBin*fNRpBin*fNBadChDistBin*fNpid;\r
228 \r
229 fRealOmega0 =new TH2F*[ndim];\r
230 fMixAOmega0 =new TH2F*[ndim];\r
231 fMixBOmega0 =new TH2F*[ndim];\r
232 fMixCOmega0 =new TH2F*[ndim];\r
233 \r
234 fRealOmega1 =new TH2F*[ndim];\r
235 fMixAOmega1 =new TH2F*[ndim];\r
236 fMixBOmega1 =new TH2F*[ndim];\r
237 fMixCOmega1 =new TH2F*[ndim];\r
238 \r
239 fRealOmega2 =new TH2F*[ndim];\r
240 fMixAOmega2 =new TH2F*[ndim];\r
241 fMixBOmega2 =new TH2F*[ndim];\r
242 fMixCOmega2 =new TH2F*[ndim];\r
243 \r
244 for(Int_t i=0;i<fNVtxZBin;i++){\r
245 for(Int_t j=0;j<fNCentBin;j++){\r
246 for(Int_t k=0;k<fNRpBin;k++){ //at event level\r
247 Int_t idim=i*fNCentBin*fNRpBin+j*fNRpBin+k;\r
248 for(Int_t ipid=0;ipid<fNpid;ipid++){\r
249 for(Int_t idist=0;idist<fNBadChDistBin;idist++){ //at particle level\r
250 \r
251 Int_t index=idim*fNpid*fNBadChDistBin+ipid*fNBadChDistBin+idist;\r
252 \r
253 snprintf(key,buffersize,"RealToPi0Gamma_Vz%dC%dRp%dPid%dDist%d",i,j,k,ipid,idist);\r
254 snprintf(title,buffersize, "%s Real Pi0GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%d",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,idist);\r
255 fhEtalon->Clone(key);\r
256 fRealOmega0[index]=(TH2F*)fhEtalon->Clone(key) ;\r
257 fRealOmega0[index]->SetName(key) ;\r
258 fRealOmega0[index]->SetTitle(title);\r
259 outputContainer->Add(fRealOmega0[index]);\r
260 \r
261 snprintf(key,buffersize,"MixAToPi0Gamma_Vz%dC%dRp%dPid%dDist%d",i,j,k,ipid,idist);\r
262 snprintf(title,buffersize, "%s MixA Pi0GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%d",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,idist);\r
263 fhEtalon->Clone(key);\r
264 fMixAOmega0[index]=(TH2F*)fhEtalon->Clone(key) ;\r
265 fMixAOmega0[index]->SetName(key) ;\r
266 fMixAOmega0[index]->SetTitle(title);\r
267 outputContainer->Add(fMixAOmega0[index]);\r
268 \r
269 snprintf(key,buffersize,"MixBToPi0Gamma_Vz%dC%dRp%dPid%dDist%d",i,j,k,ipid,idist);\r
270 snprintf(title,buffersize, "%s MixB Pi0GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%d",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,idist);\r
271 fhEtalon->Clone(key);\r
272 fMixBOmega0[index]=(TH2F*)fhEtalon->Clone(key) ;\r
273 fMixBOmega0[index]->SetName(key) ;\r
274 fMixBOmega0[index]->SetTitle(title);\r
275 outputContainer->Add(fMixBOmega0[index]);\r
276 \r
277 snprintf(key,buffersize,"MixCToPi0Gamma_Vz%dC%dRp%dPid%dDist%d",i,j,k,ipid,idist);\r
278 snprintf(title,buffersize, "%s MixC Pi0GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%d",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,idist);\r
279 fhEtalon->Clone(key);\r
280 fMixCOmega0[index]=(TH2F*)fhEtalon->Clone(key) ;\r
281 fMixCOmega0[index]->SetName(key) ;\r
282 fMixCOmega0[index]->SetTitle(title);\r
283 outputContainer->Add(fMixCOmega0[index]);\r
284 \r
285 snprintf(key,buffersize,"RealToPi0Gamma1_Vz%dC%dRp%dPid%dDist%d",i,j,k,ipid,idist);\r
286 snprintf(title,buffersize, "%s Real Pi0(A<0.7)GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%d",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,idist);\r
287 fhEtalon->Clone(key);\r
288 fRealOmega1[index]=(TH2F*)fhEtalon->Clone(key) ;\r
289 fRealOmega1[index]->SetName(key) ;\r
290 fRealOmega1[index]->SetTitle(title);\r
291 outputContainer->Add(fRealOmega1[index]);\r
292 \r
293 snprintf(key,buffersize,"MixAToPi0Gamma1_Vz%dC%dRp%dPid%dDist%d",i,j,k,ipid,idist);\r
294 snprintf(title,buffersize, "%s MixA Pi0(A<0.7)GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%d",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,idist);\r
295 fhEtalon->Clone(key);\r
296 fMixAOmega1[index]=(TH2F*)fhEtalon->Clone(key) ;\r
297 fMixAOmega1[index]->SetName(key) ;\r
298 fMixAOmega1[index]->SetTitle(title);\r
299 outputContainer->Add(fMixAOmega1[index]);\r
300 \r
301 snprintf(key,buffersize,"MixBToPi0Gamma1_Vz%dC%dRp%dPid%dDist%d",i,j,k,ipid,idist);\r
302 snprintf(title,buffersize, "%s MixB Pi0(A<0.7)GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%d",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,idist);\r
303 fhEtalon->Clone(key);\r
304 fMixBOmega1[index]=(TH2F*)fhEtalon->Clone(key) ;\r
305 fMixBOmega1[index]->SetName(key) ;\r
306 fMixBOmega1[index]->SetTitle(title);\r
307 outputContainer->Add(fMixBOmega1[index]);\r
308 \r
309 snprintf(key,buffersize,"MixCToPi0Gamma1_Vz%dC%dRp%dPid%dDist%d",i,j,k,ipid,idist);\r
310 snprintf(title,buffersize, "%s MixC Pi0(A<0.7)GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%d",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,idist);\r
311 fhEtalon->Clone(key);\r
312 fMixCOmega1[index]=(TH2F*)fhEtalon->Clone(key) ;\r
313 fMixCOmega1[index]->SetName(key) ;\r
314 fMixCOmega1[index]->SetTitle(title);\r
315 outputContainer->Add(fMixCOmega1[index]);\r
316 \r
317 snprintf(key,buffersize,"RealToPi0Gamma2_Vz%dC%dRp%dPid%dDist%d",i,j,k,ipid,idist);\r
318 snprintf(title,buffersize, "%s Real Pi0(A<0.8)GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%d",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,idist);\r
319 fhEtalon->Clone(key);\r
320 fRealOmega2[index]=(TH2F*)fhEtalon->Clone(key) ;\r
321 fRealOmega2[index]->SetName(key) ;\r
322 fRealOmega2[index]->SetTitle(title);\r
323 outputContainer->Add(fRealOmega2[index]);\r
324 \r
325 snprintf(key,buffersize,"MixAToPi0Gamma2_Vz%dC%dRp%dPid%dDist%d",i,j,k,ipid,idist);\r
326 snprintf(title,buffersize, "%s MixA Pi0(A<0.8)GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%d",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,idist);\r
327 fhEtalon->Clone(key);\r
328 fMixAOmega2[index]=(TH2F*)fhEtalon->Clone(key) ;\r
329 fMixAOmega2[index]->SetName(key) ;\r
330 fMixAOmega2[index]->SetTitle(title);\r
331 outputContainer->Add(fMixAOmega2[index]);\r
332 \r
333 snprintf(key,buffersize,"MixBToPi0Gamma2_Vz%dC%dRp%dPid%dDist%d",i,j,k,ipid,idist);\r
334 snprintf(title,buffersize, "%s MixB Pi0(A<0.8)GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%d",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,idist);\r
335 fhEtalon->Clone(key);\r
336 fMixBOmega2[index]=(TH2F*)fhEtalon->Clone(key) ;\r
337 fMixBOmega2[index]->SetName(key) ;\r
338 fMixBOmega2[index]->SetTitle(title);\r
339 outputContainer->Add(fMixBOmega2[index]);\r
340 \r
341 snprintf(key,buffersize,"MixCToPi0Gamma2_Vz%dC%dRp%dPid%dDist%d",i,j,k,ipid,idist);\r
342 snprintf(title,buffersize, "%s MixC Pi0(A<0.8)GammaIVM vz_%2.1f_ct_%2.1f_Rp_%2.1f_pid_%d_dist_%d",detector,fVtxZCut[i],fCent[j],fRp[k],ipid,idist);\r
343 fhEtalon->Clone(key);\r
344 fMixCOmega2[index]=(TH2F*)fhEtalon->Clone(key) ;\r
345 fMixCOmega2[index]->SetName(key) ;\r
346 fMixCOmega2[index]->SetTitle(title);\r
347 outputContainer->Add(fMixCOmega2[index]);\r
348 }\r
349 }\r
350 }\r
351 } \r
352 }\r
353 \r
354 if(IsDataMC()){\r
355 snprintf(key,buffersize, "%sOmegaPri",detector);\r
356 snprintf(title,buffersize,"primary #omega in %s",detector);\r
4007d0a2 357 fhOmegaPriPt=new TH1F(key, title,nptbins,ptmin,ptmax);\r
0852acd0 358 fhOmegaPriPt->GetXaxis()->SetTitle("P_{T}");\r
359 fhOmegaPriPt->GetYaxis()->SetTitle("dN/P_{T}");\r
360 outputContainer->Add(fhOmegaPriPt);\r
f8006433 361 }\r
362 \r
363 delete fhEtalon;\r
364 return outputContainer;\r
0852acd0 365}\r
366\r
367//______________________________________________________________________________\r
368void AliAnaOmegaToPi0Gamma::Print(const Option_t * /*opt*/) const\r
369{\r
370 //Print some relevant parameters set in the analysis\r
371 printf("**** Print %s %s ****\n", GetName(), GetTitle() ) ;\r
372 AliAnaPartCorrBaseClass::Print(" ");\r
373 printf("Omega->pi0+gamma->3gamma\n");\r
374 printf("Cuts at event level: \n");\r
375 printf("Bins of vertex Z: %d \n", fNVtxZBin);\r
376 printf("Bins of centrality: %d \n",fNCentBin);\r
377 printf("Bins of Reaction plane: %d\n",fNRpBin);\r
378 printf("Cuts at AOD particle level:\n");\r
379 printf("Number of PID: %d \n", fNpid);\r
380 printf("Number of DistToBadChannel cuts: %d\n", fNBadChDistBin);\r
381 printf("number of events buffer to be mixed: %d\n",fNmaxMixEv);\r
382} \r
383\r
384//______________________________________________________________________________\r
385void AliAnaOmegaToPi0Gamma::MakeAnalysisFillHistograms() \r
386{\r
f8006433 387 //fill the MC AOD if needed first\r
388 //-----------\r
389 //need to be further implemented\r
390 AliStack * stack = 0x0;\r
391 // TParticle * primary = 0x0;\r
392 TClonesArray * mcparticles0 = 0x0;\r
393 //TClonesArray * mcparticles1 = 0x0;\r
394 AliAODMCParticle * aodprimary = 0x0;\r
395 Int_t pdg=0;\r
396 Double_t pt=0;\r
397 Double_t eta=0;\r
398 \r
399 if(IsDataMC()){\r
400 if(GetReader()->ReadStack()){\r
401 stack = GetMCStack() ;\r
402 if(!stack){\r
0852acd0 403 printf("AliAnaAcceptance::MakeAnalysisFillHistograms() - There is no stack!\n");\r
f8006433 404 }\r
405 else{\r
0852acd0 406 for(Int_t i=0 ; i<stack->GetNtrack(); i++){\r
f8006433 407 TParticle * prim = stack->Particle(i) ;\r
408 pdg = prim->GetPdgCode() ;\r
409 eta=prim->Eta();\r
410 pt=prim->Pt();\r
411 if(TMath::Abs(eta)<0.5) {\r
412 if(pdg==223) fhOmegaPriPt->Fill(pt);\r
413 }\r
0852acd0 414 }\r
415 }\r
f8006433 416 }\r
417 else if(GetReader()->ReadAODMCParticles()){\r
418 //Get the list of MC particles\r
419 mcparticles0 = GetReader()->GetAODMCParticles(0);\r
420 if(!mcparticles0 ) {\r
421 if(GetDebug() > 0) printf("AliAnaAcceptance::MakeAnalysisFillHistograms() - Standard MCParticles not available!\n");\r
422 }\r
423 // if(GetReader()->GetSecondInputAODTree()){\r
424 // mcparticles1 = GetReader()->GetAODMCParticles(1);\r
425 // if(!mcparticles1 && GetDebug() > 0) {\r
426 // printf("AliAnaAcceptance::MakeAnalysisFillHistograms() - Second input MCParticles not available!\n");\r
427 // }\r
428 // }\r
429 else{\r
430 for(Int_t i=0;i<mcparticles0->GetEntries();i++){\r
431 aodprimary =(AliAODMCParticle*)mcparticles0->At(i);\r
432 pdg = aodprimary->GetPdgCode() ;\r
433 eta=aodprimary->Eta();\r
434 pt=aodprimary->Pt();\r
435 if(TMath::Abs(eta)<0.5) {\r
436 if(pdg==223) fhOmegaPriPt->Fill(pt);\r
437 }\r
438 \r
439 }\r
440 }//mcparticles0 exists\r
441 }//AOD MC Particles\r
442 }// is data and MC\r
0852acd0 443 \r
f8006433 444 \r
445 //process event from AOD brach \r
446 //extract pi0, eta and omega analysis\r
447 Int_t iRun=(GetReader()->GetInputEvent())->GetRunNumber() ;\r
448 if(IsBadRun(iRun)) return ; \r
449 \r
450 //vertex z\r
451 Double_t vert[]={0,0,0} ;\r
452 GetVertex(vert);\r
453 Int_t curEventBin =0;\r
454 \r
455 Int_t ivtxzbin=(Int_t)TMath::Abs(vert[2])/10;\r
456 if(ivtxzbin>=fNVtxZBin)return;\r
457 \r
458 //centrality\r
459 Int_t icentbin=0;\r
460 \r
461 //reaction plane\r
462 Int_t irpbin=0;\r
463 \r
464 if(ivtxzbin==-1) return; \r
465 curEventBin = ivtxzbin*fNCentBin*fNRpBin + icentbin*fNRpBin + irpbin;\r
edff0964 466 TClonesArray *aodGamma = (TClonesArray*) GetAODBranch(fInputAODGammaName); //photon array\r
f8006433 467 // TClonesArray *aodGamma = (TClonesArray *) GetReader()->GetOutputEvent()->FindListObject(fInputAODGammaName); //photon array\r
468 Int_t nphotons =0;\r
469 if(aodGamma) nphotons= aodGamma->GetEntries(); \r
470 else return;\r
471 \r
472 fInputAODPi0 = (TClonesArray*)GetInputAODBranch(); //pi0 array\r
473 Int_t npi0s = 0;\r
474 if(fInputAODPi0) npi0s= fInputAODPi0 ->GetEntries();\r
475 else return;\r
476 \r
477 if(nphotons<3 || npi0s<1)return; //for pi0, eta and omega->pi0+gamma->3gamma reconstruction\r
478 \r
479 //reconstruction of omega(782)->pi0+gamma->3gamma\r
480 //loop for pi0 and photon\r
481 if(GetDebug() > 0) printf("omega->pi0+gamma->3gamma invariant mass analysis ! This event have %d photons and %d pi0 \n", nphotons, npi0s);\r
482 for(Int_t i=0;i<npi0s;i++){\r
483 AliAODPWG4Particle * pi0 = (AliAODPWG4Particle*) (fInputAODPi0->At(i)) ; //pi0\r
484 TLorentzVector vpi0(pi0->Px(),pi0->Py(),pi0->Pz(),pi0->E());\r
485 Int_t lab1=pi0->GetCaloLabel(0); // photon1 from pi0 decay\r
486 Int_t lab2=pi0->GetCaloLabel(1); // photon2 from pi0 decay\r
487 //for omega->pi0+gamma, it needs at least three photons per event\r
488 //Get the two decay photons from pi0\r
489 AliAODPWG4Particle * photon1 =0;\r
490 AliAODPWG4Particle * photon2 =0;\r
491 for(Int_t d1=0;d1<nphotons;d1++){\r
492 for(Int_t d2=0;d2<nphotons;d2++){\r
493 AliAODPWG4Particle * dp1 = (AliAODPWG4Particle*) (aodGamma->At(d1));\r
494 AliAODPWG4Particle * dp2 = (AliAODPWG4Particle*) (aodGamma->At(d2));\r
495 Int_t dlab1=dp1->GetCaloLabel(0);\r
496 Int_t dlab2=dp2->GetCaloLabel(0);\r
497 if(dlab1==lab1 && dlab2==lab2){\r
498 photon1=dp1;\r
499 photon2=dp2;\r
500 }\r
501 else continue;\r
502 }\r
0852acd0 503 }\r
0852acd0 504 //caculate the asy and dist of the two photon from pi0 decay\r
505 TLorentzVector dph1(photon1->Px(),photon1->Py(),photon1->Pz(),photon1->E());\r
506 TLorentzVector dph2(photon2->Px(),photon2->Py(),photon2->Pz(),photon2->E());\r
f8006433 507 \r
0852acd0 508 Double_t pi0asy= TMath::Abs(dph1.E()-dph2.E())/(dph1.E()+dph2.E());\r
f8006433 509 // Double_t phi1=dph1.Phi();\r
510 // Double_t phi2=dph2.Phi();\r
511 // Double_t eta1=dph1.Eta();\r
512 // Double_t eta2=dph2.Eta();\r
513 // Double_t pi0dist=TMath::Sqrt((phi1-phi2)*(phi1-phi2)+(eta1-eta2)*(eta1-eta2));\r
d5b2a3f1 514 \r
f8006433 515 if(pi0->GetPdg()==111 && nphotons>2 && npi0s\r
516 && TMath::Abs(vpi0.M()-fPi0Mass)<fPi0MassWindow) { //pi0 candidates\r
517 \r
518 //avoid the double counting\r
519 Int_t * dc1= new Int_t[nphotons];\r
520 Int_t * dc2= new Int_t[nphotons];\r
521 Int_t index1=0;\r
522 Int_t index2=0;\r
523 for(Int_t k=0;k<i;k++){\r
524 AliAODPWG4Particle * p3=(AliAODPWG4Particle*)(fInputAODPi0->At(k));\r
525 Int_t lab4=p3->GetCaloLabel(0);\r
526 Int_t lab5=p3->GetCaloLabel(1);\r
527 if(lab1==lab4){ dc1[index1]=lab5; index1++; }\r
528 if(lab2==lab5){ dc2[index2]=lab4; index2++; }\r
529 }\r
530 \r
531 \r
532 //loop the pi0 with third gamma\r
533 for(Int_t j=0;j<nphotons;j++){\r
534 AliAODPWG4Particle *photon3 = (AliAODPWG4Particle*) (aodGamma->At(j));\r
535 TLorentzVector dph3(photon3->Px(),photon3->Py(),photon3->Pz(),photon3->E());\r
536 Int_t lab3=photon3->GetCaloLabel(0);\r
537 Double_t pi0gammapt=(vpi0+dph3).Pt();\r
538 Double_t pi0gammamass=(vpi0+dph3).M();\r
539 Double_t pi0OverOmegaPtRatio =vpi0.Pt()/pi0gammapt; \r
540 Double_t gammaOverOmegaPtRatio= dph3.Pt()/pi0gammapt;\r
541 \r
542 //pi0, gamma pt cut \r
543 if(pi0OverOmegaPtRatio>fPi0OverOmegaPtCut || \r
544 gammaOverOmegaPtRatio<fGammaOverOmegaPtCut) continue;\r
545 \r
546 for(Int_t l=0;l<index1;l++) if(lab3==dc1[l]) lab3=-1;\r
547 for(Int_t l=0;l<index2;l++) if(lab3==dc2[l]) lab3=-1;\r
548 \r
549 if(lab3>0 && lab3!=lab1 && lab3!=lab2){\r
0852acd0 550 for(Int_t ipid=0;ipid<fNpid;ipid++){\r
f8006433 551 for(Int_t idist=0;idist<fNBadChDistBin;idist++){\r
552 Int_t index=curEventBin*fNpid*fNBadChDistBin+ipid*fNBadChDistBin+idist;\r
553 if(photon1->IsPIDOK(ipid,AliCaloPID::kPhoton) &&\r
554 photon2->IsPIDOK(ipid,AliCaloPID::kPhoton) && \r
555 photon3->IsPIDOK(ipid,AliCaloPID::kPhoton) &&\r
556 photon1->DistToBad()>=idist &&\r
557 photon2->DistToBad()>=idist &&\r
558 photon3->DistToBad()>=idist ){\r
559 //fill the histograms\r
560 if(GetDebug() > 2) printf("Real: index %d pt %2.3f mass %2.3f \n", index, pi0gammapt, pi0gammamass);\r
561 fRealOmega0[index]->Fill(pi0gammapt,pi0gammamass);\r
562 if(pi0asy<0.7) fRealOmega1[index]->Fill(pi0gammapt,pi0gammamass);\r
563 if(pi0asy<0.8) fRealOmega2[index]->Fill(pi0gammapt,pi0gammamass);\r
564 }\r
565 }\r
0852acd0 566 }\r
f8006433 567 }\r
568 } \r
569 delete []dc1;\r
570 delete []dc2;\r
571 if(GetDebug() > 0) printf("MixA: (r1_event1+r2_event1)+r3_event2 \n");\r
572 //-------------------------\r
573 //background analysis\r
574 //three background\r
575 // --A (r1_event1+r2_event1)+r3_event2\r
576 Int_t nMixed = fEventsList[curEventBin]->GetSize();\r
577 for(Int_t im=0;im<nMixed;im++){\r
578 TClonesArray* ev2= (TClonesArray*) (fEventsList[curEventBin]->At(im));\r
579 for(Int_t mix1=0;mix1<ev2->GetEntries();mix1++){\r
580 AliAODPWG4Particle *mix1ph = (AliAODPWG4Particle*) (ev2->At(mix1)); \r
581 TLorentzVector vmixph(mix1ph->Px(),mix1ph->Py(),mix1ph->Pz(),mix1ph->E());\r
582 Double_t pi0gammapt=(vpi0+vmixph).Pt();\r
583 Double_t pi0gammamass=(vpi0+vmixph).M();\r
584 Double_t pi0OverOmegaPtRatio =vpi0.Pt()/pi0gammapt;\r
585 Double_t gammaOverOmegaPtRatio= vmixph.Pt()/pi0gammapt;\r
586 \r
587 //pi0, gamma pt cut \r
588 if(pi0OverOmegaPtRatio>fPi0OverOmegaPtCut || \r
589 gammaOverOmegaPtRatio<fGammaOverOmegaPtCut) continue;\r
590 \r
0852acd0 591 for(Int_t ipid=0;ipid<fNpid;ipid++){\r
f8006433 592 for(Int_t idist=0;idist<fNBadChDistBin;idist++){\r
593 Int_t index=curEventBin*fNpid*fNBadChDistBin+ipid*fNBadChDistBin+idist;\r
594 if(photon1->IsPIDOK(ipid,AliCaloPID::kPhoton)&&\r
595 photon2->IsPIDOK(ipid,AliCaloPID::kPhoton)&&\r
596 mix1ph->IsPIDOK(ipid,AliCaloPID::kPhoton) &&\r
597 photon1->DistToBad()>=idist &&\r
598 photon2->DistToBad()>=idist &&\r
599 mix1ph->DistToBad()>=idist ){\r
600 if(GetDebug() > 2) printf("MixA: index %d pt %2.3f mass %2.3f \n",index, pi0gammapt, pi0gammamass);\r
601 //fill the histograms\r
602 fMixAOmega0[index]->Fill(pi0gammapt,pi0gammamass);\r
603 if(pi0asy<0.7)fMixAOmega1[index]->Fill(pi0gammapt,pi0gammamass);\r
604 if(pi0asy<0.8)fMixAOmega2[index]->Fill(pi0gammapt,pi0gammamass);\r
605 //printf("mix A %d %2.2f \n", index, pi0gammamass);\r
606 \r
607 }\r
608 }\r
609 }\r
610 }\r
611 }\r
612 }\r
613 }\r
614 \r
615 //\r
616 // --B (r1_event1+r2_event2)+r3_event2\r
617 //\r
618 if(GetDebug() >0)printf("MixB: (r1_event1+r2_event2)+r3_event2 \n");\r
619 for(Int_t i=0;i<nphotons;i++){\r
620 AliAODPWG4Particle *ph1 = (AliAODPWG4Particle*) (aodGamma->At(i)); \r
621 TLorentzVector vph1(ph1->Px(),ph1->Py(),ph1->Pz(),ph1->E());\r
622 \r
623 Int_t nMixed = fEventsList[curEventBin]->GetSize();\r
624 for(Int_t ie=0;ie<nMixed;ie++){\r
625 TClonesArray* ev2= (TClonesArray*) (fEventsList[curEventBin]->At(ie));\r
626 for(Int_t mix1=0;mix1<ev2->GetEntries();mix1++){\r
627 AliAODPWG4Particle *ph2 = (AliAODPWG4Particle*) (ev2->At(mix1));\r
628 TLorentzVector vph2(ph2->Px(),ph2->Py(),ph2->Pz(),ph2->E());\r
629 Double_t pi0asy = TMath::Abs(vph1.E()-vph2.E())/(vph1.E()+vph2.E()); \r
630 Double_t pi0mass=(vph1+vph2).M();\r
631 \r
632 if(TMath::Abs(pi0mass-fPi0Mass)<fPi0MassWindow){//for pi0 selection\r
633 for(Int_t mix2=(mix1+1);mix2<ev2->GetEntries();mix2++){\r
634 AliAODPWG4Particle *ph3 = (AliAODPWG4Particle*) (ev2->At(mix2));\r
635 TLorentzVector vph3(ph3->Px(),ph3->Py(),ph3->Pz(),ph3->E());\r
636 \r
637 Double_t pi0gammapt=(vph1+vph2+vph3).Pt();\r
638 Double_t pi0gammamass=(vph1+vph2+vph3).M(); \r
639 Double_t pi0OverOmegaPtRatio =(vph1+vph2).Pt()/pi0gammapt;\r
640 Double_t gammaOverOmegaPtRatio= vph3.Pt()/pi0gammapt;\r
641 //pi0, gamma pt cut \r
642 if(pi0OverOmegaPtRatio>fPi0OverOmegaPtCut ||\r
643 gammaOverOmegaPtRatio<fGammaOverOmegaPtCut) continue;\r
644 \r
645 for(Int_t ipid=0;ipid<fNpid;ipid++){\r
646 for(Int_t idist=0;idist<fNBadChDistBin;idist++){\r
647 Int_t index=curEventBin*fNpid*fNBadChDistBin+ipid*fNBadChDistBin+idist;\r
648 if(ph1->IsPIDOK(ipid,AliCaloPID::kPhoton) &&\r
0852acd0 649 ph2->IsPIDOK(ipid,AliCaloPID::kPhoton) &&\r
f8006433 650 ph3->IsPIDOK(ipid,AliCaloPID::kPhoton) &&\r
651 ph1->DistToBad()>=idist &&\r
652 ph2->DistToBad()>=idist &&\r
653 ph3->DistToBad()>=idist ){\r
654 if(GetDebug() > 2) printf("MixB: index %d pt %2.3f mass %2.3f \n", index, pi0gammapt, pi0gammamass);\r
655 //fill histograms\r
656 fMixBOmega0[index]->Fill(pi0gammapt,pi0gammamass);\r
657 if(pi0asy<0.7) fMixBOmega1[index]->Fill(pi0gammapt,pi0gammamass);\r
658 if(pi0asy<0.8) fMixBOmega2[index]->Fill(pi0gammapt,pi0gammamass);\r
659 //printf("mix B %d %2.2f \n", index, pi0gammamass);\r
660 }\r
661 } \r
662 }\r
663 }\r
664 \r
665 //\r
0852acd0 666 // --C (r1_event1+r2_event2)+r3_event3\r
f8006433 667 //\r
668 if(GetDebug() >0)printf("MixC: (r1_event1+r2_event2)+r3_event3\n");\r
669 for(Int_t je=(ie+1);je<nMixed;je++){\r
670 TClonesArray* ev3= (TClonesArray*) (fEventsList[curEventBin]->At(je));\r
671 for(Int_t mix3=0;mix3<ev3->GetEntries();mix3++){\r
672 AliAODPWG4Particle *ph3 = (AliAODPWG4Particle*) (ev3->At(mix3));\r
673 TLorentzVector vph3(ph3->Px(),ph3->Py(),ph3->Pz(),ph3->E());\r
674 \r
675 Double_t pi0gammapt=(vph1+vph2+vph3).Pt();\r
676 Double_t pi0gammamass=(vph1+vph2+vph3).M();\r
677 Double_t pi0OverOmegaPtRatio =(vph1+vph2).Pt()/pi0gammapt;\r
678 Double_t gammaOverOmegaPtRatio= vph3.Pt()/pi0gammapt;\r
679 //pi0, gamma pt cut \r
680 if(pi0OverOmegaPtRatio>fPi0OverOmegaPtCut ||\r
681 gammaOverOmegaPtRatio<fGammaOverOmegaPtCut) continue;\r
682 \r
683 for(Int_t ipid=0;ipid<fNpid;ipid++){\r
684 for(Int_t idist=0;idist<fNBadChDistBin;idist++){\r
685 Int_t index=curEventBin*fNpid*fNBadChDistBin+ipid*fNBadChDistBin+idist;\r
686 if(ph1->IsPIDOK(ipid,AliCaloPID::kPhoton) &&\r
687 ph2->IsPIDOK(ipid,AliCaloPID::kPhoton) &&\r
688 ph3->IsPIDOK(ipid,AliCaloPID::kPhoton) &&\r
689 ph1->DistToBad()>=idist &&\r
690 ph2->DistToBad()>=idist &&\r
691 ph3->DistToBad()>=idist ){\r
692 if(GetDebug() > 2) printf("MixC: index %d pt %2.3f mass %2.3f \n", index, pi0gammapt, pi0gammamass);\r
693 //fill histograms\r
694 fMixCOmega0[index]->Fill(pi0gammapt,pi0gammamass);\r
695 if(pi0asy<0.7) fMixCOmega1[index]->Fill(pi0gammapt,pi0gammamass);\r
696 if(pi0asy<0.8) fMixCOmega2[index]->Fill(pi0gammapt,pi0gammamass);\r
697 //printf("mix C %d %2.2f \n", index, pi0gammamass);\r
698 }\r
699 }\r
700 }\r
701 }\r
702 }\r
703 } //for pi0 selecton \r
704 }\r
705 }\r
706 }\r
707 \r
708 \r
709 //event buffer \r
710 TClonesArray *currentEvent = new TClonesArray(*aodGamma);\r
711 if(currentEvent->GetEntriesFast()>0){\r
712 fEventsList[curEventBin]->AddFirst(currentEvent) ;\r
713 currentEvent=0 ; \r
714 if(fEventsList[curEventBin]->GetSize()>=fNmaxMixEv) {\r
715 TClonesArray * tmp = (TClonesArray*) (fEventsList[curEventBin]->Last()) ;\r
716 fEventsList[curEventBin]->RemoveLast() ;\r
717 delete tmp ;\r
718 }\r
719 }\r
720 else{ \r
721 delete currentEvent ;\r
722 currentEvent=0 ;\r
723 }\r
724 \r
0852acd0 725}\r
726\r
727//______________________________________________________________________________\r
728void AliAnaOmegaToPi0Gamma::ReadHistograms(TList * outputList)\r
729{\r
730 //read the histograms \r
731 //for the finalization of the terminate analysis\r
732\r
d5b2a3f1 733 Int_t index = outputList->IndexOf(outputList->FindObject(GetAddedHistogramsStringToName()+"RealToPi0Gamma_Vz0C0Rp0Pid0Dist0"));\r
0852acd0 734\r
735 Int_t ndim=fNVtxZBin*fNCentBin*fNRpBin*fNBadChDistBin*fNpid;\r
736\r
d5b2a3f1 737 if(!fRealOmega0) fRealOmega0 =new TH2F*[ndim];\r
738 if(!fMixAOmega0) fMixAOmega0 =new TH2F*[ndim];\r
739 if(!fMixBOmega0) fMixBOmega0 =new TH2F*[ndim];\r
740 if(!fMixCOmega0) fMixCOmega0 =new TH2F*[ndim];\r
0852acd0 741\r
742 if(!fRealOmega1) fRealOmega1 =new TH2F*[ndim];\r
743 if(!fMixAOmega1) fMixAOmega1 =new TH2F*[ndim];\r
744 if(!fMixBOmega1) fMixBOmega1 =new TH2F*[ndim];\r
745 if(!fMixCOmega1) fMixCOmega1 =new TH2F*[ndim];\r
746\r
747 if(!fRealOmega2) fRealOmega2 =new TH2F*[ndim];\r
748 if(!fMixAOmega2) fMixAOmega2 =new TH2F*[ndim];\r
749 if(!fMixBOmega2) fMixBOmega2 =new TH2F*[ndim];\r
750 if(!fMixCOmega2) fMixCOmega2 =new TH2F*[ndim];\r
751\r
752 for(Int_t i=0;i<fNVtxZBin;i++){\r
753 for(Int_t j=0;j<fNCentBin;j++){\r
754 for(Int_t k=0;k<fNRpBin;k++){ //at event level\r
755 Int_t idim=i*fNCentBin*fNRpBin+j*fNRpBin+k;\r
756 for(Int_t ipid=0;ipid<fNpid;ipid++){ \r
757 for(Int_t idist=0;idist<fNBadChDistBin;idist++){ //at particle\r
758 Int_t ind=idim*fNpid*fNBadChDistBin+ipid*fNBadChDistBin+idist;\r
d5b2a3f1 759 fRealOmega0[ind]= (TH2F*) outputList->At(index++);\r
760 fMixAOmega0[ind]= (TH2F*) outputList->At(index++);\r
761 fMixBOmega0[ind]= (TH2F*) outputList->At(index++);\r
762 fMixCOmega0[ind]= (TH2F*) outputList->At(index++);\r
0852acd0 763\r
764 fRealOmega1[ind]= (TH2F*) outputList->At(index++);\r
765 fMixAOmega1[ind]= (TH2F*) outputList->At(index++);\r
766 fMixBOmega1[ind]= (TH2F*) outputList->At(index++);\r
767 fMixCOmega1[ind]= (TH2F*) outputList->At(index++);\r
768\r
769 fRealOmega2[ind]= (TH2F*) outputList->At(index++);\r
770 fMixAOmega2[ind]= (TH2F*) outputList->At(index++);\r
771 fMixBOmega2[ind]= (TH2F*) outputList->At(index++);\r
772 fMixCOmega2[ind]= (TH2F*) outputList->At(index++);\r
773 \r
774 \r
775 }\r
776 }\r
777 }\r
778 }\r
779 }\r
780 \r
781 if(IsDataMC()){\r
782 fhOmegaPriPt = (TH1F*) outputList->At(index++);\r
783 }\r
784\r
785}\r
786\r
787//______________________________________________________________________________\r
788void AliAnaOmegaToPi0Gamma::Terminate(TList * outputList) \r
789{\r
790// //Do some calculations and plots from the final histograms.\r
791 if(GetDebug() >= 0) printf("AliAnaOmegaToPi0Gamma::Terminate() \n");\r
792 ReadHistograms(outputList);\r
f8006433 793 const Int_t buffersize = 255;\r
794 char cvs1[buffersize]; \r
164a1d84 795 snprintf(cvs1, buffersize, "Neutral_%s_IVM",fInputAODGammaName.Data());\r
0852acd0 796\r
797 TCanvas * cvsIVM = new TCanvas(cvs1, cvs1, 400, 10, 600, 700) ;\r
798 cvsIVM->Divide(2, 2);\r
799\r
800 cvsIVM->cd(1);\r
f8006433 801 char dec[buffersize];\r
802 snprintf(dec,buffersize,"h2Real_%s",fInputAODGammaName.Data());\r
d5b2a3f1 803 TH2F * h2Real= (TH2F*)fRealOmega0[0]->Clone(dec);\r
0852acd0 804 h2Real->GetXaxis()->SetRangeUser(4,6);\r
805 TH1F * hRealOmega = (TH1F*) h2Real->ProjectionY();\r
806 hRealOmega->SetTitle("RealPi0Gamma 4<pt<6");\r
807 hRealOmega->SetLineColor(2);\r
808 hRealOmega->Draw();\r
809\r
810 cvsIVM->cd(2);\r
f8006433 811 snprintf(dec,buffersize,"hMixA_%s",fInputAODGammaName.Data());\r
d5b2a3f1 812 TH2F *h2MixA= (TH2F*)fMixAOmega0[0]->Clone(dec);\r
0852acd0 813 h2MixA->GetXaxis()->SetRangeUser(4,6);\r
814 TH1F * hMixAOmega = (TH1F*) h2MixA->ProjectionY();\r
815 hMixAOmega->SetTitle("MixA 4<pt<6");\r
816 hMixAOmega->SetLineColor(2);\r
817 hMixAOmega->Draw();\r
818\r
819 cvsIVM->cd(3);\r
f8006433 820 snprintf(dec,buffersize,"hMixB_%s",fInputAODGammaName.Data());\r
d5b2a3f1 821 TH2F * h2MixB= (TH2F*)fMixBOmega0[0]->Clone(dec);\r
0852acd0 822 h2MixB->GetXaxis()->SetRangeUser(4,6);\r
823 TH1F * hMixBOmega = (TH1F*) h2MixB->ProjectionY();\r
824 hMixBOmega->SetTitle("MixB 4<pt<6");\r
825 hMixBOmega->SetLineColor(2);\r
826 hMixBOmega->Draw();\r
827\r
828 cvsIVM->cd(4);\r
f8006433 829 snprintf(dec,buffersize,"hMixC_%s",fInputAODGammaName.Data());\r
d5b2a3f1 830 TH2F *h2MixC= (TH2F*)fMixCOmega0[0]->Clone(dec);\r
0852acd0 831 h2MixC->GetXaxis()->SetRangeUser(4,6);\r
832 TH1F * hMixCOmega = (TH1F*) h2MixC->ProjectionY();\r
833 hMixCOmega->SetTitle("MixC 4<pt<6");\r
834 hMixCOmega->SetLineColor(2);\r
835 hMixCOmega->Draw();\r
836\r
f8006433 837 char eps[buffersize];\r
838 snprintf(eps,buffersize,"CVS_%s_IVM.eps",fInputAODGammaName.Data());\r
0852acd0 839 cvsIVM->Print(eps);\r
840 cvsIVM->Modified();\r
841 \r
842}\r