PlotSpectra.h File Reference

#include "StandardRecord/Proxy/SRProxy.h"
#include "CAFAna/Core/Binning.h"
#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Core/SpectrumLoader.h"
#include "CAFAna/Core/EventList.h"
#include "CAFAna/Core/Utilities.h"
#include "CAFAna/Cuts/Cuts.h"
#include "CAFAna/Vars/Vars.h"
#include "3FlavorAna/Vars/NueVarsExtra.h"
#include "3FlavorAna/Vars/NumuVarsExtra.h"
#include "3FlavorAna/Cuts/NueCuts2020.h"
#include "3FlavorAna/Cuts/NumuCuts2020.h"
#include "NuXAna/Vars/NusVars.h"
#include "NuXAna/Cuts/NusCuts20.h"
#include "NuXAna/Cuts/NusCuts.h"
#include "NuXAna/macros/Nus20/Selection/CutFlow.h"
#include "NuXAna/macros/Nus20/Selection/CutVariables.h"
#include "NuXAna/macros/Nus20/Selection/FOMUtilities.h"
#include "NuXAna/macros/Nus20/Selection/DecorrelateMatrix.h"
#include "TFile.h"
#include "TString.h"
#include "TH1.h"
#include "TCanvas.h"
#include "TPad.h"
#include "THStack.h"
#include "TLegend.h"
#include "TBox.h"
#include "TLatex.h"
#include "TStyle.h"
#include "/nashome/a/alister1/NOvAPlotting.h"

Go to the source code of this file.

Enumerations
enum	Detector { kNear, kFar }
	helper plotting scripts More...

Functions
void	LoopHistograms (TFile *spectraFile, Detector thisDet)
void	StyleHistograms (TH1D *h, Detector det, bool isBeamData, bool isCosmic, bool isTotal)
void	SetupCanvas (TPad *topPad, TPad *botPad, TPad *thisPad)
bool	PlotData (std::string name)
void	StyleLegend (TLegend *leg)
void	DrawBoxes (TH1D *histMCTotal, std::vector< double > cutValues, CutSide thisCutSide)
void	DrawRatio (TPad *thisPad, TH1D *mc, TH1D *dat)
void	GetFractionalUncertainties (TFile *covmx, std::vector< float > *nd, std::vector< float > *fd)
void	InitialiseHistograms (TFile *spectraFile, Detector thisDet, int iCut, int iVar, std::string dirName, std::string histPathPre, std::string histPath)
void	InitialiseFOMHists (Detector thisDet, int iCut, int iVar)
void	DrawFOM (std::set< FomType >, int iVar)
void	PrintSelectedEvents (int iCut, int iVar)
void	DrawDataMC (Detector thisDet, int iCut, int iVar, bool isLogPlot, bool isCosmic, bool isAreaNormalised)
void	DrawMCOnly (Detector thisDet, int iCut, int iVar, bool isLogPlot, bool isCosmic)
void	DrawIndependentSpectra (Detector thisDet, int iCut, int iVar)
TLegend *	GetStandardLegend ()

Variables
std::vector< float >	ndFractionalUncertainties = {}
std::vector< float >	fdFractionalUncertainties = {}
std::string	dirName
std::vector< VarContainer >	vars
std::vector< CutContainer >	cuts
std::vector< CutContainer >	subcuts
TH1D *	histTotal = new TH1D()
TH1D *	histMCTotal = new TH1D()
	total sig + bg More...
TH1D *	histMCSignal = new TH1D()
	total mc only More...
TH1D *	histMCSignalNoCut = new TH1D()
	total mc signal More...
TH1D *	histCosmic = new TH1D()
	total mc signal (no cut) More...
TH1D *	histMCBg = new TH1D()
	total cosmic More...
TH1D *	histData = new TH1D()
	mc-only bg (no cosmic) More...
TH1D *	EffTimesPurHist = new TH1D()
	total data More...
TH1D *	SOverSBHist = new TH1D()
TH1D *	SOverSSBHist = new TH1D()
TH1D *	BinnedSOverSSBSystHist = new TH1D()
TH2D *	histMCTotalEnergy = nullptr
TH2D *	histMCSignalEnergy = nullptr
bool	isAreaNorm = true
TCanvas *	cRatio = new TCanvas("cRatio", "cRatio", 600, 600)
TPad *	topPad = new TPad("topPad", "", 0.005, 0.3, 0.995, 0.995)
TPad *	botPad = new TPad("botPad", "", 0.005, 0.005, 0.995, 0.3)
TCanvas *	cNoRatio = new TCanvas("cNoRatio", "cNoRatio", 600, 500)
TPad *	thisPad = new TPad("thisPad", "", 0.005, 0.005, 0.995, 0.995)

Enumeration Type Documentation

enum Detector

helper plotting scripts

Enumerator
kNear
kFar

Definition at line 39 of file PlotSpectra.h.

              {
  kNear,
  kFar
};

Function Documentation

void DrawBoxes	(	TH1D *	histMCTotal,
		std::vector< double >	cutValues,
		CutSide	thisCutSide
	)

Definition at line 233 of file PlotSpectra.C.

References ana::kBoth, kHigh, and kLow.

