Functions
twodvtxcontplotter.C File Reference
#include "TLine.h"
#include "Rtypes.h"
#include "TCanvas.h"
#include "TFile.h"
#include "TGraph.h"
#include "TH1.h"
#include "TH2.h"
#include "TMath.h"
#include "TGaxis.h"
#include "TMultiGraph.h"
#include "TLegend.h"
#include "TLegendEntry.h"
#include "TLatex.h"
#include "TStyle.h"
#include "THStack.h"
#include "TPaveText.h"
#include "TList.h"
#include "TAttLine.h"
#include "TAttMarker.h"
#include <cmath>
#include <iostream>
#include <vector>
#include <list>
#include <sstream>
#include <string>
#include <fstream>
#include <iomanip>
#include "CAFAna/Core/Binning.h"
#include "CAFAna/Cuts/Cuts.h"
#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Core/SpectrumLoader.h"
#include "CAFAna/Vars/Vars.h"
#include "CAFAna/Core/Utilities.h"
#include "CAFAna/Analysis/Plots.h"
#include "CAFAna/Cuts/SpillCuts.h"
#include "3FlavorAna/Cuts/NumuCuts.h"
#include "CAFAna/Cuts/TruthCuts.h"
#include "3FlavorAna/Vars/NumuVars.h"
#include "CAFAna/Core/SpectrumLoaderBase.h"
#include "CAFAna/Analysis/Style.h"
#include "TColor.h"
#include "Utilities/rootlogon.C"

Go to the source code of this file.

Functions

void NDSimulation ()
 
void twodvtxcontplotter ()
 

Function Documentation

void NDSimulation ( )

Definition at line 60 of file twodvtxcontplotter.C.

References prelim.

Referenced by twodvtxcontplotter().

61 {
62  TLatex* prelim = new TLatex(.9, .95, "NOvA Simulation");
63  prelim->SetTextColor(kGray+1);
64  prelim->SetNDC();
65  prelim->SetTextSize(2/30.);
66  prelim->SetTextAlign(32);
67  TGaxis::SetMaxDigits(3);
68  prelim->Draw();
69 }
TLatex * prelim
Definition: Xsec_final.C:133
void twodvtxcontplotter ( )

Read the Spectrum Root file /////

Vertex

Vertex

Vertex

Vertex

Vertex

Vertex

How to Scale Histogram////

All Interactions////

Plot 1 vertex x-dir ///

Plot 2 vertex y-dir ///

Plot 3 vertex z-dir ///

Plot 4 Containment x-dir ///

Plot 5 Containment y-dir ///

Plot 6 Containment z-dir ///

CC Non pi0 Interactions////

Plot 1 vertex x-dir ///

Plot 2 vertex y-dir ///

Plot 3 vertex z-dir ///

Plot 4 Containment x-dir ///

Plot 5 Containment y-dir ///

Plot 6 Containment z-dir ///

NC Non pi0 Interactions////

Plot 1 vertex x-dir ///

Plot 2 vertex y-dir ///

Plot 3 vertex z-dir ///

Plot 4 Containment x-dir ///

Plot 5 Containment y-dir ///

Plot 6 Containment z-dir ///

CC Pi0 Interactions////

Plot 1 vertex x-dir ///

Plot 2 vertex y-dir ///

Plot 3 vertex z-dir ///

Plot 4 Containment x-dir ///

Plot 5 Containment y-dir ///

Plot 6 Containment z-dir ///

NC Pi0 Interactions////

Plot 1 vertex x-dir ///

Plot 2 vertex y-dir ///

Plot 3 vertex z-dir ///

Plot 4 Containment x-dir ///

Plot 5 Containment y-dir ///

Plot 6 Containment z-dir ///

NC Pi0 low energy Interactions////

Plot 1 vertex x-dir ///

Plot 2 vertex y-dir ///

Plot 3 vertex z-dir ///

Plot 4 Containment x-dir ///

Plot 5 Containment y-dir ///

Plot 6 Containment z-dir ///

Definition at line 71 of file twodvtxcontplotter.C.

References CurrentDate, CurrentFile, DateUp::date, Date, f1, File, ana::Spectrum::LoadFrom(), NDSimulation(), outDir2, and string.

