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);
91 TFile *
f1=
new TFile (
"/nova/app/users/acedeno/tag_releasesS18-06-14/NDAna/ncpi0_semi_inc_png_cvn/new2pngvariable_prod4_03GeV_Cuts_1CVN_07_2CVN_09_WITHDATA_GenieLike.root",
"READ");
103 auto sndNumiMasssdataNoCuts =*
Spectrum::LoadFrom(f1->GetDirectory(
"ndNumiMassdataNoCuts"));
126 auto sndNumiMassdata2CVN1=*
Spectrum::LoadFrom(f1->GetDirectory(
"ndNumiMassdata2CVN"));
217 auto sndMC3PPrelCVN2Cuts =*
Spectrum::LoadFrom(f1->GetDirectory(
"ndMC3PPrelCVN2Cuts"));
218 auto sndNCNonPi0PrelCVN2Cuts =*
Spectrum::LoadFrom(f1->GetDirectory(
"ndNCNonPi0PrelCVN2Cuts"));
219 auto sndCCPi0PrelCVN2Cuts =*
Spectrum::LoadFrom(f1->GetDirectory(
"ndCCPi0PrelCVN2Cuts"));
220 auto sndNCPi0PrelCVN2Cuts =*
Spectrum::LoadFrom(f1->GetDirectory(
"ndNCPi0PrelCVN2Cuts"));
221 auto sndNCPi0LPrelCVN2Cuts =*
Spectrum::LoadFrom(f1->GetDirectory(
"ndNCPi0LPrelCVN2Cuts"));
222 auto sndCCNonPi0PrelCVN2Cuts =*
Spectrum::LoadFrom(f1->GetDirectory(
"ndCCNonPi0PrelCVN2Cuts"));
243 const double myPOT = 8.12e20;
252 auto hndNoCuts = sndNoCuts.ToTH1(myPOT);
253 auto hndNoCutCCNonPi0 = sndNoCutCCNonPi0.ToTH1(myPOT);
254 auto hndNoCutNCNonPi0 = sndNoCutNCNonPi0.ToTH1(myPOT);
255 auto hndNoCutCCPi0 = sndNoCutCCPi0.ToTH1(myPOT);
256 auto hndNoCutNCPi0 = sndNoCutNCPi0.ToTH1(myPOT);
257 auto hndNoCutNCPi0L = sndNoCutNCPi0L.ToTH1(myPOT);
258 auto hndNumiMasssdataNoCuts = sndNumiMasssdataNoCuts.ToTH1(myPOT);
261 auto hndMC = sndMC.ToTH1(myPOT);
262 auto hndMC3P = sndMC3P.ToTH1(myPOT);
263 auto hndNCNonPi0 = sndNCNonPi0.ToTH1(myPOT);
264 auto hndCCPi0 = sndCCPi0.ToTH1(myPOT);
265 auto hndNCPi0 = sndNCPi0.ToTH1(myPOT);
266 auto hndNCPi0L = sndNCPi0L.ToTH1(myPOT);
267 auto hndCCNonPi0 = sndCCNonPi0.ToTH1(myPOT);
268 auto hndNumiMassdata2CVN = sndNumiMassdata2CVN.ToTH1(myPOT);
271 auto hndMC1 = sndMC.ToTH1(myPOT);
272 auto hndMC3P1 = sndMC3P.ToTH1(myPOT);
273 auto hndNCNonPi01 = sndNCNonPi0.ToTH1(myPOT);
274 auto hndCCPi01 = sndCCPi0.ToTH1(myPOT);
275 auto hndNCPi01 = sndNCPi0.ToTH1(myPOT);
276 auto hndNCPi0L1 = sndNCPi0L.ToTH1(myPOT);
277 auto hndCCNonPi01 = sndCCNonPi0.ToTH1(myPOT);
278 auto hndNumiMassdata2CVN1 = sndNumiMassdata2CVN1.ToTH1(myPOT);
281 auto hndMC1Cut = sndMC1Cut.ToTH1(myPOT);
282 auto hndMC3P1Cut = sndMC3P1Cut.ToTH1(myPOT);
283 auto hndNCNonPi01Cut = sndNCNonPi01Cut.ToTH1(myPOT);
284 auto hndCCPi01Cut = sndCCPi01Cut.ToTH1(myPOT);
285 auto hndNCPi01Cut = sndNCPi01Cut.ToTH1(myPOT);
286 auto hndNCPi0L1Cut = sndNCPi0L1Cut.ToTH1(myPOT);
287 auto hndCCNonPi01Cut = sndCCNonPi01Cut.ToTH1(myPOT);
288 auto hndNumiMassdataCVN = sndNumiMassdataCVN.ToTH1(myPOT);
291 auto hndMC1Cut2 = sndMC1Cut2.ToTH1(myPOT);
292 auto hndMC3P1Cut2 = sndMC3P1Cut2.ToTH1(myPOT);
293 auto hndNCNonPi01Cut2 = sndNCNonPi01Cut2.ToTH1(myPOT);
294 auto hndCCPi01Cut2 = sndCCPi01Cut2.ToTH1(myPOT);
295 auto hndNCPi01Cut2 = sndNCPi01Cut2.ToTH1(myPOT);
296 auto hndNCPi0L1Cut2 = sndNCPi0L1Cut2.ToTH1(myPOT);
297 auto hndCCNonPi01Cut2 = sndCCNonPi01Cut2.ToTH1(myPOT);
298 auto hndNumiMassdataCVN2 = sndNumiMassdataCVN2.ToTH1(myPOT);
302 auto hndCVNMC3P = sndCVNMC3P.ToTH1(myPOT);
303 auto hndCVNNCNonPi0 = sndCVNNCNonPi0.ToTH1(myPOT);
304 auto hndCVNCCPi0 = sndCVNCCPi0.ToTH1(myPOT);
305 auto hndCVNNCPi0 = sndCVNNCPi0.ToTH1(myPOT);
306 auto hndCVNNCPi0L = sndCVNNCPi0L.ToTH1(myPOT);
307 auto hndCVNCCNonPi0 = sndCVNCCNonPi0.ToTH1(myPOT);
308 auto hndNumiCVN2PData = sndNumiCVN2PData.ToTH1(myPOT);
311 auto hndCVNMC3PR = sndCVNMC3PR.ToTH1(myPOT);
312 auto hndCVNNCNonPi0R = sndCVNNCNonPi0R.ToTH1(myPOT);
313 auto hndCVNCCPi0R = sndCVNCCPi0R.ToTH1(myPOT);
314 auto hndCVNNCPi0R = sndCVNNCPi0R.ToTH1(myPOT);
315 auto hndCVNNCPi0LR = sndCVNNCPi0LR.ToTH1(myPOT);
316 auto hndCVNCCNonPi0R = sndCVNCCNonPi0R.ToTH1(myPOT);
317 auto hndNumiCVN2PRData = sndNumiCVN2PRData.ToTH1(myPOT);
320 auto hnd2CVNMC3PR = snd2CVNMC3PR.ToTH1(myPOT);
321 auto hnd2CVNNCNonPi0R = snd2CVNNCNonPi0R.ToTH1(myPOT);
322 auto hnd2CVNCCPi0R = snd2CVNCCPi0R.ToTH1(myPOT);
323 auto hnd2CVNNCPi0R = snd2CVNNCPi0R.ToTH1(myPOT);
324 auto hnd2CVNNCPi0LR = snd2CVNNCPi0LR.ToTH1(myPOT);
325 auto hnd2CVNCCNonPi0R = snd2CVNCCNonPi0R.ToTH1(myPOT);
326 auto hndNumi2CVN2PRData = sndNumi2CVN2PRData.ToTH1(myPOT);
330 auto hndMC3PRem = sndMC3PRem.ToTH1(myPOT);
331 auto hndNCNonPi0Rem = sndNCNonPi0Rem.ToTH1(myPOT);
332 auto hndCCPi0Rem = sndCCPi0Rem.ToTH1(myPOT);
333 auto hndNCPi0Rem = sndNCPi0Rem.ToTH1(myPOT);
334 auto hndNCPi0LRem = sndNCPi0LRem.ToTH1(myPOT);
335 auto hndCCNonPi0Rem = sndCCNonPi0Rem.ToTH1(myPOT);
336 auto hndNumiRemIDData = sndNumiRemIDData.ToTH1(myPOT);
340 auto hndMC3PPrel = sndMC3PPrel.ToTH1(myPOT);
341 auto hndNCNonPi0Prel = sndNCNonPi0Prel.ToTH1(myPOT);
342 auto hndCCPi0Prel = sndCCPi0Prel.ToTH1(myPOT);
343 auto hndNCPi0Prel = sndNCPi0Prel.ToTH1(myPOT);
344 auto hndNCPi0LPrel = sndNCPi0LPrel.ToTH1(myPOT);
345 auto hndCCNonPi0Prel = sndCCNonPi0Prel.ToTH1(myPOT);
349 auto hndMC3PPrelCVN = sndMC3PPrelCVN.ToTH1(myPOT);
350 auto hndNCNonPi0PrelCVN = sndNCNonPi0PrelCVN.ToTH1(myPOT);
351 auto hndCCPi0PrelCVN = sndCCPi0PrelCVN.ToTH1(myPOT);
352 auto hndNCPi0PrelCVN = sndNCPi0PrelCVN.ToTH1(myPOT);
353 auto hndNCPi0LPrelCVN = sndNCPi0LPrelCVN.ToTH1(myPOT);
354 auto hndCCNonPi0PrelCVN = sndCCNonPi0PrelCVN.ToTH1(myPOT);
358 auto hndMC3PPrelCVN2Cuts = sndMC3PPrelCVN2Cuts.ToTH1(myPOT);
359 auto hndNCNonPi0PrelCVN2Cuts = sndNCNonPi0PrelCVN2Cuts.ToTH1(myPOT);
360 auto hndCCPi0PrelCVN2Cuts = sndCCPi0PrelCVN2Cuts.ToTH1(myPOT);
361 auto hndNCPi0PrelCVN2Cuts = sndNCPi0PrelCVN2Cuts.ToTH1(myPOT);
362 auto hndNCPi0LPrelCVN2Cuts = sndNCPi0LPrelCVN2Cuts.ToTH1(myPOT);
363 auto hndCCNonPi0PrelCVN2Cuts = sndCCNonPi0PrelCVN2Cuts.ToTH1(myPOT);
367 auto hndMC3PPrelR = sndMC3PPrelR.ToTH1(myPOT);
368 auto hndNCNonPi0PrelR = sndNCNonPi0PrelR.ToTH1(myPOT);
369 auto hndCCPi0PrelR = sndCCPi0PrelR.ToTH1(myPOT);
370 auto hndNCPi0PrelR = sndNCPi0PrelR.ToTH1(myPOT);
371 auto hndNCPi0LPrelR = sndNCPi0LPrelR.ToTH1(myPOT);
372 auto hndCCNonPi0PrelR = sndCCNonPi0PrelR.ToTH1(myPOT);
376 auto hnd2dplot2prong = snd2dplot2prong.ToTH2(myPOT);
378 TString
outDir2 =
"/nova/app/users/acedeno/tag_releasesS18-06-14/NDAna/ncpi0_semi_inc_png_cvn/NewPNGCVAnalysis";
389 {
"/plot1_2Pr_03GeV_Cuts_1CVN_07_2CVN_09_NOCUTS_Pi0Mass"},
390 {
"/plot2_2Pr_03GeV_Cuts_1CVN_07_2CVN_09_Selectionection"},
391 {
"/plot3_2Pr_03GeV_Cuts_1CVN_07_2CVN_09_StackedInvariantMass2Cuts"},
392 {
"/plot3b_2Pr_03GeV_Cuts_1CVN_07_StackedInvariantMass1Cut"},
393 {
"/plot4_2Pr_03GeV_Cuts_{RemID+Presel+1CVN_07}-CVN-Analysis-{Var-2P-PhotonID}"},
394 {
"/plot5_2Pr_03GeV_Cuts_{RemID+Selection}-CVN-Analysis-{Var-1P-PhotonID}"},
395 {
"/plot6_2Pr_03GeV_Cuts_{RemID+Selection}-CVN-Analysis-{Var-2P-PhotonID}"},
396 {
"/plot7_2Pr_03GeV_Cuts_{Selection}-CVN-Analysis-{Var-RemID}"},
397 {
"/plot8_2Pr_03GeV_Cuts_{Selection}-TrueEnergy-Analysis-{Var-TruePi0En}"},
398 {
"/plot9_2Pr_03GeV_Cuts_{RemID+Selection}-TrueEnergy-Analysis-{Var-TruePi0En}"},
399 {
"/plot10_2Pr_03GeV_Cuts_{RemID+Selection+1CVN_07}-TrueEnergy-Analysis-{Var-TruePi0En}"},
400 {
"/plot10b_2Pr_03GeV_Cuts_{RemID+Selection+1CVN_07+2CVN_09}-TrueEnergy-Analysis-{Var-TruePi0En}"},
401 {
"/plot11_2Pr_2dplot_03GeV_Cuts_1CVN_07_2CVN_09"},
402 {
"/plot13_2Pr_03GeV_Cuts_1CVN_07_2CVN_09_UnstackedInvariantMass2Cuts"},
403 {
"/plot13b_2Pr_03GeV_Cuts_1CVN_07_UnstackedInvariantMass1Cut"},
404 {
"/plot11b_2Pr_2dplot_LogVersion_03GeV_Cuts_1CVN_07_2CVN_09"},
410 {
"_07--18_newPng_GenieLike"},
419 TCanvas *c1_ndNoCuts =
new TCanvas(
"c1_ndNoCuts",
"ndNoCuts",800,500);
422 hndNoCuts->SetLineColor(
kRed+2);
427 hndNoCutCCPi0->SetLineColor(kMagenta+2);
432 hndNoCutCCNonPi0->SetLineColor(kAzure+2);
437 hndNoCutNCPi0L->SetLineColor(
kOrange+2);
442 hndNoCutNCPi0->SetLineColor(
kGreen+2);
447 hndNoCutNCNonPi0->SetLineColor(kCyan-6);
453 hndNumiMasssdataNoCuts->SetMarkerStyle(2);
454 hndNumiMasssdataNoCuts->SetMarkerSize(3);
456 hndNoCuts->Draw(
"hist E1");
457 hndNoCutCCPi0->Draw(
"hist E1 same");
458 hndNoCutCCNonPi0->Draw(
"hist E1 same");
459 hndNoCutNCPi0L->Draw(
"hist E1 same");
460 hndNoCutNCPi0->Draw(
"hist E1 same");
461 hndNoCutNCNonPi0->Draw(
"hist E1 same");
464 hndNoCuts->GetYaxis()->SetTitle(
"Events/8e20 POT");
465 hndNoCuts->SetTitleOffset(0.95,
"Y");
468 TLegend * leg0 =
new TLegend(0.45,0.26,0.55,0.5);
469 leg0->SetTextFont(42);
470 leg0->SetTextSize(0.05);
472 leg0->AddEntry(hndNoCuts,
"No Cuts",
"l");
473 leg0->AddEntry(hndNoCutCCPi0,
"CCPi0",
"l");
474 leg0->AddEntry(hndNoCutCCNonPi0,
"CCNonPi0",
"l");
475 leg0->AddEntry(hndNoCutNCPi0L,
"NC w/Pi0 < 0.3 E",
"l");
476 leg0->AddEntry(hndNoCutNCPi0,
"NC w/Pi0 > 0.3 E",
"l");
477 leg0->AddEntry(hndNoCutNCNonPi0,
"NCNonPi0",
"l");
482 c1_ndNoCuts->SaveAs(outDir2+CurrentFile[0].
name.c_str()+CurrentDate[0].
date.c_str()+
".png");
489 TCanvas *c1_ndMC =
new TCanvas(
"c1_ndMC",
"ndMC",800,500);
492 hndMC->SetLineColor(
kBlue);
496 hndMC3P->SetLineColor(
kRed);
501 hndMC->Draw(
"hist E1");
502 hndMC3P->Draw(
"hist E1 same");
503 hndMC->GetYaxis()->SetTitle(
"Events/8e20 POT");
504 hndMC->SetTitleOffset(0.95,
"Y");
507 TLegend *
leg =
new TLegend(0.55,0.7,0.65,0.8);
508 leg->SetTextFont(42);
509 leg->SetTextSize(0.05);
510 leg->AddEntry(hndMC,
"Fid.+ Cont. (PreSel)",
"l");
511 leg->AddEntry(hndMC3P,
"Selection",
"l");
515 c1_ndMC->SaveAs(outDir2+CurrentFile[1].
name.c_str()+CurrentDate[0].
date.c_str()+
".png");
522 TCanvas *cs2 =
new TCanvas(
"cs2",
"cs2",10,10,855,500);
523 THStack *hs2 =
new THStack(
"hs2",
" ");
535 hndNCPi0L->SetLineColor(kBlack);
536 hndNCPi0L->SetFillColor(
kOrange);
537 hndNCPi0L->SetLineWidth(1);
538 hndNCPi0L->SetMarkerStyle(21);
543 hndCCPi0->SetLineColor(kBlack);
544 hndCCPi0->SetFillColor(kMagenta);
545 hndCCPi0->SetLineWidth(1);
546 hndCCPi0->SetMarkerStyle(21);
547 hndCCPi0->SetMarkerColor(kMagenta);
551 hndCCNonPi0->SetLineColor(kBlack);
552 hndCCNonPi0->SetFillColor(kAzure);
553 hndCCNonPi0->SetLineWidth(1);
554 hndCCNonPi0->SetMarkerStyle(21);
555 hndCCNonPi0->SetMarkerColor(kAzure);
556 hs2->Add(hndCCNonPi0);
559 hndNCNonPi0->SetLineColor(kBlack);
560 hndNCNonPi0->SetFillColor(kCyan);
561 hndNCNonPi0->SetLineWidth(1);
562 hndNCNonPi0->SetMarkerStyle(21);
563 hndNCNonPi0->SetMarkerColor(kCyan);
564 hs2->Add(hndNCNonPi0);
567 hndNCPi0->SetLineColor(kBlack);
568 hndNCPi0->SetFillColor(
kGreen);
569 hndNCPi0->SetLineWidth(1);
570 hndNCPi0->SetMarkerStyle(21);
571 hndNCPi0->SetMarkerColor(
kGreen);
575 TText
T; T.SetTextFont(42); T.SetTextAlign(21);
576 hs2->Draw(
"hist"); T.DrawTextNDC(.5,.95,
"");
577 hs2->GetYaxis()->SetTitle(
"Events/ 8e20 POT");
578 hs2->GetXaxis()->SetTitle(
"#pi^{0} Reconstructed Mass [GeV]");
580 hs2->GetXaxis()->SetRangeUser(0.,0.4);
584 hndMC3P->SetLineColor(
kRed);
585 hndMC3P->SetLineWidth(3);
587 hndNumiMassdata2CVN->SetMarkerSize(3);
589 hndMC3P->Draw(
"same");
591 TLegend * leg2 =
new TLegend(0.62,0.5,0.73,0.82);
593 leg2->SetTextFont(42);
594 leg2->SetTextSize(0.05);
595 leg2->SetFillStyle(0);
597 leg2->AddEntry(hndMC3P,
"Selection",
"kFullDotSmall");
598 leg2->AddEntry(hndCCPi0,
"CCPi0",
"kFullDotSmall");
599 leg2->AddEntry(hndCCNonPi0,
"CCNonPi0",
"kFullDotSmall");
600 leg2->AddEntry(hndNCPi0L,
"NC w/Pi0 < 0.3 E",
"kFullDotSmall");
601 leg2->AddEntry(hndNCPi0,
"NC w/Pi0 > 0.3 E",
"kFullDotSmall");
602 leg2->AddEntry(hndNCNonPi0,
"NCNonPi0",
"kFullDotSmall");
608 cs2->SaveAs(outDir2+CurrentFile[2].
name.c_str()+CurrentDate[0].
date.c_str()+
".png");
614 TCanvas *cs2b =
new TCanvas(
"cs2b",
"cs2b",10,10,855,500);
615 THStack *hs2b =
new THStack(
"hs2b",
" ");
627 hndNCPi0L1Cut->SetLineColor(kBlack);
628 hndNCPi0L1Cut->SetFillColor(
kOrange);
629 hndNCPi0L1Cut->SetLineWidth(1);
630 hndNCPi0L1Cut->SetMarkerStyle(21);
631 hndNCPi0L1Cut->SetMarkerColor(
kOrange);
632 hs2b->Add(hndNCPi0L1Cut);
635 hndCCPi01Cut->SetLineColor(kBlack);
636 hndCCPi01Cut->SetFillColor(kMagenta);
637 hndCCPi01Cut->SetLineWidth(1);
638 hndCCPi01Cut->SetMarkerStyle(21);
639 hndCCPi01Cut->SetMarkerColor(kMagenta);
641 hs2b->Add(hndCCPi01Cut);
643 hndCCNonPi01Cut->SetLineColor(kBlack);
644 hndCCNonPi01Cut->SetFillColor(kAzure);
645 hndCCNonPi01Cut->SetLineWidth(1);
646 hndCCNonPi01Cut->SetMarkerStyle(21);
647 hndCCNonPi01Cut->SetMarkerColor(kAzure);
648 hs2b->Add(hndCCNonPi01Cut);
651 hndNCNonPi01Cut->SetLineColor(kBlack);
652 hndNCNonPi01Cut->SetFillColor(kCyan);
653 hndNCNonPi01Cut->SetLineWidth(1);
654 hndNCNonPi01Cut->SetMarkerStyle(21);
655 hndNCNonPi01Cut->SetMarkerColor(kCyan);
656 hs2b->Add(hndNCNonPi01Cut);
659 hndNCPi01Cut->SetLineColor(kBlack);
660 hndNCPi01Cut->SetFillColor(
kGreen);
661 hndNCPi01Cut->SetLineWidth(1);
662 hndNCPi01Cut->SetMarkerStyle(21);
663 hndNCPi01Cut->SetMarkerColor(
kGreen);
664 hs2b->Add(hndNCPi01Cut);
667 TText T99; T99.SetTextFont(42); T99.SetTextAlign(21);
668 hs2b->Draw(
"hist"); T99.DrawTextNDC(.5,.95,
"");
669 hs2b->GetYaxis()->SetTitle(
"Events/ 8e20 POT");
670 hs2b->GetXaxis()->SetTitle(
"#pi^{0} Reconstructed Mass [GeV]");
671 hs2b->SetMaximum(2500);
672 hs2b->GetXaxis()->SetRangeUser(0.,0.4);
676 hndMC3P1Cut->SetLineColor(
kRed);
677 hndMC3P1Cut->SetLineWidth(3);
679 hndNumiMassdataCVN->SetMarkerSize(3);
681 hndMC3P1Cut->Draw(
"same");
683 TLegend * leg2b =
new TLegend(0.62,0.5,0.73,0.82);
685 leg2b->SetTextFont(42);
686 leg2b->SetTextSize(0.05);
687 leg2b->SetFillStyle(0);
689 leg2b->AddEntry(hndMC3P1Cut,
"Selection",
"kFullDotSmall");
690 leg2b->AddEntry(hndCCPi01Cut,
"CCPi0",
"kFullDotSmall");
691 leg2b->AddEntry(hndCCNonPi01Cut,
"CCNonPi0",
"kFullDotSmall");
692 leg2b->AddEntry(hndNCPi0L1Cut,
"NC w/Pi0 < 0.3 E",
"kFullDotSmall");
693 leg2b->AddEntry(hndNCPi01Cut,
"NC w/Pi0 > 0.3 E",
"kFullDotSmall");
694 leg2b->AddEntry(hndNCNonPi01Cut,
"NCNonPi0",
"kFullDotSmall");
699 cs2b->SaveAs(outDir2+CurrentFile[3].
name.c_str()+CurrentDate[0].
date.c_str()+
".png");
704 TCanvas *cs3 =
new TCanvas(
"cs3",
"cs3",10,10,855,500);
705 THStack *hs3 =
new THStack(
"hs3",
"");
717 hndCVNNCPi0L->SetLineColor(kBlack);
718 hndCVNNCPi0L->SetFillColor(
kOrange);
719 hndCVNNCPi0L->SetLineWidth(1);;
720 hndCVNNCPi0L->SetMarkerStyle(21);
721 hndCVNNCPi0L->SetMarkerColor(
kOrange);
722 hs3->Add(hndCVNNCPi0L);
725 hndCVNCCPi0->SetLineColor(kBlack);
726 hndCVNCCPi0->SetFillColor(kMagenta);
727 hndCVNCCPi0->SetLineWidth(1);;
728 hndCVNCCPi0->SetMarkerStyle(21);
729 hndCVNCCPi0->SetMarkerColor(kMagenta);
731 hs3->Add(hndCVNCCPi0);
733 hndCVNCCNonPi0->SetLineColor(kBlack);
734 hndCVNCCNonPi0->SetFillColor(kAzure);
735 hndCVNCCNonPi0->SetLineWidth(1);;
736 hndCVNCCNonPi0->SetMarkerStyle(21);
737 hndCVNCCNonPi0->SetMarkerColor(kAzure);
738 hs3->Add(hndCVNCCNonPi0);
741 hndCVNNCNonPi0->SetLineColor(kBlack);
742 hndCVNNCNonPi0->SetFillColor(kCyan);
743 hndCVNNCNonPi0->SetLineWidth(1);;
744 hndCVNNCNonPi0->SetMarkerStyle(21);
745 hndCVNNCNonPi0->SetMarkerColor(kCyan);
746 hs3->Add(hndCVNNCNonPi0);
749 hndCVNNCPi0->SetLineColor(kBlack);
750 hndCVNNCPi0->SetFillColor(
kGreen);
751 hndCVNNCPi0->SetLineWidth(1);;
752 hndCVNNCPi0->SetMarkerStyle(21);
753 hndCVNNCPi0->SetMarkerColor(
kGreen);
754 hs3->Add(hndCVNNCPi0);
757 TText T2; T2.SetTextFont(42); T2.SetTextAlign(21);
758 hs3->Draw(
"hist"); T2.DrawTextNDC(.5,.95,
"");
759 hs3->GetYaxis()->SetTitle(
"Events/ 8e20 POT");
760 hs3->GetXaxis()->SetTitle(
"Prong 2 CVN Gamma ID");
762 hs3->GetXaxis()->SetRangeUser(0,1);
766 hndCVNMC3P->SetLineColor(
kRed);
767 hndCVNMC3P->SetLineWidth(3);
769 hndNumiCVN2PData->SetMarkerSize(3);
771 hndCVNMC3P->Draw(
"same");
773 TLegend * leg3 =
new TLegend(0.55,0.5,0.65,0.82);
775 leg3->SetTextFont(42);
776 leg3->SetTextSize(0.05);
777 leg3->SetFillStyle(0);
779 leg3->AddEntry(hndCVNMC3P,
"Selection",
"kFullDotSmall");
780 leg3->AddEntry(hndCVNCCPi0,
"CCPi0",
"kFullDotSmall");
781 leg3->AddEntry(hndCVNCCNonPi0,
"CCNonPi0",
"kFullDotSmall");
782 leg3->AddEntry(hndCVNNCPi0L,
"NC w/Pi0 < 0.3 E",
"kFullDotSmall");
783 leg3->AddEntry(hndCVNNCPi0,
"NC w/Pi0 > 0.3 E",
"kFullDotSmall");
784 leg3->AddEntry(hndCVNNCNonPi0,
"NCNonPi0",
"kFullDotSmall");
789 cs3->SaveAs(outDir2+CurrentFile[4].
name.c_str()+CurrentDate[0].
date.c_str()+
".png");
795 TCanvas *cs4 =
new TCanvas(
"cs4",
"cs4",10,10,855,500);
796 THStack *hs4 =
new THStack(
"hs4",
"");
808 hndCVNNCPi0LR->SetLineColor(kBlack);
809 hndCVNNCPi0LR->SetFillColor(
kOrange);
810 hndCVNNCPi0LR->SetLineWidth(1);;
811 hndCVNNCPi0LR->SetMarkerStyle(21);
812 hndCVNNCPi0LR->SetMarkerColor(
kOrange);
813 hs4->Add(hndCVNNCPi0LR);
816 hndCVNCCPi0R->SetLineColor(kBlack);
817 hndCVNCCPi0R->SetFillColor(kMagenta);
818 hndCVNCCPi0R->SetLineWidth(1);;
819 hndCVNCCPi0R->SetMarkerStyle(21);
820 hndCVNCCPi0R->SetMarkerColor(kMagenta);
822 hs4->Add(hndCVNCCPi0R);
824 hndCVNCCNonPi0R->SetLineColor(kBlack);
825 hndCVNCCNonPi0R->SetFillColor(kAzure);
826 hndCVNCCNonPi0R->SetLineWidth(1);;
827 hndCVNCCNonPi0R->SetMarkerStyle(21);
828 hndCVNCCNonPi0R->SetMarkerColor(kAzure);
829 hs4->Add(hndCVNCCNonPi0R);
832 hndCVNNCNonPi0R->SetLineColor(kBlack);
833 hndCVNNCNonPi0R->SetFillColor(kCyan);
834 hndCVNNCNonPi0R->SetLineWidth(1);;
835 hndCVNNCNonPi0R->SetMarkerStyle(21);
836 hndCVNNCNonPi0R->SetMarkerColor(kCyan);
837 hs4->Add(hndCVNNCNonPi0R);
840 hndCVNNCPi0R->SetLineColor(kBlack);
841 hndCVNNCPi0R->SetFillColor(
kGreen);
842 hndCVNNCPi0R->SetLineWidth(1);;
843 hndCVNNCPi0R->SetMarkerStyle(21);
844 hndCVNNCPi0R->SetMarkerColor(
kGreen);
845 hs4->Add(hndCVNNCPi0R);
848 TText T4; T4.SetTextFont(42); T4.SetTextAlign(21);
849 hs4->Draw(
"hist"); T4.DrawTextNDC(.5,.95,
"");
850 hs4->GetYaxis()->SetTitle(
"Events/ 8e20 POT");
851 hs4->GetXaxis()->SetTitle(
"Prong CVN Gamma ID");
852 hs4->SetMaximum(10000);
853 hs4->GetXaxis()->SetRangeUser(0,1);
857 hndCVNMC3PR->SetLineColor(
kRed);
858 hndCVNMC3PR->SetLineWidth(3);
860 hndNumiCVN2PRData->SetMarkerSize(3);
863 hndCVNMC3PR->Draw(
"same");
865 TLegend * leg4 =
new TLegend(0.55,0.5,0.65,0.82);
867 leg4->SetTextFont(42);
868 leg4->SetTextSize(0.05);
869 leg4->SetFillStyle(0);
871 leg4->AddEntry(hndCVNMC3PR,
"Selection",
"kFullDotSmall");
872 leg4->AddEntry(hndCVNCCPi0R,
"CCPi0",
"kFullDotSmall");
873 leg4->AddEntry(hndCVNCCNonPi0R,
"CCNonPi0",
"kFullDotSmall");
874 leg4->AddEntry(hndCVNNCPi0LR,
"NC w/Pi0 < 0.3 E",
"kFullDotSmall");
875 leg4->AddEntry(hndCVNNCPi0R,
"NC w/Pi0 > 0.3 E",
"kFullDotSmall");
876 leg4->AddEntry(hndCVNNCNonPi0R,
"NCNonPi0",
"kFullDotSmall");
881 cs4->SaveAs(outDir2+CurrentFile[5].
name.c_str()+CurrentDate[0].
date.c_str()+
".png");
888 TCanvas *cs5 =
new TCanvas(
"cs5",
"cs5",10,10,855,500);
889 THStack *hs5 =
new THStack(
"hs5",
"");
901 hnd2CVNNCPi0LR->SetLineColor(kBlack);
902 hnd2CVNNCPi0LR->SetFillColor(
kOrange);
903 hnd2CVNNCPi0LR->SetLineWidth(1);;
904 hnd2CVNNCPi0LR->SetMarkerStyle(21);
905 hnd2CVNNCPi0LR->SetMarkerColor(
kOrange);
906 hs5->Add(hnd2CVNNCPi0LR);
909 hnd2CVNCCPi0R->SetLineColor(kBlack);
910 hnd2CVNCCPi0R->SetFillColor(kMagenta);
911 hnd2CVNCCPi0R->SetLineWidth(1);;
912 hnd2CVNCCPi0R->SetMarkerStyle(21);
913 hnd2CVNCCPi0R->SetMarkerColor(kMagenta);
915 hs5->Add(hnd2CVNCCPi0R);
917 hnd2CVNCCNonPi0R->SetLineColor(kBlack);
918 hnd2CVNCCNonPi0R->SetFillColor(kAzure);
919 hnd2CVNCCNonPi0R->SetLineWidth(1);;
920 hnd2CVNCCNonPi0R->SetMarkerStyle(21);
921 hnd2CVNCCNonPi0R->SetMarkerColor(kAzure);
922 hs5->Add(hnd2CVNCCNonPi0R);
925 hnd2CVNNCNonPi0R->SetLineColor(kBlack);
926 hnd2CVNNCNonPi0R->SetFillColor(kCyan);
927 hnd2CVNNCNonPi0R->SetLineWidth(1);;
928 hnd2CVNNCNonPi0R->SetMarkerStyle(21);
929 hnd2CVNNCNonPi0R->SetMarkerColor(kCyan);
930 hs5->Add(hnd2CVNNCNonPi0R);
933 hnd2CVNNCPi0R->SetLineColor(kBlack);
934 hnd2CVNNCPi0R->SetFillColor(
kGreen);
935 hnd2CVNNCPi0R->SetLineWidth(1);;
936 hnd2CVNNCPi0R->SetMarkerStyle(21);
937 hnd2CVNNCPi0R->SetMarkerColor(
kGreen);
938 hs5->Add(hnd2CVNNCPi0R);
941 TText T5; T5.SetTextFont(42); T5.SetTextAlign(21);
942 hs5->Draw(
"hist"); T5.DrawTextNDC(.5,.95,
"");
943 hs5->GetYaxis()->SetTitle(
"Events/ 8e20 POT");
944 hs5->GetXaxis()->SetTitle(
"Prong 2 CVN Gamma ID");
945 hs5->SetMaximum(14500);
946 hs5->GetXaxis()->SetRangeUser(0,1);
950 hnd2CVNMC3PR->SetLineColor(
kRed);
951 hnd2CVNMC3PR->SetLineWidth(3);
953 hndNumi2CVN2PRData->SetMarkerSize(3);
955 hnd2CVNMC3PR->Draw(
"same");
957 TLegend * leg5 =
new TLegend(0.55,0.5,0.65,0.82);
959 leg5->SetTextFont(42);
960 leg5->SetTextSize(0.05);
961 leg5->SetFillStyle(0);
963 leg5->AddEntry(hnd2CVNMC3PR,
"Selection",
"kFullDotSmall");
964 leg5->AddEntry(hnd2CVNCCPi0R,
"CCPi0",
"kFullDotSmall");
965 leg5->AddEntry(hnd2CVNCCNonPi0R,
"CCNonPi0",
"kFullDotSmall");
966 leg5->AddEntry(hnd2CVNNCPi0LR,
"NC w/Pi0 < 0.3 E",
"kFullDotSmall");
967 leg5->AddEntry(hnd2CVNNCPi0R,
"NC w/Pi0 > 0.3 E",
"kFullDotSmall");
968 leg5->AddEntry(hnd2CVNNCNonPi0R,
"NCNonPi0",
"kFullDotSmall");
973 cs5->SaveAs(outDir2+CurrentFile[6].
name.c_str()+CurrentDate[0].
date.c_str()+
".png");
980 TCanvas *cs6 =
new TCanvas(
"cs6",
"cs6",10,10,855,500);
981 THStack *hs6 =
new THStack(
"hs6",
"");
993 hndNCPi0LRem->SetLineColor(kBlack);
994 hndNCPi0LRem->SetFillColor(
kOrange);
995 hndNCPi0LRem->SetLineWidth(1);;
996 hndNCPi0LRem->SetMarkerStyle(21);
997 hndNCPi0LRem->SetMarkerColor(
kOrange);
998 hs6->Add(hndNCPi0LRem);
1001 hndCCPi0Rem->SetLineColor(kBlack);
1002 hndCCPi0Rem->SetFillColor(kMagenta);
1003 hndCCPi0Rem->SetLineWidth(1);;
1004 hndCCPi0Rem->SetMarkerStyle(21);
1005 hndCCPi0Rem->SetMarkerColor(kMagenta);
1007 hs6->Add(hndCCPi0Rem);
1009 hndCCNonPi0Rem->SetLineColor(kBlack);
1010 hndCCNonPi0Rem->SetFillColor(kAzure);
1011 hndCCNonPi0Rem->SetLineWidth(1);;
1012 hndCCNonPi0Rem->SetMarkerStyle(21);
1013 hndCCNonPi0Rem->SetMarkerColor(kAzure);
1014 hs6->Add(hndCCNonPi0Rem);
1017 hndNCNonPi0Rem->SetLineColor(kBlack);
1018 hndNCNonPi0Rem->SetFillColor(kCyan);
1019 hndNCNonPi0Rem->SetLineWidth(1);;
1020 hndNCNonPi0Rem->SetMarkerStyle(21);
1021 hndNCNonPi0Rem->SetMarkerColor(kCyan);
1022 hs6->Add(hndNCNonPi0Rem);
1025 hndNCPi0Rem->SetLineColor(kBlack);
1026 hndNCPi0Rem->SetFillColor(
kGreen);
1027 hndNCPi0Rem->SetLineWidth(1);;
1028 hndNCPi0Rem->SetMarkerStyle(21);
1029 hndNCPi0Rem->SetMarkerColor(
kGreen);
1030 hs6->Add(hndNCPi0Rem);
1033 TText T6; T6.SetTextFont(42); T6.SetTextAlign(21);
1034 hs6->Draw(
"hist"); T.DrawTextNDC(.5,.95,
"");
1035 hs6->GetYaxis()->SetTitle(
"Events/ 8e20 POT");
1036 hs6->GetXaxis()->SetTitle(
"RemID");
1037 hs6->SetMaximum(17000);
1038 hs6->GetXaxis()->SetRangeUser(0,1);
1042 hndMC3PRem->SetLineColor(
kRed);
1043 hndMC3PRem->SetLineWidth(3);
1045 hndNumiRemIDData->SetMarkerSize(10);
1048 hndMC3PRem->Draw(
"same");
1050 TLegend * leg6 =
new TLegend(0.55,0.5,0.65,0.82);
1052 leg6->SetTextFont(42);
1053 leg6->SetTextSize(0.05);
1054 leg6->SetFillStyle(0);
1056 leg6->AddEntry(hndMC3PRem,
"Selection",
"kFullDotSmall");
1057 leg6->AddEntry(hndCCPi0Rem,
"CCPi0",
"kFullDotSmall");
1058 leg6->AddEntry(hndCCNonPi0Rem,
"CCNonPi0",
"kFullDotSmall");
1059 leg6->AddEntry(hndNCPi0LRem,
"NC w/Pi0 < 0.3 E",
"kFullDotSmall");
1060 leg6->AddEntry(hndNCPi0Rem,
"NC w/Pi0 > 0.3 E",
"kFullDotSmall");
1061 leg6->AddEntry(hndNCNonPi0Rem,
"NCNonPi0",
"kFullDotSmall");
1066 cs6->SaveAs(outDir2+CurrentFile[7].
name.c_str()+CurrentDate[0].
date.c_str()+
".png");
1072 TCanvas *cs7 =
new TCanvas(
"cs7",
"cs7",10,10,855,500);
1073 THStack *hs7 =
new THStack(
"hs7",
"");
1085 hndNCPi0Prel->SetLineColor(kBlack);
1086 hndNCPi0Prel->SetFillColor(
kGreen);
1087 hndNCPi0Prel->SetLineWidth(1);;
1088 hndNCPi0Prel->SetMarkerStyle(21);
1089 hndNCPi0Prel->SetMarkerColor(
kGreen);
1090 hs7->Add(hndNCPi0Prel);
1093 hndNCPi0LPrel->SetLineColor(kBlack);
1094 hndNCPi0LPrel->SetFillColor(
kOrange);
1095 hndNCPi0LPrel->SetLineWidth(1);;
1096 hndNCPi0LPrel->SetMarkerStyle(21);
1097 hndNCPi0LPrel->SetMarkerColor(
kOrange);
1098 hs7->Add(hndNCPi0LPrel);
1102 hndNCNonPi0Prel->SetLineColor(kBlack);
1103 hndNCNonPi0Prel->SetFillColor(kCyan);
1104 hndNCNonPi0Prel->SetLineWidth(1);;
1105 hndNCNonPi0Prel->SetMarkerStyle(21);
1106 hndNCNonPi0Prel->SetMarkerColor(kCyan);
1107 hs7->Add(hndNCNonPi0Prel);
1110 hndCCPi0Prel->SetLineColor(kBlack);
1111 hndCCPi0Prel->SetFillColor(kMagenta);
1112 hndCCPi0Prel->SetLineWidth(1);;
1113 hndCCPi0Prel->SetMarkerStyle(21);
1114 hndCCPi0Prel->SetMarkerColor(kMagenta);
1116 hs7->Add(hndCCPi0Prel);
1119 hndCCNonPi0Prel->SetLineColor(kBlack);
1120 hndCCNonPi0Prel->SetFillColor(kAzure);
1121 hndCCNonPi0Prel->SetLineWidth(1);;
1122 hndCCNonPi0Prel->SetMarkerStyle(21);
1123 hndCCNonPi0Prel->SetMarkerColor(kAzure);
1124 hs7->Add(hndCCNonPi0Prel);
1127 TText T7; T7.SetTextFont(42); T7.SetTextAlign(21);
1128 hs7->Draw(
"hist"); T7.DrawTextNDC(.5,.95,
"");
1129 hs7->GetYaxis()->SetTitle(
"Events/ 8e20 POT");
1130 hs7->GetXaxis()->SetTitle(
"#pi^{0} True Energy [GeV]");
1131 hs7->SetMaximum(12000);
1136 hndMC3PPrel->SetLineColor(
kRed);
1137 hndMC3PPrel->SetLineWidth(3);
1138 hndMC3PPrel->Draw(
"same");
1139 TLegend * leg7 =
new TLegend(0.55,0.5,0.65,0.82);
1141 leg7->SetTextFont(42);
1142 leg7->SetTextSize(0.05);
1143 leg7->SetFillStyle(0);
1144 leg7->AddEntry(hndMC3PPrel,
"Selection",
"kFullDotSmall");
1148 leg7->AddEntry(hndNCPi0Prel,
"NC w/Pi0 > 0.3 E",
"kFullDotSmall");
1149 leg7->AddEntry(hndNCPi0LPrel,
"NC w/Pi0 < 0.3 E",
"kFullDotSmall");
1150 leg7->AddEntry(hndNCNonPi0Prel,
"NCNonPi0",
"kFullDotSmall");
1151 leg7->AddEntry(hndCCPi0Prel,
"CCPi0",
"kFullDotSmall");
1152 leg7->AddEntry(hndCCNonPi0Prel,
"CCNonPi0",
"kFullDotSmall");
1157 cs7->SaveAs(outDir2+CurrentFile[8].
name.c_str()+CurrentDate[0].
date.c_str()+
".png");
1164 TCanvas *cs8 =
new TCanvas(
"cs8",
"cs8",10,10,855,500);
1165 THStack *hs8 =
new THStack(
"hs8",
"");
1177 hndNCPi0PrelR->SetLineColor(kBlack);
1178 hndNCPi0PrelR->SetFillColor(
kGreen);
1179 hndNCPi0PrelR->SetLineWidth(1);;
1180 hndNCPi0PrelR->SetMarkerStyle(21);
1181 hndNCPi0PrelR->SetMarkerColor(
kGreen);
1182 hs8->Add(hndNCPi0PrelR);
1185 hndNCPi0LPrelR->SetLineColor(kBlack);
1186 hndNCPi0LPrelR->SetFillColor(
kOrange);
1187 hndNCPi0LPrelR->SetLineWidth(1);;
1188 hndNCPi0LPrelR->SetMarkerStyle(21);
1189 hndNCPi0LPrelR->SetMarkerColor(
kOrange);
1190 hs8->Add(hndNCPi0LPrelR);
1194 hndNCNonPi0PrelR->SetLineColor(kBlack);
1195 hndNCNonPi0PrelR->SetFillColor(kCyan);
1196 hndNCNonPi0PrelR->SetLineWidth(1);;
1197 hndNCNonPi0PrelR->SetMarkerStyle(21);
1198 hndNCNonPi0PrelR->SetMarkerColor(kCyan);
1199 hs8->Add(hndNCNonPi0PrelR);
1202 hndCCPi0PrelR->SetLineColor(kBlack);
1203 hndCCPi0PrelR->SetFillColor(kMagenta);
1204 hndCCPi0PrelR->SetLineWidth(1);;
1205 hndCCPi0PrelR->SetMarkerStyle(21);
1206 hndCCPi0PrelR->SetMarkerColor(kMagenta);
1208 hs8->Add(hndCCPi0PrelR);
1211 hndCCNonPi0PrelR->SetLineColor(kBlack);
1212 hndCCNonPi0PrelR->SetFillColor(kAzure);
1213 hndCCNonPi0PrelR->SetLineWidth(1);;
1214 hndCCNonPi0PrelR->SetMarkerStyle(21);
1215 hndCCNonPi0PrelR->SetMarkerColor(kAzure);
1216 hs8->Add(hndCCNonPi0PrelR);
1219 TText T8; T8.SetTextFont(42); T8.SetTextAlign(21);
1220 hs8->Draw(
"hist"); T8.DrawTextNDC(.5,.95,
"");
1221 hs8->GetYaxis()->SetTitle(
"Events/ 8e20 POT");
1222 hs8->GetXaxis()->SetTitle(
"#pi^{0} True Energy [GeV]");
1223 hs8->SetMaximum(9800);
1228 hndMC3PPrelR->SetLineColor(
kRed);
1229 hndMC3PPrelR->SetLineWidth(1);;
1230 hndMC3PPrelR->Draw(
"same");
1231 TLegend * leg8 =
new TLegend(0.55,0.5,0.65,0.82);
1233 leg8->SetTextFont(42);
1234 leg8->SetTextSize(0.05);
1235 leg8->SetFillStyle(0);
1236 leg8->AddEntry(hndMC3PPrel,
"Selection",
"kFullDotSmall");
1240 leg8->AddEntry(hndNCPi0PrelR,
"NC w/Pi0 > 0.3 E",
"kFullDotSmall");
1241 leg8->AddEntry(hndNCPi0LPrelR,
"NC w/Pi0 < 0.3 E",
"kFullDotSmall");
1242 leg8->AddEntry(hndNCNonPi0PrelR,
"NCNonPi0",
"kFullDotSmall");
1243 leg8->AddEntry(hndCCPi0PrelR,
"CCPi0",
"kFullDotSmall");
1244 leg8->AddEntry(hndCCNonPi0PrelR,
"CCNonPi0",
"kFullDotSmall");
1249 cs8->SaveAs(outDir2+CurrentFile[9].
name.c_str()+CurrentDate[0].
date.c_str()+
".png");
1255 TCanvas *cs9 =
new TCanvas(
"cs9",
"cs9",10,10,855,500);
1256 THStack *hs9 =
new THStack(
"hs9",
"");
1268 hndNCPi0PrelCVN->SetLineColor(kBlack);
1269 hndNCPi0PrelCVN->SetFillColor(
kGreen);
1270 hndNCPi0PrelCVN->SetLineWidth(1);;
1271 hndNCPi0PrelCVN->SetMarkerStyle(21);
1272 hndNCPi0PrelCVN->SetMarkerColor(
kGreen);
1273 hs9->Add(hndNCPi0PrelCVN);
1276 hndNCPi0LPrelCVN->SetLineColor(kBlack);
1277 hndNCPi0LPrelCVN->SetFillColor(
kOrange);
1278 hndNCPi0LPrelCVN->SetLineWidth(1);;
1279 hndNCPi0LPrelCVN->SetMarkerStyle(21);
1280 hndNCPi0LPrelCVN->SetMarkerColor(
kOrange);
1281 hs9->Add(hndNCPi0LPrelCVN);
1285 hndNCNonPi0PrelCVN->SetLineColor(kBlack);
1286 hndNCNonPi0PrelCVN->SetFillColor(kCyan);
1287 hndNCNonPi0PrelCVN->SetLineWidth(1);;
1288 hndNCNonPi0PrelCVN->SetMarkerStyle(21);
1289 hndNCNonPi0PrelCVN->SetMarkerColor(kCyan);
1290 hs9->Add(hndNCNonPi0PrelCVN);
1293 hndCCPi0PrelCVN->SetLineColor(kBlack);
1294 hndCCPi0PrelCVN->SetFillColor(kMagenta);
1295 hndCCPi0PrelCVN->SetLineWidth(1);;
1296 hndCCPi0PrelCVN->SetMarkerStyle(21);
1297 hndCCPi0PrelCVN->SetMarkerColor(kMagenta);
1299 hs9->Add(hndCCPi0PrelCVN);
1302 hndCCNonPi0PrelCVN->SetLineColor(kBlack);
1303 hndCCNonPi0PrelCVN->SetFillColor(kAzure);
1304 hndCCNonPi0PrelCVN->SetLineWidth(1);;
1305 hndCCNonPi0PrelCVN->SetMarkerStyle(21);
1306 hndCCNonPi0PrelCVN->SetMarkerColor(kAzure);
1307 hs9->Add(hndCCNonPi0PrelCVN);
1310 TText T9; T9.SetTextFont(42); T9.SetTextAlign(21);
1311 hs9->Draw(
"hist"); T9.DrawTextNDC(.5,.95,
"");
1312 hs9->GetYaxis()->SetTitle(
"Events/ 8e20 POT");
1313 hs9->GetXaxis()->SetTitle(
"#pi^{0} True Energy [GeV]");
1314 hs9->SetMaximum(4000);
1319 hndMC3PPrelCVN->SetLineColor(
kRed);
1320 hndMC3PPrelCVN->SetLineWidth(3);
1321 hndMC3PPrelCVN->Draw(
"same");
1322 TLegend * leg9 =
new TLegend(0.55,0.5,0.65,0.82);
1324 leg9->SetTextFont(42);
1325 leg9->SetTextSize(0.05);
1326 leg9->SetFillStyle(0);
1327 leg9->AddEntry(hndMC3PPrelCVN,
"Selection",
"kFullDotSmall");
1331 leg9->AddEntry(hndNCPi0PrelCVN,
"NC w/Pi0 > 0.3 E",
"kFullDotSmall");
1332 leg9->AddEntry(hndNCPi0LPrelCVN,
"NC w/Pi0 < 0.3 E",
"kFullDotSmall");
1333 leg9->AddEntry(hndNCNonPi0PrelCVN,
"NCNonPi0",
"kFullDotSmall");
1334 leg9->AddEntry(hndCCPi0PrelCVN,
"CCPi0",
"kFullDotSmall");
1335 leg9->AddEntry(hndCCNonPi0PrelCVN,
"CCNonPi0",
"kFullDotSmall");
1340 cs9->SaveAs(outDir2+CurrentFile[10].
name.c_str()+CurrentDate[0].
date.c_str()+
".png");
1345 TCanvas *cs9b =
new TCanvas(
"cs9b",
"cs9b",10,10,855,500);
1346 THStack *hs9b =
new THStack(
"hs9b",
"");
1358 hndNCPi0PrelCVN2Cuts->SetLineColor(kBlack);
1359 hndNCPi0PrelCVN2Cuts->SetFillColor(
kGreen);
1360 hndNCPi0PrelCVN2Cuts->SetLineWidth(1);;
1361 hndNCPi0PrelCVN2Cuts->SetMarkerStyle(21);
1362 hndNCPi0PrelCVN2Cuts->SetMarkerColor(
kGreen);
1363 hs9b->Add(hndNCPi0PrelCVN2Cuts);
1366 hndNCPi0LPrelCVN2Cuts->SetLineColor(kBlack);
1367 hndNCPi0LPrelCVN2Cuts->SetFillColor(
kOrange);
1368 hndNCPi0LPrelCVN2Cuts->SetLineWidth(1);;
1369 hndNCPi0LPrelCVN2Cuts->SetMarkerStyle(21);
1370 hndNCPi0LPrelCVN2Cuts->SetMarkerColor(
kOrange);
1371 hs9b->Add(hndNCPi0LPrelCVN2Cuts);
1375 hndNCNonPi0PrelCVN2Cuts->SetLineColor(kBlack);
1376 hndNCNonPi0PrelCVN2Cuts->SetFillColor(kCyan);
1377 hndNCNonPi0PrelCVN2Cuts->SetLineWidth(1);;
1378 hndNCNonPi0PrelCVN2Cuts->SetMarkerStyle(21);
1379 hndNCNonPi0PrelCVN2Cuts->SetMarkerColor(kCyan);
1380 hs9b->Add(hndNCNonPi0PrelCVN2Cuts);
1383 hndCCPi0PrelCVN2Cuts->SetLineColor(kBlack);
1384 hndCCPi0PrelCVN2Cuts->SetFillColor(kMagenta);
1385 hndCCPi0PrelCVN2Cuts->SetLineWidth(1);;
1386 hndCCPi0PrelCVN2Cuts->SetMarkerStyle(21);
1387 hndCCPi0PrelCVN2Cuts->SetMarkerColor(kMagenta);
1389 hs9b->Add(hndCCPi0PrelCVN2Cuts);
1392 hndCCNonPi0PrelCVN2Cuts->SetLineColor(kBlack);
1393 hndCCNonPi0PrelCVN2Cuts->SetFillColor(kAzure);
1394 hndCCNonPi0PrelCVN2Cuts->SetLineWidth(1);;
1395 hndCCNonPi0PrelCVN2Cuts->SetMarkerStyle(21);
1396 hndCCNonPi0PrelCVN2Cuts->SetMarkerColor(kAzure);
1397 hs9b->Add(hndCCNonPi0PrelCVN2Cuts);
1400 TText T91; T91.SetTextFont(42); T91.SetTextAlign(21);
1401 hs9b->Draw(
"hist"); T9.DrawTextNDC(.5,.95,
"");
1402 hs9b->GetYaxis()->SetTitle(
"Events/ 8e20 POT");
1403 hs9b->GetXaxis()->SetTitle(
"#pi^{0} True Energy [GeV]");
1404 hs9b->SetMaximum(1100);
1409 hndMC3PPrelCVN2Cuts->SetLineColor(
kRed);
1410 hndMC3PPrelCVN2Cuts->SetLineWidth(3);
1412 hndMC3PPrelCVN2Cuts->Draw(
"same");
1413 TLegend * leg9b =
new TLegend(0.55,0.5,0.65,0.82);
1415 leg9b->SetTextFont(42);
1416 leg9b->SetTextSize(0.05);
1417 leg9b->SetFillStyle(0);
1418 leg9b->AddEntry(hndMC3PPrelCVN2Cuts,
"Selection",
"kFullDotSmall");
1422 leg9b->AddEntry(hndNCPi0PrelCVN2Cuts,
"NC w/Pi0 > 0.3 E",
"kFullDotSmall");
1423 leg9b->AddEntry(hndNCPi0LPrelCVN2Cuts,
"NC w/Pi0 < 0.3 E",
"kFullDotSmall");
1424 leg9b->AddEntry(hndNCNonPi0PrelCVN2Cuts,
"NCNonPi0",
"kFullDotSmall");
1425 leg9b->AddEntry(hndCCPi0PrelCVN2Cuts,
"CCPi0",
"kFullDotSmall");
1426 leg9b->AddEntry(hndCCNonPi0PrelCVN2Cuts,
"CCNonPi0",
"kFullDotSmall");
1428 leg9b->Draw(
"same");
1431 cs9b->SaveAs(outDir2+CurrentFile[11].
name.c_str()+CurrentDate[0].
date.c_str()+
".png");
1435 TCanvas *cs10 =
new TCanvas(
"cs14",
"cs14",10,10,855,500);
1437 hnd2dplot2prong->Draw(
"colz");
1438 hnd2dplot2prong->GetYaxis()->SetTitle(
"#pi^{0} True Energy [GeV]");
1439 hnd2dplot2prong->GetXaxis()->SetTitle(
"#pi^{0} Reconstructed Mass [GeV]");
1440 hnd2dplot2prong->SetTitleOffset(0.9,
"Y");
1444 cs10->SaveAs(outDir2+CurrentFile[12].
name.c_str()+CurrentDate[0].
date.c_str()+
".png");
1449 TCanvas *cs10b =
new TCanvas(
"cs10b",
"cs10b",10,10,855,500);
1452 hnd2dplot2prong->Draw(
"colz");
1453 hnd2dplot2prong->GetYaxis()->SetTitle(
"#pi^{0} True Energy [GeV]");
1454 hnd2dplot2prong->GetXaxis()->SetTitle(
"#pi^{0} Reconstructed Mass [GeV]");
1455 hnd2dplot2prong->SetTitleOffset(0.9,
"Y");
1459 cs10b->SaveAs(outDir2+CurrentFile[15].
name.c_str()+CurrentDate[0].
date.c_str()+
".png");
1465 TCanvas *cs13 =
new TCanvas(
"cs13",
"cs13",10,10,855,500);
1468 hndMC3P1->SetLineColor(
kRed);
1473 hndCCPi01->SetLineColor(kMagenta+2);
1478 hndCCNonPi01->SetLineColor(kAzure+2);
1483 hndNCPi0L1->SetLineColor(
kOrange+2);
1488 hndNCPi01->SetLineColor(
kGreen+2);
1494 hndNCNonPi01->SetLineColor(kCyan-6);
1499 hndNumiMassdata2CVN1->SetMarkerSize(10);
1502 hndMC3P1->Draw(
"hist E1");
1504 hndCCPi01->Draw(
"hist E1 same");
1505 hndCCNonPi01->Draw(
"hist E1 same");
1506 hndNCPi0L1->Draw(
"hist E1 same");
1507 hndNCPi01->Draw(
"hist E1 same");
1508 hndNCNonPi01->Draw(
"hist E1 same");
1509 hndMC3P1->GetYaxis()->SetTitle(
"Events/8e20 POT");
1510 hndMC3P1->GetXaxis()->SetRangeUser(0.,0.4);
1511 hndMC3P1->SetTitleOffset(0.95,
"Y");
1514 TLegend * leg13 =
new TLegend(0.62,0.5,0.73,0.82);
1515 leg13->SetTextFont(42);
1516 leg13->SetTextSize(0.05);
1518 leg13->AddEntry(hndMC3P1,
"Selection",
"l");
1519 leg13->AddEntry(hndCCPi01,
"CCPi0",
"l");
1520 leg13->AddEntry(hndCCNonPi01,
"CCNonPi0",
"l");
1521 leg13->AddEntry(hndNCPi0L1,
"NC w/Pi0 < 0.3 E",
"l");
1522 leg13->AddEntry(hndNCPi01,
"NC w/Pi0 > 0.3 E",
"l");
1523 leg13->AddEntry(hndNCNonPi01,
"NCNonPi0",
"l");
1525 leg13->Draw(
"same");
1528 cs13->SaveAs(outDir2+CurrentFile[13].
name.c_str()+CurrentDate[0].
date.c_str()+
".png");
1534 TCanvas *cs13b =
new TCanvas(
"cs13b",
"cs13b",10,10,855,500);
1537 hndMC3P1Cut2->SetLineColor(
kRed);
1542 hndCCPi01Cut2->SetLineColor(kMagenta+2);
1547 hndCCNonPi01Cut2->SetLineColor(kAzure+2);
1552 hndNCPi0L1Cut2->SetLineColor(
kOrange+2);
1557 hndNCPi01Cut2->SetLineColor(
kGreen+2);
1563 hndNCNonPi01Cut2->SetLineColor(kCyan-6);
1568 hndNumiMassdataCVN2->SetMarkerSize(3);
1571 hndMC3P1Cut2->Draw(
"hist E1");
1573 hndCCPi01Cut2->Draw(
"hist E1 same");
1574 hndCCNonPi01Cut2->Draw(
"hist E1 same");
1575 hndNCPi0L1Cut2->Draw(
"hist E1 same");
1576 hndNCPi01Cut2->Draw(
"hist E1 same");
1577 hndNCNonPi01Cut2->Draw(
"hist E1 same");
1578 hndMC3P1Cut2->GetYaxis()->SetTitle(
"Events/8e20 POT");
1579 hndMC3P1Cut2->GetXaxis()->SetRangeUser(0.,0.4);
1580 hndMC3P1Cut2->SetTitleOffset(0.95,
"Y");
1583 TLegend * leg13b =
new TLegend(0.62,0.5,0.73,0.82);
1584 leg13b->SetTextFont(42);
1585 leg13b->SetTextSize(0.05);
1587 leg13b->AddEntry(hndMC3P1Cut2,
"Selection",
"l");
1588 leg13b->AddEntry(hndCCPi01Cut2,
"CCPi0",
"l");
1589 leg13b->AddEntry(hndCCNonPi01Cut2,
"CCNonPi0",
"l");
1590 leg13b->AddEntry(hndNCPi0L1Cut2,
"NC w/Pi0 < 0.3 E",
"l");
1591 leg13b->AddEntry(hndNCPi01Cut2,
"NC w/Pi0 > 0.3 E",
"l");
1592 leg13b->AddEntry(hndNCNonPi01Cut2,
"NCNonPi0",
"l");
1594 leg13b->Draw(
"same");
1597 cs13b->SaveAs(outDir2+CurrentFile[14].
name.c_str()+CurrentDate[0].
date.c_str()+
".png");
void plotHist_SpectrumCVNID()
::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)