Referenced by DrawDataMC(), DrawIndependentSpectra(), and DrawMCOnly().

                                                                                   {

  TBox* box  = new TBox(-5000, -5000, -5000, -5000);
  TBox* box2 = new TBox(-5000, -5000, -5000, -5000); 

  box->SetY1(0);
  box->SetY2(histMCTotal->GetMaximum());
  if (thisCutSide == kLow){
    box->SetX1(histMCTotal->GetBinLowEdge(1));
    box->SetX2(theseCuts.at(0));
  }
  if (thisCutSide == kHigh){
    box->SetX1(theseCuts.at(0));
    box->SetX2(histMCTotal->GetBinLowEdge(histMCTotal->GetNbinsX()+1));
  }
  if (thisCutSide == kBoth){
    box->SetX1(histMCTotal->GetBinLowEdge(1));
    box->SetX2(theseCuts.at(0));
    box2->SetX1(theseCuts.at(1));
    box2->SetX2(histMCTotal->GetBinLowEdge(histMCTotal->GetNbinsX()+1));
    box2->SetLineWidth(0);
    box2->SetFillColorAlpha(kGray, 0.4);
    box2->SetY1(0);
    box2->SetY2(histMCTotal->GetMaximum());
    box2->Draw("same");
  }
  box->SetLineWidth(0);
  box->SetFillColorAlpha(kGray, 0.4);
  if (theseCuts.at(0) != 0 || theseCuts.size() > 1)
    box->Draw("same");
  
  gPad->RedrawAxis();
  gPad->Modified();

}

void DrawDataMC	(	Detector	thisDet,
		int	iCut,
		int	iVar,
		bool	isLogPlot,
		bool	isCosmic,
		bool	isAreaNormalised
	)

Definition at line 437 of file PlotSpectra.C.

References botPad, cRatio, update_sam_good_runs_metadata::cuts, dirName, DrawBoxes(), DrawRatio(), GetStandardLegend(), histCosmic, histData, histMCBg, histMCSignal, histMCTotal, histTotal, kNear, MECModelEnuComparisons::leg, scale, Scale(), string, topPad, and vars.

Referenced by LoopHistograms().

                                                                     {
 
  // need to define this inside function. THStacks are awkward
  THStack* histSigandBgStack = new THStack("histSigandBgStack", "histSigandBgStack");
  histSigandBgStack->Add(histMCSignal);
  if (histCosmic != nullptr)
    histSigandBgStack->Add(histCosmic);
  histSigandBgStack->Add(histMCBg);

  TCanvas* thisCanvas = cRatio;
  topPad->cd();

  std::vector<double> theseCuts = vars.at(iVar).cuts;
  CutSide thisCutSide           = vars.at(iVar).cutside;

  double maxVal;
  std::string plotType = "";
  if (isLogPlot){
    topPad->SetLogy();
    plotType += "Log";
    maxVal = 10000;
  }
  else{
    topPad->SetLogy(0);
    maxVal=1.4;
  }
  if (isAreaNormalised){
    plotType += "AreaNorm";
  }
  
  histMCTotal->GetYaxis()->SetRangeUser(0.0001, histMCTotal->GetMaximum()*maxVal);
  histTotal  ->GetYaxis()->SetRangeUser(0.0001, histTotal->GetMaximum()*maxVal);

  // if area normalised, we don't want to scale the data, but scale the mc to same integral
  int tB = histData->GetNbinsX()+1;
  double scale = isAreaNormalised ? histData->Integral(0,tB)/histMCTotal->Integral(0,tB) : 1./1e4;
  if (thisDet == kNear) {
    histMCTotal       -> Scale(scale);
    histMCSignal      -> Scale(scale);
    //histSigandBgStack -> Scale(scale); // no Scale function for THStack, but 
                                         // probably ok since don't want to 
                                         // area norm, and only exist in FD
    histTotal         -> Scale(scale);
    if (!isAreaNormalised)
      histData   ->Scale(scale);
  }
    
  if (isCosmic){
    histTotal  ->GetYaxis()->SetRangeUser(0.0001, histTotal->GetMaximum()*maxVal);
    histTotal->Draw("hist");
    histSigandBgStack->Draw("hist same");
    histTotal->Draw("hist same");
    DrawBoxes(histTotal, theseCuts, thisCutSide);
  }
  else{
    histMCTotal->GetYaxis()->SetRangeUser(0.0001, histMCTotal->GetMaximum()*maxVal);
    histMCTotal  ->Draw("hist");
    histMCSignal ->Draw("hist same");
    histMCTotal  ->Draw("hist same");
    DrawBoxes(histMCTotal, theseCuts, thisCutSide);
  }

  histData->Draw("p same");
  
  TLegend* leg = GetStandardLegend();
  leg->Draw("same");

  DrawRatio(botPad, histTotal, histData);

  thisCanvas->SaveAs((dirName+vars.at(iVar).name+cuts.at(iCut).name+plotType+"DataMC.pdf").c_str());
  thisCanvas->SaveAs((dirName+vars.at(iVar).name+cuts.at(iCut).name+plotType+"DataMC.png").c_str());
  
}

void DrawFOM	(	std::set< FomType >	,
		int	iVar
	)

Definition at line 352 of file PlotSpectra.C.

References EffTimesPurHist, kEffTimesPurity, kSOverSB, kSOverSSB, SOverSBHist, SOverSSBHist, StyleLegend(), and vars.

Referenced by LoopHistograms().

                                                {

  if (fomTypes.size() == 0) return;

  if ((vars.at(iVar).name.find("CVN") != std::string::npos) ||
                vars.at(iVar).name.find("HitsPerPlane") != std::string::npos ||
          vars.at(iVar).name.find("NSliceHits") != std::string::npos ||
      vars.at(iVar).name.find("PartPtp")  != std::string::npos){

    TLegend* fomleg = new TLegend(0.5, 0.7, 0.85, 0.89);
    if (fomTypes.count(kEffTimesPurity)){
          EffTimesPurHist->Draw("same hist");
          fomleg->AddEntry(EffTimesPurHist  , "Eff. #times Pur.");
    }
    if (fomTypes.count(kSOverSB)){
        SOverSBHist->Draw("same hist");
          fomleg->AddEntry(SOverSBHist      , "S / #sqrt{B}"    );
    }
    if (fomTypes.count(kSOverSSB)){
        SOverSSBHist->Draw("same hist");
          fomleg->AddEntry(SOverSSBHist     , "S / #sqrt{S+B+#sigma^{2}}");
    }
  
        StyleLegend(fomleg);
  
        fomleg->Draw("same");
  }
}

void DrawIndependentSpectra	(	Detector	thisDet,
		int	iCut,
		int	iVar
	)

Definition at line 583 of file PlotSpectra.C.

References cNoRatio, update_sam_good_runs_metadata::cuts, dirName, DrawBoxes(), histCosmic, histMCBg, histMCSignal, histMCTotal, histTotal, kRed, std::max(), StyleLegend(), thisPad, and vars.

Referenced by LoopHistograms().

{
  TCanvas* thisCanvas = cNoRatio;
  thisPad->cd();
  thisPad->SetLogy(0);

  std::vector<double> theseCuts = vars.at(iVar).cuts;
  CutSide thisCutSide           = vars.at(iVar).cutside;

  // annoyingly the standard offsets don't work
  histTotal->GetYaxis()->SetTitleOffset(1.1);
  histMCTotal->GetYaxis()->SetTitleOffset(1.1);
        histMCSignal->GetYaxis()->SetTitleOffset(1.1);

  histMCSignal->SetFillStyle(0);
  histMCBg    ->SetFillStyle(0);
  histCosmic  ->SetFillStyle(0);
  histMCSignal->Scale(1./histMCSignal->Integral());
  histMCBg    ->Scale(1./histMCBg->Integral());
  histMCBg    ->SetLineColor(kRed+1);
  histCosmic  ->Scale(1./histCosmic->Integral());
  histMCSignal->GetYaxis()->SetRangeUser(0.0001, std::max({histMCSignal->GetMaximum(), histCosmic->GetMaximum(), histMCBg->GetMaximum()})*1.4);
  histMCSignal->Draw("hist");
  histMCBg    ->Draw("hist same");
  histCosmic  ->Draw("hist same");
  
  DrawBoxes(histMCSignal, theseCuts, thisCutSide);
  
  TLegend *areaLeg = new TLegend(0.15, 0.70, 0.5, 0.85);
  StyleLegend(areaLeg);
  areaLeg->AddEntry(histMCSignal, "True NC Events", "l");
  areaLeg->AddEntry(histMCBg    , "MC Background" , "l");
  if (histCosmic != nullptr)
        areaLeg->AddEntry(histCosmic  , "Cosmic Events" , "l");
  areaLeg->Draw("same");

  thisCanvas->SaveAs((dirName+vars.at(iVar).name+cuts.at(iCut).name+"IndependentSpectra.pdf").c_str());
  thisCanvas->SaveAs((dirName+vars.at(iVar).name+cuts.at(iCut).name+"IndependentSpectra.png").c_str());
}

void DrawMCOnly	(	Detector	thisDet,
		int	iCut,
		int	iVar,
		bool	isLogPlot,
		bool	isCosmic
	)

Definition at line 513 of file PlotSpectra.C.

References cNoRatio, update_sam_good_runs_metadata::cuts, dirName, DrawBoxes(), GetStandardLegend(), histCosmic, histMCBg, histMCSignal, histMCTotal, histTotal, kNear, MECModelEnuComparisons::leg, scale, string, thisPad, and vars.

