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

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-08-04/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::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 IExtrap *	GetExtrap () 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< PredictionExtrap >	LoadFrom (TDirectory *dir)

Protected Attributes
std::unique_ptr< IExtrap >	fExtrap

Private Member Functions
template<typename U , typename T >
U	_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.

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

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

   {
     U 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(), nue_decomp_scales(), and ana::PredictionExtendToPeripheral::PredictionExtendToPeripheral().

   {
     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(), nue_decomp_scales(), and ana::PredictionExtendToPeripheral::PredictionExtendToPeripheral().

   {
     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(), and ana::PredictionExtendToPeripheral::PredictionExtendToPeripheral().

   {
     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(), and ana::PredictionExtendToPeripheral::PredictionExtendToPeripheral().

   {
     return fExtrap->NCTotalComponent();
   }

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

       {
         // Implementing this function is optional. If you don't implement it,
         // this default implementation will be used, which signals to callers
         // that your Prediction doesn't implement this feature.
         dchi.clear();
       }

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

inline

Definition at line 45 of file PredictionExtrap.h.

References fExtrap.

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

45 {return fExtrap.get();}

ana::PredictionExtrap::fExtrap

std::unique_ptr< IExtrap > fExtrap

Definition: PredictionExtrap.h:48

std::unique_ptr< PredictionExtrap > ana::PredictionExtrap::LoadFrom ( TDirectory * dir )

static

Definition at line 151 of file PredictionExtrap.cxx.

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

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

   {
     assert(dir->GetDirectory("extrap"));
     return std::unique_ptr<PredictionExtrap>(new PredictionExtrap(ana::LoadFrom<IExtrap>(dir->GetDirectory("extrap"))));
   }

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_Data_PosComp(), FD_plots(), GetExtrap(), GetNuePredictionFuture(), getTimePeakPlots(), ana::NumuAnalysis::MakeExperiment(), MakePeriCutPlot(), MichelDecompTest(), modularextrap_demo_nue(), modularextrap_demo_numu(), numu_demo_split_plot(), 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);
   }

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

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

Spectrum ana::PredictionExtrap::PredictComponent	(	osc::IOscCalculator *	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 calc, febshutoff_auto::curr, and canMan::sign().

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<Spectrum>(calc, flav, curr, sign);
   }

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

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

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::PredictionCombinePeriods, and ana::PredictionAddRock.

Definition at line 116 of file IPrediction.cxx.

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

Referenced by PlotSysts(), ana::PredictionAddRock::PredictComponentSyst(), 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);
   }

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 130 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::IOscCalculator *	calc,
		const SystShifts &	syst
	)		const

virtualinherited

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

Definition at line 86 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(), fill_col(), ana::GetBFSystBands(), ana::GetSystBands(), joint_fit_2017_make_fc_slice(), joint_fit_2017_make_fc_surf(), ana::SingleSampleExperiment::LogLikelihood(), 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);
   }

SpectrumStan ana::IPrediction::PredictSyst	(	osc::IOscCalculatorStan *	calc,
		const SystShifts &	syst
	)		const

virtualinherited

Reimplemented in ana::PredictionInterp.

Definition at line 96 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 70 of file IPrediction.cxx.

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

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

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

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

   {
     TDirectory* tmp = gDirectory;
 
     dir->cd();
 
     TObjString("PredictionExtrap").Write("type");
 
     fExtrap->SaveTo(dir->mkdir("extrap"));
 
     tmp->cd();
   }

Member Data Documentation

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

protected

Definition at line 48 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-08-04/CAFAna/Prediction/PredictionExtrap.h
/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-08-04/CAFAna/Prediction/PredictionExtrap.cxx

Public Member Functions

Static Public Member Functions

Protected Attributes

Private Member Functions

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation