PlotXSec.C File Reference

#include "CAFAna/Core/Utilities.h"
#include "CAFAna/Analysis/Plots.h"
#include "NDAna/numucc_inc/Multiverse.h"
#include "TH1.h"
#include "TH2.h"
#include "TH3.h"
#include "TFile.h"
#include "TAxis.h"
#include "TCanvas.h"
#include "TLegend.h"
#include <string>
#include <vector>
#include <map>
#include <iostream>

Go to the source code of this file.

Functions
const TAxis	angleAxis (1, 0.5, 1.0)
const TAxis	muoneAxis (1, 0.5, 2.5)
const TAxis	uncertAxis (1,-0.2, 0.2)
const TAxis	enuAxis (1, 0.0, 5.0)
const TAxis	q2Axis (1, 0.0, 2.8)
void	FixAxisRangeQ2 (TH1 *hist)
void	FixAxisRangeENu (TH1 *hist)
void	FixAxisRange2D (TH2 *hist)
void	AddUpDwUncerts (TH1 *destHist, TH1 *upHist, TH1 *dwHist)
void	AddUpDwUncerts2D (TH2 *destHist, TH2 *upHist, TH2 *dwHist)
TH3D *	Load3DXSec (string syst)
TH2D *	XSec2D (string syst)
TH1D *	XSecENu (string syst)
TH1D *	XSecQ2 (string syst)
TH2D *	XSec2DUncert (string syst)
TH1D *	XSecENuUncert (string syst)
TH1D *	XSecQ2Uncert (string syst)
void	LoadNominal (string syst)
void	LoadSyst (string syst)
void	LoadUncerts (string syst)
void	LoadStatUncert (string nom)
void	AddQuadrature (TH1 *dest, TH1 *hist)
void	AddQuadrature (TH1 *dest, vector< TH1 * > hists)
void	AddQuadrature2D (TH2 *dest, TH2 *hist)
void	AddQuadrature2D (TH2 *dest, vector< TH2 * > hists)
TH1 *	GetMaxUpDown (TH1 *upHist, TH1 *dwHist)
TH2 *	GetMaxUpDown2D (TH2 *upHist, TH2 *dwHist)
TH1 *	GetAbsolute (TH1 *hist)
TH2 *	GetAbsolute2D (TH2 *hist)
TH1 *	GetSliceHist (TH2 *hist, int ibinx)
TH2 *	GetTotalBeamUncert2D ()
TH1 *	GetTotalBeamUncertENu ()
TH1 *	GetTotalBeamUncertQ2 ()
void	LimitTMuBins (TH1 *hist, int ibiny)
void	MakeSystematicPlots (vector< string > systs, string systFilename)
void	PlotXSec (string filename, string outputFolder)

Variables
map< string, string >	systTitles
vector< string >	standardSysts
map< string, pair< string, string > >	standardSystsUpDown
const map< string, string >	otherNomSysts
vector< string >	beamSysts
map< string, Color_t >	systColors
map< int, double >	maxTMuBins
const string	var3DName = "EAvail"
string	uncertLabel = "_uncert"
const vector< string >	vars3D
TH3D *	nom3DXSec
TH2D *	nom2DXSec
string	defaultNominal = "cv"
map< string, TH3D * >	nominalXSecs3D
map< string, TH2D * >	nominalXSecs2D
map< string, TH1D * >	nominalXSecsENu
map< string, TH1D * >	nominalXSecsQ2
map< string, TH3D * >	all3DXSecs
map< string, TH2D * >	xSecs2D
map< string, TH1D * >	xSecsENu
map< string, TH1D * >	xSecsQ2
map< string, TH2D * >	uncerts2D
map< string, TH1D * >	uncertsENu
map< string, TH1D * >	uncertsQ2
TFile *	inputFile
TFile *	outFile
TCanvas *	c1
const unsigned int	N_GENIE = 1000
const unsigned int	N_PPFX = 100

Function Documentation

void AddQuadrature	(	TH1 *	dest,
		TH1 *	hist
	)

Definition at line 437 of file PlotXSec.C.

References ana::assert(), make_syst_table_plots::ibin, cet::pow(), and std::sqrt().

Referenced by AddQuadrature(), and MakeSystematicPlots().

{
  assert(dest && hist && dest->GetNbinsX() == hist->GetNbinsX());
  for(int ibin = 0; ibin <= hist->GetNbinsX() + 1; ++ibin)
    dest->SetBinContent(ibin, sqrt(pow(dest->GetBinContent(ibin), 2) + pow(hist->GetBinContent(ibin), 2)));
}

void AddQuadrature	(	TH1 *	dest,
		vector< TH1 * >	hists
	)

Definition at line 444 of file PlotXSec.C.

References AddQuadrature(), ana::assert(), and analysePickle::hist.

{
  assert(!hists.empty());
  for (TH1 * hist : hists)
    AddQuadrature(dest, hist);
}

void AddQuadrature2D	(	TH2 *	dest,
		TH2 *	hist
	)

Definition at line 451 of file PlotXSec.C.

References ana::assert(), cet::pow(), and std::sqrt().

Referenced by AddQuadrature2D(), and MakeSystematicPlots().

{
  assert(dest && hist && dest->GetNbinsX() == hist->GetNbinsX() && dest->GetNbinsY() == hist->GetNbinsY());
  for(int ibinx = 0; ibinx <= dest->GetNbinsX() + 1; ++ibinx)
    for(int ibiny = 0; ibiny <= dest->GetNbinsY() + 1; ++ibiny)
      dest->SetBinContent(ibinx, ibiny, sqrt(pow(dest->GetBinContent(ibinx, ibiny), 2) + pow(hist->GetBinContent(ibinx, ibiny), 2)));
  // for(int ibinx = 0; ibinx <= dest->GetNbinsX() + 1; ++ibinx)
  //   for (int ibiny = 0; ibiny <= dest->GetNbinsY() + 1; ++ibiny)
  //   {
  //     double addErr = hist->GetBinContent(ibinx, ibiny);
  //     double currentErr = dest->GetBinContent(ibinx, ibiny);
  //     dest->SetBinContent(ibinx, ibiny, sqrt( pow(currentErr, 2) + pow(addErr, 2)));
  //   }
}

void AddQuadrature2D	(	TH2 *	dest,
		vector< TH2 * >	hists
	)

Definition at line 466 of file PlotXSec.C.

References AddQuadrature2D(), ana::assert(), and analysePickle::hist.

{
  assert(!hists.empty());
  for (TH2 * hist : hists)
    AddQuadrature2D(dest, hist);
}

void AddUpDwUncerts	(	TH1 *	destHist,
		TH1 *	upHist,
		TH1 *	dwHist
	)

Definition at line 177 of file PlotXSec.C.

References abs(), ana::assert(), cet::sqlite::max(), cet::pow(), and std::sqrt().

Referenced by GetTotalBeamUncertENu(), and GetTotalBeamUncertQ2().

{
  assert(destHist->GetNbinsX() == upHist->GetNbinsX() && destHist->GetNbinsX() == dwHist->GetNbinsX());
  for(int ibinx = 0; ibinx <= destHist->GetNbinsX() + 1; ++ibinx)
  {
    double upErr = abs(upHist->GetBinContent(ibinx));
    double dwErr = abs(dwHist->GetBinContent(ibinx));
    double addErr = max(upErr, dwErr);
    double currentErr = destHist->GetBinContent(ibinx);
    destHist->SetBinContent(ibinx, sqrt( pow(currentErr, 2) + pow(addErr, 2)));
  }
}

