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

Take the output of an extrapolation and oscillate it as required. More...

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-02-16/CAFAna/Prediction/PredictionExtrap.h"

Inheritance diagram for ana::PredictionExtrap:
ana::IPrediction ana::PredictionNoExtrap ana::PredictionSterile

Public Member Functions

 PredictionExtrap (std::unique_ptr< IExtrap > &&extrap)
 Takes ownership of extrap. More...
 
 PredictionExtrap ()=delete
 
Spectrum Predict (osc::IOscCalculator *calc) const override
 
SpectrumStan Predict (osc::IOscCalculatorStan *calc) const override
 
Spectrum PredictComponent (osc::IOscCalculator *calc, Flavors::Flavors_t flav, Current::Current_t curr, Sign::Sign_t sign) const override
 
SpectrumStan PredictComponent (osc::IOscCalculatorStan *calc, Flavors::Flavors_t flav, Current::Current_t curr, Sign::Sign_t sign) const override
 
OscillatableSpectrum ComponentCC (int from, int to) const override
 
Spectrum ComponentNCTotal () const override
 
Spectrum ComponentNC () const override
 
Spectrum ComponentNCAnti () const override
 
virtual void SaveTo (TDirectory *dir) const override
 
const IExtrapGetExtrap () const
 
virtual Spectrum PredictUnoscillated () const
 
virtual Spectrum PredictSyst (osc::IOscCalculator *calc, const SystShifts &syst) const
 
virtual SpectrumStan PredictSyst (osc::IOscCalculatorStan *calc, const SystShifts &syst) const
 
virtual Spectrum PredictComponentSyst (osc::IOscCalculator *calc, const SystShifts &syst, Flavors::Flavors_t flav, Current::Current_t curr, Sign::Sign_t sign) const
 
virtual SpectrumStan PredictComponentSyst (osc::IOscCalculatorStan *calc, const SystShifts &syst, Flavors::Flavors_t flav, Current::Current_t curr, Sign::Sign_t sign) const
 
virtual void Derivative (osc::IOscCalculator *calc, const SystShifts &shift, double pot, std::unordered_map< const ISyst *, std::vector< double >> &dchi) const
 

Static Public Member Functions

static std::unique_ptr< PredictionExtrapLoadFrom (TDirectory *dir)
 

Protected Attributes

std::unique_ptr< IExtrapfExtrap
 

Private Member Functions

template<typename U , typename T >
_PredictComponent (osc::_IOscCalculator< T > *calc, Flavors::Flavors_t flav, Current::Current_t curr, Sign::Sign_t sign) const
 Templated helper function called by the non-templated versions. More...
 

Detailed Description

Take the output of an extrapolation and oscillate it as required.

Definition at line 12 of file PredictionExtrap.h.

Constructor & Destructor Documentation

ana::PredictionExtrap::PredictionExtrap ( std::unique_ptr< IExtrap > &&  extrap)

Takes ownership of extrap.

Definition at line 14 of file PredictionExtrap.cxx.

15  : fExtrap(std::move(extrap))
16  {
17  }
std::unique_ptr< IExtrap > fExtrap
ana::PredictionExtrap::PredictionExtrap ( )
delete

Referenced by LoadFrom().

Member Function Documentation

template<typename U , typename T >
U ana::PredictionExtrap::_PredictComponent ( osc::_IOscCalculator< T > *  calc,
Flavors::Flavors_t  flav,
Current::Current_t  curr,
Sign::Sign_t  sign 
) const
private

Templated helper function called by the non-templated versions.

Definition at line 40 of file PredictionExtrap.cxx.

References fExtrap, ana::Flavors::kAll, ana::Sign::kAntiNu, ana::Current::kCC, ana::Current::kNC, ana::Sign::kNu, ana::Flavors::kNuEToNuE, ana::Flavors::kNuEToNuMu, ana::Flavors::kNuEToNuTau, ana::Flavors::kNuMuToNuE, ana::Flavors::kNuMuToNuMu, ana::Flavors::kNuMuToNuTau, and runNovaSAM::ret.

