12 #include "TDirectory.h" 15 #include "TObjString.h" 54 NDExtrap extrap(nonswap, swap, axis, selection, shiftMC, weight);
149 const Cut& selection,
174 TDirectory*
tmp = gDirectory;
176 dir = dir->mkdir(name.c_str());
179 TObjString(
"NDExtrap").Write(
"type");
185 fEEComp->SaveTo(dir,
"EEComp");
187 fEMComp->SaveTo(dir,
"EMComp");
189 fMMComp->SaveTo(dir,
"MMComp");
191 fMEComp->SaveTo(dir,
"MEComp");
209 dir = dir->GetDirectory(name.c_str());
216 ret->fEEComp = ana::LoadFrom<NDOscillatableSpectrum>(
dir,
"EEComp");
218 ret->fAEAEComp = ana::LoadFrom<NDOscillatableSpectrum>(
dir,
"AEAEComp");
220 ret->fEMComp = ana::LoadFrom<NDOscillatableSpectrum>(
dir,
"EMComp");
222 ret->fAEAMComp = ana::LoadFrom<NDOscillatableSpectrum>(
dir,
"AEAMComp");
224 ret->fMMComp = ana::LoadFrom<NDOscillatableSpectrum>(
dir,
"MMComp");
226 ret->fAMAMComp = ana::LoadFrom<NDOscillatableSpectrum>(
dir,
"AMAMComp");
228 ret->fMEComp = ana::LoadFrom<NDOscillatableSpectrum>(
dir,
"MEComp");
230 ret->fAMAEComp = ana::LoadFrom<NDOscillatableSpectrum>(
dir,
"AMAEComp");
232 ret->fENCComp = ana::LoadFrom<NDOscillatableSpectrum>(
dir,
"ENCComp");
234 ret->fMNCComp = ana::LoadFrom<NDOscillatableSpectrum>(
dir,
"MNCComp");
236 ret->fAENCComp = ana::LoadFrom<NDOscillatableSpectrum>(
dir,
"AENCComp");
238 ret->fAMNCComp = ana::LoadFrom<NDOscillatableSpectrum>(
dir,
"AMNCComp");
240 ret->fNCCompSpec = ana::LoadFrom<Spectrum>(
dir,
"NCComp");
Near Detector underground.
const Cut kIsNumuCC(CCFlavSel(14, 14))
Select CC .
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 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 L/E information, allowing it to be oscillated