make_DiF_plots.C

Go to the documentation of this file.

/////////////////////////////////////////////////
//       DISCLAIMER


//   I couldn't find an easy way to automate the yaxis upper limit in a way that 
//   allowed for the legend to not get in the way because of all the rebinning
//   and normalization.
//   Therefore, unfortunately, you kind of have to iteratively go through and manually 
//   set the y axis for each plot using the structure created at L87

/////////////////////////////////////////////////

#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/Spectrum.h"
#include "CAFAna/Analysis/Plots.h"
#include "CAFAna/Analysis/Style.h"
#include "TColor.h"

using namespace ana;

void myPreliminary()
{
  TLatex* prelim = new TLatex(.9, .95, "NOvA Preliminary");
  prelim->SetTextColor(kBlue);
  prelim->SetNDC();
  prelim->SetTextSize(2/30.);
  prelim->SetTextAlign(32);
  prelim->Draw();
}

//the DrawBeamLabel in Plots.h works for the other plots
//but not quite for the efficiency plots
void DrawBeamLabelEff(bool isFHC)
{
  TPaveText *pText = new TPaveText (0.1, 0.94, 0.16, 0.98,"NDC");
  pText->SetFillStyle(0);
  pText->SetLineColor(0);
  pText->SetLineWidth(0);
  pText->SetBorderSize(1);
  pText->SetTextColor(kGray+2);
  //pText->SetNDC();
  if (isFHC) pText->AddText("#nu-beam");
  else       pText->AddText("#bar{#nu}-beam"); 
  pText->SetTextSize(2/40.);  
  pText->SetTextAlign(11);
  pText->Draw();
}


//same list from the Step 2 macro, but different structure as 
//we want to automate different things now
//rebinAmt = 0 will mean no rebinning
struct Plot{const char* savename; const char* XaxisLabel; int rebinAmt; int legSide; int yAxisLimit;};
const std::vector<Plot> plots = {
  {"Shw_Angle",       "Cos(#theta_{beam})",      0, 0, 7000},
  {"NHits_in_slice",  "Hits in Slice",           4, 1, 8000},
  {"NHits_in_shower", "Hits in Shower",          4, 1, 8000},
  {"Shw_E",           "Energy in Shower (GeV)",  8, 1, 6000},
  {"E_per_Hit",       "GeV / Hit",               0, 1, 0},
  {"CVNe",            "CVN #nu_{e} selector",    4, 0, 11000},
  {"Shw_Width",       "Shower Width (cm)",       2, 1, 8000},
  {"Shw_Length",      "Shower Length (cm)",      2, 1, 8000},
  {"Shw_Start X",     "Shower Start X",          32, 0, 3000},
  {"Shw_Start Y",     "Shower Start Y",          32, 0, 0},
  {"Shw_Start Z",     "Shower Start Z",          50, 0, 1500},
  {"Vtx_X",           "Vertex Location X",       32, 0, 3000},
  {"Vtx_Y",           "Vertex Location Y",       32, 0, 0},
  {"Vtx_Z",           "Vertex Location Z",       50, 0, 1500},
  {"Distance_Top",    "Distance to Top (cm)",    100, 1, 4000},
  {"Distance_Bottom", "Distance to Bottom (cm)", 100, 0, 2500},
  {"Distance_East",   "Distance to East Face",   16, 1, 5000},
  {"Distance_West",   "Distance to West Face",   16, 1, 5000},
  {"Distance_Front",  "Distance to Front",       50, 0, 1500},
  {"Distance_Back",   "Distance to Back",        100, 0, 2800}
};

void makePlot(TCanvas *c1, TH1 *hCosmicsData, TH1 *hCRYMC, TH1 *hGENIEFluxswap, int sel, int legSide, bool isFHC);
void makeWgtdPlot(TCanvas *c1, TH1 *hCosmicsData, TH1 *hCRYMC, TH1 *hGENIEFluxswap, int sel, int legSide, bool isFHC);
void makePlot(TCanvas *c1, TH1 *hCosmicsData, TH1 *hCRYMC, int sel, int legSide);
void makeWgtdPlot(TCanvas *c1, TH1 *hCosmicsData, TH1 *hCRYMC, int sel, int legSide);
void makeDataCryRatio(TCanvas *c1, TH1 *hCosmicsData, TH1 *hCRYMC, int sel, int legSide, bool isFHC, bool isWgtd);
void makeEffPlot(TCanvas *c1,  TH1 *hCosmicsDataCore, TH1 *hCosmicsDataPresel, TH1 *hCRYMCCore, TH1 *hCRYMCPresel, int legSide, bool isFHC);
void makeWgtdEffPlot(TCanvas *c1,  TH1 *hCosmicsDataCore, TH1 *hCosmicsDataPresel, TH1 *hCRYMCCore, TH1 *hCRYMCPresel, int legSide, bool isFHC);
void makeEffPlotNoRatio(TCanvas *c1,  TH1 *hCosmicsDataCEff, TH1 *hCRYMCCoreEff, int legSide, bool isFHC, bool isWgtd);


