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

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

Inheritance diagram for ana::PredictionExtrap:

Public Member Functions
	PredictionExtrap (std::unique_ptr< IExtrap > &&extrap)
	Takes ownership of extrap. More...

	PredictionExtrap ()=delete

Spectrum	Predict (osc::IOscCalc *calc) const override

Spectrum	Predict (osc::IOscCalcStan *calc) const override

Spectrum	PredictComponent (osc::IOscCalc *calc, Flavors::Flavors_t flav, Current::Current_t curr, Sign::Sign_t sign) const override

Spectrum	PredictComponent (osc::IOscCalcStan *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 std::string &name) const override

const IExtrap *	GetExtrap () const

virtual Spectrum	PredictUnoscillated () const

virtual Spectrum	PredictSyst (osc::IOscCalc *calc, const SystShifts &syst) const

virtual Spectrum	PredictSyst (osc::IOscCalcStan *calc, const SystShifts &syst) 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

Static Public Member Functions
static std::unique_ptr< PredictionExtrap >	LoadFrom (TDirectory *dir, const std::string &name)

Protected Member Functions
template<typename T >
Spectrum	_PredictComponent (osc::_IOscCalc< 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...

Protected Attributes
std::unique_ptr< IExtrap >	fExtrap

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.

     : fExtrap(std::move(extrap))
   {
   }

ana::PredictionExtrap::PredictionExtrap ( )

delete

Referenced by LoadFrom().

Member Function Documentation

template<typename T >

Spectrum ana::PredictionExtrap::_PredictComponent	(	osc::_IOscCalc< T > *	calc,
		Flavors::Flavors_t	flav,
		Current::Current_t	curr,
		Sign::Sign_t	sign
	)		const

protected

Templated helper function called by the non-templated versions.

Definition at line 40 of file PredictionExtrap.cxx.

References ana::assert(), ana::Spectrum::Clear(), 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.

Referenced by GetExtrap(), and PredictComponent().

   {
     Spectrum ret = fExtrap->NCTotalComponent(); // Get binning
     ret.Clear();
 
     if(curr & Current::kCC){
       if(flav & Flavors::kNuEToNuE    && sign & Sign::kNu)     ret += fExtrap->NueSurvComponent().    Oscillated(calc, +12, +12);
       if(flav & Flavors::kNuEToNuE    && sign & Sign::kAntiNu) ret += fExtrap->AntiNueSurvComponent().Oscillated(calc, -12, -12);
 
       if(flav & Flavors::kNuEToNuMu   && sign & Sign::kNu)     ret += fExtrap->NumuAppComponent().    Oscillated(calc, +12, +14);
       if(flav & Flavors::kNuEToNuMu   && sign & Sign::kAntiNu) ret += fExtrap->AntiNumuAppComponent().Oscillated(calc, -12, -14);
 
       if(flav & Flavors::kNuEToNuTau  && sign & Sign::kNu)     ret += fExtrap->TauFromEComponent().    Oscillated(calc, +12, +16);
       if(flav & Flavors::kNuEToNuTau  && sign & Sign::kAntiNu) ret += fExtrap->AntiTauFromEComponent().Oscillated(calc, -12, -16);
 
       if(flav & Flavors::kNuMuToNuE   && sign & Sign::kNu)     ret += fExtrap->NueAppComponent().    Oscillated(calc, +14, +12);
       if(flav & Flavors::kNuMuToNuE   && sign & Sign::kAntiNu) ret += fExtrap->AntiNueAppComponent().Oscillated(calc, -14, -12);
 
       if(flav & Flavors::kNuMuToNuMu  && sign & Sign::kNu)     ret += fExtrap->NumuSurvComponent().    Oscillated(calc, +14, +14);
       if(flav & Flavors::kNuMuToNuMu  && sign & Sign::kAntiNu) ret += fExtrap->AntiNumuSurvComponent().Oscillated(calc, -14, -14);
 
       if(flav & Flavors::kNuMuToNuTau && sign & Sign::kNu)     ret += fExtrap->TauFromMuComponent().    Oscillated(calc, +14, +16);
       if(flav & Flavors::kNuMuToNuTau && sign & Sign::kAntiNu) ret += fExtrap->AntiTauFromMuComponent().Oscillated(calc, -14, -16);
     }
     if(curr & Current::kNC){
     assert(flav == Flavors::kAll); // Don't know how to calculate anything else
 
       if(sign & Sign::kNu)     ret += fExtrap->NCComponent();
       if(sign & Sign::kAntiNu) ret += fExtrap->NCAntiComponent();
     }
 
     return ret;
   }

OscillatableSpectrum ana::PredictionExtrap::ComponentCC	(	int	from,
		int	to
	)		const

overridevirtual

Reimplemented from ana::IPrediction.

Definition at line 98 of file PredictionExtrap.cxx.

References ana::assert(), and fExtrap.

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

   {
     if(from == +12 && to == +12) return fExtrap->NueSurvComponent();
     if(from == -12 && to == -12) return fExtrap->AntiNueSurvComponent();
 
     if(from == +12 && to == +14) return fExtrap->NumuAppComponent();
     if(from == -12 && to == -14) return fExtrap->AntiNumuAppComponent();
 
     if(from == +12 && to == +16) return fExtrap->TauFromEComponent();
     if(from == -12 && to == -16) return fExtrap->AntiTauFromEComponent();
 
     if(from == +14 && to == +12) return fExtrap->NueAppComponent();
     if(from == -14 && to == -12) return fExtrap->AntiNueAppComponent();
 
     if(from == +14 && to == +14) return fExtrap->NumuSurvComponent();
     if(from == -14 && to == -14) return fExtrap->AntiNumuSurvComponent();
 
     if(from == +14 && to == +16) return fExtrap->TauFromMuComponent();
     if(from == -14 && to == -16) return fExtrap->AntiTauFromMuComponent();
 
     assert(0 && "Not reached");
   }

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 nue_decomp_scales().

   {
     return fExtrap->NCComponent();
   }

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().

   {
     return fExtrap->NCAntiComponent();
   }

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().

   {
     return fExtrap->NCTotalComponent();
   }

const IExtrap* ana::PredictionExtrap::GetExtrap ( ) const

inline

Definition at line 43 of file PredictionExtrap.h.

References _PredictComponent(), and fExtrap.

Referenced by ana::GetDecomp(), and make_files_for_decomp_plots().

43 {return fExtrap.get();}

ana::PredictionExtrap::fExtrap

std::unique_ptr< IExtrap > fExtrap

Definition: PredictionExtrap.h:53

std::unique_ptr< PredictionExtrap > ana::PredictionExtrap::LoadFrom	(	TDirectory *	dir,
		const std::string &	name
	)

static

Definition at line 155 of file PredictionExtrap.cxx.

References ana::assert(), dir, ana::LoadFrom< IExtrap >(), PredictionExtrap(), and runNovaSAM::ret.

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

   {
     dir = dir->GetDirectory(name.c_str()); // switch to subdir
     assert(dir);
 
     std::unique_ptr<PredictionExtrap> ret(new PredictionExtrap(ana::LoadFrom<IExtrap>(dir, "extrap")));
 
     delete dir;
 
     return ret;
   }

Spectrum ana::PredictionExtrap::Predict ( osc::IOscCalc * 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_Data_PosComp(), FD_plots(), GetExtrap(), GetNuePredictionFuture(), getTimePeakPlots(), MakePeriCutPlot(), MichelDecompTest(), modularextrap_demo_nue(), modularextrap_demo_numu(), numu_sig_nonmax(), PeripheralCuts(), plot_nue_xsec_pred(), PlotNus17Prediction(), ana::PredictionSterile::PredictUnoscillated(), PredRatioToNom(), resolution2018(), saveContours_addExpt(), saveContours_simple(), sterile_demo(), template_basic(), test_ana(), test_genie_systs(), test_nueextrapsyst(), and test_saloaders().

   {
     return PredictComponent(calc,
                             Flavors::kAll,
                             Current::kBoth,
                             Sign::kBoth);
   }

Spectrum ana::PredictionExtrap::Predict ( osc::IOscCalcStan * 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().

   {
     return PredictComponent(calc,
                             Flavors::kAll,
                             Current::kBoth,
                             Sign::kBoth);
   }

Spectrum ana::PredictionExtrap::PredictComponent	(	osc::IOscCalc *	calc,
		Flavors::Flavors_t	flav,
		Current::Current_t	curr,
		Sign::Sign_t	sign
	)		const

overridevirtual

Implements ana::IPrediction.

Reimplemented in ana::PredictionSterile.

Definition at line 79 of file PredictionExtrap.cxx.

References _PredictComponent().

Referenced by BackPredict(), CompareMCCompPrediction(), CVNCuts(), GetExtrap(), getTimePeakPlots(), MakePeriCutPlot(), MCCompPredictionTable(), MichelDecompTest(), nue_pid_effs(), nue_pid_effs_miniprod(), nue_pid_effs_paper_numu_neweff(), PeripheralCuts(), Predict(), ana::PredictionSterile::PredictComponent(), ana::PredictionAddRock::RockComponent(), and test_saloaders().

   {
     return _PredictComponent(calc, flav, curr, sign);
   }

Spectrum ana::PredictionExtrap::PredictComponent	(	osc::IOscCalcStan *	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 _PredictComponent().

   {
     return _PredictComponent(calc, flav, curr, sign);
   }

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().

   {
     assert(syst.IsNominal() && "This Prediction doesn't support PredictSyst(). Did you just mean Predict()?");
 
     // Default implementation: no treatment of systematics
     return PredictComponent(calc, flav, curr, sign);
   }

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().

   {
     assert(syst.IsNominal() && "This Prediction doesn't support PredictSyst(). Did you just mean Predict()?");
 
     // Default implementation: no treatment of systematics
     return PredictComponent(calc, flav, curr, sign);
   }

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().

   {
     assert(syst.IsNominal() && "This Prediction doesn't support PredictSyst(). Did you just mean Predict()?");
 
     // Default implementation: no treatment of systematics
     return Predict(calc);
   }

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().

   {
     assert(syst.IsNominal() && "This Prediction doesn't support PredictSyst(). Did you just mean Predict()?");
 
     // Default implementation: no treatment of systematics
     return Predict(calc);
   }

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().

   {
     // Default implementation
     osc::NoOscillations noosc;
     return Predict(&noosc);
   }

void ana::PredictionExtrap::SaveTo	(	TDirectory *	dir,
		const std::string &	name
	)		const

overridevirtual

Reimplemented from ana::IPrediction.

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

Definition at line 137 of file PredictionExtrap.cxx.

References dir, fExtrap, and tmp.

Referenced by make_nue_filesyst_pred(), 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_numu(), and ana::PredictionSterile::SaveTo().

   {
     TDirectory* tmp = gDirectory;
 
     dir = dir->mkdir(name.c_str()); // switch to subdir
     dir->cd();
 
     TObjString("PredictionExtrap").Write("type");
 
     fExtrap->SaveTo(dir, "extrap");
 
     dir->Write();
     delete dir;
 
     tmp->cd();
   }

Member Data Documentation

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

protected

Definition at line 53 of file PredictionExtrap.h.

Referenced by _PredictComponent(), ComponentCC(), ComponentNC(), ComponentNCAnti(), ComponentNCTotal(), GetExtrap(), ana::PredictionSterile::PredictComponent(), SaveTo(), and ana::PredictionNoExtrap::SaveTo().

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

/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-09-26/CAFAna/Prediction/PredictionExtrap.h
/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-09-26/CAFAna/Prediction/PredictionExtrap.cxx

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation