#include <stdio.h>
#include <stdlib.h>
#include <string>
#include "Utilities/rootlogon.C"
#include "CAFAna/Cuts/SpillCuts.h"
#include "CAFAna/Core/Utilities.h"
#include "CAFAna/Core/SpectrumLoader.h"
#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Cuts/TimingCuts.h"
#include "CAFAna/Analysis/Plots.h"
#include "CAFAna/Analysis/Style.h"
#include "CAFAna/Core/LoadFromFile.h"
#include "CAFAna/Core/Ratio.h"
#include "CAFAna/XSec/GenieMultiverseSyst.h"
#include "TH1D.h"
#include "TCanvas.h"
#include "TGraph.h"
#include "TMultiGraph.h"
#include "TLine.h"
#include "TLegend.h"
#include "TPad.h"
#include "TGaxis.h"
Go to the source code of this file.
Classes | |
struct | FileUp |
struct | DateUp |
Functions | |
void | NDSimulation () |
void | PlotWithSystErrorBandWidth (TH1 *&nom, std::vector< TH1 * > &ups, std::vector< TH1 * > &dns, int col, int errCol, float headroom, bool newaxis) |
void | PlotDeltaSigmaWithSigma (TH1 *&nom, TH1 *&nomsel, std::vector< TH1 * > &ups, std::vector< TH1 * > &dns, std::vector< TH1 * > &upssel, std::vector< TH1 * > &dnssel, TH1 *&heff, std::vector< TH1 * > &effups, std::vector< TH1 * > &effdns, TString outname1, TString outname2, TString outname3, TString outname4, TString outname5, TString outname6, TString outname7, TString outname8, TString outname9, Bool_t printme=false) |
void | GetFluxError (const Spectrum *nom, std::vector< ana::Spectrum * > univs, std::vector< ana::Spectrum * > &supdown) |
void | plot_uncertainty () |
Variables | |
const double | pot = 8.09e20 |
TString | outDir2 = "/nova/app/users/acedeno/tag_releasesS18-06-14/NDAna/ncpi0_semi_inc_png_cvn/NewDataSetSyst" |
const int | File =13 |
const FileUp | CurrentFile [File] |
const int | Date = 1 |
const DateUp | CurrentDate [Date] |
void GetFluxError | ( | const Spectrum * | nom, |
std::vector< ana::Spectrum * > | univs, | ||
std::vector< ana::Spectrum * > & | supdown | ||
) |
Definition at line 428 of file plot_uncertainty.C.
References make_syst_table_plots::ibin, ana::Spectrum::Livetime(), extractScale::mean, pot, cet::pow(), std::sqrt(), and ana::Spectrum::ToTH1().
Referenced by plot_uncertainty().
void NDSimulation | ( | ) |
Definition at line 41 of file plot_uncertainty.C.
References prelim.
Referenced by plot_uncertainty(), and PlotDeltaSigmaWithSigma().
void plot_uncertainty | ( | ) |
***********************************************************************************************************************************************************************************************////
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
***********************************************************************************************************************************************************************************************////
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
***********************************************************************************************************************************************************************************************////
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
***********************************************************************************************************************************************************************************************////
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
TH1* truesignal_ckv = ckv[2]->ToTH1(pot); ///comment this if just caliip or down is needed for plots
***********************************************************************************************************************************************************************************************////
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
*************************************************************************************************************************************************************************************///// /Total PreSel : kVtxInFiducial && kContainment && k2Prongs && kPi0RemIDND && kNpCVN1gammaID;
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment lightup[3]->ToTH1(pot), ckv[3]->ToTH1(pot), calib_shape[3]->ToTH1(pot), fluxup[3]->ToTH1(pot), xsecup[3]->ToTH1(pot) if plotting just calib.
comment lightup[3]->ToTH1(pot), ckv[3]->ToTH1(pot), calib_shape[3]->ToTH1(pot), fluxup[3]->ToTH1(pot), xsecup[3]->ToTH1(pot) if plotting just calib.
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
comment this if just caliip or down is needed for plots
Definition at line 457 of file plot_uncertainty.C.
References om::cout, cut_labels, DateUp::date, allTimeWatchdog::endl, MakeMiniprodValidationCuts::f, fin, plot_validation_datamc::fname, GetFluxError(), ana::GenieMultiverseSpectra::kBandFromNominal, kRed, ana::Spectrum::LoadFrom(), ana::GenieMultiverseSpectra::LowerSigma(), getGoodRuns4SAM::n, NDSimulation(), outDir2, PlotDeltaSigmaWithSigma(), PlotWithSystErrorBandWidth(), pot, PandAna.Demos.pi0_spectra::pur, runNovaSAM::release, std::sqrt(), string, art::to_string(), and ana::GenieMultiverseSpectra::UpperSigma().
void PlotDeltaSigmaWithSigma | ( | TH1 *& | nom, |
TH1 *& | nomsel, | ||
std::vector< TH1 * > & | ups, | ||
std::vector< TH1 * > & | dns, | ||
std::vector< TH1 * > & | upssel, | ||
std::vector< TH1 * > & | dnssel, | ||
TH1 *& | heff, | ||
std::vector< TH1 * > & | effups, | ||
std::vector< TH1 * > & | effdns, | ||
TString | outname1, | ||
TString | outname2, | ||
TString | outname3, | ||
TString | outname4, | ||
TString | outname5, | ||
TString | outname6, | ||
TString | outname7, | ||
TString | outname8, | ||
TString | outname9, | ||
Bool_t | printme = false |
||
) |
uncertainty on Selected events
sigma
bkg syst
bkg Nsel
Numerator bkg Nsel
Denominator bkg Nsel
Sig Nsel
Definition at line 140 of file plot_uncertainty.C.
References demo5::c1, demo5::c2, chisquared::c3, chisquared::c4, om::cout, allTimeWatchdog::endl, kRed, NDSimulation(), outDir, std::sqrt(), and submit_syst::y.
Referenced by plot_uncertainty().
void PlotWithSystErrorBandWidth | ( | TH1 *& | nom, |
std::vector< TH1 * > & | ups, | ||
std::vector< TH1 * > & | dns, | ||
int | col, | ||
int | errCol, | ||
float | headroom, | ||
bool | newaxis | ||
) |
Definition at line 82 of file plot_uncertainty.C.
References om::cout, allTimeWatchdog::endl, MECModelEnuComparisons::g, hi(), ana::kTotalMCColor, ana::kTotalMCErrorBandColor, lo(), std::sqrt(), std::swap(), w, and submit_syst::y.
Referenced by plot_uncertainty().
Definition at line 77 of file plot_uncertainty.C.
Referenced by bdtstudyplotter(), plotHist_SpectrumCVNID(), resolutionplotter(), and twodvtxcontplotter().
Definition at line 62 of file plot_uncertainty.C.
Referenced by bdtstudyplotter(), plotHist_SpectrumCVNID(), resolutionplotter(), and twodvtxcontplotter().
const int Date = 1 |
Definition at line 76 of file plot_uncertainty.C.
Referenced by bdtstudyplotter(), plotHist_SpectrumCVNID(), resolutionplotter(), and twodvtxcontplotter().
const int File =13 |
Definition at line 61 of file plot_uncertainty.C.
Referenced by bdtstudyplotter(), FEBLoop(), hdf5::HDF5Maker::InitializeOutfile(), plotHist_SpectrumCVNID(), resolutionplotter(), and twodvtxcontplotter().
TString outDir2 = "/nova/app/users/acedeno/tag_releasesS18-06-14/NDAna/ncpi0_semi_inc_png_cvn/NewDataSetSyst" |
Definition at line 53 of file plot_uncertainty.C.
Referenced by bdtstudyplotter(), plot_uncertainty(), plotHist_SpectrumCVNID(), resolutionplotter(), and twodvtxcontplotter().
const double pot = 8.09e20 |
Definition at line 37 of file plot_uncertainty.C.
Referenced by GetFluxError(), and plot_uncertainty().