void make_eff_plots_areaNorm(std::string beam = "FHC")
{



  //get input file opened
  std::string inFile = "/nova/app/users/amhall/DiF_Cafana_4-6-2020/3FlavorAna/Ana2020/MRstudies/MRDiF_FD_"+beam+"_Step2_v2_clean.root";
  TFile *f1 = new TFile(inFile.c_str());
  f1->cd();

  //get just a couple of the spectra to set our livetime and POT before getting the rest
  auto sForLivetime=*Spectrum::LoadFrom(f1, "sDataPresel_V0");
  auto sForPOT=*Spectrum::LoadFrom(f1, "sGENIEPresel_V0");

  double mylivetime=sForLivetime.Livetime();
  double myPOT=sForPOT.POT();
  std::cout<<"myPOT:"<<myPOT<<" mylivetime:"<<mylivetime<<std::endl;

  std::vector<TH1D*> hDataPresel, hCRYMCPresel, hGENIEPresel, hData, hCRYMC, hGENIE;
  std::vector<TH1D*> hDataPreselWtd, hCRYMCPreselWtd, hDataWtd, hCRYMCWtd;
  std::vector<TH1D*> hDataEff, hDataWtdEff, hCRYMCEff, hCRYMCWtdEff;
  std::vector<TH1D*> hDataPreselEff, hDataPreselWtdEff, hCRYMCPreselEff, hCRYMCPreselWtdEff;

  char dirname[64];  
  for( unsigned int i =0; i < plots.size(); i++)
    {
      //load in Data histos
      
      sprintf(dirname, "sDataPresel_V%d", i);
      std::unique_ptr<Spectrum> sDataPresel = Spectrum::LoadFrom(f1, dirname);
      hDataPresel.push_back(sDataPresel->ToTH1(mylivetime, kLivetime));
      
      sprintf(dirname, "sData_V%d", i);
      std::unique_ptr<Spectrum> sData = Spectrum::LoadFrom(f1, dirname);
      hData.push_back(sData->ToTH1(mylivetime, kLivetime));

      sprintf(dirname, "sDataPreselWtd_V%d", i);
      std::unique_ptr<Spectrum> sDataPreselWtd = Spectrum::LoadFrom(f1, dirname);
      hDataPreselWtd.push_back(sDataPreselWtd->ToTH1(mylivetime, kLivetime));

      sprintf(dirname, "sDataWtd_V%d", i);
      std::unique_ptr<Spectrum> sDataWtd = Spectrum::LoadFrom(f1, dirname);
      hDataWtd.push_back(sDataWtd->ToTH1(mylivetime, kLivetime));

      //and the CRYMC histos

      sprintf(dirname, "sCRYMCPresel_V%d", i);
      std::unique_ptr<Spectrum> sCRYMCPresel = Spectrum::LoadFrom(f1, dirname);
      hCRYMCPresel.push_back(sCRYMCPresel->ToTH1(mylivetime, kLivetime));

      sprintf(dirname, "sCRYMC_V%d", i);
      std::unique_ptr<Spectrum> sCRYMC = Spectrum::LoadFrom(f1, dirname);
      hCRYMC.push_back(sCRYMC->ToTH1(mylivetime, kLivetime));

      sprintf(dirname, "sCRYMCPreselWtd_V%d", i);
      std::unique_ptr<Spectrum> sCRYMCPreselWtd = Spectrum::LoadFrom(f1, dirname);
      hCRYMCPreselWtd.push_back(sCRYMCPreselWtd->ToTH1(mylivetime, kLivetime));

      sprintf(dirname, "sCRYMCWtd_V%d", i);
      std::unique_ptr<Spectrum> sCRYMCWtd = Spectrum::LoadFrom(f1, dirname);
      hCRYMCWtd.push_back(sCRYMCWtd->ToTH1(mylivetime, kLivetime));

      //and the GENIE histos
      
      sprintf(dirname, "sGENIEPresel_V%d", i);
      std::unique_ptr<Spectrum> sGENIEPresel = Spectrum::LoadFrom(f1, dirname);
      hGENIEPresel.push_back(sGENIEPresel->ToTH1(myPOT));

      sprintf(dirname, "sGENIE_V%d", i);
      std::unique_ptr<Spectrum> sGENIE = Spectrum::LoadFrom(f1, dirname);
      hGENIE.push_back(sGENIE->ToTH1(myPOT));

      //and get copies of the data and CRY selected histos to make efficiencies with
      hDataEff.push_back(sData->ToTH1(mylivetime, kLivetime));
      hDataWtdEff.push_back(sDataWtd->ToTH1(mylivetime, kLivetime));
      hDataPreselEff.push_back(sDataPresel->ToTH1(mylivetime, kLivetime));
      hDataPreselWtdEff.push_back(sDataPreselWtd->ToTH1(mylivetime, kLivetime));

      hCRYMCEff.push_back(sCRYMC->ToTH1(mylivetime, kLivetime));
      hCRYMCWtdEff.push_back(sCRYMCWtd->ToTH1(mylivetime, kLivetime));
      hCRYMCPreselEff.push_back(sCRYMCPresel->ToTH1(mylivetime, kLivetime));
      hCRYMCPreselWtdEff.push_back(sCRYMCPreselWtd->ToTH1(mylivetime, kLivetime));
     
    }

  //kNormFactor = hDataPresel[0]->Integral(); //will area normalize to the cosmic data preselected number of events
 
 
  std::cout<<"NO NORM FACTORS HERE \n";
  std::cout<<"===================================== Before reweighting================================================"<<std::endl;
  std::cout<<"             Cosmic Data "<<"           CRYMC         "<<"                GENIE            "<<std::endl;
  std::cout<<"========================================================================================================="<<std::endl;
  std::cout<<"CorePresel      :"<<hDataPresel[0]->Integral()<<"               "<<hCRYMCPresel[0]->Integral()<<"              "<< hGENIEPresel[0]->Integral()<<std::endl;
  std::cout<<"Core            :"<<hData[0]->Integral()<<"               "<<hCRYMC[0]->Integral()<<"              "<< hGENIE[0]->Integral()<<std::endl;

  std::cout<<"NO NORM FACTORS HERE \n";
  std::cout<<"===================================== After reweighting================================================"<<std::endl;
  std::cout<<"             Cosmic Data "<<"           CRYMC         "<<"                GENIE            "<<std::endl;
  std::cout<<"========================================================================================================="<<std::endl;
  std::cout<<"CorePresel      :"<<hDataPreselWtd[0]->Integral()<<"               "<<hCRYMCPreselWtd[0]->Integral()<<"              "<< hGENIEPresel[0]->Integral()<<std::endl;
  std::cout<<"Core            :"<<hDataWtd[0]->Integral()<<"               "<<hCRYMCWtd[0]->Integral()<<"              "<< hGENIE[0]->Integral()<<std::endl;

  
  //make any messy variable specific adjustments that have to happen before the plot making-------------------------------------
 
  for(unsigned int var=0; var<hData.size(); var++)
    {
      //take care of any and all rebinning first
      if(plots[var].rebinAmt != 0)
        {
          hDataPresel[var]->Rebin(plots[var].rebinAmt);
          hCRYMCPresel[var]->Rebin(plots[var].rebinAmt);
          hGENIEPresel[var]->Rebin(plots[var].rebinAmt);
          hData[var]->Rebin(plots[var].rebinAmt);
          hCRYMC[var]->Rebin(plots[var].rebinAmt);
          hGENIE[var]->Rebin(plots[var].rebinAmt);

          hDataPreselWtd[var]->Rebin(plots[var].rebinAmt);
          hCRYMCPreselWtd[var]->Rebin(plots[var].rebinAmt);
          hDataWtd[var]->Rebin(plots[var].rebinAmt);
          hCRYMCWtd[var]->Rebin(plots[var].rebinAmt);

          hDataPreselEff[var]->Rebin(plots[var].rebinAmt);
          hCRYMCPreselEff[var]->Rebin(plots[var].rebinAmt);
          hDataEff[var]->Rebin(plots[var].rebinAmt);
          hCRYMCEff[var]->Rebin(plots[var].rebinAmt);

          hDataPreselWtdEff[var]->Rebin(plots[var].rebinAmt);
          hCRYMCPreselWtdEff[var]->Rebin(plots[var].rebinAmt);
          hDataWtdEff[var]->Rebin(plots[var].rebinAmt);
          hCRYMCWtdEff[var]->Rebin(plots[var].rebinAmt);
        }
    }

  //even with our angle cut at the previous stage, a few events still snuck into the last bin by being on the bin edge
  //so gonna clean those out
  int lastBin = hData[0]->GetXaxis()->GetLast();
  hData[0]->SetBinContent(lastBin, 0);
  hDataPresel[0]->SetBinContent(lastBin, 0);
  hDataEff[0]->SetBinContent(lastBin, 0);
  hDataWtd[0]->SetBinContent(lastBin, 0);
  hDataPreselWtd[0]->SetBinContent(lastBin, 0);
  hDataWtdEff[0]->SetBinContent(lastBin, 0);

  hData[0]->SetBinError(lastBin, 0);
  hDataPresel[0]->SetBinError(lastBin, 0);
  hDataEff[0]->SetBinError(lastBin, 0);
  hDataWtd[0]->SetBinError(lastBin, 0);
  hDataPreselWtd[0]->SetBinError(lastBin, 0);
  hDataWtdEff[0]->SetBinError(lastBin, 0);

  hCRYMC[0]->GetXaxis()->SetRangeUser(0.4, 1.0);
  hCRYMCPresel[0]->GetXaxis()->SetRangeUser(0.4, 1.0);
  hCRYMCEff[0]->GetXaxis()->SetRangeUser(0.4, 1.0);
  hCRYMCWtd[0]->GetXaxis()->SetRangeUser(0.4, 1.0);
  hCRYMCPreselWtd[0]->GetXaxis()->SetRangeUser(0.4, 1.0);
  hCRYMCWtdEff[0]->GetXaxis()->SetRangeUser(0.4, 1.0);

  hData[0]->GetXaxis()->SetRangeUser(0.4, 1.0);
  hDataPresel[0]->GetXaxis()->SetRangeUser(0.4, 1.0);
  hDataEff[0]->GetXaxis()->SetRangeUser(0.4, 1.0);
  hDataWtd[0]->GetXaxis()->SetRangeUser(0.4, 1.0);
  hDataPreselWtd[0]->GetXaxis()->SetRangeUser(0.4, 1.0);
  hDataWtdEff[0]->GetXaxis()->SetRangeUser(0.4, 1.0);

  hCRYMC[3]->GetXaxis()->SetRangeUser(0, 4.0);
  hCRYMCPresel[3]->GetXaxis()->SetRangeUser(0, 4.0);
  hCRYMCEff[3]->GetXaxis()->SetRangeUser(0, 4.0);
  hCRYMCWtd[3]->GetXaxis()->SetRangeUser(0, 4.0);
  hCRYMCPreselWtd[3]->GetXaxis()->SetRangeUser(0, 4.0);
  hCRYMCWtdEff[3]->GetXaxis()->SetRangeUser(0, 4.0);

  hData[3]->GetXaxis()->SetRangeUser(0, 4.0);
  hDataPresel[3]->GetXaxis()->SetRangeUser(0, 4.0);
  hDataEff[3]->GetXaxis()->SetRangeUser(0, 4.0);
  hDataWtd[3]->GetXaxis()->SetRangeUser(0, 4.0);
  hDataPreselWtd[3]->GetXaxis()->SetRangeUser(0, 4.0);
  hDataWtdEff[3]->GetXaxis()->SetRangeUser(0, 4.0);

  //this is just because the last energy bin only had data events in it
  //so clearing it out makes the eff. plots look nicer
  //as there isn't an eff. for the CRY in this bin
  lastBin = hDataEff[3]->GetXaxis()->GetLast();
  hDataEff[3]->SetBinContent(lastBin, 0);
  hDataWtdEff[3]->SetBinContent(lastBin, 0);
  hDataEff[3]->SetBinError(lastBin, 0);
  hDataWtdEff[3]->SetBinError(lastBin, 0);

  hCRYMC[4]->GetXaxis()->SetRangeUser(0, 0.1);
  hCRYMCPresel[4]->GetXaxis()->SetRangeUser(0, 0.1);
  hCRYMCEff[4]->GetXaxis()->SetRangeUser(0, 0.1);
  hCRYMCWtd[4]->GetXaxis()->SetRangeUser(0, 0.1);
  hCRYMCPreselWtd[4]->GetXaxis()->SetRangeUser(0, 0.1);
  hCRYMCWtdEff[4]->GetXaxis()->SetRangeUser(0, 0.1);
  
  hCRYMC[5]->GetXaxis()->SetRangeUser(0.3, 1);
  hCRYMCPresel[5]->GetXaxis()->SetRangeUser(0.3, 1.0);
  hCRYMCEff[5]->GetXaxis()->SetRangeUser(0.3, 1.0);
  hCRYMCWtd[5]->GetXaxis()->SetRangeUser(0.3, 1.0);
  hCRYMCPreselWtd[5]->GetXaxis()->SetRangeUser(0.3, 1.0);
  hCRYMCWtdEff[5]->GetXaxis()->SetRangeUser(0.3, 1.0);

  hData[4]->GetXaxis()->SetRangeUser(0, 0.1);
  hDataPresel[4]->GetXaxis()->SetRangeUser(0, 0.1);
  hDataEff[4]->GetXaxis()->SetRangeUser(0, 0.1);
  hDataWtd[4]->GetXaxis()->SetRangeUser(0, 0.1);
  hDataPreselWtd[4]->GetXaxis()->SetRangeUser(0, 0.1);
  hDataWtdEff[4]->GetXaxis()->SetRangeUser(0, 0.1);
  
  hData[5]->GetXaxis()->SetRangeUser(0.3, 1);
  hDataPresel[5]->GetXaxis()->SetRangeUser(0.3, 1.0);
  hDataEff[5]->GetXaxis()->SetRangeUser(0.3, 1.0);
  hDataWtd[5]->GetXaxis()->SetRangeUser(0.3, 1.0);
  hDataPreselWtd[5]->GetXaxis()->SetRangeUser(0.3, 1.0);
  hDataWtdEff[5]->GetXaxis()->SetRangeUser(0.3, 1.0);
  
 
  //and let's make some plots!

  char saveName[128];
  for(unsigned int var=0; var < 6; var ++)//var<hData.size(); var++)
    {

       //unweighted distributions--------------------------------------------------
      TCanvas *c1 = new TCanvas("c1", "presel distributions", 800,600);
      hCRYMCPresel[var]->GetXaxis()->SetTitle(plots[var].XaxisLabel);
      makePlot(c1, &*hDataPresel[var], &*hCRYMCPresel[var], &*hGENIEPresel[var], 0, plots[var].legSide, beam=="FHC" ? true : false);
      DrawBeamLabel(beam=="FHC" ? true : false);

      TCanvas *c2 = new TCanvas("c2", "distributions", 800,600);
      hCRYMC[var]->GetXaxis()->SetTitle(plots[var].XaxisLabel);
      makePlot(c2, &*hData[var], &*hCRYMC[var], &*hGENIE[var], 1, plots[var].legSide, beam=="FHC" ? true : false);
      DrawBeamLabel(beam=="FHC" ? true : false);

      if(plots[var].yAxisLimit != 0) {
        hCRYMCPresel[var]->GetYaxis()->SetRangeUser(0, plots[var].yAxisLimit);
        hCRYMC[var]->GetYaxis()->SetRangeUser(0, plots[var].yAxisLimit);
      }

      //weighted distributions-------------------------------------------------------
      TCanvas *c3 = new TCanvas("c3", "wgtd presel distributions", 800,600);
      hCRYMCPreselWtd[var]->GetXaxis()->SetTitle(plots[var].XaxisLabel);
      makeWgtdPlot(c3, &*hDataPreselWtd[var], &*hCRYMCPreselWtd[var], &*hGENIEPresel[var], 0, plots[var].legSide, beam=="FHC" ? true : false);
      DrawBeamLabel(beam=="FHC" ? true : false);

      TCanvas *c4 = new TCanvas("c4", "wgtd distributions", 800,600);
      hCRYMCWtd[var]->GetXaxis()->SetTitle(plots[var].XaxisLabel);
      makeWgtdPlot(c4, &*hDataWtd[var], &*hCRYMCWtd[var], &*hGENIE[var], 1, plots[var].legSide, beam=="FHC" ? true : false);
      DrawBeamLabel(beam=="FHC" ? true : false);

      if(plots[var].yAxisLimit != 0) {
        hCRYMCPreselWtd[var]->GetYaxis()->SetRangeUser(0, plots[var].yAxisLimit);
        hCRYMCWtd[var]->GetYaxis()->SetRangeUser(0, plots[var].yAxisLimit);
      }

      //unweighted efficiencies------------------------------------------------------
      TCanvas *c5 =new TCanvas("c5","Eff",800,800);
      //DrawBeamLabelEff(beam=="FHC" ? true : false);
      hDataEff[var]->GetXaxis()->SetTitle(plots[var].XaxisLabel);

      if (var == 3 && beam != "FHC") makeEffPlot(c5, &*hDataEff[var], &*hDataPreselEff[var], &*hCRYMCEff[var], &*hCRYMCPreselEff[var], 0, false);
      else makeEffPlot(c5, &*hDataEff[var], &*hDataPreselEff[var], &*hCRYMCEff[var], &*hCRYMCPreselEff[var], plots[var].legSide, beam=="FHC" ? true : false);
      hCRYMCEff[3]->GetYaxis()->SetRangeUser(0,1.2);

      //weighted efficiencies-------------------------------------------------------------
      TCanvas *c6 =new TCanvas("c6","Weighted Eff",800,800);
      //DrawBeamLabelEff(beam=="FHC" ? true : false);
      hDataWtdEff[var]->GetXaxis()->SetTitle(plots[var].XaxisLabel);

      if (var == 3 && beam != "FHC") makeWgtdEffPlot(c6, &*hDataWtdEff[var], &*hDataPreselWtdEff[var], &*hCRYMCWtdEff[var], &*hCRYMCPreselWtdEff[var], 0, false);
      else makeWgtdEffPlot(c6, &*hDataWtdEff[var], &*hDataPreselWtdEff[var], &*hCRYMCWtdEff[var], &*hCRYMCPreselWtdEff[var], plots[var].legSide, beam=="FHC" ? true : false);
      hCRYMCWtdEff[3]->GetYaxis()->SetRangeUser(0, 2.0);

      //and repeat the distributions without the GENIE sample
       //unweighted distributions--------------------------------------------------
      TCanvas *c7 = new TCanvas("c7", "presel distributions", 800,600);
      hCRYMCPresel[var]->GetXaxis()->SetTitle(plots[var].XaxisLabel);
      makePlot(c7, &*hDataPresel[var], &*hCRYMCPresel[var], 0, plots[var].legSide);
      DrawBeamLabel(beam=="FHC" ? true : false);

      TCanvas *c8 = new TCanvas("c8", "distributions", 800,600);
      hCRYMC[var]->GetXaxis()->SetTitle(plots[var].XaxisLabel);
      makePlot(c8, &*hData[var], &*hCRYMC[var], 1, plots[var].legSide);
      DrawBeamLabel(beam=="FHC" ? true : false);

      if(plots[var].yAxisLimit != 0) {
        hCRYMCPresel[var]->GetYaxis()->SetRangeUser(0, plots[var].yAxisLimit);
        hCRYMC[var]->GetYaxis()->SetRangeUser(0, plots[var].yAxisLimit);
      }

      //weighted distributions-------------------------------------------------------
      TCanvas *c9 = new TCanvas("c9", "wgtd presel distributions", 800,600);
      hCRYMCPreselWtd[var]->GetXaxis()->SetTitle(plots[var].XaxisLabel);
      makeWgtdPlot(c9, &*hDataPreselWtd[var], &*hCRYMCPreselWtd[var], 0, plots[var].legSide);
      DrawBeamLabel(beam=="FHC" ? true : false);

      TCanvas *c10 = new TCanvas("c10", "wgtd distributions", 800,600);
      hCRYMCWtd[var]->GetXaxis()->SetTitle(plots[var].XaxisLabel);
      makeWgtdPlot(c10, &*hDataWtd[var], &*hCRYMCWtd[var], 1, plots[var].legSide);
      DrawBeamLabel(beam=="FHC" ? true : false);

      if(plots[var].yAxisLimit != 0) {
        hCRYMCPreselWtd[var]->GetYaxis()->SetRangeUser(0, plots[var].yAxisLimit);
        hCRYMCWtd[var]->GetYaxis()->SetRangeUser(0, plots[var].yAxisLimit);
      }


      //and save all the canvases------------------------------------------------------------
      sprintf(saveName, "EffPlots/%s/%s_Presel_%s.png", beam.c_str(), plots[var].savename, beam.c_str());
      c1->SaveAs(saveName);
      sprintf(saveName, "EffPlots/%s/%s_%s.png", beam.c_str(), plots[var].savename, beam.c_str());
      c2->SaveAs(saveName);

      sprintf(saveName, "EffPlots/%s/%s_Presel_Wtd_%s.png", beam.c_str(), plots[var].savename, beam.c_str());
      c3->SaveAs(saveName);
      sprintf(saveName, "EffPlots/%s/%s_Wtd_%s.png", beam.c_str(), plots[var].savename, beam.c_str());
      c4->SaveAs(saveName);

      sprintf(saveName, "EffPlots/%s/%s_Eff_%s.png", beam.c_str(), plots[var].savename, beam.c_str());
      c5->SaveAs(saveName);
      sprintf(saveName, "EffPlots/%s/%s_Eff_Wtd_%s.png", beam.c_str(), plots[var].savename, beam.c_str());
      c6->SaveAs(saveName);

      sprintf(saveName, "EffPlots/%s/%s_Presel_%s_noGENIE.png", beam.c_str(), plots[var].savename, beam.c_str());
      c7->SaveAs(saveName);
      sprintf(saveName, "EffPlots/%s/%s_%s_noGENIE.png", beam.c_str(), plots[var].savename, beam.c_str());
      c8->SaveAs(saveName);

      sprintf(saveName, "EffPlots/%s/%s_Presel_Wtd_%s_noGENIE.png", beam.c_str(), plots[var].savename, beam.c_str());
      c9->SaveAs(saveName);
      sprintf(saveName, "EffPlots/%s/%s_Wtd_%s_noGENIE.png", beam.c_str(), plots[var].savename, beam.c_str());
      c10->SaveAs(saveName);

      // make data/CRY comparison for the CVN spectra
      if(var == 5){
        //unweighted
        TCanvas *c11 =new TCanvas("c11","Eff",800,800); 
        hDataPresel[var]->GetXaxis()->SetTitle(plots[var].XaxisLabel);
        makeDataCryRatio(c11, &*hDataPresel[var], &*hCRYMCPresel[var], 0, plots[var].legSide, beam=="FHC" ? true : false, false);
        sprintf(saveName, "EffPlots/%s/%s_%s_Ratio.png", beam.c_str(), plots[var].savename, beam.c_str());
        c11->SaveAs(saveName);

        //weighted
        TCanvas *c12 =new TCanvas("c12","Eff",800,800); 
        hDataPreselWtd[var]->GetXaxis()->SetTitle(plots[var].XaxisLabel);
        makeDataCryRatio(c12, &*hDataPreselWtd[var], &*hCRYMCPreselWtd[var], 0, plots[var].legSide, beam=="FHC" ? true : false, true);
        sprintf(saveName, "EffPlots/%s/%s_Wtd_%s_Ratio.png", beam.c_str(), plots[var].savename, beam.c_str());
        c12->SaveAs(saveName);
      }
      //------------
      //and make the efficiency plots without the ratio on the bottom
      //unweighted efficiencies------------------------------------------------------
      TCanvas *c13 =new TCanvas("c13","Eff",800,600);
      hCRYMCEff[var]->GetXaxis()->SetTitle(plots[var].XaxisLabel);
      if (var == 3 && beam != "FHC") makeEffPlotNoRatio(c13, &*hDataEff[var], &*hCRYMCEff[var], 0, false, false);
      else makeEffPlotNoRatio(c13, &*hDataEff[var], &*hCRYMCEff[var], plots[var].legSide, beam=="FHC" ? true : false, false);
      hCRYMCEff[3]->GetYaxis()->SetRangeUser(0,1.2);

      //weighted efficiencies-------------------------------------------------------------
      TCanvas *c14 =new TCanvas("c14","Weighted Eff",800,600);
      hCRYMCWtdEff[var]->GetXaxis()->SetTitle(plots[var].XaxisLabel);
      if (var == 3 && beam != "FHC") makeEffPlotNoRatio(c14, &*hDataWtdEff[var], &*hCRYMCWtdEff[var], 0, false, true);
      else makeEffPlotNoRatio(c14, &*hDataWtdEff[var], &*hCRYMCWtdEff[var], plots[var].legSide, beam=="FHC" ? true : false, true);
      //hCRYMCWtdEff[3]->GetYaxis()->SetRangeUser(0, 2.0);

      sprintf(saveName, "EffPlots/%s/%s_Eff_%s_noRatio.png", beam.c_str(), plots[var].savename, beam.c_str());
      c13->SaveAs(saveName);
      sprintf(saveName, "EffPlots/%s/%s_Eff_Wtd_%s_noRatio.png", beam.c_str(), plots[var].savename, beam.c_str());
      c14->SaveAs(saveName);

      
    }


}