Referenced by LoopHistograms().

                                              {

  // need to define this inside function. THStacks are awkward
  THStack* histSigandBgStack = new THStack("histSigandBgStack", "histSigandBgStack");
  histSigandBgStack->Add(histMCSignal);
  if (histCosmic != nullptr)
    histSigandBgStack->Add(histCosmic);
  histSigandBgStack->Add(histMCBg);

  TCanvas* thisCanvas = cNoRatio;
  thisPad->cd();

  std::vector<double> theseCuts = vars.at(iVar).cuts;
  CutSide thisCutSide           = vars.at(iVar).cutside;

  // annoyingly the standard offsets don't work
  histTotal->GetYaxis()->SetTitleOffset(1.1);
  histMCTotal->GetYaxis()->SetTitleOffset(1.1);
        histMCSignal->GetYaxis()->SetTitleOffset(1.1);

  double maxVal;
  std::string plotType = "";
  if (isLogPlot){
    thisPad->SetLogy();
    plotType += "Log";
    maxVal = 10000;
  }
  else{
    thisPad->SetLogy(0);
    maxVal=1.4;
  }
  if (!isCosmic){
    plotType += "NoCosmics";
  }

  // if area normalised, we don't want to scale the data, but scale the mc to same integral
  double scale = 1./1e4;
  if (thisDet == kNear) {
    histMCTotal ->Scale(scale);
    histMCSignal->Scale(scale);
    histTotal   ->Scale(scale); 
  }
   
  if (isCosmic){
    histTotal  ->GetYaxis()->SetRangeUser(0.0001, histTotal->GetMaximum()*maxVal);
    histTotal->Draw("hist");
    histSigandBgStack->Draw("hist same");
    histTotal->Draw("hist same");
    DrawBoxes(histTotal, theseCuts, thisCutSide);
  }
  else{
    histMCTotal->GetYaxis()->SetRangeUser(0.0001, histMCTotal->GetMaximum()*maxVal);
    histMCTotal  ->Draw("hist");
    histMCSignal ->Draw("hist same");
    histMCTotal  ->Draw("hist same");
    DrawBoxes(histMCTotal, theseCuts, thisCutSide);
  }

  TLegend* leg = GetStandardLegend();
  if (isCosmic)
    leg->AddEntry(histCosmic, "Cosmic Events", "f");
  leg->Draw("same");

  thisCanvas->SaveAs((dirName+vars.at(iVar).name+cuts.at(iCut).name+plotType+"MCSpectra.pdf").c_str());
  thisCanvas->SaveAs((dirName+vars.at(iVar).name+cuts.at(iCut).name+plotType+"MCSpectra.png").c_str());

}

void DrawRatio	(	TPad *	thisPad,
		TH1D *	mc,
		TH1D *	dat
	)

Definition at line 270 of file PlotSpectra.C.

Referenced by DrawDataMC().

                                                   {

  thisPad->cd();
  thisPad->SetGridy();
  TH1D* rat = (TH1D*)data->Clone("rat");
  rat->Divide(mc);
  rat->GetYaxis()->SetRangeUser(0.7, 1.29);
  rat->GetYaxis()->SetNdivisions(505);
  rat->GetYaxis()->SetTitle("Data/Sim");
  rat->GetYaxis()->SetTitleOffset(1.3);
  rat->GetXaxis()->SetTitleOffset(3.0);

  rat->Draw("p");

}

void GetFractionalUncertainties	(	TFile *	covmx,
		std::vector< float > *	nd,
		std::vector< float > *	fd
	)

Definition at line 287 of file PlotSpectra.C.

References MECModelEnuComparisons::i, and std::sqrt().

                                                                                                                                                     {

        int ndBins = 20;
        int fdBins = 13;

        TMatrixD* cov = (TMatrixD*)covmx->Get("fraccovmx");
        //TMatrixD* cov = (TMatrixD*)covmx->Get("decorrerr");
        
        for (int i =0; i < cov->GetNrows(); ++i){
                if (i < ndBins){
                        nd->push_back(std::sqrt((*cov)(i,i)));
                }
                else {
                        fd->push_back(std::sqrt((*cov)(i,i)));
                }
        }
}

TLegend* GetStandardLegend

(

)

Definition at line 623 of file PlotSpectra.C.

References histData, histMCSignal, histTotal, MECModelEnuComparisons::leg, and StyleLegend().

Referenced by DrawDataMC(), and DrawMCOnly().

                            {
  TLegend* leg = new TLegend(0.15, 0.7, 0.5, 0.89);
  leg->AddEntry(histTotal , "All Events", "l");
  leg->AddEntry(histMCSignal, "True NC Events", "f");
  if (histData != nullptr)
    leg->AddEntry(histData, "Data", "p");
  StyleLegend(leg);
  return leg;
}

void InitialiseFOMHists	(	Detector	thisDet,
		int	iCut,
		int	iVar
	)

Definition at line 308 of file PlotSpectra.C.

References BinnedSOverSSBSystHist, update_sam_good_runs_metadata::cuts, dirName, EffTimesPurHist, fdFractionalUncertainties, getFOMHist(), histMCSignal, histMCSignalEnergy, histMCSignalNoCut, histMCTotal, histMCTotalEnergy, histTotal, kBinnedSOverSSB, kEffTimesPurity, kNear, kSOverSB, kSOverSSB, ndFractionalUncertainties, SOverSBHist, SOverSSBHist, StyleFOMs(), and vars.

