neutKEsyst_plot.C File Reference

#include "TGaxis.h"
#include "TCanvas.h"
#include "TFile.h"
#include "TH1D.h"
#include "TLatex.h"
#include "TLegend.h"
#include "TStyle.h"
#include "TSystem.h"
#include "TF1.h"
#include "TH2.h"
#include "TGraph.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)
void	NormByBinWidth (TH1D *hist)
void	Print1D (TH1D *hist, std::string plotName)
void	MakeComparisonPlots (std::vector< TH1D * > TheHists, std::string plotName, std::vector< std::string > Labels, bool var)
void	MakeComparisonPlotsWithRatio (std::vector< TH1D * > TheHists, std::string plotName, std::vector< std::string > Labels, bool var)
TH1D *	MakeRatio (TH1D *num, TH1D *denom, int Col, std::string FType, bool var)
void	FitGauss (TH1D *h)
void	neutKEsyst_plot ()

Variables
const std::vector< std::string >	GENIEStr = { "All" }
const unsigned int	nGENIE = 1
std::vector< std::string >	ND_CutNames = {"ND_2018full"}
const unsigned int	nND_CutNames = 1
unsigned int	FCut_ND = nND_CutNames - 1
std::vector< std::string >	FD_CutNames = {"FD_2018full"}
const unsigned int	nFD_CutNames = 1
unsigned int	FCut_FD = nFD_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< double >	POT = {kAna2018FHCPOT, kAna2018RHCPOT}
std::vector< std::string >	DetNames = {"ND", "FD"}
const unsigned int	nDet = 2

Function Documentation

void FitGauss

(

TH1D *

h

)

Definition at line 475 of file neutKEsyst_plot.C.

References om::cout, allTimeWatchdog::endl, and extractScale::mean.

Referenced by MakeComparisonPlots(), and MakeComparisonPlotsWithRatio().

{
  double mean, bin1, bin2, fwhm;
  mean  = h->GetMean();
  bin1  = h->FindFirstBinAbove(h->GetMaximum()/2);
  bin2  = h->FindLastBinAbove(h->GetMaximum()/2);
  fwhm  = h->GetBinCenter(bin2 ) - h->GetBinCenter(bin1 );

  int ilo = 1;
  int ihi = h->GetNbinsX();

  int imx = h->GetMaximumBin();
  double mx = h->GetBinContent(imx);

  int ileft;
  for (ileft=imx; ileft>=ilo; --ileft){
    if (h->GetBinContent(ileft)<0.5*mx) {
      ++ileft;
      break;
    }
  }

  int iright;
  for (iright=imx; iright<=ihi; ++iright) {
    if (h->GetBinContent(iright)<0.5*mx) {
      break;
    }
  }

  TF1 *myFit = new TF1("myFit","gaus", h->GetBinLowEdge(ileft), h->GetBinLowEdge(iright));
  h->Fit(myFit,"R q 0");
  myFit->SetLineColor(h->GetLineColor());
  std::cout<<"Mean: " << myFit->GetParameter(1) << "+/-" << myFit->GetParError(1) <<std::endl;
  return;
}

void MakeComparisonPlots	(	std::vector< TH1D * >	TheHists,
		std::string	plotName,
		std::vector< std::string >	Labels,
		bool	var
	)

Definition at line 257 of file neutKEsyst_plot.C.

