PlotSpectra.C File Reference

#include "NuXAna/macros/Nus20/Selection/PlotSpectra.h"

Go to the source code of this file.

Functions
void	PlotSpectra (TString spectraFileName, TString opt)
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 > theseCuts, CutSide thisCutSide)
void	DrawRatio (TPad *thisPad, TH1D *mc, TH1D *data)
void	GetFractionalUncertainties (TFile *covmx, std::vector< float > *nd, std::vector< float > *fd)
void	InitialiseFOMHists (Detector thisDet, int iVar, int iCut)
void	DrawFOM (std::set< FomType > fomTypes, int iVar)
void	PrintSelectedEvents (int iCut, int iVar)
void	InitialiseHistograms (TFile *spectraFile, Detector thisDet, int iCut, int iVar, std::string dirName, std::string histPathPre, std::string histPath)
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 ()

Function Documentation

void DrawBoxes	(	TH1D *	histMCTotal,
		std::vector< double >	theseCuts,
		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 >	fomTypes,
		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 *	data
	)

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	iVar,
		int	iCut
	)

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 PlotSpectra	(	TString	spectraFileName,
		TString	opt
	)

Definition at line 3 of file PlotSpectra.C.

References botPad, kFar, kNear, LoopHistograms(), SetupCanvas(), thisPad, and topPad.

                                                      {

  TFile* spectraFile = new TFile(spectraFileName, "read");

  SetupCanvas(topPad, botPad, thisPad);

  if (opt.Contains("nd", TString::kIgnoreCase))
    LoopHistograms(spectraFile, kNear);
  if (opt.Contains("fd", TString::kIgnoreCase))
    LoopHistograms(spectraFile, kFar);

}

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);
}