Referenced by LoopHistograms().

                                                             {
  // setup FOM histograms
        EffTimesPurHist  = (TH1D*)getFOMHist("EffTimesPur"+dirName+vars.at(iVar).name+cuts.at(iCut).name,
                                        vars.at(iVar).cutside,
                                                                                                                                                   FomType::kEffTimesPurity,
                                                                                                                                                         histMCTotal,
                                                                                                                                                         histMCSignal,
                                                                                                                                                         histMCSignalNoCut);

        SOverSBHist      = (TH1D*)getFOMHist("SOverSB"+dirName+vars.at(iVar).name+cuts.at(iCut).name,
                                                                                                                                                         vars.at(iVar).cutside,
                                                                                                                                                         FomType::kSOverSB,
                                                                                                                                                         histMCTotal,
                                                                                                                                                         histMCSignal,
                                                                                                                                                         histMCSignalNoCut);


        BinnedSOverSSBSystHist = (TH1D*)getFOMHist("BinnedSOverSSBSyst"+dirName+vars.at(iVar).name+cuts.at(iCut).name,
                                                                                                                                                                           vars.at(iVar).cutside,
                                                                                                                                                                           FomType::kBinnedSOverSSB,
                                                                                                                                                                           histMCTotal,
                                                                                                                                                                           histMCSignal,
                                                                                                                                                                           histMCSignalNoCut,
                                                                                                                                                                           histMCTotalEnergy,
                                                                                                                                                                           histMCSignalEnergy,
                                                                                                                                                                           thisDet == kNear ? ndFractionalUncertainties : fdFractionalUncertainties);

        SOverSSBHist          = (TH1D*)getFOMHist("SOverSSB"+dirName+vars.at(iVar).name+cuts.at(iCut).name,
                                                                                                                                                                          vars.at(iVar).cutside,
                                                                                                                                                                          FomType::kSOverSSB,
                                                                                                                                                                        vars.at(iVar).name.find("CosRej") != std::string::npos ? histTotal : histMCTotal,
                                                                                                                                                                          histMCSignal,
                                                                                                                                                                          histMCSignalNoCut,
                                                                                                                                                                          histMCTotalEnergy,
                                                                                                                                                                          histMCSignalEnergy,
                                                                                                                                                                                {});

        // style histograms
        StyleFOMs(EffTimesPurHist, SOverSBHist, BinnedSOverSSBSystHist, SOverSSBHist, 
             vars.at(iVar).name.find("CosRej") != std::string::npos ? histTotal->GetMaximum(): histMCTotal->GetMaximum());

}

void InitialiseHistograms	(	TFile *	spectraFile,
		Detector	thisDet,
		int	iCut,
		int	iVar,
		std::string	dirName,
		std::string	histPathPre,
		std::string	histPath
	)

Definition at line 396 of file PlotSpectra.C.

References update_sam_good_runs_metadata::cuts, histCosmic, histData, histMCBg, histMCSignal, histMCSignalEnergy, histMCSignalNoCut, histMCTotal, histMCTotalEnergy, histTotal, string, StyleHistograms(), and vars.

Referenced by LoopHistograms().

                                                                                       {
  // get histograms
  histMCTotal       = (TH1D*)spectraFile->Get((histPath+"TotalMC"          ).c_str());
  histMCSignal      = (TH1D*)spectraFile->Get((histPath+"SignalMC"         ).c_str());
  histMCSignalNoCut = (TH1D*)spectraFile->Get((histPathPre+"PreselSignalMC").c_str());
  histCosmic        = (TH1D*)spectraFile->Get((histPath+"TotalCosmic"      ).c_str());
  histData          = (TH1D*)spectraFile->Get((histPath+"TotalData").c_str());
  
  histMCBg = (TH1D*)histMCTotal->Clone("histMCBg");
  histMCBg->Add(histMCSignal, -1);

  histTotal = (TH1D*)histMCTotal->Clone("histTotal");
  if (histCosmic  != nullptr)
    histTotal->Add(histCosmic);

        // this pulls out TH2 of variables correlated with calorimetric energy
        std::string simEnergyString = "kCaloE"          +
                                                                                                                                vars.at(iVar).name +
                                                                                                                                cuts.at(iCut).name;

  histMCTotalEnergy = nullptr;
        TDirectoryFile* thisDir = (TDirectoryFile*)spectraFile->Get(dirName.c_str());
        if (thisDir->GetListOfKeys()->Contains((simEnergyString+"TotalMC").c_str())){
                histMCTotalEnergy = (TH2D*)spectraFile->Get((dirName+"/"+simEnergyString+"TotalMC").c_str());
        }
  histMCSignalEnergy = nullptr;
        if (thisDir->GetListOfKeys()->Contains((simEnergyString+"SignalMC").c_str())){
                histMCSignalEnergy = (TH2D*)spectraFile->Get((dirName+"/"+simEnergyString+"SignalMC").c_str());
        }

  StyleHistograms(histTotal   , thisDet, false, false, true );
  StyleHistograms(histMCTotal , thisDet, false, false, true );
  StyleHistograms(histMCSignal, thisDet, false, false, false);
  if (histCosmic != nullptr)
    StyleHistograms(histCosmic, thisDet, false, true, false);
  if (histData != nullptr)
    StyleHistograms(histData, thisDet, true, false, false);
}

void LoopHistograms	(	TFile *	spectraFile,
		Detector	thisDet
	)

Definition at line 16 of file PlotSpectra.C.

References om::cout, update_sam_good_runs_metadata::cuts, cutsFD, cutsND, dirName, DrawDataMC(), DrawFOM(), DrawIndependentSpectra(), DrawMCOnly(), allTimeWatchdog::endl, fdFractionalUncertainties, GetDecorrelatedUncertainty(), histCosmic, histMCSignal, histMCSignalNoCut, histMCTotal, MECModelEnuComparisons::i, InitialiseFOMHists(), InitialiseHistograms(), makeTrainCVSamples::int, kFar, kNear, ndFractionalUncertainties, PlotData(), PrintSelectedEvents(), string, subcuts, subCutsFD, subCutsND, vars, varsFD, and varsND.

