24 #include "TPaveText.h"
25 #include "TDirectory.h"
30 static const Ssiz_t
NTxt = 3;
31 static const Ssiz_t
NPad = 2;
32 static const Ssiz_t
NPar = 3;
33 static const Ssiz_t
NVtx = 4;
43 static const TString
SInDef =
"input/embed_only/final_merge/sPhenixG4_forPtCheck_embedScanOn_embedOnly.pt2040n100pim.d8m5y2023.root";
44 static const TString
SOutDef =
"varyDeltaPtCut.withPtDependentDeltaPtCut.pt2040n100pim.d10m5y2023.root";
51 gErrorIgnoreLevel = kError;
52 cout <<
"\n Beginning delta-pt extractor script..." << endl;
55 const UInt_t nInttTrkMin = 1;
56 const UInt_t nMVtxTrkMin = 2;
57 const UInt_t nTpcTrkMin = 35;
58 const Double_t qualTrkMax = 10.;
59 const Double_t vzTrkMax = 10.;
60 const Double_t ptTrkMin = 0.1;
61 const Double_t ptDeltaMax[
NDPtCuts] = {0.5, 0.25, 0.1, 0.05, 0.03, 0.02, 0.01};
62 const Double_t ptDeltaSig[
NSigCuts] = {1., 1.5, 2., 2.5, 3.};
63 const Double_t normRange[
NRange] = {0.2, 1.2};
66 const TString sPtTrueBase =
"PtTrue";
67 const TString sPtRecoBase =
"PtReco";
68 const TString sPtFracBase =
"PtFrac";
69 const TString sPtDeltaBase =
"DeltaPt";
70 const TString sPtTrkTruBase =
"PtTrkTruth";
71 const TString sPtProjBase =
"DeltaPtProj";
72 const TString sRejCutBase =
"Reject_flatDPtCut";
73 const TString sRejSigBase =
"Reject_sigmaCut";
74 const TString sMuHiBase =
"MeanPlusSigma";
75 const TString sMuLoBase =
"MeanMinusSigma";
76 const TString sSigBase =
"ProjectionSigma";
77 const TString sMuBase =
"ProjectionMean";
78 const TString sEffBase =
"Efficiency";
79 const TString sProjSuffix[
NProj] = {
"_pt05",
"_pt1",
"_pt2",
"_pt5",
"_pt10",
"_pt20",
"_pt30",
"_pt40"};
80 const TString sDPtSuffix[
NDPtCuts] = {
"_dPt50",
"_dPt25",
"_dPt10",
"_dPt05",
"_dPt03",
"_dPt02",
"_dPt01"};
81 const TString sSigSuffix[
NSigCuts] = {
"_sigDPt1",
"_sidDpt15",
"_sigDPt2",
"_sigDPt25",
"_sigDPt3"};
84 const UInt_t iCutToDraw =
NDPtCuts - 3;
85 const UInt_t iSigToDraw =
NSigCuts - 3;
86 const UInt_t nEffRebin = 5;
87 const Bool_t doEffRebin =
true;
88 const TString sBeforeTitle =
"Before #Deltap_{T}/p_{T} cuts";
89 const TString sAfterCutTitle =
"After #Deltap_{T}/p_{T} < 0.03 cut";
90 const TString sAfterSigTitle =
"After 2 #times #sigma(#Deltap_{T}/p_{T}) cut";
93 const TString sTitle =
"";
94 const TString sCounts =
"counts";
95 const TString sPtTrueAxis =
"p_{T}^{true} [GeV/c]";
96 const TString sPtRecoAxis =
"p_{T}^{reco} [GeV/c]";
97 const TString sPtFracAxis =
"p_{T}^{reco} / p_{T}^{true}";
98 const TString sPtDeltaAxis =
"#Deltap_{T} / p_{T}^{reco}";
99 const TString sDeltaCutAxis =
"max #Deltap_{T} / p_{T}^{reco}";
100 const TString sSigmaCutAxis =
"n #times #sigma(#Deltap_{T} / p_{T}^{reco})";
101 const TString sSigProjAxis =
"#sigma(#Deltap_{T} / p_{T}^{reco})";
102 const TString sMuProjAxis =
"#mu(#Deltap_{T} / p_{T}^{reco}) #pm (n #times #sigma(#Deltap_{T} / p_{T}^{reco}))";
103 const TString sRejectAxis =
"rejection factor";
104 const TString sEffAxis =
"#epsilon_{trk}";
107 const Double_t ptProj[
NProj] = {0.5, 1., 2., 5., 10., 20., 30., 40.};
108 const Double_t sigHiGuess[
NPar] = {1., -1., 1.};
109 const Double_t sigLoGuess[
NPar] = {1., -1., 1.};
110 const Double_t deltaFitRange[
NRange] = {0., 0.1};
111 const Double_t ptFitRange[
NRange] = {0.5, 40.};
114 const UInt_t fColTrue(923);
115 const UInt_t fColPure(923);
116 const UInt_t fColTrk(809);
117 const UInt_t fMarTrue(20);
118 const UInt_t fMarPure(20);
119 const UInt_t fMarTrk(46);
120 const UInt_t fColProj[
NProj] = {799, 633, 899, 617, 879, 859, 839, 819};
121 const UInt_t fMarProj[
NProj] = {20, 22, 23, 21, 33, 34, 47, 20};
122 const UInt_t fColCut[
NDPtCuts] = {899, 909, 879, 889, 859, 869, 839};
123 const UInt_t fMarCut[
NDPtCuts] = {24, 26, 32, 25, 27, 28, 30};
124 const Float_t rPtRange[
NRange] = {0., 60.};
125 const Float_t rFracRange[
NRange] = {0., 4.};
126 const Float_t rDeltaRange[
NRange] = {0., 0.1};
129 const UInt_t fColFit[
NProj] = {803, 636, 893, 620, 883, 863, 843, 813};
130 const UInt_t fColSigFit[
NSigCuts] = {893, 903, 873, 883, 863};
131 const UInt_t fColSig[
NSigCuts] = {899, 909, 879, 889, 859};
132 const UInt_t fMarSig[
NSigCuts] = {24, 26, 32, 25, 27};
135 const TString sLegTrue(
"truth");
136 const TString sLegTrack(
"tracks (w/ cuts)");
137 const TString sLegMu(
"Mean #Deltap_{T} / p_{T}^{reco} (n = 0)");
138 const TString sInfo[
NTxt] = {
139 "#bf{#it{sPHENIX}} Simulation",
140 "100 #pi^{-}/event, p_{T} #in (20, 40) GeV/c",
143 const TString sLegProj[
NProj] = {
144 "p_{T}^{reco} = 0.5 GeV/c",
145 "p_{T}^{reco} = 1 GeV/c",
146 "p_{T}^{reco} = 2 GeV/c",
147 "p_{T}^{reco} = 5 GeV/c",
148 "p_{T}^{reco} = 10 GeV/c",
149 "p_{T}^{reco} = 20 GeV/c",
150 "p_{T}^{reco} = 30 GeV/c",
151 "p_{T}^{reco} = 40 GeV/c"
154 "#Deltap_{T} / p_{T}^{reco} < 0.5",
155 "#Deltap_{T} / p_{T}^{reco} < 0.25",
156 "#Deltap_{T} / p_{T}^{reco} < 0.1",
157 "#Deltap_{T} / p_{T}^{reco} < 0.05",
158 "#Deltap_{T} / p_{T}^{reco} < 0.03",
159 "#Deltap_{T} / p_{T}^{reco} < 0.02",
160 "#Deltap_{T} / p_{T}^{reco} < 0.01"
162 const TString sLegProjSig[
NSigCuts] = {
169 const TString sLegDelta[
NSigCuts] = {
170 "1 #times #sigma(#Deltap_{T} / p_{T}^{reco})",
171 "1.5 #times #sigma(#Deltap_{T} / p_{T}^{reco})",
172 "2 #times #sigma(#Deltap_{T} / p_{T}^{reco})",
173 "2.5 #times #sigma(#Deltap_{T} / p_{T}^{reco})",
174 "3 #times #sigma(#Deltap_{T} / p_{T}^{reco})"
177 "#Deltap_{T} / p_{T}^{reco} #in 1 #times sigma(#Deltap_{T} / p_{T}^{reco})",
178 "#Deltap_{T} / p_{T}^{reco} #in 1.5 #times sigma(#Deltap_{T} / p_{T}^{reco})",
179 "#Deltap_{T} / p_{T}^{reco} #in 2 #times sigma(#Deltap_{T} / p_{T}^{reco})",
180 "#Deltap_{T} / p_{T}^{reco} #in 2.5 #times sigma(#Deltap_{T} / p_{T}^{reco})",
181 "#Deltap_{T} / p_{T}^{reco} #in 3 #times sigma(#Deltap_{T} / p_{T}^{reco})"
183 const TString sTrkCuts[
NTrkCuts] = {
185 "N_{hit}^{intt} #geq 1",
186 "N_{hit}^{mvtx} > 2",
187 "N_{hit}^{tpc} > 35",
188 "p_{T}^{reco} > 0.1 GeV/c",
193 TFile *fOutput =
new TFile(sOutput.Data(),
"recreate");
194 TFile *fInput =
new TFile(sInput.Data(),
"read");
195 if (!fInput || !fOutput) {
196 cerr <<
"PANIC: couldn't open a file!\n"
197 <<
" fInput = " << fInput <<
"\n"
198 <<
" fOutput = " << fOutput <<
"\n"
202 cout <<
" Opened files." << endl;
205 TNtuple *ntTrack = (TNtuple*) fInput ->
Get(
SInTrack.Data());
206 TNtuple *ntTruth = (TNtuple*) fInput ->
Get(
SInTruth.Data());
207 if (!ntTrack || !ntTruth) {
208 cerr <<
"PANIC: couldn't grab aninput tuple!\n"
209 <<
" ntTrack = " << ntTrack <<
"\n"
210 <<
" ntTruth = " << ntTruth <<
"\n"
214 cout <<
" Grabbed input tuples." << endl;
220 Float_t trk_crossing;
228 Float_t trk_deltaeta;
229 Float_t trk_deltaphi;
248 Float_t trk_vertexID;
253 Float_t trk_dca2dsigma;
255 Float_t trk_dca3dxysigma;
257 Float_t trk_dca3dzsigma;
261 Float_t trk_gtrackID;
289 Float_t trk_gprimary;
290 Float_t trk_nfromtruth;
292 Float_t trk_ntrumaps;
293 Float_t trk_ntruintt;
296 Float_t trk_ntrutpc1;
297 Float_t trk_ntrutpc11;
298 Float_t trk_ntrutpc2;
299 Float_t trk_ntrutpc3;
300 Float_t trk_layersfromtruth;
301 Float_t trk_nhittpcall;
302 Float_t trk_nhittpcin;
303 Float_t trk_nhittpcmid;
304 Float_t trk_nhittpcout;
305 Float_t trk_nclusall;
306 Float_t trk_nclustpc;
307 Float_t trk_nclusintt;
308 Float_t trk_nclusmaps;
309 Float_t trk_nclusmms;
314 Float_t tru_gntracks;
315 Float_t tru_gtrackID;
351 Float_t tru_gprimary;
360 Float_t tru_deltaeta;
361 Float_t tru_deltaphi;
380 Float_t tru_vertexID;
385 Float_t tru_dca2dsigma;
387 Float_t tru_dca3dxysigma;
389 Float_t tru_dca3dzsigma;
393 Float_t tru_nfromtruth;
395 Float_t tru_ntrumaps;
396 Float_t tru_ntruintt;
399 Float_t tru_ntrutpc1;
400 Float_t tru_ntrutpc11;
401 Float_t tru_ntrutpc2;
402 Float_t tru_ntrutpc3;
403 Float_t tru_layersfromtruth;
404 Float_t tru_nhittpcall;
405 Float_t tru_nhittpcin;
406 Float_t tru_nhittpcmid;
407 Float_t tru_nhittpcout;
408 Float_t tru_nclusall;
409 Float_t tru_nclustpc;
410 Float_t tru_nclusintt;
411 Float_t tru_nclusmaps;
412 Float_t tru_nclusmms;
611 cout <<
" Set track tuple branches." << endl;
620 TH1D *hPtDeltaProj[
NProj];
632 TH2D *hPtDeltaVsFrac;
633 TH2D *hPtDeltaVsTrue;
634 TH2D *hPtDeltaVsTrack;
635 TH2D *hPtTrueVsTrack;
646 const UInt_t nPtBins(1000);
647 const UInt_t nFracBins(1000);
648 const UInt_t nDeltaBins(5000);
649 const Float_t rPtBins[
NRange] = {0., 100.};
650 const Float_t rFracBins[
NRange] = {0., 10.};
651 const Float_t rDeltaBins[
NRange] = {0., 5.};
654 TString sPtTruth(
"h");
655 TString sPtDelta(
"h");
656 TString sPtTrack(
"h");
657 TString sPtFrac(
"h");
658 TString sPtTrkTru(
"h");
659 sPtTruth.Append(sPtTrueBase.Data());
660 sPtDelta.Append(sPtDeltaBase.Data());
661 sPtTrack.Append(sPtRecoBase.Data());
662 sPtFrac.Append(sPtFracBase.Data());
663 sPtTrkTru.Append(sPtTrkTruBase.Data());
665 TString sPtDeltaVsFrac(
"h");
666 TString sPtDeltaVsTrue(
"h");
667 TString sPtDeltaVsTrack(
"h");
668 TString sPtTrueVsTrack(
"h");
669 sPtDeltaVsFrac.Append(sPtDeltaBase.Data());
670 sPtDeltaVsTrue.Append(sPtDeltaBase.Data());
671 sPtDeltaVsTrack.Append(sPtDeltaBase.Data());
672 sPtTrueVsTrack.Append(sPtTrueBase.Data());
673 sPtDeltaVsFrac.Append(
"Vs");
674 sPtDeltaVsTrue.Append(
"Vs");
675 sPtDeltaVsTrack.Append(
"Vs");
676 sPtTrueVsTrack.Append(
"Vs");
677 sPtDeltaVsFrac.Append(sPtFracBase.Data());
678 sPtDeltaVsTrue.Append(sPtTrueBase.Data());
679 sPtDeltaVsTrack.Append(sPtRecoBase.Data());
680 sPtTrueVsTrack.Append(sPtRecoBase.Data());
683 TString sPtProj[
NProj];
684 for (Ssiz_t iProj = 0; iProj <
NProj; iProj++) {
685 sPtProj[iProj] =
"h";
686 sPtProj[iProj].Append(sPtProjBase.Data());
687 sPtProj[iProj].Append(sProjSuffix[iProj].
Data());
695 TString sPtDeltaVsFracCut[
NDPtCuts];
696 TString sPtDeltaVsTrueCut[
NDPtCuts];
697 TString sPtDeltaVsTrackCut[
NDPtCuts];
698 TString sPtTrueVsTrackCut[
NDPtCuts];
699 for (Ssiz_t iCut = 0; iCut <
NDPtCuts; iCut++) {
700 sPtDeltaCut[iCut] =
"h";
701 sPtTrackCut[iCut] =
"h";
702 sPtFracCut[iCut] =
"h";
703 sPtTrkTruCut[iCut] =
"h";
704 sPtDeltaCut[iCut].Append(sPtDeltaBase.Data());
705 sPtTrackCut[iCut].Append(sPtRecoBase.Data());
706 sPtFracCut[iCut].Append(sPtFracBase.Data());
707 sPtTrkTruCut[iCut].Append(sPtTrkTruBase.Data());
708 sPtDeltaCut[iCut].Append(sDPtSuffix[iCut].
Data());
709 sPtTrackCut[iCut].Append(sDPtSuffix[iCut].
Data());
710 sPtFracCut[iCut].Append(sDPtSuffix[iCut].
Data());
711 sPtTrkTruCut[iCut].Append(sDPtSuffix[iCut].
Data());
713 sPtDeltaVsFracCut[iCut] =
"h";
714 sPtDeltaVsTrueCut[iCut] =
"h";
715 sPtDeltaVsTrackCut[iCut] =
"h";
716 sPtTrueVsTrackCut[iCut] =
"h";
717 sPtDeltaVsFracCut[iCut].Append(sPtDeltaBase.Data());
718 sPtDeltaVsFracCut[iCut].Append(sPtDeltaBase.Data());
719 sPtDeltaVsTrueCut[iCut].Append(sPtDeltaBase.Data());
720 sPtDeltaVsTrackCut[iCut].Append(sPtDeltaBase.Data());
721 sPtTrueVsTrackCut[iCut].Append(sPtTrueBase.Data());
722 sPtDeltaVsFracCut[iCut].Append(
"Vs");
723 sPtDeltaVsTrueCut[iCut].Append(
"Vs");
724 sPtDeltaVsTrackCut[iCut].Append(
"Vs");
725 sPtTrueVsTrackCut[iCut].Append(
"Vs");
726 sPtDeltaVsFracCut[iCut].Append(sPtFracBase.Data());
727 sPtDeltaVsTrueCut[iCut].Append(sPtTrueBase.Data());
728 sPtDeltaVsTrackCut[iCut].Append(sPtRecoBase.Data());
729 sPtTrueVsTrackCut[iCut].Append(sPtRecoBase.Data());
730 sPtDeltaVsFracCut[iCut].Append(sDPtSuffix[iCut].
Data());
731 sPtDeltaVsTrueCut[iCut].Append(sDPtSuffix[iCut].
Data());
732 sPtDeltaVsTrackCut[iCut].Append(sDPtSuffix[iCut].
Data());
733 sPtTrueVsTrackCut[iCut].Append(sDPtSuffix[iCut].
Data());
741 TString sPtDeltaVsFracSig[
NSigCuts];
742 TString sPtDeltaVsTrueSig[
NSigCuts];
743 TString sPtDeltaVsTrackSig[
NSigCuts];
744 TString sPtTrueVsTrackSig[
NSigCuts];
745 for (Ssiz_t iSig = 0; iSig <
NSigCuts; iSig++) {
746 sPtDeltaSig[iSig] =
"h";
747 sPtTrackSig[iSig] =
"h";
748 sPtFracSig[iSig] =
"h";
749 sPtTrkTruSig[iSig] =
"h";
750 sPtDeltaSig[iSig].Append(sPtDeltaBase.Data());
751 sPtTrackSig[iSig].Append(sPtRecoBase.Data());
752 sPtFracSig[iSig].Append(sPtFracBase.Data());
753 sPtTrkTruSig[iSig].Append(sPtTrkTruBase.Data());
754 sPtDeltaSig[iSig].Append(sSigSuffix[iSig].
Data());
755 sPtTrackSig[iSig].Append(sSigSuffix[iSig].
Data());
756 sPtFracSig[iSig].Append(sSigSuffix[iSig].
Data());
757 sPtTrkTruSig[iSig].Append(sSigSuffix[iSig].
Data());
759 sPtDeltaVsFracSig[iSig] =
"h";
760 sPtDeltaVsTrueSig[iSig] =
"h";
761 sPtDeltaVsTrackSig[iSig] =
"h";
762 sPtTrueVsTrackSig[iSig] =
"h";
763 sPtDeltaVsFracSig[iSig].Append(sPtDeltaBase.Data());
764 sPtDeltaVsFracSig[iSig].Append(sPtDeltaBase.Data());
765 sPtDeltaVsTrueSig[iSig].Append(sPtDeltaBase.Data());
766 sPtDeltaVsTrackSig[iSig].Append(sPtDeltaBase.Data());
767 sPtTrueVsTrackSig[iSig].Append(sPtTrueBase.Data());
768 sPtDeltaVsFracSig[iSig].Append(
"Vs");
769 sPtDeltaVsTrueSig[iSig].Append(
"Vs");
770 sPtDeltaVsTrackSig[iSig].Append(
"Vs");
771 sPtTrueVsTrackSig[iSig].Append(
"Vs");
772 sPtDeltaVsFracSig[iSig].Append(sPtFracBase.Data());
773 sPtDeltaVsTrueSig[iSig].Append(sPtTrueBase.Data());
774 sPtDeltaVsTrackSig[iSig].Append(sPtRecoBase.Data());
775 sPtTrueVsTrackSig[iSig].Append(sPtRecoBase.Data());
776 sPtDeltaVsFracSig[iSig].Append(sSigSuffix[iSig].
Data());
777 sPtDeltaVsTrueSig[iSig].Append(sSigSuffix[iSig].
Data());
778 sPtDeltaVsTrackSig[iSig].Append(sSigSuffix[iSig].
Data());
779 sPtTrueVsTrackSig[iSig].Append(sSigSuffix[iSig].
Data());
783 hPtTruth =
new TH1D(sPtTruth.Data(),
"", nPtBins, rPtBins[0], rPtBins[1]);
784 hPtDelta =
new TH1D(sPtDelta.Data(),
"", nDeltaBins, rDeltaBins[0], rDeltaBins[1]);
785 hPtTrack =
new TH1D(sPtTrack.Data(),
"", nPtBins, rPtBins[0], rPtBins[1]);
786 hPtFrac =
new TH1D(sPtFrac.Data(),
"", nFracBins, rFracBins[0], rFracBins[1]);
787 hPtTrkTru =
new TH1D(sPtTrkTru.Data(),
"", nPtBins, rPtBins[0], rPtBins[1]);
792 hPtTrkTru -> Sumw2();
794 hPtDeltaVsFrac =
new TH2D(sPtDeltaVsFrac.Data(),
"", nFracBins, rFracBins[0], rFracBins[1], nDeltaBins, rDeltaBins[0], rDeltaBins[1]);
795 hPtDeltaVsTrue =
new TH2D(sPtDeltaVsTrue.Data(),
"", nPtBins, rPtBins[0], rPtBins[1], nDeltaBins, rDeltaBins[0], rDeltaBins[1]);
796 hPtDeltaVsTrack =
new TH2D(sPtDeltaVsTrack.Data(),
"", nPtBins, rPtBins[0], rPtBins[1], nDeltaBins, rDeltaBins[0], rDeltaBins[1]);
797 hPtTrueVsTrack =
new TH2D(sPtTrueVsTrack.Data(),
"", nPtBins, rPtBins[0], rPtBins[1], nPtBins, rPtBins[0], rPtBins[1]);
798 hPtDeltaVsFrac -> Sumw2();
799 hPtDeltaVsTrue -> Sumw2();
800 hPtDeltaVsTrack -> Sumw2();
801 hPtTrueVsTrack -> Sumw2();
804 for (Ssiz_t iProj = 0; iProj <
NProj; iProj++) {
805 hPtDeltaProj[iProj] =
new TH1D(sPtProj[iProj].
Data(),
"", nDeltaBins, rDeltaBins[0], rDeltaBins[1]);
806 hPtDeltaProj[iProj] -> Sumw2();
810 for (Ssiz_t iCut = 0; iCut <
NDPtCuts; iCut++) {
811 hPtDeltaCut[iCut] =
new TH1D(sPtDeltaCut[iCut].
Data(),
"", nDeltaBins, rDeltaBins[0], rDeltaBins[1]);
812 hPtTrackCut[iCut] =
new TH1D(sPtTrackCut[iCut].
Data(),
"", nPtBins, rPtBins[0], rPtBins[1]);
813 hPtFracCut[iCut] =
new TH1D(sPtFracCut[iCut].
Data(),
"", nFracBins, rFracBins[0], rFracBins[1]);
814 hPtTrkTruCut[iCut] =
new TH1D(sPtTrkTruCut[iCut].
Data(),
"", nPtBins, rPtBins[0], rPtBins[1]);
815 hPtDeltaCut[iCut] -> Sumw2();
816 hPtTrackCut[iCut] -> Sumw2();
817 hPtFracCut[iCut] -> Sumw2();
818 hPtTrkTruCut[iCut] -> Sumw2();
820 hPtDeltaVsFracCut[iCut] =
new TH2D(sPtDeltaVsFracCut[iCut].
Data(),
"", nFracBins, rFracBins[0], rFracBins[1], nDeltaBins, rDeltaBins[0], rDeltaBins[1]);
821 hPtDeltaVsTrueCut[iCut] =
new TH2D(sPtDeltaVsTrueCut[iCut].
Data(),
"", nPtBins, rPtBins[0], rPtBins[1], nDeltaBins, rDeltaBins[0], rDeltaBins[1]);
822 hPtDeltaVsTrackCut[iCut] =
new TH2D(sPtDeltaVsTrackCut[iCut].
Data(),
"", nPtBins, rPtBins[0], rPtBins[1], nDeltaBins, rDeltaBins[0], rDeltaBins[1]);
823 hPtTrueVsTrackCut[iCut] =
new TH2D(sPtTrueVsTrackCut[iCut].
Data(),
"", nPtBins, rPtBins[0], rPtBins[1], nPtBins, rPtBins[0], rPtBins[1]);
824 hPtDeltaVsFracCut[iCut] -> Sumw2();
825 hPtDeltaVsTrueCut[iCut] -> Sumw2();
826 hPtDeltaVsTrackCut[iCut] -> Sumw2();
827 hPtTrueVsTrackCut[iCut] -> Sumw2();
831 for (Ssiz_t iSig = 0; iSig <
NSigCuts; iSig++) {
832 hPtDeltaSig[iSig] =
new TH1D(sPtDeltaSig[iSig].
Data(),
"", nDeltaBins, rDeltaBins[0], rDeltaBins[1]);
833 hPtTrackSig[iSig] =
new TH1D(sPtTrackSig[iSig].
Data(),
"", nPtBins, rPtBins[0], rPtBins[1]);
834 hPtFracSig[iSig] =
new TH1D(sPtFracSig[iSig].
Data(),
"", nFracBins, rFracBins[0], rFracBins[1]);
835 hPtTrkTruSig[iSig] =
new TH1D(sPtTrkTruSig[iSig].
Data(),
"", nPtBins, rPtBins[0], rPtBins[1]);
836 hPtDeltaSig[iSig] -> Sumw2();
837 hPtTrackSig[iSig] -> Sumw2();
838 hPtFracSig[iSig] -> Sumw2();
839 hPtTrkTruSig[iSig] -> Sumw2();
841 hPtDeltaVsFracSig[iSig] =
new TH2D(sPtDeltaVsFracSig[iSig].
Data(),
"", nFracBins, rFracBins[0], rFracBins[1], nDeltaBins, rDeltaBins[0], rDeltaBins[1]);
842 hPtDeltaVsTrueSig[iSig] =
new TH2D(sPtDeltaVsTrueSig[iSig].
Data(),
"", nPtBins, rPtBins[0], rPtBins[1], nDeltaBins, rDeltaBins[0], rDeltaBins[1]);
843 hPtDeltaVsTrackSig[iSig] =
new TH2D(sPtDeltaVsTrackSig[iSig].
Data(),
"", nPtBins, rPtBins[0], rPtBins[1], nDeltaBins, rDeltaBins[0], rDeltaBins[1]);
844 hPtTrueVsTrackSig[iSig] =
new TH2D(sPtTrueVsTrackSig[iSig].
Data(),
"", nPtBins, rPtBins[0], rPtBins[1], nPtBins, rPtBins[0], rPtBins[1]);
845 hPtDeltaVsFracSig[iSig] -> Sumw2();
846 hPtDeltaVsTrueSig[iSig] -> Sumw2();
847 hPtDeltaVsTrackSig[iSig] -> Sumw2();
848 hPtTrueVsTrackSig[iSig] -> Sumw2();
852 const Long64_t nTrks = ntTrack -> GetEntries();
853 const Long64_t nTrus = ntTruth -> GetEntries();
854 cout <<
" Beginning tuple loops: " << nTrks <<
" reco. tracks and " << nTrus <<
" particles to process" << endl;
857 Double_t muProj[
NProj];
858 Double_t sigProj[
NProj];
861 for (Ssiz_t iProj = 0; iProj <
NProj; iProj++) {
865 for (Ssiz_t iSig = 0; iSig <
NSigCuts; iSig++) {
866 for (Ssiz_t iProj = 0; iProj <
NProj; iProj++) {
867 muHiProj[iSig][iProj] = 0.;
868 muLoProj[iSig][iProj] = 0.;
879 for (Ssiz_t iCut = 0; iCut <
NDPtCuts; iCut++) {
884 for (Ssiz_t iSig = 0; iSig <
NSigCuts; iSig++) {
889 cout <<
" First loop over reco. tracks:" << endl;
892 Long64_t nBytesTrk = 0;
893 for (Long64_t iTrk = 0; iTrk < nTrks; iTrk++) {
896 const Long64_t bytesTrk = ntTrack ->
GetEntry(iTrk);
898 cerr <<
"WARNING: something wrong with track #" << iTrk <<
"! Aborting loop!" << endl;
901 nBytesTrk += bytesTrk;
904 const Long64_t iProgTrk = iTrk + 1;
912 const Double_t ptFrac = trk_pt / trk_gpt;
913 const Double_t ptDelta = trk_deltapt / trk_pt;
916 const Bool_t isInZVtxCut = (abs(trk_vz) < vzTrkMax);
917 const Bool_t isInInttCut = (trk_nintt >= nInttTrkMin);
918 const Bool_t isInMVtxCut = (trk_nlmaps > nMVtxTrkMin);
919 const Bool_t isInTpcCut = (trk_ntpc > nTpcTrkMin);
920 const Bool_t isInPtCut = (trk_pt > ptTrkMin);
921 const Bool_t isInQualCut = (trk_quality < qualTrkMax);
922 const Bool_t isGoodTrk = (isInZVtxCut && isInInttCut && isInMVtxCut && isInTpcCut && isInPtCut && isInQualCut);
923 if (!isGoodTrk)
continue;
926 cout <<
"CHECK intt = " << trk_nintt <<
", mvtx = " << trk_nlmaps <<
", tpc = " << trk_ntpc <<
", quality = " << trk_quality <<
", ptTrue = " << trk_gpt
927 <<
", dcaXY = " << trk_dca3dxy <<
", dcaZ = " << trk_dca3dz << endl;
931 hPtDelta ->
Fill(ptDelta);
932 hPtTrack ->
Fill(trk_pt);
933 hPtFrac ->
Fill(ptFrac);
934 hPtTrkTru ->
Fill(trk_gpt);
935 hPtDeltaVsFrac ->
Fill(ptFrac, ptDelta);
936 hPtDeltaVsTrue ->
Fill(trk_gpt, ptDelta);
937 hPtDeltaVsTrack ->
Fill(trk_pt, ptDelta);
938 hPtTrueVsTrack ->
Fill(trk_pt, trk_gpt);
941 const Bool_t isNormalTrk = ((ptFrac > normRange[0]) && (ptFrac < normRange[1]));
942 for (Ssiz_t iCut = 0; iCut <
NDPtCuts; iCut++) {
943 const Bool_t isInDeltaPtCut = (ptDelta < ptDeltaMax[iCut]);
944 if (isInDeltaPtCut) {
947 hPtDeltaCut[iCut] ->
Fill(ptDelta);
948 hPtTrackCut[iCut] ->
Fill(trk_pt);
949 hPtFracCut[iCut] ->
Fill(ptFrac);
950 hPtTrkTruCut[iCut] ->
Fill(trk_gpt);
951 hPtDeltaVsFracCut[iCut] ->
Fill(ptFrac, ptDelta);
952 hPtDeltaVsTrueCut[iCut] ->
Fill(trk_gpt, ptDelta);
953 hPtDeltaVsTrackCut[iCut] ->
Fill(trk_pt, ptDelta);
954 hPtTrueVsTrackCut[iCut] ->
Fill(trk_pt, trk_gpt);
965 cout <<
" First loop over reco. tracks finished!" << endl;
968 for (Ssiz_t iCut = 0; iCut <
NDPtCuts; iCut++) {
969 rejCut[iCut] = (Double_t) nNormCut[iCut] / (Double_t) nWeirdCut[iCut];
971 cout <<
" Calculated flat delta-pt rejection factors." << endl;
974 TString sFitProj[
NProj];
975 for (Ssiz_t iProj = 0; iProj <
NProj; iProj++) {
976 sFitProj[iProj] =
"f";
977 sFitProj[iProj].Append(sPtProjBase.Data());
978 sFitProj[iProj].Append(sProjSuffix[iProj].
Data());
982 const UInt_t fWidFit = 2;
983 const UInt_t fLinFit = 1;
985 TF1 *fPtDeltaProj[
NProj];
986 for (Ssiz_t iProj = 0; iProj <
NProj; iProj++) {
989 const UInt_t iBinProj = hPtDeltaVsTrack ->
GetXaxis() -> FindBin(ptProj[iProj]);
990 hPtDeltaProj[iProj] = hPtDeltaVsTrack ->
ProjectionY(sPtProj[iProj], iBinProj, iBinProj,
"");
993 const Float_t ampGuess = hPtDeltaProj[iProj] -> GetMaximum();
994 const Float_t muGuess = hPtDeltaProj[iProj] -> GetMean();
995 const Float_t sigGuess = hPtDeltaProj[iProj] -> GetRMS();
998 fPtDeltaProj[iProj] =
new TF1(sFitProj[iProj].
Data(),
"gaus", deltaFitRange[0], deltaFitRange[1]);
1005 hPtDeltaProj[iProj] -> Fit(sFitProj[iProj].
Data(),
"R");
1008 muProj[iProj] = fPtDeltaProj[iProj] ->
GetParameter(1);
1009 sigProj[iProj] = fPtDeltaProj[iProj] ->
GetParameter(2);
1010 for (Ssiz_t iSig = 0; iSig <
NSigCuts; iSig++) {
1011 muHiProj[iSig][iProj] = muProj[iProj] + (ptDeltaSig[iSig] * sigProj[iProj]);
1012 muLoProj[iSig][iProj] = muProj[iProj] - (ptDeltaSig[iSig] * sigProj[iProj]);
1015 cout <<
" Obtained delta-pt projections, fits, and sigmas." << endl;
1018 TString sMuProj(
"gr");
1019 TString sSigProj(
"gr");
1020 sMuProj.Append(sMuBase.Data());
1021 sSigProj.Append(sSigBase.Data());
1027 for (Ssiz_t iSig = 0; iSig <
NSigCuts; iSig++) {
1028 sGrMuHiProj[iSig] =
"gr";
1029 sGrMuLoProj[iSig] =
"gr";
1030 sFnMuHiProj[iSig] =
"f";
1031 sFnMuLoProj[iSig] =
"f";
1032 sGrMuHiProj[iSig].Append(sMuHiBase.Data());
1033 sGrMuLoProj[iSig].Append(sMuLoBase.Data());
1034 sFnMuHiProj[iSig].Append(sMuHiBase.Data());
1035 sFnMuLoProj[iSig].Append(sMuLoBase.Data());
1036 sGrMuHiProj[iSig].Append(sSigSuffix[iSig].
Data());
1037 sGrMuLoProj[iSig].Append(sSigSuffix[iSig].
Data());
1038 sFnMuHiProj[iSig].Append(sSigSuffix[iSig].
Data());
1039 sFnMuLoProj[iSig].Append(sSigSuffix[iSig].
Data());
1043 TGraph *grMuProj =
new TGraph(NProj, ptProj, muProj);
1044 TGraph *grSigProj =
new TGraph(NProj, ptProj, sigProj);
1046 grSigProj ->
SetName(sSigProj);
1052 for (Ssiz_t iSig = 0; iSig <
NSigCuts; iSig++) {
1055 grMuHiProj[iSig] =
new TGraph(NProj, ptProj, muHiProj[iSig]);
1056 grMuLoProj[iSig] =
new TGraph(NProj, ptProj, muLoProj[iSig]);
1057 grMuHiProj[iSig] ->
SetName(sGrMuHiProj[iSig].
Data());
1058 grMuLoProj[iSig] ->
SetName(sGrMuLoProj[iSig].
Data());
1061 fMuHiProj[iSig] =
new TF1(sFnMuHiProj[iSig].
Data(),
"pol2", rPtRange[0], rPtRange[1]);
1062 fMuLoProj[iSig] =
new TF1(sFnMuLoProj[iSig].
Data(),
"pol2", rPtRange[0], rPtRange[1]);
1077 grMuHiProj[iSig] -> Fit(sFnMuHiProj[iSig].
Data(),
"",
"", ptFitRange[0], ptFitRange[1]);
1078 grMuLoProj[iSig] -> Fit(sFnMuLoProj[iSig].
Data(),
"",
"", ptFitRange[0], ptFitRange[1]);
1080 cout <<
" Created and fit sigma graphs.\n"
1081 <<
" Second loop over reco. tracks:"
1086 for (Long64_t iTrk = 0; iTrk < nTrks; iTrk++) {
1089 const Long64_t bytesTrk = ntTrack ->
GetEntry(iTrk);
1090 if (bytesTrk < 0.) {
1091 cerr <<
"WARNING: something wrong with track #" << iTrk <<
"! Aborting loop!" << endl;
1094 nBytesTrk += bytesTrk;
1097 const Long64_t iProgTrk = iTrk + 1;
1098 if (iProgTrk == nTrks) {
1099 cout <<
" Processing track " << iProgTrk <<
"/" << nTrks <<
"..." << endl;
1101 cout <<
" Processing track " << iProgTrk <<
"/" << nTrks <<
"...\r" <<
flush;
1105 const Double_t ptFrac = trk_pt / trk_gpt;
1106 const Double_t ptDelta = trk_deltapt / trk_pt;
1109 const Bool_t isInZVtxCut = (abs(trk_vz) < vzTrkMax);
1110 const Bool_t isInInttCut = (trk_nintt >= nInttTrkMin);
1111 const Bool_t isInMVtxCut = (trk_nlmaps > nMVtxTrkMin);
1112 const Bool_t isInTpcCut = (trk_ntpc > nTpcTrkMin);
1113 const Bool_t isInPtCut = (trk_pt > ptTrkMin);
1114 const Bool_t isInQualCut = (trk_quality < qualTrkMax);
1115 const Bool_t isGoodTrk = (isInZVtxCut && isInInttCut && isInMVtxCut && isInTpcCut && isInPtCut && isInQualCut);
1116 if (!isGoodTrk)
continue;
1119 const Bool_t isNormalTrk = ((ptFrac > normRange[0]) && (ptFrac < normRange[1]));
1120 for (Ssiz_t iSig = 0; iSig <
NSigCuts; iSig++) {
1123 const Float_t ptDeltaMin = fMuLoProj[iSig] -> Eval(trk_pt);
1124 const Float_t ptDeltaMax = fMuHiProj[iSig] -> Eval(trk_pt);
1126 const Bool_t isInDeltaPtSigma = ((ptDelta >= ptDeltaMin) && (ptDelta <= ptDeltaMax));
1127 if (isInDeltaPtSigma) {
1130 hPtDeltaSig[iSig] ->
Fill(ptDelta);
1131 hPtTrackSig[iSig] ->
Fill(trk_pt);
1132 hPtFracSig[iSig] ->
Fill(ptFrac);
1133 hPtTrkTruSig[iSig] ->
Fill(trk_gpt);
1134 hPtDeltaVsFracSig[iSig] ->
Fill(ptFrac, ptDelta);
1135 hPtDeltaVsTrueSig[iSig] ->
Fill(trk_gpt, ptDelta);
1136 hPtDeltaVsTrackSig[iSig] ->
Fill(trk_pt, ptDelta);
1137 hPtTrueVsTrackSig[iSig] ->
Fill(trk_pt, trk_gpt);
1148 cout <<
" Second loop over reco. tracks finished!" << endl;
1151 for (Ssiz_t iSig = 0; iSig <
NSigCuts; iSig++) {
1152 rejSig[iSig] = (Double_t) nNormSig[iSig] / (Double_t) nWeirdSig[iSig];
1154 cout <<
" Calculated pt-depdendent delta-pt rejection factors.\n"
1155 <<
" Loop over particles:"
1159 Long64_t nBytesTru = 0;
1160 for (Long64_t iTru = 0; iTru < nTrus; iTru++) {
1163 const Long64_t bytesTru = ntTruth ->
GetEntry(iTru);
1164 if (bytesTru < 0.) {
1165 cerr <<
"WARNING: something wrong with particle #" << iTru <<
"! Aborting loop!" << endl;
1168 nBytesTru += bytesTru;
1171 const Long64_t iProgTru = iTru + 1;
1172 if (iProgTru == nTrus) {
1173 cout <<
" Processing particle " << iProgTru <<
"/" << nTrus <<
"..." << endl;
1175 cout <<
" Processing particle" << iProgTru <<
"/" << nTrus <<
"...\r" <<
flush;
1179 const Bool_t isPrimary = (tru_gprimary == 1);
1181 hPtTruth ->
Fill(tru_gpt);
1184 cout <<
" Loop over particles finished!" << endl;
1187 cout <<
" Finished tuple loops! Rejection factors:" << endl;
1190 cout <<
" Flat delta-pt cuts" << endl;
1191 for (Ssiz_t iCut = 0; iCut <
NDPtCuts; iCut++) {
1192 cout <<
" n(Norm, Weird) = (" << nNormCut[iCut] <<
", " << nWeirdCut[iCut] <<
"), rejection = " << rejCut[iCut] << endl;
1196 cout <<
" Pt-dependent delta-pt cuts" << endl;
1197 for (Ssiz_t iSig = 0; iSig <
NSigCuts; iSig++) {
1198 cout <<
" n(Norm, Weird) = (" << nNormSig[iSig] <<
", " << nWeirdSig[iSig] <<
"), rejection = " << rejSig[iSig] << endl;
1202 TString sRejCut(
"gr");
1203 TString sRejSig(
"gr");
1204 sRejCut.Append(sRejCutBase.Data());
1205 sRejSig.Append(sRejSigBase.Data());
1207 TGraph *grRejCut =
new TGraph(NDPtCuts, ptDeltaMax, rejCut);
1208 TGraph *grRejSig =
new TGraph(NSigCuts, ptDeltaSig, rejSig);
1209 grRejCut ->
SetName(sRejCut.Data());
1210 grRejSig ->
SetName(sRejSig.Data());
1211 cout <<
" Made rejection factor graph." << endl;
1215 hPtTruth -> Rebin(nEffRebin);
1216 hPtTrkTru -> Rebin(nEffRebin);
1217 for (Ssiz_t iCut = 0; iCut <
NDPtCuts; iCut++) {
1218 hPtTrkTruCut[iCut] -> Rebin(nEffRebin);
1220 for (Ssiz_t iSig = 0; iSig <
NSigCuts; iSig++) {
1221 hPtTrkTruSig[iSig] -> Rebin(nEffRebin);
1223 cout <<
" Rebinned efficiency histograms." << endl;
1228 sEff.Append(sEffBase.Data());
1232 for (Ssiz_t iCut = 0; iCut <
NDPtCuts; iCut++) {
1233 sEffCut[iCut] =
"h";
1234 sEffCut[iCut].Append(sEffBase.Data());
1235 sEffCut[iCut].Append(sDPtSuffix[iCut].
Data());
1240 for (Ssiz_t iSig = 0; iSig <
NSigCuts; iSig++) {
1241 sEffSig[iSig] =
"h";
1242 sEffSig[iSig].Append(sEffBase.Data());
1243 sEffSig[iSig].Append(sSigSuffix[iSig].
Data());
1246 hEff = (TH1D*) hPtTruth -> Clone();
1248 hEff ->
Reset(
"ICES");
1249 hEff -> Divide(hPtTrkTru, hPtTruth, 1., 1.);
1252 for (Ssiz_t iCut = 0; iCut <
NDPtCuts; iCut++) {
1253 hEffCut[iCut] = (TH1D*) hPtTruth -> Clone();
1255 hEffCut[iCut] ->
Reset(
"ICES");
1256 hEffCut[iCut] -> Divide(hPtTrkTruCut[iCut], hPtTruth, 1., 1.);
1260 for (Ssiz_t iSig = 0; iSig <
NSigCuts; iSig++) {
1261 hEffSig[iSig] = (TH1D*) hPtTruth -> Clone();
1263 hEffSig[iSig] ->
Reset(
"ICES");
1264 hEffSig[iSig] -> Divide(hPtTrkTruSig[iSig], hPtTruth, 1., 1.);
1266 cout <<
" Calculated efficiencies." << endl;
1269 const UInt_t fFil(0);
1270 const UInt_t fLin(1);
1271 const UInt_t fWid(1);
1272 const UInt_t fTxt(42);
1273 const UInt_t fAln(12);
1274 const UInt_t fCnt(1);
1275 const Float_t fLab[
NPad] = {0.074, 0.04};
1276 const Float_t fTit[
NPad] = {0.074, 0.04};
1277 const Float_t fOffX[
NPad] = {1.1, 1.};
1278 const Float_t fOffY[
NPad] = {0.7, 1.3};
1279 const Float_t fOffZ[
NPad] = {1.1, 1.1};
1287 grRejCut ->
SetTitle(sTitle.Data());
1288 grRejCut ->
GetXaxis() -> SetRangeUser(rDeltaRange[0], rDeltaRange[1]);
1290 grRejCut ->
GetXaxis() -> SetTitleFont(fTxt);
1291 grRejCut ->
GetXaxis() -> SetTitleSize(fTit[1]);
1292 grRejCut ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
1293 grRejCut ->
GetXaxis() -> SetLabelFont(fTxt);
1294 grRejCut ->
GetXaxis() -> SetLabelSize(fLab[1]);
1295 grRejCut ->
GetXaxis() -> CenterTitle(fCnt);
1297 grRejCut ->
GetYaxis() -> SetTitleFont(fTxt);
1298 grRejCut ->
GetYaxis() -> SetTitleSize(fTit[1]);
1299 grRejCut ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
1300 grRejCut ->
GetYaxis() -> SetLabelFont(fTxt);
1301 grRejCut ->
GetYaxis() -> SetLabelSize(fLab[1]);
1302 grRejCut ->
GetYaxis() -> CenterTitle(fCnt);
1310 grRejSig ->
SetTitle(sTitle.Data());
1311 grRejSig ->
GetXaxis() -> SetRangeUser(rDeltaRange[0], rDeltaRange[1]);
1313 grRejSig ->
GetXaxis() -> SetTitleFont(fTxt);
1314 grRejSig ->
GetXaxis() -> SetTitleSize(fTit[1]);
1315 grRejSig ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
1316 grRejSig ->
GetXaxis() -> SetLabelFont(fTxt);
1317 grRejSig ->
GetXaxis() -> SetLabelSize(fLab[1]);
1318 grRejSig ->
GetXaxis() -> CenterTitle(fCnt);
1320 grRejSig ->
GetYaxis() -> SetTitleFont(fTxt);
1321 grRejSig ->
GetYaxis() -> SetTitleSize(fTit[1]);
1322 grRejSig ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
1323 grRejSig ->
GetYaxis() -> SetLabelFont(fTxt);
1324 grRejSig ->
GetYaxis() -> SetLabelSize(fLab[1]);
1325 grRejSig ->
GetYaxis() -> CenterTitle(fCnt);
1334 hEff -> SetTitleFont(fTxt);
1335 hEff ->
GetXaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
1337 hEff ->
GetXaxis() -> SetTitleFont(fTxt);
1338 hEff ->
GetXaxis() -> SetTitleSize(fTit[0]);
1339 hEff ->
GetXaxis() -> SetTitleOffset(fOffX[0]);
1340 hEff ->
GetXaxis() -> SetLabelFont(fTxt);
1341 hEff ->
GetXaxis() -> SetLabelSize(fLab[0]);
1342 hEff ->
GetXaxis() -> CenterTitle(fCnt);
1344 hEff ->
GetYaxis() -> SetTitleFont(fTxt);
1345 hEff ->
GetYaxis() -> SetTitleSize(fTit[0]);
1346 hEff ->
GetYaxis() -> SetTitleOffset(fOffY[0]);
1347 hEff ->
GetYaxis() -> SetLabelFont(fTxt);
1348 hEff ->
GetYaxis() -> SetLabelSize(fLab[0]);
1349 hEff ->
GetYaxis() -> CenterTitle(fCnt);
1357 hPtTruth ->
SetTitle(sTitle.Data());
1358 hPtTruth -> SetTitleFont(fTxt);
1359 hPtTruth ->
GetXaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
1361 hPtTruth ->
GetXaxis() -> SetTitleFont(fTxt);
1362 hPtTruth ->
GetXaxis() -> SetTitleSize(fTit[1]);
1363 hPtTruth ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
1364 hPtTruth ->
GetXaxis() -> SetLabelFont(fTxt);
1365 hPtTruth ->
GetXaxis() -> SetLabelSize(fLab[1]);
1366 hPtTruth ->
GetXaxis() -> CenterTitle(fCnt);
1368 hPtTruth ->
GetYaxis() -> SetTitleFont(fTxt);
1369 hPtTruth ->
GetYaxis() -> SetTitleSize(fTit[1]);
1370 hPtTruth ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
1371 hPtTruth ->
GetYaxis() -> SetLabelFont(fTxt);
1372 hPtTruth ->
GetYaxis() -> SetLabelSize(fLab[1]);
1373 hPtTruth ->
GetYaxis() -> CenterTitle(fCnt);
1381 hPtDelta ->
SetTitle(sTitle.Data());
1382 hPtDelta -> SetTitleFont(fTxt);
1383 hPtDelta ->
GetXaxis() -> SetRangeUser(rDeltaRange[0], rDeltaRange[1]);
1385 hPtDelta ->
GetXaxis() -> SetTitleFont(fTxt);
1386 hPtDelta ->
GetXaxis() -> SetTitleSize(fTit[1]);
1387 hPtDelta ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
1388 hPtDelta ->
GetXaxis() -> SetLabelFont(fTxt);
1389 hPtDelta ->
GetXaxis() -> SetLabelSize(fLab[1]);
1390 hPtDelta ->
GetXaxis() -> CenterTitle(fCnt);
1392 hPtDelta ->
GetYaxis() -> SetTitleFont(fTxt);
1393 hPtDelta ->
GetYaxis() -> SetTitleSize(fTit[1]);
1394 hPtDelta ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
1395 hPtDelta ->
GetYaxis() -> SetLabelFont(fTxt);
1396 hPtDelta ->
GetYaxis() -> SetLabelSize(fLab[1]);
1397 hPtDelta ->
GetYaxis() -> CenterTitle(fCnt);
1405 hPtTrack ->
SetTitle(sTitle.Data());
1406 hPtTrack -> SetTitleFont(fTxt);
1407 hPtTrack ->
GetXaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
1409 hPtTrack ->
GetXaxis() -> SetTitleFont(fTxt);
1410 hPtTrack ->
GetXaxis() -> SetTitleSize(fTit[1]);
1411 hPtTrack ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
1412 hPtTrack ->
GetXaxis() -> SetLabelFont(fTxt);
1413 hPtTrack ->
GetXaxis() -> SetLabelSize(fLab[1]);
1414 hPtTrack ->
GetXaxis() -> CenterTitle(fCnt);
1416 hPtTrack ->
GetYaxis() -> SetTitleFont(fTxt);
1417 hPtTrack ->
GetYaxis() -> SetTitleSize(fTit[1]);
1418 hPtTrack ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
1419 hPtTrack ->
GetYaxis() -> SetLabelFont(fTxt);
1420 hPtTrack ->
GetYaxis() -> SetLabelSize(fLab[1]);
1421 hPtTrack ->
GetYaxis() -> CenterTitle(fCnt);
1429 hPtFrac ->
SetTitle(sTitle.Data());
1430 hPtFrac -> SetTitleFont(fTxt);
1431 hPtFrac ->
GetXaxis() -> SetRangeUser(rFracRange[0], rFracRange[1]);
1433 hPtFrac ->
GetXaxis() -> SetTitleFont(fTxt);
1434 hPtFrac ->
GetXaxis() -> SetTitleSize(fTit[1]);
1435 hPtFrac ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
1436 hPtFrac ->
GetXaxis() -> SetLabelFont(fTxt);
1437 hPtFrac ->
GetXaxis() -> SetLabelSize(fLab[1]);
1438 hPtFrac ->
GetXaxis() -> CenterTitle(fCnt);
1440 hPtFrac ->
GetYaxis() -> SetTitleFont(fTxt);
1441 hPtFrac ->
GetYaxis() -> SetTitleSize(fTit[1]);
1442 hPtFrac ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
1443 hPtFrac ->
GetYaxis() -> SetLabelFont(fTxt);
1444 hPtFrac ->
GetYaxis() -> SetLabelSize(fLab[1]);
1445 hPtFrac ->
GetYaxis() -> CenterTitle(fCnt);
1453 hPtTrkTru ->
SetTitle(sTitle.Data());
1454 hPtTrkTru -> SetTitleFont(fTxt);
1455 hPtTrkTru ->
GetXaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
1457 hPtTrkTru ->
GetXaxis() -> SetTitleFont(fTxt);
1458 hPtTrkTru ->
GetXaxis() -> SetTitleSize(fTit[1]);
1459 hPtTrkTru ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
1460 hPtTrkTru ->
GetXaxis() -> SetLabelFont(fTxt);
1461 hPtTrkTru ->
GetXaxis() -> SetLabelSize(fLab[1]);
1462 hPtTrkTru ->
GetXaxis() -> CenterTitle(fCnt);
1464 hPtTrkTru ->
GetYaxis() -> SetTitleFont(fTxt);
1465 hPtTrkTru ->
GetYaxis() -> SetTitleSize(fTit[1]);
1466 hPtTrkTru ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
1467 hPtTrkTru ->
GetYaxis() -> SetLabelFont(fTxt);
1468 hPtTrkTru ->
GetYaxis() -> SetLabelSize(fLab[1]);
1469 hPtTrkTru ->
GetYaxis() -> CenterTitle(fCnt);
1477 hPtDeltaVsFrac ->
SetTitle(sTitle.Data());
1478 hPtDeltaVsFrac -> SetTitleFont(fTxt);
1479 hPtDeltaVsFrac ->
GetXaxis() -> SetRangeUser(rFracRange[0], rFracRange[1]);
1481 hPtDeltaVsFrac ->
GetXaxis() -> SetTitleFont(fTxt);
1482 hPtDeltaVsFrac ->
GetXaxis() -> SetTitleSize(fTit[1]);
1483 hPtDeltaVsFrac ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
1484 hPtDeltaVsFrac ->
GetXaxis() -> SetLabelFont(fTxt);
1485 hPtDeltaVsFrac ->
GetXaxis() -> SetLabelSize(fLab[1]);
1486 hPtDeltaVsFrac ->
GetXaxis() -> CenterTitle(fCnt);
1487 hPtDeltaVsFrac ->
GetYaxis() -> SetRangeUser(rDeltaRange[0], rDeltaRange[1]);
1489 hPtDeltaVsFrac ->
GetYaxis() -> SetTitleFont(fTxt);
1490 hPtDeltaVsFrac ->
GetYaxis() -> SetTitleSize(fTit[1]);
1491 hPtDeltaVsFrac ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
1492 hPtDeltaVsFrac ->
GetYaxis() -> SetLabelFont(fTxt);
1493 hPtDeltaVsFrac ->
GetYaxis() -> SetLabelSize(fLab[1]);
1494 hPtDeltaVsFrac ->
GetYaxis() -> CenterTitle(fCnt);
1496 hPtDeltaVsFrac ->
GetZaxis() -> SetTitleFont(fTxt);
1497 hPtDeltaVsFrac ->
GetZaxis() -> SetTitleSize(fTit[1]);
1498 hPtDeltaVsFrac ->
GetZaxis() -> SetTitleOffset(fOffZ[1]);
1499 hPtDeltaVsFrac ->
GetZaxis() -> SetLabelFont(fTxt);
1500 hPtDeltaVsFrac ->
GetZaxis() -> SetLabelSize(fLab[1]);
1501 hPtDeltaVsFrac ->
GetZaxis() -> CenterTitle(fCnt);
1509 hPtDeltaVsTrue ->
SetTitle(sTitle.Data());
1510 hPtDeltaVsTrue -> SetTitleFont(fTxt);
1511 hPtDeltaVsTrue ->
GetXaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
1513 hPtDeltaVsTrue ->
GetXaxis() -> SetTitleFont(fTxt);
1514 hPtDeltaVsTrue ->
GetXaxis() -> SetTitleSize(fTit[1]);
1515 hPtDeltaVsTrue ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
1516 hPtDeltaVsTrue ->
GetXaxis() -> SetLabelFont(fTxt);
1517 hPtDeltaVsTrue ->
GetXaxis() -> SetLabelSize(fLab[1]);
1518 hPtDeltaVsTrue ->
GetXaxis() -> CenterTitle(fCnt);
1519 hPtDeltaVsTrue ->
GetYaxis() -> SetRangeUser(rDeltaRange[0], rDeltaRange[1]);
1521 hPtDeltaVsTrue ->
GetYaxis() -> SetTitleFont(fTxt);
1522 hPtDeltaVsTrue ->
GetYaxis() -> SetTitleSize(fTit[1]);
1523 hPtDeltaVsTrue ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
1524 hPtDeltaVsTrue ->
GetYaxis() -> SetLabelFont(fTxt);
1525 hPtDeltaVsTrue ->
GetYaxis() -> SetLabelSize(fLab[1]);
1526 hPtDeltaVsTrue ->
GetYaxis() -> CenterTitle(fCnt);
1528 hPtDeltaVsTrue ->
GetZaxis() -> SetTitleFont(fTxt);
1529 hPtDeltaVsTrue ->
GetZaxis() -> SetTitleSize(fTit[1]);
1530 hPtDeltaVsTrue ->
GetZaxis() -> SetTitleOffset(fOffZ[1]);
1531 hPtDeltaVsTrue ->
GetZaxis() -> SetLabelFont(fTxt);
1532 hPtDeltaVsTrue ->
GetZaxis() -> SetLabelSize(fLab[1]);
1533 hPtDeltaVsTrue ->
GetZaxis() -> CenterTitle(fCnt);
1541 hPtDeltaVsTrack ->
SetTitle(sTitle.Data());
1542 hPtDeltaVsTrack -> SetTitleFont(fTxt);
1543 hPtDeltaVsTrack ->
GetXaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
1545 hPtDeltaVsTrack ->
GetXaxis() -> SetTitleFont(fTxt);
1546 hPtDeltaVsTrack ->
GetXaxis() -> SetTitleSize(fTit[1]);
1547 hPtDeltaVsTrack ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
1548 hPtDeltaVsTrack ->
GetXaxis() -> SetLabelFont(fTxt);
1549 hPtDeltaVsTrack ->
GetXaxis() -> SetLabelSize(fLab[1]);
1550 hPtDeltaVsTrack ->
GetXaxis() -> CenterTitle(fCnt);
1551 hPtDeltaVsTrack ->
GetYaxis() -> SetRangeUser(rDeltaRange[0], rDeltaRange[1]);
1553 hPtDeltaVsTrack ->
GetYaxis() -> SetTitleFont(fTxt);
1554 hPtDeltaVsTrack ->
GetYaxis() -> SetTitleSize(fTit[1]);
1555 hPtDeltaVsTrack ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
1556 hPtDeltaVsTrack ->
GetYaxis() -> SetLabelFont(fTxt);
1557 hPtDeltaVsTrack ->
GetYaxis() -> SetLabelSize(fLab[1]);
1558 hPtDeltaVsTrack ->
GetYaxis() -> CenterTitle(fCnt);
1560 hPtDeltaVsTrack ->
GetZaxis() -> SetTitleFont(fTxt);
1561 hPtDeltaVsTrack ->
GetZaxis() -> SetTitleSize(fTit[1]);
1562 hPtDeltaVsTrack ->
GetZaxis() -> SetTitleOffset(fOffZ[1]);
1563 hPtDeltaVsTrack ->
GetZaxis() -> SetLabelFont(fTxt);
1564 hPtDeltaVsTrack ->
GetZaxis() -> SetLabelSize(fLab[1]);
1565 hPtDeltaVsTrack ->
GetZaxis() -> CenterTitle(fCnt);
1573 hPtTrueVsTrack ->
SetTitle(sTitle.Data());
1574 hPtTrueVsTrack -> SetTitleFont(fTxt);
1575 hPtTrueVsTrack ->
GetXaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
1577 hPtTrueVsTrack ->
GetXaxis() -> SetTitleFont(fTxt);
1578 hPtTrueVsTrack ->
GetXaxis() -> SetTitleSize(fTit[1]);
1579 hPtTrueVsTrack ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
1580 hPtTrueVsTrack ->
GetXaxis() -> SetLabelFont(fTxt);
1581 hPtTrueVsTrack ->
GetXaxis() -> SetLabelSize(fLab[1]);
1582 hPtTrueVsTrack ->
GetXaxis() -> CenterTitle(fCnt);
1583 hPtTrueVsTrack ->
GetYaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
1585 hPtTrueVsTrack ->
GetYaxis() -> SetTitleFont(fTxt);
1586 hPtTrueVsTrack ->
GetYaxis() -> SetTitleSize(fTit[1]);
1587 hPtTrueVsTrack ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
1588 hPtTrueVsTrack ->
GetYaxis() -> SetLabelFont(fTxt);
1589 hPtTrueVsTrack ->
GetYaxis() -> SetLabelSize(fLab[1]);
1590 hPtTrueVsTrack ->
GetYaxis() -> CenterTitle(fCnt);
1592 hPtTrueVsTrack ->
GetZaxis() -> SetTitleFont(fTxt);
1593 hPtTrueVsTrack ->
GetZaxis() -> SetTitleSize(fTit[1]);
1594 hPtTrueVsTrack ->
GetZaxis() -> SetTitleOffset(fOffZ[1]);
1595 hPtTrueVsTrack ->
GetZaxis() -> SetLabelFont(fTxt);
1596 hPtTrueVsTrack ->
GetZaxis() -> SetLabelSize(fLab[1]);
1597 hPtTrueVsTrack ->
GetZaxis() -> CenterTitle(fCnt);
1607 grMuProj ->
SetTitle(sTitle.Data());
1608 grMuProj ->
GetXaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
1610 grMuProj ->
GetXaxis() -> SetTitleFont(fTxt);
1611 grMuProj ->
GetXaxis() -> SetTitleSize(fTit[1]);
1612 grMuProj ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
1613 grMuProj ->
GetXaxis() -> SetLabelFont(fTxt);
1614 grMuProj ->
GetXaxis() -> SetLabelSize(fLab[1]);
1615 grMuProj ->
GetXaxis() -> CenterTitle(fCnt);
1617 grMuProj ->
GetYaxis() -> SetTitleFont(fTxt);
1618 grMuProj ->
GetYaxis() -> SetTitleSize(fTit[1]);
1619 grMuProj ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
1620 grMuProj ->
GetYaxis() -> SetLabelFont(fTxt);
1621 grMuProj ->
GetYaxis() -> SetLabelSize(fLab[1]);
1622 grMuProj ->
GetYaxis() -> CenterTitle(fCnt);
1630 grSigProj ->
SetTitle(sTitle.Data());
1631 grSigProj ->
GetXaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
1633 grSigProj ->
GetXaxis() -> SetTitleFont(fTxt);
1634 grSigProj ->
GetXaxis() -> SetTitleSize(fTit[1]);
1635 grSigProj ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
1636 grSigProj ->
GetXaxis() -> SetLabelFont(fTxt);
1637 grSigProj ->
GetXaxis() -> SetLabelSize(fLab[1]);
1638 grSigProj ->
GetXaxis() -> CenterTitle(fCnt);
1640 grSigProj ->
GetYaxis() -> SetTitleFont(fTxt);
1641 grSigProj ->
GetYaxis() -> SetTitleSize(fTit[1]);
1642 grSigProj ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
1643 grSigProj ->
GetYaxis() -> SetLabelFont(fTxt);
1644 grSigProj ->
GetYaxis() -> SetLabelSize(fLab[1]);
1645 grSigProj ->
GetYaxis() -> CenterTitle(fCnt);
1646 for (Ssiz_t iProj = 0; iProj <
NProj; iProj++) {
1654 hPtDeltaProj[iProj] ->
SetTitle(sTitle.Data());
1655 hPtDeltaProj[iProj] -> SetTitleFont(fTxt);
1656 hPtDeltaProj[iProj] ->
GetXaxis() -> SetRangeUser(rDeltaRange[0], rDeltaRange[1]);
1658 hPtDeltaProj[iProj] ->
GetXaxis() -> SetTitleFont(fTxt);
1659 hPtDeltaProj[iProj] ->
GetXaxis() -> SetTitleSize(fTit[1]);
1660 hPtDeltaProj[iProj] ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
1661 hPtDeltaProj[iProj] ->
GetXaxis() -> SetLabelFont(fTxt);
1662 hPtDeltaProj[iProj] ->
GetXaxis() -> SetLabelSize(fLab[1]);
1663 hPtDeltaProj[iProj] ->
GetXaxis() -> CenterTitle(fCnt);
1665 hPtDeltaProj[iProj] ->
GetYaxis() -> SetTitleFont(fTxt);
1666 hPtDeltaProj[iProj] ->
GetYaxis() -> SetTitleSize(fTit[1]);
1667 hPtDeltaProj[iProj] ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
1668 hPtDeltaProj[iProj] ->
GetYaxis() -> SetLabelFont(fTxt);
1669 hPtDeltaProj[iProj] ->
GetYaxis() -> SetLabelSize(fLab[1]);
1670 hPtDeltaProj[iProj] ->
GetYaxis() -> CenterTitle(fCnt);
1672 for (Ssiz_t iSig = 0; iSig <
NSigCuts; iSig++) {
1680 grMuHiProj[iSig] ->
SetTitle(sTitle.Data());
1681 grMuHiProj[iSig] ->
GetXaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
1683 grMuHiProj[iSig] ->
GetXaxis() -> SetTitleFont(fTxt);
1684 grMuHiProj[iSig] ->
GetXaxis() -> SetTitleSize(fTit[1]);
1685 grMuHiProj[iSig] ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
1686 grMuHiProj[iSig] ->
GetXaxis() -> SetLabelFont(fTxt);
1687 grMuHiProj[iSig] ->
GetXaxis() -> SetLabelSize(fLab[1]);
1688 grMuHiProj[iSig] ->
GetXaxis() -> CenterTitle(fCnt);
1690 grMuHiProj[iSig] ->
GetYaxis() -> SetTitleFont(fTxt);
1691 grMuHiProj[iSig] ->
GetYaxis() -> SetTitleSize(fTit[1]);
1692 grMuHiProj[iSig] ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
1693 grMuHiProj[iSig] ->
GetYaxis() -> SetLabelFont(fTxt);
1694 grMuHiProj[iSig] ->
GetYaxis() -> SetLabelSize(fLab[1]);
1695 grMuHiProj[iSig] ->
GetYaxis() -> CenterTitle(fCnt);
1703 grMuLoProj[iSig] ->
SetTitle(sTitle.Data());
1704 grMuLoProj[iSig] ->
GetXaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
1706 grMuLoProj[iSig] ->
GetXaxis() -> SetTitleFont(fTxt);
1707 grMuLoProj[iSig] ->
GetXaxis() -> SetTitleSize(fTit[1]);
1708 grMuLoProj[iSig] ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
1709 grMuLoProj[iSig] ->
GetXaxis() -> SetLabelFont(fTxt);
1710 grMuLoProj[iSig] ->
GetXaxis() -> SetLabelSize(fLab[1]);
1711 grMuLoProj[iSig] ->
GetXaxis() -> CenterTitle(fCnt);
1713 grMuLoProj[iSig] ->
GetYaxis() -> SetTitleFont(fTxt);
1714 grMuLoProj[iSig] ->
GetYaxis() -> SetTitleSize(fTit[1]);
1715 grMuLoProj[iSig] ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
1716 grMuLoProj[iSig] ->
GetYaxis() -> SetLabelFont(fTxt);
1717 grMuLoProj[iSig] ->
GetYaxis() -> SetLabelSize(fLab[1]);
1718 grMuLoProj[iSig] ->
GetYaxis() -> CenterTitle(fCnt);
1722 for (Ssiz_t iCut = 0; iCut <
NDPtCuts; iCut++) {
1730 hEffCut[iCut] ->
SetTitle(sTitle.Data());
1731 hEffCut[iCut] -> SetTitleFont(fTxt);
1732 hEffCut[iCut] ->
GetXaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
1734 hEffCut[iCut] ->
GetXaxis() -> SetTitleFont(fTxt);
1735 hEffCut[iCut] ->
GetXaxis() -> SetTitleSize(fTit[0]);
1736 hEffCut[iCut] ->
GetXaxis() -> SetTitleOffset(fOffX[0]);
1737 hEffCut[iCut] ->
GetXaxis() -> SetLabelFont(fTxt);
1738 hEffCut[iCut] ->
GetXaxis() -> SetLabelSize(fLab[0]);
1739 hEffCut[iCut] ->
GetXaxis() -> CenterTitle(fCnt);
1741 hEffCut[iCut] ->
GetYaxis() -> SetTitleFont(fTxt);
1742 hEffCut[iCut] ->
GetYaxis() -> SetTitleSize(fTit[0]);
1743 hEffCut[iCut] ->
GetYaxis() -> SetTitleOffset(fOffY[0]);
1744 hEffCut[iCut] ->
GetYaxis() -> SetLabelFont(fTxt);
1745 hEffCut[iCut] ->
GetYaxis() -> SetLabelSize(fLab[0]);
1746 hEffCut[iCut] ->
GetYaxis() -> CenterTitle(fCnt);
1754 hPtDeltaCut[iCut] ->
SetTitle(sTitle.Data());
1755 hPtDeltaCut[iCut] -> SetTitleFont(fTxt);
1756 hPtDeltaCut[iCut] ->
GetXaxis() -> SetRangeUser(rDeltaRange[0], rDeltaRange[1]);
1758 hPtDeltaCut[iCut] ->
GetXaxis() -> SetTitleFont(fTxt);
1759 hPtDeltaCut[iCut] ->
GetXaxis() -> SetTitleSize(fTit[1]);
1760 hPtDeltaCut[iCut] ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
1761 hPtDeltaCut[iCut] ->
GetXaxis() -> SetLabelFont(fTxt);
1762 hPtDeltaCut[iCut] ->
GetXaxis() -> SetLabelSize(fLab[1]);
1763 hPtDeltaCut[iCut] ->
GetXaxis() -> CenterTitle(fCnt);
1765 hPtDeltaCut[iCut] ->
GetYaxis() -> SetTitleFont(fTxt);
1766 hPtDeltaCut[iCut] ->
GetYaxis() -> SetTitleSize(fTit[1]);
1767 hPtDeltaCut[iCut] ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
1768 hPtDeltaCut[iCut] ->
GetYaxis() -> SetLabelFont(fTxt);
1769 hPtDeltaCut[iCut] ->
GetYaxis() -> SetLabelSize(fLab[1]);
1770 hPtDeltaCut[iCut] ->
GetYaxis() -> CenterTitle(fCnt);
1778 hPtTrackCut[iCut] ->
SetTitle(sTitle.Data());
1779 hPtTrackCut[iCut] -> SetTitleFont(fTxt);
1780 hPtTrackCut[iCut] ->
GetXaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
1782 hPtTrackCut[iCut] ->
GetXaxis() -> SetTitleFont(fTxt);
1783 hPtTrackCut[iCut] ->
GetXaxis() -> SetTitleSize(fTit[1]);
1784 hPtTrackCut[iCut] ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
1785 hPtTrackCut[iCut] ->
GetXaxis() -> SetLabelFont(fTxt);
1786 hPtTrackCut[iCut] ->
GetXaxis() -> SetLabelSize(fLab[1]);
1787 hPtTrackCut[iCut] ->
GetXaxis() -> CenterTitle(fCnt);
1789 hPtTrackCut[iCut] ->
GetYaxis() -> SetTitleFont(fTxt);
1790 hPtTrackCut[iCut] ->
GetYaxis() -> SetTitleSize(fTit[1]);
1791 hPtTrackCut[iCut] ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
1792 hPtTrackCut[iCut] ->
GetYaxis() -> SetLabelFont(fTxt);
1793 hPtTrackCut[iCut] ->
GetYaxis() -> SetLabelSize(fLab[1]);
1794 hPtTrackCut[iCut] ->
GetYaxis() -> CenterTitle(fCnt);
1802 hPtFracCut[iCut] ->
SetTitle(sTitle.Data());
1803 hPtFracCut[iCut] -> SetTitleFont(fTxt);
1804 hPtFracCut[iCut] ->
GetXaxis() -> SetRangeUser(rFracRange[0], rFracRange[1]);
1806 hPtFracCut[iCut] ->
GetXaxis() -> SetTitleFont(fTxt);
1807 hPtFracCut[iCut] ->
GetXaxis() -> SetTitleSize(fTit[1]);
1808 hPtFracCut[iCut] ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
1809 hPtFracCut[iCut] ->
GetXaxis() -> SetLabelFont(fTxt);
1810 hPtFracCut[iCut] ->
GetXaxis() -> SetLabelSize(fLab[1]);
1811 hPtFracCut[iCut] ->
GetXaxis() -> CenterTitle(fCnt);
1813 hPtFracCut[iCut] ->
GetYaxis() -> SetTitleFont(fTxt);
1814 hPtFracCut[iCut] ->
GetYaxis() -> SetTitleSize(fTit[1]);
1815 hPtFracCut[iCut] ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
1816 hPtFracCut[iCut] ->
GetYaxis() -> SetLabelFont(fTxt);
1817 hPtFracCut[iCut] ->
GetYaxis() -> SetLabelSize(fLab[1]);
1818 hPtFracCut[iCut] ->
GetYaxis() -> CenterTitle(fCnt);
1826 hPtTrkTruCut[iCut] ->
SetTitle(sTitle.Data());
1827 hPtTrkTruCut[iCut] -> SetTitleFont(fTxt);
1828 hPtTrkTruCut[iCut] ->
GetXaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
1830 hPtTrkTruCut[iCut] ->
GetXaxis() -> SetTitleFont(fTxt);
1831 hPtTrkTruCut[iCut] ->
GetXaxis() -> SetTitleSize(fTit[1]);
1832 hPtTrkTruCut[iCut] ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
1833 hPtTrkTruCut[iCut] ->
GetXaxis() -> SetLabelFont(fTxt);
1834 hPtTrkTruCut[iCut] ->
GetXaxis() -> SetLabelSize(fLab[1]);
1835 hPtTrkTruCut[iCut] ->
GetXaxis() -> CenterTitle(fCnt);
1837 hPtTrkTruCut[iCut] ->
GetYaxis() -> SetTitleFont(fTxt);
1838 hPtTrkTruCut[iCut] ->
GetYaxis() -> SetTitleSize(fTit[1]);
1839 hPtTrkTruCut[iCut] ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
1840 hPtTrkTruCut[iCut] ->
GetYaxis() -> SetLabelFont(fTxt);
1841 hPtTrkTruCut[iCut] ->
GetYaxis() -> SetLabelSize(fLab[1]);
1842 hPtTrkTruCut[iCut] ->
GetYaxis() -> CenterTitle(fCnt);
1845 hPtDeltaVsFracCut[iCut] ->
SetFillColor(fColCut[iCut]);
1847 hPtDeltaVsFracCut[iCut] ->
SetLineColor(fColCut[iCut]);
1850 hPtDeltaVsFracCut[iCut] ->
SetTitle(sTitle.Data());
1851 hPtDeltaVsFracCut[iCut] -> SetTitleFont(fTxt);
1852 hPtDeltaVsFracCut[iCut] ->
GetXaxis() -> SetRangeUser(rFracRange[0], rFracRange[1]);
1854 hPtDeltaVsFracCut[iCut] ->
GetXaxis() -> SetTitleFont(fTxt);
1855 hPtDeltaVsFracCut[iCut] ->
GetXaxis() -> SetTitleSize(fTit[1]);
1856 hPtDeltaVsFracCut[iCut] ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
1857 hPtDeltaVsFracCut[iCut] ->
GetXaxis() -> SetLabelFont(fTxt);
1858 hPtDeltaVsFracCut[iCut] ->
GetXaxis() -> SetLabelSize(fLab[1]);
1859 hPtDeltaVsFracCut[iCut] ->
GetXaxis() -> CenterTitle(fCnt);
1860 hPtDeltaVsFracCut[iCut] ->
GetYaxis() -> SetRangeUser(rDeltaRange[0], rDeltaRange[1]);
1862 hPtDeltaVsFracCut[iCut] ->
GetYaxis() -> SetTitleFont(fTxt);
1863 hPtDeltaVsFracCut[iCut] ->
GetYaxis() -> SetTitleSize(fTit[1]);
1864 hPtDeltaVsFracCut[iCut] ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
1865 hPtDeltaVsFracCut[iCut] ->
GetYaxis() -> SetLabelFont(fTxt);
1866 hPtDeltaVsFracCut[iCut] ->
GetYaxis() -> SetLabelSize(fLab[1]);
1867 hPtDeltaVsFracCut[iCut] ->
GetYaxis() -> CenterTitle(fCnt);
1869 hPtDeltaVsFracCut[iCut] ->
GetZaxis() -> SetTitleFont(fTxt);
1870 hPtDeltaVsFracCut[iCut] ->
GetZaxis() -> SetTitleSize(fTit[1]);
1871 hPtDeltaVsFracCut[iCut] ->
GetZaxis() -> SetTitleOffset(fOffZ[1]);
1872 hPtDeltaVsFracCut[iCut] ->
GetZaxis() -> SetLabelFont(fTxt);
1873 hPtDeltaVsFracCut[iCut] ->
GetZaxis() -> SetLabelSize(fLab[1]);
1874 hPtDeltaVsFracCut[iCut] ->
GetZaxis() -> CenterTitle(fCnt);
1882 hPtDeltaVsTrueCut[iCut] ->
SetTitle(sTitle.Data());
1883 hPtDeltaVsTrueCut[iCut] -> SetTitleFont(fTxt);
1884 hPtDeltaVsTrueCut[iCut] ->
GetXaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
1886 hPtDeltaVsTrueCut[iCut] ->
GetXaxis() -> SetTitleFont(fTxt);
1887 hPtDeltaVsTrueCut[iCut] ->
GetXaxis() -> SetTitleSize(fTit[1]);
1888 hPtDeltaVsTrueCut[iCut] ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
1889 hPtDeltaVsTrueCut[iCut] ->
GetXaxis() -> SetLabelFont(fTxt);
1890 hPtDeltaVsTrueCut[iCut] ->
GetXaxis() -> SetLabelSize(fLab[1]);
1891 hPtDeltaVsTrueCut[iCut] ->
GetXaxis() -> CenterTitle(fCnt);
1892 hPtDeltaVsTrueCut[iCut] ->
GetYaxis() -> SetRangeUser(rDeltaRange[0], rDeltaRange[1]);
1894 hPtDeltaVsTrueCut[iCut] ->
GetYaxis() -> SetTitleFont(fTxt);
1895 hPtDeltaVsTrueCut[iCut] ->
GetYaxis() -> SetTitleSize(fTit[1]);
1896 hPtDeltaVsTrueCut[iCut] ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
1897 hPtDeltaVsTrueCut[iCut] ->
GetYaxis() -> SetLabelFont(fTxt);
1898 hPtDeltaVsTrueCut[iCut] ->
GetYaxis() -> SetLabelSize(fLab[1]);
1899 hPtDeltaVsTrueCut[iCut] ->
GetYaxis() -> CenterTitle(fCnt);
1901 hPtDeltaVsTrueCut[iCut] ->
GetZaxis() -> SetTitleFont(fTxt);
1902 hPtDeltaVsTrueCut[iCut] ->
GetZaxis() -> SetTitleSize(fTit[1]);
1903 hPtDeltaVsTrueCut[iCut] ->
GetZaxis() -> SetTitleOffset(fOffZ[1]);
1904 hPtDeltaVsTrueCut[iCut] ->
GetZaxis() -> SetLabelFont(fTxt);
1905 hPtDeltaVsTrueCut[iCut] ->
GetZaxis() -> SetLabelSize(fLab[1]);
1906 hPtDeltaVsTrueCut[iCut] ->
GetZaxis() -> CenterTitle(fCnt);
1914 hPtDeltaVsTrackCut[iCut] ->
SetTitle(sTitle.Data());
1915 hPtDeltaVsTrackCut[iCut] -> SetTitleFont(fTxt);
1916 hPtDeltaVsTrackCut[iCut] ->
GetXaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
1918 hPtDeltaVsTrackCut[iCut] ->
GetXaxis() -> SetTitleFont(fTxt);
1919 hPtDeltaVsTrackCut[iCut] ->
GetXaxis() -> SetTitleSize(fTit[1]);
1920 hPtDeltaVsTrackCut[iCut] ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
1921 hPtDeltaVsTrackCut[iCut] ->
GetXaxis() -> SetLabelFont(fTxt);
1922 hPtDeltaVsTrackCut[iCut] ->
GetXaxis() -> SetLabelSize(fLab[1]);
1923 hPtDeltaVsTrackCut[iCut] ->
GetXaxis() -> CenterTitle(fCnt);
1924 hPtDeltaVsTrackCut[iCut] ->
GetYaxis() -> SetRangeUser(rDeltaRange[0], rDeltaRange[1]);
1926 hPtDeltaVsTrackCut[iCut] ->
GetYaxis() -> SetTitleFont(fTxt);
1927 hPtDeltaVsTrackCut[iCut] ->
GetYaxis() -> SetTitleSize(fTit[1]);
1928 hPtDeltaVsTrackCut[iCut] ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
1929 hPtDeltaVsTrackCut[iCut] ->
GetYaxis() -> SetLabelFont(fTxt);
1930 hPtDeltaVsTrackCut[iCut] ->
GetYaxis() -> SetLabelSize(fLab[1]);
1931 hPtDeltaVsTrackCut[iCut] ->
GetYaxis() -> CenterTitle(fCnt);
1933 hPtDeltaVsTrackCut[iCut] ->
GetZaxis() -> SetTitleFont(fTxt);
1934 hPtDeltaVsTrackCut[iCut] ->
GetZaxis() -> SetTitleSize(fTit[1]);
1935 hPtDeltaVsTrackCut[iCut] ->
GetZaxis() -> SetTitleOffset(fOffZ[1]);
1936 hPtDeltaVsTrackCut[iCut] ->
GetZaxis() -> SetLabelFont(fTxt);
1937 hPtDeltaVsTrackCut[iCut] ->
GetZaxis() -> SetLabelSize(fLab[1]);
1938 hPtDeltaVsTrackCut[iCut] ->
GetZaxis() -> CenterTitle(fCnt);
1941 hPtTrueVsTrackCut[iCut] ->
SetFillColor(fColCut[iCut]);
1943 hPtTrueVsTrackCut[iCut] ->
SetLineColor(fColCut[iCut]);
1946 hPtTrueVsTrackCut[iCut] ->
SetTitle(sTitle.Data());
1947 hPtTrueVsTrackCut[iCut] -> SetTitleFont(fTxt);
1948 hPtTrueVsTrackCut[iCut] ->
GetXaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
1950 hPtTrueVsTrackCut[iCut] ->
GetXaxis() -> SetTitleFont(fTxt);
1951 hPtTrueVsTrackCut[iCut] ->
GetXaxis() -> SetTitleSize(fTit[1]);
1952 hPtTrueVsTrackCut[iCut] ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
1953 hPtTrueVsTrackCut[iCut] ->
GetXaxis() -> SetLabelFont(fTxt);
1954 hPtTrueVsTrackCut[iCut] ->
GetXaxis() -> SetLabelSize(fLab[1]);
1955 hPtTrueVsTrackCut[iCut] ->
GetXaxis() -> CenterTitle(fCnt);
1956 hPtTrueVsTrackCut[iCut] ->
GetYaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
1958 hPtTrueVsTrackCut[iCut] ->
GetYaxis() -> SetTitleFont(fTxt);
1959 hPtTrueVsTrackCut[iCut] ->
GetYaxis() -> SetTitleSize(fTit[1]);
1960 hPtTrueVsTrackCut[iCut] ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
1961 hPtTrueVsTrackCut[iCut] ->
GetYaxis() -> SetLabelFont(fTxt);
1962 hPtTrueVsTrackCut[iCut] ->
GetYaxis() -> SetLabelSize(fLab[1]);
1963 hPtTrueVsTrackCut[iCut] ->
GetYaxis() -> CenterTitle(fCnt);
1965 hPtTrueVsTrackCut[iCut] ->
GetZaxis() -> SetTitleFont(fTxt);
1966 hPtTrueVsTrackCut[iCut] ->
GetZaxis() -> SetTitleSize(fTit[1]);
1967 hPtTrueVsTrackCut[iCut] ->
GetZaxis() -> SetTitleOffset(fOffZ[1]);
1968 hPtTrueVsTrackCut[iCut] ->
GetZaxis() -> SetLabelFont(fTxt);
1969 hPtTrueVsTrackCut[iCut] ->
GetZaxis() -> SetLabelSize(fLab[1]);
1970 hPtTrueVsTrackCut[iCut] ->
GetZaxis() -> CenterTitle(fCnt);
1974 for (Ssiz_t iSig = 0; iSig <
NSigCuts; iSig++) {
1982 hEffSig[iSig] ->
SetTitle(sTitle.Data());
1983 hEffSig[iSig] -> SetTitleFont(fTxt);
1984 hEffSig[iSig] ->
GetXaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
1986 hEffSig[iSig] ->
GetXaxis() -> SetTitleFont(fTxt);
1987 hEffSig[iSig] ->
GetXaxis() -> SetTitleSize(fTit[0]);
1988 hEffSig[iSig] ->
GetXaxis() -> SetTitleOffset(fOffX[0]);
1989 hEffSig[iSig] ->
GetXaxis() -> SetLabelFont(fTxt);
1990 hEffSig[iSig] ->
GetXaxis() -> SetLabelSize(fLab[0]);
1991 hEffSig[iSig] ->
GetXaxis() -> CenterTitle(fCnt);
1993 hEffSig[iSig] ->
GetYaxis() -> SetTitleFont(fTxt);
1994 hEffSig[iSig] ->
GetYaxis() -> SetTitleSize(fTit[0]);
1995 hEffSig[iSig] ->
GetYaxis() -> SetTitleOffset(fOffY[0]);
1996 hEffSig[iSig] ->
GetYaxis() -> SetLabelFont(fTxt);
1997 hEffSig[iSig] ->
GetYaxis() -> SetLabelSize(fLab[0]);
1998 hEffSig[iSig] ->
GetYaxis() -> CenterTitle(fCnt);
2006 hPtDeltaSig[iSig] ->
SetTitle(sTitle.Data());
2007 hPtDeltaSig[iSig] -> SetTitleFont(fTxt);
2008 hPtDeltaSig[iSig] ->
GetXaxis() -> SetRangeUser(rDeltaRange[0], rDeltaRange[1]);
2010 hPtDeltaSig[iSig] ->
GetXaxis() -> SetTitleFont(fTxt);
2011 hPtDeltaSig[iSig] ->
GetXaxis() -> SetTitleSize(fTit[1]);
2012 hPtDeltaSig[iSig] ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
2013 hPtDeltaSig[iSig] ->
GetXaxis() -> SetLabelFont(fTxt);
2014 hPtDeltaSig[iSig] ->
GetXaxis() -> SetLabelSize(fLab[1]);
2015 hPtDeltaSig[iSig] ->
GetXaxis() -> CenterTitle(fCnt);
2017 hPtDeltaSig[iSig] ->
GetYaxis() -> SetTitleFont(fTxt);
2018 hPtDeltaSig[iSig] ->
GetYaxis() -> SetTitleSize(fTit[1]);
2019 hPtDeltaSig[iSig] ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
2020 hPtDeltaSig[iSig] ->
GetYaxis() -> SetLabelFont(fTxt);
2021 hPtDeltaSig[iSig] ->
GetYaxis() -> SetLabelSize(fLab[1]);
2022 hPtDeltaSig[iSig] ->
GetYaxis() -> CenterTitle(fCnt);
2030 hPtTrackSig[iSig] ->
SetTitle(sTitle.Data());
2031 hPtTrackSig[iSig] -> SetTitleFont(fTxt);
2032 hPtTrackSig[iSig] ->
GetXaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
2034 hPtTrackSig[iSig] ->
GetXaxis() -> SetTitleFont(fTxt);
2035 hPtTrackSig[iSig] ->
GetXaxis() -> SetTitleSize(fTit[1]);
2036 hPtTrackSig[iSig] ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
2037 hPtTrackSig[iSig] ->
GetXaxis() -> SetLabelFont(fTxt);
2038 hPtTrackSig[iSig] ->
GetXaxis() -> SetLabelSize(fLab[1]);
2039 hPtTrackSig[iSig] ->
GetXaxis() -> CenterTitle(fCnt);
2041 hPtTrackSig[iSig] ->
GetYaxis() -> SetTitleFont(fTxt);
2042 hPtTrackSig[iSig] ->
GetYaxis() -> SetTitleSize(fTit[1]);
2043 hPtTrackSig[iSig] ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
2044 hPtTrackSig[iSig] ->
GetYaxis() -> SetLabelFont(fTxt);
2045 hPtTrackSig[iSig] ->
GetYaxis() -> SetLabelSize(fLab[1]);
2046 hPtTrackSig[iSig] ->
GetYaxis() -> CenterTitle(fCnt);
2054 hPtFracSig[iSig] ->
SetTitle(sTitle.Data());
2055 hPtFracSig[iSig] -> SetTitleFont(fTxt);
2056 hPtFracSig[iSig] ->
GetXaxis() -> SetRangeUser(rFracRange[0], rFracRange[1]);
2058 hPtFracSig[iSig] ->
GetXaxis() -> SetTitleFont(fTxt);
2059 hPtFracSig[iSig] ->
GetXaxis() -> SetTitleSize(fTit[1]);
2060 hPtFracSig[iSig] ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
2061 hPtFracSig[iSig] ->
GetXaxis() -> SetLabelFont(fTxt);
2062 hPtFracSig[iSig] ->
GetXaxis() -> SetLabelSize(fLab[1]);
2063 hPtFracSig[iSig] ->
GetXaxis() -> CenterTitle(fCnt);
2065 hPtFracSig[iSig] ->
GetYaxis() -> SetTitleFont(fTxt);
2066 hPtFracSig[iSig] ->
GetYaxis() -> SetTitleSize(fTit[1]);
2067 hPtFracSig[iSig] ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
2068 hPtFracSig[iSig] ->
GetYaxis() -> SetLabelFont(fTxt);
2069 hPtFracSig[iSig] ->
GetYaxis() -> SetLabelSize(fLab[1]);
2070 hPtFracSig[iSig] ->
GetYaxis() -> CenterTitle(fCnt);
2078 hPtTrkTruSig[iSig] ->
SetTitle(sTitle.Data());
2079 hPtTrkTruSig[iSig] -> SetTitleFont(fTxt);
2080 hPtTrkTruSig[iSig] ->
GetXaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
2082 hPtTrkTruSig[iSig] ->
GetXaxis() -> SetTitleFont(fTxt);
2083 hPtTrkTruSig[iSig] ->
GetXaxis() -> SetTitleSize(fTit[1]);
2084 hPtTrkTruSig[iSig] ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
2085 hPtTrkTruSig[iSig] ->
GetXaxis() -> SetLabelFont(fTxt);
2086 hPtTrkTruSig[iSig] ->
GetXaxis() -> SetLabelSize(fLab[1]);
2087 hPtTrkTruSig[iSig] ->
GetXaxis() -> CenterTitle(fCnt);
2089 hPtTrkTruSig[iSig] ->
GetYaxis() -> SetTitleFont(fTxt);
2090 hPtTrkTruSig[iSig] ->
GetYaxis() -> SetTitleSize(fTit[1]);
2091 hPtTrkTruSig[iSig] ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
2092 hPtTrkTruSig[iSig] ->
GetYaxis() -> SetLabelFont(fTxt);
2093 hPtTrkTruSig[iSig] ->
GetYaxis() -> SetLabelSize(fLab[1]);
2094 hPtTrkTruSig[iSig] ->
GetYaxis() -> CenterTitle(fCnt);
2097 hPtDeltaVsFracSig[iSig] ->
SetFillColor(fColSig[iSig]);
2099 hPtDeltaVsFracSig[iSig] ->
SetLineColor(fColSig[iSig]);
2102 hPtDeltaVsFracSig[iSig] ->
SetTitle(sTitle.Data());
2103 hPtDeltaVsFracSig[iSig] -> SetTitleFont(fTxt);
2104 hPtDeltaVsFracSig[iSig] ->
GetXaxis() -> SetRangeUser(rFracRange[0], rFracRange[1]);
2106 hPtDeltaVsFracSig[iSig] ->
GetXaxis() -> SetTitleFont(fTxt);
2107 hPtDeltaVsFracSig[iSig] ->
GetXaxis() -> SetTitleSize(fTit[1]);
2108 hPtDeltaVsFracSig[iSig] ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
2109 hPtDeltaVsFracSig[iSig] ->
GetXaxis() -> SetLabelFont(fTxt);
2110 hPtDeltaVsFracSig[iSig] ->
GetXaxis() -> SetLabelSize(fLab[1]);
2111 hPtDeltaVsFracSig[iSig] ->
GetXaxis() -> CenterTitle(fCnt);
2112 hPtDeltaVsFracSig[iSig] ->
GetYaxis() -> SetRangeUser(rDeltaRange[0], rDeltaRange[1]);
2114 hPtDeltaVsFracSig[iSig] ->
GetYaxis() -> SetTitleFont(fTxt);
2115 hPtDeltaVsFracSig[iSig] ->
GetYaxis() -> SetTitleSize(fTit[1]);
2116 hPtDeltaVsFracSig[iSig] ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
2117 hPtDeltaVsFracSig[iSig] ->
GetYaxis() -> SetLabelFont(fTxt);
2118 hPtDeltaVsFracSig[iSig] ->
GetYaxis() -> SetLabelSize(fLab[1]);
2119 hPtDeltaVsFracSig[iSig] ->
GetYaxis() -> CenterTitle(fCnt);
2121 hPtDeltaVsFracSig[iSig] ->
GetZaxis() -> SetTitleFont(fTxt);
2122 hPtDeltaVsFracSig[iSig] ->
GetZaxis() -> SetTitleSize(fTit[1]);
2123 hPtDeltaVsFracSig[iSig] ->
GetZaxis() -> SetTitleOffset(fOffZ[1]);
2124 hPtDeltaVsFracSig[iSig] ->
GetZaxis() -> SetLabelFont(fTxt);
2125 hPtDeltaVsFracSig[iSig] ->
GetZaxis() -> SetLabelSize(fLab[1]);
2126 hPtDeltaVsFracSig[iSig] ->
GetZaxis() -> CenterTitle(fCnt);
2134 hPtDeltaVsTrueSig[iSig] ->
SetTitle(sTitle.Data());
2135 hPtDeltaVsTrueSig[iSig] -> SetTitleFont(fTxt);
2136 hPtDeltaVsTrueSig[iSig] ->
GetXaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
2138 hPtDeltaVsTrueSig[iSig] ->
GetXaxis() -> SetTitleFont(fTxt);
2139 hPtDeltaVsTrueSig[iSig] ->
GetXaxis() -> SetTitleSize(fTit[1]);
2140 hPtDeltaVsTrueSig[iSig] ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
2141 hPtDeltaVsTrueSig[iSig] ->
GetXaxis() -> SetLabelFont(fTxt);
2142 hPtDeltaVsTrueSig[iSig] ->
GetXaxis() -> SetLabelSize(fLab[1]);
2143 hPtDeltaVsTrueSig[iSig] ->
GetXaxis() -> CenterTitle(fCnt);
2144 hPtDeltaVsTrueSig[iSig] ->
GetYaxis() -> SetRangeUser(rDeltaRange[0], rDeltaRange[1]);
2146 hPtDeltaVsTrueSig[iSig] ->
GetYaxis() -> SetTitleFont(fTxt);
2147 hPtDeltaVsTrueSig[iSig] ->
GetYaxis() -> SetTitleSize(fTit[1]);
2148 hPtDeltaVsTrueSig[iSig] ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
2149 hPtDeltaVsTrueSig[iSig] ->
GetYaxis() -> SetLabelFont(fTxt);
2150 hPtDeltaVsTrueSig[iSig] ->
GetYaxis() -> SetLabelSize(fLab[1]);
2151 hPtDeltaVsTrueSig[iSig] ->
GetYaxis() -> CenterTitle(fCnt);
2153 hPtDeltaVsTrueSig[iSig] ->
GetZaxis() -> SetTitleFont(fTxt);
2154 hPtDeltaVsTrueSig[iSig] ->
GetZaxis() -> SetTitleSize(fTit[1]);
2155 hPtDeltaVsTrueSig[iSig] ->
GetZaxis() -> SetTitleOffset(fOffZ[1]);
2156 hPtDeltaVsTrueSig[iSig] ->
GetZaxis() -> SetLabelFont(fTxt);
2157 hPtDeltaVsTrueSig[iSig] ->
GetZaxis() -> SetLabelSize(fLab[1]);
2158 hPtDeltaVsTrueSig[iSig] ->
GetZaxis() -> CenterTitle(fCnt);
2166 hPtDeltaVsTrackSig[iSig] ->
SetTitle(sTitle.Data());
2167 hPtDeltaVsTrackSig[iSig] -> SetTitleFont(fTxt);
2168 hPtDeltaVsTrackSig[iSig] ->
GetXaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
2170 hPtDeltaVsTrackSig[iSig] ->
GetXaxis() -> SetTitleFont(fTxt);
2171 hPtDeltaVsTrackSig[iSig] ->
GetXaxis() -> SetTitleSize(fTit[1]);
2172 hPtDeltaVsTrackSig[iSig] ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
2173 hPtDeltaVsTrackSig[iSig] ->
GetXaxis() -> SetLabelFont(fTxt);
2174 hPtDeltaVsTrackSig[iSig] ->
GetXaxis() -> SetLabelSize(fLab[1]);
2175 hPtDeltaVsTrackSig[iSig] ->
GetXaxis() -> CenterTitle(fCnt);
2176 hPtDeltaVsTrackSig[iSig] ->
GetYaxis() -> SetRangeUser(rDeltaRange[0], rDeltaRange[1]);
2178 hPtDeltaVsTrackSig[iSig] ->
GetYaxis() -> SetTitleFont(fTxt);
2179 hPtDeltaVsTrackSig[iSig] ->
GetYaxis() -> SetTitleSize(fTit[1]);
2180 hPtDeltaVsTrackSig[iSig] ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
2181 hPtDeltaVsTrackSig[iSig] ->
GetYaxis() -> SetLabelFont(fTxt);
2182 hPtDeltaVsTrackSig[iSig] ->
GetYaxis() -> SetLabelSize(fLab[1]);
2183 hPtDeltaVsTrackSig[iSig] ->
GetYaxis() -> CenterTitle(fCnt);
2185 hPtDeltaVsTrackSig[iSig] ->
GetZaxis() -> SetTitleFont(fTxt);
2186 hPtDeltaVsTrackSig[iSig] ->
GetZaxis() -> SetTitleSize(fTit[1]);
2187 hPtDeltaVsTrackSig[iSig] ->
GetZaxis() -> SetTitleOffset(fOffZ[1]);
2188 hPtDeltaVsTrackSig[iSig] ->
GetZaxis() -> SetLabelFont(fTxt);
2189 hPtDeltaVsTrackSig[iSig] ->
GetZaxis() -> SetLabelSize(fLab[1]);
2190 hPtDeltaVsTrackSig[iSig] ->
GetZaxis() -> CenterTitle(fCnt);
2193 hPtTrueVsTrackSig[iSig] ->
SetFillColor(fColSig[iSig]);
2195 hPtTrueVsTrackSig[iSig] ->
SetLineColor(fColSig[iSig]);
2198 hPtTrueVsTrackSig[iSig] ->
SetTitle(sTitle.Data());
2199 hPtTrueVsTrackSig[iSig] -> SetTitleFont(fTxt);
2200 hPtTrueVsTrackSig[iSig] ->
GetXaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
2202 hPtTrueVsTrackSig[iSig] ->
GetXaxis() -> SetTitleFont(fTxt);
2203 hPtTrueVsTrackSig[iSig] ->
GetXaxis() -> SetTitleSize(fTit[1]);
2204 hPtTrueVsTrackSig[iSig] ->
GetXaxis() -> SetTitleOffset(fOffX[1]);
2205 hPtTrueVsTrackSig[iSig] ->
GetXaxis() -> SetLabelFont(fTxt);
2206 hPtTrueVsTrackSig[iSig] ->
GetXaxis() -> SetLabelSize(fLab[1]);
2207 hPtTrueVsTrackSig[iSig] ->
GetXaxis() -> CenterTitle(fCnt);
2208 hPtTrueVsTrackSig[iSig] ->
GetYaxis() -> SetRangeUser(rPtRange[0], rPtRange[1]);
2210 hPtTrueVsTrackSig[iSig] ->
GetYaxis() -> SetTitleFont(fTxt);
2211 hPtTrueVsTrackSig[iSig] ->
GetYaxis() -> SetTitleSize(fTit[1]);
2212 hPtTrueVsTrackSig[iSig] ->
GetYaxis() -> SetTitleOffset(fOffY[1]);
2213 hPtTrueVsTrackSig[iSig] ->
GetYaxis() -> SetLabelFont(fTxt);
2214 hPtTrueVsTrackSig[iSig] ->
GetYaxis() -> SetLabelSize(fLab[1]);
2215 hPtTrueVsTrackSig[iSig] ->
GetYaxis() -> CenterTitle(fCnt);
2217 hPtTrueVsTrackSig[iSig] ->
GetZaxis() -> SetTitleFont(fTxt);
2218 hPtTrueVsTrackSig[iSig] ->
GetZaxis() -> SetTitleSize(fTit[1]);
2219 hPtTrueVsTrackSig[iSig] ->
GetZaxis() -> SetTitleOffset(fOffZ[1]);
2220 hPtTrueVsTrackSig[iSig] ->
GetZaxis() -> SetLabelFont(fTxt);
2221 hPtTrueVsTrackSig[iSig] ->
GetZaxis() -> SetLabelSize(fLab[1]);
2222 hPtTrueVsTrackSig[iSig] ->
GetZaxis() -> CenterTitle(fCnt);
2224 cout <<
" Set styles." << endl;
2227 const UInt_t fColLe = 0;
2228 const UInt_t fFilLe = 0;
2229 const UInt_t fLinLe = 0;
2230 const Float_t yObjLe = 0.1 + ((NDPtCuts + 2) * 0.05);
2231 const Float_t yObjMu = 0.1 + ((NSigCuts + 1) * 0.05);
2232 const Float_t yObjDel = 0.1 + (NSigCuts * 0.05);
2233 const Float_t yObjPro = 0.1 + (NProj * 0.05);
2234 const Float_t yObjSig = 0.1 + ((NSigCuts + 2) * 0.05);
2235 const Float_t fLegXY[
NVtx] = {0.1, 0.1, 0.3, yObjLe};
2236 const Float_t fLegMuXY[
NVtx] = {0.1, 0.1, 0.3, yObjMu};
2237 const Float_t fLegDelXY[
NVtx] = {0.1, 0.1, 0.3, yObjDel};
2238 const Float_t fLegProXY[
NVtx] = {0.1, 0.1, 0.3, yObjPro};
2239 const Float_t fLegSigXY[
NVtx] = {0.1, 0.1, 0.3, yObjSig};
2241 TLegend *
leg =
new TLegend(fLegXY[0], fLegXY[1], fLegXY[2], fLegXY[3]);
2247 leg -> SetTextAlign(fAln);
2248 leg ->
AddEntry(hPtTruth, sLegTrue.Data(),
"pf");
2249 leg ->
AddEntry(hPtTrkTru, sLegTrack.Data(),
"pf");
2250 for (Ssiz_t iCut = 0; iCut <
NDPtCuts; iCut++) {
2251 leg ->
AddEntry(hPtTrkTruCut[iCut], sLegCut[iCut].
Data(),
"pf");
2254 TLegend *legMu =
new TLegend(fLegMuXY[0], fLegMuXY[1], fLegMuXY[2], fLegMuXY[3]);
2260 legMu -> SetTextAlign(fAln);
2261 legMu ->
AddEntry(grMuProj, sLegMu.Data(),
"p");
2262 for (Ssiz_t iSig = 0; iSig <
NSigCuts; iSig++) {
2263 legMu ->
AddEntry(grMuHiProj[iSig], sLegProjSig[iSig].
Data(),
"p");
2266 TLegend *legDel =
new TLegend(fLegDelXY[0], fLegDelXY[1], fLegDelXY[2], fLegDelXY[3]);
2272 legDel -> SetTextAlign(fAln);
2273 for (Ssiz_t iSig = 0; iSig <
NSigCuts; iSig++) {
2274 legDel ->
AddEntry(fMuHiProj[iSig], sLegDelta[iSig].
Data(),
"l");
2277 TLegend *legPro =
new TLegend(fLegProXY[0], fLegProXY[1], fLegProXY[2], fLegProXY[3]);
2283 legPro -> SetTextAlign(fAln);
2284 for (Ssiz_t iProj = 0; iProj <
NProj; iProj++) {
2285 legPro ->
AddEntry(hPtDeltaProj[iProj], sLegProj[iProj].
Data(),
"pf");
2288 TLegend *legSig =
new TLegend(fLegSigXY[0], fLegSigXY[1], fLegSigXY[2], fLegSigXY[3]);
2294 legSig -> SetTextAlign(fAln);
2295 legSig ->
AddEntry(hPtTruth, sLegTrue.Data(),
"pf");
2296 legSig ->
AddEntry(hPtTrkTru, sLegTrack.Data(),
"pf");
2297 for (Ssiz_t iSig = 0; iSig <
NSigCuts; iSig++) {
2298 legSig ->
AddEntry(hPtTrkTruSig[iSig], sLegSig[iSig].
Data(),
"pf");
2300 cout <<
" Made legends." << endl;
2303 const UInt_t fColInf = 0;
2304 const UInt_t fFilInf = 0;
2305 const UInt_t fLinInf = 0;
2306 const Float_t yObjInf = 0.1 + (
NTxt * 0.05);
2307 const Float_t yObjCut = 0.1 + (
NTrkCuts * 0.05);
2308 const Float_t fInfXY[
NVtx] = {0.3, 0.1, 0.5, yObjInf};
2309 const Float_t fCutXY[
NVtx] = {0.5, 0.1, 0.7, yObjCut};
2311 TPaveText *info =
new TPaveText(fInfXY[0], fInfXY[1], fInfXY[2], fInfXY[3],
"NDC NB");
2317 info -> SetTextAlign(fAln);
2318 for (Ssiz_t iTxt = 0; iTxt <
NTxt; iTxt++) {
2319 info -> AddText(sInfo[iTxt].
Data());
2322 TPaveText *
cuts =
new TPaveText(fCutXY[0], fCutXY[1], fCutXY[2], fCutXY[3],
"NDC NB");
2328 cuts -> SetTextAlign(fAln);
2329 for (Ssiz_t iTrkCut = 0; iTrkCut <
NTrkCuts; iTrkCut++) {
2330 cuts -> AddText(sTrkCuts[iTrkCut].
Data());
2332 cout <<
" Made text." << endl;
2335 const UInt_t fColLi = 1;
2336 const UInt_t fLinLi = 9;
2337 const UInt_t fWidLi = 1;
2338 const Float_t fLinXY[
NVtx] = {rPtRange[0], 1., rPtRange[1], 1.};
2340 TLine *
line =
new TLine(fLinXY[0], fLinXY[1], fLinXY[2], fLinXY[3]);
2344 cout <<
" Made line." << endl;
2347 const UInt_t
width(750);
2348 const UInt_t width2D(1500);
2349 const UInt_t height(950);
2350 const UInt_t heightNR(750);
2351 const UInt_t fMode(0);
2352 const UInt_t fBord(2);
2353 const UInt_t fGrid(0);
2354 const UInt_t fTick(1);
2355 const UInt_t fLogX(0);
2356 const UInt_t fLogY1(1);
2357 const UInt_t fLogY2(1);
2358 const UInt_t fLogYNR(0);
2359 const UInt_t fFrame(0);
2360 const Float_t fMarginL(0.15);
2361 const Float_t fMarginR(0.02);
2362 const Float_t fMarginR2D(0.15);
2363 const Float_t fMarginT1(0.005);
2364 const Float_t fMarginT2(0.02);
2365 const Float_t fMarginTNR(0.02);
2366 const Float_t fMarginB1(0.25);
2367 const Float_t fMarginB2(0.005);
2368 const Float_t fMarginBNR(0.15);
2369 const Float_t fEffXY[
NVtx] = {0., 0., 1., 0.35};
2370 const Float_t fTrksXY[
NVtx] = {0., 0.35, 1., 1.};
2371 const Float_t fTwoDimXY[
NVtx] = {0., 0., 0.5, 1.};
2372 const Float_t fProjectXY[
NVtx] = {0.5, 0., 1., 1.};
2373 const Float_t fBeforeDPtXY[
NVtx] = {0., 0., 0.5, 1.};
2374 const Float_t fAfterDPtXY[
NVtx] = {0.5, 0., 1., 1.};
2376 TCanvas *cEffCut =
new TCanvas(
"cEfficiency_FlatCut",
"", width, height);
2377 TPad *pEffCut =
new TPad(
"pEffCut",
"", fEffXY[0], fEffXY[1], fEffXY[2], fEffXY[3]);
2378 TPad *pTrkCut =
new TPad(
"pTrkCut",
"", fTrksXY[0], fTrksXY[1], fTrksXY[2], fTrksXY[3]);
2379 cEffCut -> SetGrid(fGrid, fGrid);
2380 cEffCut -> SetTicks(fTick, fTick);
2383 pEffCut -> SetGrid(fGrid, fGrid);
2384 pEffCut -> SetTicks(fTick, fTick);
2385 pEffCut -> SetLogx(fLogX);
2394 pTrkCut -> SetGrid(fGrid, fGrid);
2395 pTrkCut -> SetTicks(fTick, fTick);
2396 pTrkCut -> SetLogx(fLogX);
2410 for (Ssiz_t iCut = 0; iCut <
NDPtCuts; iCut++) {
2411 hEffCut[iCut] ->
Draw(
"SAME");
2416 hPtTrkTru ->
Draw(
"SAME");
2417 for (Ssiz_t iCut = 0; iCut <
NDPtCuts; iCut++) {
2418 hPtTrkTruCut[iCut] ->
Draw(
"SAME");
2427 TCanvas *cEffSig =
new TCanvas(
"cEfficiency_SigmaCut",
"", width, height);
2428 TPad *pEffSig =
new TPad(
"pEffSig",
"", fEffXY[0], fEffXY[1], fEffXY[2], fEffXY[3]);
2429 TPad *pTrkSig =
new TPad(
"pTrkSig",
"", fTrksXY[0], fTrksXY[1], fTrksXY[2], fTrksXY[3]);
2430 cEffSig -> SetGrid(fGrid, fGrid);
2431 cEffSig -> SetTicks(fTick, fTick);
2434 pEffSig -> SetGrid(fGrid, fGrid);
2435 pEffSig -> SetTicks(fTick, fTick);
2436 pEffSig -> SetLogx(fLogX);
2445 pTrkSig -> SetGrid(fGrid, fGrid);
2446 pTrkSig -> SetTicks(fTick, fTick);
2447 pTrkSig -> SetLogx(fLogX);
2461 for (Ssiz_t iSig = 0; iSig <
NSigCuts; iSig++) {
2462 hEffSig[iSig] ->
Draw(
"SAME");
2467 hPtTrkTru ->
Draw(
"SAME");
2468 for (Ssiz_t iSig = 0; iSig <
NSigCuts; iSig++) {
2469 hPtTrkTruSig[iSig] ->
Draw(
"SAME");
2478 TCanvas *cRejCut =
new TCanvas(
"cReject_FlatCut",
"", width, heightNR);
2479 cRejCut -> SetGrid(fGrid, fGrid);
2480 cRejCut -> SetTicks(fTick, fTick);
2488 cRejCut -> SetLogx(fLogX);
2491 grRejCut ->
Draw(
"ALP");
2498 TCanvas *cRejSig =
new TCanvas(
"cReject_SigmaCut",
"", width, heightNR);
2499 cRejSig -> SetGrid(fGrid, fGrid);
2500 cRejSig -> SetTicks(fTick, fTick);
2508 cRejSig -> SetLogx(fLogX);
2511 grRejSig ->
Draw(
"ALP");
2518 TCanvas *cPtTruVsTrkCut =
new TCanvas(
"cPtTruthVsReco_FlatCut",
"", width2D, heightNR);
2519 TPad *pBeforeCut =
new TPad(
"pBeforeCut",
"", fBeforeDPtXY[0], fBeforeDPtXY[1], fBeforeDPtXY[2], fBeforeDPtXY[3]);
2520 TPad *pAfterCut =
new TPad(
"pAfterCut",
"", fAfterDPtXY[0], fAfterDPtXY[1], fAfterDPtXY[2], fAfterDPtXY[3]);
2521 cPtTruVsTrkCut -> SetGrid(fGrid, fGrid);
2522 cPtTruVsTrkCut -> SetTicks(fTick, fTick);
2525 pBeforeCut -> SetGrid(fGrid, fGrid);
2526 pBeforeCut -> SetTicks(fTick, fTick);
2527 pBeforeCut -> SetLogx(fLogX);
2528 pBeforeCut ->
SetLogy(fLogYNR);
2535 pAfterCut -> SetGrid(fGrid, fGrid);
2536 pAfterCut -> SetTicks(fTick, fTick);
2537 pAfterCut -> SetLogx(fLogX);
2538 pAfterCut ->
SetLogy(fLogYNR);
2545 cPtTruVsTrkCut ->
cd();
2546 pBeforeCut ->
Draw();
2547 pAfterCut ->
Draw();
2549 hPtTrueVsTrack ->
SetTitle(sBeforeTitle.Data());
2550 hPtTrueVsTrack ->
Draw(
"colz");
2553 hPtTrueVsTrackCut[iCutToDraw] ->
SetTitle(sAfterCutTitle.Data());
2554 hPtTrueVsTrackCut[iCutToDraw] ->
Draw(
"colz");
2557 cPtTruVsTrkCut ->
Write();
2558 cPtTruVsTrkCut ->
Close();
2560 TCanvas *cPtTruVsTrkSig =
new TCanvas(
"cPtTruthVsReco_SigmaCut",
"", width2D, heightNR);
2561 TPad *pBeforeSig =
new TPad(
"pBeforeSig",
"", fBeforeDPtXY[0], fBeforeDPtXY[1], fBeforeDPtXY[2], fBeforeDPtXY[3]);
2562 TPad *pAfterSig =
new TPad(
"pAfterSig",
"", fAfterDPtXY[0], fAfterDPtXY[1], fAfterDPtXY[2], fAfterDPtXY[3]);
2563 cPtTruVsTrkSig -> SetGrid(fGrid, fGrid);
2564 cPtTruVsTrkSig -> SetTicks(fTick, fTick);
2567 pBeforeSig -> SetGrid(fGrid, fGrid);
2568 pBeforeSig -> SetTicks(fTick, fTick);
2569 pBeforeSig -> SetLogx(fLogX);
2570 pBeforeSig ->
SetLogy(fLogYNR);
2577 pAfterSig -> SetGrid(fGrid, fGrid);
2578 pAfterSig -> SetTicks(fTick, fTick);
2579 pAfterSig -> SetLogx(fLogX);
2580 pAfterSig ->
SetLogy(fLogYNR);
2587 cPtTruVsTrkSig ->
cd();
2588 pBeforeSig ->
Draw();
2589 pAfterSig ->
Draw();
2591 hPtTrueVsTrack ->
Draw(
"colz");
2594 hPtTrueVsTrackSig[iSigToDraw] ->
SetTitle(sAfterSigTitle.Data());
2595 hPtTrueVsTrackSig[iSigToDraw] ->
Draw(
"colz");
2598 cPtTruVsTrkSig ->
Write();
2599 cPtTruVsTrkSig ->
Close();
2601 TCanvas *cPtDelVsTrk =
new TCanvas(
"cPtDeltaVsTrack",
"", width2D, heightNR);
2602 TPad *pTwoDim =
new TPad(
"pTwoDim",
"", fTwoDimXY[0], fTwoDimXY[1], fTwoDimXY[2], fTwoDimXY[3]);
2603 TPad *pProject =
new TPad(
"pProjections",
"", fProjectXY[0], fProjectXY[1], fProjectXY[2], fProjectXY[3]);
2604 cPtDelVsTrk -> SetGrid(fGrid, fGrid);
2605 cPtDelVsTrk -> SetTicks(fTick, fTick);
2608 pTwoDim -> SetGrid(fGrid, fGrid);
2609 pTwoDim -> SetTicks(fTick, fTick);
2610 pTwoDim -> SetLogx(fLogX);
2619 pProject -> SetGrid(fGrid, fGrid);
2620 pProject -> SetTicks(fTick, fTick);
2621 pProject -> SetLogx(fLogX);
2630 cPtDelVsTrk ->
cd();
2634 hPtDeltaVsTrack ->
Draw(
"colz");
2635 for (Ssiz_t iSig = 0; iSig <
NSigCuts; iSig++) {
2636 fMuHiProj[iSig] ->
Draw(
"same");
2637 fMuLoProj[iSig] ->
Draw(
"same");
2642 hPtDeltaProj[0] ->
Draw();
2643 fPtDeltaProj[0] ->
Draw(
"same");
2644 for (Ssiz_t iProj = 1; iProj <
NProj; iProj++) {
2645 hPtDeltaProj[iProj] ->
Draw(
"same");
2646 fPtDeltaProj[iProj] ->
Draw(
"same");
2651 cPtDelVsTrk ->
Write();
2652 cPtDelVsTrk ->
Close();
2654 TCanvas *cDeltaPt =
new TCanvas(
"cDeltaPt",
"", width, heightNR);
2655 cDeltaPt -> SetGrid(fGrid, fGrid);
2656 cDeltaPt -> SetTicks(fTick, fTick);
2664 cDeltaPt -> SetLogx(fLogX);
2671 cDeltaPt ->
Write();
2672 cDeltaPt ->
Close();
2674 TCanvas *cMuProj =
new TCanvas(
"cMuDeltaPt",
"", width, heightNR);
2675 cMuProj -> SetGrid(fGrid, fGrid);
2676 cMuProj -> SetTicks(fTick, fTick);
2684 cMuProj -> SetLogx(fLogX);
2687 grMuProj ->
Draw(
"ALP");
2688 for (Ssiz_t iSig = 0; iSig <
NSigCuts; iSig++) {
2689 grMuHiProj[iSig] ->
Draw(
"LP");
2690 grMuLoProj[iSig] ->
Draw(
"LP");
2699 TCanvas *cSigProj =
new TCanvas(
"cSigmaDeltaPt",
"", width, heightNR);
2700 cSigProj -> SetGrid(fGrid, fGrid);
2701 cSigProj -> SetTicks(fTick, fTick);
2709 cSigProj -> SetLogx(fLogX);
2712 grSigProj ->
Draw(
"ALP");
2716 cSigProj ->
Write();
2717 cSigProj ->
Close();
2718 cout <<
" Made plots." << endl;
2721 TDirectory *dNoCut = (TDirectory*) fOutput -> mkdir(
"NoCuts");
2722 TDirectory *dFlatCut = (TDirectory*) fOutput -> mkdir(
"FlatCuts");
2723 TDirectory *dSigmaCut = (TDirectory*) fOutput -> mkdir(
"SigmaCuts");
2724 TDirectory *dProject = (TDirectory*) fOutput -> mkdir(
"Projections");
2730 hPtTruth ->
Write();
2731 hPtDelta ->
Write();
2732 hPtTrack ->
Write();
2734 hPtTrkTru ->
Write();
2735 hPtDeltaVsFrac ->
Write();
2736 hPtDeltaVsTrue ->
Write();
2737 hPtDeltaVsTrack ->
Write();
2738 hPtTrueVsTrack ->
Write();
2742 grRejCut ->
Write();
2743 for (Ssiz_t iCut = 0; iCut <
NDPtCuts; iCut++) {
2744 hEffCut[iCut] ->
Write();
2745 hPtDeltaCut[iCut] ->
Write();
2746 hPtTrackCut[iCut] ->
Write();
2747 hPtFracCut[iCut] ->
Write();
2748 hPtTrkTruCut[iCut] ->
Write();
2749 hPtDeltaVsFracCut[iCut] ->
Write();
2750 hPtDeltaVsTrueCut[iCut] ->
Write();
2751 hPtDeltaVsTrackCut[iCut] ->
Write();
2752 hPtTrueVsTrackCut[iCut] ->
Write();
2757 grRejSig ->
Write();
2758 for (Ssiz_t iSig = 0; iSig <
NSigCuts; iSig++) {
2759 hEffSig[iSig] ->
Write();
2760 hPtDeltaSig[iSig] ->
Write();
2761 hPtTrackSig[iSig] ->
Write();
2762 hPtFracSig[iSig] ->
Write();
2763 hPtTrkTruSig[iSig] ->
Write();
2764 hPtDeltaVsFracSig[iSig] ->
Write();
2765 hPtDeltaVsTrueSig[iSig] ->
Write();
2766 hPtDeltaVsTrackSig[iSig] ->
Write();
2767 hPtTrueVsTrackSig[iSig] ->
Write();
2772 grMuProj ->
Write();
2773 grSigProj ->
Write();
2774 for (Ssiz_t iProj = 0; iProj <
NProj; iProj++) {
2775 hPtDeltaProj[iProj] ->
Write();
2776 fPtDeltaProj[iProj] ->
Write();
2778 for (Ssiz_t iSig = 0; iSig <
NSigCuts; iSig++) {
2779 fMuHiProj[iSig] ->
Write();
2780 fMuLoProj[iSig] ->
Write();
2781 grMuHiProj[iSig] ->
Write();
2782 grMuLoProj[iSig] ->
Write();
2784 cout <<
" Saved histograms." << endl;
2791 cout <<
" Finished delta-pt extractor script!\n" << endl;