plot_ehade_quantbound.C

Go to the documentation of this file.

#include "includes.h"

using namespace ana;


void plot_ehade_quantbound(bool fhc = true){


  const int nBins= 19;
  const int nQuant= 4;

  std::string horn;
  std::string name;

  //  int color_pal = kCool;
  int color_pal = kFuchsia;
  Style_t line = kSolid;
  Color_t color;


  if(fhc){
    horn = "fhc", name = "FHC";
    color = kBlue-7;
  }
  if(!fhc){
    horn = "rhc", name = "RHC";
    color = kMagenta+1;
  }

  std::string filename = "/nova/ana/nu_mu_ana/Ana2018/Quantiles/quantiles__" + horn + "_full__numu2018.root";
  TFile* infile = TFile::Open(pnfs2xrootd(filename).c_str());
  TH2* quanthist = (TH2*)infile->Get("FDSpec2D");
  std::vector< std::vector<double> > QuantBounds = GetQuantileBins(quanthist, kNumuCCOptimisedAxis, kHadEFracAxis, nQuant, 1);


  std::vector< TH1* > hist_bound;
  for(int thisQuant = 0; thisQuant < nQuant; thisQuant++){
    TH1* h_temp = quanthist->ProjectionX();
    h_temp -> Clear();
    for(int bin = 1; bin <= nBins; bin++){
      double bound = QuantBounds[bin][thisQuant];
      h_temp -> SetBinContent(bin, bound);
    }
    hist_bound.push_back(h_temp);
  }


  gStyle->SetPalette(color_pal);
  gStyle->SetPadRightMargin(0.11);

  new TCanvas("canvas");
  gPad->SetFillStyle(0);
  quanthist -> Scale(0.1,"width");
  quanthist -> Draw("colz");
  quanthist->GetYaxis()->SetTitle("E_{had} / E_{#nu}");
  quanthist->GetXaxis()->CenterTitle();
  quanthist->GetYaxis()->CenterTitle();
  for(int thisQuant = 1; thisQuant < nQuant; thisQuant++){
    hist_bound[thisQuant] -> SetLineColor(color);
    hist_bound[thisQuant] -> Draw("histsame][");
  }
  TLegend *legend = new TLegend(0.15,0.70,0.30,0.85);
  legend->AddEntry(hist_bound[1], Form("%s", name.c_str()),"l");
  legend->SetTextSize(0.05); //no border for legend                                                                                                                                
  legend->SetBorderSize(0); //no border for legend                                                                                                                                 
  legend->SetFillStyle(0);  //fill colour is transparent                                                                                                                           
  legend->Draw();
  gPad->Update();

  TString outname = "quantiles__" + horn + "_full__numu2018";
  gPad-> Print(outname + ".pdf");
  gPad-> Print(outname + ".png");

}