Referenced by PlotSpectra().

                                                         {

        // load 
        TFile* covmx = new TFile("/pnfs/nova/persistent/analysis/nux/nus20/selection/covmx.root", "read");
        TMatrixD* cov = (TMatrixD*)covmx->Get("fraccovmx");

  fdFractionalUncertainties = {};
  ndFractionalUncertainties = {};
        
        int ndBins = 20;
        int fdBins = 13;
        std::vector<double> fulldecorrelatedunc = 
                GetDecorrelatedUncertainty(cov,
                                                                                                   ndBins,
                                                                                                                         fdBins);

        // separate into nd and fd
        for (size_t i = 0; i < fulldecorrelatedunc.size(); ++i){
                if ((int)i < ndBins)
                        ndFractionalUncertainties.push_back(fulldecorrelatedunc[i]);
                else
                        fdFractionalUncertainties.push_back(fulldecorrelatedunc[i]);
        }

        if ((int)ndFractionalUncertainties.size() != ndBins ||
                        (int)fdFractionalUncertainties.size() != fdBins)
                throw std::logic_error("FractionalUncertainties vector is incorrect size");

  if (thisDet == kNear){
    dirName   = "ND";
    vars    = varsND;
    cuts    = cutsND;
    subcuts = subCutsND;
  }
  else if (thisDet == kFar){
    dirName   = "FD";
    vars    = varsFD;
    cuts    = cutsFD;
    subcuts = subCutsFD;
  }
  else {
    throw std::logic_error("options for detector are nd or fd");
  }

  histMCTotal      ->Sumw2();
  histMCSignal     ->Sumw2();
  histMCSignalNoCut->Sumw2();
  histCosmic       ->Sumw2();

  // loop over every cut level and variable we want to plot
  for (size_t iCut = 0; iCut < cuts.size(); ++iCut){

    for (size_t iVar = 0; iVar < vars.size(); ++iVar){

                        std::cout << dirName << " " << vars.at(iVar).name << " " << cuts.at(iCut).name << std::endl;
      std::string histPathPre = dirName + "/" + vars.at(iVar).name;
      std::string histPath    = histPathPre + cuts.at(iCut).name;

      // figure out whether we want to plot data for this histogram
      bool isPlotData = PlotData(histPath);

      InitialiseHistograms(spectraFile, thisDet, iCut, iVar, dirName, histPathPre, histPath);

      // print POT-normalised simulation information
      PrintSelectedEvents(iCut, iVar);
     
      // define boxes
      InitialiseFOMHists(thisDet, iCut, iVar);
      DrawFOM(std::set<FomType>(), iVar); // need a vector of FomTypes if you want to draw these

      std::cout << "isPlotData: " << isPlotData << std::endl;
      // plot histograms
      if (isPlotData)
        DrawDataMC(thisDet, iCut, iVar, false, false, true); // area normalised right now
      DrawMCOnly(thisDet, iCut, iVar, false, false);
      if (thisDet == kFar){
        DrawMCOnly(thisDet, iCut, iVar, false, true);
        DrawMCOnly(thisDet, iCut, iVar, true, true);
        DrawIndependentSpectra(thisDet, iCut, iVar);
      }
    }
  }
}

bool PlotData

(

std::string

name

)

Definition at line 207 of file PlotSpectra.C.

Referenced by LoopHistograms().

                             {

  if (name.find("FD") != std::string::npos)
    return false;
  else if (
    (name.find("CVN") == std::string::npos &&
     name.find("CaloE") == std::string::npos) ||
    (name.find("CVN") != std::string::npos &&
      name.find("Fiducial") != std::string::npos) ||
    (name.find("CVN") != std::string::npos &&
      name.find("FullSelection") != std::string::npos) ||
    (name.find("CaloE") != std::string::npos &&
      name.find("FullSelection") != std::string::npos))
    return true;
  else return false;
}

void PrintSelectedEvents	(	int	iCut,
		int	iVar
	)

Definition at line 382 of file PlotSpectra.C.

References om::cout, update_sam_good_runs_metadata::cuts, allTimeWatchdog::endl, histCosmic, histMCBg, histMCSignal, histTotal, and vars.

Referenced by LoopHistograms().

{
  std::cout << "---------------------------------" << std::endl;
  std::cout << " Printing selected events" << std::endl;
  std::cout << " Stage/Var: " << cuts.at(iCut).name << "/" << vars.at(iVar).name << std::endl;
  std::cout << " -- Total Events selected: " << histTotal   ->Integral(0, histTotal->GetNbinsX()+1) << std::endl; 
  std::cout << " -- Total Signal selected: " << histMCSignal->Integral(0, histTotal->GetNbinsX()+1) << std::endl; 
  std::cout << " -- Total mc bgs selected: " << histMCBg    ->Integral(0, histTotal->GetNbinsX()+1) << std::endl;
  if (histCosmic != nullptr)
    std::cout << " -- Total Cosmic selected: " << histCosmic->Integral(0, histTotal->GetNbinsX()+1) << std::endl; 
  std::cout << "---------------------------------" << std::endl;
}

void SetupCanvas	(	TPad *	topPad,
		TPad *	botPad,
		TPad *	thisPad
	)

Definition at line 186 of file PlotSpectra.C.

