EnergyCont_plot.C File Reference

#include "3FlavorAna/Ana2018/numu/SupportingStudies/EnergyContainment/EnergyCont_help.h"
#include "TCanvas.h"
#include "TFile.h"
#include "TH1D.h"
#include "TLatex.h"
#include "TLegend.h"
#include "TStyle.h"
#include "TSystem.h"
#include "CAFAna/Analysis/Exposures.h"
#include "CAFAna/Analysis/Plots.h"
#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Core/SpectrumLoader.h"
#include <iostream>
#include <sstream>
#include <fstream>
#include <vector>
#include <string>
#include <TROOT.h>

Go to the source code of this file.

Classes
struct	MySpectra

Functions
TH1D *	SpecToHist (TFile *MyF, std::string SpecName, double POT, std::string AxisTitle)
TH2 *	Spec2DtoHist (TFile *MyF, std::string SpecName, double POT, std::string AxisTitle)
void	MakeComparisonPlots (std::vector< TH1D * > TheHists, std::string PlotName, std::vector< std::string > Labels)
void	Print2D (TH2 *hist, std::string PlotName)
void	Print0binFrac (TH1 *hist, std::string PlotName)
void	EnergyCont_plot ()
void	GetPercentages (TH1 *hist, std::string PlotName)

Variables
const std::vector< std::string >	GENIEStr = { "All" }
const unsigned int	nGENIE = 1
std::vector< std::string >	FD_CutNames = { "FD_Contain", "FD_PID2018", "FD_CosRej" }
const unsigned int	nFD_CutNames = 3
std::vector< std::string >	ND_CutNames = { "ND_Contain", "ND_PID2018" }
const unsigned int	nND_CutNames = 2
unsigned int	FCut_FD = nFD_CutNames - 1
unsigned int	FCut_ND = nND_CutNames - 1
std::vector< std::string >	QuantNames = { "AllQuant", "Quant1", "Quant2", "Quant3", "Quant4" }
const unsigned int	nQuantNames = 5
std::vector< std::string >	HCNames = {"fhc", "rhc"}
const unsigned int	nHC = 2
std::vector< std::string >	DetNames = {"ND", "FD"}

Function Documentation

void EnergyCont_plot

(

)

Definition at line 253 of file EnergyCont_plot.C.