44  {
45  U ret = fExtrap->NCTotalComponent(); // Get binning
46  ret.Clear();
47 
48  if(curr & Current::kCC){
49  if(flav & Flavors::kNuEToNuE && sign & Sign::kNu) ret += fExtrap->NueSurvComponent(). Oscillated(calc, +12, +12);
50  if(flav & Flavors::kNuEToNuE && sign & Sign::kAntiNu) ret += fExtrap->AntiNueSurvComponent().Oscillated(calc, -12, -12);
51 
52  if(flav & Flavors::kNuEToNuMu && sign & Sign::kNu) ret += fExtrap->NumuAppComponent(). Oscillated(calc, +12, +14);
53  if(flav & Flavors::kNuEToNuMu && sign & Sign::kAntiNu) ret += fExtrap->AntiNumuAppComponent().Oscillated(calc, -12, -14);
54 
55  if(flav & Flavors::kNuEToNuTau && sign & Sign::kNu) ret += fExtrap->TauFromEComponent(). Oscillated(calc, +12, +16);
56  if(flav & Flavors::kNuEToNuTau && sign & Sign::kAntiNu) ret += fExtrap->AntiTauFromEComponent().Oscillated(calc, -12, -16);
57 
58  if(flav & Flavors::kNuMuToNuE && sign & Sign::kNu) ret += fExtrap->NueAppComponent(). Oscillated(calc, +14, +12);
59  if(flav & Flavors::kNuMuToNuE && sign & Sign::kAntiNu) ret += fExtrap->AntiNueAppComponent().Oscillated(calc, -14, -12);
60 
61  if(flav & Flavors::kNuMuToNuMu && sign & Sign::kNu) ret += fExtrap->NumuSurvComponent(). Oscillated(calc, +14, +14);
62  if(flav & Flavors::kNuMuToNuMu && sign & Sign::kAntiNu) ret += fExtrap->AntiNumuSurvComponent().Oscillated(calc, -14, -14);
63 
64  if(flav & Flavors::kNuMuToNuTau && sign & Sign::kNu) ret += fExtrap->TauFromMuComponent(). Oscillated(calc, +14, +16);
65  if(flav & Flavors::kNuMuToNuTau && sign & Sign::kAntiNu) ret += fExtrap->AntiTauFromMuComponent().Oscillated(calc, -14, -16);
66  }
67  if(curr & Current::kNC){
68  assert(flav == Flavors::kAll); // Don't know how to calculate anything else
69 
70  if(sign & Sign::kNu) ret += fExtrap->NCComponent();
71  if(sign & Sign::kAntiNu) ret += fExtrap->NCAntiComponent();
72  }
73 
74  return ret;
75  }
Antineutrinos-only.
Definition: IPrediction.h:50
(&#39; appearance&#39;)
Definition: IPrediction.h:18
(&#39;beam &#39;)
Definition: IPrediction.h:15
Charged-current interactions.
Definition: IPrediction.h:39
Neutrinos-only.
Definition: IPrediction.h:49
def sign(x)
Definition: canMan.py:204
(&#39; survival&#39;)
Definition: IPrediction.h:19
std::unique_ptr< IExtrap > fExtrap
Neutral-current interactions.
Definition: IPrediction.h:40
All neutrinos, any flavor.
Definition: IPrediction.h:26
(&#39; appearance&#39;)
Definition: IPrediction.h:16
OscillatableSpectrum ana::PredictionExtrap::ComponentCC ( int  from,
int  to 
) const
overridevirtual

Reimplemented from ana::IPrediction.

Definition at line 98 of file PredictionExtrap.cxx.

References fExtrap.

Referenced by ana::PredictionAddRock::ComponentCC(), and ana::PredictionExtendToPeripheral::PredictionExtendToPeripheral().

99  {
100  if(from == +12 && to == +12) return fExtrap->NueSurvComponent();
101  if(from == -12 && to == -12) return fExtrap->AntiNueSurvComponent();
102 
103  if(from == +12 && to == +14) return fExtrap->NumuAppComponent();
104  if(from == -12 && to == -14) return fExtrap->AntiNumuAppComponent();
105 
106  if(from == +12 && to == +16) return fExtrap->TauFromEComponent();
107  if(from == -12 && to == -16) return fExtrap->AntiTauFromEComponent();
108 
109  if(from == +14 && to == +12) return fExtrap->NueAppComponent();
110  if(from == -14 && to == -12) return fExtrap->AntiNueAppComponent();
111 
112  if(from == +14 && to == +14) return fExtrap->NumuSurvComponent();
113  if(from == -14 && to == -14) return fExtrap->AntiNumuSurvComponent();
114 
115  if(from == +14 && to == +16) return fExtrap->TauFromMuComponent();
116  if(from == -14 && to == -16) return fExtrap->AntiTauFromMuComponent();
117 
118  assert(0 && "Not reached");
119  }
std::unique_ptr< IExtrap > fExtrap
Spectrum ana::PredictionExtrap::ComponentNC ( ) const
overridevirtual

Reimplemented from ana::IPrediction.

Definition at line 127 of file PredictionExtrap.cxx.

References fExtrap.

Referenced by ana::PredictionAddRock::ComponentNC(), and ana::PredictionExtendToPeripheral::PredictionExtendToPeripheral().

128  {
129  return fExtrap->NCComponent();
130  }
std::unique_ptr< IExtrap > fExtrap
Spectrum ana::PredictionExtrap::ComponentNCAnti ( ) const
overridevirtual

Reimplemented from ana::IPrediction.

Definition at line 131 of file PredictionExtrap.cxx.

References fExtrap.

Referenced by ana::PredictionAddRock::ComponentNCAnti(), and ana::PredictionExtendToPeripheral::PredictionExtendToPeripheral().

132  {
133  return fExtrap->NCAntiComponent();
134  }
std::unique_ptr< IExtrap > fExtrap
Spectrum ana::PredictionExtrap::ComponentNCTotal ( ) const
overridevirtual

Reimplemented from ana::IPrediction.

Definition at line 123 of file PredictionExtrap.cxx.

References fExtrap.

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

124  {
125  return fExtrap->NCTotalComponent();
126  }
std::unique_ptr< IExtrap > fExtrap
virtual void ana::IPrediction::Derivative ( osc::IOscCalculator calc,
const SystShifts shift,
double  pot,
std::unordered_map< const ISyst *, std::vector< double >> &  dchi 
) const
inlinevirtualinherited

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

Definition at line 92 of file IPrediction.h.

Referenced by ana::PredictionAddRock::Derivative(), and ana::SingleSampleExperiment::Derivative().

96  {
97  // Implementing this function is optional. If you don't implement it,
98  // this default implementation will be used, which signals to callers
99  // that your Prediction doesn't implement this feature.
100  dchi.clear();
101  }
const IExtrap* ana::PredictionExtrap::GetExtrap ( ) const
inline
std::unique_ptr< PredictionExtrap > ana::PredictionExtrap::LoadFrom ( TDirectory *  dir)
static

Definition at line 151 of file PredictionExtrap.cxx.

References ana::LoadFrom< IExtrap >(), and PredictionExtrap().

Referenced by BlessedPlotsAna(), ana::ComparePredictionsPeripheralFromVector(), ana::PredictionSterile::LoadFrom(), ana::LoadFrom< IPrediction >(), ana::LoadMaps(), make_files_for_decomp_plots(), PlotNus17Prediction(), plots_interactionMode(), resolution2018(), and test_stanfit_statsonly().

152  {
153  assert(dir->GetDirectory("extrap"));
154  return std::unique_ptr<PredictionExtrap>(new PredictionExtrap(ana::LoadFrom<IExtrap>(dir->GetDirectory("extrap"))));
155  }
std::unique_ptr< IExtrap > LoadFrom< IExtrap >(TDirectory *dir)
Definition: IExtrap.cxx:20
TDirectory * dir
Definition: macro.C:5
Spectrum ana::PredictionExtrap::Predict ( osc::IOscCalculator calc) const
overridevirtual

Implements ana::IPrediction.

Definition at line 20 of file PredictionExtrap.cxx.

References ana::Flavors::kAll, ana::Current::kBoth, ana::Sign::kBoth, and PredictComponent().

Referenced by BackPredict(), BlessedPlotsLoad(), caf_numu_sensitivity(), caf_numu_sensitivity_no_tau(), cc(), ComparePredictions(), CVNCuts(), demo5(), demo_CPT(), demo_flat(), drawPlot(), drawSystsShiftingNDdata_updatedAna(), FD_plots(), GetExtrap(), GetNuePredictionFuture(), getTimePeakPlots(), ana::NumuAnalysis::MakeExperiment(), MakePeriCutPlot(), MichelDecompTest(), modularextrap_demo_nue(), modularextrap_demo_numu(), nue_ana_basic(), numu_demo_split_plot(), numu_sig_nonmax(), PeripheralCuts(), plot_nue_xsec_pred(), PlotNus17Prediction(), ana::PredictionSterile::PredictUnoscillated(), PredRatioToNom(), resolution2018(), rhc_sensitivity(), saveContours_addExpt(), saveContours_simple(), sterile_demo(), template_basic(), test_ana(), test_genie_systs(), test_nueextrapsyst(), and test_saloaders().

21  {
22  return PredictComponent(calc,
25  Sign::kBoth);
26  }
Spectrum PredictComponent(osc::IOscCalculator *calc, Flavors::Flavors_t flav, Current::Current_t curr, Sign::Sign_t sign) const override
Interactions of both types.
Definition: IPrediction.h:42
Both neutrinos and antineutrinos.
Definition: IPrediction.h:52
All neutrinos, any flavor.
Definition: IPrediction.h:26
SpectrumStan ana::PredictionExtrap::Predict ( osc::IOscCalculatorStan calc) const
overridevirtual

Reimplemented from ana::IPrediction.

Definition at line 29 of file PredictionExtrap.cxx.

References ana::Flavors::kAll, ana::Current::kBoth, ana::Sign::kBoth, and PredictComponent().

30  {
31  return PredictComponent(calc,
34  Sign::kBoth);
35  }
Spectrum PredictComponent(osc::IOscCalculator *calc, Flavors::Flavors_t flav, Current::Current_t curr, Sign::Sign_t sign) const override
Interactions of both types.
Definition: IPrediction.h:42
Both neutrinos and antineutrinos.
Definition: IPrediction.h:52
All neutrinos, any flavor.
Definition: IPrediction.h:26
Spectrum ana::PredictionExtrap::PredictComponent ( osc::IOscCalculator calc,
Flavors::Flavors_t  flav,
Current::Current_t  curr,
Sign::Sign_t  sign 
) const
overridevirtual
SpectrumStan ana::PredictionExtrap::PredictComponent ( osc::IOscCalculatorStan calc,
Flavors::Flavors_t  flav,
Current::Current_t  curr,
Sign::Sign_t  sign 
) const
overridevirtual

Reimplemented from ana::IPrediction.

Definition at line 89 of file PredictionExtrap.cxx.

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

93  {
94  return _PredictComponent<SpectrumStan>(calc, flav, curr, sign);
95  }
osc::OscCalculatorDumb calc
def sign(x)
Definition: canMan.py:204
Spectrum ana::IPrediction::PredictComponentSyst ( osc::IOscCalculator calc,
const SystShifts syst,
Flavors::Flavors_t  flav,
Current::Current_t  curr,
Sign::Sign_t  sign 
) const
virtualinherited

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

Definition at line 128 of file IPrediction.cxx.

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

Referenced by PlotSysts(), ana::PredictionAddRock::PredictComponentSyst(), ana::PredictionCombinePeriods::PredictComponentSyst(), ana::PredictionNumuFAHadE::PredictComponentSyst(), systematics_extrap_comp_from_pred_interp(), systematics_summary_from_pred_interp(), and systematics_table_from_pred_interp().

133  {
134  assert(syst.IsNominal() && "This Prediction doesn't support PredictSyst(). Did you just mean Predict()?");
135 
136  // Default implementation: no treatment of systematics
137  return PredictComponent(calc, flav, curr, sign);
138  }
virtual Spectrum PredictComponent(osc::IOscCalculator *calc, Flavors::Flavors_t flav, Current::Current_t curr, Sign::Sign_t sign) const =0
def sign(x)
Definition: canMan.py:204
string syst
Definition: plotSysts.py:176
SpectrumStan ana::IPrediction::PredictComponentSyst ( osc::IOscCalculatorStan calc,
const SystShifts syst,
Flavors::Flavors_t  flav,
Current::Current_t  curr,
Sign::Sign_t  sign 
) const
virtualinherited

Reimplemented in ana::PredictionInterp.

Definition at line 142 of file IPrediction.cxx.

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

147  {
148  assert(syst.IsNominal() && "This Prediction doesn't support PredictSyst(). Did you just mean Predict()?");
149 
150  // Default implementation: no treatment of systematics
151  return PredictComponent(calc, flav, curr, sign);
152  }
virtual Spectrum PredictComponent(osc::IOscCalculator *calc, Flavors::Flavors_t flav, Current::Current_t curr, Sign::Sign_t sign) const =0
def sign(x)
Definition: canMan.py:204
string syst
Definition: plotSysts.py:176
Spectrum ana::IPrediction::PredictSyst ( osc::IOscCalculator calc,
const SystShifts syst 
) const
virtualinherited
SpectrumStan ana::IPrediction::PredictSyst ( osc::IOscCalculatorStan calc,
const SystShifts syst 
) const
virtualinherited

Reimplemented in ana::PredictionInterp.

Definition at line 108 of file IPrediction.cxx.

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

110  {
111  assert(syst.IsNominal() && "This Prediction doesn't support PredictSyst(). Did you just mean Predict()?");
112 
113  // Default implementation: no treatment of systematics
114  return Predict(calc);
115  }
virtual Spectrum Predict(osc::IOscCalculator *calc) const =0
string syst
Definition: plotSysts.py:176
Spectrum ana::IPrediction::PredictUnoscillated ( ) const
virtualinherited

Reimplemented in ana::PredictionSterile.

Definition at line 82 of file IPrediction.cxx.

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

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

83  {
84  // Default implementation
86  return Predict(&noosc);
87  }
Pass neutrinos through unchanged.
virtual Spectrum Predict(osc::IOscCalculator *calc) const =0
osc::NoOscillations noosc
Definition: NueSystFuncs.h:44
void ana::PredictionExtrap::SaveTo ( TDirectory *  dir) const
overridevirtual

Reimplemented from ana::IPrediction.

Reimplemented in ana::PredictionNoExtrap, and ana::PredictionSterile.

Definition at line 137 of file PredictionExtrap.cxx.

References fExtrap, and tmp.

Referenced by make_nue_filesyst_pred(), make_nueFDprediction(), make_nueFDprediction_kinematics(), make_nueFDprediction_kinematics_FHC(), make_nueFDprediction_kinematics_FHC_REW(), make_nueFDprediction_kinematics_REW(), make_nueFDprediction_kinematics_RHC(), make_nueFDprediction_kinematics_RHC_REW(), modularextrap_demo_nue(), modularextrap_demo_numu(), and ana::PredictionSterile::SaveTo().

138  {
139  TDirectory* tmp = gDirectory;
140 
141  dir->cd();
142 
143  TObjString("PredictionExtrap").Write("type");
144 
145  fExtrap->SaveTo(dir->mkdir("extrap"));
146 
147  tmp->cd();
148  }
Float_t tmp
Definition: plot.C:36
std::unique_ptr< IExtrap > fExtrap
TDirectory * dir
Definition: macro.C:5

Member Data Documentation

std::unique_ptr<IExtrap> ana::PredictionExtrap::fExtrap
protected

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