19 #include "TMultiGraph.h" 21 #include "TLegendEntry.h" 25 #include "TPaveText.h" 29 #include "TAttMarker.h" 41 #include "CAFAna/Core/Binning.h" 62 TLatex*
prelim =
new TLatex(.9, .95,
"NOvA Simulation");
63 prelim->SetTextColor(kGray+1);
65 prelim->SetTextSize(2/30.);
66 prelim->SetTextAlign(32);
67 TGaxis::SetMaxDigits(3);
73 TLatex*
prelim =
new TLatex(.9, .95,
"NOvA Preliminary");
74 prelim->SetTextColor(
kBlue);
76 prelim->SetTextSize(2/30.);
77 prelim->SetTextAlign(32);
78 TGaxis::SetMaxDigits(3);
87 TFile *
f1=
new TFile (
"/nova/app/users/acedeno/tag_releasesS18-06-14/NDAna/ncpi0_semi_inc_png_cvn/BDTSpectrum_FlatFluxdatasetslh15.root",
"READ");
104 const double myPOT = 8.12e20;
108 auto hndMC3PRem = sndMC3PRem.ToTH1(myPOT);
109 auto hndNCNonPi0Rem = sndNCNonPi0Rem.ToTH1(myPOT);
110 auto hndCCPi0Rem = sndCCPi0Rem.ToTH1(myPOT);
111 auto hndNCPi0Rem = sndNCPi0Rem.ToTH1(myPOT);
112 auto hndNCPi0LRem = sndNCPi0LRem.ToTH1(myPOT);
113 auto hndCCNonPi0Rem = sndCCNonPi0Rem.ToTH1(myPOT);
116 TString
outDir2 =
"/nova/app/users/acedeno/tag_releasesS18-06-14/NDAna/ncpi0_semi_inc_png_cvn/BDTGStudy";
127 {
"/BDTG_Distribution"},
128 {
"/BDTG_Distribution_stacked"},
135 {
"_07-03-18_BDTG_FlatFlux_slh_15"},
143 TCanvas *cs6 =
new TCanvas(
"cs6",
"cs6",10,10,855,500);
146 hndNCPi0LRem->SetLineColor(
kOrange);
147 hndNCPi0LRem->SetLineWidth(1);
149 hndCCPi0Rem->SetLineColor(kMagenta);
150 hndCCPi0Rem->SetLineWidth(1);
153 hndCCNonPi0Rem->SetLineColor(kAzure);
154 hndCCNonPi0Rem->SetLineWidth(1);
157 hndNCNonPi0Rem->SetLineColor(kCyan);
158 hndNCNonPi0Rem->SetLineWidth(1);
161 hndNCPi0Rem->SetLineColor(
kGreen);
162 hndNCPi0Rem->SetLineWidth(1);
165 hndMC3PRem->SetLineColor(
kRed);
166 hndMC3PRem->SetLineWidth(1);
172 hndMC3PRem->Draw(
"hist E1");
173 hndCCPi0Rem->Draw(
"hist E1 same");
174 hndCCNonPi0Rem->Draw(
"hist E1 same");
175 hndNCPi0LRem->Draw(
"hist E1 same");
176 hndNCPi0Rem->Draw(
"hist E1 same");
177 hndNCNonPi0Rem->Draw(
"hist E1 same");
178 hndMC3PRem->SetMaximum(6000);
179 hndMC3PRem->GetXaxis()->SetRangeUser(-1.,1.);
180 hndMC3PRem->GetYaxis()->SetTitle(
"Events/ 8e20 POT");
183 TLegend * leg6 =
new TLegend(0.55,0.5,0.65,0.82);
185 leg6->SetTextFont(42);
186 leg6->SetTextSize(0.05);
187 leg6->SetFillStyle(0);
189 leg6->AddEntry(hndMC3PRem,
"Selection",
"kFullDotSmall");
190 leg6->AddEntry(hndCCPi0Rem,
"CCPi0",
"kFullDotSmall");
191 leg6->AddEntry(hndCCNonPi0Rem,
"CCNonPi0",
"kFullDotSmall");
192 leg6->AddEntry(hndNCPi0LRem,
"NC w/Pi0 < 0.3 E",
"kFullDotSmall");
193 leg6->AddEntry(hndNCPi0Rem,
"NC w/Pi0 > 0.3 E",
"kFullDotSmall");
194 leg6->AddEntry(hndNCNonPi0Rem,
"NCNonPi0",
"kFullDotSmall");
199 cs6->SaveAs(outDir2+CurrentFile[0].
name.c_str()+CurrentDate[0].
date.c_str()+
".png");
205 TCanvas *cs4 =
new TCanvas(
"cs4",
"cs4",10,10,855,500);
206 THStack *hs4 =
new THStack(
"hs4",
"");
208 hndNCPi0LRem->SetLineColor(kBlack);
209 hndNCPi0LRem->SetFillColor(
kOrange);
210 hndNCPi0LRem->SetLineWidth(1);;
211 hndNCPi0LRem->SetMarkerStyle(21);
212 hndNCPi0LRem->SetMarkerColor(
kOrange);
213 hs4->Add(hndNCPi0LRem);
216 hndCCPi0Rem->SetLineColor(kBlack);
217 hndCCPi0Rem->SetFillColor(kMagenta);
218 hndCCPi0Rem->SetLineWidth(1);;
219 hndCCPi0Rem->SetMarkerStyle(21);
220 hndCCPi0Rem->SetMarkerColor(kMagenta);
222 hs4->Add(hndCCPi0Rem);
224 hndCCNonPi0Rem->SetLineColor(kBlack);
225 hndCCNonPi0Rem->SetFillColor(kAzure);
226 hndCCNonPi0Rem->SetLineWidth(1);;
227 hndCCNonPi0Rem->SetMarkerStyle(21);
228 hndCCNonPi0Rem->SetMarkerColor(kAzure);
229 hs4->Add(hndCCNonPi0Rem);
232 hndNCNonPi0Rem->SetLineColor(kBlack);
233 hndNCNonPi0Rem->SetFillColor(kCyan);
234 hndNCNonPi0Rem->SetLineWidth(1);;
235 hndNCNonPi0Rem->SetMarkerStyle(21);
236 hndNCNonPi0Rem->SetMarkerColor(kCyan);
237 hs4->Add(hndNCNonPi0Rem);
240 hndNCPi0Rem->SetLineColor(kBlack);
241 hndNCPi0Rem->SetFillColor(
kGreen);
242 hndNCPi0Rem->SetLineWidth(1);;
243 hndNCPi0Rem->SetMarkerStyle(21);
244 hndNCPi0Rem->SetMarkerColor(
kGreen);
245 hs4->Add(hndNCPi0Rem);
248 TText T4; T4.SetTextFont(42); T4.SetTextAlign(21);
249 hs4->Draw(
"hist"); T4.DrawTextNDC(.5,.95,
"");
250 hs4->GetYaxis()->SetTitle(
"Events/ 8e20 POT");
251 hs4->GetXaxis()->SetTitle(
"BDTG");
252 hs4->SetMaximum(6000);
253 hs4->GetXaxis()->SetRangeUser(-1.,1.);
257 hndMC3PRem->SetLineColor(
kRed);
258 hndMC3PRem->SetLineWidth(3);
263 hndMC3PRem->Draw(
"same");
265 TLegend * leg4 =
new TLegend(0.55,0.5,0.65,0.82);
267 leg4->SetTextFont(42);
268 leg4->SetTextSize(0.05);
269 leg4->SetFillStyle(0);
271 leg4->AddEntry(hndMC3PRem,
"Selection",
"kFullDotSmall");
272 leg4->AddEntry(hndCCPi0Rem,
"CCPi0",
"kFullDotSmall");
273 leg4->AddEntry(hndCCNonPi0Rem,
"CCNonPi0",
"kFullDotSmall");
274 leg4->AddEntry(hndNCPi0LRem,
"NC w/Pi0 < 0.3 E",
"kFullDotSmall");
275 leg4->AddEntry(hndNCPi0Rem,
"NC w/Pi0 > 0.3 E",
"kFullDotSmall");
276 leg4->AddEntry(hndNCNonPi0Rem,
"NCNonPi0",
"kFullDotSmall");
281 cs4->SaveAs(outDir2+CurrentFile[1].
name.c_str()+CurrentDate[0].
date.c_str()+
".png");
::xsd::cxx::tree::date< char, simple_type > date
Cuts and Vars for the 2020 FD DiF Study.
const DateUp CurrentDate[Date]
const FileUp CurrentFile[File]
static std::unique_ptr< Spectrum > LoadFrom(TDirectory *dir, const std::string &name)