void AddUpDwUncerts2D	(	TH2 *	destHist,
		TH2 *	upHist,
		TH2 *	dwHist
	)

Definition at line 190 of file PlotXSec.C.

References abs(), ana::assert(), cet::sqlite::max(), cet::pow(), and std::sqrt().

Referenced by GetTotalBeamUncert2D().

{
  assert(destHist->GetNbinsX() == upHist->GetNbinsX() && destHist->GetNbinsX() == dwHist->GetNbinsX());
  assert(destHist->GetNbinsY() == upHist->GetNbinsY() && destHist->GetNbinsY() == dwHist->GetNbinsY());
  for(int ibinx = 0; ibinx <= destHist->GetNbinsX() + 1; ++ibinx)
    for (int ibiny = 0; ibiny <= destHist->GetNbinsY() + 1; ++ibiny)
    {
      double upErr = abs(upHist->GetBinContent(ibinx, ibiny));
      double dwErr = abs(dwHist->GetBinContent(ibinx, ibiny));
      double addErr = max(upErr, dwErr);
      double currentErr = destHist->GetBinContent(ibinx, ibiny);
      destHist->SetBinContent(ibinx, ibiny, sqrt( pow(currentErr, 2) + pow(addErr, 2)));
    }
}

const TAxis angleAxis	(	1	,
		0.	5,
		1.	0
	)

Referenced by FixAxisRange2D().

const TAxis enuAxis	(	1	,
		0.	0,
		5.	0
	)

Referenced by FixAxisRangeENu().

void FixAxisRange2D

(

TH2 *

hist

)

Definition at line 168 of file PlotXSec.C.

References angleAxis(), and muoneAxis().

Referenced by XSec2D(), and XSec2DUncert().

{
  hist->GetXaxis()->SetRangeUser(angleAxis.GetXmin(), angleAxis.GetXmax());
  hist->GetYaxis()->SetRangeUser(muoneAxis.GetXmin(), muoneAxis.GetXmax());
  // string histName = hist->GetName();
  // if (histName.find(uncertLabel) != string::npos)
  //   hist->GetZaxis()->SetRangeUser(uncertAxis.GetXmin(), uncertAxis.GetXmax());
}

void FixAxisRangeENu

(

TH1 *

hist

)

Definition at line 159 of file PlotXSec.C.

References enuAxis().

Referenced by XSecENu(), and XSecENuUncert().

{
  hist->GetXaxis()->SetRangeUser(enuAxis.GetXmin(), enuAxis.GetXmax());
  // string histName = hist->GetName();
  // if (histName.find(uncertLabel) != string::npos)
  //   hist->GetYaxis()->SetRangeUser(uncertAxis.GetXmin(), uncertAxis.GetXmax());
}

void FixAxisRangeQ2

(

TH1 *

hist

)

Definition at line 151 of file PlotXSec.C.

References q2Axis().

Referenced by XSecQ2(), and XSecQ2Uncert().

{
  hist->GetXaxis()->SetRangeUser(q2Axis.GetXmin(), q2Axis.GetXmax());
  // string histName = hist->GetName();
  // if (histName.find(uncertLabel) != string::npos)
    // hist->GetYaxis()->SetRangeUser(uncertAxis.GetXmin(), uncertAxis.GetXmax());
}

TH1* GetAbsolute

(

TH1 *

hist

)

Definition at line 493 of file PlotXSec.C.

References abs(), make_syst_table_plots::ibin, and ana::UniqueName().

Referenced by GetSliceHist(), and MakeSystematicPlots().

{
  TH1 * histAbs = (TH1*) hist->Clone(ana::UniqueName().c_str());
  for(int ibin = 0; ibin <= hist->GetNbinsX(); ++ibin)
    histAbs->SetBinContent(ibin, abs(hist->GetBinContent(ibin)));
  return histAbs;
}

TH2* GetAbsolute2D

(

TH2 *

hist

)

Definition at line 502 of file PlotXSec.C.

References abs(), and ana::UniqueName().

Referenced by MakeSystematicPlots().

{
  TH2 * histAbs = (TH2*) hist->Clone(ana::UniqueName().c_str());
  for(int ibinx = 0; ibinx <= hist->GetNbinsX(); ++ibinx)
    for(int ibiny = 0; ibiny <= hist->GetNbinsY(); ++ibiny)
      histAbs->SetBinContent(ibinx, ibiny, abs(hist->GetBinContent(ibinx, ibiny)));
  return histAbs;
}

TH1* GetMaxUpDown	(	TH1 *	upHist,
		TH1 *	dwHist
	)

Definition at line 473 of file PlotXSec.C.

References abs(), make_syst_table_plots::ibin, cet::sqlite::max(), and ana::UniqueName().

Referenced by MakeSystematicPlots().

{
  TH1 * upDwError = (TH1*) upHist->Clone(ana::UniqueName().c_str());
  upDwError->Reset();
  for (int ibin = 0; ibin <= upHist->GetNbinsX() + 1; ++ibin)
    upDwError->SetBinContent(ibin, max(abs(upHist->GetBinContent(ibin)), abs(dwHist->GetBinContent(ibin))));
  return upDwError;
}

TH2* GetMaxUpDown2D	(	TH2 *	upHist,
		TH2 *	dwHist
	)

Definition at line 482 of file PlotXSec.C.

References abs(), cet::sqlite::max(), and ana::UniqueName().

Referenced by MakeSystematicPlots().

{
  TH2 * upDwError = (TH2*) upHist->Clone(ana::UniqueName().c_str());
  upDwError->Reset();
  for (int ibinx = 0; ibinx <= upHist->GetNbinsX() + 1; ++ibinx)
    for (int ibiny = 0; ibiny <= upHist->GetNbinsY() + 1; ++ibiny)
      upDwError->SetBinContent(ibinx, ibiny, max(abs(upHist->GetBinContent(ibinx, ibiny)), abs(dwHist->GetBinContent(ibinx, ibiny))));
  return upDwError;
}

TH1* GetSliceHist	(	TH2 *	hist,
		int	ibinx
	)

Definition at line 512 of file PlotXSec.C.

References GetAbsolute(), string, art::to_string(), and ana::UniqueName().

Referenced by MakeSystematicPlots().

{
  // TH1 * histSlice = hist->ProjectionY((hist->GetName() + string("_") + to_string(ibinx)).c_str(), ibinx, ibinx);
  TH1 * histSlice = hist->ProjectionY((ana::UniqueName()).c_str(), ibinx, ibinx);
  TH1 * absSlice  = GetAbsolute(histSlice);
  absSlice->SetName((histSlice->GetName() + string("_abs")).c_str());
  absSlice->SetTitle((to_string(histSlice->GetXaxis()->GetXmin()) + " < cos(#theta_{#mu}) < " + to_string(histSlice->GetXaxis()->GetXmax())).c_str());
  return absSlice;
}

TH2* GetTotalBeamUncert2D

(

)

Definition at line 523 of file PlotXSec.C.

References AddUpDwUncerts2D(), beamSysts, plot_validation_datamc::Clone(), nominalXSecs2D, and uncerts2D.

Referenced by MakeSystematicPlots().

{
  TH2D * totalBeamUncert = (TH2D*) nominalXSecs2D["cv"]->Clone("totalBeamUncert2D");
  totalBeamUncert->Reset();
  for (string beamSyst : beamSysts)
  {
    TH2 * upSyst = uncerts2D[beamSyst + "_Up"];
    TH2 * dwSyst = uncerts2D[beamSyst + "_Dw"];
    AddUpDwUncerts2D(totalBeamUncert, upSyst, dwSyst);
  }
  uncerts2D["totalBeam"] = totalBeamUncert;
  return totalBeamUncert;
}

TH1* GetTotalBeamUncertENu

(

)

Definition at line 537 of file PlotXSec.C.

References AddUpDwUncerts(), beamSysts, plot_validation_datamc::Clone(), nominalXSecsENu, and uncertsENu.

Referenced by MakeSystematicPlots().

{
  TH1D * totalBeamUncert = (TH1D*) nominalXSecsENu["cv"]->Clone("totalBeamUncertENu");
  totalBeamUncert->Reset();
  for (string beamSyst : beamSysts)
  {
    TH1 * upSyst = uncertsENu[beamSyst + "_Up"];
    TH1 * dwSyst = uncertsENu[beamSyst + "_Dw"];
    AddUpDwUncerts(totalBeamUncert, upSyst, dwSyst);
  }
  uncertsENu["totalBeam"] = totalBeamUncert;
  return totalBeamUncert;
}

TH1* GetTotalBeamUncertQ2

(

)

Definition at line 551 of file PlotXSec.C.

References AddUpDwUncerts(), beamSysts, plot_validation_datamc::Clone(), nominalXSecsQ2, and uncertsQ2.

Referenced by MakeSystematicPlots().

{
  TH1D * totalBeamUncert = (TH1D*) nominalXSecsQ2["cv"]->Clone("totalBeamUncertQ2");
  totalBeamUncert->Reset();
  for (string beamSyst : beamSysts)
  {
    TH1 * upSyst = uncertsQ2[beamSyst + "_Up"];
    TH1 * dwSyst = uncertsQ2[beamSyst + "_Dw"];
    AddUpDwUncerts(totalBeamUncert, upSyst, dwSyst);
  }
  uncertsQ2["totalBeam"] = totalBeamUncert;
  return totalBeamUncert;
}

void LimitTMuBins	(	TH1 *	hist,
		int	ibiny
	)

Definition at line 566 of file PlotXSec.C.

References make_syst_table_plots::ibin, and maxTMuBins.

Referenced by MakeSystematicPlots().

{
  if (maxTMuBins.find(ibiny) != maxTMuBins.end())
    for (int ibin = 1; ibin <= hist->GetNbinsX(); ++ibin)
      if (hist->GetXaxis()->GetBinCenter(ibin) >= maxTMuBins.at(ibiny)){
        hist->SetBinContent(ibin, 0);
        hist->SetBinError(ibin, 0);
      }
}

TH3D* Load3DXSec

(

string

syst

)

Definition at line 206 of file PlotXSec.C.

References all3DXSecs, inputFile, systTitles, and var3DName.

Referenced by LoadNominal(), LoadSyst(), and XSec2D().

{
  if (all3DXSecs.find(syst) != all3DXSecs.end())
    return all3DXSecs.at(syst);

  TH3D* xsec3D = (TH3D*) inputFile->Get((syst + "/" + var3DName + "/xsec_" + syst + "_" + var3DName).c_str());
  all3DXSecs[syst] = xsec3D;
  if (!xsec3D)
    return NULL;
  xsec3D->SetTitle((systTitles[syst] + " XSec").c_str());
  xsec3D->Write();
  return xsec3D;
}

void LoadNominal

(

string

syst

)

Definition at line 348 of file PlotXSec.C.

References dir, Load3DXSec(), nominalXSecs2D, nominalXSecs3D, nominalXSecsENu, nominalXSecsQ2, outFile, XSec2D(), XSecENu(), and XSecQ2().

Referenced by PlotXSec().

{
  if (nominalXSecs3D.find(syst) != nominalXSecs3D.end())
    return;

  TDirectory * dir = outFile->GetDirectory(syst.c_str());
  if (!dir)
    dir = outFile->mkdir(syst.c_str());
  dir->cd();

  TH3D * nomXSec    = Load3DXSec(syst);
  TH2D * nomXSec2D  = XSec2D(syst);
  TH1D * nomXSecENu = XSecENu(syst);
  TH1D * nomXSecQ2  = XSecQ2(syst);

  nominalXSecs3D[syst]  = nomXSec;
  nominalXSecs2D[syst]  = nomXSec2D;
  nominalXSecsENu[syst] = nomXSecENu;
  nominalXSecsQ2 [syst] = nomXSecQ2;

  outFile->cd();
}

void LoadStatUncert

(

string

nom

)

Definition at line 398 of file PlotXSec.C.

References plot_validation_datamc::Clone(), nominalXSecs2D, nominalXSecsENu, nominalXSecsQ2, uncerts2D, uncertsENu, and uncertsQ2.

Referenced by PlotXSec().

{
  if (uncerts2D. find("stat") != uncerts2D.end())
    return; // Already loaded

  // Create Stat Uncerts
  TH2D * statUncert2D = NULL;
  if (nominalXSecs2D.find(nom) != nominalXSecs2D.end() && nominalXSecs2D[nom])
    statUncert2D = (TH2D*) nominalXSecs2D[nom]->Clone("stat_uncert_2D");
  for (int ibinx = 0; statUncert2D && ibinx <= statUncert2D->GetNbinsX() + 1; ++ibinx)
    for (int ibiny = 0; statUncert2D && ibiny <= statUncert2D->GetNbinsY() + 1; ++ibiny){
      double content = statUncert2D->GetBinContent(ibinx, ibiny);
      double uncert  = statUncert2D->GetBinError  (ibinx, ibiny);
      statUncert2D->SetTitle("Stat Uncert;#cos{#theta_{#mu}};T_{#mu};;");
      statUncert2D->SetBinContent(ibinx, ibiny,  content == 0.0 ? 0.0 : uncert / content);
      statUncert2D->SetBinError(ibinx, ibiny, 0);
    }
  TH1D * statUncertENu = (TH1D*) nominalXSecsENu[nom]->Clone("stat_uncert_ENu");
  for (int ibinx = 0; ibinx <= statUncertENu->GetNbinsX() + 1; ++ibinx){
    double content = statUncertENu->GetBinContent(ibinx);
    double uncert  = statUncertENu->GetBinError(ibinx);
    statUncertENu->SetTitle("Statistical Uncertainty;E_{nu} (GeV);;;");
    statUncertENu->SetBinContent(ibinx, content == 0.0 ? 0.0 : uncert / content);
    statUncertENu->SetBinError(ibinx, 0);
  }
  TH1D * statUncertQ2 = (TH1D*) nominalXSecsQ2[nom]->Clone("stat_uncert_Q2");
  for (int ibinx = 0; ibinx <= statUncertQ2->GetNbinsX() + 1; ++ibinx){
    double content = statUncertQ2->GetBinContent(ibinx);
    double uncert  = statUncertQ2->GetBinError(ibinx);
    statUncertQ2->SetTitle("Statistical Uncertainty;Q^{2} (GeV^{2});;;");
    statUncertQ2->SetBinContent(ibinx, content == 0.0 ? 0.0 : uncert / content);
    statUncertQ2->SetBinError(ibinx, 0);
  }
  if (statUncert2D)
    uncerts2D ["stat"] = statUncert2D;
  uncertsENu["stat"] = statUncertENu;
  uncertsQ2 ["stat"] = statUncertQ2;
}

void LoadSyst

(

string

syst

)

Definition at line 371 of file PlotXSec.C.

References dir, Load3DXSec(), outFile, XSec2D(), XSecENu(), and XSecQ2().

Referenced by PlotXSec().

{
  TDirectory * dir = outFile->GetDirectory(syst.c_str());
  if (!dir)
    dir = outFile->mkdir(syst.c_str());
  dir->cd();
  TH3D * xSec3D  = Load3DXSec(syst);
  TH2D * xSec2D  = XSec2D(syst);
  TH1D * xSecENu = XSecENu(syst);
  TH1D * xSecQ2  = XSecQ2(syst);

  outFile->cd();
}

void LoadUncerts

(

string

syst

)

Definition at line 385 of file PlotXSec.C.

References dir, outFile, XSec2DUncert(), XSecENuUncert(), and XSecQ2Uncert().

Referenced by PlotXSec().

{
  TDirectory * dir = outFile->GetDirectory(syst.c_str());
  if (!dir)
    dir = outFile->mkdir(syst.c_str());
  dir->cd();
  TH2 * uncert2D  = XSec2DUncert(syst);
  TH1 * uncertENu = XSecENuUncert(syst);
  TH1 * uncertQ2  = XSecQ2Uncert(syst);

  outFile->cd();
}

void MakeSystematicPlots	(	vector< string >	systs,
		string	systFilename
	)

Definition at line 579 of file PlotXSec.C.

References AddQuadrature(), AddQuadrature2D(), ana::AutoPlaceLegend(), plot_validation_datamc::Clone(), om::cout, allTimeWatchdog::endl, GetAbsolute(), GetAbsolute2D(), GetMaxUpDown(), GetMaxUpDown2D(), ana::Multiverse::GetMinusOneSigmaShift(), ana::Multiverse::GetPlusOneSigmaShift(), GetSliceHist(), GetTotalBeamUncert2D(), GetTotalBeamUncertENu(), GetTotalBeamUncertQ2(), analysePickle::hist, hists, MECModelEnuComparisons::leg, LimitTMuBins(), N_GENIE, N_PPFX, nominalXSecs2D, nominalXSecsENu, nominalXSecsQ2, push_back(), standardSysts, standardSystsUpDown, systColors, systTitles, art::to_string(), uncerts2D, uncertsENu, uncertsQ2, xSecs2D, xSecsENu, and xSecsQ2.

Referenced by PlotXSec().