References MySpectra::AllUncontain_FD, MySpectra::AllUncontain_ND, cut, DetNames, MySpectra::ElectronEFrac_FD, MySpectra::ElectronEFrac_ND, MySpectra::ElectronEFracvsTrueE_FD, MySpectra::ElectronEFracvsTrueE_ND, MySpectra::ElectronEnergy_FD, MySpectra::ElectronEnergy_ND, MySpectra::ElecUncontain_FD, MySpectra::ElecUncontain_ND, FCut_FD, FCut_ND, FD_CutNames, compareCafs::files, GENIEStr, HCNames, ana::kAna2018FHCPOT, ana::kAna2018RHCPOT, MySpectra::KaonEFrac_FD, MySpectra::KaonEFrac_ND, MySpectra::KaonEFracvsTrueE_FD, MySpectra::KaonEFracvsTrueE_ND, MySpectra::KaonEnergy_FD, MySpectra::KaonEnergy_ND, MySpectra::KaonUncontain_FD, MySpectra::KaonUncontain_ND, MakeComparisonPlots(), MySpectra::MiscEnergy_FD, MySpectra::MiscEnergy_ND, MySpectra::MiscUncontain_FD, MySpectra::MiscUncontain_ND, MySpectra::MuonEFrac_FD, MySpectra::MuonEFrac_ND, MySpectra::MuonEFracvsTrueE_FD, MySpectra::MuonEFracvsTrueE_ND, MySpectra::MuonEnergy_FD, MySpectra::MuonEnergy_ND, MySpectra::MuonUncontain_FD, MySpectra::MuonUncontain_ND, ND_CutNames, MySpectra::NeutrinoEnergy_FD, MySpectra::NeutrinoEnergy_ND, MySpectra::NeutronEFrac_FD, MySpectra::NeutronEFrac_ND, MySpectra::NeutronEFracvsTrueE_FD, MySpectra::NeutronEFracvsTrueE_ND, MySpectra::NeutronEnergy_FD, MySpectra::NeutronEnergy_ND, MySpectra::NeutUncontain_FD, MySpectra::NeutUncontain_ND, nFD_CutNames, nGENIE, nHC, nND_CutNames, MySpectra::NNeut_FD, MySpectra::NNeut_ND, MySpectra::NPi0s_FD, MySpectra::NPi0s_ND, MySpectra::NPion_FD, MySpectra::NPion_ND, MySpectra::NProt_FD, MySpectra::NProt_ND, nQuantNames, MySpectra::PhotonEFrac_FD, MySpectra::PhotonEFrac_ND, MySpectra::PhotonEFracvsTrueE_FD, MySpectra::PhotonEFracvsTrueE_ND, MySpectra::PhotonEnergy_FD, MySpectra::PhotonEnergy_ND, MySpectra::PhotUncontain_FD, MySpectra::PhotUncontain_ND, MySpectra::Pi0EFrac_FD, MySpectra::Pi0EFrac_ND, MySpectra::Pi0EFracvsTrueE_FD, MySpectra::Pi0EFracvsTrueE_ND, MySpectra::Pi0Energy_FD, MySpectra::Pi0Energy_ND, MySpectra::Pi0sUncontain_FD, MySpectra::Pi0sUncontain_ND, MySpectra::PionEFrac_FD, MySpectra::PionEFrac_ND, MySpectra::PionEFracvsTrueE_FD, MySpectra::PionEFracvsTrueE_ND, MySpectra::PionEnergy_FD, MySpectra::PionEnergy_ND, MySpectra::PionUncontain_FD, MySpectra::PionUncontain_ND, POT, Print2D(), MySpectra::ProtonEFrac_FD, MySpectra::ProtonEFrac_ND, MySpectra::ProtonEFracvsTrueE_FD, MySpectra::ProtonEFracvsTrueE_ND, MySpectra::ProtonEnergy_FD, MySpectra::ProtonEnergy_ND, MySpectra::ProtUncontain_FD, MySpectra::ProtUncontain_ND, QuantNames, Spec2DtoHist(), SpecToHist(), string, MySpectra::UncontElecEF_FD, MySpectra::UncontElecEF_ND, MySpectra::UncontElecEFT_FD, MySpectra::UncontElecEFT_ND, MySpectra::UncontElecEFTvsTrueE_FD, MySpectra::UncontElecEFTvsTrueE_ND, MySpectra::UncontElecEFvsTrueE_FD, MySpectra::UncontElecEFvsTrueE_ND, MySpectra::UncontKaonEF_FD, MySpectra::UncontKaonEF_ND, MySpectra::UncontKaonEFT_FD, MySpectra::UncontKaonEFT_ND, MySpectra::UncontKaonEFTvsTrueE_FD, MySpectra::UncontKaonEFTvsTrueE_ND, MySpectra::UncontKaonEFvsTrueE_FD, MySpectra::UncontKaonEFvsTrueE_ND, MySpectra::UncontMuonEF_FD, MySpectra::UncontMuonEF_ND, MySpectra::UncontMuonEFT_FD, MySpectra::UncontMuonEFT_ND, MySpectra::UncontMuonEFTvsTrueE_FD, MySpectra::UncontMuonEFTvsTrueE_ND, MySpectra::UncontMuonEFvsTrueE_FD, MySpectra::UncontMuonEFvsTrueE_ND, MySpectra::UncontNeutEF_FD, MySpectra::UncontNeutEF_ND, MySpectra::UncontNeutEFT_FD, MySpectra::UncontNeutEFT_ND, MySpectra::UncontNeutEFTvsTrueE_FD, MySpectra::UncontNeutEFTvsTrueE_ND, MySpectra::UncontNeutEFvsTrueE_FD, MySpectra::UncontNeutEFvsTrueE_ND, MySpectra::UncontPhotEF_FD, MySpectra::UncontPhotEF_ND, MySpectra::UncontPhotEFT_FD, MySpectra::UncontPhotEFT_ND, MySpectra::UncontPhotEFTvsTrueE_FD, MySpectra::UncontPhotEFTvsTrueE_ND, MySpectra::UncontPhotEFvsTrueE_FD, MySpectra::UncontPhotEFvsTrueE_ND, MySpectra::UncontPi0sEF_FD, MySpectra::UncontPi0sEF_ND, MySpectra::UncontPi0sEFT_FD, MySpectra::UncontPi0sEFT_ND, MySpectra::UncontPi0sEFTvsTrueE_FD, MySpectra::UncontPi0sEFTvsTrueE_ND, MySpectra::UncontPi0sEFvsTrueE_FD, MySpectra::UncontPi0sEFvsTrueE_ND, MySpectra::UncontPionEF_FD, MySpectra::UncontPionEF_ND, MySpectra::UncontPionEFT_FD, MySpectra::UncontPionEFT_ND, MySpectra::UncontPionEFTvsTrueE_FD, MySpectra::UncontPionEFTvsTrueE_ND, MySpectra::UncontPionEFvsTrueE_FD, MySpectra::UncontPionEFvsTrueE_ND, MySpectra::UncontProtEF_FD, MySpectra::UncontProtEF_ND, MySpectra::UncontProtEFT_FD, MySpectra::UncontProtEFT_ND, MySpectra::UncontProtEFTvsTrueE_FD, MySpectra::UncontProtEFTvsTrueE_ND, MySpectra::UncontProtEFvsTrueE_FD, MySpectra::UncontProtEFvsTrueE_ND, MySpectra::UncontTotEFT_FD, MySpectra::UncontTotEFT_ND, MySpectra::UncontTotEFTvsTrueE_FD, MySpectra::UncontTotEFTvsTrueE_ND, MySpectra::UncontTotEvsX_FD, MySpectra::UncontTotEvsX_ND, MySpectra::UncontTotEvsY_FD, and MySpectra::UncontTotEvsY_ND.

                       {
  // ---- First off, lets set the styles...
  gStyle->SetOptStat(0);
  gROOT->SetStyle("novaStyle");

  // --- Make my iteration of MySpectra.
  MySpectra sp;
  // --- Loop through my input files...
  for (unsigned int HC=0; HC<nHC; ++HC) {
    // --- Open my files.
//    std::string FDName = "final/SpectrumFiles_FD_"+HCNames[HC]+".root"; TFile *FDFile = TFile::Open( FDName.c_str() );
//    std::string NDName = "final/SpectrumFiles_ND_"+HCNames[HC]+".root"; TFile *NDFile = TFile::Open( NDName.c_str() );

    std::string FDName = "SpectrumFiles_FD_"+HCNames[HC]+"_Full_finerbinning.root"; TFile *FDFile = TFile::Open( FDName.c_str() );
    std::string NDName = "SpectrumFiles_ND_"+HCNames[HC]+"_Full_finerbinning.root"; TFile *NDFile = TFile::Open( NDName.c_str() );
    std::vector<TFile*> files = {FDFile, NDFile};
    // --- Determine the POT
    double POT = kAna2018FHCPOT;
    if (HCNames[HC] == "rhc") POT = kAna2018RHCPOT;
    // Loop through GENIE interaction types.
    for (unsigned int gen=0; gen<nGENIE; ++gen) {
      // Loop through the quantile cuts.
      for (unsigned int quant=0; quant<nQuantNames; ++quant) {
        // --- Loop through the FD cuts first...
        for(unsigned int cut=0; cut<nFD_CutNames; ++cut) {
          // Figure out what the end of my spectra names is.
          std::string Comb = "_"+GENIEStr[gen]+"_"+FD_CutNames[cut]+"_"+QuantNames[quant];
          // Load my spectra
          // multiplicity
          sp.NNeut_FD[gen][cut][quant][HC] = SpecToHist(FDFile, "NumOfNeut"+Comb, POT, ";Num. Neutrons;Events");
          sp.NProt_FD[gen][cut][quant][HC] = SpecToHist(FDFile, "NumOfProt"+Comb, POT, ";Num. Protons;Events" );
          sp.NPion_FD[gen][cut][quant][HC] = SpecToHist(FDFile, "NumOfPion"+Comb, POT, ";Num. Pions;Events"   );
          sp.NPi0s_FD[gen][cut][quant][HC] = SpecToHist(FDFile, "NumOfPi0s"+Comb, POT, ";Num. Pi0s;Events"    );

          // energy
          sp.NeutrinoEnergy_FD[gen][cut][quant][HC] = SpecToHist(FDFile, "NeutrinoE"+Comb, POT, ";True neutrino E (GeV);Events"   );
          sp.MuonEnergy_FD    [gen][cut][quant][HC] = SpecToHist(FDFile, "MuonE"    +Comb, POT, ";Total E in muons (GeV);Events"    );
          sp.NeutronEnergy_FD [gen][cut][quant][HC] = SpecToHist(FDFile, "NeutronE" +Comb, POT, ";Total E in neutrons (GeV);Events" );
          sp.ProtonEnergy_FD  [gen][cut][quant][HC] = SpecToHist(FDFile, "ProtonE"  +Comb, POT, ";Total E in protons (GeV);Events"  );
          sp.PionEnergy_FD    [gen][cut][quant][HC] = SpecToHist(FDFile, "PionE"    +Comb, POT, ";Total E in pions (GeV);Events"    );
          sp.Pi0Energy_FD     [gen][cut][quant][HC] = SpecToHist(FDFile, "Pi0E"     +Comb, POT, ";Total E in pi0s (GeV);Events"     );
          sp.PhotonEnergy_FD  [gen][cut][quant][HC] = SpecToHist(FDFile, "PhotonE"  +Comb, POT, ";Total E in photons (GeV);Events"  );
          sp.KaonEnergy_FD    [gen][cut][quant][HC] = SpecToHist(FDFile, "KaonE"    +Comb, POT, ";Total E in kaons (GeV);Events"    );
          sp.ElectronEnergy_FD[gen][cut][quant][HC] = SpecToHist(FDFile, "ElectronE"+Comb, POT, ";Total E in electrons (GeV);Events");
          sp.MiscEnergy_FD    [gen][cut][quant][HC] = SpecToHist(FDFile, "MiscE"    +Comb, POT, ";Total E in misc (GeV);Events"     );

          // uncontained energy
          sp.AllUncontain_FD  [gen][cut][quant][HC] = SpecToHist(FDFile, "UncontEn" +Comb, POT, ";Total uncontained E (GeV);Events"     );
          sp.MuonUncontain_FD [gen][cut][quant][HC] = SpecToHist(FDFile, "UncontMu" +Comb, POT, ";Uncont. E from muons (GeV);Events"    );
          sp.NeutUncontain_FD [gen][cut][quant][HC] = SpecToHist(FDFile, "UncontNeu"+Comb, POT, ";Uncont. E from neutrons (GeV);Events" );
          sp.ProtUncontain_FD [gen][cut][quant][HC] = SpecToHist(FDFile, "UncontPro"+Comb, POT, ";Uncont. E from protons (GeV);Events"  );
          sp.PionUncontain_FD [gen][cut][quant][HC] = SpecToHist(FDFile, "UncontPi" +Comb, POT, ";Uncont. E from pions (GeV);Events"    );
          sp.Pi0sUncontain_FD [gen][cut][quant][HC] = SpecToHist(FDFile, "UncontPi0"+Comb, POT, ";Uncont. E from pi0s (GeV);Events"     );
          sp.PhotUncontain_FD [gen][cut][quant][HC] = SpecToHist(FDFile, "UncontPho"+Comb, POT, ";Uncont. E from photons (GeV);Events"  );
          sp.KaonUncontain_FD [gen][cut][quant][HC] = SpecToHist(FDFile, "UncontKa" +Comb, POT, ";Uncont. E from kaons (GeV);Events"    );
          sp.ElecUncontain_FD [gen][cut][quant][HC] = SpecToHist(FDFile, "UncontEle"+Comb, POT, ";Uncont. E from electrons (GeV);Events");
          sp.MiscUncontain_FD [gen][cut][quant][HC] = SpecToHist(FDFile, "UncontMis"+Comb, POT, ";Uncont. E from misc. (GeV);Events"    );

          // energy fraction
          sp.MuonEFrac_FD    [gen][cut][quant][HC] = SpecToHist(FDFile, "MuonEFrac"    +Comb, POT, ";Muon E/true #nu E; Events"    );
          sp.NeutronEFrac_FD [gen][cut][quant][HC] = SpecToHist(FDFile, "NeutronEFrac" +Comb, POT, ";Neutron E/true #nu E; Events" );
          sp.ProtonEFrac_FD  [gen][cut][quant][HC] = SpecToHist(FDFile, "ProtonEFrac"  +Comb, POT, ";Proton E/true #nu E; Events"  );
          sp.PionEFrac_FD    [gen][cut][quant][HC] = SpecToHist(FDFile, "PionEFrac"    +Comb, POT, ";Pion E/true #nu E; Events"    );
          sp.Pi0EFrac_FD     [gen][cut][quant][HC] = SpecToHist(FDFile, "Pi0EFrac"     +Comb, POT, ";Pi0 E/true #nu E; Events"     );
          sp.PhotonEFrac_FD  [gen][cut][quant][HC] = SpecToHist(FDFile, "PhotonEFrac"  +Comb, POT, ";Photon E/true #nu E; Events"  );
          sp.KaonEFrac_FD    [gen][cut][quant][HC] = SpecToHist(FDFile, "KaonEFrac"    +Comb, POT, ";Kaon E/true #nu E; Events"    ); 
          sp.ElectronEFrac_FD[gen][cut][quant][HC] = SpecToHist(FDFile, "ElectronEFrac"+Comb, POT, ";Electon E/true #nu E; Events" );

          // uncont. E / tot. uncont. E
          sp.UncontMuonEF_FD[gen][cut][quant][HC] = SpecToHist(FDFile, "UncontMuonEF"+Comb, POT, ";Uncont. muon E/total uncont. E;Events"    );
          sp.UncontNeutEF_FD[gen][cut][quant][HC] = SpecToHist(FDFile, "UncontNeutEF"+Comb, POT, ";Uncont. neutron E/total uncont. E;Events" );
          sp.UncontProtEF_FD[gen][cut][quant][HC] = SpecToHist(FDFile, "UncontProtEF"+Comb, POT, ";Uncont. proton E/total uncont. E;Events"  );
          sp.UncontPionEF_FD[gen][cut][quant][HC] = SpecToHist(FDFile, "UncontPionEF"+Comb, POT, ";Uncont. pion E/total uncont. E;Events"    );
          sp.UncontPi0sEF_FD[gen][cut][quant][HC] = SpecToHist(FDFile, "UncontPi0sEF"+Comb, POT, ";Uncont. pi0 E/total uncont. E;Events"     );
          sp.UncontPhotEF_FD[gen][cut][quant][HC] = SpecToHist(FDFile, "UncontPhotEF"+Comb, POT, ";Uncont. photon E/total uncont. E;Events"  );
          sp.UncontKaonEF_FD[gen][cut][quant][HC] = SpecToHist(FDFile, "UncontKaonEF"+Comb, POT, ";Uncont. kaon E/total uncont. E;Events"    );
          sp.UncontElecEF_FD[gen][cut][quant][HC] = SpecToHist(FDFile, "UncontElecEF"+Comb, POT, ";Uncont. electron E/total uncont. E;Events");

          // uncont. E / true nu E
          sp.UncontTotEFT_FD[gen][cut][quant][HC]  = SpecToHist(FDFile, "UncontTotEFT" +Comb, POT, ";Total uncont. E/true #nu E;Events"   );
          sp.UncontMuonEFT_FD[gen][cut][quant][HC] = SpecToHist(FDFile, "UncontMuonEFT"+Comb, POT, ";Uncont. muon E/true #nu E;Events"    );
          sp.UncontNeutEFT_FD[gen][cut][quant][HC] = SpecToHist(FDFile, "UncontNeutEFT"+Comb, POT, ";Uncont. neutron E/true #nu E;Events" );
          sp.UncontProtEFT_FD[gen][cut][quant][HC] = SpecToHist(FDFile, "UncontProtEFT"+Comb, POT, ";Uncont. proton E/true #nu E;Events"  );
          sp.UncontPionEFT_FD[gen][cut][quant][HC] = SpecToHist(FDFile, "UncontPionEFT"+Comb, POT, ";Uncont. pion E/true #nu E;Events"    );
          sp.UncontPi0sEFT_FD[gen][cut][quant][HC] = SpecToHist(FDFile, "UncontPi0sEFT"+Comb, POT, ";Uncont. pi0 E/true #nu E;Events"     );
          sp.UncontPhotEFT_FD[gen][cut][quant][HC] = SpecToHist(FDFile, "UncontPhotEFT"+Comb, POT, ";Uncont. photon E/true #nu E;Events"  );
          sp.UncontKaonEFT_FD[gen][cut][quant][HC] = SpecToHist(FDFile, "UncontKaonEFT"+Comb, POT, ";Uncont. kaon E/true #nu E;Events"    );
          sp.UncontElecEFT_FD[gen][cut][quant][HC] = SpecToHist(FDFile, "UncontElecEFT"+Comb, POT, ";Uncont. electron E/true #nu E;Events");

          // uncont. E vs. vertex x, y
          sp.UncontTotEvsX_FD[gen][cut][quant][HC] = Spec2DtoHist(FDFile, "UncontTotEvsX"+Comb, POT, ";x (cm); Uncont. E (GeV)");
          sp.UncontTotEvsY_FD[gen][cut][quant][HC] = Spec2DtoHist(FDFile, "UncontTotEvsY"+Comb, POT, ";y (cm); Uncont. E (GeV)");

          // energy fraction vs. true nu E
          sp.MuonEFracvsTrueE_FD    [gen][cut][quant][HC] = 
                                  Spec2DtoHist(FDFile, "MuonEFracvsTrueE"    +Comb, POT, ";True #nu E (GeV);Muon E/true #nu E"    );
          sp.NeutronEFracvsTrueE_FD [gen][cut][quant][HC] = 
                                  Spec2DtoHist(FDFile, "NeutronEFracvsTrueE" +Comb, POT, ";True #nu E (GeV);Neutron E/true #nu E" );
          sp.ProtonEFracvsTrueE_FD  [gen][cut][quant][HC] = 
                                  Spec2DtoHist(FDFile, "ProtonEFracvsTrueE"  +Comb, POT, ";True #nu E (GeV);Proton E/true #nu E"  );
          sp.PionEFracvsTrueE_FD    [gen][cut][quant][HC] = 
                                  Spec2DtoHist(FDFile, "PionEFracvsTrueE"    +Comb, POT, ";True #nu E (GeV);Pion E/true #nu E"    );
          sp.Pi0EFracvsTrueE_FD     [gen][cut][quant][HC] = 
                                  Spec2DtoHist(FDFile, "Pi0EFracvsTrueE"     +Comb, POT, ";True #nu E (GeV);Pi0 E/true #nu E"     );
          sp.PhotonEFracvsTrueE_FD  [gen][cut][quant][HC] = 
                                  Spec2DtoHist(FDFile, "PhotonEFracvsTrueE"  +Comb, POT, ";True #nu E (GeV);Photon E/true #nu E"  );
          sp.KaonEFracvsTrueE_FD    [gen][cut][quant][HC] = 
                                  Spec2DtoHist(FDFile, "KaonEFracvsTrueE"    +Comb, POT, ";True #nu E (GeV);Kaon E/true #nu E"    );
          sp.ElectronEFracvsTrueE_FD[gen][cut][quant][HC] = 
                                  Spec2DtoHist(FDFile, "ElectronEFracvsTrueE"+Comb, POT, ";True #nu E (GeV);Electron E/true #nu E");

          // uncont. E / tot. uncont. E vs. true nu E
          sp.UncontMuonEFvsTrueE_FD[gen][cut][quant][HC] = 
                                 Spec2DtoHist(FDFile, "UncontMuonEFvsTrueE"+Comb, POT, ";True #nu E (GeV);Muon uncont. E/total uncont. E"   );
          sp.UncontNeutEFvsTrueE_FD[gen][cut][quant][HC] = 
                                 Spec2DtoHist(FDFile, "UncontNeutEFvsTrueE"+Comb, POT, ";True #nu E (GeV);Neutron uncont. E/total uncont. E");
          sp.UncontProtEFvsTrueE_FD[gen][cut][quant][HC] = 
                                 Spec2DtoHist(FDFile, "UncontProtEFvsTrueE"+Comb, POT, ";True #nu E (GeV);Proton uncont. E/total uncont. E" );
          sp.UncontPionEFvsTrueE_FD[gen][cut][quant][HC] = 
                                 Spec2DtoHist(FDFile, "UncontPionEFvsTrueE"+Comb, POT, ";True #nu E (GeV);Pion uncont. E/total uncont. E"   );
          sp.UncontPi0sEFvsTrueE_FD[gen][cut][quant][HC] = 
                                 Spec2DtoHist(FDFile, "UncontPi0sEFvsTrueE"+Comb, POT, ";True #nu E (GeV);Pi0 uncont. E/total uncont. E"    ); 
          sp.UncontPhotEFvsTrueE_FD[gen][cut][quant][HC] = 
                                 Spec2DtoHist(FDFile, "UncontPhotEFvsTrueE"+Comb, POT, ";True #nu E (GeV);Photon uncont. E/total uncont. E" );
          sp.UncontKaonEFvsTrueE_FD[gen][cut][quant][HC] = 
                                Spec2DtoHist(FDFile, "UncontKaonEFvsTrueE"+Comb, POT, ";True #nu E (GeV);Kaon uncont. E/total uncont. E"    );
          sp.UncontElecEFvsTrueE_FD[gen][cut][quant][HC] = 
                                Spec2DtoHist(FDFile, "UncontElecEFvsTrueE"+Comb, POT, ";True #nu E (GeV);Electron uncont. E/total uncont. E");

          // uncont. E / true nu E vs. true nu E
          sp.UncontTotEFTvsTrueE_FD [gen][cut][quant][HC] = 
                                Spec2DtoHist(FDFile, "UncontTotEFTvsTrueE"+Comb, POT, ";True #nu E (GeV);Total uncont. E/true #nu E"    );
          sp.UncontMuonEFTvsTrueE_FD[gen][cut][quant][HC] = 
                                 Spec2DtoHist(FDFile, "UncontMuonEFTvsTrueE"+Comb, POT, ";True #nu E (GeV);Muon uncont. E/true #nu E"    );
          sp.UncontNeutEFTvsTrueE_FD[gen][cut][quant][HC] = 
                                 Spec2DtoHist(FDFile, "UncontNeutEFTvsTrueE"+Comb, POT, ";True #nu E (GeV);Neutron uncont. E/true #nu E" );
          sp.UncontProtEFTvsTrueE_FD[gen][cut][quant][HC] = 
                                 Spec2DtoHist(FDFile, "UncontProtEFTvsTrueE"+Comb, POT, ";True #nu E (GeV);Proton uncont. E/true #nu E"  );
          sp.UncontPionEFTvsTrueE_FD[gen][cut][quant][HC] = 
                                 Spec2DtoHist(FDFile, "UncontPionEFTvsTrueE"+Comb, POT, ";True #nu E (GeV);Pion uncont. E/true #nu E"    );
          sp.UncontPi0sEFTvsTrueE_FD[gen][cut][quant][HC] = 
                                 Spec2DtoHist(FDFile, "UncontPi0sEFTvsTrueE"+Comb, POT, ";True #nu E (GeV);Pi0 uncont. E/true #nu E"     );
          sp.UncontPhotEFTvsTrueE_FD[gen][cut][quant][HC] = 
                                 Spec2DtoHist(FDFile, "UncontPhotEFTvsTrueE"+Comb, POT, ";True #nu E (GeV);Photon uncont. E/true #nu E"  );
          sp.UncontKaonEFTvsTrueE_FD[gen][cut][quant][HC] = 
                                 Spec2DtoHist(FDFile, "UncontKaonEFTvsTrueE"+Comb, POT, ";True #nu E (GeV);Kaon uncont. E/true #nu E"    );
          sp.UncontElecEFTvsTrueE_FD[gen][cut][quant][HC] = 
                                 Spec2DtoHist(FDFile, "UncontElecEFTvsTrueE"+Comb, POT, ";True #nu E (GeV);Electron uncont. E/true #nu E");

        } // nFD_CutNames

        // --- Now, loop through my ND cuts...
        for(unsigned int cut=0; cut<nND_CutNames; ++cut) {
          // Figure out what the end of my spectra names is.
          std::string Comb = "_"+GENIEStr[gen]+"_"+ND_CutNames[cut]+"_"+QuantNames[quant];
          // Load my spectra
          // Multiplicity
          sp.NNeut_ND[gen][cut][quant][HC] = SpecToHist(NDFile, "NumOfNeut"+Comb, POT, ";Num. Neutrons;Events");
          sp.NProt_ND[gen][cut][quant][HC] = SpecToHist(NDFile, "NumOfProt"+Comb, POT, ";Num. Protons;Events" );
          sp.NPion_ND[gen][cut][quant][HC] = SpecToHist(NDFile, "NumOfPion"+Comb, POT, ";Num. Pions;Events"   );
          sp.NPi0s_ND[gen][cut][quant][HC] = SpecToHist(NDFile, "NumOfPi0s"+Comb, POT, ";Num. Pi0s;Events"    );

          // energy
          sp.NeutrinoEnergy_ND[gen][cut][quant][HC] = SpecToHist(NDFile, "NeutrinoE"+Comb, POT, ";Recon. neutrino E (GeV);Events"   );
          sp.MuonEnergy_ND    [gen][cut][quant][HC] = SpecToHist(NDFile, "MuonE"    +Comb, POT, ";Total E in muons (GeV);Events"    );
          sp.NeutronEnergy_ND [gen][cut][quant][HC] = SpecToHist(NDFile, "NeutronE" +Comb, POT, ";Total E in neutrons (GeV);Events" );
          sp.ProtonEnergy_ND  [gen][cut][quant][HC] = SpecToHist(NDFile, "ProtonE"  +Comb, POT, ";Total E in protons (GeV);Events"  );
          sp.PionEnergy_ND    [gen][cut][quant][HC] = SpecToHist(NDFile, "PionE"    +Comb, POT, ";Total E in pions (GeV);Events"    );
          sp.Pi0Energy_ND     [gen][cut][quant][HC] = SpecToHist(NDFile, "Pi0E"     +Comb, POT, ";Total E in pi0s (GeV);Events"     );
          sp.PhotonEnergy_ND  [gen][cut][quant][HC] = SpecToHist(NDFile, "PhotonE"  +Comb, POT, ";Total E in photons (GeV);Events"  );
          sp.KaonEnergy_ND    [gen][cut][quant][HC] = SpecToHist(NDFile, "KaonE"    +Comb, POT, ";Total E in kaons (GeV);Events"    );
          sp.ElectronEnergy_ND[gen][cut][quant][HC] = SpecToHist(NDFile, "ElectronE"+Comb, POT, ";Total E in electrons (GeV);Events");
          sp.MiscEnergy_ND    [gen][cut][quant][HC] = SpecToHist(NDFile, "MiscE"    +Comb, POT, ";Total E in misc (GeV);Events"     );

          // uncontained energy
          sp.AllUncontain_ND  [gen][cut][quant][HC] = SpecToHist(NDFile, "UncontEn" +Comb, POT, ";Total uncontained E (GeV);Events"     );
          sp.MuonUncontain_ND [gen][cut][quant][HC] = SpecToHist(NDFile, "UncontMu" +Comb, POT, ";Uncont. E from muons (GeV);Events"    );
          sp.NeutUncontain_ND [gen][cut][quant][HC] = SpecToHist(NDFile, "UncontNeu"+Comb, POT, ";Uncont. E from neutrons (GeV);Events" );
          sp.ProtUncontain_ND [gen][cut][quant][HC] = SpecToHist(NDFile, "UncontPro"+Comb, POT, ";Uncont. E from protons (GeV);Events"  );
          sp.PionUncontain_ND [gen][cut][quant][HC] = SpecToHist(NDFile, "UncontPi" +Comb, POT, ";Uncont. E from pions (GeV);Events"    );
          sp.Pi0sUncontain_ND [gen][cut][quant][HC] = SpecToHist(NDFile, "UncontPi0"+Comb, POT, ";Uncont. E from pi0s (GeV);Events"     );
          sp.PhotUncontain_ND [gen][cut][quant][HC] = SpecToHist(NDFile, "UncontPho"+Comb, POT, ";Uncont. E from photons (GeV);Events"  );
          sp.KaonUncontain_ND [gen][cut][quant][HC] = SpecToHist(NDFile, "UncontKa" +Comb, POT, ";Uncont. E from kaons (GeV);Events"    );
          sp.ElecUncontain_ND [gen][cut][quant][HC] = SpecToHist(NDFile, "UncontEle"+Comb, POT, ";Uncont. E from electrons (GeV);Events");
          sp.MiscUncontain_ND [gen][cut][quant][HC] = SpecToHist(NDFile, "UncontMis"+Comb, POT, ";Uncont. E from misc. (GeV);Events"    );

          // energy fraction
          sp.MuonEFrac_ND    [gen][cut][quant][HC] = SpecToHist(NDFile, "MuonEFrac"    +Comb, POT, ";Muon E/true #nu E; Events"    );
          sp.NeutronEFrac_ND [gen][cut][quant][HC] = SpecToHist(NDFile, "NeutronEFrac" +Comb, POT, ";Neutron E/true #nu E; Events" );
          sp.ProtonEFrac_ND  [gen][cut][quant][HC] = SpecToHist(NDFile, "ProtonEFrac"  +Comb, POT, ";Proton E/true #nu E; Events"  );
          sp.PionEFrac_ND    [gen][cut][quant][HC] = SpecToHist(NDFile, "PionEFrac"    +Comb, POT, ";Pion E/true #nu E; Events"    );
          sp.Pi0EFrac_ND     [gen][cut][quant][HC] = SpecToHist(NDFile, "Pi0EFrac"     +Comb, POT, ";Pi0 E/true #nu E; Events"     );
          sp.PhotonEFrac_ND  [gen][cut][quant][HC] = SpecToHist(NDFile, "PhotonEFrac"  +Comb, POT, ";Photon E/true #nu E; Events"  );
          sp.KaonEFrac_ND    [gen][cut][quant][HC] = SpecToHist(NDFile, "KaonEFrac"    +Comb, POT, ";Kaon E/true #nu E; Events"    ); 
          sp.ElectronEFrac_ND[gen][cut][quant][HC] = SpecToHist(NDFile, "ElectronEFrac"+Comb, POT, ";Electon E/true #nu E; Events" );

          // uncont. E / tot. uncont. E
          sp.UncontMuonEF_ND[gen][cut][quant][HC] = SpecToHist(NDFile, "UncontMuonEF"+Comb, POT, ";Uncont. muon E/total uncont. E;Events"    );
          sp.UncontNeutEF_ND[gen][cut][quant][HC] = SpecToHist(NDFile, "UncontNeutEF"+Comb, POT, ";Uncont. neutron E/total uncont. E;Events" );
          sp.UncontProtEF_ND[gen][cut][quant][HC] = SpecToHist(NDFile, "UncontProtEF"+Comb, POT, ";Uncont. proton E/total uncont. E;Events"  );
          sp.UncontPionEF_ND[gen][cut][quant][HC] = SpecToHist(NDFile, "UncontPionEF"+Comb, POT, ";Uncont. pion E/total uncont. E;Events"    );
          sp.UncontPi0sEF_ND[gen][cut][quant][HC] = SpecToHist(NDFile, "UncontPi0sEF"+Comb, POT, ";Uncont. pi0 E/total uncont. E;Events"     );
          sp.UncontPhotEF_ND[gen][cut][quant][HC] = SpecToHist(NDFile, "UncontPhotEF"+Comb, POT, ";Uncont. photon E/total uncont. E;Events"  );
          sp.UncontKaonEF_ND[gen][cut][quant][HC] = SpecToHist(NDFile, "UncontKaonEF"+Comb, POT, ";Uncont. kaon E/total uncont. E;Events"    );
          sp.UncontElecEF_ND[gen][cut][quant][HC] = SpecToHist(NDFile, "UncontElecEF"+Comb, POT, ";Uncont. electron E/total uncont. E;Events");

          // uncont. E / true nu E
          sp.UncontTotEFT_ND[gen][cut][quant][HC]  = SpecToHist(NDFile, "UncontTotEFT" +Comb, POT, ";Total uncont. E/true #nu E;Events"   );
          sp.UncontMuonEFT_ND[gen][cut][quant][HC] = SpecToHist(NDFile, "UncontMuonEFT"+Comb, POT, ";Uncont. muon E/true #nu E;Events"    );
          sp.UncontNeutEFT_ND[gen][cut][quant][HC] = SpecToHist(NDFile, "UncontNeutEFT"+Comb, POT, ";Uncont. neutron E/true #nu E;Events" );
          sp.UncontProtEFT_ND[gen][cut][quant][HC] = SpecToHist(NDFile, "UncontProtEFT"+Comb, POT, ";Uncont. proton E/true #nu E;Events"  );
          sp.UncontPionEFT_ND[gen][cut][quant][HC] = SpecToHist(NDFile, "UncontPionEFT"+Comb, POT, ";Uncont. pion E/true #nu E;Events"    );
          sp.UncontPi0sEFT_ND[gen][cut][quant][HC] = SpecToHist(NDFile, "UncontPi0sEFT"+Comb, POT, ";Uncont. pi0 E/true #nu E;Events"     );
          sp.UncontPhotEFT_ND[gen][cut][quant][HC] = SpecToHist(NDFile, "UncontPhotEFT"+Comb, POT, ";Uncont. photon E/true #nu E;Events"  );
          sp.UncontKaonEFT_ND[gen][cut][quant][HC] = SpecToHist(NDFile, "UncontKaonEFT"+Comb, POT, ";Uncont. kaon E/true #nu E;Events"    );
          sp.UncontElecEFT_ND[gen][cut][quant][HC] = SpecToHist(NDFile, "UncontElecEFT"+Comb, POT, ";Uncont. electron E/true #nu E;Events");

          // uncont. E vs. vertex x, y
          sp.UncontTotEvsX_ND[gen][cut][quant][HC] = Spec2DtoHist(NDFile, "UncontTotEvsX"+Comb, POT, ";x (cm); Uncont. E (GeV)");
          sp.UncontTotEvsY_ND[gen][cut][quant][HC] = Spec2DtoHist(NDFile, "UncontTotEvsY"+Comb, POT, ";y (cm); Uncont. E (GeV)");

          // energy fraction vs. true nu E
          sp.MuonEFracvsTrueE_ND    [gen][cut][quant][HC] = 
                                  Spec2DtoHist(NDFile, "MuonEFracvsTrueE"    +Comb, POT, ";True #nu E (GeV);Muon E/true #nu E"    );
          sp.NeutronEFracvsTrueE_ND [gen][cut][quant][HC] = 
                                  Spec2DtoHist(NDFile, "NeutronEFracvsTrueE" +Comb, POT, ";True #nu E (GeV);Neutron E/true #nu E" );
          sp.ProtonEFracvsTrueE_ND  [gen][cut][quant][HC] = 
                                  Spec2DtoHist(NDFile, "ProtonEFracvsTrueE"  +Comb, POT, ";True #nu E (GeV);Proton E/true #nu E"  );
          sp.PionEFracvsTrueE_ND    [gen][cut][quant][HC] = 
                                  Spec2DtoHist(NDFile, "PionEFracvsTrueE"    +Comb, POT, ";True #nu E (GeV);Pion E/true #nu E"    );
          sp.Pi0EFracvsTrueE_ND     [gen][cut][quant][HC] = 
                                  Spec2DtoHist(NDFile, "Pi0EFracvsTrueE"     +Comb, POT, ";True #nu E (GeV);Pi0 E/true #nu E"     );
          sp.PhotonEFracvsTrueE_ND  [gen][cut][quant][HC] = 
                                  Spec2DtoHist(NDFile, "PhotonEFracvsTrueE"  +Comb, POT, ";True #nu E (GeV);Photon E/true #nu E"  );
          sp.KaonEFracvsTrueE_ND    [gen][cut][quant][HC] = 
                                  Spec2DtoHist(NDFile, "KaonEFracvsTrueE"    +Comb, POT, ";True #nu E (GeV);Kaon E/true #nu E"    );
          sp.ElectronEFracvsTrueE_ND[gen][cut][quant][HC] = 
                                  Spec2DtoHist(NDFile, "ElectronEFracvsTrueE"+Comb, POT, ";True #nu E (GeV);Electron E/true #nu E");

          // uncont. E / tot. uncont. E vs. true nu E
          sp.UncontMuonEFvsTrueE_ND[gen][cut][quant][HC] = 
                                 Spec2DtoHist(NDFile, "UncontMuonEFvsTrueE"+Comb, POT, ";True #nu E (GeV);Muon uncont. E/total uncont. E"   );
          sp.UncontNeutEFvsTrueE_ND[gen][cut][quant][HC] = 
                                 Spec2DtoHist(NDFile, "UncontNeutEFvsTrueE"+Comb, POT, ";True #nu E (GeV);Neutron uncont. E/total uncont. E");
          sp.UncontProtEFvsTrueE_ND[gen][cut][quant][HC] = 
                                 Spec2DtoHist(NDFile, "UncontProtEFvsTrueE"+Comb, POT, ";True #nu E (GeV);Proton uncont. E/total uncont. E" );
          sp.UncontPionEFvsTrueE_ND[gen][cut][quant][HC] = 
                                 Spec2DtoHist(NDFile, "UncontPionEFvsTrueE"+Comb, POT, ";True #nu E (GeV);Pion uncont. E/total uncont. E"   );
          sp.UncontPi0sEFvsTrueE_ND[gen][cut][quant][HC] = 
                                 Spec2DtoHist(NDFile, "UncontPi0sEFvsTrueE"+Comb, POT, ";True #nu E (GeV);Pi0 uncont. E/total uncont. E"    ); 
          sp.UncontPhotEFvsTrueE_ND[gen][cut][quant][HC] = 
                                 Spec2DtoHist(NDFile, "UncontPhotEFvsTrueE"+Comb, POT, ";True #nu E (GeV);Photon uncont. E/total uncont. E" );
          sp.UncontKaonEFvsTrueE_ND[gen][cut][quant][HC] = 
                                Spec2DtoHist(NDFile, "UncontKaonEFvsTrueE"+Comb, POT, ";True #nu E (GeV);Kaon uncont. E/total uncont. E"    );
          sp.UncontElecEFvsTrueE_ND[gen][cut][quant][HC] = 
                                Spec2DtoHist(NDFile, "UncontElecEFvsTrueE"+Comb, POT, ";True #nu E (GeV);Electron uncont. E/total uncont. E");

          // uncont. E / true nu E vs. true nu E 
          sp.UncontTotEFTvsTrueE_ND [gen][cut][quant][HC] = 
                                Spec2DtoHist(NDFile, "UncontTotEFTvsTrueE"+Comb, POT, ";True #nu E (GeV);Total uncont. E/true #nu E"    );
          sp.UncontMuonEFTvsTrueE_ND[gen][cut][quant][HC] = 
                                 Spec2DtoHist(NDFile, "UncontMuonEFTvsTrueE"+Comb, POT, ";True #nu E (GeV);Muon uncont. E/true #nu E"    );
          sp.UncontNeutEFTvsTrueE_ND[gen][cut][quant][HC] = 
                                 Spec2DtoHist(NDFile, "UncontNeutEFTvsTrueE"+Comb, POT, ";True #nu E (GeV);Neutron uncont. E/true #nu E" );
          sp.UncontProtEFTvsTrueE_ND[gen][cut][quant][HC] = 
                                 Spec2DtoHist(NDFile, "UncontProtEFTvsTrueE"+Comb, POT, ";True #nu E (GeV);Proton uncont. E/true #nu E"  );
          sp.UncontPionEFTvsTrueE_ND[gen][cut][quant][HC] = 
                                 Spec2DtoHist(NDFile, "UncontPionEFTvsTrueE"+Comb, POT, ";True #nu E (GeV);Pion uncont. E/true #nu E"    );
          sp.UncontPi0sEFTvsTrueE_ND[gen][cut][quant][HC] = 
                                 Spec2DtoHist(NDFile, "UncontPi0sEFTvsTrueE"+Comb, POT, ";True #nu E (GeV);Pi0 uncont. E/true #nu E"     );
          sp.UncontPhotEFTvsTrueE_ND[gen][cut][quant][HC] = 
                                 Spec2DtoHist(NDFile, "UncontPhotEFTvsTrueE"+Comb, POT, ";True #nu E (GeV);Photon uncont. E/true #nu E"  );
          sp.UncontKaonEFTvsTrueE_ND[gen][cut][quant][HC] = 
                                 Spec2DtoHist(NDFile, "UncontKaonEFTvsTrueE"+Comb, POT, ";True #nu E (GeV);Kaon uncont. E/true #nu E"    );
          sp.UncontElecEFTvsTrueE_ND[gen][cut][quant][HC] = 
                                 Spec2DtoHist(NDFile, "UncontElecEFTvsTrueE"+Comb, POT, ";True #nu E (GeV);Electron uncont. E/true #nu E");
        } // nND_CutNames
      } // nQuantNames
    } // nGENIE
  } // files

  // ...................................................... Make my combined plots ......................................................
  // *************** Want to plot {HC + GENIE + Quantile} as a function of Cut Tiers ***************
  for (unsigned int HC=0; HC<nHC; ++HC){
    for (unsigned int gen=0; gen<nGENIE; ++gen) {
      for (unsigned int quant=0; quant<nQuantNames; ++quant) {
        // --- Now I can make the combination string...
        std::string Comb = "_"+GENIEStr[gen]+"_"+QuantNames[quant]+"_"+HCNames[HC];
        // --- FD First.
        // energy 
        MakeComparisonPlots( {sp.NeutrinoEnergy_FD[gen][0][quant][HC], sp.NeutrinoEnergy_FD[gen][1][quant][HC], sp.NeutrinoEnergy_FD[gen][2][quant][HC]}, "NeutrinoEnergy_FD"+Comb, FD_CutNames );
        MakeComparisonPlots( {sp.MuonEnergy_FD    [gen][0][quant][HC], sp.MuonEnergy_FD    [gen][1][quant][HC], sp.MuonEnergy_FD    [gen][2][quant][HC]}, "MuonEnergy_FD"    +Comb, FD_CutNames );
        MakeComparisonPlots( {sp.NeutronEnergy_FD [gen][0][quant][HC], sp.NeutronEnergy_FD [gen][1][quant][HC], sp.NeutronEnergy_FD [gen][2][quant][HC]}, "NeutronEnergy_FD" +Comb, FD_CutNames );
        MakeComparisonPlots( {sp.ProtonEnergy_FD  [gen][0][quant][HC], sp.ProtonEnergy_FD  [gen][1][quant][HC], sp.ProtonEnergy_FD  [gen][2][quant][HC]}, "ProtonEnergy_FD"  +Comb, FD_CutNames );
        MakeComparisonPlots( {sp.PionEnergy_FD    [gen][0][quant][HC], sp.PionEnergy_FD    [gen][1][quant][HC], sp.PionEnergy_FD    [gen][2][quant][HC]}, "PionEnergy_FD"  +Comb, FD_CutNames );
        MakeComparisonPlots( {sp.Pi0Energy_FD     [gen][0][quant][HC], sp.Pi0Energy_FD     [gen][1][quant][HC], sp.Pi0Energy_FD     [gen][2][quant][HC]}, "Pi0Energy_FD"  +Comb, FD_CutNames );
        MakeComparisonPlots( {sp.PhotonEnergy_FD  [gen][0][quant][HC], sp.PhotonEnergy_FD  [gen][1][quant][HC], sp.PhotonEnergy_FD  [gen][2][quant][HC]}, "PhotonEnergy_FD"  +Comb, FD_CutNames );
        MakeComparisonPlots( {sp.KaonEnergy_FD    [gen][0][quant][HC], sp.KaonEnergy_FD    [gen][1][quant][HC], sp.KaonEnergy_FD    [gen][2][quant][HC]}, "KaonEnergy_FD"  +Comb, FD_CutNames );
        MakeComparisonPlots( {sp.ElectronEnergy_FD[gen][0][quant][HC], sp.ElectronEnergy_FD[gen][1][quant][HC], sp.ElectronEnergy_FD[gen][2][quant][HC]}, "ElectronEnergy_FD"  +Comb, FD_CutNames );
        MakeComparisonPlots( {sp.MiscEnergy_FD    [gen][0][quant][HC], sp.MiscEnergy_FD    [gen][1][quant][HC], sp.MiscEnergy_FD    [gen][2][quant][HC]}, "MiscEnergy_FD"  +Comb, FD_CutNames );

        // uncontained energy 
        MakeComparisonPlots( {sp.AllUncontain_FD  [gen][0][quant][HC], sp.AllUncontain_FD  [gen][1][quant][HC], sp.AllUncontain_FD  [gen][2][quant][HC]}, "AllUncontain_FD"  +Comb, FD_CutNames );
        MakeComparisonPlots( {sp.MuonUncontain_FD [gen][0][quant][HC], sp.MuonUncontain_FD [gen][1][quant][HC], sp.MuonUncontain_FD [gen][2][quant][HC]}, "MuonUncontain_FD" +Comb, FD_CutNames );
        MakeComparisonPlots( {sp.NeutUncontain_FD [gen][0][quant][HC], sp.NeutUncontain_FD [gen][1][quant][HC], sp.NeutUncontain_FD [gen][2][quant][HC]}, "NeutUncontain_FD" +Comb, FD_CutNames );
        MakeComparisonPlots( {sp.ProtUncontain_FD [gen][0][quant][HC], sp.ProtUncontain_FD [gen][1][quant][HC], sp.ProtUncontain_FD [gen][2][quant][HC]}, "ProtUncontain_FD" +Comb, FD_CutNames );
        MakeComparisonPlots( {sp.PionUncontain_FD [gen][0][quant][HC], sp.PionUncontain_FD [gen][1][quant][HC], sp.PionUncontain_FD [gen][2][quant][HC]}, "PionUncontain_FD" +Comb, FD_CutNames );
        MakeComparisonPlots( {sp.Pi0sUncontain_FD [gen][0][quant][HC], sp.Pi0sUncontain_FD [gen][1][quant][HC], sp.Pi0sUncontain_FD [gen][2][quant][HC]}, "Pi0sUncontain_FD" +Comb, FD_CutNames );
        MakeComparisonPlots( {sp.PhotUncontain_FD [gen][0][quant][HC], sp.PhotUncontain_FD [gen][1][quant][HC], sp.PhotUncontain_FD [gen][2][quant][HC]}, "PhotUncontain_FD" +Comb, FD_CutNames );
        MakeComparisonPlots( {sp.KaonUncontain_FD [gen][0][quant][HC], sp.KaonUncontain_FD [gen][1][quant][HC], sp.KaonUncontain_FD [gen][2][quant][HC]}, "KaonUncontain_FD" +Comb, FD_CutNames );
        MakeComparisonPlots( {sp.ElecUncontain_FD [gen][0][quant][HC], sp.ElecUncontain_FD [gen][1][quant][HC], sp.ElecUncontain_FD [gen][2][quant][HC]}, "ElecUncontain_FD" +Comb, FD_CutNames );
        MakeComparisonPlots( {sp.MiscUncontain_FD [gen][0][quant][HC], sp.MiscUncontain_FD [gen][1][quant][HC], sp.MiscUncontain_FD [gen][2][quant][HC]}, "MiscUncontain_FD" +Comb, FD_CutNames );

        // energy fraction
        MakeComparisonPlots( {sp.MuonEFrac_FD    [gen][0][quant][HC], sp.MuonEFrac_FD    [gen][1][quant][HC], sp.MuonEFrac_FD [gen][2][quant][HC]}, "MuonEFrac_FD" +Comb, FD_CutNames      );
        MakeComparisonPlots( {sp.NeutronEFrac_FD [gen][0][quant][HC], sp.NeutronEFrac_FD [gen][1][quant][HC], sp.NeutronEFrac_FD [gen][2][quant][HC]}, "NeutronEFrac_FD" +Comb, FD_CutNames);
        MakeComparisonPlots( {sp.ProtonEFrac_FD  [gen][0][quant][HC], sp.ProtonEFrac_FD  [gen][1][quant][HC], sp.ProtonEFrac_FD  [gen][2][quant][HC]}, "ProtonEFrac_FD" +Comb, FD_CutNames );
        MakeComparisonPlots( {sp.PionEFrac_FD    [gen][0][quant][HC], sp.PionEFrac_FD    [gen][1][quant][HC], sp.PionEFrac_FD    [gen][2][quant][HC]}, "PionEFrac_FD" +Comb, FD_CutNames   );
        MakeComparisonPlots( {sp.Pi0EFrac_FD     [gen][0][quant][HC], sp.Pi0EFrac_FD     [gen][1][quant][HC], sp.Pi0EFrac_FD     [gen][2][quant][HC]}, "Pi0EFrac_FD" +Comb, FD_CutNames    );
        MakeComparisonPlots( {sp.PhotonEFrac_FD  [gen][0][quant][HC], sp.PhotonEFrac_FD  [gen][1][quant][HC], sp.PhotonEFrac_FD  [gen][2][quant][HC]}, "PhotonEFrac_FD" +Comb, FD_CutNames );
        MakeComparisonPlots( {sp.KaonEFrac_FD    [gen][0][quant][HC], sp.KaonEFrac_FD    [gen][1][quant][HC], sp.KaonEFrac_FD    [gen][2][quant][HC]}, "KaonEFrac_FD" +Comb, FD_CutNames   );
        MakeComparisonPlots( {sp.ElectronEFrac_FD[gen][0][quant][HC], sp.ElectronEFrac_FD[gen][1][quant][HC], sp.ElectronEFrac_FD[gen][2][quant][HC]}, "ElectronEFrac_FD"+Comb, FD_CutNames);

        // uncont. E / tot. uncont. E
        MakeComparisonPlots( {sp.UncontMuonEF_FD[gen][0][quant][HC], sp.UncontMuonEF_FD[gen][1][quant][HC], sp.UncontMuonEF_FD[gen][2][quant][HC]}, "UncontMuonEF_FD"+Comb, FD_CutNames);
        MakeComparisonPlots( {sp.UncontNeutEF_FD[gen][0][quant][HC], sp.UncontNeutEF_FD[gen][1][quant][HC], sp.UncontNeutEF_FD[gen][2][quant][HC]}, "UncontNeutEF_FD"+Comb, FD_CutNames);
        MakeComparisonPlots( {sp.UncontProtEF_FD[gen][0][quant][HC], sp.UncontProtEF_FD[gen][1][quant][HC], sp.UncontProtEF_FD[gen][2][quant][HC]}, "UncontProtEF_FD"+Comb, FD_CutNames);
        MakeComparisonPlots( {sp.UncontPionEF_FD[gen][0][quant][HC], sp.UncontPionEF_FD[gen][1][quant][HC], sp.UncontPionEF_FD[gen][2][quant][HC]}, "UncontPionEF_FD"+Comb, FD_CutNames);
        MakeComparisonPlots( {sp.UncontPi0sEF_FD[gen][0][quant][HC], sp.UncontPi0sEF_FD[gen][1][quant][HC], sp.UncontPi0sEF_FD[gen][2][quant][HC]}, "UncontPi0sEF_FD"+Comb, FD_CutNames);
        MakeComparisonPlots( {sp.UncontPhotEF_FD[gen][0][quant][HC], sp.UncontPhotEF_FD[gen][1][quant][HC], sp.UncontPhotEF_FD[gen][2][quant][HC]}, "UncontPhotEF_FD"+Comb, FD_CutNames);
        MakeComparisonPlots( {sp.UncontKaonEF_FD[gen][0][quant][HC], sp.UncontKaonEF_FD[gen][1][quant][HC], sp.UncontKaonEF_FD[gen][2][quant][HC]}, "UncontKaonEF_FD"+Comb, FD_CutNames);
        MakeComparisonPlots( {sp.UncontElecEF_FD[gen][0][quant][HC], sp.UncontElecEF_FD[gen][1][quant][HC], sp.UncontElecEF_FD[gen][2][quant][HC]}, "UncontElecEF_FD"+Comb, FD_CutNames);

        // uncont. E / true nu E
        MakeComparisonPlots( {sp.UncontTotEFT_FD[gen][0][quant][HC], sp.UncontTotEFT_FD[gen][1][quant][HC], sp.UncontTotEFT_FD[gen][2][quant][HC]}, "UncontTotEFT_FD"+Comb, FD_CutNames);
        MakeComparisonPlots( {sp.UncontMuonEFT_FD[gen][0][quant][HC], sp.UncontMuonEFT_FD[gen][1][quant][HC], sp.UncontMuonEFT_FD[gen][2][quant][HC]}, "UncontMuonEFT_FD"+Comb, FD_CutNames);
        MakeComparisonPlots( {sp.UncontNeutEFT_FD[gen][0][quant][HC], sp.UncontNeutEFT_FD[gen][1][quant][HC], sp.UncontNeutEFT_FD[gen][2][quant][HC]}, "UncontNeutEFT_FD"+Comb, FD_CutNames);
        MakeComparisonPlots( {sp.UncontProtEFT_FD[gen][0][quant][HC], sp.UncontProtEFT_FD[gen][1][quant][HC], sp.UncontProtEFT_FD[gen][2][quant][HC]}, "UncontProtEFT_FD"+Comb, FD_CutNames);
        MakeComparisonPlots( {sp.UncontPionEFT_FD[gen][0][quant][HC], sp.UncontPionEFT_FD[gen][1][quant][HC], sp.UncontPionEFT_FD[gen][2][quant][HC]}, "UncontPionEFT_FD"+Comb, FD_CutNames);
        MakeComparisonPlots( {sp.UncontPi0sEFT_FD[gen][0][quant][HC], sp.UncontPi0sEFT_FD[gen][1][quant][HC], sp.UncontPi0sEFT_FD[gen][2][quant][HC]}, "UncontPi0sEFT_FD"+Comb, FD_CutNames);
        MakeComparisonPlots( {sp.UncontPhotEFT_FD[gen][0][quant][HC], sp.UncontPhotEFT_FD[gen][1][quant][HC], sp.UncontPhotEFT_FD[gen][2][quant][HC]}, "UncontPhotEFT_FD"+Comb, FD_CutNames);
        MakeComparisonPlots( {sp.UncontKaonEFT_FD[gen][0][quant][HC], sp.UncontKaonEFT_FD[gen][1][quant][HC], sp.UncontKaonEFT_FD[gen][2][quant][HC]}, "UncontKaonEFT_FD"+Comb, FD_CutNames);
        MakeComparisonPlots( {sp.UncontElecEFT_FD[gen][0][quant][HC], sp.UncontElecEFT_FD[gen][1][quant][HC], sp.UncontElecEFT_FD[gen][2][quant][HC]}, "UncontElecEFT_FD"+Comb, FD_CutNames);

        // --- ND Second.
        // energy 
        MakeComparisonPlots( {sp.NeutrinoEnergy_ND[gen][0][quant][HC], sp.NeutrinoEnergy_ND[gen][1][quant][HC]}, "NeutrinoEnergy_ND"+Comb, ND_CutNames );
        MakeComparisonPlots( {sp.MuonEnergy_ND    [gen][0][quant][HC], sp.MuonEnergy_ND    [gen][1][quant][HC]}, "MuonEnergy_ND"    +Comb, ND_CutNames );
        MakeComparisonPlots( {sp.NeutronEnergy_ND [gen][0][quant][HC], sp.NeutronEnergy_ND [gen][1][quant][HC]}, "NeutronEnergy_ND" +Comb, ND_CutNames );
        MakeComparisonPlots( {sp.ProtonEnergy_ND  [gen][0][quant][HC], sp.ProtonEnergy_ND  [gen][1][quant][HC]}, "ProtonEnergy_ND"  +Comb, ND_CutNames );
        MakeComparisonPlots( {sp.PionEnergy_ND    [gen][0][quant][HC], sp.PionEnergy_ND    [gen][1][quant][HC]}, "PionEnergy_ND"    +Comb, ND_CutNames );
        MakeComparisonPlots( {sp.Pi0Energy_ND     [gen][0][quant][HC], sp.Pi0Energy_ND     [gen][1][quant][HC]}, "Pi0Energy_ND"     +Comb, ND_CutNames );
        MakeComparisonPlots( {sp.PhotonEnergy_ND  [gen][0][quant][HC], sp.PhotonEnergy_ND  [gen][1][quant][HC]}, "PhotonEnergy_ND"  +Comb, ND_CutNames );
        MakeComparisonPlots( {sp.KaonEnergy_ND    [gen][0][quant][HC], sp.KaonEnergy_ND    [gen][1][quant][HC]}, "KaonEnergy_ND"    +Comb, ND_CutNames );
        MakeComparisonPlots( {sp.ElectronEnergy_ND[gen][0][quant][HC], sp.ElectronEnergy_ND[gen][1][quant][HC]}, "ElectronEnergy_ND"+Comb, ND_CutNames );
        MakeComparisonPlots( {sp.MiscEnergy_ND    [gen][0][quant][HC], sp.MiscEnergy_ND    [gen][1][quant][HC]}, "MiscEnergy_ND"    +Comb, ND_CutNames );

        // uncontained energy 
        MakeComparisonPlots( {sp.AllUncontain_ND  [gen][0][quant][HC], sp.AllUncontain_ND  [gen][1][quant][HC]}, "AllUncontain_ND"  +Comb, ND_CutNames );
        MakeComparisonPlots( {sp.MuonUncontain_ND [gen][0][quant][HC], sp.MuonUncontain_ND [gen][1][quant][HC]}, "MuonUncontain_ND" +Comb, ND_CutNames );
        MakeComparisonPlots( {sp.NeutUncontain_ND [gen][0][quant][HC], sp.NeutUncontain_ND [gen][1][quant][HC]}, "NeutUncontain_ND" +Comb, ND_CutNames );
        MakeComparisonPlots( {sp.ProtUncontain_ND [gen][0][quant][HC], sp.ProtUncontain_ND [gen][1][quant][HC]}, "ProtUncontain_ND" +Comb, ND_CutNames );
        MakeComparisonPlots( {sp.PionUncontain_ND [gen][0][quant][HC], sp.PionUncontain_ND [gen][1][quant][HC]}, "PionUncontain_ND" +Comb, ND_CutNames );
        MakeComparisonPlots( {sp.Pi0sUncontain_ND [gen][0][quant][HC], sp.Pi0sUncontain_ND [gen][1][quant][HC]}, "Pi0sUncontain_ND" +Comb, ND_CutNames );
        MakeComparisonPlots( {sp.PhotUncontain_ND [gen][0][quant][HC], sp.PhotUncontain_ND [gen][1][quant][HC]}, "PhotUncontain_ND" +Comb, ND_CutNames );
        MakeComparisonPlots( {sp.KaonUncontain_ND [gen][0][quant][HC], sp.KaonUncontain_ND [gen][1][quant][HC]}, "KaonUncontain_ND" +Comb, ND_CutNames );
        MakeComparisonPlots( {sp.ElecUncontain_ND [gen][0][quant][HC], sp.ElecUncontain_ND [gen][1][quant][HC]}, "ElecUncontain_ND" +Comb, ND_CutNames );
        MakeComparisonPlots( {sp.MiscUncontain_ND [gen][0][quant][HC], sp.MiscUncontain_ND [gen][1][quant][HC]}, "MiscUncontain_ND" +Comb, ND_CutNames );

        // energy fraction
        MakeComparisonPlots( {sp.MuonEFrac_ND    [gen][0][quant][HC], sp.MuonEFrac_ND    [gen][1][quant][HC]}, "MuonEFrac_ND"    +Comb, ND_CutNames);
        MakeComparisonPlots( {sp.NeutronEFrac_ND [gen][0][quant][HC], sp.NeutronEFrac_ND [gen][1][quant][HC]}, "NeutronEFrac_ND" +Comb, ND_CutNames);
        MakeComparisonPlots( {sp.ProtonEFrac_ND  [gen][0][quant][HC], sp.ProtonEFrac_ND  [gen][1][quant][HC]}, "ProtonEFrac_ND"  +Comb, ND_CutNames);
        MakeComparisonPlots( {sp.PionEFrac_ND    [gen][0][quant][HC], sp.PionEFrac_ND    [gen][1][quant][HC]}, "PionEFrac_ND"    +Comb, ND_CutNames);
        MakeComparisonPlots( {sp.Pi0EFrac_ND     [gen][0][quant][HC], sp.Pi0EFrac_ND     [gen][1][quant][HC]}, "Pi0EFrac_ND"     +Comb, ND_CutNames);
        MakeComparisonPlots( {sp.PhotonEFrac_ND  [gen][0][quant][HC], sp.PhotonEFrac_ND  [gen][1][quant][HC]}, "PhotonEFrac_ND"  +Comb, ND_CutNames);
        MakeComparisonPlots( {sp.KaonEFrac_ND    [gen][0][quant][HC], sp.KaonEFrac_ND    [gen][1][quant][HC]}, "KaonEFrac_ND"    +Comb, ND_CutNames);
        MakeComparisonPlots( {sp.ElectronEFrac_ND[gen][0][quant][HC], sp.ElectronEFrac_ND[gen][1][quant][HC]}, "ElectronEFrac_ND"+Comb, ND_CutNames);

        // uncont. E / tot. uncont. E
        MakeComparisonPlots( {sp.UncontMuonEF_ND[gen][0][quant][HC], sp.UncontMuonEF_ND[gen][1][quant][HC]}, "UncontMuonEF_ND"+Comb, ND_CutNames);
        MakeComparisonPlots( {sp.UncontNeutEF_ND[gen][0][quant][HC], sp.UncontNeutEF_ND[gen][1][quant][HC]}, "UncontNeutEF_ND"+Comb, ND_CutNames);
        MakeComparisonPlots( {sp.UncontProtEF_ND[gen][0][quant][HC], sp.UncontProtEF_ND[gen][1][quant][HC]}, "UncontProtEF_ND"+Comb, ND_CutNames);
        MakeComparisonPlots( {sp.UncontPionEF_ND[gen][0][quant][HC], sp.UncontPionEF_ND[gen][1][quant][HC]}, "UncontPionEF_ND"+Comb, ND_CutNames);
        MakeComparisonPlots( {sp.UncontPi0sEF_ND[gen][0][quant][HC], sp.UncontPi0sEF_ND[gen][1][quant][HC]}, "UncontPi0sEF_ND"+Comb, ND_CutNames);
        MakeComparisonPlots( {sp.UncontPhotEF_ND[gen][0][quant][HC], sp.UncontPhotEF_ND[gen][1][quant][HC]}, "UncontPhotEF_ND"+Comb, ND_CutNames);
        MakeComparisonPlots( {sp.UncontKaonEF_ND[gen][0][quant][HC], sp.UncontKaonEF_ND[gen][1][quant][HC]}, "UncontKaonEF_ND"+Comb, ND_CutNames);
        MakeComparisonPlots( {sp.UncontElecEF_ND[gen][0][quant][HC], sp.UncontElecEF_ND[gen][1][quant][HC]}, "UncontElecEF_ND"+Comb, ND_CutNames);

        // uncont. E / true nu E 
        MakeComparisonPlots( {sp.UncontTotEFT_ND [gen][0][quant][HC], sp.UncontTotEFT_ND [gen][1][quant][HC]}, "UncontTotEFT_ND" +Comb, ND_CutNames);
        MakeComparisonPlots( {sp.UncontMuonEFT_ND[gen][0][quant][HC], sp.UncontMuonEFT_ND[gen][1][quant][HC]}, "UncontMuonEFT_ND"+Comb, ND_CutNames);
        MakeComparisonPlots( {sp.UncontNeutEFT_ND[gen][0][quant][HC], sp.UncontNeutEFT_ND[gen][1][quant][HC]}, "UncontNeutEFT_ND"+Comb, ND_CutNames);
        MakeComparisonPlots( {sp.UncontProtEFT_ND[gen][0][quant][HC], sp.UncontProtEFT_ND[gen][1][quant][HC]}, "UncontProtEFT_ND"+Comb, ND_CutNames);
        MakeComparisonPlots( {sp.UncontPionEFT_ND[gen][0][quant][HC], sp.UncontPionEFT_ND[gen][1][quant][HC]}, "UncontPionEFT_ND"+Comb, ND_CutNames);
        MakeComparisonPlots( {sp.UncontPi0sEFT_ND[gen][0][quant][HC], sp.UncontPi0sEFT_ND[gen][1][quant][HC]}, "UncontPi0sEFT_ND"+Comb, ND_CutNames);
        MakeComparisonPlots( {sp.UncontPhotEFT_ND[gen][0][quant][HC], sp.UncontPhotEFT_ND[gen][1][quant][HC]}, "UncontPhotEFT_ND"+Comb, ND_CutNames);
        MakeComparisonPlots( {sp.UncontKaonEFT_ND[gen][0][quant][HC], sp.UncontKaonEFT_ND[gen][1][quant][HC]}, "UncontKaonEFT_ND"+Comb, ND_CutNames);
        MakeComparisonPlots( {sp.UncontElecEFT_ND[gen][0][quant][HC], sp.UncontElecEFT_ND[gen][1][quant][HC]}, "UncontElecEFT_ND"+Comb, ND_CutNames);

        // *************** Now look at plots where we only care about the final selection cut  ***************
        // total uncontained energy vs. vertex position
        Print2D(sp.UncontTotEvsX_FD[gen][FCut_FD][quant][HC], "UncontTotEvsX_FD"+Comb );
        Print2D(sp.UncontTotEvsY_FD[gen][FCut_FD][quant][HC], "UncontTotEvsY_FD"+Comb );

        // EF vs. true nuE
        Print2D(sp.MuonEFracvsTrueE_FD    [gen][FCut_FD][quant][HC], "MuonEFracvsTrueE_FD"    +Comb ); 
        Print2D(sp.NeutronEFracvsTrueE_FD [gen][FCut_FD][quant][HC], "NeutronEFracvsTrueE_FD" +Comb ); 
        Print2D(sp.ProtonEFracvsTrueE_FD  [gen][FCut_FD][quant][HC], "ProtonEFracvsTrueE_FD"  +Comb ); 
        Print2D(sp.PionEFracvsTrueE_FD    [gen][FCut_FD][quant][HC], "PionEFracvsTrueE_FD"    +Comb ); 
        Print2D(sp.Pi0EFracvsTrueE_FD     [gen][FCut_FD][quant][HC], "Pi0EFracvsTrueE_FD"     +Comb ); 
        Print2D(sp.PhotonEFracvsTrueE_FD  [gen][FCut_FD][quant][HC], "PhotonEFracvsTrueE_FD"  +Comb ); 
        Print2D(sp.KaonEFracvsTrueE_FD    [gen][FCut_FD][quant][HC], "KaonEFracvsTrueE_FD"    +Comb ); 
        Print2D(sp.ElectronEFracvsTrueE_FD[gen][FCut_FD][quant][HC], "ElectronEFracvsTrueE_FD"+Comb ); 

        // uncont EF vs. true nuE
        Print2D(sp.UncontMuonEFvsTrueE_FD[gen][FCut_FD][quant][HC], "UncontMuonEFvsTrueE_FD"+Comb ); 
        Print2D(sp.UncontNeutEFvsTrueE_FD[gen][FCut_FD][quant][HC], "UncontNeutEFvsTrueE_FD"+Comb ); 
        Print2D(sp.UncontProtEFvsTrueE_FD[gen][FCut_FD][quant][HC], "UncontProtEFvsTrueE_FD"+Comb ); 
        Print2D(sp.UncontPionEFvsTrueE_FD[gen][FCut_FD][quant][HC], "UncontPionEFvsTrueE_FD"+Comb ); 
        Print2D(sp.UncontPi0sEFvsTrueE_FD[gen][FCut_FD][quant][HC], "UncontPi0sEFvsTrueE_FD"+Comb ); 
        Print2D(sp.UncontPhotEFvsTrueE_FD[gen][FCut_FD][quant][HC], "UncontPhotEFvsTrueE_FD"+Comb ); 
        Print2D(sp.UncontKaonEFvsTrueE_FD[gen][FCut_FD][quant][HC], "UncontKaonEFvsTrueE_FD"+Comb ); 
        Print2D(sp.UncontElecEFvsTrueE_FD[gen][FCut_FD][quant][HC], "UncontElecEFvsTrueE_FD"+Comb ); 

        // uncont EFT vs. true nuE
        Print2D(sp.UncontTotEFTvsTrueE_FD [gen][FCut_FD][quant][HC], "UncontTotEFTvsTrueE_FD"+Comb ); 
        Print2D(sp.UncontMuonEFTvsTrueE_FD[gen][FCut_FD][quant][HC], "UncontMuonEFTvsTrueE_FD"+Comb ); 
        Print2D(sp.UncontNeutEFTvsTrueE_FD[gen][FCut_FD][quant][HC], "UncontNeutEFTvsTrueE_FD"+Comb ); 
        Print2D(sp.UncontProtEFTvsTrueE_FD[gen][FCut_FD][quant][HC], "UncontProtEFTvsTrueE_FD"+Comb ); 
        Print2D(sp.UncontPionEFTvsTrueE_FD[gen][FCut_FD][quant][HC], "UncontPionEFTvsTrueE_FD"+Comb ); 
        Print2D(sp.UncontPi0sEFTvsTrueE_FD[gen][FCut_FD][quant][HC], "UncontPi0sEFTvsTrueE_FD"+Comb ); 
        Print2D(sp.UncontPhotEFTvsTrueE_FD[gen][FCut_FD][quant][HC], "UncontPhotEFTvsTrueE_FD"+Comb ); 
        Print2D(sp.UncontKaonEFTvsTrueE_FD[gen][FCut_FD][quant][HC], "UncontKaonEFTvsTrueE_FD"+Comb ); 
        Print2D(sp.UncontElecEFTvsTrueE_FD[gen][FCut_FD][quant][HC], "UncontElecEFTvsTrueE_FD"+Comb ); 


        // --- Pull the ND plots.
        // total uncontained energy vs. vertex position
        Print2D(sp.UncontTotEvsX_ND[gen][FCut_ND][quant][HC], "UncontTotEvsX_ND"+Comb );
        Print2D(sp.UncontTotEvsY_ND[gen][FCut_ND][quant][HC], "UncontTotEvsY_ND"+Comb );

        // EF vs. true nuE
        Print2D(sp.MuonEFracvsTrueE_ND    [gen][FCut_ND][quant][HC], "MuonEFracvsTrueE_ND"    +Comb ); 
        Print2D(sp.NeutronEFracvsTrueE_ND [gen][FCut_ND][quant][HC], "NeutronEFracvsTrueE_ND" +Comb ); 
        Print2D(sp.ProtonEFracvsTrueE_ND  [gen][FCut_ND][quant][HC], "ProtonEFracvsTrueE_ND"  +Comb ); 
        Print2D(sp.PionEFracvsTrueE_ND    [gen][FCut_ND][quant][HC], "PionEFracvsTrueE_ND"    +Comb ); 
        Print2D(sp.Pi0EFracvsTrueE_ND     [gen][FCut_ND][quant][HC], "Pi0EFracvsTrueE_ND"     +Comb ); 
        Print2D(sp.PhotonEFracvsTrueE_ND  [gen][FCut_ND][quant][HC], "PhotonEFracvsTrueE_ND"  +Comb ); 
        Print2D(sp.KaonEFracvsTrueE_ND    [gen][FCut_ND][quant][HC], "KaonEFracvsTrueE_ND"    +Comb ); 
        Print2D(sp.ElectronEFracvsTrueE_ND[gen][FCut_ND][quant][HC], "ElectronEFracvsTrueE_ND"+Comb ); 

        // uncont EF vs. true nuE
        Print2D(sp.UncontMuonEFvsTrueE_ND[gen][FCut_ND][quant][HC], "UncontMuonEFvsTrueE_ND"+Comb ); 
        Print2D(sp.UncontNeutEFvsTrueE_ND[gen][FCut_ND][quant][HC], "UncontNeutEFvsTrueE_ND"+Comb ); 
        Print2D(sp.UncontProtEFvsTrueE_ND[gen][FCut_ND][quant][HC], "UncontProtEFvsTrueE_ND"+Comb ); 
        Print2D(sp.UncontPionEFvsTrueE_ND[gen][FCut_ND][quant][HC], "UncontPionEFvsTrueE_ND"+Comb ); 
        Print2D(sp.UncontPi0sEFvsTrueE_ND[gen][FCut_ND][quant][HC], "UncontPi0sEFvsTrueE_ND"+Comb ); 
        Print2D(sp.UncontPhotEFvsTrueE_ND[gen][FCut_ND][quant][HC], "UncontPhotEFvsTrueE_ND"+Comb ); 
        Print2D(sp.UncontKaonEFvsTrueE_ND[gen][FCut_ND][quant][HC], "UncontKaonEFvsTrueE_ND"+Comb ); 
        Print2D(sp.UncontElecEFvsTrueE_ND[gen][FCut_ND][quant][HC], "UncontElecEFvsTrueE_ND"+Comb ); 

        //  uncont EFT vs. true nuE
        Print2D(sp.UncontTotEFTvsTrueE_ND [gen][FCut_ND][quant][HC], "UncontTotEFTvsTrueE_ND"+Comb ); 
        Print2D(sp.UncontMuonEFTvsTrueE_ND[gen][FCut_ND][quant][HC], "UncontMuonEFTvsTrueE_ND"+Comb ); 
        Print2D(sp.UncontNeutEFTvsTrueE_ND[gen][FCut_ND][quant][HC], "UncontNeutEFTvsTrueE_ND"+Comb ); 
        Print2D(sp.UncontProtEFTvsTrueE_ND[gen][FCut_ND][quant][HC], "UncontProtEFTvsTrueE_ND"+Comb ); 
        Print2D(sp.UncontPionEFTvsTrueE_ND[gen][FCut_ND][quant][HC], "UncontPionEFTvsTrueE_ND"+Comb ); 
        Print2D(sp.UncontPi0sEFTvsTrueE_ND[gen][FCut_ND][quant][HC], "UncontPi0sEFTvsTrueE_ND"+Comb ); 
        Print2D(sp.UncontPhotEFTvsTrueE_ND[gen][FCut_ND][quant][HC], "UncontPhotEFTvsTrueE_ND"+Comb ); 
        Print2D(sp.UncontKaonEFTvsTrueE_ND[gen][FCut_ND][quant][HC], "UncontKaonEFTvsTrueE_ND"+Comb ); 
        Print2D(sp.UncontElecEFTvsTrueE_ND[gen][FCut_ND][quant][HC], "UncontElecEFTvsTrueE_ND"+Comb ); 

        // --- compare ND and FD -- these will only make sense if area normalized
        // --- uncontained energy

        MakeComparisonPlots({sp.AllUncontain_ND [gen][FCut_ND][quant][HC], sp.AllUncontain_FD [gen][FCut_FD][quant][HC]}, "AllUncontain_NDFD"+Comb, DetNames);
        MakeComparisonPlots({sp.MuonUncontain_ND[gen][FCut_ND][quant][HC], sp.MuonUncontain_FD[gen][FCut_FD][quant][HC]}, "MuonUncontain_NDFD"+Comb, DetNames); 
        MakeComparisonPlots({sp.NeutUncontain_ND[gen][FCut_ND][quant][HC], sp.NeutUncontain_FD[gen][FCut_FD][quant][HC]}, "NeutUncontain_NDFD"+Comb, DetNames); 
        MakeComparisonPlots({sp.ProtUncontain_ND[gen][FCut_ND][quant][HC], sp.ProtUncontain_FD[gen][FCut_FD][quant][HC]}, "ProtUncontain_ND[FD"+Comb, DetNames); 
        MakeComparisonPlots({sp.PionUncontain_ND[gen][FCut_ND][quant][HC], sp.PionUncontain_FD[gen][FCut_FD][quant][HC]}, "PionUncontain_NDFD"+Comb, DetNames); 
        MakeComparisonPlots({sp.Pi0sUncontain_ND[gen][FCut_ND][quant][HC], sp.Pi0sUncontain_FD[gen][FCut_FD][quant][HC]}, "Pi0sUncontain_NDFD"+Comb, DetNames); 
        MakeComparisonPlots({sp.PhotUncontain_ND[gen][FCut_ND][quant][HC], sp.PhotUncontain_FD[gen][FCut_FD][quant][HC]}, "PhotUncontain_NDFD"+Comb, DetNames); 
        MakeComparisonPlots({sp.KaonUncontain_ND[gen][FCut_ND][quant][HC], sp.KaonUncontain_FD[gen][FCut_FD][quant][HC]}, "KaonUncontain_NDFD"+Comb, DetNames); 
        MakeComparisonPlots({sp.ElecUncontain_ND[gen][FCut_ND][quant][HC], sp.ElecUncontain_FD[gen][FCut_FD][quant][HC]}, "ElecUncontain_NDFD"+Comb, DetNames); 
        MakeComparisonPlots({sp.MiscUncontain_ND[gen][FCut_ND][quant][HC], sp.MiscUncontain_FD[gen][FCut_FD][quant][HC]}, "MiscUncontain_NDFD"+Comb, DetNames); 

        // --- energy / true neutrino energy 
        MakeComparisonPlots({sp.MuonEFrac_ND    [gen][FCut_ND][quant][HC], sp.MuonEFrac_FD    [gen][FCut_FD][quant][HC]}, "MuonEFrac_NDFD"+Comb, DetNames);
        MakeComparisonPlots({sp.NeutronEFrac_ND [gen][FCut_ND][quant][HC], sp.NeutronEFrac_FD [gen][FCut_FD][quant][HC]}, "NeutronEFrac_NDFD"+Comb, DetNames);
        MakeComparisonPlots({sp.ProtonEFrac_ND  [gen][FCut_ND][quant][HC], sp.ProtonEFrac_FD  [gen][FCut_FD][quant][HC]}, "ProtonEFrac_NDFD"+Comb, DetNames);
        MakeComparisonPlots({sp.PionEFrac_ND    [gen][FCut_ND][quant][HC], sp.PionEFrac_FD    [gen][FCut_FD][quant][HC]}, "PionEFrac_NDFD"+Comb, DetNames);
        MakeComparisonPlots({sp.Pi0EFrac_ND     [gen][FCut_ND][quant][HC], sp.Pi0EFrac_FD     [gen][FCut_FD][quant][HC]}, "Pi0EFrac_NDFD"+Comb, DetNames);
        MakeComparisonPlots({sp.PhotonEFrac_ND  [gen][FCut_ND][quant][HC], sp.PhotonEFrac_FD  [gen][FCut_FD][quant][HC]}, "PhotonEFrac_NDFD"+Comb, DetNames);
        MakeComparisonPlots({sp.KaonEFrac_ND    [gen][FCut_ND][quant][HC], sp.KaonEFrac_FD    [gen][FCut_FD][quant][HC]}, "KaonEFrac_NDFD"+Comb, DetNames);
        MakeComparisonPlots({sp.ElectronEFrac_ND[gen][FCut_ND][quant][HC], sp.ElectronEFrac_FD[gen][FCut_FD][quant][HC]}, "ElectronEFrac_NDFD"+Comb, DetNames);
                                                                                                                       
        // --- uncontained energy / true neutrino energy                                                              
        MakeComparisonPlots({sp.UncontTotEFT_ND [gen][FCut_ND][quant][HC], sp.UncontTotEFT_FD [gen][FCut_FD][quant][HC]}, "UncontTotEFT_NDFD"+Comb, DetNames);
        MakeComparisonPlots({sp.UncontMuonEFT_ND[gen][FCut_ND][quant][HC], sp.UncontMuonEFT_FD[gen][FCut_FD][quant][HC]}, "UncontMuonEFT_NDFD"+Comb, DetNames);
        MakeComparisonPlots({sp.UncontNeutEFT_ND[gen][FCut_ND][quant][HC], sp.UncontNeutEFT_FD[gen][FCut_FD][quant][HC]}, "UncontNeutEFT_NDFD"+Comb, DetNames);
        MakeComparisonPlots({sp.UncontProtEFT_ND[gen][FCut_ND][quant][HC], sp.UncontProtEFT_FD[gen][FCut_FD][quant][HC]}, "UncontProtEFT_NDFD"+Comb, DetNames);
        MakeComparisonPlots({sp.UncontPionEFT_ND[gen][FCut_ND][quant][HC], sp.UncontPionEFT_FD[gen][FCut_FD][quant][HC]}, "UncontPionEFT_NDFD"+Comb, DetNames);
        MakeComparisonPlots({sp.UncontPi0sEFT_ND[gen][FCut_ND][quant][HC], sp.UncontPi0sEFT_FD[gen][FCut_FD][quant][HC]}, "UncontPi0sEFT_NDFD"+Comb, DetNames);
        MakeComparisonPlots({sp.UncontPhotEFT_ND[gen][FCut_ND][quant][HC], sp.UncontPhotEFT_FD[gen][FCut_FD][quant][HC]}, "UncontPhotEFT_NDFD"+Comb, DetNames);
        MakeComparisonPlots({sp.UncontKaonEFT_ND[gen][FCut_ND][quant][HC], sp.UncontKaonEFT_FD[gen][FCut_FD][quant][HC]}, "UncontKaonEFT_NDFD"+Comb, DetNames);
        MakeComparisonPlots({sp.UncontElecEFT_ND[gen][FCut_ND][quant][HC], sp.UncontElecEFT_FD[gen][FCut_FD][quant][HC]}, "UncontElecEFT_NDFD"+Comb, DetNames);
                                                                                                                       
      } // Loop over quants.
    } // Loop through GENIE interaction types.
  }// FHC/RHC

  // --- FHC/RHC comparisons
  for(unsigned int gen=0; gen<nGENIE; ++gen) {
    for (unsigned int quant=0; quant<nQuantNames; ++quant) {
      std::string Comb = "_"+GENIEStr[gen]+"_"+QuantNames[quant];
      // --- FD 
      // --- uncontained energy
      MakeComparisonPlots({sp.AllUncontain_FD[gen][FCut_FD][quant][0] , sp.AllUncontain_FD[gen][FCut_FD][quant][1]} , "AllUncontain_FD_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.MuonUncontain_FD[gen][FCut_FD][quant][0], sp.MuonUncontain_FD[gen][FCut_FD][quant][1]}, "MuonUncontain_FD_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.NeutUncontain_FD[gen][FCut_FD][quant][0], sp.NeutUncontain_FD[gen][FCut_FD][quant][1]}, "NeutUncontain_FD_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.ProtUncontain_FD[gen][FCut_FD][quant][0], sp.ProtUncontain_FD[gen][FCut_FD][quant][1]}, "ProtUncontain_FD_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.PionUncontain_FD[gen][FCut_FD][quant][0], sp.PionUncontain_FD[gen][FCut_FD][quant][1]}, "PionUncontain_FD_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.Pi0sUncontain_FD[gen][FCut_FD][quant][0], sp.Pi0sUncontain_FD[gen][FCut_FD][quant][1]}, "Pi0sUncontain_FD_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.PhotUncontain_FD[gen][FCut_FD][quant][0], sp.PhotUncontain_FD[gen][FCut_FD][quant][1]}, "PhotUncontain_FD_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.KaonUncontain_FD[gen][FCut_FD][quant][0], sp.KaonUncontain_FD[gen][FCut_FD][quant][1]}, "KaonUncontain_FD_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.ElecUncontain_FD[gen][FCut_FD][quant][0], sp.ElecUncontain_FD[gen][FCut_FD][quant][1]}, "ElecUncontain_FD_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.MiscUncontain_FD[gen][FCut_FD][quant][0], sp.MiscUncontain_FD[gen][FCut_FD][quant][1]}, "MiscUncontain_FD_FHCRHC"+Comb, HCNames);

      // --- energy / true neutrino energy
      MakeComparisonPlots({sp.MuonEFrac_FD    [gen][FCut_FD][quant][0], sp.MuonEFrac_FD    [gen][FCut_FD][quant][1]}, "MuonEFrac_FD_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.NeutronEFrac_FD [gen][FCut_FD][quant][0], sp.NeutronEFrac_FD [gen][FCut_FD][quant][1]}, "NeutronEFrac_FD_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.ProtonEFrac_FD  [gen][FCut_FD][quant][0], sp.ProtonEFrac_FD  [gen][FCut_FD][quant][1]}, "ProtonEFrac_FD_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.PionEFrac_FD    [gen][FCut_FD][quant][0], sp.PionEFrac_FD    [gen][FCut_FD][quant][1]}, "PionEFrac_FD_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.Pi0EFrac_FD     [gen][FCut_FD][quant][0], sp.Pi0EFrac_FD     [gen][FCut_FD][quant][1]}, "Pi0EFrac_FD_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.PhotonEFrac_FD  [gen][FCut_FD][quant][0], sp.PhotonEFrac_FD  [gen][FCut_FD][quant][1]}, "PhotonEFrac_FD_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.KaonEFrac_FD    [gen][FCut_FD][quant][0], sp.KaonEFrac_FD    [gen][FCut_FD][quant][1]}, "KaonEFrac_FD_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.ElectronEFrac_FD[gen][FCut_FD][quant][0], sp.ElectronEFrac_FD[gen][FCut_FD][quant][1]}, "ElectronEFrac_FD_FHCRHC"+Comb, HCNames);

      // --- uncontained energy / true neutrino energy
      MakeComparisonPlots({sp.UncontTotEFT_FD [gen][FCut_FD][quant][0], sp.UncontTotEFT_FD [gen][FCut_FD][quant][1]}, "UncontTotEFT_FD_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.UncontMuonEFT_FD[gen][FCut_FD][quant][0], sp.UncontMuonEFT_FD[gen][FCut_FD][quant][1]}, "UncontMuonEFT_FD_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.UncontNeutEFT_FD[gen][FCut_FD][quant][0], sp.UncontNeutEFT_FD[gen][FCut_FD][quant][1]}, "UncontNeutEFT_FD_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.UncontProtEFT_FD[gen][FCut_FD][quant][0], sp.UncontProtEFT_FD[gen][FCut_FD][quant][1]}, "UncontProtEFT_FD_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.UncontPionEFT_FD[gen][FCut_FD][quant][0], sp.UncontPionEFT_FD[gen][FCut_FD][quant][1]}, "UncontPionEFT_FD_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.UncontPi0sEFT_FD[gen][FCut_FD][quant][0], sp.UncontPi0sEFT_FD[gen][FCut_FD][quant][1]}, "UncontPi0sEFT_FD_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.UncontPhotEFT_FD[gen][FCut_FD][quant][0], sp.UncontPhotEFT_FD[gen][FCut_FD][quant][1]}, "UncontPhotEFT_FD_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.UncontKaonEFT_FD[gen][FCut_FD][quant][0], sp.UncontKaonEFT_FD[gen][FCut_FD][quant][1]}, "UncontKaonEFT_FD_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.UncontElecEFT_FD[gen][FCut_FD][quant][0], sp.UncontElecEFT_FD[gen][FCut_FD][quant][1]}, "UncontElecEFT_FD_FHCRHC"+Comb, HCNames);

      // --- ND
      // --- uncontained energy
      MakeComparisonPlots({sp.AllUncontain_ND[gen][FCut_ND][quant][0] , sp.AllUncontain_ND[gen][FCut_ND][quant][1]} , "AllUncontain_ND_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.MuonUncontain_ND[gen][FCut_ND][quant][0], sp.MuonUncontain_ND[gen][FCut_ND][quant][1]}, "MuonUncontain_ND_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.NeutUncontain_ND[gen][FCut_ND][quant][0], sp.NeutUncontain_ND[gen][FCut_ND][quant][1]}, "NeutUncontain_ND_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.ProtUncontain_ND[gen][FCut_ND][quant][0], sp.ProtUncontain_ND[gen][FCut_ND][quant][1]}, "ProtUncontain_ND_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.PionUncontain_ND[gen][FCut_ND][quant][0], sp.PionUncontain_ND[gen][FCut_ND][quant][1]}, "PionUncontain_ND_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.Pi0sUncontain_ND[gen][FCut_ND][quant][0], sp.Pi0sUncontain_ND[gen][FCut_ND][quant][1]}, "Pi0sUncontain_ND_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.PhotUncontain_ND[gen][FCut_ND][quant][0], sp.PhotUncontain_ND[gen][FCut_ND][quant][1]}, "PhotUncontain_ND_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.KaonUncontain_ND[gen][FCut_ND][quant][0], sp.KaonUncontain_ND[gen][FCut_ND][quant][1]}, "KaonUncontain_ND_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.ElecUncontain_ND[gen][FCut_ND][quant][0], sp.ElecUncontain_ND[gen][FCut_ND][quant][1]}, "ElecUncontain_ND_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.MiscUncontain_ND[gen][FCut_ND][quant][0], sp.MiscUncontain_ND[gen][FCut_ND][quant][1]}, "MiscUncontain_ND_FHCRHC"+Comb, HCNames);

      // --- energy / true neutrino energy
      MakeComparisonPlots({sp.MuonEFrac_ND    [gen][FCut_ND][quant][0], sp.MuonEFrac_ND    [gen][FCut_ND][quant][1]}, "MuonEFrac_ND_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.NeutronEFrac_ND [gen][FCut_ND][quant][0], sp.NeutronEFrac_ND [gen][FCut_ND][quant][1]}, "NeutronEFrac_ND_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.ProtonEFrac_ND  [gen][FCut_ND][quant][0], sp.ProtonEFrac_ND  [gen][FCut_ND][quant][1]}, "ProtonEFrac_ND_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.PionEFrac_ND    [gen][FCut_ND][quant][0], sp.PionEFrac_ND    [gen][FCut_ND][quant][1]}, "PionEFrac_ND_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.Pi0EFrac_ND     [gen][FCut_ND][quant][0], sp.Pi0EFrac_ND     [gen][FCut_ND][quant][1]}, "Pi0EFrac_ND_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.PhotonEFrac_ND  [gen][FCut_ND][quant][0], sp.PhotonEFrac_ND  [gen][FCut_ND][quant][1]}, "PhotonEFrac_ND_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.KaonEFrac_ND    [gen][FCut_ND][quant][0], sp.KaonEFrac_ND    [gen][FCut_ND][quant][1]}, "KaonEFrac_ND_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.ElectronEFrac_ND[gen][FCut_ND][quant][0], sp.ElectronEFrac_ND[gen][FCut_ND][quant][1]}, "ElectronEFrac_ND_FHCRHC"+Comb, HCNames);

      // --- uncontained energy / true neutrino energy
      MakeComparisonPlots({sp.UncontTotEFT_ND [gen][FCut_ND][quant][0], sp.UncontTotEFT_ND [gen][FCut_ND][quant][1]}, "UncontTotEFT_ND_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.UncontMuonEFT_ND[gen][FCut_ND][quant][0], sp.UncontMuonEFT_ND[gen][FCut_ND][quant][1]}, "UncontMuonEFT_ND_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.UncontNeutEFT_ND[gen][FCut_ND][quant][0], sp.UncontNeutEFT_ND[gen][FCut_ND][quant][1]}, "UncontNeutEFT_ND_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.UncontProtEFT_ND[gen][FCut_ND][quant][0], sp.UncontProtEFT_ND[gen][FCut_ND][quant][1]}, "UncontProtEFT_ND_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.UncontPionEFT_ND[gen][FCut_ND][quant][0], sp.UncontPionEFT_ND[gen][FCut_ND][quant][1]}, "UncontPionEFT_ND_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.UncontPi0sEFT_ND[gen][FCut_ND][quant][0], sp.UncontPi0sEFT_ND[gen][FCut_ND][quant][1]}, "UncontPi0sEFT_ND_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.UncontPhotEFT_ND[gen][FCut_ND][quant][0], sp.UncontPhotEFT_ND[gen][FCut_ND][quant][1]}, "UncontPhotEFT_ND_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.UncontKaonEFT_ND[gen][FCut_ND][quant][0], sp.UncontKaonEFT_ND[gen][FCut_ND][quant][1]}, "UncontKaonEFT_ND_FHCRHC"+Comb, HCNames);
      MakeComparisonPlots({sp.UncontElecEFT_ND[gen][FCut_ND][quant][0], sp.UncontElecEFT_ND[gen][FCut_ND][quant][1]}, "UncontElecEFT_ND_FHCRHC"+Comb, HCNames);

    } // quant
  } // GENIE

