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

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

  TFile *f1= new TFile ("/nova/app/users/acedeno/tag_releasesS18-06-14/NDAna/ncpi0_semi_inc_png_cvn/BDTSpectrum_FlatFluxdatasetslh15.root","READ");
  f1->cd(); 


  
  ///load reco mass 2 prong CVN cuts spectrums///// (stacked) /////
  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"));
 

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

  ///load histogram /////
  
  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);
  

   TString outDir2 = "/nova/app/users/acedeno/tag_releasesS18-06-14/NDAna/ncpi0_semi_inc_png_cvn/BDTGStudy";
  //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 =2;
  const FileUp CurrentFile[File] = {
    {"/BDTG_Distribution"},
    {"/BDTG_Distribution_stacked"},
    
    
  };

   const int Date = 1;
   const DateUp CurrentDate[Date] = {
     {"_07-03-18_BDTG_FlatFlux_slh_15"},
  };


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

/// Plot 6 ///
  //gStyle->SetTitleY(0.9);
   TCanvas *cs6 = new TCanvas("cs6","cs6",10,10,855,500);
   cs6->cd();
  
   hndNCPi0LRem->SetLineColor(kOrange);
   hndNCPi0LRem->SetLineWidth(1);

   hndCCPi0Rem->SetLineColor(kMagenta);
   hndCCPi0Rem->SetLineWidth(1);

   
   hndCCNonPi0Rem->SetLineColor(kAzure);
   hndCCNonPi0Rem->SetLineWidth(1);
  
   
   hndNCNonPi0Rem->SetLineColor(kCyan);
   hndNCNonPi0Rem->SetLineWidth(1);
 
   
   hndNCPi0Rem->SetLineColor(kGreen);
   hndNCPi0Rem->SetLineWidth(1);

  
   hndMC3PRem->SetLineColor(kRed);
   hndMC3PRem->SetLineWidth(1);
   //hndNumiRemIDData->SetMarkerStyle(2);
   //hndNumiRemIDData->SetMarkerSize(10);
 


   hndMC3PRem->Draw("hist E1");
   hndCCPi0Rem->Draw("hist E1 same");
   hndCCNonPi0Rem->Draw("hist E1 same");
   hndNCPi0LRem->Draw("hist E1 same");
   hndNCPi0Rem->Draw("hist E1 same");
   hndNCNonPi0Rem->Draw("hist E1 same");
   hndMC3PRem->SetMaximum(6000);
   hndMC3PRem->GetXaxis()->SetRangeUser(-1.,1.);
   hndMC3PRem->GetYaxis()->SetTitle("Events/ 8e20 POT");

   //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[0].name.c_str()+CurrentDate[0].date.c_str()+".png");

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

    /// Plot 4 ///
 
   TCanvas *cs4 = new TCanvas("cs4","cs4",10,10,855,500);
   THStack *hs4 = new THStack("hs4",""); 

   hndNCPi0LRem->SetLineColor(kBlack);
   hndNCPi0LRem->SetFillColor(kOrange);
   hndNCPi0LRem->SetLineWidth(1);;
   hndNCPi0LRem->SetMarkerStyle(21);
   hndNCPi0LRem->SetMarkerColor(kOrange);
   hs4->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);

   hs4->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);
   hs4->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);
   hs4->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);
   hs4->Add(hndNCPi0Rem);
   
   //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("BDTG");
   hs4->SetMaximum(6000);
   hs4->GetXaxis()->SetRangeUser(-1.,1.);
   gPad->Modified();
   cs4->cd();
   //cs2->SetTopMargin(0.18);
   hndMC3PRem->SetLineColor(kRed);
   hndMC3PRem->SetLineWidth(3);
   //hndNumiCVN2PRData->SetMarkerStyle(2);
   //hndNumiCVN2PRData->SetMarkerSize(3);


   hndMC3PRem->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(hndMC3PRem,"Selection","kFullDotSmall");
   leg4->AddEntry(hndCCPi0Rem,"CCPi0","kFullDotSmall");
   leg4->AddEntry(hndCCNonPi0Rem,"CCNonPi0","kFullDotSmall");
   leg4->AddEntry(hndNCPi0LRem,"NC w/Pi0 < 0.3 E","kFullDotSmall");
   leg4->AddEntry(hndNCPi0Rem,"NC w/Pi0 > 0.3 E","kFullDotSmall");
   leg4->AddEntry(hndNCNonPi0Rem,"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[1].name.c_str()+CurrentDate[0].date.c_str()+".png");
    /*********************************************************************************************************************/

    f1->Close();

}
#endif