plotHist_SpectrumCVNID.C

Go to the documentation of this file.

/* Gets ND Data/MC Numu varaibles spectra and converts them into POT normalized histograms.
 */
#ifdef __CINT__
void plotHist_SpectrumCVNID()
{
  std::cout<<" Sorry, you must run in compiled mode" <<std::endl;
}
#else
#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 "TGaxis.h"
#include "TAttLine.h"
#include "TAttMarker.h"

#include <cmath>
#include <iostream>
#include <vector>
#include <list>
#include <sstream>
#include <string>
#include <sstream>
#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"



using namespace ana;
void NDSimulation()
{
  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 NDSimulation1()
{
  TLatex* prelim = new TLatex(.9, .95, "NOvA Preliminary");
  prelim->SetTextColor(kBlue);
  prelim->SetNDC();
  prelim->SetTextSize(2/30.);
  prelim->SetTextAlign(32);
   TGaxis::SetMaxDigits(3);
  prelim->Draw();
}

void plotHist_SpectrumCVNID()
{



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

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

  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");
  f1->cd(); 
 
  //////////INVARIANT MASSSS//////////////////////

  //load w/ "NoCUTS" (only 2 prong) pi0 mass//////
  auto sndNoCuts=*Spectrum::LoadFrom(f1->GetDirectory("ndNoCuts"));
  auto sndNoCutCCNonPi0=*Spectrum::LoadFrom(f1->GetDirectory("ndNoCutCCNonPi0"));
  auto sndNoCutNCNonPi0=*Spectrum::LoadFrom(f1->GetDirectory("ndNoCutNCNonPi0"));
  auto sndNoCutCCPi0=*Spectrum::LoadFrom(f1->GetDirectory("ndNoCutCCPi0"));
  auto sndNoCutNCPi0=*Spectrum::LoadFrom(f1->GetDirectory("ndNoCutNCPi0"));
  auto sndNoCutNCPi0L=*Spectrum::LoadFrom(f1->GetDirectory("ndNoCutNCPi0L"));
  auto sndNumiMasssdataNoCuts =*Spectrum::LoadFrom(f1->GetDirectory("ndNumiMassdataNoCuts"));



  ///load reco mass 2 prong CVN cuts spectrums///// (stacked) /////
  auto sndMC=*Spectrum::LoadFrom(f1->GetDirectory("ndMC"));
  auto sndMC3P=*Spectrum::LoadFrom(f1->GetDirectory("ndMC3P"));
  auto sndNCNonPi0=*Spectrum::LoadFrom(f1->GetDirectory("ndNCNonPi0"));
  auto sndNumiMassdata2CVN=*Spectrum::LoadFrom(f1->GetDirectory("ndNumiMassdata2CVN"));
  auto sndNCPi0=*Spectrum::LoadFrom(f1->GetDirectory("ndNCPi0"));
  auto sndCCPi0=*Spectrum::LoadFrom(f1->GetDirectory("ndCCPi0"));
  auto sndNCPi0L=*Spectrum::LoadFrom(f1->GetDirectory("ndNCPi0L"));
  auto sndCCNonPi0=*Spectrum::LoadFrom(f1->GetDirectory("ndCCNonPi0"));
 
 ///load reco mass 2 prong CVN cuts spectrums/////      (Unstacked)/////
  auto sndMC1=*Spectrum::LoadFrom(f1->GetDirectory("ndMC"));
  auto sndMC3P1=*Spectrum::LoadFrom(f1->GetDirectory("ndMC3P"));
  auto sndNCNonPi01=*Spectrum::LoadFrom(f1->GetDirectory("ndNCNonPi0"));
  //auto sndNCMCSig=*Spectrum::LoadFrom(f1->GetDirectory("ndNCMCSig"));
  auto sndNCPi01=*Spectrum::LoadFrom(f1->GetDirectory("ndNCPi0"));
  auto sndCCPi01=*Spectrum::LoadFrom(f1->GetDirectory("ndCCPi0"));
  auto sndNCPi0L1=*Spectrum::LoadFrom(f1->GetDirectory("ndNCPi0L"));
  auto sndCCNonPi01=*Spectrum::LoadFrom(f1->GetDirectory("ndCCNonPi0"));
  auto sndNumiMassdata2CVN1=*Spectrum::LoadFrom(f1->GetDirectory("ndNumiMassdata2CVN"));



  ///load reco mass 1 prong CVN cut spectrums/////     (stacked)/////
  auto sndMC1Cut=*Spectrum::LoadFrom(f1->GetDirectory("ndMC1Cut"));
  auto sndMC3P1Cut=*Spectrum::LoadFrom(f1->GetDirectory("ndMC3P1Cut"));
  auto sndNCNonPi01Cut=*Spectrum::LoadFrom(f1->GetDirectory("ndNCNonPi01Cut"));
  //auto sndNCMCSig=*Spectrum::LoadFrom(f1->GetDirectory("ndNCMCSig"));
  auto sndNCPi01Cut=*Spectrum::LoadFrom(f1->GetDirectory("ndNCPi01Cut"));
  auto sndCCPi01Cut=*Spectrum::LoadFrom(f1->GetDirectory("ndCCPi01Cut"));
  auto sndNCPi0L1Cut=*Spectrum::LoadFrom(f1->GetDirectory("ndNCPi0L1Cut"));
  auto sndCCNonPi01Cut=*Spectrum::LoadFrom(f1->GetDirectory("ndCCNonPi01Cut"));
  auto sndNumiMassdataCVN=*Spectrum::LoadFrom(f1->GetDirectory("ndNumiMassdataCVN"));


  ///load reco mass 1 prong CVN cut spectrums/////          (unstacked) /////
  auto sndMC1Cut2=*Spectrum::LoadFrom(f1->GetDirectory("ndMC1Cut"));
  auto sndMC3P1Cut2=*Spectrum::LoadFrom(f1->GetDirectory("ndMC3P1Cut"));
  auto sndNCNonPi01Cut2=*Spectrum::LoadFrom(f1->GetDirectory("ndNCNonPi01Cut"));
  //auto sndNCMCSig=*Spectrum::LoadFrom(f1->GetDirectory("ndNCMCSig"));
  auto sndNCPi01Cut2=*Spectrum::LoadFrom(f1->GetDirectory("ndNCPi01Cut"));
  auto sndCCPi01Cut2=*Spectrum::LoadFrom(f1->GetDirectory("ndCCPi01Cut"));
  auto sndNCPi0L1Cut2=*Spectrum::LoadFrom(f1->GetDirectory("ndNCPi0L1Cut"));
  auto sndCCNonPi01Cut2=*Spectrum::LoadFrom(f1->GetDirectory("ndCCNonPi01Cut"));
  auto sndNumiMassdataCVN2=*Spectrum::LoadFrom(f1->GetDirectory("ndNumiMassdataCVN"));
  


////////CVN STUDY////////////////
  
  //load /w/Prel + remID + Prong 1 CVN cuts & variable (2prongs photonID spectrums////
 
  auto sndCVNMC3P=*Spectrum::LoadFrom(f1->GetDirectory("ndCVNMC3P"));
  auto sndCVNNCNonPi0=*Spectrum::LoadFrom(f1->GetDirectory("ndCVNNCNonPi0"));
  auto sndCVNNCPi0=*Spectrum::LoadFrom(f1->GetDirectory("ndCVNNCPi0"));
  auto sndCVNCCPi0=*Spectrum::LoadFrom(f1->GetDirectory("ndCVNCCPi0"));
  auto sndCVNNCPi0L=*Spectrum::LoadFrom(f1->GetDirectory("ndCVNNCPi0L"));
  auto sndCVNCCNonPi0=*Spectrum::LoadFrom(f1->GetDirectory("ndCVNCCNonPi0"));
  auto sndNumiCVN2PData=*Spectrum::LoadFrom(f1->GetDirectory("ndNumiCVN2PData"));
  //load 1 prongCVN variable & Presel + remid cut spectrums////

  auto sndCVNMC3PR=*Spectrum::LoadFrom(f1->GetDirectory("ndCVNMC3PR"));
  auto sndCVNNCNonPi0R=*Spectrum::LoadFrom(f1->GetDirectory("ndCVNNCNonPi0R"));
  auto sndCVNNCPi0R=*Spectrum::LoadFrom(f1->GetDirectory("ndCVNNCPi0R"));
  auto sndCVNCCPi0R=*Spectrum::LoadFrom(f1->GetDirectory("ndCVNCCPi0R"));
  auto sndCVNNCPi0LR=*Spectrum::LoadFrom(f1->GetDirectory("ndCVNNCPi0LR"));
  auto sndCVNCCNonPi0R=*Spectrum::LoadFrom(f1->GetDirectory("ndCVNCCNonPi0R"));
  auto sndNumiCVN2PRData=*Spectrum::LoadFrom(f1->GetDirectory("ndNumiCVN2PRData"));
  
  //////2 prong CVN variable & Presel + remid cut ///////////
  auto snd2CVNMC3PR=*Spectrum::LoadFrom(f1->GetDirectory("nd2CVNMC3PR"));
  auto snd2CVNNCNonPi0R=*Spectrum::LoadFrom(f1->GetDirectory("nd2CVNNCNonPi0R"));
  auto snd2CVNNCPi0R=*Spectrum::LoadFrom(f1->GetDirectory("nd2CVNNCPi0R"));
  auto snd2CVNCCPi0R=*Spectrum::LoadFrom(f1->GetDirectory("nd2CVNCCPi0R"));
  auto snd2CVNNCPi0LR=*Spectrum::LoadFrom(f1->GetDirectory("nd2CVNNCPi0LR"));
  auto snd2CVNCCNonPi0R=*Spectrum::LoadFrom(f1->GetDirectory("nd2CVNCCNonPi0R"));
  auto sndNumi2CVN2PRData=*Spectrum::LoadFrom(f1->GetDirectory("ndNumi2CVN2PRData"));
  
  /// load RemID variable & w/Presel + 2 prong cuts  spectrums /////

  auto sndMC3PRem=*Spectrum::LoadFrom(f1->GetDirectory("ndMC3PRem"));
  auto sndNCNonPi0Rem=*Spectrum::LoadFrom(f1->GetDirectory("ndNCNonPi0Rem"));
  auto sndNCPi0Rem=*Spectrum::LoadFrom(f1->GetDirectory("ndNCPi0Rem"));
  auto sndCCPi0Rem=*Spectrum::LoadFrom(f1->GetDirectory("ndCCPi0Rem"));
  auto sndNCPi0LRem=*Spectrum::LoadFrom(f1->GetDirectory("ndNCPi0LRem"));
  auto sndCCNonPi0Rem=*Spectrum::LoadFrom(f1->GetDirectory("ndCCNonPi0Rem"));
  auto sndNumiRemIDData=*Spectrum::LoadFrom(f1->GetDirectory("ndNumiRemIDData"));

  ///////TRUE ENERGY STUDY//////////////////

  ////True Highest Pi0 Energy w/Prel + 2 Prong//////

  auto sndMC3PPrel=*Spectrum::LoadFrom(f1->GetDirectory("ndMC3PPrel"));
  auto sndNCNonPi0Prel=*Spectrum::LoadFrom(f1->GetDirectory("ndNCNonPi0Prel"));
  auto sndCCPi0Prel=*Spectrum::LoadFrom(f1->GetDirectory("ndCCPi0Prel"));
  auto sndNCPi0Prel=*Spectrum::LoadFrom(f1->GetDirectory("ndNCPi0Prel"));
  auto sndNCPi0LPrel=*Spectrum::LoadFrom(f1->GetDirectory("ndNCPi0LPrel"));
  auto sndCCNonPi0Prel=*Spectrum::LoadFrom(f1->GetDirectory("ndCCNonPi0Prel"));

  ////True Highest Pi0 Energy w/Prel + RemId < 0.36 + 1  Prong CVN Gamma ID > 0.7 + 2 prong  //////

  auto sndMC3PPrelCVN =*Spectrum::LoadFrom(f1->GetDirectory("ndMC3PPrelCVN"));
  auto sndNCNonPi0PrelCVN =*Spectrum::LoadFrom(f1->GetDirectory("ndNCNonPi0PrelCVN"));
  auto sndCCPi0PrelCVN =*Spectrum::LoadFrom(f1->GetDirectory("ndCCPi0PrelCVN"));
  auto sndNCPi0PrelCVN =*Spectrum::LoadFrom(f1->GetDirectory("ndNCPi0PrelCVN"));
  auto sndNCPi0LPrelCVN =*Spectrum::LoadFrom(f1->GetDirectory("ndNCPi0LPrelCVN"));
  auto sndCCNonPi0PrelCVN =*Spectrum::LoadFrom(f1->GetDirectory("ndCCNonPi0PrelCVN"));
  
  ////True Highest Pi0 Energy w/Prel + RemId < 0.36 + 1 & 2 Prong CVN Gamma ID > 0.7 + 2 prong  //////

  auto sndMC3PPrelCVN2Cuts =*Spectrum::LoadFrom(f1->GetDirectory("ndMC3PPrelCVN2Cuts"));
  auto sndNCNonPi0PrelCVN2Cuts =*Spectrum::LoadFrom(f1->GetDirectory("ndNCNonPi0PrelCVN2Cuts"));
  auto sndCCPi0PrelCVN2Cuts =*Spectrum::LoadFrom(f1->GetDirectory("ndCCPi0PrelCVN2Cuts"));
  auto sndNCPi0PrelCVN2Cuts =*Spectrum::LoadFrom(f1->GetDirectory("ndNCPi0PrelCVN2Cuts"));
  auto sndNCPi0LPrelCVN2Cuts =*Spectrum::LoadFrom(f1->GetDirectory("ndNCPi0LPrelCVN2Cuts"));
  auto sndCCNonPi0PrelCVN2Cuts =*Spectrum::LoadFrom(f1->GetDirectory("ndCCNonPi0PrelCVN2Cuts"));

  ////True Highest Pi0 Energy w/Prel + RemID + 2 Prong//////

  auto sndMC3PPrelR =*Spectrum::LoadFrom(f1->GetDirectory("ndMC3PPrelR"));
  auto sndNCNonPi0PrelR =*Spectrum::LoadFrom(f1->GetDirectory("ndNCNonPi0PrelR"));
  auto sndCCPi0PrelR =*Spectrum::LoadFrom(f1->GetDirectory("ndCCPi0PrelR"));
  auto sndNCPi0PrelR =*Spectrum::LoadFrom(f1->GetDirectory("ndNCPi0PrelR"));
  auto sndNCPi0LPrelR =*Spectrum::LoadFrom(f1->GetDirectory("ndNCPi0LPrelR"));
  auto sndCCNonPi0PrelR =*Spectrum::LoadFrom(f1->GetDirectory("ndCCNonPi0PrelR"));




 // ////2D Plots 2 prong & 3 prong///////
  auto snd2dplot2prong=*Spectrum::LoadFrom(f1->GetDirectory("nd2dplot2prong"));
/*********************************************************************************************************************/

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

  /////////////////////INVARIANT MASS////////////////////////

  //load w/ "NoCUTS" (only 2 prong) pi0 mass//////
  auto hndNoCuts = sndNoCuts.ToTH1(myPOT);
  auto hndNoCutCCNonPi0 = sndNoCutCCNonPi0.ToTH1(myPOT);
  auto hndNoCutNCNonPi0 = sndNoCutNCNonPi0.ToTH1(myPOT);
  auto hndNoCutCCPi0 = sndNoCutCCPi0.ToTH1(myPOT);
  auto hndNoCutNCPi0 = sndNoCutNCPi0.ToTH1(myPOT);
  auto hndNoCutNCPi0L = sndNoCutNCPi0L.ToTH1(myPOT);
  auto hndNumiMasssdataNoCuts = sndNumiMasssdataNoCuts.ToTH1(myPOT);

  ///load reco mass spectrums 2 prong CVN cuts (stacked) /////
  auto hndMC = sndMC.ToTH1(myPOT);
  auto hndMC3P = sndMC3P.ToTH1(myPOT);
  auto hndNCNonPi0 = sndNCNonPi0.ToTH1(myPOT);
  auto hndCCPi0 = sndCCPi0.ToTH1(myPOT);
  auto hndNCPi0 = sndNCPi0.ToTH1(myPOT);
  auto hndNCPi0L = sndNCPi0L.ToTH1(myPOT);
  auto hndCCNonPi0 = sndCCNonPi0.ToTH1(myPOT);
  auto hndNumiMassdata2CVN = sndNumiMassdata2CVN.ToTH1(myPOT);

  ///load reco mass spectrums 2 prong CVN cuts (unstacked) /////
  auto hndMC1 = sndMC.ToTH1(myPOT);
  auto hndMC3P1 = sndMC3P.ToTH1(myPOT);
  auto hndNCNonPi01 = sndNCNonPi0.ToTH1(myPOT);
  auto hndCCPi01 = sndCCPi0.ToTH1(myPOT);
  auto hndNCPi01 = sndNCPi0.ToTH1(myPOT);
  auto hndNCPi0L1 = sndNCPi0L.ToTH1(myPOT);
  auto hndCCNonPi01 = sndCCNonPi0.ToTH1(myPOT);
  auto hndNumiMassdata2CVN1 = sndNumiMassdata2CVN1.ToTH1(myPOT);

  ///load reco mass spectrums 1 prong CVN cut (Stacked) /////
  auto hndMC1Cut = sndMC1Cut.ToTH1(myPOT);
  auto hndMC3P1Cut = sndMC3P1Cut.ToTH1(myPOT);
  auto hndNCNonPi01Cut = sndNCNonPi01Cut.ToTH1(myPOT);
  auto hndCCPi01Cut = sndCCPi01Cut.ToTH1(myPOT);
  auto hndNCPi01Cut = sndNCPi01Cut.ToTH1(myPOT);
  auto hndNCPi0L1Cut = sndNCPi0L1Cut.ToTH1(myPOT);
  auto hndCCNonPi01Cut = sndCCNonPi01Cut.ToTH1(myPOT);
  auto hndNumiMassdataCVN = sndNumiMassdataCVN.ToTH1(myPOT);

  ///load reco mass spectrums 1 prong CVN cut (Unstacked)/////
  auto hndMC1Cut2 = sndMC1Cut2.ToTH1(myPOT);
  auto hndMC3P1Cut2 = sndMC3P1Cut2.ToTH1(myPOT);
  auto hndNCNonPi01Cut2 = sndNCNonPi01Cut2.ToTH1(myPOT);
  auto hndCCPi01Cut2 = sndCCPi01Cut2.ToTH1(myPOT);
  auto hndNCPi01Cut2 = sndNCPi01Cut2.ToTH1(myPOT);
  auto hndNCPi0L1Cut2 = sndNCPi0L1Cut2.ToTH1(myPOT);
  auto hndCCNonPi01Cut2 = sndCCNonPi01Cut2.ToTH1(myPOT);
  auto hndNumiMassdataCVN2 = sndNumiMassdataCVN2.ToTH1(myPOT);
////////CVN STUDY////////////////

   //load /w/Prel + remID + Prong 1 CVN cuts & 2prong photon ID variable (2prongs photonID spectrums)////
  auto hndCVNMC3P = sndCVNMC3P.ToTH1(myPOT);
  auto hndCVNNCNonPi0 = sndCVNNCNonPi0.ToTH1(myPOT);
  auto hndCVNCCPi0 = sndCVNCCPi0.ToTH1(myPOT);
  auto hndCVNNCPi0 = sndCVNNCPi0.ToTH1(myPOT);
  auto hndCVNNCPi0L = sndCVNNCPi0L.ToTH1(myPOT);
  auto hndCVNCCNonPi0 = sndCVNCCNonPi0.ToTH1(myPOT);
  auto hndNumiCVN2PData = sndNumiCVN2PData.ToTH1(myPOT); 

  //load 1 prongCVN variable & Presel + remid cut spectrums////
  auto hndCVNMC3PR = sndCVNMC3PR.ToTH1(myPOT);
  auto hndCVNNCNonPi0R = sndCVNNCNonPi0R.ToTH1(myPOT);
  auto hndCVNCCPi0R = sndCVNCCPi0R.ToTH1(myPOT);
  auto hndCVNNCPi0R = sndCVNNCPi0R.ToTH1(myPOT);
  auto hndCVNNCPi0LR = sndCVNNCPi0LR.ToTH1(myPOT);
  auto hndCVNCCNonPi0R = sndCVNCCNonPi0R.ToTH1(myPOT);
  auto hndNumiCVN2PRData = sndNumiCVN2PRData.ToTH1(myPOT);
  
  //////2 prong CVN variable & Presel + remid cut ///////////
  auto hnd2CVNMC3PR = snd2CVNMC3PR.ToTH1(myPOT);
  auto hnd2CVNNCNonPi0R = snd2CVNNCNonPi0R.ToTH1(myPOT);
  auto hnd2CVNCCPi0R = snd2CVNCCPi0R.ToTH1(myPOT);
  auto hnd2CVNNCPi0R = snd2CVNNCPi0R.ToTH1(myPOT);
  auto hnd2CVNNCPi0LR = snd2CVNNCPi0LR.ToTH1(myPOT);
  auto hnd2CVNCCNonPi0R = snd2CVNCCNonPi0R.ToTH1(myPOT);
  auto hndNumi2CVN2PRData = sndNumi2CVN2PRData.ToTH1(myPOT);

  ////RemID//////

  auto hndMC3PRem = sndMC3PRem.ToTH1(myPOT);
  auto hndNCNonPi0Rem = sndNCNonPi0Rem.ToTH1(myPOT);
  auto hndCCPi0Rem = sndCCPi0Rem.ToTH1(myPOT);
  auto hndNCPi0Rem = sndNCPi0Rem.ToTH1(myPOT);
  auto hndNCPi0LRem = sndNCPi0LRem.ToTH1(myPOT);
  auto hndCCNonPi0Rem = sndCCNonPi0Rem.ToTH1(myPOT);
  auto hndNumiRemIDData = sndNumiRemIDData.ToTH1(myPOT);

  ////True Highest Pi0 Energy w/Prel + 2 Prong//////

  auto hndMC3PPrel = sndMC3PPrel.ToTH1(myPOT);
  auto hndNCNonPi0Prel = sndNCNonPi0Prel.ToTH1(myPOT);
  auto hndCCPi0Prel = sndCCPi0Prel.ToTH1(myPOT);
  auto hndNCPi0Prel = sndNCPi0Prel.ToTH1(myPOT);
  auto hndNCPi0LPrel = sndNCPi0LPrel.ToTH1(myPOT);
  auto hndCCNonPi0Prel = sndCCNonPi0Prel.ToTH1(myPOT);

  ////True Highest Pi0 Energy w/Prel + RemId < 0.36 + 1  Prong CVN Gamma ID > 0.7 + 2 prong  //////

  auto hndMC3PPrelCVN = sndMC3PPrelCVN.ToTH1(myPOT);
  auto hndNCNonPi0PrelCVN = sndNCNonPi0PrelCVN.ToTH1(myPOT);
  auto hndCCPi0PrelCVN = sndCCPi0PrelCVN.ToTH1(myPOT);
  auto hndNCPi0PrelCVN = sndNCPi0PrelCVN.ToTH1(myPOT);
  auto hndNCPi0LPrelCVN = sndNCPi0LPrelCVN.ToTH1(myPOT);
  auto hndCCNonPi0PrelCVN = sndCCNonPi0PrelCVN.ToTH1(myPOT);

  ////True Highest Pi0 Energy w/Prel + RemId < 0.36 + 1 & 2 Prong CVN Gamma ID > 0.7 + 2 prong  //////

  auto hndMC3PPrelCVN2Cuts = sndMC3PPrelCVN2Cuts.ToTH1(myPOT);
  auto hndNCNonPi0PrelCVN2Cuts = sndNCNonPi0PrelCVN2Cuts.ToTH1(myPOT);
  auto hndCCPi0PrelCVN2Cuts = sndCCPi0PrelCVN2Cuts.ToTH1(myPOT);
  auto hndNCPi0PrelCVN2Cuts = sndNCPi0PrelCVN2Cuts.ToTH1(myPOT);
  auto hndNCPi0LPrelCVN2Cuts = sndNCPi0LPrelCVN2Cuts.ToTH1(myPOT);
  auto hndCCNonPi0PrelCVN2Cuts = sndCCNonPi0PrelCVN2Cuts.ToTH1(myPOT);

  ////True Highest Pi0 Energy w/Prel + RemID + 2 Prong//////

  auto hndMC3PPrelR = sndMC3PPrelR.ToTH1(myPOT);
  auto hndNCNonPi0PrelR = sndNCNonPi0PrelR.ToTH1(myPOT);
  auto hndCCPi0PrelR = sndCCPi0PrelR.ToTH1(myPOT);
  auto hndNCPi0PrelR = sndNCPi0PrelR.ToTH1(myPOT);
  auto hndNCPi0LPrelR = sndNCPi0LPrelR.ToTH1(myPOT);
  auto hndCCNonPi0PrelR = sndCCNonPi0PrelR.ToTH1(myPOT);

  ////2D Plots 2 prong & 3 prong///////

  auto hnd2dplot2prong = snd2dplot2prong.ToTH2(myPOT);
//****************************Updated Date and string path for saving files **********************************************************///
  TString outDir2 = "/nova/app/users/acedeno/tag_releasesS18-06-14/NDAna/ncpi0_semi_inc_png_cvn/NewPNGCVAnalysis";
  //TString outDir2 = "/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/NcPi0Plots/March_PCVN_Analysis/";
  struct FileUp{
    std::string name;
  };
  struct DateUp{
    std::string date;
  };
 
  const int File =18;
  const FileUp CurrentFile[File] = {
    {"/plot1_2Pr_03GeV_Cuts_1CVN_07_2CVN_09_NOCUTS_Pi0Mass"},
    {"/plot2_2Pr_03GeV_Cuts_1CVN_07_2CVN_09_Selectionection"},
    {"/plot3_2Pr_03GeV_Cuts_1CVN_07_2CVN_09_StackedInvariantMass2Cuts"},
    {"/plot3b_2Pr_03GeV_Cuts_1CVN_07_StackedInvariantMass1Cut"},
    {"/plot4_2Pr_03GeV_Cuts_{RemID+Presel+1CVN_07}-CVN-Analysis-{Var-2P-PhotonID}"},
    {"/plot5_2Pr_03GeV_Cuts_{RemID+Selection}-CVN-Analysis-{Var-1P-PhotonID}"},
    {"/plot6_2Pr_03GeV_Cuts_{RemID+Selection}-CVN-Analysis-{Var-2P-PhotonID}"},
    {"/plot7_2Pr_03GeV_Cuts_{Selection}-CVN-Analysis-{Var-RemID}"},
    {"/plot8_2Pr_03GeV_Cuts_{Selection}-TrueEnergy-Analysis-{Var-TruePi0En}"},
    {"/plot9_2Pr_03GeV_Cuts_{RemID+Selection}-TrueEnergy-Analysis-{Var-TruePi0En}"},
    {"/plot10_2Pr_03GeV_Cuts_{RemID+Selection+1CVN_07}-TrueEnergy-Analysis-{Var-TruePi0En}"},
    {"/plot10b_2Pr_03GeV_Cuts_{RemID+Selection+1CVN_07+2CVN_09}-TrueEnergy-Analysis-{Var-TruePi0En}"},
    {"/plot11_2Pr_2dplot_03GeV_Cuts_1CVN_07_2CVN_09"},
    {"/plot13_2Pr_03GeV_Cuts_1CVN_07_2CVN_09_UnstackedInvariantMass2Cuts"},
    {"/plot13b_2Pr_03GeV_Cuts_1CVN_07_UnstackedInvariantMass1Cut"},
    {"/plot11b_2Pr_2dplot_LogVersion_03GeV_Cuts_1CVN_07_2CVN_09"},
    
  };

   const int Date = 1;
   const DateUp CurrentDate[Date] = {
     {"_07--18_newPng_GenieLike"},
  };
   //newPng1at07Png2at09_flat
//------------------------------------------------------------

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

  /// Plot  ////

  TCanvas *c1_ndNoCuts =new TCanvas("c1_ndNoCuts","ndNoCuts",800,500);
   c1_ndNoCuts->cd();
  
   hndNoCuts->SetLineColor(kRed+2);
   // hndNoCuts->SetLineWidth(3);
   // hndNoCuts->SetMarkerStyle(21);
   // hndNoCuts->SetMarkerColor(kBlack);

   hndNoCutCCPi0->SetLineColor(kMagenta+2);
   // hndNoCutCCPi0->SetLineWidth(1);
   //hndNoCutCCPi0->SetMarkerStyle(21);
   //hndNoCutCCPi0->SetMarkerColor(kMagenta);

   hndNoCutCCNonPi0->SetLineColor(kAzure+2);
   // hndNoCutCCNonPi0->SetLineWidth(3);
   // hndNoCutCCNonPi0->SetMarkerStyle(21);
   // hndNoCutCCNonPi0->SetMarkerColor(kAzure);

   hndNoCutNCPi0L->SetLineColor(kOrange+2);
   // hndNoCutNCPi0L->SetLineWidth(3);
   // hndNoCutNCPi0L->SetMarkerStyle(21);
   // hndNoCutNCPi0L->SetMarkerColor(kOrange);

   hndNoCutNCPi0->SetLineColor(kGreen+2);
   // hndNoCutNCPi0->SetLineWidth(3);
   // hndNoCutNCPi0->SetMarkerStyle(21);
   // hndNoCutNCPi0->SetMarkerColor(kGreen);

   hndNoCutNCNonPi0->SetLineColor(kCyan-6);
   // hndNoCutNCNonPi0->SetLineWidth(3);
   // hndNoCutNCNonPi0->SetMarkerStyle(21);
   // hndNoCutNCNonPi0->SetMarkerColor(kCyan);

   // hndNumiMasssdataNoCuts->SetLineColor(kBlack+2);
   hndNumiMasssdataNoCuts->SetMarkerStyle(2);
   hndNumiMasssdataNoCuts->SetMarkerSize(3);

   hndNoCuts->Draw("hist E1");
   hndNoCutCCPi0->Draw("hist E1 same");
   hndNoCutCCNonPi0->Draw("hist E1 same");
   hndNoCutNCPi0L->Draw("hist E1 same");
   hndNoCutNCPi0->Draw("hist E1 same");
   hndNoCutNCNonPi0->Draw("hist E1 same");
   //hndNumiMasssdataNoCuts->Draw("hist E1 same");
   //hndNoCuts->SetMaximum(300000);
   hndNoCuts->GetYaxis()->SetTitle("Events/8e20 POT");
   hndNoCuts->SetTitleOffset(0.95,"Y");
   // hndNoCuts->SetTitle("2 Prong Events No Cuts");

   TLegend * leg0 = new TLegend(0.45,0.26,0.55,0.5);
   leg0->SetTextFont(42);
   leg0->SetTextSize(0.05);
   //leg0->AddEntry(hndNumiMasssdataNoCuts,"Data","l");
   leg0->AddEntry(hndNoCuts,"No Cuts","l");
   leg0->AddEntry(hndNoCutCCPi0,"CCPi0","l");
   leg0->AddEntry(hndNoCutCCNonPi0,"CCNonPi0","l");
   leg0->AddEntry(hndNoCutNCPi0L,"NC w/Pi0 < 0.3 E","l");
   leg0->AddEntry(hndNoCutNCPi0,"NC w/Pi0 > 0.3 E","l");
   leg0->AddEntry(hndNoCutNCNonPi0,"NCNonPi0","l");

   leg0->Draw("same");
    NDSimulation();
    //c1_ndNoCuts->SaveAs("/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/NcPi0Plots/plot1_2Pr_0_3_CVN07_11-02-17.pdf");
    c1_ndNoCuts->SaveAs(outDir2+CurrentFile[0].name.c_str()+CurrentDate[0].date.c_str()+".png");

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


  /// Plot 1 ////

  TCanvas *c1_ndMC =new TCanvas("c1_ndMC","ndMC",800,500);
   c1_ndMC->cd();
  
   hndMC->SetLineColor(kBlue);
   // hndMC->SetLineWidth(3);
   // hndMC->SetMarkerStyle(21);
   // hndMC->SetMarkerColor(kBlue);
   hndMC3P->SetLineColor(kRed);
   // hndMC3P->SetLineWidth(3);
   //hndMC3P->SetFillColor(kRed);
   // hndMC3P->SetMarkerStyle(21);
   // hndMC3P->SetMarkerColor(kRed);
   hndMC->Draw("hist E1"); 
   hndMC3P->Draw("hist E1 same");
   hndMC->GetYaxis()->SetTitle("Events/8e20 POT");
   hndMC->SetTitleOffset(0.95,"Y");
   // hndMC->SetTitle("2 Prong Event Preselection");
   //gPad->Modified();
   TLegend * leg = new TLegend(0.55,0.7,0.65,0.8);
   leg->SetTextFont(42);
   leg->SetTextSize(0.05);
   leg->AddEntry(hndMC,"Fid.+ Cont. (PreSel)","l");
   leg->AddEntry(hndMC3P,"Selection","l");
   leg->Draw("same");
    NDSimulation();
    //c1_ndMC->SaveAs("/nova/app/users/acedeno/tag_releaseS17-05-31/NDAna/ncpi0/PngCVNID/NcPi0Plots/plot2_2Pr_0_3_CVN07__07-29-17.pdf");
    c1_ndMC->SaveAs(outDir2+CurrentFile[1].name.c_str()+CurrentDate[0].date.c_str()+".png");

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

/// Plot 2 ///
 
   //gStyle->SetTitleY(0.9);
   TCanvas *cs2 = new TCanvas("cs2","cs2",10,10,855,500);
   THStack *hs2 = new THStack("hs2"," ");
 //TCanvas *c2_ndMC =new TCanvas("c2_ndMC","ndMC",800,500);
   //TH1D *c2_ndNCMCSig = new TH1D("c2_ndNCMCSig","ndNCMCSig",100,-4,4);
   //hXndNumiMassdata->SetLineColor(kBlack);
   // hndNCMCSig->SetFillColor(kGray);
   // hndNCMCSig->SetLineWidth(3);
   // hndNCMCSig->SetMarkerStyle(21);
   // hndNCMCSig->SetMarkerColor(kGray);
   
   // hs2->Add(hndNCMCSig);
  
   //TH1D *c2_ndCCMCSig = new TH1D("c2_ndCCMCSig","ndCCMCSig",100,-4,4);
   hndNCPi0L->SetLineColor(kBlack);
   hndNCPi0L->SetFillColor(kOrange);
   hndNCPi0L->SetLineWidth(1);
   hndNCPi0L->SetMarkerStyle(21);
   //hndNCPi0L->SetMarkerColor(kOrange);
   hs2->Add(hndNCPi0L);

   //TH1D *c2_hndCCPi0 = new TH1D("c2_ndCCPi0","ndCCPi0",100,-4,4);
   hndCCPi0->SetLineColor(kBlack);
   hndCCPi0->SetFillColor(kMagenta);
   hndCCPi0->SetLineWidth(1);
   hndCCPi0->SetMarkerStyle(21);
   hndCCPi0->SetMarkerColor(kMagenta);

   hs2->Add(hndCCPi0);
   //TH1D *c2_hndCCNonPi0 = new TH1D("c2_hndCCNonPi0","hndCCNonPi0",100,-4,4);
   hndCCNonPi0->SetLineColor(kBlack);
   hndCCNonPi0->SetFillColor(kAzure);
   hndCCNonPi0->SetLineWidth(1);
   hndCCNonPi0->SetMarkerStyle(21);
   hndCCNonPi0->SetMarkerColor(kAzure);
   hs2->Add(hndCCNonPi0);
   
   //TH1D *c2_hndNCPi0 = new TH1D("c2_ndNCPi0","ndNCPi0",100,-4,4);
   hndNCNonPi0->SetLineColor(kBlack);
   hndNCNonPi0->SetFillColor(kCyan);
   hndNCNonPi0->SetLineWidth(1);
   hndNCNonPi0->SetMarkerStyle(21);
   hndNCNonPi0->SetMarkerColor(kCyan);
   hs2->Add(hndNCNonPi0);
   
   //TH1D *c2_hndNCPi0L = new TH1D("c2_ndNCPi0L","ndNCPi0L",100,-4,4);
   hndNCPi0->SetLineColor(kBlack);
   hndNCPi0->SetFillColor(kGreen);
   hndNCPi0->SetLineWidth(1);
   hndNCPi0->SetMarkerStyle(21);
   hndNCPi0->SetMarkerColor(kGreen);
   hs2->Add(hndNCPi0);

   //TCanvas *cs2 = new TCanvas("cs2","cs2",10,10,855,500);
   TText T; T.SetTextFont(42); T.SetTextAlign(21);
   hs2->Draw("hist"); T.DrawTextNDC(.5,.95,"");
   hs2->GetYaxis()->SetTitle("Events/ 8e20 POT");
   hs2->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   //xhs2->SetMaximum(1300);
   hs2->GetXaxis()->SetRangeUser(0.,0.4);
   gPad->Modified();
   cs2->cd();
   //cs2->SetTopMargin(0.18);
   hndMC3P->SetLineColor(kRed);
   hndMC3P->SetLineWidth(3);
   //hndNumiMassdata2CVN->SetMarkerStyle(2);
   hndNumiMassdata2CVN->SetMarkerSize(3);

   hndMC3P->Draw("same");
   //hndNumiMassdata2CVN->Draw("same");
   TLegend * leg2 = new TLegend(0.62,0.5,0.73,0.82);
   //leg2->SetTextAlign(40);
   leg2->SetTextFont(42);
   leg2->SetTextSize(0.05);
   leg2->SetFillStyle(0);
   //leg2->AddEntry(hndNumiMassdata2CVN,"Data","kFullDotSmall");
   leg2->AddEntry(hndMC3P,"Selection","kFullDotSmall");
   leg2->AddEntry(hndCCPi0,"CCPi0","kFullDotSmall");
   leg2->AddEntry(hndCCNonPi0,"CCNonPi0","kFullDotSmall");
   leg2->AddEntry(hndNCPi0L,"NC w/Pi0 < 0.3 E","kFullDotSmall");
   leg2->AddEntry(hndNCPi0,"NC w/Pi0 > 0.3 E","kFullDotSmall");
   leg2->AddEntry(hndNCNonPi0,"NCNonPi0","kFullDotSmall");
  
   
   leg2->Draw("same"); 
    NDSimulation();
    //cs2->SaveAs("/nova/app/users/acedeno/tag_releas9-24/NDAna/ncpi0_Systematics/NcPi0Plots/plot3_2Pr_0_3_CVN07__11-02-17.pdf");
    cs2->SaveAs(outDir2+CurrentFile[2].name.c_str()+CurrentDate[0].date.c_str()+".png");

/*********************************************************************************************************************/
/// Plot 2b ///
 
   //gStyle->SetTitleY(0.9);
   TCanvas *cs2b = new TCanvas("cs2b","cs2b",10,10,855,500);
   THStack *hs2b = new THStack("hs2b"," ");
 //TCanvas *c2_ndMC =new TCanvas("c2_ndMC","ndMC",800,500);
   //TH1D *c2_ndNCMCSig = new TH1D("c2_ndNCMCSig","ndNCMCSig",100,-4,4);
   // hndNCMCSig->SetLineColor(kBlack);
   // hndNCMCSig->SetFillColor(kGray);
   // hndNCMCSig->SetLineWidth(3);
   // hndNCMCSig->SetMarkerStyle(21);
   // hndNCMCSig->SetMarkerColor(kGray);
   
   // hs2->Add(hndNCMCSig);
  
   //TH1D *c2_ndCCMCSig = new TH1D("c2_ndCCMCSig","ndCCMCSig",100,-4,4);
   hndNCPi0L1Cut->SetLineColor(kBlack);
   hndNCPi0L1Cut->SetFillColor(kOrange);
   hndNCPi0L1Cut->SetLineWidth(1);
   hndNCPi0L1Cut->SetMarkerStyle(21);
   hndNCPi0L1Cut->SetMarkerColor(kOrange);
   hs2b->Add(hndNCPi0L1Cut);

   //TH1D *c2_hndCCPi0 = new TH1D("c2_ndCCPi0","ndCCPi0",100,-4,4);
   hndCCPi01Cut->SetLineColor(kBlack);
   hndCCPi01Cut->SetFillColor(kMagenta);
   hndCCPi01Cut->SetLineWidth(1);
   hndCCPi01Cut->SetMarkerStyle(21);
   hndCCPi01Cut->SetMarkerColor(kMagenta);

   hs2b->Add(hndCCPi01Cut);
   //TH1D *c2_hndCCNonPi0 = new TH1D("c2_hndCCNonPi0","hndCCNonPi0",100,-4,4);
   hndCCNonPi01Cut->SetLineColor(kBlack);
   hndCCNonPi01Cut->SetFillColor(kAzure);
   hndCCNonPi01Cut->SetLineWidth(1);
   hndCCNonPi01Cut->SetMarkerStyle(21);
   hndCCNonPi01Cut->SetMarkerColor(kAzure);
   hs2b->Add(hndCCNonPi01Cut);
   
   //TH1D *c2_hndNCPi0 = new TH1D("c2_ndNCPi0","ndNCPi0",100,-4,4);
   hndNCNonPi01Cut->SetLineColor(kBlack);
   hndNCNonPi01Cut->SetFillColor(kCyan);
   hndNCNonPi01Cut->SetLineWidth(1);
   hndNCNonPi01Cut->SetMarkerStyle(21);
   hndNCNonPi01Cut->SetMarkerColor(kCyan);
   hs2b->Add(hndNCNonPi01Cut);
   
   //TH1D *c2_hndNCPi0L = new TH1D("c2_ndNCPi0L","ndNCPi0L",100,-4,4);
   hndNCPi01Cut->SetLineColor(kBlack);
   hndNCPi01Cut->SetFillColor(kGreen);
   hndNCPi01Cut->SetLineWidth(1);
   hndNCPi01Cut->SetMarkerStyle(21);
   hndNCPi01Cut->SetMarkerColor(kGreen);
   hs2b->Add(hndNCPi01Cut);

   //TCanvas *cs2 = new TCanvas("cs2","cs2",10,10,855,500);
   TText T99; T99.SetTextFont(42); T99.SetTextAlign(21);
   hs2b->Draw("hist"); T99.DrawTextNDC(.5,.95,"");
   hs2b->GetYaxis()->SetTitle("Events/ 8e20 POT");
   hs2b->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hs2b->SetMaximum(2500);
   hs2b->GetXaxis()->SetRangeUser(0.,0.4);
   gPad->Modified();
   cs2b->cd();
   //cs2->SetTopMargin(0.18);
   hndMC3P1Cut->SetLineColor(kRed);
   hndMC3P1Cut->SetLineWidth(3);
   //hndNumiMassdataCVN->SetMarkerStyle(2);
   hndNumiMassdataCVN->SetMarkerSize(3);

   hndMC3P1Cut->Draw("same");
   //hndNumiMassdataCVN->Draw("same");
   TLegend * leg2b = new TLegend(0.62,0.5,0.73,0.82);
   //leg2->SetTextAlign(40);
   leg2b->SetTextFont(42);
   leg2b->SetTextSize(0.05);
   leg2b->SetFillStyle(0);
   //leg2b->AddEntry(hndNumiMassdataCVN,"Data","kFullStar");
   leg2b->AddEntry(hndMC3P1Cut,"Selection","kFullDotSmall");
   leg2b->AddEntry(hndCCPi01Cut,"CCPi0","kFullDotSmall");
   leg2b->AddEntry(hndCCNonPi01Cut,"CCNonPi0","kFullDotSmall");
   leg2b->AddEntry(hndNCPi0L1Cut,"NC w/Pi0 < 0.3 E","kFullDotSmall");
   leg2b->AddEntry(hndNCPi01Cut,"NC w/Pi0 > 0.3 E","kFullDotSmall");
   leg2b->AddEntry(hndNCNonPi01Cut,"NCNonPi0","kFullDotSmall");
   
   leg2b->Draw("same"); 
    NDSimulation();
    //cs2->SaveAs("/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/NcPi0Plots/plot3_2Pr_0_3_CVN07__11-02-17.pdf");
    cs2b->SaveAs(outDir2+CurrentFile[3].name.c_str()+CurrentDate[0].date.c_str()+".png");

/*********************************************************************************************************************/
/// Plot 3 ///
  //gStyle->SetTitleY(0.9);
   TCanvas *cs3 = new TCanvas("cs3","cs3",10,10,855,500);
   THStack *hs3 = new THStack("hs3","");
 //TCanvas *c2_ndMC =new TCanvas("c2_ndMC","ndMC",800,500);
   //TH1D *c2_ndNCMCSig = new TH1D("c2_ndNCMCSig","ndNCMCSig",100,-4,4);
   // hndNCMCSig->SetLineColor(kBlack);
   // hndNCMCSig->SetFillColor(kGray);
   // hndNCMCSig->SetLineWidth(3);
   // hndNCMCSig->SetMarkerStyle(21);
   // hndNCMCSig->SetMarkerColor(kGray);
   
   // hs2->Add(hndNCMCSig);
  
   //TH1D *c2_ndCCMCSig = new TH1D("c2_ndCCMCSig","ndCCMCSig",100,-4,4);
   hndCVNNCPi0L->SetLineColor(kBlack);
   hndCVNNCPi0L->SetFillColor(kOrange);
   hndCVNNCPi0L->SetLineWidth(1);;
   hndCVNNCPi0L->SetMarkerStyle(21);
   hndCVNNCPi0L->SetMarkerColor(kOrange);
   hs3->Add(hndCVNNCPi0L);

   //TH1D *c2_hndCCPi0 = new TH1D("c2_ndCCPi0","ndCCPi0",100,-4,4);
   hndCVNCCPi0->SetLineColor(kBlack);
   hndCVNCCPi0->SetFillColor(kMagenta);
   hndCVNCCPi0->SetLineWidth(1);;
   hndCVNCCPi0->SetMarkerStyle(21);
   hndCVNCCPi0->SetMarkerColor(kMagenta);

   hs3->Add(hndCVNCCPi0);
   //TH1D *c2_hndCCNonPi0 = new TH1D("c2_hndCCNonPi0","hndCCNonPi0",100,-4,4);
   hndCVNCCNonPi0->SetLineColor(kBlack);
   hndCVNCCNonPi0->SetFillColor(kAzure);
   hndCVNCCNonPi0->SetLineWidth(1);;
   hndCVNCCNonPi0->SetMarkerStyle(21);
   hndCVNCCNonPi0->SetMarkerColor(kAzure);
   hs3->Add(hndCVNCCNonPi0);
   
   //TH1D *c2_hndNCPi0 = new TH1D("c2_ndNCPi0","ndNCPi0",100,-4,4);
   hndCVNNCNonPi0->SetLineColor(kBlack);
   hndCVNNCNonPi0->SetFillColor(kCyan);
   hndCVNNCNonPi0->SetLineWidth(1);;
   hndCVNNCNonPi0->SetMarkerStyle(21);
   hndCVNNCNonPi0->SetMarkerColor(kCyan);
   hs3->Add(hndCVNNCNonPi0);
   
   //TH1D *c2_hndNCPi0L = new TH1D("c2_ndNCPi0L","ndNCPi0L",100,-4,4);
   hndCVNNCPi0->SetLineColor(kBlack);
   hndCVNNCPi0->SetFillColor(kGreen);
   hndCVNNCPi0->SetLineWidth(1);;
   hndCVNNCPi0->SetMarkerStyle(21);
   hndCVNNCPi0->SetMarkerColor(kGreen);
   hs3->Add(hndCVNNCPi0);

   //TCanvas *cs2 = new TCanvas("cs2","cs2",10,10,855,500);
   TText T2; T2.SetTextFont(42); T2.SetTextAlign(21);
   hs3->Draw("hist"); T2.DrawTextNDC(.5,.95,"");
   hs3->GetYaxis()->SetTitle("Events/ 8e20 POT");
   hs3->GetXaxis()->SetTitle("Prong 2 CVN Gamma ID");
   //hs3->SetMaximum(6000);
   hs3->GetXaxis()->SetRangeUser(0,1);
   gPad->Modified();
   cs3->cd();
   //cs2->SetTopMargin(0.18);
   hndCVNMC3P->SetLineColor(kRed);
   hndCVNMC3P->SetLineWidth(3);
   //hndNumiCVN2PData->SetMarkerStyle(2);
   hndNumiCVN2PData->SetMarkerSize(3);

   hndCVNMC3P->Draw("same");
   //hndNumiCVN2PData->Draw("same");
   TLegend * leg3 = new TLegend(0.55,0.5,0.65,0.82);
   //leg2->SetTextAlign(40);
   leg3->SetTextFont(42);
   leg3->SetTextSize(0.05);
   leg3->SetFillStyle(0);
   //leg3->AddEntry(hndNumiCVN2PData,"Data","kFullDotSmall");
   leg3->AddEntry(hndCVNMC3P,"Selection","kFullDotSmall");
   leg3->AddEntry(hndCVNCCPi0,"CCPi0","kFullDotSmall");
   leg3->AddEntry(hndCVNCCNonPi0,"CCNonPi0","kFullDotSmall");
   leg3->AddEntry(hndCVNNCPi0L,"NC w/Pi0 < 0.3 E","kFullDotSmall");
   leg3->AddEntry(hndCVNNCPi0,"NC w/Pi0 > 0.3 E","kFullDotSmall");
   leg3->AddEntry(hndCVNNCNonPi0,"NCNonPi0","kFullDotSmall");
   
   leg3->Draw("same"); 
    NDSimulation();
    //cs3->SaveAs("/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/NcPi0Plots/plot4_2Pr_0_3_CVN07__11-02-17.pdf");
    cs3->SaveAs(outDir2+CurrentFile[4].name.c_str()+CurrentDate[0].date.c_str()+".png");

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

/// Plot 4 ///
  //gStyle->SetTitleY(0.9);
   TCanvas *cs4 = new TCanvas("cs4","cs4",10,10,855,500);
   THStack *hs4 = new THStack("hs4","");
 //TCanvas *c2_ndMC =new TCanvas("c2_ndMC","ndMC",800,500);
   //TH1D *c2_ndNCMCSig = new TH1D("c2_ndNCMCSig","ndNCMCSig",100,-4,4);
   // hndNCMCSig->SetLineColor(kBlack);
   // hndNCMCSig->SetFillColor(kGray);
   // hndNCMCSig->SetLineWidth(3);
   // hndNCMCSig->SetMarkerStyle(21);
   // hndNCMCSig->SetMarkerColor(kGray);
   
   // hs2->Add(hndNCMCSig);
  
   //TH1D *c2_ndCCMCSig = new TH1D("c2_ndCCMCSig","ndCCMCSig",100,-4,4);
   hndCVNNCPi0LR->SetLineColor(kBlack);
   hndCVNNCPi0LR->SetFillColor(kOrange);
   hndCVNNCPi0LR->SetLineWidth(1);;
   hndCVNNCPi0LR->SetMarkerStyle(21);
   hndCVNNCPi0LR->SetMarkerColor(kOrange);
   hs4->Add(hndCVNNCPi0LR);

   //TH1D *c2_hndCCPi0 = new TH1D("c2_ndCCPi0","ndCCPi0",100,-4,4);
   hndCVNCCPi0R->SetLineColor(kBlack);
   hndCVNCCPi0R->SetFillColor(kMagenta);
   hndCVNCCPi0R->SetLineWidth(1);;
   hndCVNCCPi0R->SetMarkerStyle(21);
   hndCVNCCPi0R->SetMarkerColor(kMagenta);

   hs4->Add(hndCVNCCPi0R);
   //TH1D *c2_hndCCNonPi0 = new TH1D("c2_hndCCNonPi0","hndCCNonPi0",100,-4,4);
   hndCVNCCNonPi0R->SetLineColor(kBlack);
   hndCVNCCNonPi0R->SetFillColor(kAzure);
   hndCVNCCNonPi0R->SetLineWidth(1);;
   hndCVNCCNonPi0R->SetMarkerStyle(21);
   hndCVNCCNonPi0R->SetMarkerColor(kAzure);
   hs4->Add(hndCVNCCNonPi0R);
   
   //TH1D *c2_hndNCPi0 = new TH1D("c2_ndNCPi0","ndNCPi0",100,-4,4);
   hndCVNNCNonPi0R->SetLineColor(kBlack);
   hndCVNNCNonPi0R->SetFillColor(kCyan);
   hndCVNNCNonPi0R->SetLineWidth(1);;
   hndCVNNCNonPi0R->SetMarkerStyle(21);
   hndCVNNCNonPi0R->SetMarkerColor(kCyan);
   hs4->Add(hndCVNNCNonPi0R);
   
   //TH1D *c2_hndNCPi0L = new TH1D("c2_ndNCPi0L","ndNCPi0L",100,-4,4);
   hndCVNNCPi0R->SetLineColor(kBlack);
   hndCVNNCPi0R->SetFillColor(kGreen);
   hndCVNNCPi0R->SetLineWidth(1);;
   hndCVNNCPi0R->SetMarkerStyle(21);
   hndCVNNCPi0R->SetMarkerColor(kGreen);
   hs4->Add(hndCVNNCPi0R);
   
   //TCanvas *cs2 = new TCanvas("cs2","cs2",10,10,855,500);
   TText T4; T4.SetTextFont(42); T4.SetTextAlign(21);
   hs4->Draw("hist"); T4.DrawTextNDC(.5,.95,"");
   hs4->GetYaxis()->SetTitle("Events/ 8e20 POT");
   hs4->GetXaxis()->SetTitle("Prong CVN Gamma ID");
   hs4->SetMaximum(10000);
   hs4->GetXaxis()->SetRangeUser(0,1);
   gPad->Modified();
   cs4->cd();
   //cs2->SetTopMargin(0.18);
   hndCVNMC3PR->SetLineColor(kRed);
   hndCVNMC3PR->SetLineWidth(3);
   //hndNumiCVN2PRData->SetMarkerStyle(2);
   hndNumiCVN2PRData->SetMarkerSize(3);


   hndCVNMC3PR->Draw("same");
   //hndNumiCVN2PRData->Draw("same");
   TLegend * leg4 = new TLegend(0.55,0.5,0.65,0.82);
   //leg2->SetTextAlign(40);
   leg4->SetTextFont(42);
   leg4->SetTextSize(0.05);
   leg4->SetFillStyle(0);
   //leg4->AddEntry(hndNumiCVN2PRData,"Data","kFullDotSmall");
   leg4->AddEntry(hndCVNMC3PR,"Selection","kFullDotSmall");
   leg4->AddEntry(hndCVNCCPi0R,"CCPi0","kFullDotSmall");
   leg4->AddEntry(hndCVNCCNonPi0R,"CCNonPi0","kFullDotSmall");
   leg4->AddEntry(hndCVNNCPi0LR,"NC w/Pi0 < 0.3 E","kFullDotSmall");
   leg4->AddEntry(hndCVNNCPi0R,"NC w/Pi0 > 0.3 E","kFullDotSmall");
   leg4->AddEntry(hndCVNNCNonPi0R,"NCNonPi0","kFullDotSmall");
   
   leg4->Draw("same"); 
    NDSimulation();
    //cs4->SaveAs("/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/NcPi0Plots/plot5_2Pr_0_3_CVN07__11-02-17.pdf");
    cs4->SaveAs(outDir2+CurrentFile[5].name.c_str()+CurrentDate[0].date.c_str()+".png");


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

/// Plot 5 ///
  //gStyle->SetTitleY(0.9);
   TCanvas *cs5 = new TCanvas("cs5","cs5",10,10,855,500);
   THStack *hs5 = new THStack("hs5","");
 //TCanvas *c2_ndMC =new TCanvas("c2_ndMC","ndMC",800,500);
   //TH1D *c2_ndNCMCSig = new TH1D("c2_ndNCMCSig","ndNCMCSig",100,-4,4);
   // hndNCMCSig->SetLineColor(kBlack);
   // hndNCMCSig->SetFillColor(kGray);
   // hndNCMCSig->SetLineWidth(3);
   // hndNCMCSig->SetMarkerStyle(21);
   // hndNCMCSig->SetMarkerColor(kGray);
   
   // hs2->Add(hndNCMCSig);
  
   //TH1D *c2_ndCCMCSig = new TH1D("c2_ndCCMCSig","ndCCMCSig",100,-4,4);
   hnd2CVNNCPi0LR->SetLineColor(kBlack);
   hnd2CVNNCPi0LR->SetFillColor(kOrange);
   hnd2CVNNCPi0LR->SetLineWidth(1);;
   hnd2CVNNCPi0LR->SetMarkerStyle(21);
   hnd2CVNNCPi0LR->SetMarkerColor(kOrange);
   hs5->Add(hnd2CVNNCPi0LR);

   //TH1D *c2_hndCCPi0 = new TH1D("c2_ndCCPi0","ndCCPi0",100,-4,4);
   hnd2CVNCCPi0R->SetLineColor(kBlack);
   hnd2CVNCCPi0R->SetFillColor(kMagenta);
   hnd2CVNCCPi0R->SetLineWidth(1);;
   hnd2CVNCCPi0R->SetMarkerStyle(21);
   hnd2CVNCCPi0R->SetMarkerColor(kMagenta);

   hs5->Add(hnd2CVNCCPi0R);
   //TH1D *c2_hndCCNonPi0 = new TH1D("c2_hndCCNonPi0","hndCCNonPi0",100,-4,4);
   hnd2CVNCCNonPi0R->SetLineColor(kBlack);
   hnd2CVNCCNonPi0R->SetFillColor(kAzure);
   hnd2CVNCCNonPi0R->SetLineWidth(1);;
   hnd2CVNCCNonPi0R->SetMarkerStyle(21);
   hnd2CVNCCNonPi0R->SetMarkerColor(kAzure);
   hs5->Add(hnd2CVNCCNonPi0R);
   
   //TH1D *c2_hndNCPi0 = new TH1D("c2_ndNCPi0","ndNCPi0",100,-4,4);
   hnd2CVNNCNonPi0R->SetLineColor(kBlack);
   hnd2CVNNCNonPi0R->SetFillColor(kCyan);
   hnd2CVNNCNonPi0R->SetLineWidth(1);;
   hnd2CVNNCNonPi0R->SetMarkerStyle(21);
   hnd2CVNNCNonPi0R->SetMarkerColor(kCyan);
   hs5->Add(hnd2CVNNCNonPi0R);
   
   //TH1D *c2_hndNCPi0L = new TH1D("c2_ndNCPi0L","ndNCPi0L",100,-4,4);
   hnd2CVNNCPi0R->SetLineColor(kBlack);
   hnd2CVNNCPi0R->SetFillColor(kGreen);
   hnd2CVNNCPi0R->SetLineWidth(1);;
   hnd2CVNNCPi0R->SetMarkerStyle(21);
   hnd2CVNNCPi0R->SetMarkerColor(kGreen);
   hs5->Add(hnd2CVNNCPi0R);

   //TCanvas *cs2 = new TCanvas("cs2","cs2",10,10,855,500);
   TText T5; T5.SetTextFont(42); T5.SetTextAlign(21);
   hs5->Draw("hist"); T5.DrawTextNDC(.5,.95,"");
   hs5->GetYaxis()->SetTitle("Events/ 8e20 POT");
   hs5->GetXaxis()->SetTitle("Prong 2 CVN Gamma ID");
   hs5->SetMaximum(14500);
   hs5->GetXaxis()->SetRangeUser(0,1);
   gPad->Modified();
   cs5->cd();
   //cs2->SetTopMargin(0.18);
   hnd2CVNMC3PR->SetLineColor(kRed);
   hnd2CVNMC3PR->SetLineWidth(3);
   //hndNumi2CVN2PRData->SetMarkerStyle(2);
   hndNumi2CVN2PRData->SetMarkerSize(3);

   hnd2CVNMC3PR->Draw("same");
   //hndNumi2CVN2PRData->Draw("same");
   TLegend * leg5 = new TLegend(0.55,0.5,0.65,0.82);
   //leg2->SetTextAlign(40);
   leg5->SetTextFont(42);
   leg5->SetTextSize(0.05);
   leg5->SetFillStyle(0);
   //leg5->AddEntry(hndNumi2CVN2PRData,"Data","kFullDotSmall");
   leg5->AddEntry(hnd2CVNMC3PR,"Selection","kFullDotSmall");
   leg5->AddEntry(hnd2CVNCCPi0R,"CCPi0","kFullDotSmall");
   leg5->AddEntry(hnd2CVNCCNonPi0R,"CCNonPi0","kFullDotSmall");
   leg5->AddEntry(hnd2CVNNCPi0LR,"NC w/Pi0 < 0.3 E","kFullDotSmall");
   leg5->AddEntry(hnd2CVNNCPi0R,"NC w/Pi0 > 0.3 E","kFullDotSmall");
   leg5->AddEntry(hnd2CVNNCNonPi0R,"NCNonPi0","kFullDotSmall");
   
   leg5->Draw("same"); 
    NDSimulation();
    //cs5->SaveAs("/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/NcPi0Plots/plot6_2Pr_0_3_CVN07__11-02-17.pdf");
    cs5->SaveAs(outDir2+CurrentFile[6].name.c_str()+CurrentDate[0].date.c_str()+".png");


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

/// Plot 6 ///
  //gStyle->SetTitleY(0.9);
   TCanvas *cs6 = new TCanvas("cs6","cs6",10,10,855,500);
   THStack *hs6 = new THStack("hs6","");
 //TCanvas *c2_ndMC =new TCanvas("c2_ndMC","ndMC",800,500);
   //TH1D *c2_ndNCMCSig = new TH1D("c2_ndNCMCSig","ndNCMCSig",100,-4,4);
   // hndNCMCSig->SetLineColor(kBlack);
   // hndNCMCSig->SetFillColor(kGray);
   // hndNCMCSig->SetLineWidth(3);
   // hndNCMCSig->SetMarkerStyle(21);
   // hndNCMCSig->SetMarkerColor(kGray);
   
   // hs2->Add(hndNCMCSig);
  
   //TH1D *c2_ndCCMCSig = new TH1D("c2_ndCCMCSig","ndCCMCSig",100,-4,4);
   hndNCPi0LRem->SetLineColor(kBlack);
   hndNCPi0LRem->SetFillColor(kOrange);
   hndNCPi0LRem->SetLineWidth(1);;
   hndNCPi0LRem->SetMarkerStyle(21);
   hndNCPi0LRem->SetMarkerColor(kOrange);
   hs6->Add(hndNCPi0LRem);

   //TH1D *c2_hndCCPi0 = new TH1D("c2_ndCCPi0","ndCCPi0",100,-4,4);
   hndCCPi0Rem->SetLineColor(kBlack);
   hndCCPi0Rem->SetFillColor(kMagenta);
   hndCCPi0Rem->SetLineWidth(1);;
   hndCCPi0Rem->SetMarkerStyle(21);
   hndCCPi0Rem->SetMarkerColor(kMagenta);

   hs6->Add(hndCCPi0Rem);
   //TH1D *c2_hndCCNonPi0 = new TH1D("c2_hndCCNonPi0","hndCCNonPi0",100,-4,4);
   hndCCNonPi0Rem->SetLineColor(kBlack);
   hndCCNonPi0Rem->SetFillColor(kAzure);
   hndCCNonPi0Rem->SetLineWidth(1);;
   hndCCNonPi0Rem->SetMarkerStyle(21);
   hndCCNonPi0Rem->SetMarkerColor(kAzure);
   hs6->Add(hndCCNonPi0Rem);
   
   //TH1D *c2_hndNCPi0 = new TH1D("c2_ndNCPi0","ndNCPi0",100,-4,4);
   hndNCNonPi0Rem->SetLineColor(kBlack);
   hndNCNonPi0Rem->SetFillColor(kCyan);
   hndNCNonPi0Rem->SetLineWidth(1);;
   hndNCNonPi0Rem->SetMarkerStyle(21);
   hndNCNonPi0Rem->SetMarkerColor(kCyan);
   hs6->Add(hndNCNonPi0Rem);
   
   //TH1D *c2_hndNCPi0L = new TH1D("c2_ndNCPi0L","ndNCPi0L",100,-4,4);
   hndNCPi0Rem->SetLineColor(kBlack);
   hndNCPi0Rem->SetFillColor(kGreen);
   hndNCPi0Rem->SetLineWidth(1);;
   hndNCPi0Rem->SetMarkerStyle(21);
   hndNCPi0Rem->SetMarkerColor(kGreen);
   hs6->Add(hndNCPi0Rem);

   //TCanvas *cs2 = new TCanvas("cs2","cs2",10,10,855,500);
   TText T6; T6.SetTextFont(42); T6.SetTextAlign(21);
   hs6->Draw("hist"); T.DrawTextNDC(.5,.95,"");
   hs6->GetYaxis()->SetTitle("Events/ 8e20 POT");
   hs6->GetXaxis()->SetTitle("RemID");
   hs6->SetMaximum(17000);
   hs6->GetXaxis()->SetRangeUser(0,1);
   gPad->Modified();
   cs6->cd();
   //cs2->SetTopMargin(0.18);
   hndMC3PRem->SetLineColor(kRed);
   hndMC3PRem->SetLineWidth(3);
   //hndNumiRemIDData->SetMarkerStyle(2);
   hndNumiRemIDData->SetMarkerSize(10);


   hndMC3PRem->Draw("same");
   //hndNumiRemIDData->Draw("same");
   TLegend * leg6 = new TLegend(0.55,0.5,0.65,0.82);
   //leg2->SetTextAlign(40);
   leg6->SetTextFont(42);
   leg6->SetTextSize(0.05);
   leg6->SetFillStyle(0);
   //leg6->AddEntry(hndNumiRemIDData,"Data","kFullDotSmall");
   leg6->AddEntry(hndMC3PRem,"Selection","kFullDotSmall");
   leg6->AddEntry(hndCCPi0Rem,"CCPi0","kFullDotSmall");
   leg6->AddEntry(hndCCNonPi0Rem,"CCNonPi0","kFullDotSmall");
   leg6->AddEntry(hndNCPi0LRem,"NC w/Pi0 < 0.3 E","kFullDotSmall");
   leg6->AddEntry(hndNCPi0Rem,"NC w/Pi0 > 0.3 E","kFullDotSmall");
   leg6->AddEntry(hndNCNonPi0Rem,"NCNonPi0","kFullDotSmall");
   
   leg6->Draw("same"); 
    NDSimulation();
    //cs6->SaveAs("/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/NcPi0Plots/plot7_2Pr_0_3_CVN07__RemID_11-02-17.pdf");
    cs6->SaveAs(outDir2+CurrentFile[7].name.c_str()+CurrentDate[0].date.c_str()+".png");

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

/// Plot 12 ///
  //gStyle->SetTitleY(0.9);
   TCanvas *cs7 = new TCanvas("cs7","cs7",10,10,855,500);
   THStack *hs7 = new THStack("hs7","");
 //TCanvas *c2_ndMC =new TCanvas("c2_ndMC","ndMC",800,500);
   //TH1D *c2_ndNCMCSig = new TH1D("c2_ndNCMCSig","ndNCMCSig",100,-4,4);
   // hndNCMCSig->SetLineColor(kBlack);
   // hndNCMCSig->SetFillColor(kGray);
   // hndNCMCSig->SetLineWidth(1);;
   // hndNCMCSig->SetMarkerStyle(21);
   // hndNCMCSig->SetMarkerColor(kGray);
   
   // hs2->Add(hndNCMCSig);

   //TH1D *c2_hndNCPi0L = new TH1D("c2_ndNCPi0L","ndNCPi0L",100,-4,4);
   hndNCPi0Prel->SetLineColor(kBlack);
   hndNCPi0Prel->SetFillColor(kGreen);
   hndNCPi0Prel->SetLineWidth(1);;
   hndNCPi0Prel->SetMarkerStyle(21);
   hndNCPi0Prel->SetMarkerColor(kGreen);
   hs7->Add(hndNCPi0Prel);
    
   //TH1D *c2_ndCCMCSig = new TH1D("c2_ndCCMCSig","ndCCMCSig",100,-4,4);
   hndNCPi0LPrel->SetLineColor(kBlack);
   hndNCPi0LPrel->SetFillColor(kOrange);
   hndNCPi0LPrel->SetLineWidth(1);;
   hndNCPi0LPrel->SetMarkerStyle(21);
   hndNCPi0LPrel->SetMarkerColor(kOrange);
   hs7->Add(hndNCPi0LPrel);

   
   //TH1D *c2_hndNCPi0 = new TH1D("c2_ndNCPi0","ndNCPi0",100,-4,4);
   hndNCNonPi0Prel->SetLineColor(kBlack);
   hndNCNonPi0Prel->SetFillColor(kCyan);
   hndNCNonPi0Prel->SetLineWidth(1);;
   hndNCNonPi0Prel->SetMarkerStyle(21);
   hndNCNonPi0Prel->SetMarkerColor(kCyan);
   hs7->Add(hndNCNonPi0Prel);
   
   //TH1D *c2_hndCCPi0 = new TH1D("c2_ndCCPi0","ndCCPi0",100,-4,4);
   hndCCPi0Prel->SetLineColor(kBlack);
   hndCCPi0Prel->SetFillColor(kMagenta);
   hndCCPi0Prel->SetLineWidth(1);;
   hndCCPi0Prel->SetMarkerStyle(21);
   hndCCPi0Prel->SetMarkerColor(kMagenta);

   hs7->Add(hndCCPi0Prel);
  
   //TH1D *c2_hndCCNonPi0 = new TH1D("c2_hndCCNonPi0","hndCCNonPi0",100,-4,4);
   hndCCNonPi0Prel->SetLineColor(kBlack);
   hndCCNonPi0Prel->SetFillColor(kAzure);
   hndCCNonPi0Prel->SetLineWidth(1);;
   hndCCNonPi0Prel->SetMarkerStyle(21);
   hndCCNonPi0Prel->SetMarkerColor(kAzure);
   hs7->Add(hndCCNonPi0Prel);

   //TCanvas *cs2 = new TCanvas("cs2","cs2",10,10,855,500);
   TText T7; T7.SetTextFont(42); T7.SetTextAlign(21);
   hs7->Draw("hist"); T7.DrawTextNDC(.5,.95,"");
   hs7->GetYaxis()->SetTitle("Events/ 8e20 POT");
   hs7->GetXaxis()->SetTitle("#pi^{0} True Energy [GeV]");
   hs7->SetMaximum(12000);
   // hs6->GetXaxis()->SetRangeUser(0.,0.3);
   gPad->Modified();
   cs7->cd();
   //cs2->SetTopMargin(0.18);
   hndMC3PPrel->SetLineColor(kRed);
   hndMC3PPrel->SetLineWidth(3);
   hndMC3PPrel->Draw("same");
   TLegend * leg7 = new TLegend(0.55,0.5,0.65,0.82);
   //leg2->SetTextAlign(40);
   leg7->SetTextFont(42);
   leg7->SetTextSize(0.05);
   leg7->SetFillStyle(0);
   leg7->AddEntry(hndMC3PPrel,"Selection","kFullDotSmall");
   // leg2->AddEntry(hndNCMCSig,"Signal","kFullDotSmall");
  
  
   leg7->AddEntry(hndNCPi0Prel,"NC w/Pi0 > 0.3 E","kFullDotSmall");
   leg7->AddEntry(hndNCPi0LPrel,"NC w/Pi0 < 0.3 E","kFullDotSmall");
   leg7->AddEntry(hndNCNonPi0Prel,"NCNonPi0","kFullDotSmall");
   leg7->AddEntry(hndCCPi0Prel,"CCPi0","kFullDotSmall");
   leg7->AddEntry(hndCCNonPi0Prel,"CCNonPi0","kFullDotSmall");
   
   leg7->Draw("same"); 
    NDSimulation();
    //cs7->SaveAs("/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/NcPi0Plots/plot8_2Pr_0_3_CVN07__TruthEn_11-02-17.pdf");
    cs7->SaveAs(outDir2+CurrentFile[8].name.c_str()+CurrentDate[0].date.c_str()+".png");


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

/// Plot 13 ///
  //gStyle->SetTitleY(0.9);
   TCanvas *cs8 = new TCanvas("cs8","cs8",10,10,855,500);
   THStack *hs8 = new THStack("hs8","");
 //TCanvas *c2_ndMC =new TCanvas("c2_ndMC","ndMC",800,500);
   //TH1D *c2_ndNCMCSig = new TH1D("c2_ndNCMCSig","ndNCMCSig",100,-4,4);
   // hndNCMCSig->SetLineColor(kBlack);
   // hndNCMCSig->SetFillColor(kGray);
   // hndNCMCSig->SetLineWidth(3);
   // hndNCMCSig->SetMarkerStyle(21);
   // hndNCMCSig->SetMarkerColor(kGray);
   
   // hs2->Add(hndNCMCSig);

   //TH1D *c2_hndNCPi0L = new TH1D("c2_ndNCPi0L","ndNCPi0L",100,-4,4);
   hndNCPi0PrelR->SetLineColor(kBlack);
   hndNCPi0PrelR->SetFillColor(kGreen);
   hndNCPi0PrelR->SetLineWidth(1);;
   hndNCPi0PrelR->SetMarkerStyle(21);
   hndNCPi0PrelR->SetMarkerColor(kGreen);
   hs8->Add(hndNCPi0PrelR);
    
   //TH1D *c2_ndCCMCSig = new TH1D("c2_ndCCMCSig","ndCCMCSig",100,-4,4);
   hndNCPi0LPrelR->SetLineColor(kBlack);
   hndNCPi0LPrelR->SetFillColor(kOrange);
   hndNCPi0LPrelR->SetLineWidth(1);;
   hndNCPi0LPrelR->SetMarkerStyle(21);
   hndNCPi0LPrelR->SetMarkerColor(kOrange);
   hs8->Add(hndNCPi0LPrelR);

   
   //TH1D *c2_hndNCPi0 = new TH1D("c2_ndNCPi0","ndNCPi0",100,-4,4);
   hndNCNonPi0PrelR->SetLineColor(kBlack);
   hndNCNonPi0PrelR->SetFillColor(kCyan);
   hndNCNonPi0PrelR->SetLineWidth(1);;
   hndNCNonPi0PrelR->SetMarkerStyle(21);
   hndNCNonPi0PrelR->SetMarkerColor(kCyan);
   hs8->Add(hndNCNonPi0PrelR);
   
   //TH1D *c2_hndCCPi0 = new TH1D("c2_ndCCPi0","ndCCPi0",100,-4,4);
   hndCCPi0PrelR->SetLineColor(kBlack);
   hndCCPi0PrelR->SetFillColor(kMagenta);
   hndCCPi0PrelR->SetLineWidth(1);;
   hndCCPi0PrelR->SetMarkerStyle(21);
   hndCCPi0PrelR->SetMarkerColor(kMagenta);

   hs8->Add(hndCCPi0PrelR);
  
   //TH1D *c2_hndCCNonPi0 = new TH1D("c2_hndCCNonPi0","hndCCNonPi0",100,-4,4);
   hndCCNonPi0PrelR->SetLineColor(kBlack);
   hndCCNonPi0PrelR->SetFillColor(kAzure);
   hndCCNonPi0PrelR->SetLineWidth(1);;
   hndCCNonPi0PrelR->SetMarkerStyle(21);
   hndCCNonPi0PrelR->SetMarkerColor(kAzure);
   hs8->Add(hndCCNonPi0PrelR);

   //TCanvas *cs2 = new TCanvas("cs2","cs2",10,10,855,500);
   TText T8; T8.SetTextFont(42); T8.SetTextAlign(21);
   hs8->Draw("hist"); T8.DrawTextNDC(.5,.95,"");
   hs8->GetYaxis()->SetTitle("Events/ 8e20 POT");
   hs8->GetXaxis()->SetTitle("#pi^{0} True Energy [GeV]");
   hs8->SetMaximum(9800);
   // hs6->GetXaxis()->SetRangeUser(0.,0.3);
   gPad->Modified();
   cs8->cd();
   //cs2->SetTopMargin(0.18);
   hndMC3PPrelR->SetLineColor(kRed);
   hndMC3PPrelR->SetLineWidth(1);;
   hndMC3PPrelR->Draw("same");
   TLegend * leg8 = new TLegend(0.55,0.5,0.65,0.82);
   //leg2->SetTextAlign(40);
   leg8->SetTextFont(42);
   leg8->SetTextSize(0.05);
   leg8->SetFillStyle(0);
   leg8->AddEntry(hndMC3PPrel,"Selection","kFullDotSmall");
   // leg2->AddEntry(hndNCMCSig,"Signal","kFullDotSmall");
  
  
   leg8->AddEntry(hndNCPi0PrelR,"NC w/Pi0 > 0.3 E","kFullDotSmall");
   leg8->AddEntry(hndNCPi0LPrelR,"NC w/Pi0 < 0.3 E","kFullDotSmall");
   leg8->AddEntry(hndNCNonPi0PrelR,"NCNonPi0","kFullDotSmall");
   leg8->AddEntry(hndCCPi0PrelR,"CCPi0","kFullDotSmall");
   leg8->AddEntry(hndCCNonPi0PrelR,"CCNonPi0","kFullDotSmall");
   
   leg8->Draw("same"); 
    NDSimulation();
    //cs8->SaveAs("/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/NcPi0Plots/plot9_2Pr_0_3_CVN07__TruthEnR_11-02-17.pdf");
    cs8->SaveAs(outDir2+CurrentFile[9].name.c_str()+CurrentDate[0].date.c_str()+".png");

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

/// Plot 14 ///
  //gStyle->SetTitleY(0.9);
   TCanvas *cs9 = new TCanvas("cs9","cs9",10,10,855,500);
   THStack *hs9 = new THStack("hs9","");
 //TCanvas *c2_ndMC =new TCanvas("c2_ndMC","ndMC",800,500);
   //TH1D *c2_ndNCMCSig = new TH1D("c2_ndNCMCSig","ndNCMCSig",100,-4,4);
   // hndNCMCSig->SetLineColor(kBlack);
   // hndNCMCSig->SetFillColor(kGray);
   // hndNCMCSig->SetLineWidth(3);
   // hndNCMCSig->SetMarkerStyle(21);
   // hndNCMCSig->SetMarkerColor(kGray);
   
   // hs2->Add(hndNCMCSig);

   //TH1D *c2_hndNCPi0L = new TH1D("c2_ndNCPi0L","ndNCPi0L",100,-4,4);
   hndNCPi0PrelCVN->SetLineColor(kBlack);
   hndNCPi0PrelCVN->SetFillColor(kGreen);
   hndNCPi0PrelCVN->SetLineWidth(1);;
   hndNCPi0PrelCVN->SetMarkerStyle(21);
   hndNCPi0PrelCVN->SetMarkerColor(kGreen);
   hs9->Add(hndNCPi0PrelCVN);
    
   //TH1D *c2_ndCCMCSig = new TH1D("c2_ndCCMCSig","ndCCMCSig",100,-4,4);
   hndNCPi0LPrelCVN->SetLineColor(kBlack);
   hndNCPi0LPrelCVN->SetFillColor(kOrange);
   hndNCPi0LPrelCVN->SetLineWidth(1);;
   hndNCPi0LPrelCVN->SetMarkerStyle(21);
   hndNCPi0LPrelCVN->SetMarkerColor(kOrange);
   hs9->Add(hndNCPi0LPrelCVN);

   
   //TH1D *c2_hndNCPi0 = new TH1D("c2_ndNCPi0","ndNCPi0",100,-4,4);
   hndNCNonPi0PrelCVN->SetLineColor(kBlack);
   hndNCNonPi0PrelCVN->SetFillColor(kCyan);
   hndNCNonPi0PrelCVN->SetLineWidth(1);;
   hndNCNonPi0PrelCVN->SetMarkerStyle(21);
   hndNCNonPi0PrelCVN->SetMarkerColor(kCyan);
   hs9->Add(hndNCNonPi0PrelCVN);
   
   //TH1D *c2_hndCCPi0 = new TH1D("c2_ndCCPi0","ndCCPi0",100,-4,4);
   hndCCPi0PrelCVN->SetLineColor(kBlack);
   hndCCPi0PrelCVN->SetFillColor(kMagenta);
   hndCCPi0PrelCVN->SetLineWidth(1);;
   hndCCPi0PrelCVN->SetMarkerStyle(21);
   hndCCPi0PrelCVN->SetMarkerColor(kMagenta);

   hs9->Add(hndCCPi0PrelCVN);
  
   //TH1D *c2_hndCCNonPi0 = new TH1D("c2_hndCCNonPi0","hndCCNonPi0",100,-4,4);
   hndCCNonPi0PrelCVN->SetLineColor(kBlack);
   hndCCNonPi0PrelCVN->SetFillColor(kAzure);
   hndCCNonPi0PrelCVN->SetLineWidth(1);;
   hndCCNonPi0PrelCVN->SetMarkerStyle(21);
   hndCCNonPi0PrelCVN->SetMarkerColor(kAzure);
   hs9->Add(hndCCNonPi0PrelCVN);

   //TCanvas *cs2 = new TCanvas("cs2","cs2",10,10,855,500);
   TText T9; T9.SetTextFont(42); T9.SetTextAlign(21);
   hs9->Draw("hist"); T9.DrawTextNDC(.5,.95,"");
   hs9->GetYaxis()->SetTitle("Events/ 8e20 POT");
   hs9->GetXaxis()->SetTitle("#pi^{0} True Energy [GeV]");
   hs9->SetMaximum(4000);
   // hs6->GetXaxis()->SetRangeUser(0.,0.3);
   gPad->Modified();
   cs9->cd();
   //cs2->SetTopMargin(0.18);
   hndMC3PPrelCVN->SetLineColor(kRed);
   hndMC3PPrelCVN->SetLineWidth(3);
   hndMC3PPrelCVN->Draw("same");
   TLegend * leg9 = new TLegend(0.55,0.5,0.65,0.82);
   //leg2->SetTextAlign(40);
   leg9->SetTextFont(42);
   leg9->SetTextSize(0.05);
   leg9->SetFillStyle(0);
   leg9->AddEntry(hndMC3PPrelCVN,"Selection","kFullDotSmall");
   // leg2->AddEntry(hndNCMCSig,"Signal","kFullDotSmall");
  
  
   leg9->AddEntry(hndNCPi0PrelCVN,"NC w/Pi0 > 0.3 E","kFullDotSmall");
   leg9->AddEntry(hndNCPi0LPrelCVN,"NC w/Pi0 < 0.3 E","kFullDotSmall");
   leg9->AddEntry(hndNCNonPi0PrelCVN,"NCNonPi0","kFullDotSmall");
   leg9->AddEntry(hndCCPi0PrelCVN,"CCPi0","kFullDotSmall");
   leg9->AddEntry(hndCCNonPi0PrelCVN,"CCNonPi0","kFullDotSmall");
   
   leg9->Draw("same"); 
    NDSimulation();
    //cs9->SaveAs("/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/NcPi0Plots/plot10_2Pr_0_3_CVN07__TruthEnCVN_11-02-17.pdf");
    cs9->SaveAs(outDir2+CurrentFile[10].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /*********************************************************************************************************************/
/// Plot 14 ///
  //gStyle->SetTitleY(0.9);
   TCanvas *cs9b = new TCanvas("cs9b","cs9b",10,10,855,500);
   THStack *hs9b = new THStack("hs9b","");
 //TCanvas *c2_ndMC =new TCanvas("c2_ndMC","ndMC",800,500);
   //TH1D *c2_ndNCMCSig = new TH1D("c2_ndNCMCSig","ndNCMCSig",100,-4,4);
   // hndNCMCSig->SetLineColor(kBlack);
   // hndNCMCSig->SetFillColor(kGray);
   // hndNCMCSig->SetLineWidth(3);
   // hndNCMCSig->SetMarkerStyle(21);
   // hndNCMCSig->SetMarkerColor(kGray);
   
   // hs2->Add(hndNCMCSig);

   //TH1D *c2_hndNCPi0L = new TH1D("c2_ndNCPi0L","ndNCPi0L",100,-4,4);
   hndNCPi0PrelCVN2Cuts->SetLineColor(kBlack);
   hndNCPi0PrelCVN2Cuts->SetFillColor(kGreen);
   hndNCPi0PrelCVN2Cuts->SetLineWidth(1);;
   hndNCPi0PrelCVN2Cuts->SetMarkerStyle(21);
   hndNCPi0PrelCVN2Cuts->SetMarkerColor(kGreen);
   hs9b->Add(hndNCPi0PrelCVN2Cuts);
    
   //TH1D *c2_ndCCMCSig = new TH1D("c2_ndCCMCSig","ndCCMCSig",100,-4,4);
   hndNCPi0LPrelCVN2Cuts->SetLineColor(kBlack);
   hndNCPi0LPrelCVN2Cuts->SetFillColor(kOrange);
   hndNCPi0LPrelCVN2Cuts->SetLineWidth(1);;
   hndNCPi0LPrelCVN2Cuts->SetMarkerStyle(21);
   hndNCPi0LPrelCVN2Cuts->SetMarkerColor(kOrange);
   hs9b->Add(hndNCPi0LPrelCVN2Cuts);

   
   //TH1D *c2_hndNCPi0 = new TH1D("c2_ndNCPi0","ndNCPi0",100,-4,4);
   hndNCNonPi0PrelCVN2Cuts->SetLineColor(kBlack);
   hndNCNonPi0PrelCVN2Cuts->SetFillColor(kCyan);
   hndNCNonPi0PrelCVN2Cuts->SetLineWidth(1);;
   hndNCNonPi0PrelCVN2Cuts->SetMarkerStyle(21);
   hndNCNonPi0PrelCVN2Cuts->SetMarkerColor(kCyan);
   hs9b->Add(hndNCNonPi0PrelCVN2Cuts);
   
   //TH1D *c2_hndCCPi0 = new TH1D("c2_ndCCPi0","ndCCPi0",100,-4,4);
   hndCCPi0PrelCVN2Cuts->SetLineColor(kBlack);
   hndCCPi0PrelCVN2Cuts->SetFillColor(kMagenta);
   hndCCPi0PrelCVN2Cuts->SetLineWidth(1);;
   hndCCPi0PrelCVN2Cuts->SetMarkerStyle(21);
   hndCCPi0PrelCVN2Cuts->SetMarkerColor(kMagenta);

   hs9b->Add(hndCCPi0PrelCVN2Cuts);
  
   //TH1D *c2_hndCCNonPi0 = new TH1D("c2_hndCCNonPi0","hndCCNonPi0",100,-4,4);
   hndCCNonPi0PrelCVN2Cuts->SetLineColor(kBlack);
   hndCCNonPi0PrelCVN2Cuts->SetFillColor(kAzure);
   hndCCNonPi0PrelCVN2Cuts->SetLineWidth(1);;
   hndCCNonPi0PrelCVN2Cuts->SetMarkerStyle(21);
   hndCCNonPi0PrelCVN2Cuts->SetMarkerColor(kAzure);
   hs9b->Add(hndCCNonPi0PrelCVN2Cuts);

   //TCanvas *cs2 = new TCanvas("cs2","cs2",10,10,855,500);
   TText T91; T91.SetTextFont(42); T91.SetTextAlign(21);
   hs9b->Draw("hist"); T9.DrawTextNDC(.5,.95,"");
   hs9b->GetYaxis()->SetTitle("Events/ 8e20 POT");
   hs9b->GetXaxis()->SetTitle("#pi^{0} True Energy [GeV]");
   hs9b->SetMaximum(1100);
   // hs6->GetXaxis()->SetRangeUser(0.,0.3);
   gPad->Modified();
   cs9b->cd();
   //cs2->SetTopMargin(0.18);
   hndMC3PPrelCVN2Cuts->SetLineColor(kRed);
   hndMC3PPrelCVN2Cuts->SetLineWidth(3);
   //hndMC3PPrelCVN2Cuts->SetMarkerStyle(2);
   hndMC3PPrelCVN2Cuts->Draw("same");
   TLegend * leg9b = new TLegend(0.55,0.5,0.65,0.82);
   //leg2->SetTextAlign(40);
   leg9b->SetTextFont(42);
   leg9b->SetTextSize(0.05);
   leg9b->SetFillStyle(0);
   leg9b->AddEntry(hndMC3PPrelCVN2Cuts,"Selection","kFullDotSmall");
   // leg2->AddEntry(hndNCMCSig,"Signal","kFullDotSmall");
  
  
   leg9b->AddEntry(hndNCPi0PrelCVN2Cuts,"NC w/Pi0 > 0.3 E","kFullDotSmall");
   leg9b->AddEntry(hndNCPi0LPrelCVN2Cuts,"NC w/Pi0 < 0.3 E","kFullDotSmall");
   leg9b->AddEntry(hndNCNonPi0PrelCVN2Cuts,"NCNonPi0","kFullDotSmall");
   leg9b->AddEntry(hndCCPi0PrelCVN2Cuts,"CCPi0","kFullDotSmall");
   leg9b->AddEntry(hndCCNonPi0PrelCVN2Cuts,"CCNonPi0","kFullDotSmall");
   
   leg9b->Draw("same"); 
    NDSimulation();
    //cs9->SaveAs("/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/NcPi0Plots/plot10_2Pr_0_3_CVN07__TruthEnCVN_11-02-17.pdf");
    cs9b->SaveAs(outDir2+CurrentFile[11].name.c_str()+CurrentDate[0].date.c_str()+".png");

 /*********************************************************************************************************************/
   /// Plot 16 ///
   TCanvas *cs10 = new TCanvas("cs14","cs14",10,10,855,500);
   cs10->cd();
   hnd2dplot2prong->Draw("colz");
   hnd2dplot2prong->GetYaxis()->SetTitle("#pi^{0} True Energy [GeV]");
   hnd2dplot2prong->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prong->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("/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/NcPi0Plots/plot11_2Pr_2dplot_0_3_CVN07__11-02-17.pdf");
   cs10->SaveAs(outDir2+CurrentFile[12].name.c_str()+CurrentDate[0].date.c_str()+".png");
  
// /*********************************************************************************************************************/
/*********************************************************************************************************************/
   /// Plot 18 ///
   TCanvas *cs10b = new TCanvas("cs10b","cs10b",10,10,855,500);
   cs10b->cd();
   gPad->SetLogz(1);
   hnd2dplot2prong->Draw("colz");
   hnd2dplot2prong->GetYaxis()->SetTitle("#pi^{0} True Energy [GeV]");
   hnd2dplot2prong->GetXaxis()->SetTitle("#pi^{0} Reconstructed Mass [GeV]");
   hnd2dplot2prong->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("/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/NcPi0Plots/plot11_2Pr_2dplot_0_3_CVN07__11-02-17.pdf");
   cs10b->SaveAs(outDir2+CurrentFile[15].name.c_str()+CurrentDate[0].date.c_str()+".png");
  
// /*********************************************************************************************************************/
/// Plot 2 ///
 
   //gStyle->SetTitleY(0.9);
   TCanvas *cs13 = new TCanvas("cs13","cs13",10,10,855,500);
   cs13->cd();

   hndMC3P1->SetLineColor(kRed);
   //hndMC3P1->SetLineWidth(3);
   // hndMC3P1->SetMarkerStyle(21);
   //hndMC3P1->SetMarkerColor(kRed);
  
   hndCCPi01->SetLineColor(kMagenta+2);
   //hndCCPi01->SetLineWidth(3);
   // hndCCPi01->SetMarkerStyle(21);
   // hndCCPi01->SetMarkerColor(kMagenta);

   hndCCNonPi01->SetLineColor(kAzure+2);
   //hndCCNonPi01->SetLineWidth(3);
   // hndCCNonPi01->SetMarkerStyle(21);
   //hndCCNonPi01->SetMarkerColor(kAzure);

   hndNCPi0L1->SetLineColor(kOrange+2);
   //hndNCPi0L1->SetLineWidth(3);
   // hndNCPi0L1->SetMarkerStyle(21);
   //hndNCPi0L1->SetMarkerColor(kOrange);
  
   hndNCPi01->SetLineColor(kGreen+2); 
   //hndNCPi01->SetLineWidth(3);
   // hndNCPi01->SetMarkerStyle(21);
   //hndNCPi01->SetMarkerColor(kGreen);
  

   hndNCNonPi01->SetLineColor(kCyan-6);
   //hndNCNonPi01->SetLineWidth(3);
   // hndNCNonPi01->SetMarkerStyle(21);
   //hndNCNonPi01->SetMarkerColor(kCyan);
   //hndNumiMassdata2CVN1->SetMarkerStyle(2);
   hndNumiMassdata2CVN1->SetMarkerSize(10);


   hndMC3P1->Draw("hist E1");
   //hndNumiMassdata2CVN1->Draw("hist E1 same");
   hndCCPi01->Draw("hist E1 same");
   hndCCNonPi01->Draw("hist E1 same");
   hndNCPi0L1->Draw("hist E1 same");
   hndNCPi01->Draw("hist E1 same");
   hndNCNonPi01->Draw("hist E1 same");
   hndMC3P1->GetYaxis()->SetTitle("Events/8e20 POT");
   hndMC3P1->GetXaxis()->SetRangeUser(0.,0.4);
   hndMC3P1->SetTitleOffset(0.95,"Y");
   // hndMC3P1->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 ");
   
   TLegend * leg13 = new TLegend(0.62,0.5,0.73,0.82);
   leg13->SetTextFont(42);
   leg13->SetTextSize(0.05);
   //leg13->AddEntry(hndNumiMassdata2CVN1,"Data","l");                        
   leg13->AddEntry(hndMC3P1,"Selection","l");
   leg13->AddEntry(hndCCPi01,"CCPi0","l");
   leg13->AddEntry(hndCCNonPi01,"CCNonPi0","l");
   leg13->AddEntry(hndNCPi0L1,"NC w/Pi0 < 0.3 E","l");
   leg13->AddEntry(hndNCPi01,"NC w/Pi0 > 0.3 E","l");
   leg13->AddEntry(hndNCNonPi01,"NCNonPi0","l");
   
   leg13->Draw("same"); 
    NDSimulation();
    //cs13->SaveAs("/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/NcPi0Plots/plot13_2Pr_0_3_CVN07__11-02-17.pdf");
   cs13->SaveAs(outDir2+CurrentFile[13].name.c_str()+CurrentDate[0].date.c_str()+".png");

/*********************************************************************************************************************/
/// Plot 2 ///
 
   //gStyle->SetTitleY(0.9);
   TCanvas *cs13b = new TCanvas("cs13b","cs13b",10,10,855,500);
   cs13b->cd();

   hndMC3P1Cut2->SetLineColor(kRed);
   //hndMC3P1Cut2->SetLineWidth(3);
   // hndMC3P1->SetMarkerStyle(21);
   //hndMC3P1->SetMarkerColor(kRed);
  
   hndCCPi01Cut2->SetLineColor(kMagenta+2);
   //hndCCPi01Cut2->SetLineWidth(3);
   // hndCCPi01->SetMarkerStyle(21);
   // hndCCPi01->SetMarkerColor(kMagenta);

   hndCCNonPi01Cut2->SetLineColor(kAzure+2);
   //hndCCNonPi01Cut2->SetLineWidth(3);
   // hndCCNonPi01->SetMarkerStyle(21);
   //hndCCNonPi01->SetMarkerColor(kAzure);

   hndNCPi0L1Cut2->SetLineColor(kOrange+2);
   //hndNCPi0L1Cut2->SetLineWidth(3);
   // hndNCPi0L1->SetMarkerStyle(21);
   //hndNCPi0L1->SetMarkerColor(kOrange);
  
   hndNCPi01Cut2->SetLineColor(kGreen+2); 
   //hndNCPi01Cut2->SetLineWidth(3);
   // hndNCPi01->SetMarkerStyle(21);
   //hndNCPi01->SetMarkerColor(kGreen);
  

   hndNCNonPi01Cut2->SetLineColor(kCyan-6);
   //hndNCNonPi01Cut2->SetLineWidth(3);
   // hndNCNonPi01->SetMarkerStyle(21);
   //hndNCNonPi01->SetMarkerColor(kCyan);
   //hndNumiMassdataCVN2->SetMarkerStyle(2);
   hndNumiMassdataCVN2->SetMarkerSize(3);


   hndMC3P1Cut2->Draw("hist E1");
   //hndNumiMassdataCVN2->Draw("hist E1 same");
   hndCCPi01Cut2->Draw("hist E1 same");
   hndCCNonPi01Cut2->Draw("hist E1 same");
   hndNCPi0L1Cut2->Draw("hist E1 same");
   hndNCPi01Cut2->Draw("hist E1 same");
   hndNCNonPi01Cut2->Draw("hist E1 same");
   hndMC3P1Cut2->GetYaxis()->SetTitle("Events/8e20 POT");
   hndMC3P1Cut2->GetXaxis()->SetRangeUser(0.,0.4);
   hndMC3P1Cut2->SetTitleOffset(0.95,"Y");
   // hndMC3P1->SetTitle("2 Prong Event Preselection + RemID < 0.36 + Prong 1 & 2 CVN Gamma ID > 0.7 ");
   
   TLegend * leg13b = new TLegend(0.62,0.5,0.73,0.82);
   leg13b->SetTextFont(42);
   leg13b->SetTextSize(0.05);
   //leg13b->AddEntry(hndNumiMassdataCVN2,"Data","l");
   leg13b->AddEntry(hndMC3P1Cut2,"Selection","l");
   leg13b->AddEntry(hndCCPi01Cut2,"CCPi0","l");
   leg13b->AddEntry(hndCCNonPi01Cut2,"CCNonPi0","l");
   leg13b->AddEntry(hndNCPi0L1Cut2,"NC w/Pi0 < 0.3 E","l");
   leg13b->AddEntry(hndNCPi01Cut2,"NC w/Pi0 > 0.3 E","l");
   leg13b->AddEntry(hndNCNonPi01Cut2,"NCNonPi0","l");
   
   leg13b->Draw("same"); 
    NDSimulation();
    //cs13->SaveAs("/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/NcPi0Plots/plot13_2Pr_0_3_CVN07__11-02-17.pdf");
   cs13b->SaveAs(outDir2+CurrentFile[14].name.c_str()+CurrentDate[0].date.c_str()+".png");
/*********************************************************************************************************************/
 f1->Close();

}
#endif