//  // print the 0bin fraction last so it doesn't get lost in the billion saved plots
//  for (unsigned int HC=0; HC<nHC; ++HC) {
//    // events in 0bin in uncontained energy / true neutrino energy histogram
//    for(unsigned int gen=0; gen<nGENIE; ++gen) {
//      // Loop through the quantile cuts.
//      for (unsigned int quant=0; quant<nQuantNames; ++quant) {
//          std::string Comb = "_"+GENIEStr[gen]+"_"+QuantNames[quant]+"_"+HCNames[HC];
//          std::cout<<"FD"<<std::endl;
//          Print0binFrac(sp.UncontTotEFT_FD [gen][FCut_FD][quant][HC], "UncontTotEFT"  +Comb);
//          Print0binFrac(sp.UncontMuonEFT_FD[gen][FCut_FD][quant][HC], "UncontMuonEFT" +Comb);
//          Print0binFrac(sp.UncontNeutEFT_FD[gen][FCut_FD][quant][HC], "UncontNeutEFT" +Comb);
//          Print0binFrac(sp.UncontProtEFT_FD[gen][FCut_FD][quant][HC], "UncontProtEFT" +Comb);
//          Print0binFrac(sp.UncontPionEFT_FD[gen][FCut_FD][quant][HC], "UncontPionEFT" +Comb);
//          Print0binFrac(sp.UncontPi0sEFT_FD[gen][FCut_FD][quant][HC], "UncontPi0sEFT" +Comb);
//          Print0binFrac(sp.UncontPhotEFT_FD[gen][FCut_FD][quant][HC], "UncontPhotEFT" +Comb);
//          Print0binFrac(sp.UncontKaonEFT_FD[gen][FCut_FD][quant][HC], "UncontKaonEFT" +Comb);
//          Print0binFrac(sp.UncontElecEFT_FD[gen][FCut_FD][quant][HC], "UncontElecEFT" +Comb);
//
//          std::cout<<"UncontND"<<std::endl;
//          Print0binFrac(sp.UncontTotEFT_ND [gen][FCut_ND][quant][HC], "UncontTotEFT"  +Comb);
//          Print0binFrac(sp.UncontMuonEFT_ND[gen][FCut_ND][quant][HC], "UncontMuonEFT" +Comb);
//          Print0binFrac(sp.UncontNeutEFT_ND[gen][FCut_ND][quant][HC], "UncontNeutEFT" +Comb);
//          Print0binFrac(sp.UncontProtEFT_ND[gen][FCut_ND][quant][HC], "UncontProtEFT" +Comb);
//          Print0binFrac(sp.UncontPionEFT_ND[gen][FCut_ND][quant][HC], "UncontPionEFT" +Comb);
//          Print0binFrac(sp.UncontPi0sEFT_ND[gen][FCut_ND][quant][HC], "UncontPi0sEFT" +Comb);
//          Print0binFrac(sp.UncontPhotEFT_ND[gen][FCut_ND][quant][HC], "UncontPhotEFT" +Comb);
//          Print0binFrac(sp.UncontKaonEFT_ND[gen][FCut_ND][quant][HC], "UncontKaonEFT" +Comb);
//          Print0binFrac(sp.UncontElecEFT_ND[gen][FCut_ND][quant][HC], "UncontElecEFT" +Comb);
//      } // Quant cut
//    } // GENIE Cuts
//  } // HC 

  return;
}

