Prediction broken down into arbitrary components whose scales can be varied independently. More...

#include <PredictionScaleComp.h>

Inheritance diagram for ana::PredictionScaleComp:

Public Member Functions
	PredictionScaleComp (SpectrumLoaderBase &loader, const HistAxis &axis, SpillCut spillcut, Cut cut, const std::vector< const SystComponentScale * > &systs, const SystShifts &shift=kNoShift, const Var &wei=kUnweighted)

	PredictionScaleComp (SpectrumLoaderBase &loader, const HistAxis &axis1, const HistAxis &axis2, SpillCut spillcut, Cut cut, const std::vector< const SystComponentScale * > &systs, const SystShifts &shift=kNoShift, const Var &wei=kUnweighted)
	Constructor to take two HistAxis's to weight 2D spectra. More...

	PredictionScaleComp (SpectrumLoaderBase &loaderNonswap, SpectrumLoaderBase &loaderNue, SpectrumLoaderBase &loaderNuTau, SpectrumLoaderBase &loaderIntrinsic, const HistAxis &axis, SpillCut spillcut, Cut cut, const std::vector< const SystComponentScale * > &systs, const SystShifts &shift=kNoShift, const Var &wei=kUnweighted)
	This is for the FD via PredictionNoExtrap. More...

virtual	~PredictionScaleComp ()

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

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

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

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

Spectrum	PredictCategory (osc::IOscCalc osc, const SystComponentScale syst) const

virtual void	SaveTo (TDirectory *dir) const override

Public Member Functions inherited from ana::IPrediction
virtual	~IPrediction ()

virtual Spectrum	PredictUnoscillated () const

virtual OscillatableSpectrum	ComponentCC (int from, int to) const

Static Public Member Functions
static std::unique_ptr < PredictionScaleComp >	LoadFrom (TDirectory *dir)

Protected Member Functions
	PredictionScaleComp (const IPrediction total, const std::vector< const IPrediction > &preds, const std::vector< const SystComponentScale * > &systs)

Protected Attributes
std::vector< const SystComponentScale * >	fSysts

std::vector< const IPrediction * >	fPreds

const IPrediction *	fTotal

Detailed Description

Prediction broken down into arbitrary components whose scales can be varied independently.

Definition at line 14 of file PredictionScaleComp.h.

Constructor & Destructor Documentation

ana::PredictionScaleComp::PredictionScaleComp	(	SpectrumLoaderBase &	loader,
		const HistAxis &	axis,
		SpillCut	spillcut,
		Cut	cut,
		const std::vector< const SystComponentScale * > &	systs,
		const SystShifts &	shift = `kNoShift`,
		const Var &	wei = `kUnweighted`
	)

Parameters

cut	Cut applied to all histograms
truthcuts	Prediction will be broken down into N components following these cuts.

Definition at line 20 of file PredictionScaleComp.cxx.

     : fSysts(systs)
   {
     Cut complementCut = kNoCut;
 
     assert(!systs.empty() && "Please give at least one systematic.");
     for(const SystComponentScale* syst: systs){
       fPreds.push_back(new PredictionNoOsc(loader, axis,
                                            spillcut, cut && syst->GetCut(), shift, wei));
     }
 
     fTotal = new PredictionNoOsc(loader, axis, spillcut, cut, shift, wei);
   }

ana::PredictionScaleComp::PredictionScaleComp	(	SpectrumLoaderBase &	loader,
		const HistAxis &	axis1,
		const HistAxis &	axis2,
		SpillCut	spillcut,
		Cut	cut,
		const std::vector< const SystComponentScale * > &	systs,
		const SystShifts &	shift = `kNoShift`,
		const Var &	wei = `kUnweighted`
	)

Constructor to take two HistAxis's to weight 2D spectra.

Definition at line 66 of file PredictionScaleComp.cxx.

   {
     assert(0 && "unimplemented");
 
     // TODO TODO TODO
     /*
     assert(truthcuts.size()>0 && "Please give at least one truth selection.");
     for(unsigned int i = 0; i < truthcuts.size(); ++i){
       fPreds.push_back(new PredictionNoOsc(loader, axis1, axis2,
                                            cut && truthcuts[i], shift, wei));
       fSysts.push_back(new DummyScaleCompSyst(i));
       fComplementCut = fComplementCut && !truthcuts[i];
     }
 
     // The idea is that if truthcuts are exhaustive, this Spectrum should wind
     // up empty
     fComplement = new Spectrum(loader, axis1, axis2,
                                cut && fComplementCut, shift, wei);
     */
   }

ana::PredictionScaleComp::PredictionScaleComp	(	SpectrumLoaderBase &	loaderNonswap,
		SpectrumLoaderBase &	loaderNue,
		SpectrumLoaderBase &	loaderNuTau,
		SpectrumLoaderBase &	loaderIntrinsic,
		const HistAxis &	axis,
		SpillCut	spillcut,
		Cut	cut,
		const std::vector< const SystComponentScale * > &	systs,
		const SystShifts &	shift = `kNoShift`,
		const Var &	wei = `kUnweighted`
	)

This is for the FD via PredictionNoExtrap.

Definition at line 42 of file PredictionScaleComp.cxx.

     : fSysts(systs)
   {
     assert(!systs.empty() && "Please give at least one systematic.");
     for(const SystComponentScale* syst: systs){
       fPreds.push_back(new PredictionNoExtrap(loaderNonswap, loaderNue, loaderNuTau, loaderIntrinsic,
                                               axis, spillcut, cut && syst->GetCut(), shift, wei));
     }
 
     fTotal = new PredictionNoExtrap(loaderNonswap, loaderNue, loaderNuTau, loaderIntrinsic,
                                     axis, spillcut, cut, shift, wei);
   }

ana::PredictionScaleComp::~PredictionScaleComp ( )

virtual

Definition at line 106 of file PredictionScaleComp.cxx.

   {
     for(const IPrediction* p: fPreds) delete p;
     delete fTotal;
   }

ana::PredictionScaleComp::PredictionScaleComp	(	const IPrediction *	total,
		const std::vector< const IPrediction * > &	preds,
		const std::vector< const SystComponentScale * > &	systs
	)

protected

Definition at line 96 of file PredictionScaleComp.cxx.

     : fSysts(systs),
       fPreds(preds),
       fTotal(total)
   {
   }

Member Function Documentation

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

static

Definition at line 176 of file PredictionScaleComp.cxx.

   {
     IPrediction* total = ana::LoadFrom<IPrediction>(dir->GetDirectory("total")).release();
 
     std::vector<const IPrediction*> preds;
     for(unsigned int i = 0; ; ++i){
       TDirectory* di = dir->GetDirectory(("pred"+std::to_string(i)).c_str());
       if(!di) break; // We got all the predictions
 
       preds.push_back(ana::LoadFrom<IPrediction>(di).release());
     }
 
     std::vector<const SystComponentScale*> systs;
     for(unsigned int i = 0; ; ++i){
       TDirectory* si = dir->GetDirectory(("syst"+std::to_string(i)).c_str());
       if(!si) break; // We got all the predictions
 
       systs.push_back(ana::LoadFrom<SystComponentScale>(si).release());
     }
 
     return std::unique_ptr<PredictionScaleComp>(new PredictionScaleComp(total, preds, systs));
   }

virtual Spectrum ana::PredictionScaleComp::Predict ( osc::IOscCalc * osc ) const

inlineoverridevirtual

Implements ana::IPrediction.

Definition at line 51 of file PredictionScaleComp.h.

     {
       return fTotal->Predict(osc);
     }

Spectrum ana::PredictionScaleComp::PredictCategory	(	osc::IOscCalc *	osc,
		const SystComponentScale *	syst
	)		const

Definition at line 143 of file PredictionScaleComp.cxx.

   {
     for(unsigned int i = 0; i < fSysts.size(); ++i){
       if(fSysts[i] == syst) return fPreds[i]->Predict(osc);
     }
 
     std::cout << "PredictionScaleComp::PredictCategory(): Unknown systematic " << syst->ShortName() << std::endl;
     abort();
   }

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

inlineoverridevirtual

Implements ana::IPrediction.

Definition at line 63 of file PredictionScaleComp.h.

     {
       return fTotal->PredictComponent(calc, flav, curr, sign);
     }

Spectrum ana::PredictionScaleComp::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 113 of file PredictionScaleComp.cxx.

   {
     SystShifts shiftClean = shift;
     for(const ISyst* s: fSysts) shiftClean.SetShift(s, 0);
 
     // Starting with the total and adding in the differences for components
     // that need to change is faster if most of the knobs are set to zero.
     Spectrum ret = fTotal->PredictSyst(calc, shiftClean);
 
     for(unsigned int i = 0; i < fPreds.size(); ++i){
       const double x = shift.GetShift(fSysts[i]);
       if(x == 0) continue; // Nominal, can skip
 
       Spectrum si = fPreds[i]->PredictComponentSyst(calc, shiftClean,
                                                     flav, curr, sign);
 
       si.Scale(pow(1+fSysts[i]->OneSigmaScale(), x) - 1);
 
       ret += si;
     }
 
     return ret;
   }

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

inlineoverridevirtual

Reimplemented from ana::IPrediction.

Definition at line 56 of file PredictionScaleComp.h.

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

void ana::PredictionScaleComp::SaveTo ( TDirectory * dir ) const

overridevirtual

Reimplemented from ana::IPrediction.

Definition at line 155 of file PredictionScaleComp.cxx.

   {
     TDirectory* tmp = gDirectory;
     dir->cd();
 
     TObjString("PredictionScaleComp").Write("type");
 
     fTotal->SaveTo(dir->mkdir("total"));
 
     for(unsigned int i = 0; i < fPreds.size(); ++i){
       fPreds[i]->SaveTo(dir->mkdir(("pred"+std::to_string(i)).c_str()));
     }
 
     for(unsigned int i = 0; i < fSysts.size(); ++i){
       fSysts[i]->SaveTo(dir->mkdir(("syst"+std::to_string(i)).c_str()));
     }
 
     tmp->cd();
   }

Member Data Documentation

std::vector<const IPrediction*> ana::PredictionScaleComp::fPreds

protected

Definition at line 89 of file PredictionScaleComp.h.

std::vector<const SystComponentScale*> ana::PredictionScaleComp::fSysts

protected

Definition at line 88 of file PredictionScaleComp.h.

const IPrediction* ana::PredictionScaleComp::fTotal

protected

Definition at line 91 of file PredictionScaleComp.h.

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

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation