Public Member Functions | Static Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
ana::PredictionNoOsc Class Reference

Prediction that wraps a simple Spectrum. More...

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-09-26/CAFAna/Prediction/PredictionNoOsc.h"

Inheritance diagram for ana::PredictionNoOsc:
ana::IPrediction

Public Member Functions

 PredictionNoOsc (SpectrumLoaderBase &loader, const HistAxis &axis, const Cut &cut, const SystShifts &shift=kNoShift, const Var &wei=kUnweighted)
 
 PredictionNoOsc (SpectrumLoaderBase &loader, const std::string &label, const Binning &bins, const Var &var, const Cut &cut, const SystShifts &shift=kNoShift, const Var &wei=kUnweighted)
 
virtual void SaveTo (TDirectory *dir, const std::string &name) const override
 
virtual Spectrum Predict (osc::IOscCalc *) const override
 
virtual Spectrum PredictComponent (osc::IOscCalc *calc, Flavors::Flavors_t flav, Current::Current_t curr, Sign::Sign_t sign) const override
 
virtual Spectrum PredictUnoscillated () const
 
virtual Spectrum Predict (osc::IOscCalcStan *calc) const
 
virtual Spectrum PredictSyst (osc::IOscCalc *calc, const SystShifts &syst) const
 
virtual Spectrum PredictSyst (osc::IOscCalcStan *calc, const SystShifts &syst) const
 
virtual Spectrum PredictComponent (osc::IOscCalcStan *calc, Flavors::Flavors_t flav, Current::Current_t curr, Sign::Sign_t sign) const
 
virtual Spectrum PredictComponentSyst (osc::IOscCalc *calc, const SystShifts &syst, Flavors::Flavors_t flav, Current::Current_t curr, Sign::Sign_t sign) const
 
virtual Spectrum PredictComponentSyst (osc::IOscCalcStan *calc, const SystShifts &syst, Flavors::Flavors_t flav, Current::Current_t curr, Sign::Sign_t sign) const
 
virtual OscillatableSpectrum ComponentCC (int from, int to) const
 
virtual Spectrum ComponentNCTotal () const
 
virtual Spectrum ComponentNC () const
 
virtual Spectrum ComponentNCAnti () const
 

Static Public Member Functions

static std::unique_ptr< PredictionNoOscLoadFrom (TDirectory *dir, const std::string &name)
 

Protected Member Functions

 PredictionNoOsc (const Spectrum &s, const Spectrum &sNC, const Spectrum &sNumu, const Spectrum &sNumubar, const Spectrum &sNue, const Spectrum &sNuebar)
 

Protected Attributes

Spectrum fSpectrum
 
Spectrum fSpectrumNC
 
Spectrum fSpectrumNumu
 
Spectrum fSpectrumNumubar
 
Spectrum fSpectrumNue
 
Spectrum fSpectrumNuebar
 

Detailed Description

Prediction that wraps a simple Spectrum.

Definition at line 11 of file PredictionNoOsc.h.

Constructor & Destructor Documentation

ana::PredictionNoOsc::PredictionNoOsc ( SpectrumLoaderBase loader,
const HistAxis axis,
const Cut cut,
const SystShifts shift = kNoShift,
const Var wei = kUnweighted 
)

Definition at line 30 of file PredictionNoOsc.cxx.

Referenced by ana::NoOscPredictionGenerator::Generate(), and LoadFrom().

35  : fSpectrum( loader, axis, cut, shift, wei),
36  fSpectrumNC( loader, axis, cut && kIsNC, shift, wei),
37  fSpectrumNumu( loader, axis, cut && !kIsNC && kIsNumuCC && !kIsAntiNu, shift, wei),
38  fSpectrumNumubar(loader, axis, cut && !kIsNC && kIsNumuCC && kIsAntiNu, shift, wei),
39  fSpectrumNue( loader, axis, cut && !kIsNC && kIsBeamNue && !kIsAntiNu, shift, wei),
40  fSpectrumNuebar( loader, axis, cut && !kIsNC && kIsBeamNue && kIsAntiNu, shift, wei)
41  {
42  }
const Cut kIsNumuCC(CCFlavSel(14, 14))
Select CC .
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?
Definition: TruthCuts.h:53
const Cut kIsBeamNue(CCFlavSel(12, 12))
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?
Definition: TruthCuts.h:8
loader
Definition: demo0.py:10
const Cut cut
Definition: exporter_fd.C:30
ana::PredictionNoOsc::PredictionNoOsc ( SpectrumLoaderBase loader,
const std::string &  label,
const Binning bins,
const Var var,
const Cut cut,
const SystShifts shift = kNoShift,
const Var wei = kUnweighted 
)