References cNoRatio, and cRatio.

Referenced by PlotSpectra().

                                                           {
  cRatio->cd();
  topPad->SetTopMargin(0.1);
  topPad->SetBottomMargin(0.005);
  topPad->SetLeftMargin(0.12);
  botPad->SetTopMargin(0.005);
  botPad->SetBottomMargin(0.35);
  botPad->SetLeftMargin(0.12);
  topPad->Draw();
  botPad->Draw();

  cNoRatio->cd();
  thisPad->SetTopMargin(0.1);
  thisPad->SetBottomMargin(0.15);
  thisPad->SetLeftMargin(0.12);
  thisPad->Draw();
}

void StyleHistograms	(	TH1D *	h,
		Detector	det,
		bool	isBeamData,
		bool	isCosmic,
		bool	isTotal
	)

Definition at line 101 of file PlotSpectra.C.

References kFar, kNear, rootlogon::kNovaFont, and kOrange.

Referenced by InitialiseHistograms().

                                                                  {

  if (isBeamData){
    h->SetMarkerStyle(20);
    h->SetMarkerSize(0.6);
  }
  else {
    if (det == kNear){
      if (isTotal){
        h->SetLineColor(kBlack);
        h->SetLineWidth(2);
      }
      else {
        h->SetFillColor(kAzure+2);
        h->SetLineColor(kAzure+2);
      }
    }
    if (det == kFar){
      if (isTotal){
        h->SetLineColor(kBlack);
        h->SetLineWidth(2);
      }
      else if(isCosmic){
        h->SetFillColor(kOrange-3);
        h->SetLineColor(kOrange-3);
      }
      else {
        h->SetFillColor(kAzure+2);
        h->SetLineColor(kAzure+2);
      }
    }
    if (det == kNear){
      h->GetYaxis()->SetTitle("10^{4} Events");
    }
  }
 
  h->GetXaxis()->CenterTitle();
  h->GetXaxis()->SetTitleFont(43);
  h->GetXaxis()->SetTitleSize(28);
  h->GetYaxis()->CenterTitle();
  h->GetYaxis()->SetTitleFont(43);
  h->GetYaxis()->SetTitleSize(28);

  h->GetXaxis()->SetLabelFont(43);
  h->GetXaxis()->SetLabelSize(28);
  h->GetYaxis()->SetLabelFont(43);
  h->GetYaxis()->SetLabelSize(28);

  h->GetYaxis()->SetMaxDigits(3);
  h->GetYaxis()->SetTitleOffset(1.3);
  h->GetYaxis()->SetNdivisions(515);

  NOvACornerLabel("Neutrino Beam");
  NOvAPreliminarySide();
  if (det == kNear)
    NOvADetector("Near Detector");
  else
    NOvADetector("Far Detector");
  NOvAPOT("14.23");
  
  // Thicker lines
  gStyle->SetHistLineWidth(2);
  gStyle->SetFrameLineWidth(2);
  gStyle->SetFuncWidth(2);
 
  // Set the number of tick marks to show
  gStyle->SetNdivisions(506, "xyz");
 
  // Set the tick mark style
  gStyle->SetPadTickX(1);
  gStyle->SetPadTickY(1);
 
  // Fonts
  const int kNovaFont = 42;
  gStyle->SetStatFont(kNovaFont);
  gStyle->SetLabelFont(kNovaFont, "xyz");
  gStyle->SetTitleFont(kNovaFont, "xyz");
  gStyle->SetTitleFont(kNovaFont, ""); // Apply same setting to plot titles
  gStyle->SetTextFont(kNovaFont);
  gStyle->SetLegendFont(kNovaFont);
 
}

void StyleLegend

(

TLegend *

leg

)

Definition at line 225 of file PlotSpectra.C.

Referenced by DrawFOM(), DrawIndependentSpectra(), and GetStandardLegend().

                              {
  leg->SetLineWidth(0);
  leg->SetFillStyle(0);
  leg->SetTextFont(43);
  leg->SetTextSize(28*0.75);
}

Variable Documentation

TH1D* BinnedSOverSSBSystHist = new TH1D()

Definition at line 66 of file PlotSpectra.h.

Referenced by InitialiseFOMHists().

TPad* botPad = new TPad("botPad", "", 0.005, 0.005, 0.995, 0.3)

Definition at line 75 of file PlotSpectra.h.

Referenced by DrawDataMC(), makeRandomUniverseHistograms(), and PlotSpectra().

TCanvas* cNoRatio = new TCanvas("cNoRatio", "cNoRatio", 600, 500)

Definition at line 77 of file PlotSpectra.h.

Referenced by DrawIndependentSpectra(), DrawMCOnly(), and SetupCanvas().

TCanvas* cRatio = new TCanvas("cRatio", "cRatio", 600, 600)

Definition at line 73 of file PlotSpectra.h.

Referenced by DrawDataMC(), and SetupCanvas().

std::vector<CutContainer> cuts

Definition at line 49 of file PlotSpectra.h.

std::string dirName

Definition at line 47 of file PlotSpectra.h.