References om::cout, allTimeWatchdog::endl, FitGauss(), make_mec_shifts_plots::GetMaximum(), analysePickle::hist, cet::sqlite::max(), SaveAs(), and gen_hdf5record::size.

                                                                                                                 {
  TCanvas* hCan = new TCanvas(plotName.c_str(), plotName.c_str());
  hCan->SetLeftMargin(0.1);
  int NBins = TheHists[0]->GetNbinsX();
  double max = 0.0;

  // --- Legend ---
  double HistLeg_X1 = 0.7, HistLeg_Y1 = 0.65, HistLeg_X2 = 0.9, HistLeg_Y2 = 0.88;
  TLegend* Leg = new TLegend(HistLeg_X1, HistLeg_Y1, HistLeg_X2, HistLeg_Y2 );
  Leg->SetFillStyle(0);
  Leg->SetBorderSize(0);
  double size = 0.9*TheHists[0]->GetXaxis()->GetTitleSize();
  Leg->SetTextSize( size );

  for (size_t hist=0; hist<TheHists.size(); ++hist) {
    TheHists[hist]->SetLineColor( 1+hist );
    if (hist==0) {
      TheHists[hist]->GetYaxis()->SetTitleOffset(0.7);
      TheHists[hist]->DrawCopy("hist");
      if(TString(plotName).Contains("NuE")){
        std::cout<<"PlotName: "<<plotName[hist]<<" Label: "<<Labels[hist]<<std::endl;
        FitGauss(TheHists[hist]);
      }
    }
    else{
      TheHists[hist] -> DrawCopy("hist same");
      if(TString(plotName).Contains("NuE")){
        std::cout<<"PlotName: "<<plotName[hist]<<" Label: "<<Labels[hist]<<std::endl;
        FitGauss(TheHists[hist]);
      }
    }
    Leg->AddEntry(TheHists[hist], Labels[hist].c_str());
    if(TheHists[hist]->GetMaximum() > max) max = TheHists[hist]->GetMaximum();
  }
  // set the new maximum
  TheHists[0]->SetMaximum(max);
  Leg->Draw();

  hCan -> SaveAs(TString("Plots/")+TString(hCan->GetName())+TString(".pdf"));
  hCan -> SaveAs(TString("Plots/")+TString(hCan->GetName())+TString(".png"));
}

void MakeComparisonPlotsWithRatio	(	std::vector< TH1D * >	TheHists,
		std::string	plotName,
		std::vector< std::string >	Labels,
		bool	var
	)

Definition at line 299 of file neutKEsyst_plot.C.

References cd(), om::cout, Draw(), allTimeWatchdog::endl, FitGauss(), make_mec_shifts_plots::GetMaximum(), analysePickle::hist, make_syst_table_plots::line, MakeRatio(), cet::sqlite::max(), allTimeWatchdog::plotName, SaveAs(), gen_hdf5record::size, string, submit_syst::x, and submit_syst::y.

