13 #include "Utilities/func/MathUtil.h" 36 LoadFromFile<PredictionExtrap>(dirName+fName1,
"predProp").
release();
39 LoadFromFile<PredictionNoExtrap>(dirName+fName1,
"fdmcNom").
release();
45 LoadFromFile<PredictionExtrap>(dirName+fName1,
"predProp").
release();
48 LoadFromFile<PredictionNoExtrap>(dirName+fName1,
"fdmcNom").
release();
69 TH1 *htot = tot.
ToTH1(pot);
71 TH1 *hbb = bb. ToTH1(pot);
73 htot->SetLineColor(
kRed);
74 hsig->SetLineColor(kPink+9);
75 hbb ->SetLineColor(kAzure);
79 TCanvas *
can =
new TCanvas();
81 hsig->Draw(
"hist,same");
82 hbb ->Draw(
"hist,same");
89 bb_merge -= sig_merge;
91 TH1 *htot_merge = tot_merge.
ToTH1(pot);
92 TH1 *hsig_merge = sig_merge.
ToTH1(pot);
93 TH1 *hbb_merge = bb_merge. ToTH1(pot);
95 htot_merge->SetLineColor(
kRed);
96 hsig_merge->SetLineColor(kPink+9);
97 hbb_merge ->SetLineColor(kAzure);
99 TCanvas *can_merge =
new TCanvas();
100 htot_merge->Draw(
"hist");
101 hsig_merge->Draw(
"hist,same");
102 hbb_merge ->Draw(
"hist,same");
105 std::cout <<
"Total: Merged - " << htot_merge->Integral(28,28) <<
" and Not Merged - " << htot->Integral(28,36) <<
std::endl;
106 std::cout <<
"Signal: Merged - " << hsig_merge->Integral(28,28) <<
" and Not Merged - " << hsig->Integral(28,36) <<
std::endl;
107 std::cout <<
"Beam Bkgd: Merged - " << hbb_merge->Integral(28,28) <<
" and Not Merged - " << hbb->Integral(28,36) <<
std::endl;
Cuts and Vars for the 2020 FD DiF Study.
TH1D * ToTH1(double exposure, Color_t col=kBlack, Style_t style=kSolid, EExposureType expotype=kPOT, EBinType bintype=kBinContent) const
Histogram made from this Spectrum, scaled to some exposure.
const FitDmSq32 kFitDmSq32
void SetValue(osc::IOscCalcAdjustable *osc, double val) const override
void SetValue(osc::IOscCalcAdjustable *osc, double val) const override
Forward to wrapped Var's SetValue()
osc::IOscCalcAdjustable * DefaultOscCalc()
Create a new calculator with default assumptions for all parameters.
Representation of a spectrum in any variable, with associated POT.
Charged-current interactions.
void test_nue2017Prediction()
void SetValue(osc::IOscCalcAdjustable *osc, double val) const override
virtual Spectrum Predict(osc::IOscCalc *calc) const override
const ConstrainedFitVarWithPrior fitSsqTh23_UniformPriorSsqTh23 & kFitSinSqTheta23
virtual Spectrum PredictComponent(osc::IOscCalc *calc, Flavors::Flavors_t flav, Current::Current_t curr, Sign::Sign_t sign) const override
const FitVarWithPrior fitDeltaInPiUnits_UniformPriordCP & kFitDeltaInPiUnits