#include "TGaxis.h"
#include "TCanvas.h"
#include "TFile.h"
#include "TH1D.h"
#include "TLatex.h"
#include "TLegend.h"
#include "TStyle.h"
#include "TSystem.h"
#include "TF1.h"
#include "TH2.h"
#include "TGraph.h"
#include "CAFAna/Analysis/Exposures.h"
#include "CAFAna/Analysis/Plots.h"
#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Core/SpectrumLoader.h"
#include <iostream>
#include <sstream>
#include <fstream>
#include <vector>
#include <string>
#include <TROOT.h>
Go to the source code of this file.
Classes | |
struct | MySpectra |
Functions | |
TH1D * | SpecToHist (TFile *MyF, std::string SpecName, double POT, std::string AxisTitle) |
void | NormByBinWidth (TH1D *hist) |
void | Print1D (TH1D *hist, std::string plotName) |
void | MakeComparisonPlots (std::vector< TH1D * > TheHists, std::string plotName, std::vector< std::string > Labels, bool var) |
void | MakeComparisonPlotsWithRatio (std::vector< TH1D * > TheHists, std::string plotName, std::vector< std::string > Labels, bool var) |
TH1D * | MakeRatio (TH1D *num, TH1D *denom, int Col, std::string FType, bool var) |
void | FitGauss (TH1D *h) |
void | neutKEsyst_plot () |
Variables | |
const std::vector< std::string > | GENIEStr = { "All" } |
const unsigned int | nGENIE = 1 |
std::vector< std::string > | ND_CutNames = {"ND_2018full"} |
const unsigned int | nND_CutNames = 1 |
unsigned int | FCut_ND = nND_CutNames - 1 |
std::vector< std::string > | FD_CutNames = {"FD_2018full"} |
const unsigned int | nFD_CutNames = 1 |
unsigned int | FCut_FD = nFD_CutNames - 1 |
std::vector< std::string > | QuantNames = { "AllQuant", "Quant1", "Quant2", "Quant3", "Quant4" } |
const unsigned int | nQuantNames = 5 |
std::vector< std::string > | HCNames = {"fhc", "rhc"} |
const unsigned int | nHC = 2 |
std::vector< double > | POT = {kAna2018FHCPOT, kAna2018RHCPOT} |
std::vector< std::string > | DetNames = {"ND", "FD"} |
const unsigned int | nDet = 2 |
void FitGauss | ( | TH1D * | h | ) |
Definition at line 475 of file neutKEsyst_plot.C.
References om::cout, allTimeWatchdog::endl, and extractScale::mean.
Referenced by MakeComparisonPlots(), and MakeComparisonPlotsWithRatio().
void MakeComparisonPlots | ( | std::vector< TH1D * > | TheHists, |
std::string | plotName, | ||
std::vector< std::string > | Labels, | ||
bool | var | ||
) |
Definition at line 257 of file neutKEsyst_plot.C.
References om::cout, allTimeWatchdog::endl, FitGauss(), make_mec_shifts_plots::GetMaximum(), analysePickle::hist, cet::sqlite::max(), SaveAs(), and gen_hdf5record::size.
void MakeComparisonPlotsWithRatio | ( | std::vector< TH1D * > | TheHists, |
std::string | plotName, | ||
std::vector< std::string > | Labels, | ||
bool | var | ||
) |
Definition at line 299 of file neutKEsyst_plot.C.
References cd(), om::cout, Draw(), allTimeWatchdog::endl, FitGauss(), make_mec_shifts_plots::GetMaximum(), analysePickle::hist, make_syst_table_plots::line, MakeRatio(), cet::sqlite::max(), allTimeWatchdog::plotName, SaveAs(), gen_hdf5record::size, string, submit_syst::x, and submit_syst::y.
Referenced by neutKEsyst_plot().
TH1D * MakeRatio | ( | TH1D * | num, |
TH1D * | denom, | ||
int | Col, | ||
std::string | FType, | ||
bool | var | ||
) |
Definition at line 423 of file neutKEsyst_plot.C.
References ana::bins, and art::to_string().
Referenced by MakeComparisonPlotsWithRatio().
void neutKEsyst_plot | ( | ) |
Definition at line 83 of file neutKEsyst_plot.C.
References fillBadChanDBTables::det, FCut_FD, FCut_ND, FD_CutNames, GENIEStr, runNovaSAM::gErrorIgnoreLevel, MySpectra::HadE_nom_FD, MySpectra::HadE_nom_ND, MySpectra::HadE_scale_m_FD, MySpectra::HadE_scale_m_ND, MySpectra::HadE_scale_p_FD, MySpectra::HadE_scale_p_ND, HCNames, MakeComparisonPlotsWithRatio(), ND_CutNames, nFD_CutNames, nGENIE, nHC, nND_CutNames, MySpectra::NNeu_FD, MySpectra::NNeu_ND, NormByBinWidth(), nQuantNames, MySpectra::NuE_nom_FD, MySpectra::NuE_nom_ND, MySpectra::NuE_scale_m_FD, MySpectra::NuE_scale_m_ND, MySpectra::NuE_scale_p_FD, MySpectra::NuE_scale_p_ND, POT, Print1D(), QuantNames, SpecToHist(), string, MySpectra::VisE_FD, MySpectra::VisE_ND, MySpectra::VisE_sum_FD, and MySpectra::VisE_sum_ND.
void NormByBinWidth | ( | TH1D * | hist | ) |
void Print1D | ( | TH1D * | hist, |
std::string | plotName | ||
) |
Definition at line 248 of file neutKEsyst_plot.C.
References SaveAs().
Referenced by neutKEsyst_plot().
TH1D* SpecToHist | ( | TFile * | MyF, |
std::string | SpecName, | ||
double | POT, | ||
std::string | AxisTitle | ||
) |
Definition at line 224 of file neutKEsyst_plot.C.
References ana::CenterTitles(), and ana::Spectrum::LoadFrom().
Referenced by neutKEsyst_plot().
std::vector<std::string> DetNames = {"ND", "FD"} |
Definition at line 48 of file neutKEsyst_plot.C.
unsigned int FCut_FD = nFD_CutNames - 1 |
Definition at line 38 of file neutKEsyst_plot.C.
Referenced by neutKEsyst_plot().
unsigned int FCut_ND = nND_CutNames - 1 |
Definition at line 36 of file neutKEsyst_plot.C.
Referenced by neutKEsyst_plot().
std::vector<std::string> FD_CutNames = {"FD_2018full"} |
Definition at line 37 of file neutKEsyst_plot.C.
Referenced by neutKEsyst_plot().
const std::vector<std::string> GENIEStr = { "All" } |
Definition at line 31 of file neutKEsyst_plot.C.
Referenced by neutKEsyst_plot().
std::vector<std::string> HCNames = {"fhc", "rhc"} |
Definition at line 44 of file neutKEsyst_plot.C.
Referenced by neutKEsyst_plot().
std::vector<std::string> ND_CutNames = {"ND_2018full"} |
Definition at line 35 of file neutKEsyst_plot.C.
Referenced by neutKEsyst_plot().
const unsigned int nDet = 2 |
Definition at line 48 of file neutKEsyst_plot.C.
Referenced by hyperon_macro(), and hyperon_nom_macro().
const unsigned int nFD_CutNames = 1 |
Definition at line 37 of file neutKEsyst_plot.C.
Referenced by neutKEsyst_plot().
const unsigned int nGENIE = 1 |
Definition at line 31 of file neutKEsyst_plot.C.
Referenced by neutKEsyst_plot().
const unsigned int nHC = 2 |
Definition at line 44 of file neutKEsyst_plot.C.
Referenced by neutKEsyst_plot().
const unsigned int nND_CutNames = 1 |
Definition at line 35 of file neutKEsyst_plot.C.
Referenced by neutKEsyst_plot().
const unsigned int nQuantNames = 5 |
Definition at line 41 of file neutKEsyst_plot.C.
Referenced by neutKEsyst_plot().
std::vector<double> POT = {kAna2018FHCPOT, kAna2018RHCPOT} |
Definition at line 45 of file neutKEsyst_plot.C.
Referenced by ana::AddErrorInQuadrature(), AddExptCorrNue(), AddExptCorrNumu(), Analyse_GetEfficiency(), Analyse_GetEfficiency_UseNEntries(), calculateWrongSignNue(), calculateWrongSignNumuQ1(), calculateWrongSignNumuQ2(), calculateWrongSignNumuQ3(), calculateWrongSignNumuQ4(), compare_cos_numi(), datamc(), demoFitContours(), demoFitSlices(), drawLongTerm(), ana::CovMxSurface::DrawSurfacePoint(), drawSystsShiftingNDdata(), drawTimePlots(), drawVsPOT(), EnergyCont_plot(), fake_future_data(), fd_plot(), ana::GenieMultiverseSpectra::FindBandBoundaries(), ana::FluxMultiverseSyst::FindSigmaBoundaries(), get_cosmic_sample(), ana::CutOptimization::GetAllSystDefs(), GetHistsFD(), GetHistsND(), Plotter::GetMCPlots(), ana::Multiverse::GetNSigmaShift(), GetNueData2018(), GetNueData2019(), ana::GetNueData2020(), GetNumuData2019(), ana::GetNumuData2020(), ana::Multiverse::GetNUniverses(), getPOT(), getTimePeakPlots(), getTotPOT(), goodness_of_fit(), hyperon_plot(), joint_fit_2018_contours(), joint_fit_2018_slices(), joint_fit_2019_bestfit(), joint_fit_2019_contours(), joint_fit_2019_slices(), joint_fit_future_bestfit_univ(), joint_fit_future_contour_univ(), ana::TrivialCrossSectionAnalysis::LoadFrom(), make1DsAndRatio(), make_fc_surfaces_2020(), make_mass_and_oct_fc_input_2019(), make_quantiles_histogram_2020(), MakeCutFlow(), makeEnergyEstimator(), makeFlatWeight(), MakeLegend(), MakeNus17CosBkgd(), makeResolutionPlots(), mre_comp_split(), multiverse_reweighting(), neutKEsyst_plot(), neutronE_plot(), ana::GenieMultiverseNormalizedSpectra::NormalizeSpectra(), nuebar_signif(), pi0_xcheck(), plot_datapredictions(), plot_pi0_xcheck(), plot_pid(), plot_predictions(), PandAna.reco_validation.prod5_pid_validation::plot_ratio(), plot_rationoosc(), ana::PlotAllRecoTrue(), ana::PlotAllSignalEstimates(), PlotStack(), Plotting_DataAndPrediction(), Plotting_OverlayStudies(), ana::SpectrumComponents::Purity(), run_joint_fit_2020_bestfit(), run_joint_fit_2020_contours(), run_joint_fit_2020_slices(), selection_story_plots(), ana::ICrossSectionAnalysis::SetPOT(), ana::SpectrumComponents::SpectrumComponents(), syst_table_fit(), ana::Multiverse::ToSpectrum(), ana::FluxMultiverseSyst::ToTH1(), ana::GenieMultiverseSpectra::ToTH1(), ana::GenieMultiverseSpectra::ToTH2(), ana::Multiverse::Transform(), Tutorial2019Fit(), Tutorial2019FitContours(), Tutorial2019FitSlices(), and ana::NDPredictionNoOsc::~NDPredictionNoOsc().
std::vector<std::string> QuantNames = { "AllQuant", "Quant1", "Quant2", "Quant3", "Quant4" } |
Definition at line 41 of file neutKEsyst_plot.C.
Referenced by neutKEsyst_plot().