Referenced by neutKEsyst_plot().

                                                                                                                          {
  // making two canvases -- one with just the comparison, and one with a ratio 

  TCanvas* hCan = new TCanvas(plotName.c_str(), plotName.c_str());
  hCan->SetLeftMargin(0.1);
  int NBins = TheHists[0]->GetNbinsX();
  double max = 0.0;

  // --- Legend ---
  double HistLeg_X1 = 0.7, HistLeg_Y1 = 0.65, HistLeg_X2 = 0.9, HistLeg_Y2 = 0.88;
  TLegend* Leg = new TLegend(HistLeg_X1, HistLeg_Y1, HistLeg_X2, HistLeg_Y2 );
  Leg->SetFillStyle(0);
  Leg->SetBorderSize(0);
  double size = 0.9*TheHists[0]->GetXaxis()->GetTitleSize();
  Leg->SetTextSize( size );

  for (size_t hist=0; hist<TheHists.size(); ++hist) {
    TheHists[hist]->SetLineColor( 1+hist );
    if (hist==0) {
      TheHists[hist]->GetYaxis()->SetTitleOffset(0.7);
      TheHists[hist]->DrawCopy("hist");
      if(TString(plotName).Contains("NuE")){
        std::cout<<"PlotName: "<<plotName<<" Label: "<<Labels[hist]<<std::endl;
        FitGauss(TheHists[hist]);
      }
    }
    else{
      TheHists[hist] -> DrawCopy("hist same");
      if(TString(plotName).Contains("NuE")){
        std::cout<<"PlotName: "<<plotName<<" Label: "<<Labels[hist]<<std::endl;
        FitGauss(TheHists[hist]);
      }
    }
    Leg->AddEntry(TheHists[hist], Labels[hist].c_str());
    if(TheHists[hist]->GetMaximum() > max) max = TheHists[hist]->GetMaximum();
  }
  // set the new maximum
  TheHists[0]->SetMaximum(max);
  Leg->Draw();

  // now we're going to make the canvas with the plot + ratio plot on it
  std::string rCname = "r"+plotName;

  // --- Get ready to make the canvas with ratios.
  TCanvas* rC = new TCanvas(rCname.c_str(), rCname.c_str());

  // Sort out the canvas.
  rC -> SetBottomMargin(0.);
  double Spl = 0.3;
  TPad* P1 = new TPad( "Temp_1", "", 0.0, Spl, 1.0, 1.0, 0 );
  TPad* P2 = new TPad( "Temp_2", "", 0.0, 0.0, 1.0, Spl, 0 );
  P2 -> SetRightMargin (.03);
  P2 -> SetTopMargin   (.00);
  P2 -> SetBottomMargin(.22);
  P2 -> SetLeftMargin  (.1);
  P2 -> Draw();
  P1 -> SetRightMargin (.03);
  P1 -> SetLeftMargin  (.1);
  P1 -> SetTopMargin   (.10);
  P1 -> SetBottomMargin(.00);
  P1 -> Draw();
  // Set some label sizes.
  double Lb1 = 0.07;
  double Lb2 = 0.10;
  // --- First, draw the ratios so cd onto Pad2
  P2 -> cd();

  TH1D* Rat = MakeRatio(TheHists[0], TheHists[1], kBlack, rCname, var); // Make the ratio.
  // Set axis ranges etc.
  Rat->GetYaxis()->SetTitleSize( Lb2 );
  Rat->GetYaxis()->SetTitleOffset(0.4);
  Rat->GetYaxis()->SetLabelSize( Lb2 );
  Rat->GetXaxis()->SetTitleSize( Lb2 );
  Rat->GetXaxis()->SetLabelSize( Lb2 );
  Rat->Draw();
  // Add a reference line
  double x[2] = {TheHists[0]->GetXaxis()->GetXmin(), TheHists[0]->GetXaxis()->GetXmax()};
  double y[2] = {1, 1};
  TGraph *line = new TGraph(2, x, y);
  line->SetLineColor(kGray);
  line->SetLineWidth(3);
  line->SetLineStyle(2);
  line->Draw("l same");

  Rat->Draw("axis same");
  P1->cd();

  // hist goes on P1  
  double Leg_X1 = 0.75, Leg_Y1 = 0.65, Leg_X2 = 0.95, Leg_Y2 = 0.88;
  TLegend* Leg2 = new TLegend(Leg_X1, Leg_Y1, Leg_X2, Leg_Y2 );
  Leg2->SetFillStyle (0);
  Leg2->SetBorderSize(0);
  Leg2->SetTextSize( 1.2*size );

  max = 0.0;
  for (size_t hist=0; hist<TheHists.size(); ++hist) {
    TheHists[hist]->SetLineColor( 1+hist );
    if (hist==0) {
      TheHists[hist]->GetYaxis()->SetTitleSize( Lb1 );
      TheHists[hist]->GetYaxis()->SetLabelSize( Lb1 );
      TheHists[hist]->GetYaxis()->SetTitleOffset( 0.7 );
      // Remove the x axis labels
      TheHists[hist]->GetXaxis()->SetLabelSize  (0 );
      TheHists[hist]->GetXaxis()->SetLabelOffset(99);
      TheHists[hist]->DrawCopy("hist");
    }
    else{
      TheHists[hist]->DrawCopy("hist same");
    }
    Leg2->AddEntry(TheHists[hist], Labels[hist].c_str());
    if(TheHists[hist]->GetMaximum() > max) max = TheHists[hist]->GetMaximum();
  }
  // set the new maximum
  TheHists[0]->SetMaximum(max);
  Leg2->Draw();
 
  hCan -> SaveAs( TString("Plots/")+TString(hCan->GetName())+TString(".pdf") );
  hCan -> SaveAs( TString("Plots/")+TString(hCan->GetName())+TString(".png") );
  rC -> SaveAs( TString("Plots/")+TString(rC->GetName())+TString(".pdf") );
  rC -> SaveAs( TString("Plots/")+TString(rC->GetName())+TString(".png") );

  return;
}