Definition at line 18 of file PredictionNoOsc.cxx.

25  : PredictionNoOsc(loader, HistAxis(label, bins, var), cut, shift, wei)
26  {
27  }
const char * label
PredictionNoOsc(SpectrumLoaderBase &loader, const HistAxis &axis, const Cut &cut, const SystShifts &shift=kNoShift, const Var &wei=kUnweighted)
loader
Definition: demo0.py:10
const Binning bins
Definition: NumuCC_CPiBin.h:8
const Cut cut
Definition: exporter_fd.C:30
GenericHistAxis< Var > HistAxis
Definition: HistAxis.h:111
ana::PredictionNoOsc::PredictionNoOsc ( const Spectrum s,
const Spectrum sNC,
const Spectrum sNumu,
const Spectrum sNumubar,
const Spectrum sNue,
const Spectrum sNuebar 
)
inlineprotected

Definition at line 48 of file PredictionNoOsc.h.

52  : fSpectrum(s),
53  fSpectrumNC(sNC),
54  fSpectrumNumu(sNumu), fSpectrumNumubar(sNumubar),
55  fSpectrumNue(sNue), fSpectrumNuebar(sNuebar)
56  {
57  }
const XML_Char * s
Definition: expat.h:262

Member Function Documentation

OscillatableSpectrum ana::IPrediction::ComponentCC ( int  from,
int  to 
) const
virtualinherited

Reimplemented in ana::PredictionExtendToPeripheral, ana::PredictionAddRock, ana::PredictionCombinePeriods, ana::PredictionExtrapSum, and ana::PredictionExtrap.

Definition at line 128 of file IPrediction.cxx.

References om::cout, and allTimeWatchdog::endl.

Referenced by ana::PredictionCombinePeriods::ComponentCC(), ana::PredictionAddRock::ComponentCC(), and ND_DataMC().

129  {
130  std::cout << "WARNING! ComponentCC is unimplemented in IPrediction" << std::endl; abort();
131  }
OStream cout
Definition: OStream.cxx:6
Spectrum ana::IPrediction::ComponentNC ( ) const
virtualinherited

Reimplemented in ana::PredictionExtendToPeripheral, ana::PredictionAddRock, ana::PredictionCombinePeriods, ana::PredictionExtrapSum, and ana::PredictionExtrap.

Definition at line 138 of file IPrediction.cxx.

References om::cout, and allTimeWatchdog::endl.

Referenced by ana::PredictionAddRock::ComponentNC().

139  {
140  std::cout << "WARNING! ComponentNC is unimplemented in IPrediction" << std::endl; abort();
141  }
OStream cout
Definition: OStream.cxx:6
Spectrum ana::IPrediction::ComponentNCAnti ( ) const
virtualinherited

Reimplemented in ana::PredictionExtendToPeripheral, ana::PredictionAddRock, ana::PredictionCombinePeriods, ana::PredictionExtrapSum, and ana::PredictionExtrap.

Definition at line 143 of file IPrediction.cxx.

References om::cout, and allTimeWatchdog::endl.

Referenced by ana::PredictionAddRock::ComponentNCAnti().

144  {
145  std::cout << "WARNING! ComponentNCAnti is unimplemented in IPrediction" << std::endl; abort();
146  }
OStream cout
Definition: OStream.cxx:6
Spectrum ana::IPrediction::ComponentNCTotal ( ) const
virtualinherited

Reimplemented in ana::PredictionExtendToPeripheral, ana::PredictionAddRock, ana::PredictionCombinePeriods, ana::PredictionExtrapSum, and ana::PredictionExtrap.

Definition at line 133 of file IPrediction.cxx.

References om::cout, and allTimeWatchdog::endl.

Referenced by ana::PredictionAddRock::ComponentNCTotal(), and ND_DataMC().

134  {
135  std::cout << "WARNING! ComponentNCTotal is unimplemented in IPrediction" << std::endl; abort();
136  }
OStream cout
Definition: OStream.cxx:6
std::unique_ptr< PredictionNoOsc > ana::PredictionNoOsc::LoadFrom ( TDirectory *  dir,
const std::string &  name 
)
static

Definition at line 103 of file PredictionNoOsc.cxx.

References ana::assert(), dir, PredictionNoOsc(), and runNovaSAM::ret.

Referenced by calculateComponentsNumu(), calculateWrongSignNue(), calculateWrongSignNumuQ1(), calculateWrongSignNumuQ2(), calculateWrongSignNumuQ3(), calculateWrongSignNumuQ4(), and ana::LoadFrom< IPrediction >().

104  {
105  dir = dir->GetDirectory(name.c_str()); // switch to subdir
106  assert(dir);
107 
109  *ana::LoadFrom<Spectrum>(dir, "spect"),
110  *ana::LoadFrom<Spectrum>(dir, "spect_nc"),
111  *ana::LoadFrom<Spectrum>(dir, "spect_numu"),
112  *ana::LoadFrom<Spectrum>(dir, "spect_numubar"),
113  *ana::LoadFrom<Spectrum>(dir, "spect_nue"),
114  *ana::LoadFrom<Spectrum>(dir, "spect_nuebar"));
115 
116  delete dir;
117 
118  // Can't use make_unique because constructor is protected
119  return std::unique_ptr<PredictionNoOsc>(ret);
120  }
const XML_Char * name
Definition: expat.h:151
PredictionNoOsc(SpectrumLoaderBase &loader, const HistAxis &axis, const Cut &cut, const SystShifts &shift=kNoShift, const Var &wei=kUnweighted)
TDirectory * dir
Definition: macro.C:5
assert(nhit_max >=nhit_nbins)
virtual Spectrum ana::PredictionNoOsc::Predict ( osc::IOscCalc ) const
inlineoverridevirtual

Implements ana::IPrediction.

Definition at line 37 of file PredictionNoOsc.h.

References calc, febshutoff_auto::curr, fSpectrum, PredictComponent(), and canMan::sign().

Referenced by PredictComponent().

38  {
39  return fSpectrum;
40  }
Spectrum ana::IPrediction::Predict ( osc::IOscCalcStan calc) const
virtualinherited

Reimplemented in ana::PredictionInterp, ana::PredictionExtendToPeripheral, ana::PredictionExtrapSum, ana::PredictionExtrap, and ana::PredictionAddRock.

Definition at line 66 of file IPrediction.cxx.

References ana::assert().

67  {
68  assert(0 && "This Prediction hasn't implemented a Stan-aware Predict()!");
69  }
assert(nhit_max >=nhit_nbins)
Spectrum ana::PredictionNoOsc::PredictComponent ( osc::IOscCalc calc,
Flavors::Flavors_t  flav,
Current::Current_t  curr,
Sign::Sign_t  sign 
) const
overridevirtual

Implements ana::IPrediction.

Definition at line 45 of file PredictionNoOsc.cxx.

References ana::assert(), ana::Spectrum::Clear(), fSpectrum, fSpectrumNC, fSpectrumNue, fSpectrumNuebar, fSpectrumNumu, fSpectrumNumubar, ana::Flavors::kAll, ana::Sign::kAntiNu, ana::Current::kBoth, ana::Sign::kBoth, ana::Current::kCC, ana::Current::kNC, ana::Flavors::kNuEToNuE, ana::Flavors::kNuMuToNuMu, Predict(), and runNovaSAM::ret.

Referenced by calculateComponentsNumu(), calculateWrongSignNue(), calculateWrongSignNumuQ1(), calculateWrongSignNumuQ2(), calculateWrongSignNumuQ3(), calculateWrongSignNumuQ4(), and Predict().