void GetPercentages	(	TH1 *	hist,
		std::string	PlotName
	)

Definition at line 1015 of file EnergyCont_plot.C.

References ana::bins, and MECModelEnuComparisons::i.

{

  double pop10 = 0;
  double pop50 = 0;
  double pop90 = 0;

  int bins = hist->GetNbinsX();
  double entriesnot0bin = hist->GetEntries()-hist->GetBinContent(1);
  double currentEntries = 0.0;

  // want to loop over all the of the events that are not in the 0 bin
  // find the energies at 10% of the events, 50%, 90%

  // skip the first bin
  for(int i = 2; i<bins; ++i){
    currentEntries += hist->GetBinContent(i); 
    // 90% of the uncontained events fall below this energy...
    if(currentEntries/entriesnot0bin >= 0.90) {
      pop90 = hist->GetXaxis()->GetBinCenter(i);
      break;
    }
    // ...and 50% of the uncontained events fall below this energy...
    else if(currentEntries/entriesnot0bin >=50){
      pop50 = hist->GetXaxis()->GetBinCenter(i);
    }

    // ...and 10% of the uncontained events fall below this energy.
    else if(currentEntries/entriesnot0bin>10){
      pop10 = hist->GetXaxis()->GetBinCenter(i);
    }
  }

}

void MakeComparisonPlots	(	std::vector< TH1D * >	TheHists,
		std::string	PlotName,
		std::vector< std::string >	Labels
	)

Definition at line 944 of file EnergyCont_plot.C.

References AddEntry(), Draw(), make_mec_shifts_plots::GetMaximum(), GetYaxis(), analysePickle::hist, cet::sqlite::max(), SaveAs(), SetBorderSize(), SetFillStyle(), SetLineColor(), and SetLogy().

Referenced by EnergyCont_plot().

                                                                                                       {
  
  // --- Make my canvas!
  TCanvas* ThisCan = new TCanvas( PlotName.c_str(), PlotName.c_str() );
  // --- Legend time!
  double HistLeg_X1 = 0.57, HistLeg_Y1 = 0.65, HistLeg_X2 = 0.88, HistLeg_Y2 = 0.88;
  TLegend* Leg = new TLegend(HistLeg_X1, HistLeg_Y1, HistLeg_X2, HistLeg_Y2 );

  double max = 0.0;
  // --- Loop through hists and plot them!
  for (size_t hist=0; hist<TheHists.size(); ++hist) {
    TheHists[hist] -> SetLineColor( 1+hist );
    // ROOT defines the plotting frame with the first histogram
    // if subsequent histograms have a larger maximum value than the first one, it will be cut off
    if(TheHists[hist]->GetMaximum() > max) max = TheHists[hist]->GetMaximum();

    if (hist==0) TheHists[hist] -> Draw("hist");
    else         TheHists[hist] -> Draw("hist same");
    Leg -> AddEntry( TheHists[hist], Labels[hist].c_str() );
  }
  Leg -> SetFillStyle(0);
  Leg -> SetBorderSize(0);
  Leg -> Draw();

  // set the new maximum
  TheHists[0]->SetMaximum(max*1.1);
  
  // --- Do I want to have a log axis?
  if ( PlotName.find("CompareUncontainEnergies") != std::string::npos) {
    for (size_t hist=0; hist<TheHists.size(); ++hist) TheHists[hist] -> GetYaxis() -> SetRangeUser( 10, TheHists[0]->GetMaximum() );
    ThisCan -> SetLogy();
  }
  

  // --- Save the canvas.
  ThisCan -> SaveAs( TString("Plots/")+TString(ThisCan->GetName())+TString(".pdf") );
  ThisCan -> SaveAs( TString("Plots/")+TString(ThisCan->GetName())+TString(".png") );

  // return
  return;
}