TH1D * MakeRatio	(	TH1D *	num,
		TH1D *	denom,
		int	Col,
		std::string	FType,
		bool	var
	)

Definition at line 423 of file neutKEsyst_plot.C.

References ana::bins, and art::to_string().

Referenced by MakeComparisonPlotsWithRatio().

                                                                              {
  // I know this function is awful but it's 1 am so this is what we get 
  if(var){
    Double_t bins[] = { 0, 0.75, 1, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2, 2.25, 2.5, 2.75, 3, 3.5, 4, 5 };
    Int_t binnum = sizeof(bins)/sizeof(Double_t) - 1;
    TH1D* ratHist = new TH1D( TString(FType)+TString("_ratio")+TString(std::to_string(Col))+TString(num->GetName()),
                              TString(";")+TString(num->GetXaxis()->GetTitle())+TString(";Ratio"),
                              binnum,
                              bins
                            );
      // Don't want any exponents?
      ratHist->GetYaxis()->SetNoExponent();
      ratHist->GetXaxis()->SetNoExponent();
      // Divide the two input histograms...
      ratHist->Divide( num , denom );
      // Set Line and Marker colors / styles.
      ratHist->SetLineColor  ( Col );
      ratHist->SetMarkerColor( Col );
      ratHist->SetMarkerStyle( 20 );
      ratHist->SetMarkerSize ( 0  );
      // Figure out what the range should be...
      ratHist->GetYaxis()->SetRangeUser( 0, 2.0 );
      // Return my shiny new ratio plot!
      return ratHist;

  }

  else{
    TH1D* ratHist = new TH1D( TString(FType)+TString("_ratio")+TString(std::to_string(Col))+TString(num->GetName()),
                              TString(";")+TString(num->GetXaxis()->GetTitle())+TString(";Ratio"),
                              num->GetNbinsX(),
                              num->GetXaxis()->GetXmin(),
                              num->GetXaxis()->GetXmax()
                            );
      // Don't want any exponents?
      ratHist->GetYaxis()->SetNoExponent();
      ratHist->GetXaxis()->SetNoExponent();
      // Divide the two input histograms...
      ratHist->Divide( num , denom );
      // Set Line and Marker colors / styles.
      ratHist->SetLineColor  ( Col );
      ratHist->SetMarkerColor( Col );
      ratHist->SetMarkerStyle( 20 );
      ratHist->SetMarkerSize ( 0  );
      // Figure out what the range should be...
      ratHist->GetYaxis()->SetRangeUser( 0, 2.0 );
      // Return my shiny new ratio plot!
      return ratHist;

    }
}

void neutKEsyst_plot

(

)

Definition at line 83 of file neutKEsyst_plot.C.

References fillBadChanDBTables::det, FCut_FD, FCut_ND, FD_CutNames, GENIEStr, runNovaSAM::gErrorIgnoreLevel, MySpectra::HadE_nom_FD, MySpectra::HadE_nom_ND, MySpectra::HadE_scale_m_FD, MySpectra::HadE_scale_m_ND, MySpectra::HadE_scale_p_FD, MySpectra::HadE_scale_p_ND, HCNames, MakeComparisonPlotsWithRatio(), ND_CutNames, nFD_CutNames, nGENIE, nHC, nND_CutNames, MySpectra::NNeu_FD, MySpectra::NNeu_ND, NormByBinWidth(), nQuantNames, MySpectra::NuE_nom_FD, MySpectra::NuE_nom_ND, MySpectra::NuE_scale_m_FD, MySpectra::NuE_scale_m_ND, MySpectra::NuE_scale_p_FD, MySpectra::NuE_scale_p_ND, POT, Print1D(), QuantNames, SpecToHist(), string, MySpectra::VisE_FD, MySpectra::VisE_ND, MySpectra::VisE_sum_FD, and MySpectra::VisE_sum_ND.

                       {
  // --- I don't want to see the print messages --- 
  gErrorIgnoreLevel = kError;

  // --- NOvA style --- 
  gStyle->SetOptStat(0);
  gROOT->SetStyle("novaStyle");
  TGaxis::SetMaxDigits(4);

  MySpectra sp;

  std::string BaseCAFFileDir="/nova/ana/users/essmith/neutronKE/syst/";
  std::string FDName, NDName;

  // ----- Spectra -----
  // pulling spectra out of files and converting to histograms
  // also want to normalize some of these plots by bin width (but not others)
  // best to do it here instead of in comparison plot code or I'll normalize
  // the first one a bunch of times by accident 

  // fhc/rhc
  for (unsigned int HC = 0; HC<nHC; ++HC){

    FDName = BaseCAFFileDir+"nKEsystBirks_FD"+HCNames[HC]+".root"; 
    TFile *FDFile = TFile::Open(FDName.c_str());

    NDName = BaseCAFFileDir+"nKEsystBirks_ND"+HCNames[HC]+".root";
    TFile *NDFile = TFile::Open(NDName.c_str());

    // interaction types
    for (unsigned int gen=0; gen<nGENIE; ++gen) {

      // quantiles
      for (unsigned int quant=0; quant<nQuantNames; ++quant) {

        // analysis cuts, ND
        for(unsigned int det=0; det<nND_CutNames; ++det) {

          std::string Comb = GENIEStr[gen]+"_"+ND_CutNames[det]+"_"+QuantNames[quant];

          sp.NuE_nom_ND         [gen][det][quant][HC] = SpecToHist(NDFile, "NuE_nom_"+Comb,
                                                        POT[HC], ";Reco #nu E (GeV);Events/0.1GeV");
          sp.HadE_nom_ND        [gen][det][quant][HC] = SpecToHist(NDFile, "HadE_nom_"+Comb, 
                                                        POT[HC], ";Reco Had E (GeV);Events");
          sp.NuE_scale_p_ND     [gen][det][quant][HC] = SpecToHist(NDFile, "NuE_scale_p_"+Comb, 
                                                        POT[HC], ";Reco #nu E (GeV);Events/0.1GeV");
          sp.HadE_scale_p_ND    [gen][det][quant][HC] = SpecToHist(NDFile, "HadE_scale_p_"+Comb, 
                                                        POT[HC], ";Reco Had E (GeV);Events");
          sp.NuE_scale_m_ND     [gen][det][quant][HC] = SpecToHist(NDFile, "NuE_scale_m_"+Comb, 
                                                        POT[HC], ";Reco #nu E (GeV);Events/0.1GeV");
          sp.HadE_scale_m_ND    [gen][det][quant][HC] = SpecToHist(NDFile, "HadE_scale_m_"+Comb, 
                                                        POT[HC], ";Reco Had E (GeV);Events");
          sp.NNeu_ND            [gen][det][quant][HC] = SpecToHist(NDFile, "NNeu_"+Comb, 
                                                        POT[HC], ";Multiplicity");
          sp.VisE_sum_ND        [gen][det][quant][HC] = SpecToHist(NDFile, "VisE_sum_"+Comb,
                                                        POT[HC], ";Summed Vis E (GeV)");
          sp.VisE_ND            [gen][det][quant][HC] = SpecToHist(NDFile, "VisE_"+Comb, 
                                                        POT[HC], ";Vis E (GeV)");

          // normalize the plots with variable binning by bin width
          NormByBinWidth(sp.NuE_nom_ND[gen][det][quant][HC]);
          NormByBinWidth(sp.NuE_scale_p_ND[gen][det][quant][HC]);
          NormByBinWidth(sp.NuE_scale_m_ND[gen][det][quant][HC]);
       }

        // analysis cuts, FD 
        for(unsigned int det=0; det<nFD_CutNames; ++det) {

          std::string Comb = GENIEStr[gen]+"_"+FD_CutNames[det]+"_"+QuantNames[quant];

          sp.NuE_nom_FD        [gen][det][quant][HC] = SpecToHist(FDFile, "NuE_nom_"+Comb,
                                                       POT[HC], ";Reco #nu E (GeV);Events/0.1GeV");
          sp.HadE_nom_FD       [gen][det][quant][HC] = SpecToHist(FDFile, "HadE_nom_"+Comb, 
                                                       POT[HC], ";Reco Had E (GeV);Events");
          sp.NuE_scale_p_FD    [gen][det][quant][HC] = SpecToHist(FDFile, "NuE_scale_p_"+Comb, 
                                                       POT[HC], ";Reco #nu E (GeV);Events/0.1GeV");
          sp.HadE_scale_p_FD   [gen][det][quant][HC] = SpecToHist(FDFile, "HadE_scale_p_"+Comb, 
                                                       POT[HC], ";Reco Had E (GeV);Events");
          sp.NuE_scale_m_FD    [gen][det][quant][HC] = SpecToHist(FDFile, "NuE_scale_m_"+Comb, 
                                                       POT[HC], ";Reco #nu E (GeV);Events/0.1GeV");
          sp.HadE_scale_m_FD   [gen][det][quant][HC] = SpecToHist(FDFile, "HadE_scale_m_"+Comb, 
                                                       POT[HC], ";Reco Had E (GeV);Events");
          sp.NNeu_FD           [gen][det][quant][HC] = SpecToHist(FDFile, "NNeu_"+Comb, 
                                                       POT[HC], ";Multiplicity");
          sp.VisE_sum_FD       [gen][det][quant][HC] = SpecToHist(FDFile, "VisE_sum_"+Comb,
                                                       POT[HC], ";Summed Vis E (GeV)");
          sp.VisE_FD           [gen][det][quant][HC] = SpecToHist(FDFile, "VisE_"+Comb, 
                                                       POT[HC], ";Vis E (GeV)");

          // normalize the plots with variable binning by bin width
          NormByBinWidth(sp.NuE_nom_FD[gen][det][quant][HC]);
          NormByBinWidth(sp.NuE_scale_p_FD[gen][det][quant][HC]);
          NormByBinWidth(sp.NuE_scale_m_FD[gen][det][quant][HC]);
        }
      }
    }
  }

  // ----- comparison plots -----
  // 

  //final cut plots
  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];
        MakeComparisonPlotsWithRatio({sp.NuE_nom_ND[gen][FCut_ND][quant][HC], sp.NuE_scale_p_ND[gen][FCut_ND][quant][HC]}, 
                            "NuE_scale_p_ND_"+Comb, {"no scale", "scale up"}, true);
        MakeComparisonPlotsWithRatio({sp.NuE_nom_ND[gen][FCut_ND][quant][HC], sp.NuE_scale_m_ND[gen][FCut_ND][quant][HC]}, 
                            "NuE_scale_m_ND_"+Comb, {"no scale", "scale down"}, true);
        MakeComparisonPlotsWithRatio({sp.HadE_nom_ND[gen][FCut_ND][quant][HC], sp.HadE_scale_p_ND[gen][FCut_ND][quant][HC]}, 
                            "HadE_scale_p_ND_"+Comb, {"no scale", "scale up"}, false);
        MakeComparisonPlotsWithRatio({sp.HadE_nom_ND[gen][FCut_ND][quant][HC], sp.HadE_scale_m_ND[gen][FCut_ND][quant][HC]}, 
                            "HadE_scale_m_ND_"+Comb, {"no scale", "scale down"}, false);
        Print1D(sp.NNeu_ND    [gen][FCut_ND][quant][HC], "NNeu_ND_"+Comb);
        Print1D(sp.VisE_sum_ND[gen][FCut_ND][quant][HC], "VisE_sum_ND_"+Comb);
        Print1D(sp.VisE_ND    [gen][FCut_ND][quant][HC], "VisE_ND_"+Comb);

        MakeComparisonPlotsWithRatio({sp.NuE_nom_FD[gen][FCut_FD][quant][HC], sp.NuE_scale_p_FD[gen][FCut_FD][quant][HC]}, 
                            "NuE_scale_p_FD_"+Comb, {"no scale", "scale up"}, true);
        MakeComparisonPlotsWithRatio({sp.NuE_nom_FD[gen][FCut_FD][quant][HC], sp.NuE_scale_m_FD[gen][FCut_FD][quant][HC]}, 
                            "NuE_scale_m_FD_"+Comb, {"no scale", "scale down"}, true);
        MakeComparisonPlotsWithRatio({sp.HadE_nom_FD[gen][FCut_FD][quant][HC], sp.HadE_scale_p_FD[gen][FCut_FD][quant][HC]}, 
                            "HadE_scale_p_FD_"+Comb, {"no scale", "scale up"}, false);
        MakeComparisonPlotsWithRatio({sp.HadE_nom_FD[gen][FCut_FD][quant][HC], sp.HadE_scale_m_FD[gen][FCut_FD][quant][HC]}, 
                            "HadE_scale_m_FD_"+Comb, {"no scale", "scale down"}, false);

        Print1D(sp.NNeu_FD    [gen][FCut_FD][quant][HC], "NNeu_FD_"+Comb);
        Print1D(sp.VisE_sum_FD[gen][FCut_FD][quant][HC], "VisE_sum_FD_"+Comb);
        Print1D(sp.VisE_FD    [gen][FCut_FD][quant][HC], "VisE_FD_"+Comb);

      } // quants
    } // int. types 
  } // hc

  return;
}