72 {
73 
74 
75 
76 /*********************************************************************************************************************/
77 
78 ///Read the Spectrum Root file /////
79 
80  TFile *f1= new TFile ("/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/2dvtxcontstudyplots2dvtxcontsutyplotsdifferentinteractions.root","READ");
81  f1->cd();
82 
83  /////All interactions///////////
84 
85  ///Vertex
86  auto snd2dplot2prongx=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongx"));
87  auto snd2dplot2prongy=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongy"));
88  auto snd2dplot2prongz=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongz"));
89 
90  //Containment
91  auto snd2dplot2prongcontX=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongcontX"));
92  auto snd2dplot2prongcontY=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongcontY"));
93  auto snd2dplot2prongcontZ=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongcontZ"));
94 
95 /*********************************************************************************************************************/
96  /////CC Non pi0 Interactions///////////
97 
98  ///Vertex
99  auto snd2dplot2prongCCNonPi0x=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongCCNonPi0x"));
100  auto snd2dplot2prongCCNonPi0y=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongCCNonPi0y"));
101  auto snd2dplot2prongCCNonPi0z=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongCCNonPi0z"));
102 
103  //Containment
104  auto snd2dplot2prongCCNonPi0contX=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongCCNonPi0contX"));
105  auto snd2dplot2prongCCNonPi0contY=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongCCNonPi0contY"));
106  auto snd2dplot2prongCCNonPi0contZ=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongCCNonPi0contZ"));
107 /*********************************************************************************************************************/
108  /////NC Non pi0 Interactions///////////
109 
110  ///Vertex
111  auto snd2dplot2prongNCNonPi0x=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCNonPi0x"));
112  auto snd2dplot2prongNCNonPi0y=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCNonPi0y"));
113  auto snd2dplot2prongNCNonPi0z=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCNonPi0z"));
114 
115  //Containment
116  auto snd2dplot2prongNCNonPi0contX=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCNonPi0contX"));
117  auto snd2dplot2prongNCNonPi0contY=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCNonPi0contY"));
118  auto snd2dplot2prongNCNonPi0contZ=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCNonPi0contZ"));
119 
120 /*********************************************************************************************************************/
121  /////CC pi0 Interactions///////////
122 
123  ///Vertex
124  auto snd2dplot2prongCCPi0x=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongCCPi0x"));
125  auto snd2dplot2prongCCPi0y=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongCCPi0y"));
126  auto snd2dplot2prongCCPi0z=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongCCPi0z"));
127 
128  //Containment
129  auto snd2dplot2prongCCPi0contX=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongCCPi0contX"));
130  auto snd2dplot2prongCCPi0contY=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongCCPi0contY"));
131  auto snd2dplot2prongCCPi0contZ=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongCCPi0contZ"));
132 /*********************************************************************************************************************/
133  /////NC pi0 Interactions///////////
134 
135  ///Vertex
136  auto snd2dplot2prongNCPi0x=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCPi0x"));
137  auto snd2dplot2prongNCPi0y=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCPi0y"));
138  auto snd2dplot2prongNCPi0z=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCPi0z"));
139 
140  //Containment
141  auto snd2dplot2prongNCPi0contX=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCPi0contX"));
142  auto snd2dplot2prongNCPi0contY=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCPi0contY"));
143  auto snd2dplot2prongNCPi0contZ=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCPi0contZ"));
144 /*********************************************************************************************************************/
145  /////NC pi0 low energy Interactions///////////
146 
147  ///Vertex
148  auto snd2dplot2prongNCPi0Lx=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCPi0Lx"));
149  auto snd2dplot2prongNCPi0Ly=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCPi0Ly"));
150  auto snd2dplot2prongNCPi0Lz=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCPi0Lz"));
151 
152  //Containment
153  auto snd2dplot2prongNCPi0LcontX=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCPi0LcontX"));
154  auto snd2dplot2prongNCPi0LcontY=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCPi0LcontY"));
155  auto snd2dplot2prongNCPi0LcontZ=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCPi0LcontZ"));
156 
157 
158 
159 
160 /*********************************************************************************************************************/
161  ///How to Scale Histogram////
162 
163  // double myPOT=sndMC.POT();
164  const double myPOT = 8e20;
165 /*********************************************************************************************************************/
166 
167  //All interactions
168 
169  //Vertex
170  auto hnd2dplot2prongx = snd2dplot2prongx.ToTH2(myPOT);
171  auto hnd2dplot2prongy = snd2dplot2prongy.ToTH2(myPOT);
172  auto hnd2dplot2prongz = snd2dplot2prongz.ToTH2(myPOT);
173 
174  //Containment
175  auto hnd2dplot2prongcontX = snd2dplot2prongcontX.ToTH2(myPOT);
176  auto hnd2dplot2prongcontY = snd2dplot2prongcontY.ToTH2(myPOT);
177  auto hnd2dplot2prongcontZ = snd2dplot2prongcontZ.ToTH2(myPOT);
178  /*********************************************************************************************************************/
179  //CC Non pi0 Interactions///////
180 
181  //Vertex
182  auto hnd2dplot2prongCCNonPi0x = snd2dplot2prongCCNonPi0x.ToTH2(myPOT);
183  auto hnd2dplot2prongCCNonPi0y = snd2dplot2prongCCNonPi0y.ToTH2(myPOT);
184  auto hnd2dplot2prongCCNonPi0z = snd2dplot2prongCCNonPi0z.ToTH2(myPOT);
185 
186  //Containment
187  auto hnd2dplot2prongCCNonPi0contX = snd2dplot2prongCCNonPi0contX.ToTH2(myPOT);
188  auto hnd2dplot2prongCCNonPi0contY = snd2dplot2prongCCNonPi0contY.ToTH2(myPOT);
189  auto hnd2dplot2prongCCNonPi0contZ = snd2dplot2prongCCNonPi0contZ.ToTH2(myPOT);
190 
191  /*********************************************************************************************************************/
192  //NC Non pi0 Interactions///////
193 
194  //Vertex
195  auto hnd2dplot2prongNCNonPi0x = snd2dplot2prongNCNonPi0x.ToTH2(myPOT);
196  auto hnd2dplot2prongNCNonPi0y = snd2dplot2prongNCNonPi0y.ToTH2(myPOT);
197  auto hnd2dplot2prongNCNonPi0z = snd2dplot2prongNCNonPi0z.ToTH2(myPOT);
198 
199  //Containment
200  auto hnd2dplot2prongNCNonPi0contX = snd2dplot2prongNCNonPi0contX.ToTH2(myPOT);
201  auto hnd2dplot2prongNCNonPi0contY = snd2dplot2prongNCNonPi0contY.ToTH2(myPOT);
202  auto hnd2dplot2prongNCNonPi0contZ = snd2dplot2prongNCNonPi0contZ.ToTH2(myPOT);
203 
204 
205 /*********************************************************************************************************************/
206 
207  //CC Pi0 Interactions//////////////////
208 
209  //Vertex
210  auto hnd2dplot2prongCCPi0x = snd2dplot2prongCCPi0x.ToTH2(myPOT);
211  auto hnd2dplot2prongCCPi0y = snd2dplot2prongCCPi0y.ToTH2(myPOT);
212  auto hnd2dplot2prongCCPi0z = snd2dplot2prongCCPi0z.ToTH2(myPOT);
213 
214  //Containment
215  auto hnd2dplot2prongCCPi0contX = snd2dplot2prongCCPi0contX.ToTH2(myPOT);
216  auto hnd2dplot2prongCCPi0contY = snd2dplot2prongCCPi0contY.ToTH2(myPOT);
217  auto hnd2dplot2prongCCPi0contZ = snd2dplot2prongCCPi0contZ.ToTH2(myPOT);
218 
219 /*********************************************************************************************************************/
220 
221  //NC Pi0 Interactions//////////////////
222 
223  //Vertex
224  auto hnd2dplot2prongNCPi0x = snd2dplot2prongNCPi0x.ToTH2(myPOT);
225  auto hnd2dplot2prongNCPi0y = snd2dplot2prongNCPi0y.ToTH2(myPOT);
226  auto hnd2dplot2prongNCPi0z = snd2dplot2prongNCPi0z.ToTH2(myPOT);
227 
228  //Containment
229  auto hnd2dplot2prongNCPi0contX = snd2dplot2prongNCPi0contX.ToTH2(myPOT);
230  auto hnd2dplot2prongNCPi0contY = snd2dplot2prongNCPi0contY.ToTH2(myPOT);
231  auto hnd2dplot2prongNCPi0contZ = snd2dplot2prongNCPi0contZ.ToTH2(myPOT);
232 
233 /*********************************************************************************************************************/
234 
235  //NC Pi0 low energy Interactions//////////////////
236 
237  //Vertex
238  auto hnd2dplot2prongNCPi0Lx = snd2dplot2prongNCPi0Lx.ToTH2(myPOT);
239  auto hnd2dplot2prongNCPi0Ly = snd2dplot2prongNCPi0Ly.ToTH2(myPOT);
240  auto hnd2dplot2prongNCPi0Lz = snd2dplot2prongNCPi0Lz.ToTH2(myPOT);
241 
242  //Containment
243  auto hnd2dplot2prongNCPi0LcontX = snd2dplot2prongNCPi0LcontX.ToTH2(myPOT);
244  auto hnd2dplot2prongNCPi0LcontY = snd2dplot2prongNCPi0LcontY.ToTH2(myPOT);
245  auto hnd2dplot2prongNCPi0LcontZ = snd2dplot2prongNCPi0LcontZ.ToTH2(myPOT);
246 /*********************************************************************************************************************/
247 //****************************Updated Date and string path for saving files **********************************************************///
248  TString outDir2 = "/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/twodimvtxcontplots";
249  struct FileUp{
251  };
252  struct DateUp{
254  };
255 
256  const int File =37;
257  const FileUp CurrentFile[File] = {
258  {"/Vertex_X_Direction_AllInteractions"},
259  {"/Vertex_Y_Direction_AllInteractions"},
260  {"/Vertex_Z_Direction_AllInteractions"},
261  {"/Containment_X_Direction_AllInteractions"},
262  {"/Containment_Y_Direction_AllInteractions"},
263  {"/Containment_Z_Direction_AllInteractions"},
264 
265  {"/Vertex_X_Direction_CCNonPi0"},
266  {"/Vertex_Y_Direction_CCNonPi0"},
267  {"/Vertex_Z_Direction_CCNonPi0"},
268  {"/Containment_X_Direction_CCNonPi0"},
269  {"/Containment_Y_Direction_CCNonPi0"},
270  {"/Containment_Z_Direction_CCNonPi0"},
271 
272  {"/Vertex_X_Direction_NCNonPi0"},
273  {"/Vertex_Y_Direction_NCNonPi0"},
274  {"/Vertex_Z_Direction_NCNonPi0"},
275  {"/Containment_X_Direction_NCNonPi0"},
276  {"/Containment_Y_Direction_NCNonPi0"},
277  {"/Containment_Z_Direction_NCNonPi0"},
278 
279  {"/Vertex_X_Direction_CCPi0"},
280  {"/Vertex_Y_Direction_CCPi0"},
281  {"/Vertex_Z_Direction_CCPi0"},
282  {"/Containment_X_Direction_CCPi0"},
283  {"/Containment_Y_Direction_CCPi0"},
284  {"/Containment_Z_Direction_CCPi0"},
285 
286  {"/Vertex_X_Direction_NCPi0HighEnergy"},
287  {"/Vertex_Y_Direction_NCPi0HighEnergy"},
288  {"/Vertex_Z_Direction_NCPi0HighEnergy"},
289  {"/Containment_X_Direction_NCPi0HighEnergy"},
290  {"/Containment_Y_Direction_NCPi0HighEnergy"},
291  {"/Containment_Z_Direction_NCPi0HighEnergy"},
292 
293  {"/Vertex_X_Direction_NCPi0_LowEnergy"},
294  {"/Vertex_Y_Direction_NCPi0_LowEnergy"},
295  {"/Vertex_Z_Direction_NCPi0_LowEnergy"},
296  {"/Containment_X_Direction_NCPi0_LowEnergy"},
297  {"/Containment_Y_Direction_NCPi0_LowEnergy"},
298  {"/Containment_Z_Direction_NCPi0_Lowenergy"},
299  };
300 
301  const int Date = 1;
302  const DateUp CurrentDate[Date] = {
303  {"_03-30-18"},
304  };
305 
306 
307 ///All Interactions////
308 
309  /// Plot 1 vertex x-dir ///
310  TCanvas *cs1 = new TCanvas("cs1","cs1",10,10,855,500);
311  cs1->cd();
312  hnd2dplot2prongx->Draw("colz");
313  hnd2dplot2prongx->GetYaxis()->SetTitle("Vertex in x-dir [cm]");
314  hnd2dplot2prongx->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
315  hnd2dplot2prongx->SetTitleOffset(0.9,"Y");
316  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
317  NDSimulation();
318  cs1->SaveAs(outDir2+CurrentFile[0].name.c_str()+CurrentDate[0].date.c_str()+".png");
319 
320  /// Plot 2 vertex y-dir ///
321  TCanvas *cs2 = new TCanvas("cs2","cs2",10,10,855,500);
322  cs2->cd();
323  hnd2dplot2prongy->Draw("colz");
324  hnd2dplot2prongy->GetYaxis()->SetTitle("Vertex in y-dir [cm]");
325  hnd2dplot2prongy->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
326  hnd2dplot2prongy->SetTitleOffset(0.9,"Y");
327  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
328  NDSimulation();
329  cs2->SaveAs(outDir2+CurrentFile[1].name.c_str()+CurrentDate[0].date.c_str()+".png");
330 
331  /// Plot 3 vertex z-dir ///
332  TCanvas *cs3 = new TCanvas("cs3","cs3",10,10,855,500);
333  cs3->cd();
334  hnd2dplot2prongz->Draw("colz");
335  hnd2dplot2prongz->GetYaxis()->SetTitle("Vertex in z-dir [cm]");
336  hnd2dplot2prongz->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
337  hnd2dplot2prongz->SetTitleOffset(0.9,"Y");
338  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
339  NDSimulation();
340  cs3->SaveAs(outDir2+CurrentFile[2].name.c_str()+CurrentDate[0].date.c_str()+".png");
341 
342  /// Plot 4 Containment x-dir ///
343  TCanvas *cs4 = new TCanvas("cs4","cs4",10,10,855,500);
344  cs4->cd();
345  hnd2dplot2prongcontX->Draw("colz");
346  hnd2dplot2prongcontX->GetYaxis()->SetTitle("Containment in x-dir [cm]");
347  hnd2dplot2prongcontX->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
348  hnd2dplot2prongcontX->SetTitleOffset(0.9,"Y");
349  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
350  NDSimulation();
351  cs4->SaveAs(outDir2+CurrentFile[3].name.c_str()+CurrentDate[0].date.c_str()+".png");
352 
353  /// Plot 5 Containment y-dir ///
354  TCanvas *cs5 = new TCanvas("cs5","cs5",10,10,855,500);
355  cs5->cd();
356  hnd2dplot2prongcontY->Draw("colz");
357  hnd2dplot2prongcontY->GetYaxis()->SetTitle("Containment in y-dir [cm]");
358  hnd2dplot2prongcontY->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
359  hnd2dplot2prongcontY->SetTitleOffset(0.9,"Y");
360  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
361  NDSimulation();
362  cs5->SaveAs(outDir2+CurrentFile[4].name.c_str()+CurrentDate[0].date.c_str()+".png");
363 
364  /// Plot 6 Containment z-dir ///
365  TCanvas *cs6 = new TCanvas("cs6","cs6",10,10,855,500);
366  cs6->cd();
367  hnd2dplot2prongcontZ->Draw("colz");
368  hnd2dplot2prongcontZ->GetYaxis()->SetTitle("Containment in z-dir [cm]");
369  hnd2dplot2prongcontZ->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
370  hnd2dplot2prongcontZ->SetTitleOffset(0.9,"Y");
371  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
372  NDSimulation();
373  cs6->SaveAs(outDir2+CurrentFile[5].name.c_str()+CurrentDate[0].date.c_str()+".png");
374 /*********************************************************************************************************************/
375 /*********************************************************************************************************************/
376 
377 
378 ///CC Non pi0 Interactions////
379 
380  /// Plot 1 vertex x-dir ///
381  TCanvas *cs7 = new TCanvas("cs7","cs7",10,10,855,500);
382  cs7->cd();
383  hnd2dplot2prongCCNonPi0x->Draw("colz");
384  hnd2dplot2prongCCNonPi0x->GetYaxis()->SetTitle("Vertex in x-dir [cm]");
385  hnd2dplot2prongCCNonPi0x->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
386  hnd2dplot2prongCCNonPi0x->SetTitleOffset(0.9,"Y");
387  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
388  NDSimulation();
389  cs7->SaveAs(outDir2+CurrentFile[6].name.c_str()+CurrentDate[0].date.c_str()+".png");
390 
391  /// Plot 2 vertex y-dir ///
392  TCanvas *cs8 = new TCanvas("cs8","cs8",10,10,855,500);
393  cs8->cd();
394  hnd2dplot2prongCCNonPi0y->Draw("colz");
395  hnd2dplot2prongCCNonPi0y->GetYaxis()->SetTitle("Vertex in y-dir [cm]");
396  hnd2dplot2prongCCNonPi0y->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
397  hnd2dplot2prongCCNonPi0y->SetTitleOffset(0.9,"Y");
398  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
399  NDSimulation();
400  cs8->SaveAs(outDir2+CurrentFile[7].name.c_str()+CurrentDate[0].date.c_str()+".png");
401 
402  /// Plot 3 vertex z-dir ///
403  TCanvas *cs55 = new TCanvas("cs55","cs55",10,10,855,500);
404  cs55->cd();
405  hnd2dplot2prongCCNonPi0z->Draw("colz");
406  hnd2dplot2prongCCNonPi0z->GetYaxis()->SetTitle("Vertex in z-dir [cm]");
407  hnd2dplot2prongCCNonPi0z->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
408  hnd2dplot2prongCCNonPi0z->SetTitleOffset(0.9,"Y");
409  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
410  NDSimulation();
411  cs55->SaveAs(outDir2+CurrentFile[8].name.c_str()+CurrentDate[0].date.c_str()+".png");
412 
413  /// Plot 4 Containment x-dir ///
414  TCanvas *cs9 = new TCanvas("cs9","cs9",10,10,855,500);
415  cs9->cd();
416  hnd2dplot2prongCCNonPi0contX->Draw("colz");
417  hnd2dplot2prongCCNonPi0contX->GetYaxis()->SetTitle("Containment in x-dir [cm]");
418  hnd2dplot2prongCCNonPi0contX->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
419  hnd2dplot2prongCCNonPi0contX->SetTitleOffset(0.9,"Y");
420  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
421  NDSimulation();
422  cs9->SaveAs(outDir2+CurrentFile[9].name.c_str()+CurrentDate[0].date.c_str()+".png");
423 
424  /// Plot 5 Containment y-dir ///
425  TCanvas *cs10 = new TCanvas("cs10","cs10",10,10,855,500);
426  cs10->cd();
427  hnd2dplot2prongCCNonPi0contY->Draw("colz");
428  hnd2dplot2prongCCNonPi0contY->GetYaxis()->SetTitle("Containment in y-dir [cm]");
429  hnd2dplot2prongCCNonPi0contY->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
430  hnd2dplot2prongCCNonPi0contY->SetTitleOffset(0.9,"Y");
431  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
432  NDSimulation();
433  cs10->SaveAs(outDir2+CurrentFile[10].name.c_str()+CurrentDate[0].date.c_str()+".png");
434 
435  /// Plot 6 Containment z-dir ///
436  TCanvas *cs11 = new TCanvas("cs11","cs11",10,10,855,500);
437  cs11->cd();
438  hnd2dplot2prongCCNonPi0contZ->Draw("colz");
439  hnd2dplot2prongCCNonPi0contZ->GetYaxis()->SetTitle("Containment in z-dir [cm]");
440  hnd2dplot2prongCCNonPi0contZ->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
441  hnd2dplot2prongCCNonPi0contZ->SetTitleOffset(0.9,"Y");
442  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
443  NDSimulation();
444  cs11->SaveAs(outDir2+CurrentFile[11].name.c_str()+CurrentDate[0].date.c_str()+".png");
445 /*********************************************************************************************************************/
446 /*********************************************************************************************************************/
447 
448 
449 ///NC Non pi0 Interactions////
450 
451  /// Plot 1 vertex x-dir ///
452  TCanvas *cs12 = new TCanvas("cs12","cs12",10,10,855,500);
453  cs12->cd();
454  hnd2dplot2prongNCNonPi0x->Draw("colz");
455  hnd2dplot2prongNCNonPi0x->GetYaxis()->SetTitle("Vertex in x-dir [cm]");
456  hnd2dplot2prongNCNonPi0x->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
457  hnd2dplot2prongNCNonPi0x->SetTitleOffset(0.9,"Y");
458  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
459  NDSimulation();
460  cs12->SaveAs(outDir2+CurrentFile[12].name.c_str()+CurrentDate[0].date.c_str()+".png");
461 
462  /// Plot 2 vertex y-dir ///
463  TCanvas *cs13 = new TCanvas("cs13","cs13",10,10,855,500);
464  cs13->cd();
465  hnd2dplot2prongNCNonPi0y->Draw("colz");
466  hnd2dplot2prongNCNonPi0y->GetYaxis()->SetTitle("Vertex in y-dir [cm]");
467  hnd2dplot2prongNCNonPi0y->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
468  hnd2dplot2prongNCNonPi0y->SetTitleOffset(0.9,"Y");
469  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
470  NDSimulation();
471  cs13->SaveAs(outDir2+CurrentFile[13].name.c_str()+CurrentDate[0].date.c_str()+".png");
472 
473  /// Plot 3 vertex z-dir ///
474  TCanvas *cs14 = new TCanvas("cs14","cs14",10,10,855,500);
475  cs14->cd();
476  hnd2dplot2prongNCNonPi0z->Draw("colz");
477  hnd2dplot2prongNCNonPi0z->GetYaxis()->SetTitle("Vertex in z-dir [cm]");
478  hnd2dplot2prongNCNonPi0z->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
479  hnd2dplot2prongNCNonPi0z->SetTitleOffset(0.9,"Y");
480  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
481  NDSimulation();
482  cs14->SaveAs(outDir2+CurrentFile[14].name.c_str()+CurrentDate[0].date.c_str()+".png");
483 
484  /// Plot 4 Containment x-dir ///
485  TCanvas *cs15 = new TCanvas("cs15","cs15",10,10,855,500);
486  cs15->cd();
487  hnd2dplot2prongNCNonPi0contX->Draw("colz");
488  hnd2dplot2prongNCNonPi0contX->GetYaxis()->SetTitle("Containment in x-dir [cm]");
489  hnd2dplot2prongNCNonPi0contX->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
490  hnd2dplot2prongNCNonPi0contX->SetTitleOffset(0.9,"Y");
491  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
492  NDSimulation();
493  cs15->SaveAs(outDir2+CurrentFile[15].name.c_str()+CurrentDate[0].date.c_str()+".png");
494 
495  /// Plot 5 Containment y-dir ///
496  TCanvas *cs16 = new TCanvas("cs16","cs16",10,10,855,500);
497  cs16->cd();
498  hnd2dplot2prongNCNonPi0contY->Draw("colz");
499  hnd2dplot2prongNCNonPi0contY->GetYaxis()->SetTitle("Containment in y-dir [cm]");
500  hnd2dplot2prongNCNonPi0contY->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
501  hnd2dplot2prongNCNonPi0contY->SetTitleOffset(0.9,"Y");
502  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
503  NDSimulation();
504  cs16->SaveAs(outDir2+CurrentFile[16].name.c_str()+CurrentDate[0].date.c_str()+".png");
505 
506  /// Plot 6 Containment z-dir ///
507  TCanvas *cs17 = new TCanvas("cs17","cs17",10,10,855,500);
508  cs17->cd();
509  hnd2dplot2prongNCNonPi0contZ->Draw("colz");
510  hnd2dplot2prongNCNonPi0contZ->GetYaxis()->SetTitle("Containment in z-dir [cm]");
511  hnd2dplot2prongNCNonPi0contZ->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
512  hnd2dplot2prongNCNonPi0contZ->SetTitleOffset(0.9,"Y");
513  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
514  NDSimulation();
515  cs17->SaveAs(outDir2+CurrentFile[17].name.c_str()+CurrentDate[0].date.c_str()+".png");
516 /*********************************************************************************************************************/
517 ///CC Pi0 Interactions////
518 
519  /// Plot 1 vertex x-dir ///
520  TCanvas *cs18 = new TCanvas("cs18","cs18",10,10,855,500);
521  cs18->cd();
522  hnd2dplot2prongCCPi0x->Draw("colz");
523  hnd2dplot2prongCCPi0x->GetYaxis()->SetTitle("Vertex in x-dir [cm]");
524  hnd2dplot2prongCCPi0x->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
525  hnd2dplot2prongCCPi0x->SetTitleOffset(0.9,"Y");
526  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
527  NDSimulation();
528  cs18->SaveAs(outDir2+CurrentFile[18].name.c_str()+CurrentDate[0].date.c_str()+".png");
529 
530  /// Plot 2 vertex y-dir ///
531  TCanvas *cs19 = new TCanvas("cs19","cs19",10,10,855,500);
532  cs19->cd();
533  hnd2dplot2prongCCPi0y->Draw("colz");
534  hnd2dplot2prongCCPi0y->GetYaxis()->SetTitle("Vertex in y-dir [cm]");
535  hnd2dplot2prongCCPi0y->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
536  hnd2dplot2prongCCPi0y->SetTitleOffset(0.9,"Y");
537  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
538  NDSimulation();
539  cs19->SaveAs(outDir2+CurrentFile[19].name.c_str()+CurrentDate[0].date.c_str()+".png");
540 
541  /// Plot 3 vertex z-dir ///
542  TCanvas *cs20 = new TCanvas("cs20","cs20",10,10,855,500);
543  cs20->cd();
544  hnd2dplot2prongCCPi0z->Draw("colz");
545  hnd2dplot2prongCCPi0z->GetYaxis()->SetTitle("Vertex in z-dir [cm]");
546  hnd2dplot2prongCCPi0z->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
547  hnd2dplot2prongCCPi0z->SetTitleOffset(0.9,"Y");
548  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
549  NDSimulation();
550  cs20->SaveAs(outDir2+CurrentFile[20].name.c_str()+CurrentDate[0].date.c_str()+".png");
551 
552  /// Plot 4 Containment x-dir ///
553  TCanvas *cs21 = new TCanvas("cs21","cs21",10,10,855,500);
554  cs21->cd();
555  hnd2dplot2prongCCPi0contX->Draw("colz");
556  hnd2dplot2prongCCPi0contX->GetYaxis()->SetTitle("Containment in x-dir [cm]");
557  hnd2dplot2prongCCPi0contX->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
558  hnd2dplot2prongCCPi0contX->SetTitleOffset(0.9,"Y");
559  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
560  NDSimulation();
561  cs21->SaveAs(outDir2+CurrentFile[21].name.c_str()+CurrentDate[0].date.c_str()+".png");
562 
563  /// Plot 5 Containment y-dir ///
564  TCanvas *cs22 = new TCanvas("cs22","cs22",10,10,855,500);
565  cs22->cd();
566  hnd2dplot2prongCCPi0contY->Draw("colz");
567  hnd2dplot2prongCCPi0contY->GetYaxis()->SetTitle("Containment in y-dir [cm]");
568  hnd2dplot2prongCCPi0contY->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
569  hnd2dplot2prongCCPi0contY->SetTitleOffset(0.9,"Y");
570  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
571  NDSimulation();
572  cs22->SaveAs(outDir2+CurrentFile[22].name.c_str()+CurrentDate[0].date.c_str()+".png");
573 
574  /// Plot 6 Containment z-dir ///
575  TCanvas *cs23 = new TCanvas("cs23","cs23",10,10,855,500);
576  cs23->cd();
577  hnd2dplot2prongCCPi0contZ->Draw("colz");
578  hnd2dplot2prongCCPi0contZ->GetYaxis()->SetTitle("Containment in z-dir [cm]");
579  hnd2dplot2prongCCPi0contZ->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
580  hnd2dplot2prongCCPi0contZ->SetTitleOffset(0.9,"Y");
581  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
582  NDSimulation();
583  cs23->SaveAs(outDir2+CurrentFile[23].name.c_str()+CurrentDate[0].date.c_str()+".png");
584 /*********************************************************************************************************************/
585 /*********************************************************************************************************************/
586 ///NC Pi0 Interactions////
587 
588  /// Plot 1 vertex x-dir ///
589  TCanvas *cs24 = new TCanvas("cs24","cs24",10,10,855,500);
590  cs24->cd();
591  hnd2dplot2prongNCPi0x->Draw("colz");
592  hnd2dplot2prongNCPi0x->GetYaxis()->SetTitle("Vertex in x-dir [cm]");
593  hnd2dplot2prongNCPi0x->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
594  hnd2dplot2prongNCPi0x->SetTitleOffset(0.9,"Y");
595  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
596  NDSimulation();
597  cs24->SaveAs(outDir2+CurrentFile[24].name.c_str()+CurrentDate[0].date.c_str()+".png");
598 
599  /// Plot 2 vertex y-dir ///
600  TCanvas *cs25 = new TCanvas("cs25","cs25",10,10,855,500);
601  cs25->cd();
602  hnd2dplot2prongNCPi0y->Draw("colz");
603  hnd2dplot2prongNCPi0y->GetYaxis()->SetTitle("Vertex in y-dir [cm]");
604  hnd2dplot2prongNCPi0y->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
605  hnd2dplot2prongNCPi0y->SetTitleOffset(0.9,"Y");
606  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
607  NDSimulation();
608  cs25->SaveAs(outDir2+CurrentFile[25].name.c_str()+CurrentDate[0].date.c_str()+".png");
609 
610  /// Plot 3 vertex z-dir ///
611  TCanvas *cs26 = new TCanvas("cs26","cs26",10,10,855,500);
612  cs26->cd();
613  hnd2dplot2prongNCPi0z->Draw("colz");
614  hnd2dplot2prongNCPi0z->GetYaxis()->SetTitle("Vertex in z-dir [cm]");
615  hnd2dplot2prongNCPi0z->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
616  hnd2dplot2prongNCPi0z->SetTitleOffset(0.9,"Y");
617  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
618  NDSimulation();
619  cs26->SaveAs(outDir2+CurrentFile[26].name.c_str()+CurrentDate[0].date.c_str()+".png");
620 
621  /// Plot 4 Containment x-dir ///
622  TCanvas *cs27 = new TCanvas("cs27","cs27",10,10,855,500);
623  cs27->cd();
624  hnd2dplot2prongNCPi0contX->Draw("colz");
625  hnd2dplot2prongNCPi0contX->GetYaxis()->SetTitle("Containment in x-dir [cm]");
626  hnd2dplot2prongNCPi0contX->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
627  hnd2dplot2prongNCPi0contX->SetTitleOffset(0.9,"Y");
628  // hnd2dplot2prNng->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
629  NDSimulation();
630  cs27->SaveAs(outDir2+CurrentFile[27].name.c_str()+CurrentDate[0].date.c_str()+".png");
631 
632  /// Plot 5 Containment y-dir ///
633  TCanvas *cs28 = new TCanvas("cs28","cs28",10,10,855,500);
634  cs28->cd();
635  hnd2dplot2prongNCPi0contY->Draw("colz");
636  hnd2dplot2prongNCPi0contY->GetYaxis()->SetTitle("Containment in y-dir [cm]");
637  hnd2dplot2prongNCPi0contY->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
638  hnd2dplot2prongNCPi0contY->SetTitleOffset(0.9,"Y");
639  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
640  NDSimulation();
641  cs28->SaveAs(outDir2+CurrentFile[28].name.c_str()+CurrentDate[0].date.c_str()+".png");
642 
643  /// Plot 6 Containment z-dir ///
644  TCanvas *cs29 = new TCanvas("cs29","cs29",10,10,855,500);
645  cs29->cd();
646  hnd2dplot2prongNCPi0contZ->Draw("colz");
647  hnd2dplot2prongNCPi0contZ->GetYaxis()->SetTitle("Containment in z-dir [cm]");
648  hnd2dplot2prongNCPi0contZ->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
649  hnd2dplot2prongNCPi0contZ->SetTitleOffset(0.9,"Y");
650  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
651  NDSimulation();
652  cs29->SaveAs(outDir2+CurrentFile[29].name.c_str()+CurrentDate[0].date.c_str()+".png");
653 /*********************************************************************************************************************/
654 /*********************************************************************************************************************/
655 ///NC Pi0 low energy Interactions////
656 
657  /// Plot 1 vertex x-dir ///
658  TCanvas *cs30 = new TCanvas("cs30","cs30",10,10,855,500);
659  cs30->cd();
660  hnd2dplot2prongNCPi0Lx->Draw("colz");
661  hnd2dplot2prongNCPi0Lx->GetYaxis()->SetTitle("Vertex in x-dir [cm]");
662  hnd2dplot2prongNCPi0Lx->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
663  hnd2dplot2prongNCPi0Lx->SetTitleOffset(0.9,"Y");
664  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
665  NDSimulation();
666  cs30->SaveAs(outDir2+CurrentFile[30].name.c_str()+CurrentDate[0].date.c_str()+".png");
667 
668  /// Plot 2 vertex y-dir ///
669  TCanvas *cs31 = new TCanvas("cs31","cs31",10,10,855,500);
670  cs31->cd();
671  hnd2dplot2prongNCPi0Ly->Draw("colz");
672  hnd2dplot2prongNCPi0Ly->GetYaxis()->SetTitle("Vertex in y-dir [cm]");
673  hnd2dplot2prongNCPi0Ly->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
674  hnd2dplot2prongNCPi0Ly->SetTitleOffset(0.9,"Y");
675  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
676  NDSimulation();
677  cs31->SaveAs(outDir2+CurrentFile[31].name.c_str()+CurrentDate[0].date.c_str()+".png");
678 
679  /// Plot 3 vertex z-dir ///
680  TCanvas *cs32 = new TCanvas("cs32","cs32",10,10,855,500);
681  cs32->cd();
682  hnd2dplot2prongNCPi0Lz->Draw("colz");
683  hnd2dplot2prongNCPi0Lz->GetYaxis()->SetTitle("Vertex in z-dir [cm]");
684  hnd2dplot2prongNCPi0Lz->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
685  hnd2dplot2prongNCPi0Lz->SetTitleOffset(0.9,"Y");
686  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
687  NDSimulation();
688  cs32->SaveAs(outDir2+CurrentFile[32].name.c_str()+CurrentDate[0].date.c_str()+".png");
689 
690  /// Plot 4 Containment x-dir ///
691  TCanvas *cs33 = new TCanvas("cs33","cs33",10,10,855,500);
692  cs33->cd();
693  hnd2dplot2prongNCPi0LcontX->Draw("colz");
694  hnd2dplot2prongNCPi0LcontX->GetYaxis()->SetTitle("Containment in x-dir [cm]");
695  hnd2dplot2prongNCPi0LcontX->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
696  hnd2dplot2prongNCPi0LcontX->SetTitleOffset(0.9,"Y");
697  // hnd2dplot2prNng->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
698  NDSimulation();
699  cs33->SaveAs(outDir2+CurrentFile[33].name.c_str()+CurrentDate[0].date.c_str()+".png");
700 
701  /// Plot 5 Containment y-dir ///
702  TCanvas *cs34 = new TCanvas("cs34","cs34",10,10,855,500);
703  cs34->cd();
704  hnd2dplot2prongNCPi0LcontY->Draw("colz");
705  hnd2dplot2prongNCPi0LcontY->GetYaxis()->SetTitle("Containment in y-dir [cm]");
706  hnd2dplot2prongNCPi0LcontY->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
707  hnd2dplot2prongNCPi0LcontY->SetTitleOffset(0.9,"Y");
708  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
709  NDSimulation();
710  cs34->SaveAs(outDir2+CurrentFile[34].name.c_str()+CurrentDate[0].date.c_str()+".png");
711 
712  /// Plot 6 Containment z-dir ///
713  TCanvas *cs35 = new TCanvas("cs35","cs35",10,10,855,500);
714  cs35->cd();
715  hnd2dplot2prongNCPi0LcontZ->Draw("colz");
716  hnd2dplot2prongNCPi0LcontZ->GetYaxis()->SetTitle("Containment in z-dir [cm]");
717  hnd2dplot2prongNCPi0LcontZ->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
718  hnd2dplot2prongNCPi0LcontZ->SetTitleOffset(0.9,"Y");
719  // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
720  NDSimulation();
721  cs35->SaveAs(outDir2+CurrentFile[35].name.c_str()+CurrentDate[0].date.c_str()+".png");
722 /*********************************************************************************************************************/
723  f1->Close();
724 
725 }
const XML_Char * name
Definition: expat.h:151
::xsd::cxx::tree::date< char, simple_type > date
Definition: Database.h:186
const DateUp CurrentDate[Date]
const int Date
std::unique_ptr< T > LoadFrom(TDirectory *dir, const std::string &label)
Definition: LoadFromFile.h:17
std::string name
const FileUp CurrentFile[File]
void NDSimulation()
Float_t f1
std::string date
TString outDir2
const int File
enum BeamMode string