Referenced by cmf::EventList::AddEvent(), ana::AddNusAna2020FileSysts(), Analyse_Data2DataComp(), Analyse_Data2DataComp_HigherEnergyCuts(), Analyse_Data2DataComp_kNumu2020ND(), Analyse_OverlayStudies(), nova::database::Table::CreateFromDB(), nova::database::Table::CreateFromXML(), nova::dbi::Table::CreateFromXML(), cmf::EventListManipulator::Deserialize(), DrawDataMC(), DrawIndependentSpectra(), DrawMCOnly(), GetBestFitPoint(), GetContours(), GetDeltaChiSqr(), GetMedianDeltaChiSqr(), GetNuePrediction2017(), GetNuePrediction2018(), GetNuePrediction2019(), GetNuePredictionFuture(), getStabilitySpectra(), InitialiseFOMHists(), LoopHistograms(), MakeCovMx(), makeMedianContours(), runCheatDecomp(), Save1DContour(), util::SaveObj(), nutools::dbi::Table::Table(), nova::dbi::Table::Table(), test_nue2017Prediction(), and test_nue2018_fitter().

TH1D* EffTimesPurHist = new TH1D()

total data

Definition at line 63 of file PlotSpectra.h.

Referenced by DrawFOM(), InitialiseFOMHists(), PlotPreSelectionPlots(), and PlotSelectionPlots().

std::vector<float> fdFractionalUncertainties = {}

Definition at line 45 of file PlotSpectra.h.

Referenced by InitialiseFOMHists(), and LoopHistograms().

TH1D* histCosmic = new TH1D()

total mc signal (no cut)

Definition at line 57 of file PlotSpectra.h.

Referenced by DrawDataMC(), DrawIndependentSpectra(), DrawMCOnly(), InitialiseHistograms(), LoopHistograms(), and PrintSelectedEvents().

TH1D* histData = new TH1D()

mc-only bg (no cosmic)

Definition at line 59 of file PlotSpectra.h.

Referenced by DrawDataMC(), GetStandardLegend(), and InitialiseHistograms().

TH1D* histMCBg = new TH1D()

total cosmic

Definition at line 58 of file PlotSpectra.h.

Referenced by DrawDataMC(), DrawIndependentSpectra(), DrawMCOnly(), InitialiseHistograms(), and PrintSelectedEvents().

TH1D* histMCSignal = new TH1D()

total mc only

Definition at line 55 of file PlotSpectra.h.

Referenced by DrawDataMC(), DrawIndependentSpectra(), DrawMCOnly(), GetStandardLegend(), InitialiseFOMHists(), InitialiseHistograms(), LoopHistograms(), and PrintSelectedEvents().

TH2D* histMCSignalEnergy = nullptr

Definition at line 69 of file PlotSpectra.h.

Referenced by InitialiseFOMHists(), and InitialiseHistograms().

TH1D* histMCSignalNoCut = new TH1D()

total mc signal

Definition at line 56 of file PlotSpectra.h.

Referenced by InitialiseFOMHists(), InitialiseHistograms(), and LoopHistograms().

TH1D* histMCTotal = new TH1D()

total sig + bg

Definition at line 54 of file PlotSpectra.h.

Referenced by DrawDataMC(), DrawIndependentSpectra(), DrawMCOnly(), InitialiseFOMHists(), InitialiseHistograms(), and LoopHistograms().

TH2D* histMCTotalEnergy = nullptr

Definition at line 68 of file PlotSpectra.h.

Referenced by InitialiseFOMHists(), and InitialiseHistograms().

TH1D* histTotal = new TH1D()

Definition at line 53 of file PlotSpectra.h.

Referenced by DrawDataMC(), DrawIndependentSpectra(), DrawMCOnly(), energyResolution(), GetStandardLegend(), InitialiseFOMHists(), InitialiseHistograms(), and PrintSelectedEvents().

bool isAreaNorm = true

Definition at line 71 of file PlotSpectra.h.

std::vector<float> ndFractionalUncertainties = {}

Definition at line 44 of file PlotSpectra.h.

Referenced by InitialiseFOMHists(), and LoopHistograms().

TH1D* SOverSBHist = new TH1D()

Definition at line 64 of file PlotSpectra.h.

Referenced by DrawFOM(), InitialiseFOMHists(), PlotPreSelectionPlots(), and PlotSelectionPlots().

TH1D* SOverSSBHist = new TH1D()

Definition at line 65 of file PlotSpectra.h.

Referenced by DrawFOM(), InitialiseFOMHists(), PlotPreSelectionPlots(), and PlotSelectionPlots().

std::vector<CutContainer> subcuts

Definition at line 50 of file PlotSpectra.h.

Referenced by addVars(), and LoopHistograms().

TPad* thisPad = new TPad("thisPad", "", 0.005, 0.005, 0.995, 0.995)

Definition at line 78 of file PlotSpectra.h.

Referenced by DrawIndependentSpectra(), DrawMCOnly(), PlotPreSelectionPlots(), PlotSelectionPlots(), and PlotSpectra().

TPad* topPad = new TPad("topPad", "", 0.005, 0.3, 0.995, 0.995)

Definition at line 74 of file PlotSpectra.h.

Referenced by DrawDataMC(), makeCalibSystRatioPlots(), makeRandomUniverseHistograms(), PlotPreSelectionPlots(), PlotSelectionPlots(), and PlotSpectra().

std::vector<VarContainer> vars

Definition at line 48 of file PlotSpectra.h.