8 #include "TObjString.h" 9 #include "TDirectory.h" 19 TDirectory*
tmp = gDirectory;
21 dir = dir->mkdir(name.c_str());
24 TObjString(
"ModularExtrap").Write(
"type");
54 TDirectory*
tmp = gDirectory;
56 fMEextrap->SavePlots( dir->mkdir(
"MEextrap"), potFD );
57 fMEAntiextrap->SavePlots( dir->mkdir(
"MEAntiextrap"), potFD );
58 fEEextrap->SavePlots( dir->mkdir(
"EEextrap"), potFD );
59 fMMextrap->SavePlots( dir->mkdir(
"MMextrap"), potFD );
68 TDirectory*
tmp = gDirectory;
70 fMEextrap->SavePlots( dir->mkdir(
"MEextrap"), potFD );
71 fMEAntiextrap->SavePlots( dir->mkdir(
"MEAntiextrap"), potFD );
72 fEEextrap->SavePlots( dir->mkdir(
"EEextrap"), potFD );
73 fEEAntiextrap->SavePlots( dir->mkdir(
"EEAntiextrap"), potFD );
74 fMMextrap->SavePlots( dir->mkdir(
"MMextrap"), potFD );
75 fMMAntiextrap->SavePlots( dir->mkdir(
"MMAntiextrap"), potFD );
76 fNCextrap->SavePlots( dir->mkdir(
"NCextrap"), potFD );
77 fNCAntiextrap->SavePlots( dir->mkdir(
"NCAntiextrap"), potFD );
85 TDirectory*
tmp = gDirectory;
87 fMMextrap->SavePlots( dir->mkdir(
"MMextrap"), potFD );
88 fMMAntiextrap->SavePlots( dir->mkdir(
"MMAntiextrap"), potFD );
96 dir = dir->GetDirectory(name.c_str());
160 farMCnonswap, axis, fdcut, shiftMC, weight,
kIsNC, farMCswap, farMCtauswap)),
162 farMCnonswap, axis, fdcut, shiftMC, weight,
kIsNC&&!
kIsAntiNu, farMCswap, farMCtauswap)),
164 farMCnonswap, axis, fdcut, shiftMC, weight,
kIsNC&&
kIsAntiNu, farMCswap, farMCtauswap)),
207 {
return fNCextrap->Return().Unoscillated();}
const Cut kIsNumuCC(CCFlavSel(14, 14))
Select CC .
std::unique_ptr< ModularExtrapComponent > LoadFrom< ModularExtrapComponent >(TDirectory *dir, const std::string &label)
Cuts and Vars for the 2020 FD DiF Study.
"Extrapolates" component by returning FD Monte Carlo.
Simple record of shifts applied to systematic parameters.
const Cut kIsAntiNu([](const caf::SRProxy *sr){if(sr->mc.nnu==0) return false;assert(sr->mc.nnu==1);return sr->mc.nu[0].pdg< 0;})
Is this truly an antineutrino?
const Cut kIsBeamNue(CCFlavSel(12, 12))
Select CC .
const Cut kIsNumuApp(CCFlavSel(14, 12))
Select CC .
Representation of a spectrum in any variable, with associated POT.
const Cut kIsTauFromE(CCFlavSel(16, 12))
Select CC .
const Cut kIsTauFromMu(CCFlavSel(16, 14))
Select CC .
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?
const Cut kIsSig(CCFlavSel(12, 14))
Select CC .
Base class for the various types of spectrum loader.
REGISTER_LOADFROM("BENDecomp", IDecomp, BENDecomp)
assert(nhit_max >=nhit_nbins)
Spectrum with true energy information, allowing it to be oscillated