void NormByBinWidth

(

TH1D *

hist

)

Definition at line 242 of file neutKEsyst_plot.C.

Referenced by neutKEsyst_plot().

                               {
  hist->Scale(0.1, "width");
  hist->SetMaximum(hist->GetMaximum()*1.1);
}

void Print1D	(	TH1D *	hist,
		std::string	plotName
	)

Definition at line 248 of file neutKEsyst_plot.C.

References SaveAs().

Referenced by neutKEsyst_plot().

                                            {
  TCanvas* hCan = new TCanvas(plotName.c_str(), plotName.c_str());
  hCan->SetLeftMargin(0.1);
  hist->DrawCopy("hist");
  hCan -> SaveAs(TString("Plots/")+TString(hCan->GetName())+TString(".pdf"));
  hCan -> SaveAs(TString("Plots/")+TString(hCan->GetName())+TString(".png"));

}

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

Definition at line 224 of file neutKEsyst_plot.C.

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

Referenced by neutKEsyst_plot().

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

  //area norm
//  MyHist->Scale(1./MyHist->Integral());
 
  // comment this out to keep the 0bin
//  MyHist->GetXaxis()->SetRange(2,MyHist->GetNbinsX());
 
  return MyHist;
}

Variable Documentation

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

Definition at line 48 of file neutKEsyst_plot.C.

unsigned int FCut_FD = nFD_CutNames - 1

Definition at line 38 of file neutKEsyst_plot.C.

Referenced by neutKEsyst_plot().

unsigned int FCut_ND = nND_CutNames - 1

Definition at line 36 of file neutKEsyst_plot.C.

Referenced by neutKEsyst_plot().

std::vector<std::string> FD_CutNames = {"FD_2018full"}

Definition at line 37 of file neutKEsyst_plot.C.

Referenced by neutKEsyst_plot().

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

Definition at line 31 of file neutKEsyst_plot.C.

Referenced by neutKEsyst_plot().

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

Definition at line 44 of file neutKEsyst_plot.C.

Referenced by neutKEsyst_plot().

std::vector<std::string> ND_CutNames = {"ND_2018full"}

Definition at line 35 of file neutKEsyst_plot.C.

Referenced by neutKEsyst_plot().

const unsigned int nDet = 2

Definition at line 48 of file neutKEsyst_plot.C.

Referenced by hyperon_macro(), and hyperon_nom_macro().

const unsigned int nFD_CutNames = 1

Definition at line 37 of file neutKEsyst_plot.C.

Referenced by neutKEsyst_plot().

const unsigned int nGENIE = 1

Definition at line 31 of file neutKEsyst_plot.C.

