14 #include "3FlavorAna/Cuts/NueCutsSecondAna.h" 40 const std::string MCfiles =
"defname: prod_caf_R17-03-01-prod3reco.d_nd_genie_nonswap_fhc_nova_v08_full_v1 with stride 4 offset 1";
44 const std::string Datafiles =
"defname: prod_caf_R17-03-01-prod3reco.d_nd_genie_nonswap_fhc_nova_v08_full_v1 with stride 4 offset 0";
54 Spectrum data_reweighted (lNDData, sideband_x_axis, sideband_y_axis,
55 kCVNCutSideband,
kNoShift, kReweight);
56 Spectrum data_tufts (lNDData, sideband_x_axis, sideband_y_axis,
57 kCVNCutSideband,
kNoShift, kTuftsWeightCC);
62 Spectrum data_sig_sb(lNDData, sideband_x_axis, sideband_y_axis,
65 Spectrum data_bkg_numu_sb(lNDData, sideband_x_axis, sideband_y_axis,
66 kCVNCutSideband && kIsNumuBkg,
68 Spectrum data_bkg_nc_sb(lNDData, sideband_x_axis, sideband_y_axis,
69 kCVNCutSideband &&
kIsNC,
71 Spectrum data_sig_anti_sb(lNDData, sideband_x_axis, sideband_y_axis,
74 Spectrum data_bkg_other_sb(lNDData, sideband_x_axis, sideband_y_axis,
75 kCVNCutSideband && kIsOther_nueccinc,
78 Spectrum data_sig_sb_weight(lNDData, sideband_x_axis, sideband_y_axis,
82 Spectrum data_bkg_numu_sb_weight(lNDData, sideband_x_axis, sideband_y_axis,
83 kCVNCutSideband && kIsNumuBkg,
85 Spectrum data_bkg_nc_sb_weight(lNDData, sideband_x_axis, sideband_y_axis,
86 kCVNCutSideband &&
kIsNC,
88 Spectrum data_sig_anti_sb_weight(lNDData, sideband_x_axis, sideband_y_axis,
91 Spectrum data_bkg_other_sb_weight(lNDData, sideband_x_axis, sideband_y_axis,
92 kCVNCutSideband && kIsOther_nueccinc,
95 Spectrum mc_sig_sb(lNDMC, sideband_x_axis, sideband_y_axis,
98 Spectrum mc_bkg_numu_sb(lNDMC, sideband_x_axis, sideband_y_axis,
99 kCVNCutSideband && kIsNumuBkg,
101 Spectrum mc_bkg_nc_sb(lNDMC, sideband_x_axis, sideband_y_axis,
102 kCVNCutSideband &&
kIsNC,
104 Spectrum mc_sig_anti_sb(lNDMC, sideband_x_axis, sideband_y_axis,
107 Spectrum mc_bkg_other_sb(lNDMC, sideband_x_axis, sideband_y_axis,
108 kCVNCutSideband && kIsOther_nueccinc,
116 TFile*
out =
new TFile(
"SidebandFittest.root",
"recreate");
119 data_reweighted.
SaveTo(out->mkdir(
"data_reweighted"));
120 data_tufts.
SaveTo(out->mkdir(
"data_tufts"));
121 data_sig_sb.
SaveTo(out->mkdir(
"data_sig_sb"));
122 data_bkg_numu_sb.
SaveTo(out->mkdir(
"data_bkg_numu_sb"));
123 data_bkg_nc_sb.
SaveTo(out->mkdir(
"data_bkg_nc_sb"));
124 data_sig_anti_sb.
SaveTo(out->mkdir(
"data_sig_anti_sb"));
125 data_bkg_other_sb.
SaveTo(out->mkdir(
"data_bkg_other_sb"));
126 mc_sig_sb.
SaveTo(out->mkdir(
"mc_sig_sb"));
127 mc_bkg_numu_sb.
SaveTo(out->mkdir(
"mc_bkg_numu_sb"));
128 mc_bkg_nc_sb.
SaveTo(out->mkdir(
"mc_bkg_nc_sb"));
129 mc_sig_anti_sb.
SaveTo(out->mkdir(
"mc_sig_anti_sb"));
130 mc_bkg_other_sb.
SaveTo(out->mkdir(
"mc_bkg_other_sb"));
131 data_sig_sb_weight.
SaveTo(out->mkdir(
"data_sig_sb_weight"));
132 data_bkg_numu_sb_weight.
SaveTo(out->mkdir(
"data_bkg_numu_sb_weight"));
133 data_bkg_nc_sb_weight.
SaveTo(out->mkdir(
"data_bkg_nc_sb_weight"));
134 data_sig_anti_sb_weight.
SaveTo(out->mkdir(
"data_sig_anti_sb_weight"));
135 data_bkg_other_sb_weight.
SaveTo(out->mkdir(
"data_bkg_other_sb_weight"));
Cuts and Vars for the 2020 FD DiF Study.
void SetSpillCut(const SpillCut &cut)
Representation of a spectrum in any variable, with associated POT.
const Cut kIsNC([](const caf::SRProxy *sr){if(sr->mc.nnu==0) return false;assert(sr->mc.nnu==1);return!sr->mc.nu[0].iscc;})
Is this a Neutral Current event?
virtual void Go() override
Load all the registered spectra.
void SaveTo(TDirectory *dir, const std::string &name) const
const SystShifts kNoShift
Collaborates with Spectrum and OscillatableSpectrum to fill spectra from CAF files.
const SpillCut kStandardSpillCuts
Apply this unless you're doing something special.