void Print0binFrac	(	TH1 *	hist,
		std::string	PlotName
	)

Definition at line 1001 of file EnergyCont_plot.C.

References om::cout, allTimeWatchdog::endl, and MECModelEnuComparisons::i.

{
  double nevts0 = hist->GetBinContent(1);
  double nevtsTot = 0;

  for(int i = 1; i < hist->GetNbinsX(); i++) {
    nevtsTot += hist->GetBinContent(i);
  }
  
  std::cout<<PlotName<<" | nevts0bin/nevtsTotal " << nevts0<<"/"<<nevtsTot<<"="<<nevts0/nevtsTot << std::endl;
}

void Print2D	(	TH2 *	hist,
		std::string	PlotName
	)

Definition at line 988 of file EnergyCont_plot.C.

References SaveAs().

Referenced by EnergyCont_plot().

{
  // --- Make my canvas!
  TCanvas* ThisCan = new TCanvas( PlotName.c_str(), PlotName.c_str() );
  hist->Draw("colz");
  ThisCan -> SaveAs( TString("Plots/")+TString(ThisCan->GetName())+TString(".pdf") );
  ThisCan -> SaveAs( TString("Plots/")+TString(ThisCan->GetName())+TString(".png") );
  return;

}

TH2 * Spec2DtoHist	(	TFile *	MyF,
		std::string	SpecName,
		double	POT,
		std::string	AxisTitle
	)

Definition at line 910 of file EnergyCont_plot.C.

References ana::CenterTitles(), MECModelEnuComparisons::i, febshutoff_auto::integral, calib::j, ana::Spectrum::LoadFrom(), and ana::ToTH2().

Referenced by EnergyCont_plot().

                                                                                     {
  std::unique_ptr<Spectrum> TempSpec = Spectrum::LoadFrom(MyF,  SpecName.c_str() ) ;
  TH2* MyHist = TempSpec -> ToTH2( POT );
  ana::CenterTitles      ( MyHist);
  MyHist->SetMinimum( 0     );
  MyHist->SetMaximum( MyHist->GetMaximum()*1.1 );
  MyHist->SetTitle  ( AxisTitle.c_str() );

  // borrowed from Michael -- thanks!!
  // Loop over the histos to normalize them in vertical strips one bin wide.
  unsigned int NXbins = 0.0;
  unsigned int NYbins = 0.0;

  // W0, X0
  NXbins = MyHist->GetNbinsX();
  NYbins = MyHist->GetNbinsY();

  for(unsigned int i = 0; i <= NXbins+1; ++i) {
    float integral = 0.0;
    for(unsigned int j = 2; j <= NYbins; ++j) // 0 bin swamping the plot... remove the 0bin 
      integral += MyHist->GetBinContent(i,j);
    if(integral == 0.0) continue;
    for(unsigned int j = 0; j <= NYbins+1; ++j)
      MyHist->SetBinContent(i,j,MyHist->GetBinContent(i,j)/integral);
  }

  // the 0 bin is swamping the plot... remove it
  MyHist->GetYaxis()->SetRange(2,NYbins);

  // since it's normalized, need to set the z-axis to go from 0-1
  MyHist->GetZaxis()->SetRangeUser(0.0, 1.0);
  return MyHist;
}

TH1D * SpecToHist	(	TFile *	MyF,
		std::string	SpecName,
		double	POT,
		std::string	AxisTitle
	)

Definition at line 876 of file EnergyCont_plot.C.

References ana::CenterTitles(), and ana::Spectrum::LoadFrom().

Referenced by EnergyCont_plot().

                                                                                    {
  std::unique_ptr<Spectrum> TempSpec = Spectrum::LoadFrom(MyF,  SpecName.c_str() ) ;
  TH1D* MyHist = TempSpec -> ToTH1( POT );
  // area norm
//  MyHist->Scale(1./MyHist->Integral());
  // artificially set the 0 bin to have no entries AFTER normalization
  // use only if you want to remove the 0 bin
  MyHist->SetBinContent(1, 0);

  ana::CenterTitles      ( MyHist);
  MyHist->SetMinimum( 0     );
  MyHist->SetMaximum( MyHist->GetMaximum()*1.1 );
  MyHist->SetTitle  ( AxisTitle.c_str() );
 
  // I want to check what my maximum bin is...
  // this is clever but if we're trying to compare two sets of things and one has the zero bin and one doesn't it becomes difficult
//  int MaxBin = MyHist->GetMaximumBin();
//  int MaxVal = MyHist->GetMaximum();
//  int SecBin = MyHist->GetBinContent(2);
//  if ( MaxBin == 1 && (MaxVal/10 > SecBin) ) {
//    double NewMax = 0;
//    for (int bin=2; bin<MyHist->GetNbinsX(); ++bin) {
//      if (MyHist->GetBinContent(bin) > NewMax) NewMax = 1.1*MyHist->GetBinContent(bin);
//    }
//    MyHist->SetMaximum( NewMax );
//
//  }

  
 
  return MyHist;
}

Variable Documentation

std::vector<std::string> DetNames = {"ND", "FD"}

Definition at line 45 of file EnergyCont_plot.C.

Referenced by EnergyCont_plot().

unsigned int FCut_FD = nFD_CutNames - 1

Definition at line 36 of file EnergyCont_plot.C.

Referenced by EnergyCont_plot().

unsigned int FCut_ND = nND_CutNames - 1

Definition at line 36 of file EnergyCont_plot.C.

Referenced by EnergyCont_plot().

std::vector<std::string> FD_CutNames = { "FD_Contain", "FD_PID2018", "FD_CosRej" }

Definition at line 34 of file EnergyCont_plot.C.

Referenced by EnergyCont_plot().

const std::vector<std::string> GENIEStr = { "All" }

Definition at line 30 of file EnergyCont_plot.C.

Referenced by EnergyCont_macro(), EnergyCont_plot(), hyperon_nom_macro(), neutKEsyst(), and neutronE_macro().

std::vector<std::string> HCNames = {"fhc", "rhc"}

Definition at line 42 of file EnergyCont_plot.C.

Referenced by EnergyCont_plot().

std::vector<std::string> ND_CutNames = { "ND_Contain", "ND_PID2018" }

Definition at line 35 of file EnergyCont_plot.C.

Referenced by EnergyCont_plot().

const unsigned int nFD_CutNames = 3

Definition at line 34 of file EnergyCont_plot.C.

Referenced by EnergyCont_plot().

const unsigned int nGENIE = 1

Definition at line 30 of file EnergyCont_plot.C.

Referenced by EnergyCont_plot(), and hyperon_nom_macro().

const unsigned int nHC = 2

Definition at line 42 of file EnergyCont_plot.C.

Referenced by EnergyCont_plot().

const unsigned int nND_CutNames = 2

Definition at line 35 of file EnergyCont_plot.C.

Referenced by EnergyCont_plot().

const unsigned int nQuantNames = 5

Definition at line 39 of file EnergyCont_plot.C.

Referenced by EnergyCont_plot().

std::vector<std::string> QuantNames = { "AllQuant", "Quant1", "Quant2", "Quant3", "Quant4" }

Definition at line 39 of file EnergyCont_plot.C.

Referenced by EnergyCont_plot().