containmentstudy.C

Go to the documentation of this file.

// Do a proper prediction, with all the components included

#ifdef __CINT__
void containmentstudy()
{
  std::cout << "Sorry, you must run in compiled mode" << std::endl;
}
#else


////////CAFAna libraries///
#include "CAFAna/Cuts/Cuts.h"
#include "CAFAna/Analysis/Plots.h"
#include "CAFAna/Core/SpectrumLoaderMockData.h"
#include "CAFAna/Prediction/PredictionNoExtrap.h"
#include "CAFAna/Experiment/SingleSampleExperiment.h"
#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Core/SpectrumLoader.h"
#include "CAFAna/Vars/Vars.h"
#include "StandardRecord/Proxy/SRProxy.h" 
////////////////////////////
 
////Plotting Libraries/////
#include "TCanvas.h"
#include "TH2.h"
#include "TH1.h"
#include "TFitResultPtr.h"
#include "TLegend.h"
//////////////////////////

///////NCPi0 Headers//////
#include "/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/ncpi0Bins.h"
#include "/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/ncpi0Cuts.h"
#include "/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/ncpi0Vars.h"
/////////////////////////

using namespace ana;
  

    
void containmentstudy()

{
        
  /*********************************************************************************************************************/

        //Enviroment variables and wildcards work

        const std::string MC = "dataset_def_name_newest_snapshot prod_caf_R17-03-01-prod3reco.d_nd_genie_nonswap_fhc_nova_v08_full_v1";
        //const std::string Numi = "prod_caf_R17-03-01-prod3reco.d_nd_numi_fhc_full_v1_goodruns";

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

        //Create loader for each 'data' type
        
        SpectrumLoader mcloader(MC);
        //SpectrumLoader dataloader(Numi);
  
  /*********************************************************************************************************************/

        Spectrum ndMC("shwlid.stop in x-dir [cm]",bins_for_vtxXY, mcloader,ContStopX , kPreSelect);
        Spectrum ndMC3P("shwlid.stop in x-dir [cm]",bins_for_vtxXY,mcloader,ContStopX,k2ProngPreSel2CVN);
        //Spectrum ndNCMCSig("Prong 3 Photon ID",PCIDbins,mcloader,kNpCVN3gamma,kPreSelSig2P);
        Spectrum ndNCPi0("shwlid.stop in x-dir [cm]",bins_for_vtxXY,mcloader,ContStopX,kIsNCPi02P); 
        Spectrum ndNCPi0L("shwlid.stop in x-dir [cm]",bins_for_vtxXY,mcloader,ContStopX,kIsNCPi0L2P);
        Spectrum ndCCPi0("shwlid.stop in x-dir [cm]",bins_for_vtxXY, mcloader, ContStopX,kIsCCPi02P);
        Spectrum ndNCNonPi0("shwlid.stop in x-dir [cm]",bins_for_vtxXY,mcloader,ContStopX , kIsNCNonPi02P);
        Spectrum ndCCNonPi0("shwlid.stopin x-dir [cm]",bins_for_vtxXY,mcloader,ContStopX , kIsCCNonPi02P);


        Spectrum ndMCy("shwlid.stop in y-dir [cm]",bins_for_vtxXY, mcloader,ContStopY , kPreSelect);
        Spectrum ndMC3Py("shwlid.stop in y-dir [cm]",bins_for_vtxXY,mcloader,ContStopY ,k2ProngPreSel2CVN);
        //Spectrum ndNCMCSig("Prong 3 Photon ID",PCIDbins,mcloader,kNpCVN3gamma,kPreSelSig2P);
        Spectrum ndNCPi0y("shwlid.stop in y-dir [cm]",bins_for_vtxXY,mcloader,ContStopY,kIsNCPi02P); 
        Spectrum ndNCPi0Ly("shwlid.stop in y-dir [cm]",bins_for_vtxXY,mcloader,ContStopY ,kIsNCPi0L2P);
        Spectrum ndCCPi0y("shwlid.stop in y-dir [cm]",bins_for_vtxXY, mcloader, ContStopY,kIsCCPi02P);
        Spectrum ndNCNonPi0y("shwlid.stop in y-dir [cm]",bins_for_vtxXY,mcloader,ContStopY , kIsNCNonPi02P);
        Spectrum ndCCNonPi0y("shwlid.stop in y-dir [cm]",bins_for_vtxXY,mcloader,ContStopY , kIsCCNonPi02P);

        Spectrum ndMCz("shwlid.stop in z-dir [cm]",bins_for_vtxZ, mcloader, ContStopZ , kPreSelect);
        Spectrum ndMC3Pz("shwlid.stopin z-dir [cm]",bins_for_vtxZ,mcloader,  ContStopZ,k2ProngPreSel2CVN);
        //Spectrum ndNCMCSig("Prong 3 Photon ID",PCIDbins,mcloader,kNpCVN3gamma,kPreSelSig2P);
        Spectrum ndNCPi0z("shwlid.stop in z-dir [cm]",bins_for_vtxZ,mcloader, ContStopZ ,kIsNCPi02P); 
        Spectrum ndNCPi0Lz("shwlid.stop in z-dir [cm]",bins_for_vtxZ,mcloader, ContStopZ ,kIsNCPi0L2P);
        Spectrum ndCCPi0z("shwlid.stopin z-dir [cm]",bins_for_vtxZ, mcloader,  ContStopZ,kIsCCPi02P);
        Spectrum ndNCNonPi0z("shwlid.stop in z-dir [cm]",bins_for_vtxZ,mcloader,  ContStopZ, kIsNCNonPi02P);
        Spectrum ndCCNonPi0z("shwlid.stop in z-dir [cm]",bins_for_vtxZ,mcloader, ContStopZ , kIsCCNonPi02P);
        
        
        /********************************************/
        mcloader.Go();
        /*******************************************/

        TFile *out = new TFile("xyzcontainmentstudy1DSpectrumsApril2018.root","RECREATE");


        ///////reco masss with &  Presel+ Remid+ prong 1 CVN cut////////
        ndMC.SaveTo(out->mkdir("ndMC"));
        ndMC3P.SaveTo(out->mkdir("ndMC3P"));
        ndNCPi0.SaveTo(out->mkdir("ndNCPi0"));
        ndCCPi0.SaveTo(out->mkdir("ndCCPi0"));
        ndNCNonPi0.SaveTo(out->mkdir("ndNCNonPi0"));
        ndNCPi0L.SaveTo(out->mkdir("ndNCPi0L"));
        ndCCNonPi0.SaveTo(out->mkdir("ndCCNonPi0"));

        ndMCy.SaveTo(out->mkdir("ndMCy"));
        ndMC3Py.SaveTo(out->mkdir("ndMC3Py"));
        ndNCPi0y.SaveTo(out->mkdir("ndNCPi0y"));
        ndCCPi0y.SaveTo(out->mkdir("ndCCPi0y"));
        ndNCNonPi0y.SaveTo(out->mkdir("ndNCNonPi0y"));
        ndNCPi0Ly.SaveTo(out->mkdir("ndNCPi0Ly"));
        ndCCNonPi0y.SaveTo(out->mkdir("ndCCNonPi0y"));


        ndMCz.SaveTo(out->mkdir("ndMCz"));
        ndMC3Pz.SaveTo(out->mkdir("ndMC3Pz"));
        ndNCPi0z.SaveTo(out->mkdir("ndNCPi0z"));
        ndCCPi0z.SaveTo(out->mkdir("ndCCPi0z"));
        ndNCNonPi0z.SaveTo(out->mkdir("ndNCNonPi0z"));
        ndNCPi0Lz.SaveTo(out->mkdir("ndNCPi0Lz"));
        ndCCNonPi0z.SaveTo(out->mkdir("ndCCNonPi0z"));


        std::cout<<"\nOutput saved to " << out->GetName()<<std::endl;

        
        out->Close();

        return;
        
}  

#endif