Referenced by neutKEsyst_plot().

const unsigned int nHC = 2

Definition at line 44 of file neutKEsyst_plot.C.

Referenced by neutKEsyst_plot().

const unsigned int nND_CutNames = 1

Definition at line 35 of file neutKEsyst_plot.C.

Referenced by neutKEsyst_plot().

const unsigned int nQuantNames = 5

Definition at line 41 of file neutKEsyst_plot.C.

Referenced by neutKEsyst_plot().

std::vector<double> POT = {kAna2018FHCPOT, kAna2018RHCPOT}

Definition at line 45 of file neutKEsyst_plot.C.

Referenced by ana::AddErrorInQuadrature(), AddExptCorrNue(), AddExptCorrNumu(), Analyse_GetEfficiency(), Analyse_GetEfficiency_UseNEntries(), calculateWrongSignNue(), calculateWrongSignNumuQ1(), calculateWrongSignNumuQ2(), calculateWrongSignNumuQ3(), calculateWrongSignNumuQ4(), compare_cos_numi(), datamc(), demoFitContours(), demoFitSlices(), drawLongTerm(), ana::CovMxSurface::DrawSurfacePoint(), drawSystsShiftingNDdata(), drawTimePlots(), drawVsPOT(), EnergyCont_plot(), fake_future_data(), fd_plot(), ana::GenieMultiverseSpectra::FindBandBoundaries(), ana::FluxMultiverseSyst::FindSigmaBoundaries(), get_cosmic_sample(), ana::CutOptimization::GetAllSystDefs(), GetHistsFD(), GetHistsND(), Plotter::GetMCPlots(), ana::Multiverse::GetNSigmaShift(), GetNueData2018(), GetNueData2019(), ana::GetNueData2020(), GetNumuData2019(), ana::GetNumuData2020(), ana::Multiverse::GetNUniverses(), getPOT(), getTimePeakPlots(), getTotPOT(), goodness_of_fit(), hyperon_plot(), joint_fit_2018_contours(), joint_fit_2018_slices(), joint_fit_2019_bestfit(), joint_fit_2019_contours(), joint_fit_2019_slices(), joint_fit_future_bestfit_univ(), joint_fit_future_contour_univ(), ana::TrivialCrossSectionAnalysis::LoadFrom(), make1DsAndRatio(), make_fc_surfaces_2020(), make_mass_and_oct_fc_input_2019(), make_quantiles_histogram_2020(), MakeCutFlow(), makeEnergyEstimator(), makeFlatWeight(), MakeLegend(), MakeNus17CosBkgd(), makeResolutionPlots(), mre_comp_split(), multiverse_reweighting(), neutKEsyst_plot(), neutronE_plot(), ana::GenieMultiverseNormalizedSpectra::NormalizeSpectra(), nuebar_signif(), pi0_xcheck(), plot_datapredictions(), plot_pi0_xcheck(), plot_pid(), plot_predictions(), PandAna.reco_validation.prod5_pid_validation::plot_ratio(), plot_rationoosc(), ana::PlotAllRecoTrue(), ana::PlotAllSignalEstimates(), PlotStack(), Plotting_DataAndPrediction(), Plotting_OverlayStudies(), ana::SpectrumComponents::Purity(), run_joint_fit_2020_bestfit(), run_joint_fit_2020_contours(), run_joint_fit_2020_slices(), selection_story_plots(), ana::ICrossSectionAnalysis::SetPOT(), ana::SpectrumComponents::SpectrumComponents(), syst_table_fit(), ana::Multiverse::ToSpectrum(), ana::FluxMultiverseSyst::ToTH1(), ana::GenieMultiverseSpectra::ToTH1(), ana::GenieMultiverseSpectra::ToTH2(), ana::Multiverse::Transform(), Tutorial2019Fit(), Tutorial2019FitContours(), Tutorial2019FitSlices(), and ana::NDPredictionNoOsc::~NDPredictionNoOsc().

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

Definition at line 41 of file neutKEsyst_plot.C.

Referenced by neutKEsyst_plot().