Sum MC predictions from different periods scaled according to data POT targets. More...

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/CAFAna/Prediction/PredictionCombinePeriods.h"

Inheritance diagram for ana::PredictionCombinePeriods:

Public Member Functions
	PredictionCombinePeriods (const std::vector< std::pair< const IPrediction *, double >> &preds)

	~PredictionCombinePeriods ()

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

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

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

Spectrum	PredictComponentSyst (osc::IOscCalc *calc, const SystShifts &syst, 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

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

unsigned int	NPredictions () const

const IPrediction *	GetPrediction (int i) const

double	GetPOT (int i) const

virtual Spectrum	PredictUnoscillated () const

virtual Spectrum	Predict (osc::IOscCalcStan *calc) 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::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< PredictionCombinePeriods >	LoadFrom (TDirectory *dir, const std::string &name)

Protected Attributes
std::vector< std::pair< const IPrediction *, double > >	fPreds

Detailed Description

Sum MC predictions from different periods scaled according to data POT targets.

Be sure you understand how putting one of these in a SingleSampleExperiment differs from putting regular predictions for the different periods into a MultiExperiment. This calculates chisq(data1+data2|mc1+mc2), whereas the MultiExperiment is chisq(data1|mc1)+chisq(data2+mc2).

Definition at line 14 of file PredictionCombinePeriods.h.

Constructor & Destructor Documentation

ana::PredictionCombinePeriods::PredictionCombinePeriods ( const std::vector< std::pair< const IPrediction *, double >> & preds )

Parameters

preds List of predictions to combine, along with the POT target for that period. Should be like: {{pred1, data1.POT()}, {pred2, data2.POT()}}

Definition at line 18 of file PredictionCombinePeriods.cxx.

References ana::assert(), and fPreds.

Referenced by LoadFrom().

     : fPreds(preds)
   {
     assert(!fPreds.empty());
   }

ana::PredictionCombinePeriods::~PredictionCombinePeriods ( )

Definition at line 25 of file PredictionCombinePeriods.cxx.

   {
     // It's pretty normal for the elements of fPred to have come from the stack
     // for(auto it: fPreds) delete it.first;
   }

Member Function Documentation

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

overridevirtual

Reimplemented from ana::IPrediction.

Definition at line 105 of file PredictionCombinePeriods.cxx.

References ana::bins, ana::IPrediction::ComponentCC(), fPreds, ana::ReweightableSpectrum::GetBinnings(), ana::ReweightableSpectrum::GetEigen(), ana::ReweightableSpectrum::GetLabels(), make_syst_table_plots::h, it, PandAna.Demos.pi0_spectra::labels, and plot_validation_datamc::pred.

   {
     // TODO this is a total copy-paste job from the above.
 
     Eigen::MatrixXd asum;
 
     std::vector<std::string> labels;
     std::vector<Binning> bins;
 
     double potsum = 0;
     double livesum = 0;
     for(auto it: fPreds){
       const IPrediction* pred = it.first;
       const double targetPOT = it.second;
       const OscillatableSpectrum s = pred->ComponentCC(from, to);
 
       Eigen::MatrixXd h = s.GetEigen(targetPOT);
       potsum += targetPOT;
       // No. FD MC is going to have zero (unknown) livetime anyway. If someone
       // needs this livetime accounting, they can come and think hard about it
       // and then reinstate it.
       // livesum += s.Livetime()*targetPOT/s.POT();
 
       if(asum.size() > 0){
         asum += h;
       }
       else{
         asum = h;
         labels = s.GetLabels();
         bins = s.GetBinnings();
       }
     }
 
     return OscillatableSpectrum(std::move(asum),
                                 HistAxis(labels, bins),
                                 potsum, livesum);
   }

Spectrum ana::PredictionCombinePeriods::ComponentNC ( ) const

overridevirtual

Reimplemented from ana::IPrediction.

Definition at line 150 of file PredictionCombinePeriods.cxx.

References ana::Flavors::kAll, ana::Current::kNC, ana::Sign::kNu, noosc, and PredictComponent().

   {
     osc::NoOscillations noosc;
     return PredictComponent(&noosc, Flavors::kAll, Current::kNC, Sign::kNu);
   }

Spectrum ana::PredictionCombinePeriods::ComponentNCAnti ( ) const

overridevirtual

Reimplemented from ana::IPrediction.

Definition at line 155 of file PredictionCombinePeriods.cxx.

References ana::Flavors::kAll, ana::Sign::kAntiNu, ana::Current::kNC, noosc, and PredictComponent().

   {
     osc::NoOscillations noosc;
     return PredictComponent(&noosc, Flavors::kAll, Current::kNC, Sign::kAntiNu);
   }

Spectrum ana::PredictionCombinePeriods::ComponentNCTotal ( ) const

overridevirtual

Reimplemented from ana::IPrediction.

Definition at line 145 of file PredictionCombinePeriods.cxx.

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

   {
     osc::NoOscillations noosc;
     return PredictComponent(&noosc, Flavors::kAll, Current::kNC, Sign::kBoth);
   }

double ana::PredictionCombinePeriods::GetPOT ( int i ) const

inline

Definition at line 57 of file PredictionCombinePeriods.h.

References fPreds, and MECModelEnuComparisons::i.

57 {return fPreds[i].second;}

MECModelEnuComparisons.i

int i

Definition: MECModelEnuComparisons.py:63

ana::PredictionCombinePeriods::fPreds

std::vector< std::pair< const IPrediction *, double > > fPreds

Definition: PredictionCombinePeriods.h:60

const IPrediction* ana::PredictionCombinePeriods::GetPrediction ( int i ) const

inline

Definition at line 56 of file PredictionCombinePeriods.h.

References fPreds, and MECModelEnuComparisons::i.

Referenced by PlotSysts().

56 {return fPreds[i].first;}

MECModelEnuComparisons.i

int i

Definition: MECModelEnuComparisons.py:63

ana::PredictionCombinePeriods::fPreds

std::vector< std::pair< const IPrediction *, double > > fPreds

Definition: PredictionCombinePeriods.h:60

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

static

Definition at line 191 of file PredictionCombinePeriods.cxx.

References ana::assert(), dir, genie::utils::style::Format(), MECModelEnuComparisons::i, ana::LoadFrom< IPrediction >(), pots, plot_validation_datamc::pred, PredictionCombinePeriods(), runNovaSAM::release, and getGoodRuns4SAM::tag.

Referenced by BlessedPlotsAnaByPeriod(), CalcRWithSysts(), FDDataMCSystBandAna(), FitParamEffectsAna(), FitSystEffectsAna(), LoadMaps(), SaveTo(), SpectrumParamEffectsAna(), and SystsOscParamAna().

   {
     dir = dir->GetDirectory(name.c_str()); // switch to subdir
     assert(dir);
 
     TObjString* tag = (TObjString*)dir->Get("type");
     assert(tag);
     assert(tag->GetString() == "PredictionCombinePeriods");
     delete tag;
 
     TVectorD* pots = (TVectorD*)dir->Get("pots");
     assert(pots);
 
     std::vector<std::pair<const IPrediction*, double>> preds;
 
     for(int i = 0; i < pots->GetNrows(); ++i){
       const IPrediction* pred = ana::LoadFrom<IPrediction>(dir, TString::Format("pred%d", i).Data()).release();
       assert(pred);
       preds.emplace_back(pred, (*pots)[i]);
     }
 
     delete dir;
 
     return std::unique_ptr<PredictionCombinePeriods>(new PredictionCombinePeriods(preds));
   }

unsigned int ana::PredictionCombinePeriods::NPredictions ( ) const

inline

Definition at line 55 of file PredictionCombinePeriods.h.

References fPreds.

55 {return fPreds.size();}

ana::PredictionCombinePeriods::fPreds

std::vector< std::pair< const IPrediction *, double > > fPreds

Definition: PredictionCombinePeriods.h:60

Spectrum ana::PredictionCombinePeriods::Predict ( osc::IOscCalc * calc ) const

overridevirtual

Implements ana::IPrediction.

Definition at line 32 of file PredictionCombinePeriods.cxx.

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

Referenced by drawSystsShiftingNDdata(), futureSig_reach_singlePOTcombo_syst(), MakeNusPrediction(), PlotStack(), PlotSysts(), saveContours_complete(), and SpectrumParamEffectsAna().

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

Spectrum ana::IPrediction::Predict ( osc::IOscCalcStan * calc ) const

virtualinherited

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

Definition at line 43 of file IPrediction.cxx.

References ana::assert().

   {
     assert(0 && "This Prediction hasn't implemented a Stan-aware Predict()!");
   }

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

overridevirtual

Implements ana::IPrediction.

Definition at line 52 of file PredictionCombinePeriods.cxx.

References ana::SystShifts::Nominal(), and PredictComponentSyst().

Referenced by CalcRWithSysts(), ComponentNC(), ComponentNCAnti(), ComponentNCTotal(), PlotStack(), PlotSysts(), Predict(), and PredictSyst().

   {
     return PredictComponentSyst(calc, SystShifts::Nominal(),
                                 flav, curr, sign);
   }

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::PredictionNoOsc, ana::PredictionExtendToPeripheral, ana::PredictionExtrapSum, ana::PredictionAddRock, and ana::PredictionExtrap.

Definition at line 70 of file IPrediction.cxx.

References ana::assert().

   {
     assert(0 && "This Prediction hasn't implemented a Stan-aware PredictComponent()!");
   }

Spectrum ana::PredictionCombinePeriods::PredictComponentSyst	(	osc::IOscCalc *	calc,
		const SystShifts &	syst,
		Flavors::Flavors_t	flav,
		Current::Current_t	curr,
		Sign::Sign_t	sign
	)		const

overridevirtual

Reimplemented from ana::IPrediction.

Definition at line 63 of file PredictionCombinePeriods.cxx.

References ana::bins, fPreds, ana::Spectrum::GetBinnings(), ana::Spectrum::GetEigen(), ana::Spectrum::GetLabels(), make_syst_table_plots::h, it, PandAna.Demos.pi0_spectra::labels, plot_validation_datamc::pred, and ana::IPrediction::PredictComponentSyst().

Referenced by CalcRWithSysts(), PredictComponent(), and PredictSyst().

   {
     Eigen::ArrayXd asum;
 
     std::vector<std::string> labels;
     std::vector<Binning> bins;
 
     double potsum = 0;
     double livesum = 0;
     for(auto it: fPreds){
       const IPrediction* pred = it.first;
       const double targetPOT = it.second;
       const Spectrum s = pred->PredictComponentSyst(calc, syst,
                                                     flav, curr, sign);
 
       Eigen::ArrayXd h = s.GetEigen(targetPOT);
       potsum += targetPOT;
       // No. FD MC is going to have zero (unknown) livetime anyway. If someone
       // needs this livetime accounting, they can come and think hard about it
       // and then reinstate it.
       // livesum += s.Livetime()*targetPOT/s.POT();
 
       if(asum.size() > 0){
         asum += h;
       }
       else{
         asum = h;
         labels = s.GetLabels();
         bins = s.GetBinnings();
       }
     }
 
     return Spectrum(std::move(asum),
                     HistAxis(labels, bins),
                     potsum, livesum);
   }

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 93 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::PredictionCombinePeriods::PredictSyst	(	osc::IOscCalc *	calc,
		const SystShifts &	syst
	)		const

overridevirtual

Reimplemented from ana::IPrediction.

Definition at line 41 of file PredictionCombinePeriods.cxx.

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

Referenced by PlotSysts().

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

Spectrum ana::IPrediction::PredictSyst	(	osc::IOscCalcStan *	calc,
		const SystShifts &	syst
	)		const

virtualinherited

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

Definition at line 59 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 33 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::PredictionCombinePeriods::SaveTo	(	TDirectory *	dir,
		const std::string &	name
	)		const

overridevirtual

Reimplemented from ana::IPrediction.

Definition at line 164 of file PredictionCombinePeriods.cxx.

References dir, genie::utils::style::Format(), fPreds, MECModelEnuComparisons::i, compare_h5_caf::idx, it, LoadFrom(), pots, moon_position_table_new3::second, and tmp.

Referenced by MakeNusPrediction().

   {
     TDirectory* tmp = gDirectory;
 
     dir = dir->mkdir(name.c_str()); // switch to subdir
     dir->cd();
 
     TObjString("PredictionCombinePeriods").Write("type");
 
     int idx = 0;
     for(auto it: fPreds)
       it.first->SaveTo(dir, TString::Format("pred%d", idx++).Data());
 
     TVectorD pots(fPreds.size());
     for(unsigned int i = 0; i < fPreds.size(); ++i)
       pots[i] = fPreds[i].second;
 
     pots.Write("pots");
 
     dir->Write();
     delete dir;
 
     tmp->cd();
   }

Member Data Documentation

std::vector<std::pair<const IPrediction*, double> > ana::PredictionCombinePeriods::fPreds

protected

Definition at line 60 of file PredictionCombinePeriods.h.

Referenced by ComponentCC(), GetPOT(), GetPrediction(), NPredictions(), PredictComponentSyst(), PredictionCombinePeriods(), and SaveTo().

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

/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/CAFAna/Prediction/PredictionCombinePeriods.h
/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/CAFAna/Prediction/PredictionCombinePeriods.cxx

Public Member Functions

Static Public Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation