-
Notifications
You must be signed in to change notification settings - Fork 0
/
CombineNeutralPionResultspp.C
3636 lines (2921 loc) · 232 KB
/
CombineNeutralPionResultspp.C
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
#include <Riostream.h>
#include "TMath.h"
#include <stdlib.h>
#include <fstream>
#include <math.h>
#include <TROOT.h>
#include <TApplication.h>
#include <TPaveLabel.h>
#include <TSystem.h>
#include <TFrame.h>
#include <TStyle.h>
#include <TString.h>
#include "TGaxis.h"
#include "TFile.h"
#include "TH1F.h"
#include "TH1D.h"
#include "TH2F.h"
#include "TF1.h"
#include "TVirtualFitter.h"
#include "TObject.h"
#include "TCanvas.h"
#include "TMultiGraph.h"
#include "TLegend.h"
#include "TDatabasePDG.h"
#include "TMinuit.h"
#include "TBenchmark.h"
#include "TRandom.h"
#include "TLatex.h"
#include "TASImage.h"
#include "TPostScript.h"
#include "TGraphErrors.h"
#include "TArrow.h"
#include "TGraphAsymmErrors.h"
#include "TGaxis.h"
#include "TMarker.h"
#include "Math/WrappedTF1.h"
#include "Math/BrentRootFinder.h"
#include "CommonHeaders/PlottingGammaConversionHistos.h"
#include "CommonHeaders/PlottingGammaConversionAdditional.h"
#include "CommonHeaders/FittingGammaConversion.h"
#include "CommonHeaders/ConversionFunctions.h"
extern TRandom* gRandom;
extern TBenchmark* gBenchmark;
extern TSystem* gSystem;
extern TMinuit* gMinuit;
void ScaleMCYield(TH1D* histoCorrectedToBeScaled, Double_t deltaRapid, Double_t scaling, Double_t nEvtMC, TString nameMeson, TString optionDalitz ){
histoCorrectedToBeScaled->Sumw2();
histoCorrectedToBeScaled->Scale(1./deltaRapid);
histoCorrectedToBeScaled->Scale(scaling);
histoCorrectedToBeScaled->Scale(1./nEvtMC);
for (Int_t i = 1; i < histoCorrectedToBeScaled->GetNbinsX()+1 ; i++){
Double_t newBinContent = histoCorrectedToBeScaled->GetBinContent(i)/histoCorrectedToBeScaled->GetBinCenter(i);
Double_t newBinError = histoCorrectedToBeScaled->GetBinError(i)/histoCorrectedToBeScaled->GetBinCenter(i);
histoCorrectedToBeScaled->SetBinContent(i,newBinContent);
histoCorrectedToBeScaled->SetBinError(i,newBinError);
}
if (nameMeson.CompareTo("Pi0") == 0 ||nameMeson.CompareTo("Pi0EtaBinning") == 0 ){
if (optionDalitz.CompareTo("kFALSE")==0){
histoCorrectedToBeScaled->Scale(1./0.98798);
} else {
histoCorrectedToBeScaled->Scale(1./0.01198);
}
}else{
if (optionDalitz.CompareTo("kFALSE")==0){
histoCorrectedToBeScaled->Scale(1./0.3931);
} else {
histoCorrectedToBeScaled->Scale(1./6.8e-5);
}
}
}
void CombineNeutralPionResultspp(TString suffix = "pdf",
TString nameFilepp1 = "Pythia",
TString nameFilepp2 = "PythiaAddSign",
TString nameFilepp3 = "Phojet",
Bool_t runDrawReweighted = kTRUE){
gROOT->Reset();
gROOT->SetStyle("Plain");
TString dateForOutput = ReturnDateStringForOutput();
TString outputDir = Form("%s/%s/CombineNeutralPionResultspp",suffix.Data(),dateForOutput.Data());
gSystem->Exec("mkdir -p "+outputDir);
gSystem->Exec(Form("cp -r %s %s/InputFilePCMPPPythia.root ",nameFilepp1.Data(),outputDir.Data()));
gSystem->Exec(Form("cp -r %s %s/InputFilePCMPPPythiaPlus.root ",nameFilepp2.Data(),outputDir.Data()));
gSystem->Exec(Form("cp -r %s %s/InputFilePCMPPPhojet.root ",nameFilepp3.Data(),outputDir.Data()));
gSystem->Exec("mkdir -p "+outputDir);
StyleSettingsThesis();
SetPlotStyle();
Color_t colorCombPP = kBlack;
Color_t colorCombMCPP = kRed+3;
Style_t markerStylePP = 33 ;
Style_t markerStylePPMC = 24 ;
Size_t markerSizePP = 2.5;
Color_t colorPi0900GeV = kRed +2;
Color_t colorPi02760GeV = kMagenta+2;
Color_t colorPi07TeV = kBlue+2;
Color_t colorPi08TeV = kGreen+2;
Color_t colorPi0900GeVBox = colorPi0900GeV-10;
Color_t colorPi02760GeVBox = colorPi02760GeV-10;
Color_t colorPi07TeVBox = colorPi07TeV-10;
Color_t colorMCPythiaPP900GeV = colorPi0900GeV-4;
Color_t colorMCPythiaPP2760GeV = colorPi02760GeV+2;
Color_t colorMCPythiaPP7TeV = colorPi07TeV+3;
Color_t colorMCPythiaPP8TeV = colorPi08TeV+2;
Color_t colorMCPhojetPP900GeV = colorPi0900GeV+2;
Color_t colorMCPhojetPP2760GeV = colorPi02760GeV-4;
Color_t colorMCPhojetPP7TeV = colorPi07TeV-3;
Color_t colorMCPhojetPP8TeV = colorPi08TeV-2;
Style_t markerStyleSpectrum7TeVMC = 24 ;
Style_t markerStyleSpectrum900GeVMC = 25 ;
Style_t markerStyleSpectrum8TeVMC = 22 ;
Style_t markerStyleSpectrum2760GeVMC = 30 ;
Style_t markerStyleSpectrum7TeV = 20 ;
Style_t markerStyleSpectrum900GeV = 21 ;
Style_t markerStyleSpectrum8TeV = 24 ;
Style_t markerStyleSpectrum2760GeV = 29 ;
Double_t xSection7TeVppINEL = 73.2*1e9;
Double_t xSection2760GeVppINEL = 62.8*1e9;
Double_t xSection900GeVppINEL = 52.5*1e9;
Style_t markerStyleMCPP8TeV = 24 ;
Style_t markerStyleMCPP7TeV = 24 ;
Style_t markerStyleMCPP900GeV = 25 ;
Style_t markerStyleMCPP2760GeV = 30 ;
Size_t markerSizePi0PP8TeV = 1.8;
Size_t markerSizePi0PP7TeV = 1.8;
Size_t markerSizePi0PP900GeV = 1.8;
Size_t markerSizePi0PP2760GeV = 2.2;
TString collisionSystemPP2760GeV = "pp #sqrt{#it{s}} = 2.76 TeV";
TString collisionSystemPP7TeV = "pp #sqrt{#it{s}} = 7 TeV";
TString collisionSystemPP8TeV = "pp #sqrt{#it{s}} = 8 TeV";
TString collisionSystemPP900GeV = "pp #sqrt{#it{s}} = 0.9 TeV";
Size_t markerSizeComparison = 0.5;
Double_t maxPtMesonEffFit = 12.;
Double_t minPtMesonEffFit = 1.2;
Int_t offsetCorrectionHighPt= 1;
TF1* fitTrueEffi = new TF1("EffiFitDummy","1 - [0]*exp([1]*x)+[1]");
fitTrueEffi->SetRange(minPtMesonEffFit,maxPtMesonEffFit);
//******************** Retrieving plots from PP file ************************
//************************** PYTHIA only *************************************
TFile* filePythia = new TFile(nameFilepp1.Data());
//********* 8 TeV *****************************************************************************************************
TDirectory* directoryPythiaPi08TeV = (TDirectory*)filePythia->Get("Pi08TeV");
TH1D* histoPi0CorrectedSpecPythia8TeV = (TH1D*)directoryPythiaPi08TeV->Get("CorrectedYieldPi0");
TGraphAsymmErrors* graphPi0CorrectedSpecSysPythia8TeV = (TGraphAsymmErrors*)directoryPythiaPi08TeV->Get("Pi0SystError");
TH1D* histoPi0MassDataPythia8TeV = (TH1D*)directoryPythiaPi08TeV->Get("MassPi0");
TH1D* histoPi0MassMCPythia8TeV = (TH1D*)directoryPythiaPi08TeV->Get("TrueMassPi0");
TH1D* histoPi0WidthDataPythia8TeV = (TH1D*)directoryPythiaPi08TeV->Get("FWHMPi0MeV");
TH1D* histoPi0WidthMCPythia8TeV = (TH1D*)directoryPythiaPi08TeV->Get("TrueFWHMPi0MeV");
TH1D* histoPi0TrueEffiPtPythia8TeV = (TH1D*)directoryPythiaPi08TeV->Get("EfficiencyPi0");
TH1D* histoPi0AcceptPtPythia8TeV = (TH1D*)directoryPythiaPi08TeV->Get("AcceptancePi0");
TH1D* histoMCPi0YieldPtPythia8TeV = (TH1D*)directoryPythiaPi08TeV->Get("Pi0_Input_Reweighted");
TH1D* histoMCPi0YieldPtPythia8TeVWOWeights = (TH1D*)directoryPythiaPi08TeV->Get("Pi0_Input");
TH1D* histoPi0WeightsPythia8TeV = (TH1D*)directoryPythiaPi08TeV->Get("Pi0_Weights");
TH1D* histoPi0RawYieldPythia8TeV = (TH1D*)directoryPythiaPi08TeV->Get("RAWYieldPerEventsPi0");
histoPi0MassDataPythia8TeV->Scale(1000.);
histoPi0MassMCPythia8TeV->Scale(1000.);
TDirectory* directoryPythiaEta8TeV = (TDirectory*)filePythia->Get("Eta_pp_8TeV");
TH1D* histoEtaCorrectedSpecPythia8TeV = (TH1D*)directoryPythiaEta8TeV->Get("CorrectedYieldEta");
TGraphAsymmErrors* graphEtaCorrectedSpecSysPythia8TeV = (TGraphAsymmErrors*)directoryPythiaEta8TeV->Get("EtaSystError");
TH1D* histoEtaMassDataPythia8TeV = (TH1D*)directoryPythiaEta8TeV->Get("MassEta");
TH1D* histoEtaMassMCPythia8TeV = (TH1D*)directoryPythiaEta8TeV->Get("TrueMassEta");
TH1D* histoEtaWidthDataPythia8TeV = (TH1D*)directoryPythiaEta8TeV->Get("FWHMEtaMeV");
TH1D* histoEtaWidthMCPythia8TeV = (TH1D*)directoryPythiaEta8TeV->Get("TrueFWHMEtaMeV");
TH1D* histoEtaTrueEffiPtPythia8TeV = (TH1D*)directoryPythiaEta8TeV->Get("EfficiencyEta");
TH1D* histoEtaAcceptPtPythia8TeV = (TH1D*)directoryPythiaEta8TeV->Get("EfficiencyEta");
TH1D* histoMCEtaYieldPtPythia8TeV = (TH1D*)directoryPythiaEta8TeV->Get("Eta_Input_Reweighted");
TH1D* histoMCEtaYieldPtPythia8TeVWOWeights = (TH1D*)directoryPythiaEta8TeV->Get("Eta_Input");
TH1D* histoEtaWeightsPythia8TeV = (TH1D*)directoryPythiaEta8TeV->Get("Eta_Weights");
TH1D* histoEtaRawYieldPythia8TeV = (TH1D*)directoryPythiaEta8TeV->Get("RAWYieldPerEventsEta");
histoEtaMassDataPythia8TeV->Scale(1000.);
histoEtaMassMCPythia8TeV->Scale(1000.);
//********* 7 TeV *****************************************************************************************************
TDirectory* directoryPythiaPi07TeV = (TDirectory*)filePythia->Get("Pi0_pp_7TeV");
TH1D* histoPi0CorrectedSpecPythia7TeV = (TH1D*)directoryPythiaPi07TeV->Get("CorrectedYieldPi0");
TGraphAsymmErrors* graphPi0CorrectedSpecSysPythia7TeV = (TGraphAsymmErrors*)directoryPythiaPi07TeV->Get("Pi0SystError");
TH1D* histoPi0MassDataPythia7TeV = (TH1D*)directoryPythiaPi07TeV->Get("MassPi0");
TH1D* histoPi0MassMCPythia7TeV = (TH1D*)directoryPythiaPi07TeV->Get("TrueMassPi0");
TH1D* histoPi0WidthDataPythia7TeV = (TH1D*)directoryPythiaPi07TeV->Get("FWHMPi0MeV");
TH1D* histoPi0WidthMCPythia7TeV = (TH1D*)directoryPythiaPi07TeV->Get("TrueFWHMPi0MeV");
TH1D* histoPi0TrueEffiPtPythia7TeV = (TH1D*)directoryPythiaPi07TeV->Get("EfficiencyPi0");
TH1D* histoPi0AcceptPtPythia7TeV = (TH1D*)directoryPythiaPi07TeV->Get("AcceptancePi0");
TH1D* histoMCPi0YieldPtPythia7TeV = (TH1D*)directoryPythiaPi07TeV->Get("Pi0_Input_Reweighted");
TH1D* histoMCPi0YieldPtPythia7TeVWOWeights = (TH1D*)directoryPythiaPi07TeV->Get("Pi0_Input");
TH1D* histoPi0WeightsPythia7TeV = (TH1D*)directoryPythiaPi07TeV->Get("Pi0_Weights");
TH1D* histoPi0RawYieldPythia7TeV = (TH1D*)directoryPythiaPi07TeV->Get("RAWYieldPerEventsPi0");
histoPi0MassDataPythia7TeV->Scale(1000.);
histoPi0MassMCPythia7TeV->Scale(1000.);
TDirectory* directoryPythiaEta7TeV = (TDirectory*)filePythia->Get("Eta_pp_7TeV");
TH1D* histoEtaCorrectedSpecPythia7TeV = (TH1D*)directoryPythiaEta7TeV->Get("CorrectedYieldEta");
TGraphAsymmErrors* graphEtaCorrectedSpecSysPythia7TeV = (TGraphAsymmErrors*)directoryPythiaEta7TeV->Get("EtaSystError");
TH1D* histoEtaMassDataPythia7TeV = (TH1D*)directoryPythiaEta7TeV->Get("MassEta");
TH1D* histoEtaMassMCPythia7TeV = (TH1D*)directoryPythiaEta7TeV->Get("TrueMassEta");
TH1D* histoEtaWidthDataPythia7TeV = (TH1D*)directoryPythiaEta7TeV->Get("FWHMEtaMeV");
TH1D* histoEtaWidthMCPythia7TeV = (TH1D*)directoryPythiaEta7TeV->Get("TrueFWHMEtaMeV");
TH1D* histoEtaTrueEffiPtPythia7TeV = (TH1D*)directoryPythiaEta7TeV->Get("EfficiencyEta");
TH1D* histoEtaAcceptPtPythia7TeV = (TH1D*)directoryPythiaEta7TeV->Get("EfficiencyEta");
TH1D* histoMCEtaYieldPtPythia7TeV = (TH1D*)directoryPythiaEta7TeV->Get("Eta_Input_Reweighted");
TH1D* histoMCEtaYieldPtPythia7TeVWOWeights = (TH1D*)directoryPythiaEta7TeV->Get("Eta_Input");
TH1D* histoEtaWeightsPythia7TeV = (TH1D*)directoryPythiaEta7TeV->Get("Eta_Weights");
TH1D* histoEtaRawYieldPythia7TeV = (TH1D*)directoryPythiaEta7TeV->Get("RAWYieldPerEventsEta");
histoEtaMassDataPythia7TeV->Scale(1000.);
histoEtaMassMCPythia7TeV->Scale(1000.);
//********* 2.76 TeV *****************************************************************************************************
TDirectory* directoryPythiaPi02760GeV = (TDirectory*)filePythia->Get("Pi0_pp_2760GeV");
TH1D* histoPi0CorrectedSpecPythia2760GeV = (TH1D*)directoryPythiaPi02760GeV->Get("CorrectedYieldPi0");
TGraphAsymmErrors* graphPi0CorrectedSpecSysPythia2760GeV = (TGraphAsymmErrors*)directoryPythiaPi02760GeV->Get("Pi0SystError");
TH1D* histoPi0MassDataPythia2760GeV = (TH1D*)directoryPythiaPi02760GeV->Get("MassPi0");
TH1D* histoPi0MassMCPythia2760GeV = (TH1D*)directoryPythiaPi02760GeV->Get("TrueMassPi0");
TH1D* histoPi0WidthDataPythia2760GeV = (TH1D*)directoryPythiaPi02760GeV->Get("FWHMPi0MeV");
TH1D* histoPi0WidthMCPythia2760GeV = (TH1D*)directoryPythiaPi02760GeV->Get("TrueFWHMPi0MeV");
TH1D* histoPi0TrueEffiPtPythia2760GeV = (TH1D*)directoryPythiaPi02760GeV->Get("EfficiencyPi0");
TH1D* histoPi0AcceptPtPythia2760GeV = (TH1D*)directoryPythiaPi02760GeV->Get("AcceptancePi0");
TH1D* histoMCPi0YieldPtPythia2760GeV = (TH1D*)directoryPythiaPi02760GeV->Get("Pi0_Input_Reweighted");
TH1D* histoMCPi0YieldPtPythia2760GeVWOWeights = (TH1D*)directoryPythiaPi02760GeV->Get("Pi0_Input");
TH1D* histoPi0WeightsPythia2760GeV = (TH1D*)directoryPythiaPi02760GeV->Get("Pi0_Weights");
TH1D* histoPi0RawYieldPythia2760GeV = (TH1D*)directoryPythiaPi02760GeV->Get("RAWYieldPerEventsPi0");
histoPi0MassDataPythia2760GeV->Scale(1000.);
histoPi0MassMCPythia2760GeV->Scale(1000.);
TDirectory* directoryPythiaEta2760GeV = (TDirectory*)filePythia->Get("Eta_pp_2760GeV");
TH1D* histoEtaCorrectedSpecPythia2760GeV = (TH1D*)directoryPythiaEta2760GeV->Get("CorrectedYieldEta");
TGraphAsymmErrors* graphEtaCorrectedSpecSysPythia2760GeV = (TGraphAsymmErrors*)directoryPythiaEta2760GeV->Get("EtaSystError");
TH1D* histoEtaMassDataPythia2760GeV = (TH1D*)directoryPythiaEta2760GeV->Get("MassEta");
TH1D* histoEtaMassMCPythia2760GeV = (TH1D*)directoryPythiaEta2760GeV->Get("TrueMassEta");
TH1D* histoEtaWidthDataPythia2760GeV = (TH1D*)directoryPythiaEta2760GeV->Get("FWHMEtaMeV");
TH1D* histoEtaWidthMCPythia2760GeV = (TH1D*)directoryPythiaEta2760GeV->Get("TrueFWHMEtaMeV");
TH1D* histoEtaTrueEffiPtPythia2760GeV = (TH1D*)directoryPythiaEta2760GeV->Get("EfficiencyEta");
TH1D* histoEtaAcceptPtPythia2760GeV = (TH1D*)directoryPythiaEta2760GeV->Get("EfficiencyEta");
TH1D* histoMCEtaYieldPtPythia2760GeV = (TH1D*)directoryPythiaEta2760GeV->Get("Eta_Input_Reweighted");
TH1D* histoMCEtaYieldPtPythia2760GeVWOWeights = (TH1D*)directoryPythiaEta2760GeV->Get("Eta_Input");
TH1D* histoEtaWeightsPythia2760GeV = (TH1D*)directoryPythiaEta2760GeV->Get("Eta_Weights");
TH1D* histoEtaRawYieldPythia2760GeV = (TH1D*)directoryPythiaEta2760GeV->Get("RAWYieldPerEventsEta");
histoEtaMassDataPythia2760GeV->Scale(1000.);
histoEtaMassMCPythia2760GeV->Scale(1000.);
//********* 0.9 TeV *****************************************************************************************************
TDirectory* directoryPythiaPi0900GeV = (TDirectory*)filePythia->Get("Pi0_pp_900GeV");
TH1D* histoPi0CorrectedSpecPythia900GeV = (TH1D*)directoryPythiaPi0900GeV->Get("CorrectedYieldPi0");
TGraphAsymmErrors* graphPi0CorrectedSpecSysPythia900GeV = (TGraphAsymmErrors*)directoryPythiaPi0900GeV->Get("Pi0SystError");
TH1D* histoPi0MassDataPythia900GeV = (TH1D*)directoryPythiaPi0900GeV->Get("MassPi0");
TH1D* histoPi0MassMCPythia900GeV = (TH1D*)directoryPythiaPi0900GeV->Get("TrueMassPi0");
TH1D* histoPi0WidthDataPythia900GeV = (TH1D*)directoryPythiaPi0900GeV->Get("FWHMPi0MeV");
TH1D* histoPi0WidthMCPythia900GeV = (TH1D*)directoryPythiaPi0900GeV->Get("TrueFWHMPi0MeV");
TH1D* histoPi0TrueEffiPtPythia900GeV = (TH1D*)directoryPythiaPi0900GeV->Get("EfficiencyPi0");
TH1D* histoPi0AcceptPtPythia900GeV = (TH1D*)directoryPythiaPi0900GeV->Get("AcceptancePi0");
TH1D* histoMCPi0YieldPtPythia900GeV = (TH1D*)directoryPythiaPi0900GeV->Get("Pi0_Input_Reweighted");
TH1D* histoMCPi0YieldPtPythia900GeVWOWeights = (TH1D*)directoryPythiaPi0900GeV->Get("Pi0_Input");
TH1D* histoPi0WeightsPythia900GeV = (TH1D*)directoryPythiaPi0900GeV->Get("Pi0_Weights");
TH1D* histoPi0RawYieldPythia900GeV = (TH1D*)directoryPythiaPi0900GeV->Get("RAWYieldPerEventsPi0");
histoPi0MassDataPythia900GeV->Scale(1000.);
histoPi0MassMCPythia900GeV->Scale(1000.);
TDirectory* directoryPythiaEta900GeV = (TDirectory*)filePythia->Get("Eta_pp_900GeV");
TH1D* histoEtaCorrectedSpecPythia900GeV = (TH1D*)directoryPythiaEta900GeV->Get("CorrectedYieldEta");
TGraphAsymmErrors* graphEtaCorrectedSpecSysPythia900GeV = (TGraphAsymmErrors*)directoryPythiaEta900GeV->Get("EtaSystError");
TH1D* histoEtaMassDataPythia900GeV = (TH1D*)directoryPythiaEta900GeV->Get("MassEta");
TH1D* histoEtaMassMCPythia900GeV = (TH1D*)directoryPythiaEta900GeV->Get("TrueMassEta");
TH1D* histoEtaWidthDataPythia900GeV = (TH1D*)directoryPythiaEta900GeV->Get("FWHMEtaMeV");
TH1D* histoEtaWidthMCPythia900GeV = (TH1D*)directoryPythiaEta900GeV->Get("TrueFWHMEtaMeV");
TH1D* histoEtaTrueEffiPtPythia900GeV = (TH1D*)directoryPythiaEta900GeV->Get("EfficiencyEta");
TH1D* histoEtaAcceptPtPythia900GeV = (TH1D*)directoryPythiaEta900GeV->Get("EfficiencyEta");
TH1D* histoMCEtaYieldPtPythia900GeV = (TH1D*)directoryPythiaEta900GeV->Get("Eta_Input_Reweighted");
TH1D* histoMCEtaYieldPtPythia900GeVWOWeights = (TH1D*)directoryPythiaEta900GeV->Get("Eta_Input");
TH1D* histoEtaWeightsPythia900GeV = (TH1D*)directoryPythiaEta900GeV->Get("Eta_Weights");
TH1D* histoEtaRawYieldPythia900GeV = (TH1D*)directoryPythiaEta900GeV->Get("RAWYieldPerEventsEta");
histoEtaMassDataPythia900GeV->Scale(1000.);
histoEtaMassMCPythia900GeV->Scale(1000.);
//*********************************************************************************************
//************************** PYTHIA plus Added Signal *************************************
TFile* filePythiaPlus = new TFile(nameFilepp2.Data());
//********* 8 TeV *****************************************************************************************************
TDirectory* directoryPythiaPlusPi08TeV = (TDirectory*)filePythiaPlus->Get("Pi0_pp_8TeV");
TH1D* histoPi0CorrectedSpecPythiaPlus8TeV = (TH1D*)directoryPythiaPlusPi08TeV->Get("CorrectedYieldPi0");
TGraphAsymmErrors* graphPi0CorrectedSpecSysPythiaPlus8TeV = (TGraphAsymmErrors*)directoryPythiaPlusPi08TeV->Get("Pi0SystError");
TH1D* histoPi0MassDataPythiaPlus8TeV = (TH1D*)directoryPythiaPlusPi08TeV->Get("MassPi0");
TH1D* histoPi0MassMCPythiaPlus8TeV = (TH1D*)directoryPythiaPlusPi08TeV->Get("TrueMassPi0");
TH1D* histoPi0WidthDataPythiaPlus8TeV = (TH1D*)directoryPythiaPlusPi08TeV->Get("FWHMPi0MeV");
TH1D* histoPi0WidthMCPythiaPlus8TeV = (TH1D*)directoryPythiaPlusPi08TeV->Get("TrueFWHMPi0MeV");
TH1D* histoPi0TrueEffiPtPythiaPlus8TeV = (TH1D*)directoryPythiaPlusPi08TeV->Get("EfficiencyPi0");
TH1D* histoPi0AcceptPtPythiaPlus8TeV = (TH1D*)directoryPythiaPlusPi08TeV->Get("AcceptancePi0");
TH1D* histoMCPi0YieldPtPythiaPlus8TeV = (TH1D*)directoryPythiaPlusPi08TeV->Get("Pi0_Input_Reweighted");
TH1D* histoMCPi0YieldPtPythiaPlus8TeVWOWeights = (TH1D*)directoryPythiaPlusPi08TeV->Get("Pi0_Input");
TH1D* histoPi0WeightsPythiaPlus8TeV = (TH1D*)directoryPythiaPlusPi08TeV->Get("Pi0_Weights");
TH1D* histoPi0RawYieldPythiaPlus8TeV = (TH1D*)directoryPythiaPlusPi08TeV->Get("RAWYieldPerEventsPi0");
TH1D* histoMCPi0YieldPtPythiaPlus8TeVAddedSig = (TH1D*)directoryPythiaPlusPi08TeV->Get("Pi0_Input_Reweighted_AddedSig");
TH1D* histoMCPi0YieldPtPythiaPlus8TeVAddedSigWOWeights = (TH1D*)directoryPythiaPlusPi08TeV->Get("Pi0_Input_AddedSig");
histoPi0MassDataPythiaPlus8TeV->Scale(1000.);
histoPi0MassMCPythiaPlus8TeV->Scale(1000.);
TDirectory* directoryPythiaPlusEta8TeV = (TDirectory*)filePythiaPlus->Get("Eta_pp_8TeV");
TH1D* histoEtaCorrectedSpecPythiaPlus8TeV = (TH1D*)directoryPythiaPlusEta8TeV->Get("CorrectedYieldEta");
TGraphAsymmErrors* graphEtaCorrectedSpecSysPythiaPlus8TeV = (TGraphAsymmErrors*)directoryPythiaPlusEta8TeV->Get("EtaSystError");
TH1D* histoEtaMassDataPythiaPlus8TeV = (TH1D*)directoryPythiaPlusEta8TeV->Get("MassEta");
TH1D* histoEtaMassMCPythiaPlus8TeV = (TH1D*)directoryPythiaPlusEta8TeV->Get("TrueMassEta");
TH1D* histoEtaWidthDataPythiaPlus8TeV = (TH1D*)directoryPythiaPlusEta8TeV->Get("FWHMEtaMeV");
TH1D* histoEtaWidthMCPythiaPlus8TeV = (TH1D*)directoryPythiaPlusEta8TeV->Get("TrueFWHMEtaMeV");
TH1D* histoEtaTrueEffiPtPythiaPlus8TeV = (TH1D*)directoryPythiaPlusEta8TeV->Get("EfficiencyEta");
TH1D* histoEtaAcceptPtPythiaPlus8TeV = (TH1D*)directoryPythiaPlusEta8TeV->Get("EfficiencyEta");
TH1D* histoMCEtaYieldPtPythiaPlus8TeV = (TH1D*)directoryPythiaPlusEta8TeV->Get("Eta_Input_Reweighted");
TH1D* histoMCEtaYieldPtPythiaPlus8TeVWOWeights = (TH1D*)directoryPythiaPlusEta8TeV->Get("Eta_Input");
TH1D* histoEtaWeightsPythiaPlus8TeV = (TH1D*)directoryPythiaPlusEta8TeV->Get("Eta_Weights");
TH1D* histoMCEtaYieldPtPythiaPlus8TeVAddedSig = (TH1D*)directoryPythiaPlusEta8TeV->Get("Eta_Input_Reweighted_AddedSig");
TH1D* histoMCEtaYieldPtPythiaPlus8TeVAddedSigWOWeights = (TH1D*)directoryPythiaPlusEta8TeV->Get("Eta_Input_AddedSig");
//TH1D* histoEtaRawYieldPythiaPlus8TeV = (TH1D*)directoryPythiaPlusEta8TeV->Get("RAWYieldPerEventsEta");
histoEtaMassDataPythiaPlus8TeV->Scale(1000.);
histoEtaMassMCPythiaPlus8TeV->Scale(1000.);
//********* 7 TeV *****************************************************************************************************
TDirectory* directoryPythiaPlusPi07TeV = (TDirectory*)filePythiaPlus->Get("Pi0_pp_7TeV");
TH1D* histoPi0CorrectedSpecPythiaPlus7TeV = (TH1D*)directoryPythiaPlusPi07TeV->Get("CorrectedYieldPi0");
TGraphAsymmErrors* graphPi0CorrectedSpecSysPythiaPlus7TeV = (TGraphAsymmErrors*)directoryPythiaPlusPi07TeV->Get("Pi0SystError");
TH1D* histoPi0MassDataPythiaPlus7TeV = (TH1D*)directoryPythiaPlusPi07TeV->Get("MassPi0");
TH1D* histoPi0MassMCPythiaPlus7TeV = (TH1D*)directoryPythiaPlusPi07TeV->Get("TrueMassPi0");
TH1D* histoPi0WidthDataPythiaPlus7TeV = (TH1D*)directoryPythiaPlusPi07TeV->Get("FWHMPi0MeV");
TH1D* histoPi0WidthMCPythiaPlus7TeV = (TH1D*)directoryPythiaPlusPi07TeV->Get("TrueFWHMPi0MeV");
TH1D* histoPi0TrueEffiPtPythiaPlus7TeV = (TH1D*)directoryPythiaPlusPi07TeV->Get("EfficiencyPi0");
TH1D* histoPi0AcceptPtPythiaPlus7TeV = (TH1D*)directoryPythiaPlusPi07TeV->Get("AcceptancePi0");
TH1D* histoMCPi0YieldPtPythiaPlus7TeV = (TH1D*)directoryPythiaPlusPi07TeV->Get("Pi0_Input_Reweighted");
TH1D* histoMCPi0YieldPtPythiaPlus7TeVWOWeights = (TH1D*)directoryPythiaPlusPi07TeV->Get("Pi0_Input");
TH1D* histoPi0WeightsPythiaPlus7TeV = (TH1D*)directoryPythiaPlusPi07TeV->Get("Pi0_Weights");
TH1D* histoMCPi0YieldPtPythiaPlus7TeVAddedSig = (TH1D*)directoryPythiaPlusPi07TeV->Get("Pi0_Input_Reweighted_AddedSig");
TH1D* histoMCPi0YieldPtPythiaPlus7TeVAddedSigWOWeights = (TH1D*)directoryPythiaPlusPi07TeV->Get("Pi0_Input_AddedSig");
TH1D* histoPi0RawYieldPythiaPlus7TeV = (TH1D*)directoryPythiaPlusPi07TeV->Get("RAWYieldPerEventsPi0");
histoPi0MassDataPythiaPlus7TeV->Scale(1000.);
histoPi0MassMCPythiaPlus7TeV->Scale(1000.);
TDirectory* directoryPythiaPlusEta7TeV = (TDirectory*)filePythiaPlus->Get("Eta_pp_7TeV");
TH1D* histoEtaCorrectedSpecPythiaPlus7TeV = (TH1D*)directoryPythiaPlusEta7TeV->Get("CorrectedYieldEta");
TGraphAsymmErrors* graphEtaCorrectedSpecSysPythiaPlus7TeV = (TGraphAsymmErrors*)directoryPythiaPlusEta7TeV->Get("EtaSystError");
TH1D* histoEtaMassDataPythiaPlus7TeV = (TH1D*)directoryPythiaPlusEta7TeV->Get("MassEta");
TH1D* histoEtaMassMCPythiaPlus7TeV = (TH1D*)directoryPythiaPlusEta7TeV->Get("TrueMassEta");
TH1D* histoEtaWidthDataPythiaPlus7TeV = (TH1D*)directoryPythiaPlusEta7TeV->Get("FWHMEtaMeV");
TH1D* histoEtaWidthMCPythiaPlus7TeV = (TH1D*)directoryPythiaPlusEta7TeV->Get("TrueFWHMEtaMeV");
TH1D* histoEtaTrueEffiPtPythiaPlus7TeV = (TH1D*)directoryPythiaPlusEta7TeV->Get("EfficiencyEta");
TH1D* histoEtaAcceptPtPythiaPlus7TeV = (TH1D*)directoryPythiaPlusEta7TeV->Get("EfficiencyEta");
TH1D* histoMCEtaYieldPtPythiaPlus7TeV = (TH1D*)directoryPythiaPlusEta7TeV->Get("Eta_Input_Reweighted");
TH1D* histoMCEtaYieldPtPythiaPlus7TeVWOWeights = (TH1D*)directoryPythiaPlusEta7TeV->Get("Eta_Input");
TH1D* histoEtaWeightsPythiaPlus7TeV = (TH1D*)directoryPythiaPlusEta7TeV->Get("Eta_Weights");
TH1D* histoMCEtaYieldPtPythiaPlus7TeVAddedSig = (TH1D*)directoryPythiaPlusEta7TeV->Get("Eta_Input_Reweighted_AddedSig");
TH1D* histoMCEtaYieldPtPythiaPlus7TeVAddedSigWOWeights = (TH1D*)directoryPythiaPlusEta7TeV->Get("Eta_Input_AddedSig");
TH1D* histoEtaRawYieldPythiaPlus7TeV = (TH1D*)directoryPythiaPlusEta7TeV->Get("RAWYieldPerEventsEta");
histoEtaMassDataPythiaPlus7TeV->Scale(1000.);
histoEtaMassMCPythiaPlus7TeV->Scale(1000.);
//********* 2.76 TeV *****************************************************************************************************
TDirectory* directoryPythiaPlusPi02760GeV = (TDirectory*)filePythiaPlus->Get("Pi0_pp_2760GeV");
TH1D* histoPi0CorrectedSpecPythiaPlus2760GeV = (TH1D*)directoryPythiaPlusPi02760GeV->Get("CorrectedYieldPi0");
TGraphAsymmErrors* graphPi0CorrectedSpecSysPythiaPlus2760GeV = (TGraphAsymmErrors*)directoryPythiaPlusPi02760GeV->Get("Pi0SystError");
TH1D* histoPi0MassDataPythiaPlus2760GeV = (TH1D*)directoryPythiaPlusPi02760GeV->Get("MassPi0");
TH1D* histoPi0MassMCPythiaPlus2760GeV = (TH1D*)directoryPythiaPlusPi02760GeV->Get("TrueMassPi0");
TH1D* histoPi0WidthDataPythiaPlus2760GeV = (TH1D*)directoryPythiaPlusPi02760GeV->Get("FWHMPi0MeV");
TH1D* histoPi0WidthMCPythiaPlus2760GeV = (TH1D*)directoryPythiaPlusPi02760GeV->Get("TrueFWHMPi0MeV");
TH1D* histoPi0TrueEffiPtPythiaPlus2760GeV = (TH1D*)directoryPythiaPlusPi02760GeV->Get("EfficiencyPi0");
TH1D* histoPi0AcceptPtPythiaPlus2760GeV = (TH1D*)directoryPythiaPlusPi02760GeV->Get("AcceptancePi0");
TH1D* histoMCPi0YieldPtPythiaPlus2760GeV = (TH1D*)directoryPythiaPlusPi02760GeV->Get("Pi0_Input_Reweighted");
TH1D* histoMCPi0YieldPtPythiaPlus2760GeVWOWeights = (TH1D*)directoryPythiaPlusPi02760GeV->Get("Pi0_Input");
TH1D* histoPi0WeightsPythiaPlus2760GeV = (TH1D*)directoryPythiaPlusPi02760GeV->Get("Pi0_Weights");
TH1D* histoMCPi0YieldPtPythiaPlus2760GeVAddedSig = (TH1D*)directoryPythiaPlusPi02760GeV->Get("Pi0_Input_Reweighted_AddedSig");
TH1D* histoMCPi0YieldPtPythiaPlus2760GeVAddedSigWOWeights = (TH1D*)directoryPythiaPlusPi02760GeV->Get("Pi0_Input_AddedSig");
TH1D* histoPi0RawYieldPythiaPlus2760GeV = (TH1D*)directoryPythiaPlusPi02760GeV->Get("RAWYieldPerEventsPi0");
histoPi0MassDataPythiaPlus2760GeV->Scale(1000.);
histoPi0MassMCPythiaPlus2760GeV->Scale(1000.);
TDirectory* directoryPythiaPlusEta2760GeV = (TDirectory*)filePythiaPlus->Get("Eta_pp_2760GeV");
TH1D* histoEtaCorrectedSpecPythiaPlus2760GeV = (TH1D*)directoryPythiaPlusEta2760GeV->Get("CorrectedYieldEta");
TGraphAsymmErrors* graphEtaCorrectedSpecSysPythiaPlus2760GeV = (TGraphAsymmErrors*)directoryPythiaPlusEta2760GeV->Get("EtaSystError");
TH1D* histoEtaMassDataPythiaPlus2760GeV = (TH1D*)directoryPythiaPlusEta2760GeV->Get("MassEta");
TH1D* histoEtaMassMCPythiaPlus2760GeV = (TH1D*)directoryPythiaPlusEta2760GeV->Get("TrueMassEta");
TH1D* histoEtaWidthDataPythiaPlus2760GeV = (TH1D*)directoryPythiaPlusEta2760GeV->Get("FWHMEtaMeV");
TH1D* histoEtaWidthMCPythiaPlus2760GeV = (TH1D*)directoryPythiaPlusEta2760GeV->Get("TrueFWHMEtaMeV");
TH1D* histoEtaTrueEffiPtPythiaPlus2760GeV = (TH1D*)directoryPythiaPlusEta2760GeV->Get("EfficiencyEta");
TH1D* histoEtaAcceptPtPythiaPlus2760GeV = (TH1D*)directoryPythiaPlusEta2760GeV->Get("EfficiencyEta");
TH1D* histoMCEtaYieldPtPythiaPlus2760GeV = (TH1D*)directoryPythiaPlusEta2760GeV->Get("Eta_Input_Reweighted");
TH1D* histoMCEtaYieldPtPythiaPlus2760GeVWOWeights = (TH1D*)directoryPythiaPlusEta2760GeV->Get("Eta_Input");
TH1D* histoEtaWeightsPythiaPlus2760GeV = (TH1D*)directoryPythiaPlusEta2760GeV->Get("Eta_Weights");
TH1D* histoMCEtaYieldPtPythiaPlus2760GeVAddedSig = (TH1D*)directoryPythiaPlusEta2760GeV->Get("Eta_Input_Reweighted_AddedSig");
TH1D* histoMCEtaYieldPtPythiaPlus2760GeVAddedSigWOWeights = (TH1D*)directoryPythiaPlusEta2760GeV->Get("Eta_Input_AddedSig");
TH1D* histoEtaRawYieldPythiaPlus2760GeV = (TH1D*)directoryPythiaPlusEta2760GeV->Get("RAWYieldPerEventsEta");
histoEtaMassDataPythiaPlus2760GeV->Scale(1000.);
histoEtaMassMCPythiaPlus2760GeV->Scale(1000.);
//********* 0.9 TeV *****************************************************************************************************
TDirectory* directoryPythiaPlusPi0900GeV = (TDirectory*)filePythiaPlus->Get("Pi0_pp_900GeV");
TH1D* histoPi0CorrectedSpecPythiaPlus900GeV = (TH1D*)directoryPythiaPlusPi0900GeV->Get("CorrectedYieldPi0");
TGraphAsymmErrors* graphPi0CorrectedSpecSysPythiaPlus900GeV = (TGraphAsymmErrors*)directoryPythiaPlusPi0900GeV->Get("Pi0SystError");
TH1D* histoPi0MassDataPythiaPlus900GeV = (TH1D*)directoryPythiaPlusPi0900GeV->Get("MassPi0");
TH1D* histoPi0MassMCPythiaPlus900GeV = (TH1D*)directoryPythiaPlusPi0900GeV->Get("TrueMassPi0");
TH1D* histoPi0WidthDataPythiaPlus900GeV = (TH1D*)directoryPythiaPlusPi0900GeV->Get("FWHMPi0MeV");
TH1D* histoPi0WidthMCPythiaPlus900GeV = (TH1D*)directoryPythiaPlusPi0900GeV->Get("TrueFWHMPi0MeV");
TH1D* histoPi0TrueEffiPtPythiaPlus900GeV = (TH1D*)directoryPythiaPlusPi0900GeV->Get("EfficiencyPi0");
TH1D* histoPi0AcceptPtPythiaPlus900GeV = (TH1D*)directoryPythiaPlusPi0900GeV->Get("AcceptancePi0");
TH1D* histoMCPi0YieldPtPythiaPlus900GeV = (TH1D*)directoryPythiaPlusPi0900GeV->Get("Pi0_Input_Reweighted");
TH1D* histoMCPi0YieldPtPythiaPlus900GeVWOWeights = (TH1D*)directoryPythiaPlusPi0900GeV->Get("Pi0_Input");
TH1D* histoPi0WeightsPythiaPlus900GeV = (TH1D*)directoryPythiaPlusPi0900GeV->Get("Pi0_Weights");
TH1D* histoMCPi0YieldPtPythiaPlus900GeVAddedSig = (TH1D*)directoryPythiaPlusPi0900GeV->Get("Pi0_Input_Reweighted_AddedSig");
TH1D* histoMCPi0YieldPtPythiaPlus900GeVAddedSigWOWeights = (TH1D*)directoryPythiaPlusPi0900GeV->Get("Pi0_Input_AddedSig");
TH1D* histoPi0RawYieldPythiaPlus900GeV = (TH1D*)directoryPythiaPlusPi0900GeV->Get("RAWYieldPerEventsPi0");
histoPi0MassDataPythiaPlus900GeV->Scale(1000.);
histoPi0MassMCPythiaPlus900GeV->Scale(1000.);
TDirectory* directoryPythiaPlusEta900GeV = (TDirectory*)filePythiaPlus->Get("Eta_pp_900GeV");
TH1D* histoEtaCorrectedSpecPythiaPlus900GeV = (TH1D*)directoryPythiaPlusEta900GeV->Get("CorrectedYieldEta");
TGraphAsymmErrors* graphEtaCorrectedSpecSysPythiaPlus900GeV = (TGraphAsymmErrors*)directoryPythiaPlusEta900GeV->Get("EtaSystError");
TH1D* histoEtaMassDataPythiaPlus900GeV = (TH1D*)directoryPythiaPlusEta900GeV->Get("MassEta");
TH1D* histoEtaMassMCPythiaPlus900GeV = (TH1D*)directoryPythiaPlusEta900GeV->Get("TrueMassEta");
TH1D* histoEtaWidthDataPythiaPlus900GeV = (TH1D*)directoryPythiaPlusEta900GeV->Get("FWHMEtaMeV");
TH1D* histoEtaWidthMCPythiaPlus900GeV = (TH1D*)directoryPythiaPlusEta900GeV->Get("TrueFWHMEtaMeV");
TH1D* histoEtaTrueEffiPtPythiaPlus900GeV = (TH1D*)directoryPythiaPlusEta900GeV->Get("EfficiencyEta");
TH1D* histoEtaAcceptPtPythiaPlus900GeV = (TH1D*)directoryPythiaPlusEta900GeV->Get("EfficiencyEta");
TH1D* histoMCEtaYieldPtPythiaPlus900GeV = (TH1D*)directoryPythiaPlusEta900GeV->Get("Eta_Input_Reweighted");
TH1D* histoMCEtaYieldPtPythiaPlus900GeVWOWeights = (TH1D*)directoryPythiaPlusEta900GeV->Get("Eta_Input");
TH1D* histoEtaWeightsPythiaPlus900GeV = (TH1D*)directoryPythiaPlusEta900GeV->Get("Eta_Weights");
TH1D* histoMCEtaYieldPtPythiaPlus900GeVAddedSig = (TH1D*)directoryPythiaPlusEta900GeV->Get("Eta_Input_Reweighted_AddedSig");
TH1D* histoMCEtaYieldPtPythiaPlus900GeVAddedSigWOWeights = (TH1D*)directoryPythiaPlusEta900GeV->Get("Eta_Input_AddedSig");
TH1D* histoEtaRawYieldPythiaPlus900GeV = (TH1D*)directoryPythiaPlusEta900GeV->Get("RAWYieldPerEventsEta");
histoEtaMassDataPythiaPlus900GeV->Scale(1000.);
histoEtaMassMCPythiaPlus900GeV->Scale(1000.);
//***************************************************************************
//************************** PHOJET *************************************
TFile* filePhojet = new TFile(nameFilepp2.Data());
//********* 8 TeV *****************************************************************************************************
TDirectory* directoryPhojetPi08TeV = (TDirectory*)filePhojet->Get("Pi0_pp_8TeV");
TH1D* histoPi0CorrectedSpecPhojet8TeV = (TH1D*)directoryPhojetPi08TeV->Get("CorrectedYieldPi0");
TGraphAsymmErrors* graphPi0CorrectedSpecSysPhojet8TeV = (TGraphAsymmErrors*)directoryPhojetPi08TeV->Get("Pi0SystError");
TH1D* histoPi0MassDataPhojet8TeV = (TH1D*)directoryPhojetPi08TeV->Get("MassPi0");
TH1D* histoPi0MassMCPhojet8TeV = (TH1D*)directoryPhojetPi08TeV->Get("TrueMassPi0");
TH1D* histoPi0WidthDataPhojet8TeV = (TH1D*)directoryPhojetPi08TeV->Get("FWHMPi0MeV");
TH1D* histoPi0WidthMCPhojet8TeV = (TH1D*)directoryPhojetPi08TeV->Get("TrueFWHMPi0MeV");
TH1D* histoPi0TrueEffiPtPhojet8TeV = (TH1D*)directoryPhojetPi08TeV->Get("EfficiencyPi0");
TH1D* histoPi0AcceptPtPhojet8TeV = (TH1D*)directoryPhojetPi08TeV->Get("AcceptancePi0");
TH1D* histoMCPi0YieldPtPhojet8TeV = (TH1D*)directoryPhojetPi08TeV->Get("Pi0_Input_Reweighted");
TH1D* histoMCPi0YieldPtPhojet8TeVWOWeights = (TH1D*)directoryPhojetPi08TeV->Get("Pi0_Input");
TH1D* histoPi0WeightsPhojet8TeV = (TH1D*)directoryPhojetPi08TeV->Get("Pi0_Weights");
TH1D* histoPi0RawYieldPhojet8TeV = (TH1D*)directoryPhojetPi08TeV->Get("RAWYieldPerEventsPi0");
histoPi0MassDataPhojet8TeV->Scale(1000.);
histoPi0MassMCPhojet8TeV->Scale(1000.);
TDirectory* directoryPhojetEta8TeV = (TDirectory*)filePhojet->Get("Eta_pp_8TeV");
TH1D* histoEtaCorrectedSpecPhojet8TeV = (TH1D*)directoryPhojetEta8TeV->Get("CorrectedYieldEta");
TGraphAsymmErrors* graphEtaCorrectedSpecSysPhojet8TeV = (TGraphAsymmErrors*)directoryPhojetEta8TeV->Get("EtaSystError");
TH1D* histoEtaMassDataPhojet8TeV = (TH1D*)directoryPhojetEta8TeV->Get("MassEta");
TH1D* histoEtaMassMCPhojet8TeV = (TH1D*)directoryPhojetEta8TeV->Get("TrueMassEta");
TH1D* histoEtaWidthDataPhojet8TeV = (TH1D*)directoryPhojetEta8TeV->Get("FWHMEtaMeV");
TH1D* histoEtaWidthMCPhojet8TeV = (TH1D*)directoryPhojetEta8TeV->Get("TrueFWHMEtaMeV");
TH1D* histoEtaTrueEffiPtPhojet8TeV = (TH1D*)directoryPhojetEta8TeV->Get("EfficiencyEta");
TH1D* histoEtaAcceptPtPhojet8TeV = (TH1D*)directoryPhojetEta8TeV->Get("EfficiencyEta");
TH1D* histoMCEtaYieldPtPhojet8TeV = (TH1D*)directoryPhojetEta8TeV->Get("Eta_Input_Reweighted");
TH1D* histoMCEtaYieldPtPhojet8TeVWOWeights = (TH1D*)directoryPhojetEta8TeV->Get("Eta_Input");
TH1D* histoEtaWeightsPhojet8TeV = (TH1D*)directoryPhojetEta8TeV->Get("Eta_Weights");
TH1D* histoEtaRawYieldPhojet8TeV = (TH1D*)directoryPhojetEta8TeV->Get("RAWYieldPerEventsEta");
histoEtaMassDataPhojet8TeV->Scale(1000.);
histoEtaMassMCPhojet8TeV->Scale(1000.);
//********* 7 TeV *****************************************************************************************************
TDirectory* directoryPhojetPi07TeV = (TDirectory*)filePhojet->Get("Pi0_pp_7TeV");
TH1D* histoPi0CorrectedSpecPhojet7TeV = (TH1D*)directoryPhojetPi07TeV->Get("CorrectedYieldPi0");
TGraphAsymmErrors* graphPi0CorrectedSpecSysPhojet7TeV = (TGraphAsymmErrors*)directoryPhojetPi07TeV->Get("Pi0SystError");
TH1D* histoPi0MassDataPhojet7TeV = (TH1D*)directoryPhojetPi07TeV->Get("MassPi0");
TH1D* histoPi0MassMCPhojet7TeV = (TH1D*)directoryPhojetPi07TeV->Get("TrueMassPi0");
TH1D* histoPi0WidthDataPhojet7TeV = (TH1D*)directoryPhojetPi07TeV->Get("FWHMPi0MeV");
TH1D* histoPi0WidthMCPhojet7TeV = (TH1D*)directoryPhojetPi07TeV->Get("TrueFWHMPi0MeV");
TH1D* histoPi0TrueEffiPtPhojet7TeV = (TH1D*)directoryPhojetPi07TeV->Get("EfficiencyPi0");
TH1D* histoPi0AcceptPtPhojet7TeV = (TH1D*)directoryPhojetPi07TeV->Get("AcceptancePi0");
TH1D* histoMCPi0YieldPtPhojet7TeV = (TH1D*)directoryPhojetPi07TeV->Get("Pi0_Input_Reweighted");
TH1D* histoMCPi0YieldPtPhojet7TeVWOWeights = (TH1D*)directoryPhojetPi07TeV->Get("Pi0_Input");
TH1D* histoPi0WeightsPhojet7TeV = (TH1D*)directoryPhojetPi07TeV->Get("Pi0_Weights");
TH1D* histoPi0RawYieldPhojet7TeV = (TH1D*)directoryPhojetPi07TeV->Get("RAWYieldPerEventsPi0");
histoPi0MassDataPhojet7TeV->Scale(1000.);
histoPi0MassMCPhojet7TeV->Scale(1000.);
TDirectory* directoryPhojetEta7TeV = (TDirectory*)filePhojet->Get("Eta_pp_7TeV");
TH1D* histoEtaCorrectedSpecPhojet7TeV = (TH1D*)directoryPhojetEta7TeV->Get("CorrectedYieldEta");
TGraphAsymmErrors* graphEtaCorrectedSpecSysPhojet7TeV = (TGraphAsymmErrors*)directoryPhojetEta7TeV->Get("EtaSystError");
TH1D* histoEtaMassDataPhojet7TeV = (TH1D*)directoryPhojetEta7TeV->Get("MassEta");
TH1D* histoEtaMassMCPhojet7TeV = (TH1D*)directoryPhojetEta7TeV->Get("TrueMassEta");
TH1D* histoEtaWidthDataPhojet7TeV = (TH1D*)directoryPhojetEta7TeV->Get("FWHMEtaMeV");
TH1D* histoEtaWidthMCPhojet7TeV = (TH1D*)directoryPhojetEta7TeV->Get("TrueFWHMEtaMeV");
TH1D* histoEtaTrueEffiPtPhojet7TeV = (TH1D*)directoryPhojetEta7TeV->Get("EfficiencyEta");
TH1D* histoEtaAcceptPtPhojet7TeV = (TH1D*)directoryPhojetEta7TeV->Get("EfficiencyEta");
TH1D* histoMCEtaYieldPtPhojet7TeV = (TH1D*)directoryPhojetEta7TeV->Get("Eta_Input_Reweighted");
TH1D* histoMCEtaYieldPtPhojet7TeVWOWeights = (TH1D*)directoryPhojetEta7TeV->Get("Eta_Input");
TH1D* histoEtaWeightsPhojet7TeV = (TH1D*)directoryPhojetEta7TeV->Get("Eta_Weights");
TH1D* histoEtaRawYieldPhojet7TeV = (TH1D*)directoryPhojetEta7TeV->Get("RAWYieldPerEventsEta");
histoEtaMassDataPhojet7TeV->Scale(1000.);
histoEtaMassMCPhojet7TeV->Scale(1000.);
//********* 2.76 TeV *****************************************************************************************************
TDirectory* directoryPhojetPi02760GeV = (TDirectory*)filePhojet->Get("Pi0_pp_2760GeV");
TH1D* histoPi0CorrectedSpecPhojet2760GeV = (TH1D*)directoryPhojetPi02760GeV->Get("CorrectedYieldPi0");
TGraphAsymmErrors* graphPi0CorrectedSpecSysPhojet2760GeV = (TGraphAsymmErrors*)directoryPhojetPi02760GeV->Get("Pi0SystError");
TH1D* histoPi0MassDataPhojet2760GeV = (TH1D*)directoryPhojetPi02760GeV->Get("MassPi0");
TH1D* histoPi0MassMCPhojet2760GeV = (TH1D*)directoryPhojetPi02760GeV->Get("TrueMassPi0");
TH1D* histoPi0WidthDataPhojet2760GeV = (TH1D*)directoryPhojetPi02760GeV->Get("FWHMPi0MeV");
TH1D* histoPi0WidthMCPhojet2760GeV = (TH1D*)directoryPhojetPi02760GeV->Get("TrueFWHMPi0MeV");
TH1D* histoPi0TrueEffiPtPhojet2760GeV = (TH1D*)directoryPhojetPi02760GeV->Get("EfficiencyPi0");
TH1D* histoPi0AcceptPtPhojet2760GeV = (TH1D*)directoryPhojetPi02760GeV->Get("AcceptancePi0");
TH1D* histoMCPi0YieldPtPhojet2760GeV = (TH1D*)directoryPhojetPi02760GeV->Get("Pi0_Input_Reweighted");
TH1D* histoMCPi0YieldPtPhojet2760GeVWOWeights = (TH1D*)directoryPhojetPi02760GeV->Get("Pi0_Input");
TH1D* histoPi0WeightsPhojet2760GeV = (TH1D*)directoryPhojetPi02760GeV->Get("Pi0_Weights");
TH1D* histoPi0RawYieldPhojet2760GeV = (TH1D*)directoryPhojetPi02760GeV->Get("RAWYieldPerEventsPi0");
histoPi0MassDataPhojet2760GeV->Scale(1000.);
histoPi0MassMCPhojet2760GeV->Scale(1000.);
TDirectory* directoryPhojetEta2760GeV = (TDirectory*)filePhojet->Get("Eta_pp_2760GeV");
TH1D* histoEtaCorrectedSpecPhojet2760GeV = (TH1D*)directoryPhojetEta2760GeV->Get("CorrectedYieldEta");
TGraphAsymmErrors* graphEtaCorrectedSpecSysPhojet2760GeV = (TGraphAsymmErrors*)directoryPhojetEta2760GeV->Get("EtaSystError");
TH1D* histoEtaMassDataPhojet2760GeV = (TH1D*)directoryPhojetEta2760GeV->Get("MassEta");
TH1D* histoEtaMassMCPhojet2760GeV = (TH1D*)directoryPhojetEta2760GeV->Get("TrueMassEta");
TH1D* histoEtaWidthDataPhojet2760GeV = (TH1D*)directoryPhojetEta2760GeV->Get("FWHMEtaMeV");
TH1D* histoEtaWidthMCPhojet2760GeV = (TH1D*)directoryPhojetEta2760GeV->Get("TrueFWHMEtaMeV");
TH1D* histoEtaTrueEffiPtPhojet2760GeV = (TH1D*)directoryPhojetEta2760GeV->Get("EfficiencyEta");
TH1D* histoEtaAcceptPtPhojet2760GeV = (TH1D*)directoryPhojetEta2760GeV->Get("EfficiencyEta");
TH1D* histoMCEtaYieldPtPhojet2760GeV = (TH1D*)directoryPhojetEta2760GeV->Get("Eta_Input_Reweighted");
TH1D* histoMCEtaYieldPtPhojet2760GeVWOWeights = (TH1D*)directoryPhojetEta2760GeV->Get("Eta_Input");
TH1D* histoEtaWeightsPhojet2760GeV = (TH1D*)directoryPhojetEta2760GeV->Get("Eta_Weights");
TH1D* histoEtaRawYieldPhojet2760GeV = (TH1D*)directoryPhojetEta2760GeV->Get("RAWYieldPerEventsEta");
histoEtaMassDataPhojet2760GeV->Scale(1000.);
histoEtaMassMCPhojet2760GeV->Scale(1000.);
//********* 0.9 TeV *****************************************************************************************************
TDirectory* directoryPhojetPi0900GeV = (TDirectory*)filePhojet->Get("Pi0_pp_900GeV");
TH1D* histoPi0CorrectedSpecPhojet900GeV = (TH1D*)directoryPhojetPi0900GeV->Get("CorrectedYieldPi0");
TGraphAsymmErrors* graphPi0CorrectedSpecSysPhojet900GeV = (TGraphAsymmErrors*)directoryPhojetPi0900GeV->Get("Pi0SystError");
TH1D* histoPi0MassDataPhojet900GeV = (TH1D*)directoryPhojetPi0900GeV->Get("MassPi0");
TH1D* histoPi0MassMCPhojet900GeV = (TH1D*)directoryPhojetPi0900GeV->Get("TrueMassPi0");
TH1D* histoPi0WidthDataPhojet900GeV = (TH1D*)directoryPhojetPi0900GeV->Get("FWHMPi0MeV");
TH1D* histoPi0WidthMCPhojet900GeV = (TH1D*)directoryPhojetPi0900GeV->Get("TrueFWHMPi0MeV");
TH1D* histoPi0TrueEffiPtPhojet900GeV = (TH1D*)directoryPhojetPi0900GeV->Get("EfficiencyPi0");
TH1D* histoPi0AcceptPtPhojet900GeV = (TH1D*)directoryPhojetPi0900GeV->Get("AcceptancePi0");
TH1D* histoMCPi0YieldPtPhojet900GeV = (TH1D*)directoryPhojetPi0900GeV->Get("Pi0_Input_Reweighted");
TH1D* histoMCPi0YieldPtPhojet900GeVWOWeights = (TH1D*)directoryPhojetPi0900GeV->Get("Pi0_Input");
TH1D* histoPi0WeightsPhojet900GeV = (TH1D*)directoryPhojetPi0900GeV->Get("Pi0_Weights");
TH1D* histoPi0RawYieldPhojet900GeV = (TH1D*)directoryPhojetPi0900GeV->Get("RAWYieldPerEventsPi0");
histoPi0MassDataPhojet900GeV->Scale(1000.);
histoPi0MassMCPhojet900GeV->Scale(1000.);
TDirectory* directoryPhojetEta900GeV = (TDirectory*)filePhojet->Get("Eta_pp_900GeV");
TH1D* histoEtaCorrectedSpecPhojet900GeV = (TH1D*)directoryPhojetEta900GeV->Get("CorrectedYieldEta");
TGraphAsymmErrors* graphEtaCorrectedSpecSysPhojet900GeV = (TGraphAsymmErrors*)directoryPhojetEta900GeV->Get("EtaSystError");
TH1D* histoEtaMassDataPhojet900GeV = (TH1D*)directoryPhojetEta900GeV->Get("MassEta");
TH1D* histoEtaMassMCPhojet900GeV = (TH1D*)directoryPhojetEta900GeV->Get("TrueMassEta");
TH1D* histoEtaWidthDataPhojet900GeV = (TH1D*)directoryPhojetEta900GeV->Get("FWHMEtaMeV");
TH1D* histoEtaWidthMCPhojet900GeV = (TH1D*)directoryPhojetEta900GeV->Get("TrueFWHMEtaMeV");
TH1D* histoEtaTrueEffiPtPhojet900GeV = (TH1D*)directoryPhojetEta900GeV->Get("EfficiencyEta");
TH1D* histoEtaAcceptPtPhojet900GeV = (TH1D*)directoryPhojetEta900GeV->Get("EfficiencyEta");
TH1D* histoMCEtaYieldPtPhojet900GeV = (TH1D*)directoryPhojetEta900GeV->Get("Eta_Input_Reweighted");
TH1D* histoMCEtaYieldPtPhojet900GeVWOWeights = (TH1D*)directoryPhojetEta900GeV->Get("Eta_Input");
TH1D* histoEtaWeightsPhojet900GeV = (TH1D*)directoryPhojetEta900GeV->Get("Eta_Weights");
TH1D* histoEtaRawYieldPhojet900GeV = (TH1D*)directoryPhojetEta900GeV->Get("RAWYieldPerEventsEta");
histoEtaMassDataPhojet900GeV->Scale(1000.);
histoEtaMassMCPhojet900GeV->Scale(1000.);
TFile* fileNeutralPionDataPP = new TFile(nameFilePPOld.Data());
//************* 7 TeV *****************
TDirectory* directoryPi07TeV = (TDirectory*)fileNeutralPionDataPP->Get("Pi07TeV");
TH1D* histoAccPi07TeV = (TH1D*)directoryPi07TeV->Get("AcceptancePi0");
TH1D* histoTrueEffPtPi07TeV = (TH1D*)directoryPi07TeV->Get("EfficiencyPi0");
TH1D* histoRawYieldPi07TeV = (TH1D*)directoryPi07TeV->Get("RAWYieldPerEventsPi0");
TDirectory* directoryEta7TeV = (TDirectory*)fileNeutralPionDataPP->Get("Eta7TeV");
TH1D* histoCorrectedYieldEta7TeV = (TH1D*)directoryEta7TeV->Get("CorrectedYieldEta");
TGraphAsymmErrors* graphCorrectedYieldSysEta7TeV = (TGraphAsymmErrors*)directoryEta7TeV->Get("EtaSystError");
TH1D* histoAccEta7TeV = (TH1D*)directoryEta7TeV->Get("AcceptanceEta");
TH1D* histoTrueEffPtEta7TeV = (TH1D*)directoryEta7TeV->Get("EfficiencyEta");
TH1D* histoRawYieldEta7TeV = (TH1D*)directoryEta7TeV->Get("RAWYieldPerEventsEta");
TH1D* histoEtaMassData7TeV = (TH1D*)directoryEta7TeV->Get("MassEta");
TH1D* histoEtaWidthData7TeV = (TH1D*)directoryEta7TeV->Get("FWHMEtaMeV");
TH1D* histoEtaMassMC7TeV = (TH1D*)directoryEta7TeV->Get("TrueMassEta");
TH1D* histoEtaWidthMC7TeV = (TH1D*)directoryEta7TeV->Get("TrueFWHMEtaMeV");
TH1D* histoRatioEtaPi07TeV= (TH1D*)directoryEta7TeV->Get("EtatoPi0RatioConversion");
TGraphAsymmErrors* graphRatioEtaPi0SystErr7TeV= (TGraphAsymmErrors*)directoryEta7TeV->Get("EtatoPi0RatioConversionSys");
histoEtaMassData7TeV->Scale(1000.);
histoEtaMassMC7TeV->Scale(1000.);
histoEtaMassData7TeV->SetBinContent(histoEtaMassData7TeV->GetNbinsX(),0);
histoEtaMassMC7TeV->SetBinContent(histoEtaMassMC7TeV->GetNbinsX(),0);
histoEtaWidthData7TeV->SetBinContent(histoEtaWidthData7TeV->GetNbinsX(),10000.);
histoEtaWidthMC7TeV->SetBinContent(histoEtaWidthMC7TeV->GetNbinsX(),10000.);
//************* 0.9 TeV *****************
TDirectory* directoryPi0900GeV = (TDirectory*)fileNeutralPionDataPP->Get("Pi0900GeV");
TH1D* histoAccPi0900GeV = (TH1D*)directoryPi0900GeV->Get("AcceptancePi0");
TH1D* histoTrueEffPtPi0900GeV = (TH1D*)directoryPi0900GeV->Get("EfficiencyPi0");
TH1D* histoRawYieldPi0900GeV = (TH1D*)directoryPi0900GeV->Get("RAWYieldPerEventsPi0");
TDirectory* directoryEta900GeV = (TDirectory*)fileNeutralPionDataPP->Get("Eta900GeV");
TH1D* histoCorrectedYieldEta900GeV = (TH1D*)directoryEta900GeV->Get("CorrectedYieldEta");
TGraphAsymmErrors* graphCorrectedYieldSysEta900GeV = (TGraphAsymmErrors*)directoryEta900GeV->Get("EtaSystError");
TH1D* histoAccEta900GeV = (TH1D*)directoryEta900GeV->Get("AcceptanceEta");
TH1D* histoTrueEffPtEta900GeV = (TH1D*)directoryEta900GeV->Get("EfficiencyEta");
TH1D* histoRawYieldEta900GeV = (TH1D*)directoryEta900GeV->Get("RAWYieldPerEventsEta");
TH1D* histoRatioEtaPi0900GeV= (TH1D*)directoryEta900GeV->Get("EtatoPi0RatioConversion");
//************* 2.76 TeV *****************
TGraphAsymmErrors* graphRatioEtaPi0SystErr900GeV= (TGraphAsymmErrors*)directoryEta900GeV->Get("EtatoPi0RatioConversionSys");
TDirectory* directoryPi02760GeV = (TDirectory*)fileNeutralPionDataPP->Get("Pi02.76TeV");
TH1D* histoAccPi02760GeV = (TH1D*)directoryPi02760GeV->Get("AcceptancePi0");
TH1D* histoTrueEffPtPi02760GeV = (TH1D*)directoryPi02760GeV->Get("EfficiencyPi0");
TH1D* histoRawYieldPi02760GeV = (TH1D*)directoryPi02760GeV->Get("RAWYieldPerEventsPi0");
TH1D* histoCorrectedYieldPi02760GeV = (TH1D*)directoryPi02760GeV->Get("CorrectedYieldPi0");
TGraphAsymmErrors* graphCorrectedYieldSysPi02760GeV = (TGraphAsymmErrors*)directoryPi02760GeV->Get("Pi0SystError");
TDirectory* directoryEta2760GeV = (TDirectory*)fileNeutralPionDataPP->Get("Eta2.76TeV");
TH1D* histoAccEta2760GeV = (TH1D*)directoryEta2760GeV->Get("AcceptanceEta");
TH1D* histoTrueEffPtEta2760GeV = (TH1D*)directoryEta2760GeV->Get("EfficiencyEta");
TH1D* histoCorrectedYieldEta2760GeV = (TH1D*)directoryEta2760GeV->Get("CorrectedYieldEta");
TGraphAsymmErrors* graphCorrectedYieldSysEta2760GeV = (TGraphAsymmErrors*)directoryEta2760GeV->Get("EtaSystError");
TH1D* histoRawYieldEta2760GeV = (TH1D*)directoryEta2760GeV->Get("RAWYieldPerEventsEta");
TH1D* histoEtaMassData2760GeV = (TH1D*)directoryEta2760GeV->Get("MassEta");
TH1D* histoEtaWidthData2760GeV = (TH1D*)directoryEta2760GeV->Get("FWHMEtaMeV");
TH1D* histoEtaMassMC2760GeV = (TH1D*)directoryEta2760GeV->Get("TrueMassEta");
TH1D* histoEtaWidthMC2760GeV = (TH1D*)directoryEta2760GeV->Get("TrueFWHMEtaMeV");
TH1D* histoRatioEtaPi02760GeV= (TH1D*)directoryEta2760GeV->Get("EtatoPi0RatioConversion");
TGraphAsymmErrors* graphRatioEtaPi0SystErr2760GeV= (TGraphAsymmErrors*)directoryEta2760GeV->Get("EtatoPi0RatioConversionSys");
histoEtaMassData2760GeV->Scale(1000.);
histoEtaMassMC2760GeV->Scale(1000.);
histoEtaMassData2760GeV->SetBinContent(histoEtaMassData2760GeV->GetNbinsX(),0);
histoEtaMassMC2760GeV->SetBinContent(histoEtaMassMC2760GeV->GetNbinsX(),0);
histoEtaWidthData2760GeV->SetBinContent(histoEtaWidthData2760GeV->GetNbinsX(),10000.);
histoEtaWidthMC2760GeV->SetBinContent(histoEtaWidthMC2760GeV->GetNbinsX(),10000.);
TH1D* histoPCMMassDataPP2760GeV = (TH1D*)directoryPi02760GeV->Get("MassPi0");
TH1D* histoPCMWidthDataPP2760GeV = (TH1D*)directoryPi02760GeV->Get("FWHMPi0MeV");
TH1D* histoPCMMassMCPP2760GeV = (TH1D*)directoryPi02760GeV->Get("TrueMassPi0");
TH1D* histoPCMWidthMCPP2760GeV = (TH1D*)directoryPi02760GeV->Get("TrueFWHMPi0MeV");
histoPCMMassDataPP2760GeV->Scale(1000.);
histoPCMMassMCPP2760GeV->Scale(1000.);
histoPCMMassDataPP2760GeV->SetBinContent(histoPCMMassDataPP2760GeV->GetNbinsX(),0);
histoPCMMassMCPP2760GeV->SetBinContent(histoPCMMassMCPP2760GeV->GetNbinsX(),0);
histoPCMWidthDataPP2760GeV->SetBinContent(histoPCMWidthDataPP2760GeV->GetNbinsX(),10000.);
histoPCMWidthMCPP2760GeV->SetBinContent(histoPCMWidthMCPP2760GeV->GetNbinsX(),10000.);
cout << "here" << endl;
TFile* filePi0Pythia = new TFile("0000011_002092570028250400000_01521065000000/8TeV/Pi0_data_GammaConvV1Correction_0000011_002092570028250400000_01521065000000.root");
TH1D* histoMCYieldPi0Pythia = (TH1D*)filePi0Pythia->Get("MCYield_Meson_oldBin");
TFile* filePi0PythiaPlus = new TFile("0000012_002092570028250400000_01521065000000/8TeV/Pi0_MC_GammaConvV1Correction_0000012_002092570028250400000_01521065000000.root");
TH1D* histoMCYieldPi0PythiaPlus = (TH1D*)filePi0PythiaPlus->Get("MCYield_Meson_oldBin");
TFile* fileEtaPythia = new TFile("0000011_002092570028250400000_01521065000000/8TeV/Eta_data_GammaConvV1Correction_0000011_002092570028250400000_01521065000000.root");
TH1D* histoMCYieldEtaPythia = (TH1D*)fileEtaPythia->Get("MCYield_Meson_oldBin");
TFile* fileEtaPythiaPlus = new TFile("0000012_002092570028250400000_01521065000000/8TeV/Eta_MC_GammaConvV1Correction_0000012_002092570028250400000_01521065000000.root");
TH1D* histoMCYieldEtaPythiaPlus = (TH1D*)fileEtaPythiaPlus->Get("MCYield_Meson_oldBin");
Width_t widthLinesBoxes = 1.4;
Width_t widthCommonFit = 2.;
Width_t widthStatErrBars = 1.5;
Width_t widthCommonErrors = 1.1;
Width_t widthCommonSpectrumBoxes = 0.99;
if (suffix.CompareTo("eps")==0){
widthLinesBoxes = 1.4;
widthCommonFit = 2.;
widthStatErrBars = 1.5;
widthCommonErrors = 1.1;
widthCommonSpectrumBoxes = 0.99;
} else {
widthLinesBoxes = 2.3;
widthCommonFit = 2.6;
widthStatErrBars = 2.6;
widthCommonErrors = 2.;
widthCommonSpectrumBoxes = 2.3;
}
cout << "Pi0 efficiency in different system - Pythia" << endl;
TCanvas* canvasEfficiencyPi0Pythia = new TCanvas("canvasEfficiencyPi0Pythia","",200,10,1350,1350); // gives the page size
DrawGammaCanvasSettings(canvasEfficiencyPi0Pythia, 0.1, 0.02, 0.035, 0.09);
TH2F * histo2DEffPi0Pythia;
histo2DEffPi0Pythia = new TH2F("histo2DEffPi0Pythia","histo2DEffPi0Pythia",1000,0,8,2000,0.e-3,3.7e-3 );
SetStyleHistoTH2ForGraphs(histo2DEffPi0Pythia, "#it{p}_{T} (GeV/#it{c})","#epsilon_{reco, #pi^{0}}",0.03,0.04, 0.03,0.04, 1.,1.);
histo2DEffPi0Pythia->Draw("copy");
DrawGammaSetMarker(histoPi0TrueEffiPtPythia8TeV, markerStyleSpectrum8TeVMC, markerSizePi0PP8TeV, colorPi08TeV, colorPi08TeV);
histoPi0TrueEffiPtPythia8TeV->DrawCopy("e1,same");
DrawGammaSetMarker(histoPi0TrueEffiPtPythia7TeV, markerStyleSpectrum7TeVMC, markerSizePi0PP7TeV, colorPi07TeV, colorPi07TeV);
histoPi0TrueEffiPtPythia7TeV->DrawCopy("e1,same");
DrawGammaSetMarker(histoPi0TrueEffiPtPythia2760GeV, markerStyleSpectrum2760GeVMC, markerSizePi0PP2760GeV, colorPi02760GeV, colorPi02760GeV);
histoPi0TrueEffiPtPythia2760GeV->DrawCopy("pe1,same");
DrawGammaSetMarker(histoPi0TrueEffiPtPythia900GeV, markerStyleSpectrum900GeVMC, markerSizePi0PP900GeV, colorPi0900GeV, colorPi0900GeV);
histoPi0TrueEffiPtPythia900GeV->DrawCopy("e1,same");
TLegend* legendEffiPi0Pythia = new TLegend(0.12,0.78,0.83,0.93);
legendEffiPi0Pythia->SetFillColor(0);
legendEffiPi0Pythia->SetLineColor(0);
legendEffiPi0Pythia->SetTextSize(0.027);
legendEffiPi0Pythia->SetNColumns(2);
legendEffiPi0Pythia->AddEntry(histoPi0TrueEffiPtPythia8TeV,"pp #sqrt{#it{s}} = 8 TeV","p");
legendEffiPi0Pythia->AddEntry(histoPi0TrueEffiPtPythia7TeV,"pp #sqrt{#it{s}} = 7 TeV","p");
legendEffiPi0Pythia->AddEntry(histoPi0TrueEffiPtPythia2760GeV,"pp #sqrt{#it{s}} = 2.76 TeV","p");
legendEffiPi0Pythia->AddEntry(histoPi0TrueEffiPtPythia900GeV,"pp #sqrt{#it{s}} = 0.9 TeV","p");
legendEffiPi0Pythia->Draw();
canvasEfficiencyPi0Pythia->SaveAs(Form("%s/EfficiencyPi0Pythia.%s",outputDir.Data(),suffix.Data()));
histo2DEffPi0Pythia->Draw("copy");
cout << "Pi0 efficiency in different system - Pythia with added signals" << endl;
TCanvas* canvasEfficiencyPi0PythiaPlus = new TCanvas("canvasEfficiencyPi0PythiaPlus","",200,10,1350,1350); // gives the page size
DrawGammaCanvasSettings(canvasEfficiencyPi0PythiaPlus, 0.1, 0.02, 0.035, 0.09);
TH2F * histo2DEffPi0PythiaPlus;
histo2DEffPi0PythiaPlus = new TH2F("histo2DEffPi0PythiaPlus","histo2DEffPi0PythiaPlus",1000,0,8,2000,0.e-3,3.7e-3 );
SetStyleHistoTH2ForGraphs(histo2DEffPi0PythiaPlus, "#it{p}_{T} (GeV/#it{c})","#epsilon_{reco, #pi^{0}}",0.03,0.04, 0.03,0.04, 1.,1.);
histo2DEffPi0PythiaPlus->Draw("copy");
// DrawGammaSetMarker(histoTrueEffPtPi07TeV, markerStyleSpectrum7TeVMC, markerSizePi0PP7TeV, colorPi07TeV, colorPi07TeV);
// histoTrueEffPtPi07TeV->DrawCopy("e1,same");
// DrawGammaSetMarker(histoTrueEffPtPi02760GeV, markerStyleSpectrum2760GeVMC, markerSizePi0PP2760GeV, colorPi02760GeV, colorPi02760GeV);
// histoTrueEffPtPi02760GeV->DrawCopy("pe1,same");
// DrawGammaSetMarker(histoTrueEffPtPi0900GeV, markerStyleSpectrum900GeVMC, markerSizePi0PP900GeV, colorPi0900GeV, colorPi0900GeV);
// histoTrueEffPtPi0900GeV->DrawCopy("e1,same");
DrawGammaSetMarker(histoPi0TrueEffiPtPythiaPlus8TeV, markerStyleSpectrum8TeVMC, markerSizePi0PP8TeV, colorPi08TeV, colorPi08TeV);
histoPi0TrueEffiPtPythiaPlus8TeV->DrawCopy("e1,same");
DrawGammaSetMarker(histoPi0TrueEffiPtPythiaPlus7TeV, markerStyleSpectrum7TeVMC, markerSizePi0PP7TeV, colorPi07TeV, colorPi07TeV);
histoPi0TrueEffiPtPythiaPlus7TeV->DrawCopy("e1,same");
DrawGammaSetMarker(histoPi0TrueEffiPtPythiaPlus2760GeV, markerStyleSpectrum2760GeVMC, markerSizePi0PP2760GeV, colorPi02760GeV, colorPi02760GeV);
histoPi0TrueEffiPtPythiaPlus2760GeV->DrawCopy("pe1,same");
DrawGammaSetMarker(histoPi0TrueEffiPtPythiaPlus900GeV, markerStyleSpectrum900GeVMC, markerSizePi0PP900GeV, colorPi0900GeV, colorPi0900GeV);
histoPi0TrueEffiPtPythiaPlus900GeV->DrawCopy("e1,same");
TLegend* legendEffiPi0PythiaPlus = new TLegend(0.12,0.78,0.83,0.93);
legendEffiPi0PythiaPlus->SetFillColor(0);
legendEffiPi0PythiaPlus->SetLineColor(0);
legendEffiPi0PythiaPlus->SetTextSize(0.027);
legendEffiPi0PythiaPlus->SetNColumns(2);
legendEffiPi0PythiaPlus->AddEntry(histoPi0TrueEffiPtPythiaPlus8TeV,"pp #sqrt{#it{s}} = 8 TeV","p");
legendEffiPi0PythiaPlus->AddEntry(histoPi0TrueEffiPtPythiaPlus7TeV,"pp #sqrt{#it{s}} = 7 TeV","p");
legendEffiPi0PythiaPlus->AddEntry(histoPi0TrueEffiPtPythiaPlus2760GeV,"pp #sqrt{#it{s}} = 2.76 TeV","p");
legendEffiPi0PythiaPlus->AddEntry(histoPi0TrueEffiPtPythiaPlus900GeV,"pp #sqrt{#it{s}} = 0.9 TeV","p");
//legendEffiPi0PythiaPlus->AddEntry(histoTrueEffPtPi07TeV,"pp #sqrt{#it{s}} = 7 TeV","p");
//legendEffiPi0PythiaPlus->AddEntry(histoTrueEffPtPi02760GeV,"pp #sqrt{#it{s}} = 2.76 TeV","p");
//legendEffiPi0PythiaPlus->AddEntry(histoTrueEffPtPi0900GeV,"pp #sqrt{#it{s}} = 0.9 TeV","p");
legendEffiPi0PythiaPlus->Draw();
canvasEfficiencyPi0PythiaPlus->SaveAs(Form("%s/EfficiencyPi0PythiaPlus.%s",outputDir.Data(),suffix.Data()));
histo2DEffPi0PythiaPlus->Draw("copy");
cout << "Pi0 efficiency in different system - Phojet" << endl;
TCanvas* canvasEfficiencyPi0Phojet = new TCanvas("canvasEfficiencyPi0Phojet","",200,10,1350,1350); // gives the page size
DrawGammaCanvasSettings(canvasEfficiencyPi0Phojet, 0.1, 0.02, 0.035, 0.09);
TH2F * histo2DEffPi0Phojet;
histo2DEffPi0Phojet = new TH2F("histo2DEffPi0Phojet","histo2DEffPi0Phojet",1000,0,8,2000,0.e-3,3.7e-3 );
SetStyleHistoTH2ForGraphs(histo2DEffPi0Phojet, "#it{p}_{T} (GeV/#it{c})","#epsilon_{reco, #pi^{0}}",0.03,0.04, 0.03,0.04, 1.,1.);
histo2DEffPi0Phojet->Draw("copy");
// DrawGammaSetMarker(histoTrueEffPtPi07TeV, markerStyleSpectrum7TeVMC, markerSizePi0PP7TeV, colorPi07TeV, colorPi07TeV);
// histoTrueEffPtPi07TeV->DrawCopy("e1,same");
// DrawGammaSetMarker(histoTrueEffPtPi02760GeV, markerStyleSpectrum2760GeVMC, markerSizePi0PP2760GeV, colorPi02760GeV, colorPi02760GeV);
// histoTrueEffPtPi02760GeV->DrawCopy("pe1,same");
// DrawGammaSetMarker(histoTrueEffPtPi0900GeV, markerStyleSpectrum900GeVMC, markerSizePi0PP900GeV, colorPi0900GeV, colorPi0900GeV);
// histoTrueEffPtPi0900GeV->DrawCopy("e1,same");
DrawGammaSetMarker(histoPi0TrueEffiPtPhojet8TeV, markerStyleSpectrum8TeVMC, markerSizePi0PP8TeV, colorPi08TeV, colorPi08TeV);
histoPi0TrueEffiPtPhojet8TeV->DrawCopy("e1,same");
DrawGammaSetMarker(histoPi0TrueEffiPtPhojet7TeV, markerStyleSpectrum7TeVMC, markerSizePi0PP7TeV, colorPi07TeV, colorPi07TeV);
histoPi0TrueEffiPtPhojet7TeV->DrawCopy("e1,same");
DrawGammaSetMarker(histoPi0TrueEffiPtPhojet2760GeV, markerStyleSpectrum2760GeVMC, markerSizePi0PP2760GeV, colorPi02760GeV, colorPi02760GeV);
histoPi0TrueEffiPtPhojet2760GeV->DrawCopy("pe1,same");
DrawGammaSetMarker(histoPi0TrueEffiPtPhojet900GeV, markerStyleSpectrum900GeVMC, markerSizePi0PP900GeV, colorPi0900GeV, colorPi0900GeV);
histoPi0TrueEffiPtPhojet900GeV->DrawCopy("e1,same");
TLegend* legendEffiPi0Phojet = new TLegend(0.12,0.78,0.83,0.93);
legendEffiPi0Phojet->SetFillColor(0);
legendEffiPi0Phojet->SetLineColor(0);
legendEffiPi0Phojet->SetTextSize(0.027);
legendEffiPi0Phojet->SetNColumns(2);
legendEffiPi0Phojet->AddEntry(histoPi0TrueEffiPtPhojet8TeV,"pp #sqrt{#it{s}} = 8 TeV","p");
legendEffiPi0Phojet->AddEntry(histoPi0TrueEffiPtPhojet7TeV,"pp #sqrt{#it{s}} = 7 TeV","p");
legendEffiPi0Phojet->AddEntry(histoPi0TrueEffiPtPhojet2760GeV,"pp #sqrt{#it{s}} = 2.76 TeV","p");
legendEffiPi0Phojet->AddEntry(histoPi0TrueEffiPtPhojet900GeV,"pp #sqrt{#it{s}} = 0.9 TeV","p");
//legendEffiPi0Phojet->AddEntry(histoTrueEffPtPi07TeV,"pp #sqrt{#it{s}} = 7 TeV","p");
//legendEffiPi0Phojet->AddEntry(histoTrueEffPtPi02760GeV,"pp #sqrt{#it{s}} = 2.76 TeV","p");
//legendEffiPi0Phojet->AddEntry(histoTrueEffPtPi0900GeV,"pp #sqrt{#it{s}} = 0.9 TeV","p");
legendEffiPi0Phojet->Draw();
canvasEfficiencyPi0Phojet->SaveAs(Form("%s/EfficiencyPi0Phojet.%s",outputDir.Data(),suffix.Data()));
histo2DEffPi0Phojet->Draw("copy");
//******************************************************************************************************************************************************
cout << "Pi0 efficiency 8TeV" << endl;
TCanvas* canvasEfficiencyPi08TeV = new TCanvas("canvasEfficiencyPi08TeV","",200,10,1350,1350); // gives the page size
DrawGammaCanvasSettings(canvasEfficiencyPi08TeV, 0.1, 0.02, 0.035, 0.09);
TH2F * histo2DEffPi08TeV;
histo2DEffPi08TeV = new TH2F("histo2DEffPi08TeV","histo2DEffPi08TeV",1000,0,8,2000,0.e-3,3.7e-3 );
SetStyleHistoTH2ForGraphs(histo2DEffPi08TeV, "#it{p}_{T} (GeV/#it{c})","#epsilon_{reco, #pi^{0}}",0.03,0.04, 0.03,0.04, 1.,1.);
histo2DEffPi08TeV->Draw("copy");
DrawGammaSetMarker(histoPi0TrueEffiPtPythia8TeV, markerStyleSpectrum8TeVMC, markerSizePi0PP8TeV, colorPi08TeV, colorPi08TeV);
histoPi0TrueEffiPtPythia8TeV->DrawCopy("e1,same");
// DrawGammaSetMarker(histoPi0TrueEffiPtPythiaPlus8TeV, markerStyleSpectrum8TeVMC+1, markerSizePi0PP8TeV, colorPi08TeV, colorPi08TeV);
// histoPi0TrueEffiPtPythiaPlus8TeV->DrawCopy("e1,same");
// DrawGammaSetMarker(histoPi0TrueEffiPtPhojet8TeV, markerStyleSpectrum8TeVMC+2, markerSizePi0PP8TeV, colorPi08TeV, colorPi08TeV);
// histoPi0TrueEffiPtPhojet8TeV->DrawCopy("e1,same");
TLegend* legendEffiPi08TeV = new TLegend(0.12,0.78,0.83,0.93);
legendEffiPi08TeV->SetFillColor(0);
legendEffiPi08TeV->SetLineColor(0);
legendEffiPi08TeV->SetTextSize(0.027);
legendEffiPi08TeV->SetNColumns(2);
legendEffiPi08TeV->AddEntry(histoPi0TrueEffiPtPythia8TeV,"pp #sqrt{#it{s}} = 8 TeV - Pythia","p");
// legendEffiPi08TeV->AddEntry(histoPi0TrueEffiPtPythiaPlus8TeV,"pp #sqrt{#it{s}} = 8 TeV - Pythia+AddSign","p");
// legendEffiPi08TeV->AddEntry(histoPi0TrueEffiPtPhojet8TeV,"pp #sqrt{#it{s}} = 8 TeV - Phojet","p");
legendEffiPi08TeV->Draw();
canvasEfficiencyPi08TeV->SaveAs(Form("%s/EfficiencyPi08TeV.%s",outputDir.Data(),suffix.Data()));
histo2DEffPi08TeV->Draw("copy");
cout << "Pi0 efficiency 7TeV" << endl;
TCanvas* canvasEfficiencyPi07TeV = new TCanvas("canvasEfficiencyPi07TeV","",200,10,1350,1350); // gives the page size
DrawGammaCanvasSettings(canvasEfficiencyPi07TeV, 0.1, 0.02, 0.035, 0.09);
TH2F * histo2DEffPi07TeV;
histo2DEffPi07TeV = new TH2F("histo2DEffPi07TeV","histo2DEffPi07TeV",1000,0,8,2000,0.e-3,3.7e-3 );
SetStyleHistoTH2ForGraphs(histo2DEffPi07TeV, "#it{p}_{T} (GeV/#it{c})","#epsilon_{reco, #pi^{0}}",0.03,0.04, 0.03,0.04, 1.,1.);
histo2DEffPi07TeV->Draw("copy");
DrawGammaSetMarker(histoPi0TrueEffiPtPythia7TeV, markerStyleSpectrum7TeVMC, markerSizePi0PP7TeV, colorPi07TeV, colorPi07TeV);
histoPi0TrueEffiPtPythia7TeV->DrawCopy("e1,same");
DrawGammaSetMarker(histoPi0TrueEffiPtPythiaPlus7TeV, markerStyleSpectrum7TeVMC+1, markerSizePi0PP7TeV, colorPi07TeV, colorPi07TeV);
histoPi0TrueEffiPtPythiaPlus7TeV->DrawCopy("e1,same");
DrawGammaSetMarker(histoPi0TrueEffiPtPhojet7TeV, markerStyleSpectrum7TeVMC+2, markerSizePi0PP7TeV, colorPi07TeV, colorPi07TeV);
histoPi0TrueEffiPtPhojet7TeV->DrawCopy("e1,same");
TLegend* legendEffiPi07TeV = new TLegend(0.12,0.78,0.83,0.93);
legendEffiPi07TeV->SetFillColor(0);
legendEffiPi07TeV->SetLineColor(0);
legendEffiPi07TeV->SetTextSize(0.027);
legendEffiPi07TeV->SetNColumns(2);
legendEffiPi07TeV->AddEntry(histoPi0TrueEffiPtPythia7TeV,"pp #sqrt{#it{s}} = 7 TeV - Pythia","p");
legendEffiPi07TeV->AddEntry(histoPi0TrueEffiPtPythiaPlus7TeV,"pp #sqrt{#it{s}} = 7 TeV - Pythia+AddSign","p");
legendEffiPi07TeV->AddEntry(histoPi0TrueEffiPtPhojet7TeV,"pp #sqrt{#it{s}} = 7 TeV - Phojet","p");
legendEffiPi07TeV->Draw();
canvasEfficiencyPi07TeV->SaveAs(Form("%s/EfficiencyPi07TeV.%s",outputDir.Data(),suffix.Data()));
histo2DEffPi07TeV->Draw("copy");
cout << "Pi0 efficiency 2760GeV" << endl;
TCanvas* canvasEfficiencyPi02760GeV = new TCanvas("canvasEfficiencyPi02760GeV","",200,10,1350,1350); // gives the page size
DrawGammaCanvasSettings(canvasEfficiencyPi02760GeV, 0.1, 0.02, 0.035, 0.09);
TH2F * histo2DEffPi02760GeV;
histo2DEffPi02760GeV = new TH2F("histo2DEffPi02760GeV","histo2DEffPi02760GeV",1000,0,8,2000,0.e-3,3.7e-3 );
SetStyleHistoTH2ForGraphs(histo2DEffPi02760GeV, "#it{p}_{T} (GeV/#it{c})","#epsilon_{reco, #pi^{0}}",0.03,0.04, 0.03,0.04, 1.,1.);
histo2DEffPi02760GeV->Draw("copy");
DrawGammaSetMarker(histoPi0TrueEffiPtPythia2760GeV, markerStyleSpectrum2760GeVMC, markerSizePi0PP2760GeV, colorPi02760GeV, colorPi02760GeV);
histoPi0TrueEffiPtPythia2760GeV->DrawCopy("e1,same");
DrawGammaSetMarker(histoPi0TrueEffiPtPythiaPlus2760GeV, markerStyleSpectrum2760GeVMC+1, markerSizePi0PP2760GeV, colorPi02760GeV, colorPi02760GeV);
histoPi0TrueEffiPtPythiaPlus2760GeV->DrawCopy("e1,same");
DrawGammaSetMarker(histoPi0TrueEffiPtPhojet2760GeV, markerStyleSpectrum2760GeVMC+2, markerSizePi0PP2760GeV, colorPi02760GeV, colorPi02760GeV);
histoPi0TrueEffiPtPhojet2760GeV->DrawCopy("e1,same");
TLegend* legendEffiPi02760GeV = new TLegend(0.12,0.78,0.83,0.93);
legendEffiPi02760GeV->SetFillColor(0);
legendEffiPi02760GeV->SetLineColor(0);
legendEffiPi02760GeV->SetTextSize(0.027);
legendEffiPi02760GeV->SetNColumns(2);
legendEffiPi02760GeV->AddEntry(histoPi0TrueEffiPtPythia2760GeV,"pp #sqrt{#it{s}} = 2.76 TeV - Pythia","p");
legendEffiPi02760GeV->AddEntry(histoPi0TrueEffiPtPythiaPlus2760GeV,"pp #sqrt{#it{s}} = 2.76 TeV - Pythia+AddSign","p");
legendEffiPi02760GeV->AddEntry(histoPi0TrueEffiPtPhojet2760GeV,"pp #sqrt{#it{s}} = 2.76 TeV - Phojet","p");
legendEffiPi02760GeV->Draw();
canvasEfficiencyPi02760GeV->SaveAs(Form("%s/EfficiencyPi02760GeV.%s",outputDir.Data(),suffix.Data()));
histo2DEffPi02760GeV->Draw("copy");
cout << "Pi0 efficiency 900GeV" << endl;
TCanvas* canvasEfficiencyPi0900GeV = new TCanvas("canvasEfficiencyPi0900GeV","",200,10,1350,1350); // gives the page size