]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWG4/GammaConv/AliGammaConversionHistograms.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
More fixes about macros
[u/mrichter/AliRoot.git] / PWG4 / GammaConv / AliGammaConversionHistograms.cxx
CommitLineData
d7d7e825 1/**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: Ana Marin, Kathrin Koch, Kenneth Aamodt *
5 * Version 1.1 *
6 * *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
15
16////////////////////////////////////////////////
17//---------------------------------------------
18// Class used to do analysis on conversion pairs
19//---------------------------------------------
20////////////////////////////////////////////////
21
22#include "AliGammaConversionHistograms.h"
23#include "TMath.h"
24#include "TObjString.h"
25#include "TMap.h"
26#include "TList.h"
27#include "TH1F.h"
28#include "TH2F.h"
29
30
31using namespace std;
32
33ClassImp(AliGammaConversionHistograms)
34
35
36AliGammaConversionHistograms::AliGammaConversionHistograms() :
37 fHistogramMap(new TMap()),
38 fNPhiIndex(0),
39 fNRIndex(0),
9640a3d1 40 fNZIndex(0),
5a34881d 41 fRBinLimits(0),
42 fZBinLimits(0),
d7d7e825 43 fMinRadius(0.),
44 fMaxRadius(0.),
45 fDeltaR(0.),
46 fMinPhi(0.),
47 fMaxPhi(0.),
48 fDeltaPhi(0.),
49 fMappingContainer(NULL),
50 fBackgroundContainer(NULL),
51 fDebugContainer(NULL),
52 fResolutionContainer(NULL),
53 fMatchContainer(NULL),
54 fESDContainer(NULL),
55 fMCContainer(NULL),
56 fTableContainer(NULL),
57 fOtherContainer(NULL)
58{
59 // see header file for documenation
60}
61
62
63AliGammaConversionHistograms::AliGammaConversionHistograms(const AliGammaConversionHistograms & original) :
64 fHistogramMap(original.fHistogramMap),
65 fNPhiIndex(original.fNPhiIndex),
66 fNRIndex(original.fNRIndex),
9640a3d1 67 fNZIndex(original.fNZIndex),
5a34881d 68 fRBinLimits(original.fRBinLimits),
69 fZBinLimits(original.fZBinLimits),
d7d7e825 70 fMinRadius(original.fMinRadius),
71 fMaxRadius(original.fMaxRadius),
72 fDeltaR(original.fDeltaR),
73 fMinPhi(original.fMinPhi),
74 fMaxPhi(original.fMaxPhi),
75 fDeltaPhi(original.fDeltaPhi),
76 fMappingContainer(original.fMappingContainer),
77 fBackgroundContainer(original.fBackgroundContainer),
78 fDebugContainer(original.fDebugContainer),
79 fResolutionContainer(original.fResolutionContainer),
80 fMatchContainer(original.fMatchContainer),
81 fESDContainer(original.fESDContainer),
82 fMCContainer(original.fMCContainer),
83 fTableContainer(original.fTableContainer),
84 fOtherContainer(original.fOtherContainer)
85{
86 //see header file for documentation
87}
88
89
90AliGammaConversionHistograms & AliGammaConversionHistograms::operator = (const AliGammaConversionHistograms & /*original*/)
91{
92 // assignment operator
93 return *this;
94}
95
96
97AliGammaConversionHistograms::~AliGammaConversionHistograms() {
98 //destructor
99
100
101}
102
103void AliGammaConversionHistograms::AddHistogram(TString histogramName, TString histogramTitle, Int_t nXBins, Double_t firstX,Double_t lastX,TString xAxisTitle, TString yAxisTitle){
104 // see header file for documentation
105 TH1F *tmp = new TH1F(histogramName, histogramTitle,nXBins,firstX,lastX);
106 tmp->GetXaxis()->SetTitle(xAxisTitle);
107 tmp->GetYaxis()->SetTitle(yAxisTitle);
108 TObjString* tobjstring = new TObjString(histogramName.Data());
109 fHistogramMap->Add((TObject*)tobjstring,(TObject*)tmp);
110}
111
112void AliGammaConversionHistograms::AddHistogram(TString histogramName, TString histogramTitle, Int_t nXBins, Double_t firstX, Double_t lastX, Int_t nYBins, Double_t firstY, Double_t lastY, TString xAxisTitle, TString yAxisTitle){
113 // see header file for documentation
114 TH2F *tmp = new TH2F(histogramName, histogramTitle,nXBins,firstX,lastX,nYBins,firstY,lastY);
115 tmp->GetXaxis()->SetTitle(xAxisTitle);
116 tmp->GetYaxis()->SetTitle(yAxisTitle);
117 TObjString *tobjstring = new TObjString(histogramName.Data());
118 fHistogramMap->Add((TObject*)tobjstring,(TObject*)tmp);
119}
120
121void AliGammaConversionHistograms::AddTable(TString tableName,TString tableTitle,Int_t nXBins,const char * axesLabel[]){
4a6157dc 122 //see header file for documentation
d7d7e825 123
4a6157dc 124 TH1F *tmp = new TH1F(tableName,tableTitle,nXBins,0,nXBins);
125 for(Int_t xbin=1; xbin<=nXBins; xbin++){
126 tmp->GetXaxis()->SetBinLabel(xbin,axesLabel[xbin-1]);
127 }
128 tmp->SetStats(0);
d7d7e825 129
4a6157dc 130 TObjString *tobjstring = new TObjString(tableName.Data());
131 fHistogramMap->Add((TObject*)tobjstring,(TObject*)tmp);
d7d7e825 132}
133
134void AliGammaConversionHistograms::FillTable(TString tableName,Double_t xValue) const {
4a6157dc 135 //see header file for documentation
136 TH1 *tmp = (TH1*)fHistogramMap->GetValue(tableName.Data());
137 if(tmp){
138 tmp->Fill(xValue);
139 }
d7d7e825 140}
141
142void AliGammaConversionHistograms::FillHistogram(TString histogramName, Double_t xValue) const{
143 //see header file for documentation
144 TH1 *tmp = (TH1*)fHistogramMap->GetValue(histogramName.Data());
145 if(tmp){
146 tmp->Fill(xValue);
147 }
148}
149
150void AliGammaConversionHistograms::FillHistogram(TString histogramName, Double_t xValue, Double_t yValue) const{
151 //see header file for documentation
152 TH1 *tmp = (TH1*)fHistogramMap->GetValue(histogramName.Data());
153 if(tmp){
154 tmp->Fill(xValue, yValue);
155 }
156}
157
158void AliGammaConversionHistograms::GetOutputContainer(TList *fOutputContainer){
159 //checking if the container is alrerady created
160
161 if(fOutputContainer == NULL){
162 cout<<"WARNING: GetOutputContainer: output container object is NULL"<<endl;
163 return;
164 }
165
166 if(fHistogramMap != NULL){
167 TIter iter(fHistogramMap);
168 TObjString *histogramName;
169 while ((histogramName = (TObjString*) iter.Next())) {
170 TString histogramString = histogramName->GetString();
171 if(histogramString.Contains("Mapping")){// means it should be put in the mapping folder
172 if(fMappingContainer == NULL){
173 fMappingContainer = new TList();
174 fMappingContainer->SetName("Mapping histograms");
175 }
176 if(fMappingContainer != NULL){
177 fMappingContainer->Add((TH1*)fHistogramMap->GetValue(histogramString.Data()));
178 }
179 }
180 else if(histogramString.Contains("Background")){// means it should be put in the background folder
181 if(fBackgroundContainer == NULL){
182 fBackgroundContainer = new TList();
183 fBackgroundContainer->SetName("Background histograms");
184 }
185 if(fBackgroundContainer != NULL){
186 fBackgroundContainer->Add((TH1*)fHistogramMap->GetValue(histogramString.Data()));
187 }
188 }
189 else if(histogramString.Contains("Debug")){// means it should be put in the debug folder
190 if(fDebugContainer == NULL){
191 fDebugContainer = new TList();
192 fDebugContainer->SetName("Debug histograms");
193 }
194 if(fDebugContainer != NULL){
195 fDebugContainer->Add((TH1*)fHistogramMap->GetValue(histogramString.Data()));
196 }
197 }
198 else if(histogramString.Contains("Resolution")){// means it should be put in the resolution folder
199 if(fResolutionContainer == NULL){
200 fResolutionContainer = new TList();
201 fResolutionContainer->SetName("Resolution histograms");
202 }
203 if(fResolutionContainer != NULL){
204 fResolutionContainer->Add((TH1*)fHistogramMap->GetValue(histogramString.Data()));
205 }
206 }
207 else if(histogramString.Contains("TrueConv")){// means it should be put in the true conv folder
208 if(fMatchContainer == NULL){
209 fMatchContainer = new TList();
210 fMatchContainer->SetName("True conversion histograms");
211 }
212 if(fMatchContainer != NULL){
213 fMatchContainer->Add((TH1*)fHistogramMap->GetValue(histogramString.Data()));
214 }
215 }
216 else if(histogramString.Contains("ESD")){// means it should be put in the ESD folder
217 if(fESDContainer == NULL){
218 fESDContainer = new TList();
219 fESDContainer->SetName("ESD histograms");
220 }
221 if(fESDContainer != NULL){
222 fESDContainer->Add((TH1*)fHistogramMap->GetValue(histogramString.Data()));
223 }
224 }
225 else if(histogramString.Contains("MC")){// means it should be put in the MC folder
226 if(fMCContainer == NULL){
227 fMCContainer = new TList();
228 fMCContainer->SetName("MC histograms");
229 }
230 if(fMCContainer != NULL){
231 fMCContainer->Add((TH1*)fHistogramMap->GetValue(histogramString.Data()));
232 }
233 }
234 else if(histogramString.Contains("Table")){// means it should be put in the Table Folder
235 if(fTableContainer == NULL){
236 fTableContainer = new TList();
237 fTableContainer->SetName("Tables");
238 }
239 if(fTableContainer != NULL){
240 fTableContainer->Add((TH1*)fHistogramMap->GetValue(histogramString.Data()));
241 }
242 }
243 else{
244 if(fOtherContainer == NULL){
245 fOtherContainer = new TList();
246 fOtherContainer->SetName("Other histograms");
247 }
248 if(fOtherContainer != NULL){
249 fOtherContainer->Add((TH1*)fHistogramMap->GetValue(histogramString.Data()));
250 }
251 }
252 histogramName = NULL;
253 } // end while
254 if(fMappingContainer != NULL){
255 fOutputContainer->Add(fMappingContainer);
256 }
257 if(fBackgroundContainer != NULL){
258 fOutputContainer->Add(fBackgroundContainer);
259 }
260 if(fDebugContainer != NULL){
261 fOutputContainer->Add(fDebugContainer);
262 }
263 if(fResolutionContainer != NULL){
264 fOutputContainer->Add(fResolutionContainer);
265 }
266 if(fMatchContainer != NULL){
267 fOutputContainer->Add(fMatchContainer);
268 }
269 if(fESDContainer != NULL){
270 fOutputContainer->Add(fESDContainer);
271 }
272 if(fMCContainer != NULL){
273 fOutputContainer->Add(fMCContainer);
274 }
275 if(fTableContainer != NULL){
276 fOutputContainer->Add(fTableContainer);
277 }
278 if(fOtherContainer != NULL){
a68437fb 279 fOutputContainer->Add(fOtherContainer);
d7d7e825 280 }
281 }
282}
283
284Int_t AliGammaConversionHistograms::GetRBin(Double_t radius) const{
285 // see header file for documentation
286 Int_t iResult=0;
9640a3d1 287// if(fDeltaR>0){
288// iResult = (Int_t)((radius - fMinRadius)/fDeltaR);
289// }
290 for(Int_t i=0;i<fNRIndex;i++){
291 // cout<<"Test-limits::"<< fRBinLimits[i]<<endl;
292 if( radius>=fRBinLimits[i] && radius<fRBinLimits[i+1] ){
293 iResult=i;
294 }
d7d7e825 295 }
296 return iResult;
297}
298
9640a3d1 299Int_t AliGammaConversionHistograms::GetZBin(Double_t zPos) const{
300 // see header file for documentation
301 Int_t iResult=0;
302
303 for(Int_t i=0;i<fNZIndex;i++){
304 // cout<<"Test-limits::"<< fZBinLimits[i]<<endl;
305 if( TMath::Abs(zPos)>=fZBinLimits[i] && TMath::Abs(zPos)<fZBinLimits[i+1] ){
306 iResult=i;
307 }
308 }
309 return iResult;
310}
311
312
313
d7d7e825 314Int_t AliGammaConversionHistograms::GetPhiBin(Double_t phi) const{
315 // see header file for documentation
316 Int_t iResult=0;
317 if(fDeltaPhi>0){
318 if(phi>TMath::Pi()){
319 phi-=2*TMath::Pi();
320 }
321 iResult = (Int_t)((phi - fMinPhi)/fDeltaPhi);
322 }
323 return iResult;
324}
325
326
327
328void AliGammaConversionHistograms::InitializeMappingValues(Int_t nPhiIndex, Int_t nRIndex, Int_t nBinsR, Double_t minRadius, Double_t maxRadius,Int_t nBinsPhi, Double_t minPhi, Double_t maxPhi){
329 // Initializing the valuse for the mapping
330
331 fNPhiIndex = nPhiIndex;
332 fNRIndex = nRIndex;
9640a3d1 333 fNZIndex = 6;
334
335 fRBinLimits= new Double_t[8];
336 fRBinLimits[0]=0.;
337 fRBinLimits[1]=12.;
338 fRBinLimits[2]=22.;
339 fRBinLimits[3]=35.;
340 fRBinLimits[4]=55.;
341 fRBinLimits[5]=72.;
342 fRBinLimits[6]=90.;
343 fRBinLimits[7]=500.;
344
345 fZBinLimits= new Double_t[7];
346 fZBinLimits[0]=0.;
347 fZBinLimits[1]=15.;
348 fZBinLimits[2]=30.;
349 fZBinLimits[3]=50.;
350 fZBinLimits[4]=100.;
351 fZBinLimits[5]=200.;
352 fZBinLimits[6]=500.;
353
354
d7d7e825 355 fMinRadius = minRadius;
356 fMaxRadius = maxRadius;
357 if(nBinsR>0 && nRIndex!=0){
358 fDeltaR = (fMaxRadius - fMinRadius)/nRIndex;
359 }
360 fMinPhi = minPhi;
361 fMaxPhi = maxPhi;
362 if(nBinsPhi>0 && nPhiIndex!=0){
363 fDeltaPhi = (fMaxPhi-fMinPhi)/nPhiIndex;
364 }
365}
366
367
368//mapping
369void AliGammaConversionHistograms::AddMappingHistograms(Int_t nPhiIndex, Int_t nRIndex,Int_t nXBins, Double_t firstX, Double_t lastX, Int_t nYBins, Double_t firstY, Double_t lastY, TString xAxisTitle, TString yAxisTitle){
370 // see header file for documentation
371
9640a3d1 372 for(Int_t phi =0; phi<fNPhiIndex;phi++){
d7d7e825 373
374 for(Int_t r =0; r<fNRIndex;r++){
375
376 // setting axis to "" changes below
e158cbc3 377 xAxisTitle="z [cm]";
378 yAxisTitle="#eta";
379
d7d7e825 380 //Creating the axis titles
e158cbc3 381 //if(xAxisTitle.Length() == 0){
382 //xAxisTitle.Form("Phi %02d",phi);
383 // }
d7d7e825 384
e158cbc3 385 //if(yAxisTitle.Length() == 0){
386 //yAxisTitle.Form("R %02d",phi);
387 // }
d7d7e825 388
389 //MC
390 TString nameMC="";
e158cbc3 391 nameMC.Form("MC_Conversion_Mapping_Phi%02d_R%02d",phi,r);
d7d7e825 392 TString titleMC="";
393 titleMC.Form("Electron-Positron MC Mapping-Phi%02d-R%02d",phi,r);
394
395 AddHistogram(nameMC, titleMC, nXBins, firstX, lastX, nYBins, firstY, lastY, xAxisTitle, yAxisTitle);
396
397 //ESD
398 TString nameESD="";
e158cbc3 399 nameESD.Form("ESD_Conversion_Mapping_Phi%02d_R%02d",phi,r);
d7d7e825 400 TString titleESD="";
401 titleESD.Form("Electron-Positron ESD Mapping-Phi%02d-R%02d",phi,r);
402
403 AddHistogram(nameESD, titleESD, nXBins, firstX, lastX, nYBins, firstY, lastY, xAxisTitle, yAxisTitle);
404 }
405 }
406
407
408 for(Int_t phi =0; phi<=nPhiIndex;phi++){
409
410 // setting axis to "" changes below
e158cbc3 411 xAxisTitle="z [cm]";
412 yAxisTitle="#eta";
413 //Creating the axis titles
414 //if(xAxisTitle.Length() == 0){
415 // xAxisTitle.Form("Phi %02d",phi);
416 //}
417 //if(yAxisTitle.Length() == 0){
418 // yAxisTitle = "Counts";
419 //}
d7d7e825 420
421 //MC
422 TString nameMC="";
e158cbc3 423 nameMC.Form("MC_Conversion_Mapping_Phi%02d",phi);
d7d7e825 424 TString titleMC="";
425 titleMC.Form("Electron-Positron MC Mapping-Phi%02d",phi);
426
9640a3d1 427 AddHistogram(nameMC, titleMC, nXBins, firstX, lastX, nYBins, firstY, lastY, xAxisTitle, yAxisTitle);
d7d7e825 428
429 //MC
430 TString nameESD="";
e158cbc3 431 nameESD.Form("ESD_Conversion_Mapping_Phi%02d",phi);
d7d7e825 432 TString titleESD="";
433 titleESD.Form("Electron-Positron ESD Mapping-Phi%02d",phi);
434
9640a3d1 435 AddHistogram(nameESD, titleESD, nXBins, firstX, lastX,nYBins, firstY, lastY, xAxisTitle, yAxisTitle);
d7d7e825 436 }
437
438
9640a3d1 439 for(Int_t r =0; r<nRIndex;r++){
d7d7e825 440
441 // setting axis to "" changes below
e158cbc3 442 xAxisTitle="#phi";
443 yAxisTitle="counts";
d7d7e825 444 //Creating the axis titles
e158cbc3 445 //if(xAxisTitle.Length() == 0){
446 // xAxisTitle.Form("R %02d",r);
447 //}
448 //if(yAxisTitle.Length() == 0){
449 // yAxisTitle = "Counts";
450 //}
d7d7e825 451
452 //MC
453 TString nameMC="";
e158cbc3 454 nameMC.Form("MC_Conversion_Mapping_R%02d",r);
d7d7e825 455 TString titleMC="";
456 titleMC.Form("Electron-Positron MC Mapping-R%02d",r);
457
9640a3d1 458 AddHistogram(nameMC, titleMC, nXBins, firstX, lastX, nYBins, firstY, lastY, xAxisTitle, yAxisTitle);
d7d7e825 459
460 //ESD
461 TString nameESD="";
e158cbc3 462 nameESD.Form("ESD_Conversion_Mapping_R%02d",r);
d7d7e825 463 TString titleESD="";
464 titleESD.Form("Electron-Positron ESD Mapping-R%02d",r);
465
9640a3d1 466 AddHistogram(nameESD, titleESD, nXBins, firstX, lastX,nYBins, firstY, lastY, xAxisTitle, yAxisTitle);
d7d7e825 467
468 //Mapping Phi in R
469 TString nameMCPhiInR="";
e158cbc3 470 nameMCPhiInR.Form("MC_Conversion_Mapping_Phi_in_R_%02d",r);
d7d7e825 471 TString titleMCPhiInR="";
472 titleMCPhiInR.Form("MC Mapping of Phi in R%02d",r);
9640a3d1 473 // AddHistogram(nameMCPhiInR, titleMCPhiInR, nXBins, firstX, lastX, xAxisTitle, yAxisTitle);
474 AddHistogram(nameMCPhiInR, titleMCPhiInR, nXBins, -TMath::Pi(), TMath::Pi(), xAxisTitle, yAxisTitle);
d7d7e825 475
9640a3d1 476
477 //Mapping Z in R
478 TString nameMCZInR="";
479 nameMCZInR.Form("MC_Conversion_Mapping_Z_in_R_%02d",r);
480 TString titleMCZInR="";
481 titleMCZInR.Form("MC Mapping of Z in R%02d",r);
482 // AddHistogram(nameMCPhiInR, titleMCPhiInR, nXBins, firstX, lastX, xAxisTitle, yAxisTitle);
483 AddHistogram(nameMCZInR, titleMCZInR, nXBins, -300, 300, xAxisTitle, yAxisTitle);
484
485
486 //Mapping Phi in R Middle Pt
487 TString nameMCMidPtPhiInR="";
488 nameMCMidPtPhiInR.Form("MC_Conversion_Mapping_MidPt_Phi_in_R_%02d",r);
489 TString titleMCMidPtPhiInR="";
490 titleMCMidPtPhiInR.Form("MC Mapping Middle Pt of Phi in R%02d",r);
491 // AddHistogram(nameMCPhiInR, titleMCPhiInR, nXBins, firstX, lastX, xAxisTitle, yAxisTitle);
492 AddHistogram(nameMCMidPtPhiInR, titleMCMidPtPhiInR, nXBins, -TMath::Pi(), TMath::Pi(), xAxisTitle, yAxisTitle);
493
494
495 //Mapping Z in R Middle Pt
496 TString nameMCMidPtZInR="";
497 nameMCMidPtZInR.Form("MC_Conversion_Mapping_MidPt_Z_in_R_%02d",r);
498 TString titleMCMidPtZInR="";
499 titleMCMidPtZInR.Form("MC Mapping Middle Pt of Z in R%02d",r);
500 // AddHistogram(nameMCPhiInR, titleMCPhiInR, nXBins, firstX, lastX, xAxisTitle, yAxisTitle);
501 AddHistogram(nameMCMidPtZInR, titleMCMidPtZInR, nXBins, -300, 300, xAxisTitle, yAxisTitle);
502
503
504
505
d7d7e825 506 //Mapping Phi in R
507 TString nameESDPhiInR="";
e158cbc3 508 nameESDPhiInR.Form("ESD_Conversion_Mapping_Phi_in_R_%02d",r);
d7d7e825 509 TString titleESDPhiInR="";
510 titleESDPhiInR.Form("ESD Mapping of Phi in R%02d",r);
9640a3d1 511 // AddHistogram(nameESDPhiInR, titleESDPhiInR, nXBins, firstX, lastX, xAxisTitle, yAxisTitle);
512 AddHistogram(nameESDPhiInR, titleESDPhiInR, nXBins, -TMath::Pi(), TMath::Pi(), xAxisTitle, yAxisTitle);
513
514 //Mapping Z in R
515 TString nameESDZInR="";
516 nameESDZInR.Form("ESD_Conversion_Mapping_Z_in_R_%02d",r);
517 TString titleESDZInR="";
518 titleESDZInR.Form("ESD Mapping of Z in R%02d",r);
519 // AddHistogram(nameESDPhiInR, titleESDPhiInR, nXBins, firstX, lastX, xAxisTitle, yAxisTitle);
520 AddHistogram(nameESDZInR, titleESDZInR, nXBins, -300, 300, xAxisTitle, yAxisTitle);
521
522 //Mapping Phi in R Middle Pt
523 TString nameESDMidPtPhiInR="";
524 nameESDMidPtPhiInR.Form("ESD_Conversion_Mapping_MidPt_Phi_in_R_%02d",r);
525 TString titleESDMidPtPhiInR="";
526 titleESDMidPtPhiInR.Form("ESD Mapping Middle Pt of Phi in R%02d",r);
527 // AddHistogram(nameESDPhiInR, titleESDPhiInR, nXBins, firstX, lastX, xAxisTitle, yAxisTitle);
528 AddHistogram(nameESDMidPtPhiInR, titleESDMidPtPhiInR, nXBins, -TMath::Pi(), TMath::Pi(), xAxisTitle, yAxisTitle);
529
530 //Mapping Z in R Middle Pt
531 TString nameESDMidPtZInR="";
532 nameESDMidPtZInR.Form("ESD_Conversion_Mapping_MidPt_Z_in_R_%02d",r);
533 TString titleESDMidPtZInR="";
534 titleESDMidPtZInR.Form("ESD Mapping Middle Pt of Z in R%02d",r);
535 // AddHistogram(nameESDPhiInR, titleESDPhiInR, nXBins, firstX, lastX, xAxisTitle, yAxisTitle);
536 AddHistogram(nameESDMidPtZInR, titleESDMidPtZInR, nXBins, -300, 300, xAxisTitle, yAxisTitle);
537
538
539
540 }
541
542
543
544 for(Int_t z =0; z<fNZIndex;z++){
545 //Mapping Phi in Z
546 TString nameMCPhiInZ="";
547 nameMCPhiInZ.Form("MC_Conversion_Mapping_Phi_in_Z_%02d",z);
548 TString titleMCPhiInZ="";
549 titleMCPhiInZ.Form("MC Mapping of Phi in Z%02d",z);
550 // AddHistogram(nameMCPhiInR, titleMCPhiInR, nXBins, firstX, lastX, xAxisTitle, yAxisTitle);
551 AddHistogram(nameMCPhiInZ, titleMCPhiInZ, nXBins, -TMath::Pi(), TMath::Pi(), xAxisTitle, yAxisTitle);
552
553
554 //Mapping R in Z
555 TString nameMCRInZ="";
556 nameMCRInZ.Form("MC_Conversion_Mapping_R_in_Z_%02d",z);
557 TString titleMCRInZ="";
558 titleMCRInZ.Form("MC Mapping of R in Z%02d",z);
559 // AddHistogram(nameMCPhiInR, titleMCPhiInR, nXBins, firstX, lastX, xAxisTitle, yAxisTitle);
560 AddHistogram(nameMCRInZ, titleMCRInZ, nXBins, fMinRadius, fMaxRadius, xAxisTitle, yAxisTitle);
561
562 //Mapping Phi in Z Middle Pt
563 TString nameMCMidPtPhiInZ="";
564 nameMCMidPtPhiInZ.Form("MC_Conversion_Mapping_MidPt_Phi_in_Z_%02d",z);
565 TString titleMCMidPtPhiInZ="";
566 titleMCMidPtPhiInZ.Form("MC Mapping Middle Pt of Phi in Z%02d",z);
567 // AddHistogram(nameMCPhiInR, titleMCPhiInR, nXBins, firstX, lastX, xAxisTitle, yAxisTitle);
568 AddHistogram(nameMCMidPtPhiInZ, titleMCMidPtPhiInZ, nXBins, -TMath::Pi(), TMath::Pi(), xAxisTitle, yAxisTitle);
569
570
571 //Mapping R in Z Middle Pt
572 TString nameMCMidPtRInZ="";
573 nameMCMidPtRInZ.Form("MC_Conversion_Mapping_MidPt_R_in_Z_%02d",z);
574 TString titleMCMidPtRInZ="";
575 titleMCMidPtRInZ.Form("MC Mapping Middle Pt of R in Z%02d",z);
576 // AddHistogram(nameMCPhiInR, titleMCPhiInR, nXBins, firstX, lastX, xAxisTitle, yAxisTitle);
577 AddHistogram(nameMCMidPtRInZ, titleMCMidPtRInZ, nXBins, fMinRadius, fMaxRadius, xAxisTitle, yAxisTitle);
578
579
580
581
582 //Mapping Phi in Z
583 TString nameESDPhiInZ="";
584 nameESDPhiInZ.Form("ESD_Conversion_Mapping_Phi_in_Z_%02d",z);
585 TString titleESDPhiInZ="";
586 titleESDPhiInZ.Form("ESD Mapping of Phi in R%02d",z);
587 // AddHistogram(nameESDPhiInR, titleESDPhiInR, nXBins, firstX, lastX, xAxisTitle, yAxisTitle);
588 AddHistogram(nameESDPhiInZ, titleESDPhiInZ, nXBins, -TMath::Pi(), TMath::Pi(), xAxisTitle, yAxisTitle);
589
590 //Mapping R in Z
591 TString nameESDRInZ="";
592 nameESDRInZ.Form("ESD_Conversion_Mapping_R_in_Z_%02d",z);
593 TString titleESDRInZ="";
594 titleESDRInZ.Form("ESD Mapping of R in Z%02d",z);
595 // AddHistogram(nameESDPhiInR, titleESDPhiInR, nXBins, firstX, lastX, xAxisTitle, yAxisTitle);
596 AddHistogram(nameESDRInZ, titleESDRInZ, nXBins, fMinRadius, fMaxRadius, xAxisTitle, yAxisTitle);
597
598
599 //Mapping Phi in Z Middle Pt
600 TString nameESDMidPtPhiInZ="";
601 nameESDMidPtPhiInZ.Form("ESD_Conversion_Mapping_MidPt_Phi_in_Z_%02d",z);
602 TString titleESDMidPtPhiInZ="";
603 titleESDMidPtPhiInZ.Form("ESD Mapping Middle Ptof Phi in R%02d",z);
604 // AddHistogram(nameESDPhiInR, titleESDPhiInR, nXBins, firstX, lastX, xAxisTitle, yAxisTitle);
605 AddHistogram(nameESDMidPtPhiInZ, titleESDMidPtPhiInZ, nXBins, -TMath::Pi(), TMath::Pi(), xAxisTitle, yAxisTitle);
606
607 //Mapping R in Z Middle Pt
608 TString nameESDMidPtRInZ="";
609 nameESDMidPtRInZ.Form("ESD_Conversion_Mapping_MidPt_R_in_Z_%02d",z);
610 TString titleESDMidPtRInZ="";
611 titleESDMidPtRInZ.Form("ESD Mapping Middle Pt of R in Z%02d",z);
612 // AddHistogram(nameESDPhiInR, titleESDPhiInR, nXBins, firstX, lastX, xAxisTitle, yAxisTitle);
613 AddHistogram(nameESDMidPtRInZ, titleESDMidPtRInZ, nXBins, fMinRadius, fMaxRadius, xAxisTitle, yAxisTitle);
614
615
616
d7d7e825 617 }
9640a3d1 618
619
620
d7d7e825 621}