{
  TH2 * totalUncert2D = NULL;
  if (nominalXSecs2D.find("cv") != nominalXSecs2D.end() && nominalXSecs2D["cv"])
    totalUncert2D = (TH2*) nominalXSecs2D["cv"] ->Clone("totalUncert2D");
  TH1 * totalUncertENu = (TH1*) nominalXSecsENu["cv"]->Clone("totalUncertENu");
  TH1 * totalUncertQ2  = (TH1*) nominalXSecsQ2["cv"] ->Clone("totalUncertQ2");

  // Switch the 
  if (totalUncert2D)
    totalUncert2D->Reset();
  totalUncertENu ->Reset();
  totalUncertQ2 ->Reset();

  map<string, TH2*> uncert2DPlots;
  map<string, TH1*> uncertENuPlots;
  map<string, TH1*> uncertQ2Plots;

  for (string syst : systs)
  {
    TH2 * uncert2D  = NULL;
    TH1 * uncertENu = NULL;
    TH1 * uncertQ2  = NULL;
    if(std::find(standardSysts.begin(), standardSysts.end(), syst) != standardSysts.end() || syst == "stat")
    {
      if (uncerts2D.find(syst) != uncerts2D.end() && uncerts2D[syst]){
        cout << "Found 2D uncert for: " << syst << endl;
        uncert2D  = GetAbsolute2D(uncerts2D[syst]);
      }
      uncertENu = GetAbsolute(uncertsENu[syst]);
      uncertQ2  = GetAbsolute(uncertsQ2[syst]);
    }
    else if (standardSystsUpDown.find(syst) != standardSystsUpDown.end())
    {
      string upSyst = standardSystsUpDown[syst].first;
      string dwSyst = standardSystsUpDown[syst].second;
      if (uncerts2D.find(upSyst) != uncerts2D.end() && uncerts2D.find(dwSyst) != uncerts2D.end() && uncerts2D[upSyst] && uncerts2D[dwSyst])
        uncert2D  = GetMaxUpDown2D(uncerts2D[upSyst], uncerts2D[dwSyst]);
      uncertENu = GetMaxUpDown(uncertsENu[upSyst], uncertsENu[dwSyst]);
      uncertQ2  = GetMaxUpDown(uncertsQ2[upSyst], uncertsQ2[dwSyst]);
    }
    else if (syst == "beam")
    {
      if (uncerts2D.find("2kA_Up") != uncerts2D.end() && uncerts2D["2kA_Up"])
        uncert2D  = GetTotalBeamUncert2D();
      uncertENu = GetTotalBeamUncertENu();
      uncertQ2  = GetTotalBeamUncertQ2();
    }
    else if (syst == "genie")
    {
      vector<TH1*> genieUnivs2D;
      vector<TH1*> genieUnivsENu;
      vector<TH1*> genieUnivsQ2;
      for (unsigned int iuniv = 0; iuniv < N_GENIE; ++iuniv)
      {
        string syst = "genie" + to_string(iuniv);
        if(xSecs2D.find(syst) != xSecs2D.end() && uncerts2D[syst])
          genieUnivs2D. push_back(xSecs2D[syst]);
        genieUnivsENu.push_back(xSecsENu[syst]);
        genieUnivsQ2. push_back(xSecsQ2[syst]);
      }
      ana::Multiverse * genieVerse2D = NULL;
      if (!genieUnivs2D.empty())
        genieVerse2D = new ana::Multiverse(genieUnivs2D,  2);
      ana::Multiverse * genieVerseENu = new ana::Multiverse(genieUnivsENu, 1);
      ana::Multiverse * genieVerseQ2  = new ana::Multiverse(genieUnivsQ2,  1);
      if (genieVerse2D)
        uncert2D = GetMaxUpDown2D(
          new TH2D(*(TH2D*)genieVerse2D->GetPlusOneSigmaShift(nominalXSecs2D["cv"]) - *nominalXSecs2D["cv"]), 
          new TH2D(*(TH2D*)genieVerse2D->GetMinusOneSigmaShift(nominalXSecs2D["cv"]) - *nominalXSecs2D["cv"]));
      uncertENu = GetMaxUpDown(
        new TH1D(*(TH1D*)genieVerseENu->GetPlusOneSigmaShift(nominalXSecsENu["cv"]) - *nominalXSecsENu["cv"]), 
        new TH1D(*(TH1D*)genieVerseENu->GetMinusOneSigmaShift(nominalXSecsENu["cv"]) - *nominalXSecsENu["cv"]));
      uncertQ2 = GetMaxUpDown(
        new TH1D(*(TH1D*)genieVerseQ2->GetPlusOneSigmaShift(nominalXSecsQ2["cv"]) - *nominalXSecsQ2["cv"]), 
        new TH1D(*(TH1D*)genieVerseQ2->GetMinusOneSigmaShift(nominalXSecsQ2["cv"]) - *nominalXSecsQ2["cv"]));

      if (uncert2D)
        uncert2D->Divide(nominalXSecs2D["cv"]);
      uncertENu->Divide(nominalXSecsENu["cv"]);
      uncertQ2->Divide(nominalXSecsQ2["cv"]);
    }
    else if (syst == "ppfx")
    {
      vector<TH1*> ppfxUnivs2D;
      vector<TH1*> ppfxUnivsENu;
      vector<TH1*> ppfxUnivsQ2;
      for (unsigned int iuniv = 0; iuniv < N_PPFX; ++iuniv)
      {
        string syst = "ppfx" + to_string(iuniv);
        if(xSecs2D.find(syst) != xSecs2D.end() && uncerts2D[syst])
          ppfxUnivs2D. push_back(xSecs2D[syst]);
        ppfxUnivsENu.push_back(xSecsENu[syst]);
        ppfxUnivsQ2. push_back(xSecsQ2[syst]);
      }
      ana::Multiverse * ppfxVerse2D = NULL;
      if (!ppfxUnivs2D.empty())
        ppfxVerse2D = new ana::Multiverse(ppfxUnivs2D,  2);
      ana::Multiverse * ppfxVerseENu = new ana::Multiverse(ppfxUnivsENu, 1);
      ana::Multiverse * ppfxVerseQ2  = new ana::Multiverse(ppfxUnivsQ2,  1);
      if (ppfxVerse2D)
        uncert2D = GetMaxUpDown2D(
          new TH2D(*(TH2D*)ppfxVerse2D->GetPlusOneSigmaShift(nominalXSecs2D["cv"]) - *nominalXSecs2D["cv"]), 
          new TH2D(*(TH2D*)ppfxVerse2D->GetMinusOneSigmaShift(nominalXSecs2D["cv"]) - *nominalXSecs2D["cv"]));
      uncertENu = GetMaxUpDown(
        new TH1D(*(TH1D*)ppfxVerseENu->GetPlusOneSigmaShift(nominalXSecsENu["cv"]) - *nominalXSecsENu["cv"]), 
        new TH1D(*(TH1D*)ppfxVerseENu->GetMinusOneSigmaShift(nominalXSecsENu["cv"]) - *nominalXSecsENu["cv"]));
      uncertQ2 = GetMaxUpDown(
        new TH1D(*(TH1D*)ppfxVerseQ2->GetPlusOneSigmaShift(nominalXSecsQ2["cv"]) - *nominalXSecsQ2["cv"]), 
        new TH1D(*(TH1D*)ppfxVerseQ2->GetMinusOneSigmaShift(nominalXSecsQ2["cv"]) - *nominalXSecsQ2["cv"]));

      if (uncert2D)
        uncert2D->Divide(nominalXSecs2D["cv"]);
      uncertENu->Divide(nominalXSecsENu["cv"]);
      uncertQ2->Divide(nominalXSecsQ2["cv"]);
    }
    if (totalUncert2D && uncert2D)
      AddQuadrature2D(totalUncert2D, uncert2D);
    AddQuadrature(totalUncertENu, uncertENu);
    AddQuadrature(totalUncertQ2, uncertQ2);
    if (uncert2D)  uncert2D ->SetLineColor(systColors[syst]);
    if (uncertENu) uncertENu->SetLineColor(systColors[syst]);
    if (uncertQ2)  uncertQ2 ->SetLineColor(systColors[syst]);
    uncert2DPlots[syst]  = uncert2D;
    uncertENuPlots[syst] = uncertENu;
    uncertQ2Plots[syst]  = uncertQ2;
  }

  // Plotting
  TCanvas * systCanv = new TCanvas();
  systCanv->SetLeftMargin(0.12);
  systCanv->SetBottomMargin(0.12);
  systCanv->SaveAs((systFilename + "[").c_str());
  TLegend * leg = new TLegend();

  // 2D Slices
  if (totalUncert2D){
    for (int ibinx = 2; ibinx <= totalUncert2D->GetNbinsX(); ++ibinx)
    {
      TH1 * totalUncertSlice = GetSliceHist(totalUncert2D, ibinx);
      LimitTMuBins(totalUncertSlice, ibinx);
      totalUncertSlice->SetTitle(("Uncertainties: " + 
        to_string(totalUncert2D->GetXaxis()->GetBinLowEdge(ibinx)) + " < cos(#theta_{#mu}) < " + 
        to_string(totalUncert2D->GetXaxis()->GetBinUpEdge(ibinx))).c_str());
      totalUncertSlice->GetYaxis()->SetRangeUser(0.0, 0.38);
      totalUncertSlice->Draw("hist");
      vector<pair<TH1*, string> > histsForLeg;
      for (string syst : systs){
        TH1 * histSlice = GetSliceHist(uncert2DPlots[syst], ibinx);
        LimitTMuBins(histSlice, ibinx);
        histSlice->SetLineColor(systColors[syst]);
        histSlice->Draw("hist same");
        histsForLeg.push_back({histSlice, systTitles[syst]});
      }
      leg = ana::AutoPlaceLegend(0.2, 0.35);
      leg->AddEntry(totalUncertENu, "Total Uncert");
      for (pair<TH1*, string> hists : histsForLeg)
        leg->AddEntry(hists.first, hists.second.c_str());
      leg->Draw();
      systCanv->SaveAs(systFilename.c_str());
      systCanv->SaveAs(("plots/SliceSyst_" + to_string(ibinx) + ".png").c_str());
      leg->Clear();
      delete leg;
    }
  }
  else{
    // ENu Plots
    totalUncertENu->SetTitle("Uncertainties - ENu");
    totalUncertENu->GetYaxis()->SetTitle("#frac{#delta #sigma}{#sigma}");
    totalUncertENu->GetYaxis()->CenterTitle();
    totalUncertENu->Draw("hist");
    for (string syst : systs){
      TH1 * hist = uncertENuPlots[syst];
      hist->SetLineColor(systColors[syst]);
      hist->Draw("hist same");
    }
    leg = ana::AutoPlaceLegend(0.2, 0.35);
    leg->AddEntry(totalUncertENu, "Total Uncert");
    for (string syst : systs)
      leg->AddEntry(uncertENuPlots[syst], systTitles[syst].c_str());
    leg->Draw();
    systCanv->SaveAs(systFilename.c_str());
    leg->Clear();
    delete leg;

    // Q2 Plots
    totalUncertQ2->SetTitle("Uncertainties - Q2");
    totalUncertQ2->GetYaxis()->SetTitle("#frac{#delta #sigma}{#sigma}");
    totalUncertQ2->GetYaxis()->CenterTitle();
    totalUncertQ2->Draw("hist");
    for (string syst : systs){
      TH1 * hist = uncertQ2Plots[syst];
      hist->SetLineColor(systColors[syst]);
      hist->Draw("hist same");
    }
    leg = ana::AutoPlaceLegend(0.2, 0.35);
    leg->AddEntry(totalUncertQ2, "Total Uncert");
    for (string syst : systs)
      leg->AddEntry(uncertQ2Plots[syst], systTitles[syst].c_str());
    leg->Draw();
    systCanv->SaveAs(systFilename.c_str());
    leg->Clear();
    delete leg;
  }

  systCanv->SaveAs((systFilename + "]").c_str());

  systCanv->Clear();
  delete systCanv;
}

const TAxis muoneAxis	(	1	,
		0.	5,
		2.	5
	)

Referenced by FixAxisRange2D().

void PlotXSec	(	string	filename,
		string	outputFolder
	)

Definition at line 796 of file PlotXSec.C.

References beamSysts, c1, make_mec_shifts_plots::GetMaximum(), analysePickle::hist, inputFile, kBlue, kRed, LoadNominal(), LoadStatUncert(), LoadSyst(), LoadUncerts(), MakeSystematicPlots(), N_GENIE, N_PPFX, nominalXSecs2D, nominalXSecsENu, nominalXSecsQ2, otherNomSysts, outFile, standardSysts, standardSystsUpDown, art::to_string(), xSecsENu, and xSecsQ2.

{
  inputFile = new TFile(filename.c_str(), "READ");
  outFile = new TFile((outputFolder + "/xsecUncertResults.root").c_str(), "RECREATE");

  c1 = new TCanvas();
  string outputPDF = outputFolder + "/allUncertainties.pdf";
  c1->SetTopMargin(0.15);
  c1->SetBottomMargin(0.15);
  c1->SetRightMargin(0.15);
  c1->SaveAs((outputPDF + "[").c_str());

  string output1D = outputFolder + "/1DXSecs.pdf";

  string doubleDiffTitle = "#frac{d^{2}#sigma}{dcos(#theta)dTmu} (cm^{2} / GeV / nucleon)";
  string uncertTitle = "#frac{#delta#sigma}{#sigma}";

  LoadNominal("cv");
  for(pair<string, string> nonNominal : otherNomSysts)
  {
    string otherNom = nonNominal.second;
    if (otherNomSysts.find(otherNom) == otherNomSysts.end())
      LoadNominal(otherNom); // Don't repeat for the same systematic
  }
  if (nominalXSecs2D["cv"]){
    nominalXSecs2D["cv"]->Draw("COLZ");
    c1->SaveAs(outputPDF.c_str());
  }
  LoadStatUncert("cv");

  // Standard XSecs
  for (string syst : standardSysts)
  {
    LoadSyst(syst);
    LoadUncerts(syst);
  }

  // Standard XSecs
  for (pair<string, pair<string, string> > systDef : standardSystsUpDown)
  {
    string systBase = systDef.first;
    string systUp = systDef.second.first;
    string systDw = systDef.second.second;

    LoadSyst(systUp);
    LoadUncerts(systUp);

    LoadSyst(systDw);
    LoadUncerts(systDw);
  }

  // Beam XSecs
  for (string beamSyst : beamSysts)
  {
    string systUp = beamSyst + "_Up";
    string systDw = beamSyst + "_Dw";

    LoadSyst(systUp);
    LoadUncerts(systUp);

    LoadSyst(systDw);
    LoadUncerts(systDw);
  }

  // Handle Genie and PPFX
  for (unsigned int igenie = 0; igenie < N_GENIE; ++igenie)
  {
    string syst = "genie" + to_string(igenie);
    LoadSyst(syst);
    LoadUncerts(syst);
  }
  for (unsigned int ippfx = 0; ippfx < N_PPFX; ++ippfx)
  {
    string syst = "ppfx" + to_string(ippfx);
    LoadSyst(syst);
    LoadUncerts(syst);
  }

  c1->SaveAs((outputPDF + "]").c_str());

  MakeSystematicPlots(
    {"muone", "angle", "neutron", "genie", "ppfx", "beam", "stat"}, 
    outputFolder + "/systematicPlots.pdf");
  c1->Clear();

  nominalXSecsENu["cv"]->GetYaxis()->SetRangeUser(0.0, nominalXSecsENu["cv"]->GetMaximum() * 1.5);
  nominalXSecsENu["cv"]->SetLineColor(kRed);
  nominalXSecsENu["cv"]->Draw("hist");
  for (unsigned int igenie = 0; igenie < N_GENIE; ++igenie){
    TH1 * hist = xSecsENu["genie" + to_string(igenie)];
    hist->SetLineColor(kBlue);
    hist->Draw("hist same");
  }
  nominalXSecsENu["cv"]->Draw("hist same");
  c1->SaveAs((outputFolder + "/enuGenie.png").c_str());

  nominalXSecsQ2["cv"]->GetYaxis()->SetRangeUser(0.0, nominalXSecsQ2["cv"]->GetMaximum() * 1.5);
  nominalXSecsQ2["cv"]->SetLineColor(kRed);
  nominalXSecsQ2["cv"]->Draw("hist");
  for (unsigned int ippfx = 0; ippfx < N_PPFX; ++ippfx){
    TH1 * hist = xSecsQ2["ppfx" + to_string(ippfx)];
    hist->SetLineColor(kBlue);
    hist->Draw("hist same");
  }
  nominalXSecsQ2["cv"]->Draw("hist same");
  c1->SaveAs((outputFolder + "/q2Genie.png").c_str());

  outFile->Close();
}

const TAxis q2Axis	(	1	,
		0.	0,
		2.	8
	)

Referenced by FixAxisRangeQ2().

const TAxis uncertAxis	(	1	,
		-0.	2,
		0.	2
	)

TH2D* XSec2D

(

string

syst

)

Definition at line 221 of file PlotXSec.C.

References FixAxisRange2D(), Load3DXSec(), systTitles, and xSecs2D.

Referenced by LoadNominal(), and LoadSyst().

{
  if (xSecs2D.find(syst) != xSecs2D.end())
    return xSecs2D.at(syst);

  TH3D * syst3DXSec = Load3DXSec(syst);
  TH2D * syst2DXSec;
  if (!syst3DXSec)
  {
    xSecs2D[syst] = NULL;
    return NULL;
  } 
  syst2DXSec = (TH2D*) syst3DXSec->Project3D("yx");
  syst2DXSec->Scale(1, "width"); // Area-normalize the xsec
  FixAxisRange2D(syst2DXSec);
  xSecs2D[syst] = syst2DXSec;
  syst2DXSec->SetName(("xsec_" + syst + "_2DProj").c_str());
  syst2DXSec->SetTitle((systTitles[syst] + " XSec").c_str());
  syst2DXSec->Write();
  return syst2DXSec;
}

TH2D* XSec2DUncert

(

string

syst

)

Definition at line 272 of file PlotXSec.C.

References ana::assert(), om::cout, allTimeWatchdog::endl, FixAxisRange2D(), nominalXSecs2D, otherNomSysts, uncerts2D, and xSecs2D.

Referenced by LoadUncerts().

{
  if (uncerts2D.find(syst) != uncerts2D.end())
    return uncerts2D.at(syst);

  assert(xSecs2D.find(syst) != xSecs2D.end());
  TH2D * systXSec = xSecs2D.at(syst);
  if (!systXSec){
    uncerts2D[syst] = NULL;
    return NULL;
  }

  TH2D * thisNom = nominalXSecs2D["cv"];
  if (otherNomSysts.find(syst) != otherNomSysts.end()){
    cout << "Using non-standard nominal : " << otherNomSysts.at(syst) << " for syst " << syst << endl;
    thisNom = nominalXSecs2D[otherNomSysts.at(syst)];
  }
  TH2D * syst2DUncert = new TH2D(*systXSec - *thisNom);
  syst2DUncert->SetName((syst + "_2d_uncert").c_str());
  syst2DUncert->SetTitle((syst + " Uncert").c_str());
  syst2DUncert->Divide(thisNom);
  FixAxisRange2D(syst2DUncert);
  // syst2DUncert->GetZaxis()->SetRangeUser(uncertAxis.GetXmin(), uncertAxis.GetXmax());
  uncerts2D[syst] = syst2DUncert;
  syst2DUncert->Write();
  return syst2DUncert;
}

TH1D* XSecENu

(

string

syst

)

Definition at line 244 of file PlotXSec.C.

References FixAxisRangeENu(), inputFile, systTitles, and xSecsENu.

Referenced by LoadNominal(), and LoadSyst().

{
  if (xSecsENu.find(syst) != xSecsENu.end())
    return xSecsENu.at(syst);

  TH1D * xSecENu = (TH1D*) inputFile->Get((syst + "/ENu/xsec_" + syst + "_ENu").c_str());
  FixAxisRangeENu(xSecENu);
  xSecsENu[syst] = xSecENu;
  xSecENu->SetTitle((systTitles[syst] + " XSec").c_str());
  xSecENu->Write();
  return xSecENu;
}

TH1D* XSecENuUncert

(

string

syst

)

Definition at line 300 of file PlotXSec.C.

References ana::assert(), om::cout, allTimeWatchdog::endl, FixAxisRangeENu(), nominalXSecsENu, otherNomSysts, uncertsENu, and xSecsENu.

Referenced by LoadUncerts().

{
  if (uncertsENu.find(syst) != uncertsENu.end())
    return uncertsENu.at(syst);

  assert(xSecsENu.find(syst) != xSecsENu.end());
  TH1D * systXSec = xSecsENu.at(syst);

  TH1D * thisNom = nominalXSecsENu["cv"];
  if (otherNomSysts.find(syst) != otherNomSysts.end()){
    cout << "Using non-standard nominal : " << otherNomSysts.at(syst) << " for syst " << syst << endl;
    thisNom = nominalXSecsENu[otherNomSysts.at(syst)];
  }
  TH1D * systUncert = new TH1D(*systXSec - *thisNom);
  systUncert->SetName((syst + "_enu_uncert").c_str());
  systUncert->SetTitle((syst + " Uncert").c_str());
  systUncert->Divide(thisNom);
  FixAxisRangeENu(systUncert);
  // systUncert->GetZaxis()->SetRangeUser(uncertAxis.GetXmin(), uncertAxis.GetXmax());
  uncertsENu[syst] = systUncert;
  systUncert->Write();
  return systUncert;
}

TH1D* XSecQ2

(

string

syst

)

Definition at line 258 of file PlotXSec.C.

References FixAxisRangeQ2(), inputFile, systTitles, and xSecsQ2.

Referenced by LoadNominal(), and LoadSyst().

{
  if (xSecsQ2.find(syst) != xSecsQ2.end())
    return xSecsQ2.at(syst);

  TH1D * xSecQ2 = (TH1D*) inputFile->Get((syst + "/Q2/xsec_" + syst + "_Q2").c_str());
  xSecQ2->Scale(1, "width");
  FixAxisRangeQ2(xSecQ2);
  xSecsQ2[syst] = xSecQ2;
  xSecQ2->SetTitle((systTitles[syst] + " XSec").c_str());
  xSecQ2->Write();
  return xSecQ2;
}

TH1D* XSecQ2Uncert

(

string

syst

)

Definition at line 324 of file PlotXSec.C.

References ana::assert(), om::cout, allTimeWatchdog::endl, FixAxisRangeQ2(), nominalXSecsQ2, otherNomSysts, uncertsQ2, and xSecsQ2.

Referenced by LoadUncerts().

{
  if (uncertsQ2.find(syst) != uncertsQ2.end())
    return uncertsQ2.at(syst);

  assert(xSecsQ2.find(syst) != xSecsQ2.end());
  TH1D * systXSec = xSecsQ2.at(syst);

  TH1D * thisNom = nominalXSecsQ2["cv"];
  if (otherNomSysts.find(syst) != otherNomSysts.end()){
    cout << "Using non-standard nominal : " << otherNomSysts.at(syst) << " for syst " << syst << endl;
    thisNom = nominalXSecsQ2[otherNomSysts.at(syst)];
  }
  TH1D * systUncert = new TH1D(*systXSec - *thisNom);
  systUncert->SetName((syst + "_q2_uncert").c_str());
  systUncert->SetTitle((syst + " Uncert").c_str());
  systUncert->Divide(thisNom);
  FixAxisRangeQ2(systUncert);
  // systUncert->GetZaxis()->SetRangeUser(uncertAxis.GetXmin(), uncertAxis.GetXmax());
  uncertsQ2[syst] = systUncert;
  systUncert->Write();
  return systUncert;
}

Variable Documentation

map<string, TH3D*> all3DXSecs

Definition at line 125 of file PlotXSec.C.

Referenced by Load3DXSec().

vector<string> beamSysts

Initial value:

{
  "2kA",
  "02mmBeamSpotSize",
  "1mmBeamShiftX",
  "1mmBeamShiftY",
  "3mmHorn1X",
  "3mmHorn1Y",
  "3mmHorn2X",
  "3mmHorn2Y",
  "7mmTargetZ",
  "MagneticFieldinDecayPipe",
  "1mmHornWater"
}

Definition at line 65 of file PlotXSec.C.

Referenced by GetTotalBeamUncert2D(), GetTotalBeamUncertENu(), GetTotalBeamUncertQ2(), and PlotXSec().

TCanvas* c1

Definition at line 138 of file PlotXSec.C.

Referenced by PlotXSec().

string defaultNominal = "cv"

Definition at line 118 of file PlotXSec.C.

TFile* inputFile

Definition at line 134 of file PlotXSec.C.

Referenced by CalculateXSec(), combineFiles(), FindUniverseSpectra(), Load3DXSec(), make_tmva_NCPi0s_2selection(), PlotXSec(), XSecENu(), and XSecQ2().

map<int, double> maxTMuBins

Initial value:

{
  {2,  1.1}, 
  {3,  1.1}, 
  {4,  1.1}, 
  {5,  1.4}, 
  {6,  1.4}, 
  {7,  1.4}, 
  {8,  1.8}, 
  {9,  1.9}, 
  {10, 2.2}, 
}

Definition at line 94 of file PlotXSec.C.

Referenced by LimitTMuBins().

const unsigned int N_GENIE = 1000

Definition at line 148 of file PlotXSec.C.

Referenced by MakeSystematicPlots(), and PlotXSec().

const unsigned int N_PPFX = 100

Definition at line 149 of file PlotXSec.C.

Referenced by MakeSystematicPlots(), and PlotXSec().

TH2D* nom2DXSec

Definition at line 115 of file PlotXSec.C.

TH3D* nom3DXSec

Definition at line 114 of file PlotXSec.C.

map<string, TH2D*> nominalXSecs2D

Definition at line 120 of file PlotXSec.C.

Referenced by GetTotalBeamUncert2D(), LoadNominal(), LoadStatUncert(), MakeSystematicPlots(), PlotXSec(), and XSec2DUncert().

map<string, TH3D*> nominalXSecs3D

Definition at line 119 of file PlotXSec.C.

Referenced by LoadNominal().

map<string, TH1D*> nominalXSecsENu

Definition at line 121 of file PlotXSec.C.

Referenced by GetTotalBeamUncertENu(), LoadNominal(), LoadStatUncert(), MakeSystematicPlots(), PlotXSec(), and XSecENuUncert().

map<string, TH1D*> nominalXSecsQ2

Definition at line 122 of file PlotXSec.C.

Referenced by GetTotalBeamUncertQ2(), LoadNominal(), LoadStatUncert(), MakeSystematicPlots(), PlotXSec(), and XSecQ2Uncert().

const map<string, string> otherNomSysts

Initial value:

{
  {"neutron",    "neutron_nom"},
  {"light",      "nom_good_seventh"},
  {"calib",      "nom_good_seventh"},
  {"calibshift", "nom_good_seventh"}
}

Definition at line 57 of file PlotXSec.C.

Referenced by PlotXSec(), XSec2DUncert(), XSecENuUncert(), and XSecQ2Uncert().

TFile* outFile

Definition at line 135 of file PlotXSec.C.

Referenced by accum_pot_equivalent_fhc(), NDPredictionHandler::AddVar(), FDPredictionHandler::AddVar(), Analyse_Data2DataComp(), Analyse_Data2DataComp_HigherEnergyCuts(), Analyse_Data2DataComp_kNumu2020ND(), AnalyzeNus18Pred(), AnalyzeNus18Systs(), AnalyzeSysts(), ApplyOscillations(), CalcCutVals(), calculateWrongSignNue(), calculateWrongSignNumuQ1(), calculateWrongSignNumuQ2(), calculateWrongSignNumuQ3(), calculateWrongSignNumuQ4(), CalculateXSec(), combineFiles(), Compare_NoExtrap(), Compare_Spectra(), ComparisonPlots_Data(), ComparisonPlots_MC(), CosmicPred(), demo2p5a(), DoTraining(), DrawExtrapSurface(), Evaluate_BDTMLP_Algs_PredNoExtrap(), Evaluate_BDTMLP_Algs_Spectra(), eventListToTextFile(), example_macro(), FillSpectra(), FillTrainingTrees(), FOMPlot(), GenerateCovMx(), get_numi_data_histogram(), LoadNominal(), LoadSyst(), LoadUncerts(), make_dst_cosrejbdttrain(), make_nue_ana2018_pot(), make_nue_ana2019_epoch7d_pot(), make_nue_ana2019_epoch8b_pot(), make_nue_ana2019_fhc_ub_pot(), make_nue_ana2019_rhc_ub_pot(), make_nue_thirdana_pot(), make_numu_thirdana_pot(), make_nus17_pot(), make_rhcpred_2017(), make_rockpred(), make_rockpred_2017(), MakeCosmicSpectra(), MakeCovMx(), MakeExtrapSurface(), MakeISysts(), MakeNumuCovMx(), MakeNus18ExtrapPred(), MakeNus18SidebandPred(), MakeNus18Systs(), makePIDCutTuning(), makePrediction(), MakePrediction(), MakePredictionNoOsc_FHC_FD(), MakePredictionNoOsc_FHC_ND(), MakePredictionNoOsc_RHC_FD(), MakePredictionNoOsc_RHC_ND(), MakeSurface(), MakeSurfaceBinningStudy(), MakeSurfaceLLTest(), mcTruthPredictions(), meanWeight_macro(), MergePredictions(), MergeSurface(), MergeSurfaces(), MichelDecompTest(), MRDiFStudy_FHC_Step2(), MRDiFStudy_RHC_Step2(), multiverse_efficiency_macro(), NoExtrap(), nuCrystalBall(), nue_michelDataMC(), NueExtrap(), NuMu2020_TrimCAFs(), NumuCosmic(), NumuExtrap(), NuSCalculateCorr(), NuSTreeMaker(), PlotCovariance(), PlotNus18Sideband(), PlotXSec(), RunCalibration(), runTwoSampleDecomp(), SaveSpectra(), SystematicComp(), test_micheldecomp(), ThrowFakeData(), TrimCAFs(), UnfoldXSec(), FDPredictionHandler::~FDPredictionHandler(), and NDPredictionHandler::~NDPredictionHandler().

vector<string> standardSysts

Initial value:

{
  
  "calibshift",
  "ckv"
}

Definition at line 41 of file PlotXSec.C.

Referenced by MakeSystematicPlots(), and PlotXSec().

map<string, pair<string, string> > standardSystsUpDown

Initial value:

{
  {"light",   {"lightup",    "lightdown"}},
  {"calib",   {"calibpos",   "calibneg"}},
  {"muone",   {"muoneup",    "muonedw"}},
  {"angle",   {"angleup",    "angledw"}},
  {"neutron", {"neutron_Up", "neutron_Dw"}}
}

Definition at line 48 of file PlotXSec.C.

Referenced by MakeSystematicPlots(), and PlotXSec().

map<string, Color_t> systColors

Initial value:

{
  {"total", kBlack},
  {"beam", kBlue},
  {"ckv", kOrange},
  {"neutron", kGreen},
  {"angle", kRed},
  {"muone", kViolet},
  {"stat", kGray},
  {"genie", kCyan},
  {"ppfx", kOrange}
}

Definition at line 80 of file PlotXSec.C.

Referenced by MakeSystematicPlots().

map<string, string> systTitles

Initial value:

{
  {"cv",         "Central Value"},
  {"light",      "Light Level"},
  {"calib",      "Calibration"},
  {"calibshift", "Calib Shape"},
  {"ckv",        "Cherenkov"},
  {"muone",      "Muon Energy"},
  {"neutron",    "Neutron Syst"},
  {"angle",      "Muon Angle"},
  {"neutron_nom","Nominal (Neutron syst)"},
  {"nom_good_seventh", "Nominal (calib / light levels)"},
  {"genie",      "XSec and FSI"},
  {"ppfx",       "Beam - HP"},
  {"stat",       "Statistical"},
  {"beam",       "Beam - Focusing"}
}

Definition at line 23 of file PlotXSec.C.

Referenced by Load3DXSec(), MakeSystematicPlots(), XSec2D(), XSecENu(), and XSecQ2().

string uncertLabel = "_uncert"

Definition at line 107 of file PlotXSec.C.

map<string, TH2D*> uncerts2D

Definition at line 129 of file PlotXSec.C.

Referenced by GetTotalBeamUncert2D(), LoadStatUncert(), MakeSystematicPlots(), and XSec2DUncert().

map<string, TH1D*> uncertsENu

Definition at line 130 of file PlotXSec.C.

Referenced by GetTotalBeamUncertENu(), LoadStatUncert(), MakeSystematicPlots(), and XSecENuUncert().

map<string, TH1D*> uncertsQ2

Definition at line 131 of file PlotXSec.C.

Referenced by GetTotalBeamUncertQ2(), LoadStatUncert(), MakeSystematicPlots(), and XSecQ2Uncert().

const string var3DName = "EAvail"

Definition at line 106 of file PlotXSec.C.

Referenced by Load3DXSec().

const vector<string> vars3D

Initial value:

{
  "EAvail"
}

Definition at line 109 of file PlotXSec.C.

map<string, TH2D*> xSecs2D

Definition at line 126 of file PlotXSec.C.

Referenced by MakeSystematicPlots(), XSec2D(), and XSec2DUncert().

map<string, TH1D*> xSecsENu

Definition at line 127 of file PlotXSec.C.

Referenced by MakeSystematicPlots(), PlotXSec(), XSecENu(), and XSecENuUncert().

map<string, TH1D*> xSecsQ2

Definition at line 128 of file PlotXSec.C.

Referenced by MakeSystematicPlots(), PlotXSec(), XSecQ2(), and XSecQ2Uncert().

---

Generated on Sun Jan 24 2021 14:15:27 for NOvA by   1.8.11