12 #include "TDirectory.h" 15 #include "TObjString.h" 54 FDExtrap extrap(nonswap, swap, tau, axis, selection, shiftMC, weight);
166 const Cut& selection,
198 TDirectory*
tmp = gDirectory;
200 dir = dir->mkdir(name.c_str());
203 TObjString(
"FDExtrap").Write(
"type");
206 fEEComp->SaveTo(dir,
"EEComp");
210 fEMComp->SaveTo(dir,
"EMComp");
214 fETComp->SaveTo(dir,
"ETComp");
218 fMEComp->SaveTo(dir,
"MEComp");
222 fMMComp->SaveTo(dir,
"MMComp");
226 fMTComp->SaveTo(dir,
"MTComp");
247 dir = dir->GetDirectory(name.c_str());
254 ret->fEEComp = ana::LoadFrom<OscillatableSpectrum>(
dir,
"EEComp");
256 ret->fAEAEComp = ana::LoadFrom<OscillatableSpectrum>(
dir,
"AEAEComp");
260 ret->fEMComp = ana::LoadFrom<OscillatableSpectrum>(
dir,
"EMComp");
262 ret->fAEAMComp = ana::LoadFrom<OscillatableSpectrum>(
dir,
"AEAMComp");
266 ret->fETComp = ana::LoadFrom<OscillatableSpectrum>(
dir,
"ETComp");
268 ret->fAEATComp = ana::LoadFrom<OscillatableSpectrum>(
dir,
"AEATComp");
272 ret->fMEComp = ana::LoadFrom<OscillatableSpectrum>(
dir,
"MEComp");
274 ret->fAMAEComp = ana::LoadFrom<OscillatableSpectrum>(
dir,
"AMAEComp");
278 ret->fMMComp = ana::LoadFrom<OscillatableSpectrum>(
dir,
"MMComp");
280 ret->fAMAMComp = ana::LoadFrom<OscillatableSpectrum>(
dir,
"AMAMComp");
284 ret->fMTComp = ana::LoadFrom<OscillatableSpectrum>(
dir,
"MTComp");
286 ret->fAMATComp = ana::LoadFrom<OscillatableSpectrum>(
dir,
"AMATComp");
290 ret->fENCComp = ana::LoadFrom<OscillatableSpectrum>(
dir,
"ENCComp");
292 ret->fMNCComp = ana::LoadFrom<OscillatableSpectrum>(
dir,
"MNCComp");
294 ret->fTNCComp = ana::LoadFrom<OscillatableSpectrum>(
dir,
"TNCComp");
296 ret->fATNCComp = ana::LoadFrom<OscillatableSpectrum>(
dir,
"ATNCComp");
298 ret->fAENCComp = ana::LoadFrom<OscillatableSpectrum>(
dir,
"AENCComp");
300 ret->fAMNCComp = ana::LoadFrom<OscillatableSpectrum>(
dir,
"AMNCComp");
302 ret->fNCCompSpec = ana::LoadFrom<Spectrum>(
dir,
"NCComp");
const Cut kIsNumuCC(CCFlavSel(14, 14))
Select CC .
Far Detector at Ash River.
Cuts and Vars for the 2020 FD DiF Study.
Simple record of shifts applied to systematic parameters.
const Cut kIsNumu([](const caf::SRProxy *sr){return(sr->mc.nnu==1 &&abs(sr->mc.nu[0].pdg)==14);})
Collection of SpectrumLoaders for many configurations.
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?
SpectrumLoaderBase & GetLoader(caf::Det_t det, DataMC datamc, DataSource src=kBeam, SwappingConfig swap=kNonSwap)
Retrieve a specific loader.
const Cut kIsSig(CCFlavSel(12, 14))
Select CC .
const Cut kIsNue([](const caf::SRProxy *sr){return(sr->mc.nnu==1 &&abs(sr->mc.nu[0].pdg)==12);})
Base class for the various types of spectrum loader.
assert(nhit_max >=nhit_nbins)
Spectrum with true energy information, allowing it to be oscillated