void makePlot(TCanvas *c1, TH1 *hCosmicsData, TH1 *hCRYMC, TH1 *hGENIEFluxswap, int sel, int legSide, bool isFHC)
{
  c1->cd();

  hCRYMC->GetXaxis()->CenterTitle();
  hCRYMC->GetYaxis()->CenterTitle();
  hCosmicsData->SetLineColor(kBlack);
  hCosmicsData->SetMarkerStyle(20);
  hCRYMC->SetFillColorAlpha(kMagenta+2, 0.5);
  hCRYMC->SetLineColor(kMagenta+2);
  hGENIEFluxswap->SetLineColor(kBlue);

  double kNormFactor = hCosmicsData->Integral();
  hGENIEFluxswap->Scale(kNormFactor/hGENIEFluxswap->Integral());
  hCRYMC->Scale(kNormFactor/hCRYMC->Integral());

  hCRYMC->GetYaxis()->SetTitle("Events");
  hCRYMC->Draw("E2");
  //this copy is literally just to be able to have the horizontal line for each bin
  auto hCRYMC2 = (TH1*)hCRYMC->Clone();
  for(int bin =0; bin <= hCRYMC2->GetXaxis()->GetLast(); bin ++) hCRYMC2->SetBinError(bin, 0.000001);
  hCRYMC2->SetFillColor(kWhite);
  hCRYMC2->SetMarkerColor(kMagenta+2);
  hCRYMC2->Draw("E hist same");
  hCosmicsData->Draw("E1 same");
  hGENIEFluxswap->Draw("hist E1 same");
  hCRYMC->GetYaxis()->SetLabelSize(0.037);
  hCRYMC->GetYaxis()->SetTitleOffset(0.95);

  TLegend * leg;
  if (legSide ==0) {
    leg = new TLegend(0.15,0.62,0.75,0.85);
    leg->SetNColumns(2);
    leg->AddEntry(hCosmicsData,"Cosmic Data","p");
    leg->AddEntry((TObject*)0, sel==0 ? "Preselected DiF Events": "Selected DiF Events", "");
    leg->AddEntry(hCRYMC,"Cosmics MC (CRY)","l");
    leg->AddEntry((TObject*)0, "Area norm. to data", "");
    leg->AddEntry(hCRYMC, "Cosmics MC Stat. Err.", "f");
    leg->AddEntry((TObject*)0, "","");
    leg->AddEntry(hGENIEFluxswap, isFHC ? "#nu_{e} Signal MC" : "#bar{#nu}_{e} Signal MC","l");
  }
  else {
    leg = new TLegend(0.15,0.62,0.85,0.85);
    leg->SetNColumns(2);
    leg->AddEntry((TObject*)0, sel==0 ? "Preselected DiF Events": "Selected DiF Events", "");
    leg->AddEntry(hCosmicsData,"Cosmic Data","p");
    leg->AddEntry((TObject*)0, "Area norm. to data", "");
    leg->AddEntry(hCRYMC,"Cosmics MC (CRY)","l");
    leg->AddEntry((TObject*)0, "","");
    leg->AddEntry(hCRYMC, "Cosmics MC Stat. Err.", "f");
    leg->AddEntry((TObject*)0, "","");
    leg->AddEntry(hGENIEFluxswap, isFHC ? "#nu_{e} Signal MC" : "#bar{#nu}_{e} Signal MC","l");
  }
  leg->SetFillStyle(0);
  leg->SetTextSize(0.04);
  leg->Draw("same");
  myPreliminary();
}

void makePlot(TCanvas *c1, TH1 *hCosmicsData, TH1 *hCRYMC, int sel, int legSide)
{
  //don't plot the GENIE histo
  c1->cd();

  hCRYMC->GetXaxis()->CenterTitle();
  hCRYMC->GetYaxis()->CenterTitle();
  hCosmicsData->SetLineColor(kBlack);
  hCosmicsData->SetMarkerStyle(20);
  hCRYMC->SetLineColor(kMagenta+2);
  hCRYMC->SetFillColorAlpha(kMagenta+2, 0.5);

  double kNormFactor = hCosmicsData->Integral();
  hCRYMC->Scale(kNormFactor/hCRYMC->Integral());

  hCRYMC->GetYaxis()->SetTitle("Events");
  hCRYMC->Draw("E2");
  //this copy is literally just to be able to have the horizontal line for each bin
  //without the vertical lines
  auto hCRYMC2 = (TH1*)hCRYMC->Clone();
  for(int bin =0; bin <= hCRYMC2->GetXaxis()->GetLast(); bin ++) hCRYMC2->SetBinError(bin, 0.000001);
  hCRYMC2->SetFillColor(kWhite);
  hCRYMC2->SetMarkerColor(kMagenta+2);
  hCRYMC2->Draw("E hist same");
  hCosmicsData->Draw("E1 same");
  hCRYMC->GetYaxis()->SetLabelSize(0.037);
  //hCRYMC->GetYaxis()->SetTitleOffset(0.95);

  TLegend * leg;
  if (legSide ==0) {
    leg = new TLegend(0.15,0.68,0.75,0.85);
    leg->SetNColumns(2);
    leg->AddEntry(hCosmicsData,"Cosmic Data","p");
    leg->AddEntry((TObject*)0, sel==0 ? "Preselected DiF Events": "Selected DiF Events", "");
    leg->AddEntry(hCRYMC,"Cosmics MC (CRY)","l");
    leg->AddEntry((TObject*)0, "Area norm. to data", "");
    leg->AddEntry(hCRYMC, "Cosmics MC Stat. Err.", "f");
  }
  else {
    leg = new TLegend(0.15,0.68,0.85,0.85);
    leg->SetNColumns(2);
    leg->AddEntry((TObject*)0, sel==0 ? "Preselected DiF Events": "Selected DiF Events", "");
    leg->AddEntry(hCosmicsData,"Cosmic Data","p");
    leg->AddEntry((TObject*)0, "Area norm. to data", "");
    leg->AddEntry(hCRYMC,"Cosmics MC (CRY)","l");
    leg->AddEntry((TObject*)0, "","");
    leg->AddEntry(hCRYMC, "Cosmics MC Stat. Err.", "f");
  }
  leg->SetFillStyle(0);
  leg->SetTextSize(0.04);
  leg->Draw("same");
  myPreliminary();
}


void makeWgtdPlot(TCanvas *c1, TH1 *hCosmicsData, TH1 *hCRYMC, TH1 *hGENIEFluxswap, int sel, int legSide, bool isFHC)
{
  c1->cd();

  hCRYMC->GetXaxis()->CenterTitle();
  hCRYMC->GetYaxis()->CenterTitle();
  hCosmicsData->SetLineColor(kBlack);
  hCosmicsData->SetMarkerStyle(20);
  hCRYMC->SetLineColor(kMagenta+2);
  hCRYMC->SetFillColorAlpha(kMagenta+2, 0.5);
  hGENIEFluxswap->SetLineColor(kBlue);

  double kNormFactor = hCosmicsData->Integral();
  hGENIEFluxswap->Scale(kNormFactor/hGENIEFluxswap->Integral());
  hCRYMC->Scale(kNormFactor/hCRYMC->Integral());

  hCRYMC->GetYaxis()->SetTitle("Events (cos#theta_{beam} weighted)");
  hCRYMC->Draw("E2");
  //this copy is literally just to be able to have the horizontal line for each bin
  //without the vertical lines
  auto hCRYMC2 = (TH1*)hCRYMC->Clone();
  for(int bin =0; bin <= hCRYMC2->GetXaxis()->GetLast(); bin ++) hCRYMC2->SetBinError(bin, 0.000001);
  hCRYMC2->SetFillColor(kWhite);
  hCRYMC2->SetMarkerColor(kMagenta+2);
  hCRYMC2->Draw("E hist same");
  hCosmicsData->Draw("E1 same");  
  //hCRYMC->Draw("E2 same");
  hGENIEFluxswap->Draw("hist E1 same");
  hCRYMC->GetYaxis()->SetLabelSize(0.037);
  hCRYMC->GetYaxis()->SetTitleOffset(1.0);
  hCRYMC->GetYaxis()->SetTitleSize(0.05);

  TLegend * leg;
  if (legSide ==0) {
    leg = new TLegend(0.15,0.62,0.75,0.85);
    leg->SetNColumns(2);
    leg->AddEntry(hCosmicsData,"Cosmic Data","p");
    leg->AddEntry((TObject*)0, sel==0 ? "Preselected DiF Events": "Selected DiF Events", "");
    leg->AddEntry(hCRYMC,"Cosmics MC (CRY)","l");
    leg->AddEntry((TObject*)0, "Area norm. to data", "");
    leg->AddEntry(hCRYMC, "Cosmics MC Stat. Err.", "f");
    leg->AddEntry((TObject*)0, "","");
    leg->AddEntry(hGENIEFluxswap, isFHC ? "#nu_{e} Signal MC" : "#bar{#nu}_{e} Signal MC","l");
  }
  else {
    leg = new TLegend(0.15,0.62,0.85,0.85);
    leg->SetNColumns(2);
    leg->AddEntry((TObject*)0, sel==0 ? "Preselected DiF Events": "Selected DiF Events", "");
    leg->AddEntry(hCosmicsData,"Cosmic Data","p");
    leg->AddEntry((TObject*)0, "Area norm. to data", "");
    leg->AddEntry(hCRYMC,"Cosmics MC (CRY)","l");
    leg->AddEntry((TObject*)0, "","");
    leg->AddEntry(hCRYMC, "Cosmics MC Stat. Err.", "f");
    leg->AddEntry((TObject*)0, "","");
    leg->AddEntry(hGENIEFluxswap, isFHC ? "#nu_{e} Signal MC" : "#bar{#nu}_{e} Signal MC","l");
  }
  leg->SetFillStyle(0);
  leg->SetTextSize(0.04);
  leg->Draw("same");
  myPreliminary();
}

void makeWgtdPlot(TCanvas *c1, TH1 *hCosmicsData, TH1 *hCRYMC, int sel, int legSide)
{
  //don't plot the GENIE histo
  c1->cd();

  hCRYMC->GetXaxis()->CenterTitle();
  hCRYMC->GetYaxis()->CenterTitle();
  hCosmicsData->SetLineColor(kBlack);
  hCosmicsData->SetMarkerStyle(20);
  hCRYMC->SetLineColor(kMagenta+2);
  hCRYMC->SetFillColorAlpha(kMagenta+2, 0.5);

  double kNormFactor = hCosmicsData->Integral();
  hCRYMC->Scale(kNormFactor/hCRYMC->Integral());
  hCRYMC->GetYaxis()->SetTitle("Events (cos#theta_{beam} weighted)");
  hCRYMC->Draw("E2");
  //this copy is literally just to be able to have the horizontal line for each bin
  //without the vertical lines
  auto hCRYMC2 = (TH1*)hCRYMC->Clone();
  for(int bin =0; bin <= hCRYMC2->GetXaxis()->GetLast(); bin ++) hCRYMC2->SetBinError(bin, 0.000001);
  hCRYMC2->SetFillColor(kWhite);
  hCRYMC2->SetMarkerColor(kMagenta+2);
  hCRYMC2->Draw("E hist same");
  hCosmicsData->Draw("E1 same");
  //hCRYMC->GetYaxis()->SetLabelSize(0.037);
  //hCRYMC->GetYaxis()->SetTitleOffset(0.95);

  TLegend * leg;
  if (legSide ==0) {
    leg = new TLegend(0.15,0.68,0.75,0.85);
    leg->SetNColumns(2);
    leg->AddEntry(hCosmicsData,"Cosmic Data","p");
    leg->AddEntry((TObject*)0, sel==0 ? "Preselected DiF Events": "Selected DiF Events", "");
    leg->AddEntry(hCRYMC,"Cosmics MC (CRY)","l");
    leg->AddEntry((TObject*)0, "Area norm. to data", "");
    leg->AddEntry(hCRYMC, "Cosmics MC Stat. Err.", "f");
  }
  else {
    leg = new TLegend(0.15,0.68,0.85,0.85);
    leg->SetNColumns(2);
    leg->AddEntry((TObject*)0, sel==0 ? "Preselected DiF Events": "Selected DiF Events", "");
    leg->AddEntry(hCosmicsData,"Cosmic Data","p");
    leg->AddEntry((TObject*)0, "Area norm. to data", "");
    leg->AddEntry(hCRYMC,"Cosmics MC (CRY)","l");
    leg->AddEntry((TObject*)0, "","");
    leg->AddEntry(hCRYMC, "Cosmics MC Stat. Err.", "f");
  }
  leg->SetFillStyle(0);
  leg->SetTextSize(0.04);
  leg->Draw("same");

  myPreliminary();
}

void makeDataCryRatio(TCanvas *c1, TH1 *hCosmicsData, TH1 *hCRYMC, int sel, int legSide, bool isFHC, bool isWgtd)
{
  c1->cd();
  TPad * pad1 = new TPad("pad1", "pad1", 0, 0.375, 1, 1.0);
  pad1->Draw();
  pad1->SetBottomMargin(0.02);
  pad1->cd();

  hCRYMC->GetXaxis()->CenterTitle();
  hCRYMC->GetYaxis()->CenterTitle();
  hCosmicsData->SetLineColor(kBlack);
  hCosmicsData->SetMarkerStyle(20);
  hCRYMC->SetLineColor(kMagenta+2);
  hCRYMC->SetFillColorAlpha(kMagenta+2, 0.5);

  double kNormFactor = hCosmicsData->Integral();
  hCRYMC->Scale(kNormFactor/hCRYMC->Integral());

  hCRYMC->GetYaxis()->SetTitle(isWgtd ? "Events (cos#theta_{beam} weighted)" : "Events");
  hCRYMC->Draw("E2");
  //this copy is literally just to be able to have the horizontal line for each bin
  //without the vertical lines
  auto hCRYMC2 = (TH1*)hCRYMC->Clone();
  for(int bin =0; bin <= hCRYMC2->GetXaxis()->GetLast(); bin ++) hCRYMC2->SetBinError(bin, 0.000001);
  hCRYMC2->SetFillColor(kWhite);
  hCRYMC2->SetMarkerColor(kMagenta+2);
  hCRYMC2->Draw("E hist same");
  hCosmicsData->Draw("E1 same");
  hCRYMC->GetYaxis()->SetLabelSize(0.037);
  hCRYMC->GetYaxis()->SetTitleOffset(0.95);
  hCRYMC->GetXaxis()->SetLabelSize(0);

  TLegend * leg;
  if (legSide ==0) leg = new TLegend(0.15,0.68,0.75,0.85);
  else leg = new TLegend(0.5,0.68,0.8,0.85);
  leg->SetNColumns(2);
  leg->AddEntry(hCosmicsData,"Cosmic Data","p");
  leg->AddEntry((TObject*)0, sel==0 ? "Preselected DiF Events": "Selected DiF Events", "");
  leg->AddEntry(hCRYMC,"Cosmics MC (CRY)","l");
  leg->AddEntry((TObject*)0, "Area norm. to data", "");
  leg->AddEntry(hCRYMC, "Cosmics MC Stat. Err.", "f");
  leg->SetTextSize(0.05);
  leg->Draw("same");
  myPreliminary();

  c1->cd();
  TPad * pad2 = new TPad("pad2", "pad2", 0, 0.002, 1, 0.375);
  pad2->SetTopMargin(0);
  pad2->SetBottomMargin(0.2);
  pad2->Draw();
  pad2->cd();
  pad2->SetGridy();

  TH1D* hRatio =(TH1D*)hCosmicsData->Clone("hRatio");

  hRatio->GetXaxis()->CenterTitle();
  hRatio->GetYaxis()->CenterTitle();
  hRatio->SetLineColor(kBlack);
  hRatio->SetMarkerStyle(20);
  hRatio->SetLineWidth(3);
  hRatio->Divide(hCRYMC);
  hRatio->Draw("E1");
  hRatio->GetXaxis()->SetLabelSize(0.075);
  hRatio->GetYaxis()->SetLabelSize(0.075);
  hRatio->GetXaxis()->SetTitleOffset(0.85);
  hRatio->GetYaxis()->SetTitleOffset(0.55);
  hRatio->GetXaxis()->SetTitleSize(0.09);
  hRatio->GetYaxis()->SetTitleSize(0.085);
  hRatio->GetYaxis()->SetRangeUser(0.,2.0);
  hRatio->GetYaxis()->SetTitle("#frac{Data}{MC}");

  hCRYMC->GetXaxis()->SetTitle("");
  c1->cd();
  DrawBeamLabelEff(isFHC);
}



void makeEffPlot(TCanvas *c1, TH1 *hCosmicsDataCore, TH1 *hCosmicsDataPresel, TH1 *hCRYMCCore, TH1 *hCRYMCPresel, int legSide, bool isFHC)
{
  c1->cd();
  TPad * pad1 = new TPad("pad1", "pad1", 0, 0.375, 1, 1.0);
  pad1->Draw();
  pad1->SetBottomMargin(0.02);
  pad1->cd();

  hCRYMCCore->GetXaxis()->CenterTitle();
  hCRYMCCore->GetYaxis()->CenterTitle();
  hCosmicsDataCore->SetLineColor(kBlack);
  hCosmicsDataCore->SetMarkerStyle(20);
  hCosmicsDataCore->SetLineWidth(2);
  hCRYMCCore->SetLineColor(kMagenta+2);
  hCRYMCCore->SetFillColorAlpha(kMagenta+2, 0.5);
  hCRYMCCore->SetLineWidth(2);  

  hCosmicsDataCore->Divide(hCosmicsDataPresel);
  hCRYMCCore->Divide(hCRYMCPresel);
  hCRYMCCore->Draw("E2");
  //this copy is literally just to be able to have the horizontal line for each bin
  //without the vertical lines
  auto hCRYMC2 = (TH1*)hCRYMCCore->Clone();
  for(int bin =0; bin <= hCRYMC2->GetXaxis()->GetLast(); bin ++) hCRYMC2->SetBinError(bin, 0.000001);
  hCRYMC2->SetFillColor(kWhite);
  hCRYMC2->SetMarkerColor(kMagenta+2);
  hCRYMC2->Draw("E hist same");
  hCosmicsDataCore->Draw("E1 same");
  
  //hCRYMCCore->GetYaxis()->SetLabelSize(0.037);
  hCRYMCCore->GetXaxis()->SetLabelSize(0);
  hCRYMCCore->GetYaxis()->SetTitle("Selection Efficiency");

  //hCRYMCCore->GetYaxis()->SetTitleOffset(0.95);
  hCRYMCCore->GetYaxis()->SetRangeUser(0,1.2*gPad->GetUymax());
  
  TLegend * leg;
  if (legSide ==0) leg = new TLegend(0.15,0.68,0.45,0.85);
  else leg = new TLegend(0.5,0.68,0.8,0.85);
  leg->AddEntry(hCosmicsDataCore,"Cosmic Data","p");
  leg->AddEntry(hCRYMCCore,"Cosmic MC (CRY)","l");
  leg->AddEntry(hCRYMCCore, "Cosmic MC Stat. Err.", "f");
  leg->SetTextSize(0.05);
  leg->Draw("same");
  myPreliminary();

  //double dif_label_y =  isWgtd ? 1.6 : 1.0;
  //double dif_label_y = hCRYMCEff->GetYaxis()->GetBinUpEdge(hCRYMCEff->GetYaxis()->GetLast());
  double dif_label_x = 0.5;//hCRYMCCore->GetXaxis()->GetBinLowEdge(2);
  TText * tt = new TText(0.9*dif_label_x, 0.6, "FD Decay-in-Flight");
  tt->SetTextSize(0.06);
  tt->Draw();

  //DrawBeamLabelEff(isFHC);

  c1->cd();
  TPad * pad2 = new TPad("pad2", "pad2", 0, 0.002, 1, 0.375);
  pad2->SetTopMargin(0);
  pad2->SetBottomMargin(0.2);
  pad2->Draw();
  pad2->cd();
  pad2->SetGridy();

  TH1D* hEffRatio =(TH1D*)hCosmicsDataCore->Clone("hEffRatio");

  hEffRatio->GetXaxis()->CenterTitle();
  hEffRatio->GetYaxis()->CenterTitle();
  hEffRatio->SetLineColor(kBlack);
  hEffRatio->SetMarkerStyle(20);
  hEffRatio->SetLineWidth(3);
  hEffRatio->Divide(hCRYMCCore);
  hEffRatio->Draw("E1");
  hEffRatio->GetXaxis()->SetLabelSize(0.075);
  hEffRatio->GetYaxis()->SetLabelSize(0.075);
  hEffRatio->GetXaxis()->SetTitleOffset(0.85);
  hEffRatio->GetYaxis()->SetTitleOffset(0.55);
  hEffRatio->GetXaxis()->SetTitleSize(0.09);
  hEffRatio->GetYaxis()->SetTitleSize(0.085);
  hEffRatio->GetYaxis()->SetRangeUser(0.5,1.5);
  hEffRatio->GetYaxis()->SetTitle("#frac{Data}{MC}");

  hCRYMCCore->GetXaxis()->SetTitle("");
  c1->cd();
  DrawBeamLabelEff(isFHC);
}


void makeWgtdEffPlot(TCanvas *c1, TH1 *hCosmicsDataCore, TH1 *hCosmicsDataPresel, TH1 *hCRYMCCore, TH1 *hCRYMCPresel, int legSide, bool isFHC)
{
  c1->cd();
  TPad * pad1 = new TPad("pad1", "pad1", 0, 0.375, 1, 1.0);
  pad1->Draw();
  pad1->SetBottomMargin(0.02);
  pad1->cd();

  hCRYMCCore->GetXaxis()->CenterTitle();
  hCRYMCCore->GetYaxis()->CenterTitle();
  hCosmicsDataCore->SetLineColor(kBlack);
  hCosmicsDataCore->SetMarkerStyle(20);
  hCosmicsDataCore->SetLineWidth(2);
  hCRYMCCore->SetLineColor(kMagenta+2);
  hCRYMCCore->SetFillColorAlpha(kMagenta+2, 0.5);
  hCRYMCCore->SetLineWidth(2);  

  hCosmicsDataCore->Divide(hCosmicsDataPresel);
  hCRYMCCore->Divide(hCRYMCPresel);
  hCRYMCCore->Draw("E2");
  //this copy is literally just to be able to have the horizontal line for each bin
  //without the vertical lines
  auto hCRYMC2 = (TH1*)hCRYMCCore->Clone();
  for(int bin =0; bin <= hCRYMC2->GetXaxis()->GetLast(); bin ++) hCRYMC2->SetBinError(bin, 0.000001);
  hCRYMC2->SetFillColor(kWhite);
  hCRYMC2->SetMarkerColor(kMagenta+2);
  hCRYMC2->SetLineWidth(2);
  hCRYMC2->Draw("E hist same");
  hCosmicsDataCore->Draw("E1 same");
  //hCRYMCCore->Draw("hist E1 same");
  
  hCRYMCCore->GetYaxis()->SetLabelSize(0.037);
  hCRYMCCore->GetXaxis()->SetLabelSize(0);
  hCRYMCCore->GetYaxis()->SetTitle("cos#theta_{beam} Weighted Efficiency");

  hCRYMCCore->GetYaxis()->SetRangeUser(0,1.2);
  
  TLegend * leg;
  if (legSide ==0) leg = new TLegend(0.15,0.68,0.45,0.85);
  else leg = new TLegend(0.5,0.68,0.8,0.85);
  leg->AddEntry(hCosmicsDataCore,"Cosmic Data","p");
  leg->AddEntry(hCRYMCCore,"Cosmic MC (CRY)","l");
  leg->AddEntry(hCRYMCCore, "Cosmic MC Stat. Err.", "f");
  leg->SetTextSize(0.05);
  leg->Draw("same");
  myPreliminary();

  double dif_label_y = 0.6;
  //double dif_label_y =  isWgtd ? 1.6 : 1.0;
  //double dif_label_y = hCRYMCEff->GetYaxis()->GetBinUpEdge(hCRYMCEff->GetYaxis()->GetLast());
  double dif_label_x = 0.5;//hCRYMCCore->GetXaxis()->GetBinLowEdge(2);
  TText * tt = new TText(0.9*dif_label_x, dif_label_y, "FD Decay-in-Flight");
  tt->SetTextSize(0.06);
  tt->Draw();

  c1->cd();
  TPad * pad2 = new TPad("pad2", "pad2", 0, 0.002, 1, 0.375);
  pad2->SetTopMargin(0);
  pad2->SetBottomMargin(0.2);
  pad2->Draw();
  pad2->cd();
  pad2->SetGridy();

  TH1D* hEffRatio =(TH1D*)hCosmicsDataCore->Clone("hEffRatio");

  hEffRatio->GetXaxis()->CenterTitle();
  hEffRatio->GetYaxis()->CenterTitle();
  hEffRatio->SetLineColor(kBlack);
  hEffRatio->SetMarkerStyle(20);
  hEffRatio->SetLineWidth(3);
  hEffRatio->Divide(hCRYMCCore);
  hEffRatio->Draw("E1");
  hEffRatio->GetXaxis()->SetLabelSize(0.075);
  hEffRatio->GetYaxis()->SetLabelSize(0.075);
  hEffRatio->GetXaxis()->SetTitleOffset(0.85);
  hEffRatio->GetYaxis()->SetTitleOffset(0.55);
  hEffRatio->GetXaxis()->SetTitleSize(0.09);
  hEffRatio->GetYaxis()->SetTitleSize(0.085);
  hEffRatio->GetYaxis()->SetRangeUser(0.5,1.5);
  hEffRatio->GetYaxis()->SetTitle("#frac{Data}{MC}");

  hCRYMCCore->GetXaxis()->SetTitle("");
  c1->cd();
  DrawBeamLabelEff(isFHC);
}


void makeEffPlotNoRatio(TCanvas *c1,  TH1 *hCosmicsDataEff, TH1 *hCRYMCEff, int legSide, bool isFHC, bool isWgtd)
{
  //should only be used after the regular eff plot functions have been used

  hCRYMCEff->GetXaxis()->SetLabelSize(0.04);
  hCRYMCEff->Draw("E2");

  //this copy is literally just to be able to have the horizontal line for each bin
  //without the vertical lines
  auto hCRYMC2 = (TH1*)hCRYMCEff->Clone();
  for(int bin =0; bin <= hCRYMC2->GetXaxis()->GetLast(); bin ++) hCRYMC2->SetBinError(bin, 0.000001);
  hCRYMC2->SetFillColor(kWhite);
  hCRYMC2->SetMarkerColor(kMagenta+2);
  hCRYMC2->Draw("E hist same");
  hCosmicsDataEff->Draw("E1 same");

  TLegend * leg;
  if (legSide ==0) leg = new TLegend(0.15,0.68,0.45,0.85);
  else leg = new TLegend(0.5,0.68,0.8,0.85);
  leg->SetTextSize(0.05);
  leg->AddEntry(hCosmicsDataEff,"Cosmic Data","p");
  leg->AddEntry(hCRYMCEff,"Cosmic MC (CRY)","l");
  leg->AddEntry(hCRYMCEff, "Cosmic MC Stat. Err.", "f");
  leg->Draw("same");
  myPreliminary();
  DrawBeamLabel(isFHC);

  double dif_label_y = 0.6;
  double dif_label_x = 0.25;
  TText * tt = new TText(dif_label_x, dif_label_y, "FD Decay-in-Flight");
  tt->SetTextSize(0.06);
  tt->Draw();

  gPad->RedrawAxis();
}