49  {
50  if(flav == Flavors::kAll &&
51  curr == Current::kBoth &&
52  sign == Sign::kBoth)
53  return Predict((osc::IOscCalc*)(nullptr)); // Faster
54 
55  if(curr & Current::kNC){
56  // We don't have NC broken down by sign or flavour
57  assert(flav & Flavors::kAll && sign & Sign::kBoth);
58  return fSpectrumNC;
59  }
60 
62 
63  using namespace Flavors;
64  using namespace Current;
65  using namespace Sign;
66 
68  ret.Clear();
69 
70  // Safe to assume by this point that it's 100% CC
71  if(flav & kNuMuToNuMu && sign & kNu) ret += fSpectrumNumu;
72  if(flav & kNuMuToNuMu && sign & kAntiNu) ret += fSpectrumNumubar;
73  if(flav & kNuEToNuE && sign & kNu) ret += fSpectrumNue;
74  if(flav & kNuEToNuE && sign & kAntiNu) ret += fSpectrumNuebar;
75 
76  return ret;
77  }
Antineutrinos-only.
Definition: IPrediction.h:50
(&#39;beam &#39;)
Definition: IPrediction.h:15
void Clear()
Definition: Spectrum.cxx:433
Charged-current interactions.
Definition: IPrediction.h:39
Interactions of both types.
Definition: IPrediction.h:42
std::vector< float > Spectrum
Definition: Constants.h:527
(&#39; survival&#39;)
Definition: IPrediction.h:19
Neutral-current interactions.
Definition: IPrediction.h:40
assert(nhit_max >=nhit_nbins)
Both neutrinos and antineutrinos.
Definition: IPrediction.h:52
virtual Spectrum Predict(osc::IOscCalc *) const override
All neutrinos, any flavor.
Definition: IPrediction.h:26
def sign(x)
Definition: canMan.py:197
Spectrum ana::IPrediction::PredictComponent ( osc::IOscCalcStan calc,
Flavors::Flavors_t  flav,
Current::Current_t  curr,
Sign::Sign_t  sign 
) const
virtualinherited

Reimplemented in ana::PredictionInterp, ana::PredictionExtendToPeripheral, ana::PredictionExtrapSum, ana::PredictionAddRock, and ana::PredictionExtrap.

Definition at line 93 of file IPrediction.cxx.

References ana::assert().

97  {
98  assert(0 && "This Prediction hasn't implemented a Stan-aware PredictComponent()!");
99  }
assert(nhit_max >=nhit_nbins)
Spectrum ana::IPrediction::PredictComponentSyst ( osc::IOscCalc calc,
const SystShifts syst,
Flavors::Flavors_t  flav,
Current::Current_t  curr,
Sign::Sign_t  sign 
) const
virtualinherited

Reimplemented in ana::PredictionInterp, ana::PredictionCombinePeriods, and ana::PredictionAddRock.

Definition at line 102 of file IPrediction.cxx.

References ana::assert(), ana::SystShifts::IsNominal(), and ana::IPrediction::PredictComponent().

Referenced by ana::PredictionAddRock::_PredictComponentSyst(), ana::CovMxManager::AddSystematic(), ana::GetBFSystBands(), ana::GetSystBands(), PlotSysts(), ana::PredictionCombinePeriods::PredictComponentSyst(), systematics_extrap_comp_from_pred_interp(), systematics_summary_from_pred_interp(), and systematics_table_from_pred_interp().

107  {
108  assert(syst.IsNominal() && "This Prediction doesn't support PredictSyst(). Did you just mean Predict()?");
109 
110  // Default implementation: no treatment of systematics
111  return PredictComponent(calc, flav, curr, sign);
112  }
virtual Spectrum PredictComponent(osc::IOscCalc *calc, Flavors::Flavors_t flav, Current::Current_t curr, Sign::Sign_t sign) const =0
assert(nhit_max >=nhit_nbins)
def sign(x)
Definition: canMan.py:197
Spectrum ana::IPrediction::PredictComponentSyst ( osc::IOscCalcStan calc,
const SystShifts syst,
Flavors::Flavors_t  flav,
Current::Current_t  curr,
Sign::Sign_t  sign 
) const
virtualinherited

Reimplemented in ana::PredictionInterp, and ana::PredictionAddRock.

Definition at line 116 of file IPrediction.cxx.

References ana::assert(), ana::SystShifts::IsNominal(), and ana::IPrediction::PredictComponent().

121  {
122  assert(syst.IsNominal() && "This Prediction doesn't support PredictSyst(). Did you just mean Predict()?");
123 
124  // Default implementation: no treatment of systematics
125  return PredictComponent(calc, flav, curr, sign);
126  }
virtual Spectrum PredictComponent(osc::IOscCalc *calc, Flavors::Flavors_t flav, Current::Current_t curr, Sign::Sign_t sign) const =0
assert(nhit_max >=nhit_nbins)
def sign(x)
Definition: canMan.py:197
Spectrum ana::IPrediction::PredictSyst ( osc::IOscCalc calc,
const SystShifts syst 
) const
virtualinherited

Reimplemented in ana::PredictionInterp, ana::PredictionScaleComp, ana::PredictionCombinePeriods, and ana::PredictionAddRock.

Definition at line 72 of file IPrediction.cxx.

References ana::assert(), ana::SystShifts::IsNominal(), and ana::IPrediction::Predict().

Referenced by BicountEllipse_dCP(), ana::CountingExperiment::ChiSq(), ana::DataMCComparison(), demo_nueNumuSysts(), DrawBins(), DrawPulls(), fill_col(), joint_fit_2017_make_fc_slice(), joint_fit_2017_make_fc_surf(), MakeFakeData(), ana::PlotWithAreaSystErrorBand(), ana::PlotWithSystErrorBand(), ana::PlotWithSystErrorBand_Quant(), ana::SingleSampleExperiment::PredHistIncCosmics(), and ana::PredictionScaleComp::~PredictionScaleComp().

74  {
75  assert(syst.IsNominal() && "This Prediction doesn't support PredictSyst(). Did you just mean Predict()?");
76 
77  // Default implementation: no treatment of systematics
78  return Predict(calc);
79  }
virtual Spectrum Predict(osc::IOscCalc *calc) const =0
assert(nhit_max >=nhit_nbins)
Spectrum ana::IPrediction::PredictSyst ( osc::IOscCalcStan calc,
const SystShifts syst 
) const
virtualinherited

Reimplemented in ana::PredictionInterp, and ana::PredictionAddRock.

Definition at line 82 of file IPrediction.cxx.

References ana::assert(), ana::SystShifts::IsNominal(), and ana::IPrediction::Predict().

84  {
85  assert(syst.IsNominal() && "This Prediction doesn't support PredictSyst(). Did you just mean Predict()?");
86 
87  // Default implementation: no treatment of systematics
88  return Predict(calc);
89  }
virtual Spectrum Predict(osc::IOscCalc *calc) const =0
assert(nhit_max >=nhit_nbins)
Spectrum ana::IPrediction::PredictUnoscillated ( ) const
virtualinherited

Reimplemented in ana::PredictionSterile.

Definition at line 56 of file IPrediction.cxx.

References noosc, and ana::IPrediction::Predict().

Referenced by cc(), demo5(), demo::DrawUpDownRatioCanvas(), efficiency(), efficiencySA(), plot_nd_data_mc(), template_basic(), and test_ana().

57  {
58  // Default implementation
60  return Predict(&noosc);
61  }
Pass neutrinos through unchanged.
Definition: IOscCalc.h:40
virtual Spectrum Predict(osc::IOscCalc *calc) const =0
osc::NoOscillations noosc
Definition: NueSystFuncs.h:44
void ana::PredictionNoOsc::SaveTo ( TDirectory *  dir,
const std::string &  name 
) const
overridevirtual

Reimplemented from ana::IPrediction.

Definition at line 80 of file PredictionNoOsc.cxx.

References dir, fSpectrum, fSpectrumNC, fSpectrumNue, fSpectrumNuebar, fSpectrumNumu, fSpectrumNumubar, ana::Spectrum::SaveTo(), and tmp.

81  {
82  TDirectory* tmp = gDirectory;
83 
84  dir = dir->mkdir(name.c_str()); // switch to subdir
85  dir->cd();
86 
87  TObjString("PredictionNoOsc").Write("type");
88 
89  fSpectrum.SaveTo(dir, "spect");
90  fSpectrumNC.SaveTo(dir, "spect_nc");
91  fSpectrumNumu.SaveTo(dir, "spect_numu");
92  fSpectrumNumubar.SaveTo(dir, "spect_numubar");
93  fSpectrumNue.SaveTo(dir, "spect_nue");
94  fSpectrumNuebar.SaveTo(dir, "spect_nuebar");
95 
96  dir->Write();
97  delete dir;
98 
99  tmp->cd();
100  }
const XML_Char * name
Definition: expat.h:151
Float_t tmp
Definition: plot.C:36
void SaveTo(TDirectory *dir, const std::string &name) const
Definition: Spectrum.cxx:578
TDirectory * dir
Definition: macro.C:5

Member Data Documentation

Spectrum ana::PredictionNoOsc::fSpectrum
protected

Definition at line 59 of file PredictionNoOsc.h.

Referenced by Predict(), PredictComponent(), and SaveTo().

Spectrum ana::PredictionNoOsc::fSpectrumNC
protected

Definition at line 61 of file PredictionNoOsc.h.

Referenced by PredictComponent(), and SaveTo().

Spectrum ana::PredictionNoOsc::fSpectrumNue
protected

Definition at line 64 of file PredictionNoOsc.h.

Referenced by PredictComponent(), and SaveTo().

Spectrum ana::PredictionNoOsc::fSpectrumNuebar
protected

Definition at line 65 of file PredictionNoOsc.h.

Referenced by PredictComponent(), and SaveTo().

Spectrum ana::PredictionNoOsc::fSpectrumNumu
protected

Definition at line 62 of file PredictionNoOsc.h.

Referenced by PredictComponent(), and SaveTo().

Spectrum ana::PredictionNoOsc::fSpectrumNumubar
protected

Definition at line 63 of file PredictionNoOsc.h.

Referenced by PredictComponent(), and SaveTo().


The documentation for this class was generated from the following files: