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().

{
  TLatex* prelim = new TLatex(.9, .95, "NOvA Simulation");
  prelim->SetTextColor(kGray+1);
  prelim->SetNDC();
  prelim->SetTextSize(2/30.);
  prelim->SetTextAlign(32);
   TGaxis::SetMaxDigits(3);
  prelim->Draw();
}

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.

{



/*********************************************************************************************************************/

///Read the Spectrum Root file /////

  TFile *f1= new TFile ("/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/2dvtxcontstudyplots2dvtxcontsutyplotsdifferentinteractions.root","READ");
  f1->cd(); 

 /////All interactions///////////

  ///Vertex
 auto snd2dplot2prongx=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongx"));
 auto snd2dplot2prongy=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongy"));
 auto snd2dplot2prongz=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongz"));

 //Containment
 auto snd2dplot2prongcontX=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongcontX"));
 auto snd2dplot2prongcontY=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongcontY"));
 auto snd2dplot2prongcontZ=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongcontZ"));

/*********************************************************************************************************************/
 /////CC Non pi0 Interactions///////////

  ///Vertex
 auto snd2dplot2prongCCNonPi0x=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongCCNonPi0x"));
 auto snd2dplot2prongCCNonPi0y=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongCCNonPi0y"));
 auto snd2dplot2prongCCNonPi0z=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongCCNonPi0z"));

 //Containment
 auto snd2dplot2prongCCNonPi0contX=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongCCNonPi0contX"));
 auto snd2dplot2prongCCNonPi0contY=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongCCNonPi0contY"));
 auto snd2dplot2prongCCNonPi0contZ=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongCCNonPi0contZ"));
/*********************************************************************************************************************/
 /////NC Non pi0 Interactions///////////

  ///Vertex
 auto snd2dplot2prongNCNonPi0x=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCNonPi0x"));
 auto snd2dplot2prongNCNonPi0y=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCNonPi0y"));
 auto snd2dplot2prongNCNonPi0z=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCNonPi0z"));

 //Containment
 auto snd2dplot2prongNCNonPi0contX=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCNonPi0contX"));
 auto snd2dplot2prongNCNonPi0contY=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCNonPi0contY"));
 auto snd2dplot2prongNCNonPi0contZ=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCNonPi0contZ"));

/*********************************************************************************************************************/
 /////CC pi0 Interactions///////////

   ///Vertex
 auto snd2dplot2prongCCPi0x=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongCCPi0x"));
 auto snd2dplot2prongCCPi0y=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongCCPi0y"));
 auto snd2dplot2prongCCPi0z=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongCCPi0z"));

 //Containment
 auto snd2dplot2prongCCPi0contX=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongCCPi0contX"));
 auto snd2dplot2prongCCPi0contY=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongCCPi0contY"));
 auto snd2dplot2prongCCPi0contZ=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongCCPi0contZ"));
/*********************************************************************************************************************/
 /////NC pi0 Interactions///////////

   ///Vertex
 auto snd2dplot2prongNCPi0x=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCPi0x"));
 auto snd2dplot2prongNCPi0y=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCPi0y"));
 auto snd2dplot2prongNCPi0z=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCPi0z"));

 //Containment
 auto snd2dplot2prongNCPi0contX=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCPi0contX"));
 auto snd2dplot2prongNCPi0contY=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCPi0contY"));
 auto snd2dplot2prongNCPi0contZ=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCPi0contZ"));
/*********************************************************************************************************************/
 /////NC pi0 low energy Interactions///////////

   ///Vertex
 auto snd2dplot2prongNCPi0Lx=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCPi0Lx"));
 auto snd2dplot2prongNCPi0Ly=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCPi0Ly"));
 auto snd2dplot2prongNCPi0Lz=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCPi0Lz"));

 //Containment
 auto snd2dplot2prongNCPi0LcontX=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCPi0LcontX"));
 auto snd2dplot2prongNCPi0LcontY=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCPi0LcontY"));
 auto snd2dplot2prongNCPi0LcontZ=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prongNCPi0LcontZ"));




/*********************************************************************************************************************/
  ///How to Scale Histogram////
 
  // double myPOT=sndMC.POT();
  const double myPOT = 8e20;
/*********************************************************************************************************************/

  //All interactions

  //Vertex
  auto hnd2dplot2prongx = snd2dplot2prongx.ToTH2(myPOT);
  auto hnd2dplot2prongy = snd2dplot2prongy.ToTH2(myPOT);
  auto hnd2dplot2prongz = snd2dplot2prongz.ToTH2(myPOT);

  //Containment
  auto hnd2dplot2prongcontX = snd2dplot2prongcontX.ToTH2(myPOT);
  auto hnd2dplot2prongcontY = snd2dplot2prongcontY.ToTH2(myPOT);
  auto hnd2dplot2prongcontZ = snd2dplot2prongcontZ.ToTH2(myPOT);
 /*********************************************************************************************************************/
  //CC Non pi0 Interactions///////

  //Vertex
  auto hnd2dplot2prongCCNonPi0x = snd2dplot2prongCCNonPi0x.ToTH2(myPOT);
  auto hnd2dplot2prongCCNonPi0y = snd2dplot2prongCCNonPi0y.ToTH2(myPOT);
  auto hnd2dplot2prongCCNonPi0z = snd2dplot2prongCCNonPi0z.ToTH2(myPOT);

  //Containment
  auto hnd2dplot2prongCCNonPi0contX = snd2dplot2prongCCNonPi0contX.ToTH2(myPOT);
  auto hnd2dplot2prongCCNonPi0contY = snd2dplot2prongCCNonPi0contY.ToTH2(myPOT);
  auto hnd2dplot2prongCCNonPi0contZ = snd2dplot2prongCCNonPi0contZ.ToTH2(myPOT);

 /*********************************************************************************************************************/
  //NC Non pi0 Interactions///////

  //Vertex
  auto hnd2dplot2prongNCNonPi0x = snd2dplot2prongNCNonPi0x.ToTH2(myPOT);
  auto hnd2dplot2prongNCNonPi0y = snd2dplot2prongNCNonPi0y.ToTH2(myPOT);
  auto hnd2dplot2prongNCNonPi0z = snd2dplot2prongNCNonPi0z.ToTH2(myPOT);

  //Containment
  auto hnd2dplot2prongNCNonPi0contX = snd2dplot2prongNCNonPi0contX.ToTH2(myPOT);
  auto hnd2dplot2prongNCNonPi0contY = snd2dplot2prongNCNonPi0contY.ToTH2(myPOT);
  auto hnd2dplot2prongNCNonPi0contZ = snd2dplot2prongNCNonPi0contZ.ToTH2(myPOT);


/*********************************************************************************************************************/

  //CC Pi0 Interactions//////////////////

  //Vertex
  auto hnd2dplot2prongCCPi0x = snd2dplot2prongCCPi0x.ToTH2(myPOT);
  auto hnd2dplot2prongCCPi0y = snd2dplot2prongCCPi0y.ToTH2(myPOT);
  auto hnd2dplot2prongCCPi0z = snd2dplot2prongCCPi0z.ToTH2(myPOT);

  //Containment
  auto hnd2dplot2prongCCPi0contX = snd2dplot2prongCCPi0contX.ToTH2(myPOT);
  auto hnd2dplot2prongCCPi0contY = snd2dplot2prongCCPi0contY.ToTH2(myPOT);
  auto hnd2dplot2prongCCPi0contZ = snd2dplot2prongCCPi0contZ.ToTH2(myPOT);

/*********************************************************************************************************************/

  //NC Pi0 Interactions//////////////////

  //Vertex
  auto hnd2dplot2prongNCPi0x = snd2dplot2prongNCPi0x.ToTH2(myPOT);
  auto hnd2dplot2prongNCPi0y = snd2dplot2prongNCPi0y.ToTH2(myPOT);
  auto hnd2dplot2prongNCPi0z = snd2dplot2prongNCPi0z.ToTH2(myPOT);

  //Containment
  auto hnd2dplot2prongNCPi0contX = snd2dplot2prongNCPi0contX.ToTH2(myPOT);
  auto hnd2dplot2prongNCPi0contY = snd2dplot2prongNCPi0contY.ToTH2(myPOT);
  auto hnd2dplot2prongNCPi0contZ = snd2dplot2prongNCPi0contZ.ToTH2(myPOT);

/*********************************************************************************************************************/

  //NC Pi0 low energy Interactions//////////////////

  //Vertex
  auto hnd2dplot2prongNCPi0Lx = snd2dplot2prongNCPi0Lx.ToTH2(myPOT);
  auto hnd2dplot2prongNCPi0Ly = snd2dplot2prongNCPi0Ly.ToTH2(myPOT);
  auto hnd2dplot2prongNCPi0Lz = snd2dplot2prongNCPi0Lz.ToTH2(myPOT);

  //Containment
  auto hnd2dplot2prongNCPi0LcontX = snd2dplot2prongNCPi0LcontX.ToTH2(myPOT);
  auto hnd2dplot2prongNCPi0LcontY = snd2dplot2prongNCPi0LcontY.ToTH2(myPOT);
  auto hnd2dplot2prongNCPi0LcontZ = snd2dplot2prongNCPi0LcontZ.ToTH2(myPOT);
/*********************************************************************************************************************/ 
//****************************Updated Date and string path for saving files **********************************************************///
  TString outDir2 = "/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/twodimvtxcontplots";
  struct FileUp{
    std::string name;
  };
  struct DateUp{
    std::string date;
  };
 
  const int File =37;
  const FileUp CurrentFile[File] = {
    {"/Vertex_X_Direction_AllInteractions"},
    {"/Vertex_Y_Direction_AllInteractions"},
    {"/Vertex_Z_Direction_AllInteractions"},
    {"/Containment_X_Direction_AllInteractions"},
    {"/Containment_Y_Direction_AllInteractions"},
    {"/Containment_Z_Direction_AllInteractions"},

    {"/Vertex_X_Direction_CCNonPi0"},
    {"/Vertex_Y_Direction_CCNonPi0"},
    {"/Vertex_Z_Direction_CCNonPi0"},
    {"/Containment_X_Direction_CCNonPi0"},
    {"/Containment_Y_Direction_CCNonPi0"},
    {"/Containment_Z_Direction_CCNonPi0"},

    {"/Vertex_X_Direction_NCNonPi0"},
    {"/Vertex_Y_Direction_NCNonPi0"},
    {"/Vertex_Z_Direction_NCNonPi0"},
    {"/Containment_X_Direction_NCNonPi0"},
    {"/Containment_Y_Direction_NCNonPi0"},
    {"/Containment_Z_Direction_NCNonPi0"},
    
    {"/Vertex_X_Direction_CCPi0"},
    {"/Vertex_Y_Direction_CCPi0"},
    {"/Vertex_Z_Direction_CCPi0"},
    {"/Containment_X_Direction_CCPi0"},
    {"/Containment_Y_Direction_CCPi0"},
    {"/Containment_Z_Direction_CCPi0"},

    {"/Vertex_X_Direction_NCPi0HighEnergy"},
    {"/Vertex_Y_Direction_NCPi0HighEnergy"},
    {"/Vertex_Z_Direction_NCPi0HighEnergy"},
    {"/Containment_X_Direction_NCPi0HighEnergy"},
    {"/Containment_Y_Direction_NCPi0HighEnergy"},
    {"/Containment_Z_Direction_NCPi0HighEnergy"},    
  
    {"/Vertex_X_Direction_NCPi0_LowEnergy"},
    {"/Vertex_Y_Direction_NCPi0_LowEnergy"},
    {"/Vertex_Z_Direction_NCPi0_LowEnergy"},
    {"/Containment_X_Direction_NCPi0_LowEnergy"},
    {"/Containment_Y_Direction_NCPi0_LowEnergy"},
    {"/Containment_Z_Direction_NCPi0_Lowenergy"},
  };

   const int Date = 1;
   const DateUp CurrentDate[Date] = {
     {"_03-30-18"},
  };


///All Interactions////

   /// Plot 1 vertex x-dir ///
   TCanvas *cs1 = new TCanvas("cs1","cs1",10,10,855,500);
   cs1->cd();
   hnd2dplot2prongx->Draw("colz");
   hnd2dplot2prongx->GetYaxis()->SetTitle("Vertex in x-dir [cm]");
   hnd2dplot2prongx->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongx->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs1->SaveAs(outDir2+CurrentFile[0].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /// Plot 2 vertex y-dir ///
   TCanvas *cs2 = new TCanvas("cs2","cs2",10,10,855,500);
   cs2->cd();
   hnd2dplot2prongy->Draw("colz");
   hnd2dplot2prongy->GetYaxis()->SetTitle("Vertex in y-dir [cm]");
   hnd2dplot2prongy->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongy->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs2->SaveAs(outDir2+CurrentFile[1].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /// Plot 3 vertex z-dir ///
   TCanvas *cs3 = new TCanvas("cs3","cs3",10,10,855,500);
   cs3->cd();
   hnd2dplot2prongz->Draw("colz");
   hnd2dplot2prongz->GetYaxis()->SetTitle("Vertex in z-dir [cm]");
   hnd2dplot2prongz->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongz->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs3->SaveAs(outDir2+CurrentFile[2].name.c_str()+CurrentDate[0].date.c_str()+".png");

     /// Plot 4 Containment x-dir ///
   TCanvas *cs4 = new TCanvas("cs4","cs4",10,10,855,500);
   cs4->cd();
   hnd2dplot2prongcontX->Draw("colz");
   hnd2dplot2prongcontX->GetYaxis()->SetTitle("Containment in x-dir [cm]");
   hnd2dplot2prongcontX->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongcontX->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs4->SaveAs(outDir2+CurrentFile[3].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /// Plot 5 Containment y-dir ///
   TCanvas *cs5 = new TCanvas("cs5","cs5",10,10,855,500);
   cs5->cd();
   hnd2dplot2prongcontY->Draw("colz");
   hnd2dplot2prongcontY->GetYaxis()->SetTitle("Containment in y-dir [cm]");
   hnd2dplot2prongcontY->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongcontY->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs5->SaveAs(outDir2+CurrentFile[4].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /// Plot 6 Containment z-dir ///
   TCanvas *cs6 = new TCanvas("cs6","cs6",10,10,855,500);
   cs6->cd();
   hnd2dplot2prongcontZ->Draw("colz");
   hnd2dplot2prongcontZ->GetYaxis()->SetTitle("Containment in z-dir [cm]");
   hnd2dplot2prongcontZ->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongcontZ->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs6->SaveAs(outDir2+CurrentFile[5].name.c_str()+CurrentDate[0].date.c_str()+".png");
/*********************************************************************************************************************/
/*********************************************************************************************************************/ 


///CC Non pi0 Interactions////

   /// Plot 1 vertex x-dir ///
   TCanvas *cs7 = new TCanvas("cs7","cs7",10,10,855,500);
   cs7->cd();
   hnd2dplot2prongCCNonPi0x->Draw("colz");
   hnd2dplot2prongCCNonPi0x->GetYaxis()->SetTitle("Vertex in x-dir [cm]");
   hnd2dplot2prongCCNonPi0x->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongCCNonPi0x->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs7->SaveAs(outDir2+CurrentFile[6].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /// Plot 2 vertex y-dir ///
   TCanvas *cs8 = new TCanvas("cs8","cs8",10,10,855,500);
   cs8->cd();
   hnd2dplot2prongCCNonPi0y->Draw("colz");
   hnd2dplot2prongCCNonPi0y->GetYaxis()->SetTitle("Vertex in y-dir [cm]");
   hnd2dplot2prongCCNonPi0y->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongCCNonPi0y->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs8->SaveAs(outDir2+CurrentFile[7].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /// Plot 3 vertex z-dir ///
   TCanvas *cs55 = new TCanvas("cs55","cs55",10,10,855,500);
   cs55->cd();
   hnd2dplot2prongCCNonPi0z->Draw("colz");
   hnd2dplot2prongCCNonPi0z->GetYaxis()->SetTitle("Vertex in z-dir [cm]");
   hnd2dplot2prongCCNonPi0z->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongCCNonPi0z->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs55->SaveAs(outDir2+CurrentFile[8].name.c_str()+CurrentDate[0].date.c_str()+".png");

     /// Plot 4 Containment x-dir ///
   TCanvas *cs9 = new TCanvas("cs9","cs9",10,10,855,500);
   cs9->cd();
   hnd2dplot2prongCCNonPi0contX->Draw("colz");
   hnd2dplot2prongCCNonPi0contX->GetYaxis()->SetTitle("Containment in x-dir [cm]");
   hnd2dplot2prongCCNonPi0contX->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongCCNonPi0contX->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs9->SaveAs(outDir2+CurrentFile[9].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /// Plot 5 Containment y-dir ///
   TCanvas *cs10 = new TCanvas("cs10","cs10",10,10,855,500);
   cs10->cd();
   hnd2dplot2prongCCNonPi0contY->Draw("colz");
   hnd2dplot2prongCCNonPi0contY->GetYaxis()->SetTitle("Containment in y-dir [cm]");
   hnd2dplot2prongCCNonPi0contY->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongCCNonPi0contY->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs10->SaveAs(outDir2+CurrentFile[10].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /// Plot 6 Containment z-dir ///
   TCanvas *cs11 = new TCanvas("cs11","cs11",10,10,855,500);
   cs11->cd();
   hnd2dplot2prongCCNonPi0contZ->Draw("colz");
   hnd2dplot2prongCCNonPi0contZ->GetYaxis()->SetTitle("Containment in z-dir [cm]");
   hnd2dplot2prongCCNonPi0contZ->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongCCNonPi0contZ->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs11->SaveAs(outDir2+CurrentFile[11].name.c_str()+CurrentDate[0].date.c_str()+".png");
/*********************************************************************************************************************/
/*********************************************************************************************************************/ 


///NC Non pi0 Interactions////

   /// Plot 1 vertex x-dir ///
   TCanvas *cs12 = new TCanvas("cs12","cs12",10,10,855,500);
   cs12->cd();
   hnd2dplot2prongNCNonPi0x->Draw("colz");
   hnd2dplot2prongNCNonPi0x->GetYaxis()->SetTitle("Vertex in x-dir [cm]");
   hnd2dplot2prongNCNonPi0x->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongNCNonPi0x->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs12->SaveAs(outDir2+CurrentFile[12].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /// Plot 2 vertex y-dir ///
   TCanvas *cs13 = new TCanvas("cs13","cs13",10,10,855,500);
   cs13->cd();
   hnd2dplot2prongNCNonPi0y->Draw("colz");
   hnd2dplot2prongNCNonPi0y->GetYaxis()->SetTitle("Vertex in y-dir [cm]");
   hnd2dplot2prongNCNonPi0y->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongNCNonPi0y->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs13->SaveAs(outDir2+CurrentFile[13].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /// Plot 3 vertex z-dir ///
   TCanvas *cs14 = new TCanvas("cs14","cs14",10,10,855,500);
   cs14->cd();
   hnd2dplot2prongNCNonPi0z->Draw("colz");
   hnd2dplot2prongNCNonPi0z->GetYaxis()->SetTitle("Vertex in z-dir [cm]");
   hnd2dplot2prongNCNonPi0z->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongNCNonPi0z->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs14->SaveAs(outDir2+CurrentFile[14].name.c_str()+CurrentDate[0].date.c_str()+".png");

     /// Plot 4 Containment x-dir ///
   TCanvas *cs15 = new TCanvas("cs15","cs15",10,10,855,500);
   cs15->cd();
   hnd2dplot2prongNCNonPi0contX->Draw("colz");
   hnd2dplot2prongNCNonPi0contX->GetYaxis()->SetTitle("Containment in x-dir [cm]");
   hnd2dplot2prongNCNonPi0contX->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongNCNonPi0contX->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs15->SaveAs(outDir2+CurrentFile[15].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /// Plot 5 Containment y-dir ///
   TCanvas *cs16 = new TCanvas("cs16","cs16",10,10,855,500);
   cs16->cd();
   hnd2dplot2prongNCNonPi0contY->Draw("colz");
   hnd2dplot2prongNCNonPi0contY->GetYaxis()->SetTitle("Containment in y-dir [cm]");
   hnd2dplot2prongNCNonPi0contY->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongNCNonPi0contY->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs16->SaveAs(outDir2+CurrentFile[16].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /// Plot 6 Containment z-dir ///
   TCanvas *cs17 = new TCanvas("cs17","cs17",10,10,855,500);
   cs17->cd();
   hnd2dplot2prongNCNonPi0contZ->Draw("colz");
   hnd2dplot2prongNCNonPi0contZ->GetYaxis()->SetTitle("Containment in z-dir [cm]");
   hnd2dplot2prongNCNonPi0contZ->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongNCNonPi0contZ->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs17->SaveAs(outDir2+CurrentFile[17].name.c_str()+CurrentDate[0].date.c_str()+".png");
/*********************************************************************************************************************/
///CC Pi0 Interactions////

   /// Plot 1 vertex x-dir ///
   TCanvas *cs18 = new TCanvas("cs18","cs18",10,10,855,500);
   cs18->cd();
   hnd2dplot2prongCCPi0x->Draw("colz");
   hnd2dplot2prongCCPi0x->GetYaxis()->SetTitle("Vertex in x-dir [cm]");
   hnd2dplot2prongCCPi0x->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongCCPi0x->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs18->SaveAs(outDir2+CurrentFile[18].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /// Plot 2 vertex y-dir ///
   TCanvas *cs19 = new TCanvas("cs19","cs19",10,10,855,500);
   cs19->cd();
   hnd2dplot2prongCCPi0y->Draw("colz");
   hnd2dplot2prongCCPi0y->GetYaxis()->SetTitle("Vertex in y-dir [cm]");
   hnd2dplot2prongCCPi0y->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongCCPi0y->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs19->SaveAs(outDir2+CurrentFile[19].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /// Plot 3 vertex z-dir ///
   TCanvas *cs20 = new TCanvas("cs20","cs20",10,10,855,500);
   cs20->cd();
   hnd2dplot2prongCCPi0z->Draw("colz");
   hnd2dplot2prongCCPi0z->GetYaxis()->SetTitle("Vertex in z-dir [cm]");
   hnd2dplot2prongCCPi0z->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongCCPi0z->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs20->SaveAs(outDir2+CurrentFile[20].name.c_str()+CurrentDate[0].date.c_str()+".png");

     /// Plot 4 Containment x-dir ///
   TCanvas *cs21 = new TCanvas("cs21","cs21",10,10,855,500);
   cs21->cd();
   hnd2dplot2prongCCPi0contX->Draw("colz");
   hnd2dplot2prongCCPi0contX->GetYaxis()->SetTitle("Containment in x-dir [cm]");
   hnd2dplot2prongCCPi0contX->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongCCPi0contX->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs21->SaveAs(outDir2+CurrentFile[21].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /// Plot 5 Containment y-dir ///
   TCanvas *cs22 = new TCanvas("cs22","cs22",10,10,855,500);
   cs22->cd();
   hnd2dplot2prongCCPi0contY->Draw("colz");
   hnd2dplot2prongCCPi0contY->GetYaxis()->SetTitle("Containment in y-dir [cm]");
   hnd2dplot2prongCCPi0contY->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongCCPi0contY->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs22->SaveAs(outDir2+CurrentFile[22].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /// Plot 6 Containment z-dir ///
   TCanvas *cs23 = new TCanvas("cs23","cs23",10,10,855,500);
   cs23->cd();
   hnd2dplot2prongCCPi0contZ->Draw("colz");
   hnd2dplot2prongCCPi0contZ->GetYaxis()->SetTitle("Containment in z-dir [cm]");
   hnd2dplot2prongCCPi0contZ->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongCCPi0contZ->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs23->SaveAs(outDir2+CurrentFile[23].name.c_str()+CurrentDate[0].date.c_str()+".png");
/*********************************************************************************************************************/
/*********************************************************************************************************************/
///NC Pi0 Interactions////

   /// Plot 1 vertex x-dir ///
   TCanvas *cs24 = new TCanvas("cs24","cs24",10,10,855,500);
   cs24->cd();
   hnd2dplot2prongNCPi0x->Draw("colz");
   hnd2dplot2prongNCPi0x->GetYaxis()->SetTitle("Vertex in x-dir [cm]");
   hnd2dplot2prongNCPi0x->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongNCPi0x->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs24->SaveAs(outDir2+CurrentFile[24].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /// Plot 2 vertex y-dir ///
   TCanvas *cs25 = new TCanvas("cs25","cs25",10,10,855,500);
   cs25->cd();
   hnd2dplot2prongNCPi0y->Draw("colz");
   hnd2dplot2prongNCPi0y->GetYaxis()->SetTitle("Vertex in y-dir [cm]");
   hnd2dplot2prongNCPi0y->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongNCPi0y->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs25->SaveAs(outDir2+CurrentFile[25].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /// Plot 3 vertex z-dir ///
   TCanvas *cs26 = new TCanvas("cs26","cs26",10,10,855,500);
   cs26->cd();
   hnd2dplot2prongNCPi0z->Draw("colz");
   hnd2dplot2prongNCPi0z->GetYaxis()->SetTitle("Vertex in z-dir [cm]");
   hnd2dplot2prongNCPi0z->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongNCPi0z->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs26->SaveAs(outDir2+CurrentFile[26].name.c_str()+CurrentDate[0].date.c_str()+".png");

     /// Plot 4 Containment x-dir ///
   TCanvas *cs27 = new TCanvas("cs27","cs27",10,10,855,500);
   cs27->cd();
   hnd2dplot2prongNCPi0contX->Draw("colz");
   hnd2dplot2prongNCPi0contX->GetYaxis()->SetTitle("Containment in x-dir [cm]");
   hnd2dplot2prongNCPi0contX->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongNCPi0contX->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prNng->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs27->SaveAs(outDir2+CurrentFile[27].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /// Plot 5 Containment y-dir ///
   TCanvas *cs28 = new TCanvas("cs28","cs28",10,10,855,500);
   cs28->cd();
   hnd2dplot2prongNCPi0contY->Draw("colz");
   hnd2dplot2prongNCPi0contY->GetYaxis()->SetTitle("Containment in y-dir [cm]");
   hnd2dplot2prongNCPi0contY->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongNCPi0contY->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs28->SaveAs(outDir2+CurrentFile[28].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /// Plot 6 Containment z-dir ///
   TCanvas *cs29 = new TCanvas("cs29","cs29",10,10,855,500);
   cs29->cd();
   hnd2dplot2prongNCPi0contZ->Draw("colz");
   hnd2dplot2prongNCPi0contZ->GetYaxis()->SetTitle("Containment in z-dir [cm]");
   hnd2dplot2prongNCPi0contZ->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongNCPi0contZ->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs29->SaveAs(outDir2+CurrentFile[29].name.c_str()+CurrentDate[0].date.c_str()+".png");
/*********************************************************************************************************************/
/*********************************************************************************************************************/
///NC Pi0 low energy Interactions////

   /// Plot 1 vertex x-dir ///
   TCanvas *cs30 = new TCanvas("cs30","cs30",10,10,855,500);
   cs30->cd();
   hnd2dplot2prongNCPi0Lx->Draw("colz");
   hnd2dplot2prongNCPi0Lx->GetYaxis()->SetTitle("Vertex in x-dir [cm]");
   hnd2dplot2prongNCPi0Lx->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongNCPi0Lx->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs30->SaveAs(outDir2+CurrentFile[30].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /// Plot 2 vertex y-dir ///
   TCanvas *cs31 = new TCanvas("cs31","cs31",10,10,855,500);
   cs31->cd();
   hnd2dplot2prongNCPi0Ly->Draw("colz");
   hnd2dplot2prongNCPi0Ly->GetYaxis()->SetTitle("Vertex in y-dir [cm]");
   hnd2dplot2prongNCPi0Ly->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongNCPi0Ly->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs31->SaveAs(outDir2+CurrentFile[31].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /// Plot 3 vertex z-dir ///
   TCanvas *cs32 = new TCanvas("cs32","cs32",10,10,855,500);
   cs32->cd();
   hnd2dplot2prongNCPi0Lz->Draw("colz");
   hnd2dplot2prongNCPi0Lz->GetYaxis()->SetTitle("Vertex in z-dir [cm]");
   hnd2dplot2prongNCPi0Lz->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongNCPi0Lz->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs32->SaveAs(outDir2+CurrentFile[32].name.c_str()+CurrentDate[0].date.c_str()+".png");

     /// Plot 4 Containment x-dir ///
   TCanvas *cs33 = new TCanvas("cs33","cs33",10,10,855,500);
   cs33->cd();
   hnd2dplot2prongNCPi0LcontX->Draw("colz");
   hnd2dplot2prongNCPi0LcontX->GetYaxis()->SetTitle("Containment in x-dir [cm]");
   hnd2dplot2prongNCPi0LcontX->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongNCPi0LcontX->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prNng->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs33->SaveAs(outDir2+CurrentFile[33].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /// Plot 5 Containment y-dir ///
   TCanvas *cs34 = new TCanvas("cs34","cs34",10,10,855,500);
   cs34->cd();
   hnd2dplot2prongNCPi0LcontY->Draw("colz");
   hnd2dplot2prongNCPi0LcontY->GetYaxis()->SetTitle("Containment in y-dir [cm]");
   hnd2dplot2prongNCPi0LcontY->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongNCPi0LcontY->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs34->SaveAs(outDir2+CurrentFile[34].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /// Plot 6 Containment z-dir ///
   TCanvas *cs35 = new TCanvas("cs35","cs35",10,10,855,500);
   cs35->cd();
   hnd2dplot2prongNCPi0LcontZ->Draw("colz");
   hnd2dplot2prongNCPi0LcontZ->GetYaxis()->SetTitle("Containment in z-dir [cm]");
   hnd2dplot2prongNCPi0LcontZ->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prongNCPi0LcontZ->SetTitleOffset(0.9,"Y");
   // hnd2dplot2prong->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 Cut");
  NDSimulation();
  cs35->SaveAs(outDir2+CurrentFile[35].name.c_str()+CurrentDate[0].date.c_str()+".png");
/*********************************************************************************************************************/
 f1->Close();

}