ana::PredictionExtendToPeripheral Class Reference

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-05-06/3FlavorAna/Prediction/PredictionExtendToPeripheral.h"

Inheritance diagram for ana::PredictionExtendToPeripheral:

Public Member Functions
	PredictionExtendToPeripheral (IPrediction *predCore, PredictionNoExtrap *predNoExtrap, bool mergePeripheral=true, int nbins=3)
	Takes ownership of predCore and predNoExtrap. More...
	PredictionExtendToPeripheral (IPrediction *predCore, PredictionNoExtrap *predNoExtrap, const ana::Binning &bins, bool mergePeripheral=true, int nbins=3)
virtual	~PredictionExtendToPeripheral ()
virtual Spectrum	Predict (osc::IOscCalc *calc) const override
virtual Spectrum	Predict (osc::IOscCalcStan *calc) const override
virtual Spectrum	PredictComponent (osc::IOscCalc *calc, Flavors::Flavors_t flav, Current::Current_t curr, Sign::Sign_t sign) const override
virtual 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
void	Reduce ()
virtual void	SaveTo (TDirectory *dir, const std::string &name) const override
	PredictionExtendToPeripheral ()=delete
IPrediction *	GetCore ()
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< PredictionExtendToPeripheral >	LoadFrom (TDirectory *dir, const std::string &name)

Protected Member Functions
OscillatableSpectrum	ReduceHelper (const OscillatableSpectrum &s) const
template<typename T >
Spectrum	_PredictComponent (osc::_IOscCalc< T > *calc, Flavors::Flavors_t flav, Current::Current_t curr, Sign::Sign_t sign) const
Spectrum	ReduceHelperNC (const Spectrum &s) const
Spectrum	MergePeripheral (const Spectrum &s) const
OscillatableSpectrum	MergePeripheralOsc (const OscillatableSpectrum &s) const
Spectrum	ExtendRecoWeight (const Spectrum &Core, const Spectrum &NoExtrap) const
OscillatableSpectrum	ExtendRecoWeightOscillatable (const OscillatableSpectrum &Core, OscillatableSpectrum NoExtrap) const

Protected Attributes
IPrediction *	fPredCore
PredictionNoExtrap *	fPredNoExtrap
Binning	fBins
bool	fMergePeripheral
int	fNPIDBin
const int	fNBinsI = fBins.NBins()
const int	fEnergyBins = fNBinsI/fNPIDBin
const int	fPeriNum = int(fEnergyBins/2)
const int	fNBinsO = fNBinsI-fPeriNum
const int	fNCore = (fNPIDBin-1)*fEnergyBins
const int	fPeriBin = fNCore+int((fEnergyBins-fPeriNum)/2)+1
const int	fReduceBin = fNBinsI - 2* fNPIDBin - fEnergyBins +1
const Binning	ReduceBinning = Binning::Simple(fReduceBin, 0, fReduceBin)
OscillatableSpectrum	fNuEToNuE
OscillatableSpectrum	fAntiNuEToAntiNuE
OscillatableSpectrum	fNuEToNuMu
OscillatableSpectrum	fAntiNuEToAntiNuMu
OscillatableSpectrum	fNuEToNuTau
OscillatableSpectrum	fAntiNuEToAntiNuTau
OscillatableSpectrum	fNuMuToNuE
OscillatableSpectrum	fAntiNuMuToAntiNuE
OscillatableSpectrum	fNuMuToNuMu
OscillatableSpectrum	fAntiNuMuToAntiNuMu
OscillatableSpectrum	fNuMuToNuTau
OscillatableSpectrum	fAntiNuMuToAntiNuTau
Spectrum	fNCTot
Spectrum	fNC
Spectrum	fNCAnti
bool	fIsMerged

Detailed Description

Combines two extrapolations, one with extrapolated predictions of the core bins (and nominal peripheral) and one with all nominal predictions, to apply the correct weights to the peripheral sample (copying the high PID bin).

Definition at line 14 of file PredictionExtendToPeripheral.h.

Constructor & Destructor Documentation

ana::PredictionExtendToPeripheral::PredictionExtendToPeripheral	(	IPrediction *	predCore,
		PredictionNoExtrap *	predNoExtrap,
		bool	mergePeripheral = true,
		int	nbins = 3
	)

Takes ownership of predCore and predNoExtrap.

Definition at line 24 of file PredictionExtendToPeripheral.cxx.

    : PredictionExtendToPeripheral(predCore, predNoExtrap, Binning::Simple(27,0,27), mergePeripheral, nbins)
  {
  }

ana::PredictionExtendToPeripheral::PredictionExtendToPeripheral	(	IPrediction *	predCore,
		PredictionNoExtrap *	predNoExtrap,
		const ana::Binning &	bins,
		bool	mergePeripheral = true,
		int	nbins = 3
	)

Definition at line 30 of file PredictionExtendToPeripheral.cxx.

References om::cout, allTimeWatchdog::endl, and fMergePeripheral.

    : fPredCore(predCore),
      fPredNoExtrap(predNoExtrap),
      fBins(bins),
      fMergePeripheral(mergePeripheral),
      fNPIDBin(nbins),
      fNuEToNuE(ExtendRecoWeightOscillatable(fPredCore->ComponentCC(12, 12),fPredNoExtrap->ComponentCC(12, 12))),
      fAntiNuEToAntiNuE(ExtendRecoWeightOscillatable(fPredCore->ComponentCC(-12, -12),fPredNoExtrap->ComponentCC(-12, -12))),
      fNuEToNuMu(ExtendRecoWeightOscillatable(fPredCore->ComponentCC(12, 14),fPredNoExtrap->ComponentCC(12, 14))),
      fAntiNuEToAntiNuMu(ExtendRecoWeightOscillatable(fPredCore->ComponentCC(-12, -14),fPredNoExtrap->ComponentCC(-12, -14))),
      fNuEToNuTau(ExtendRecoWeightOscillatable(fPredCore->ComponentCC(12, 16),fPredNoExtrap->ComponentCC(12, 16))),
      fAntiNuEToAntiNuTau(ExtendRecoWeightOscillatable(fPredCore->ComponentCC(-12, -16),fPredNoExtrap->ComponentCC(-12, -16))),
      fNuMuToNuE(ExtendRecoWeightOscillatable(fPredCore->ComponentCC(14, 12),fPredNoExtrap->ComponentCC(14, 12))),
      fAntiNuMuToAntiNuE(ExtendRecoWeightOscillatable(fPredCore->ComponentCC(-14, -12),fPredNoExtrap->ComponentCC(-14, -12))),
      fNuMuToNuMu(ExtendRecoWeightOscillatable(fPredCore->ComponentCC(14, 14),fPredNoExtrap->ComponentCC(14, 14))),
      fAntiNuMuToAntiNuMu(ExtendRecoWeightOscillatable(fPredCore->ComponentCC(-14, -14),fPredNoExtrap->ComponentCC(-14, -14))),
      fNuMuToNuTau(ExtendRecoWeightOscillatable(fPredCore->ComponentCC(14, 16),fPredNoExtrap->ComponentCC(14, 16))),
      fAntiNuMuToAntiNuTau(ExtendRecoWeightOscillatable(fPredCore->ComponentCC(-14, -16),fPredNoExtrap->ComponentCC(-14, -16))),
      fNCTot(ExtendRecoWeight(fPredCore->ComponentNCTotal(), fPredNoExtrap->ComponentNCTotal())),
      fNC(ExtendRecoWeight(fPredCore->ComponentNC(), fPredNoExtrap->ComponentNC())),
      fNCAnti(ExtendRecoWeight(fPredCore->ComponentNCAnti(), fPredNoExtrap->ComponentNCAnti())),
      fIsMerged(false)
  {

    if(!fMergePeripheral) std::cout << "WARNING: Not merging peripheral. Check inputs." << std::endl;
  }

ana::PredictionExtendToPeripheral::~PredictionExtendToPeripheral ( )

virtual

Definition at line 58 of file PredictionExtendToPeripheral.cxx.

  {
    // We don't necessarily own the predictions, leave them alone
  }

ana::PredictionExtendToPeripheral::PredictionExtendToPeripheral ( )

delete

Referenced by LoadFrom().

Member Function Documentation

template<typename T >

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

protected

Definition at line 298 of file PredictionExtendToPeripheral.cxx.

References ana::assert(), calc, ana::Spectrum::Clear(), ComponentNC(), ComponentNCAnti(), febshutoff_auto::curr, fAntiNuEToAntiNuE, fAntiNuEToAntiNuMu, fAntiNuEToAntiNuTau, fAntiNuMuToAntiNuE, fAntiNuMuToAntiNuMu, fAntiNuMuToAntiNuTau, fNCTot, fNuEToNuE, fNuEToNuMu, fNuEToNuTau, fNuMuToNuE, fNuMuToNuMu, fNuMuToNuTau, 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, ana::OscillatableSpectrum::Oscillated(), runNovaSAM::ret, and canMan::sign().

Referenced by GetCore(), and PredictComponent().

  {
    Spectrum ret = fNCTot;

    ret.Clear();

    if(curr & Current::kCC){

      if(flav & Flavors::kNuEToNuE   && sign & Sign::kNu)
        ret += fNuEToNuE.Oscillated(calc, +12, +12);
      if(flav & Flavors::kNuEToNuE   && sign & Sign::kAntiNu)
        ret += fAntiNuEToAntiNuE.Oscillated(calc, -12, -12);
      if(flav & Flavors::kNuEToNuMu  && sign & Sign::kNu)
        ret += fNuEToNuMu.Oscillated(calc, +12, +14);
      if(flav & Flavors::kNuEToNuMu  && sign & Sign::kAntiNu)
        ret += fAntiNuEToAntiNuMu.Oscillated(calc, -12, -14);
      if(flav & Flavors::kNuEToNuTau && sign & Sign::kNu)
        ret += fNuEToNuTau.Oscillated(calc, +12, +16);
      if(flav & Flavors::kNuEToNuTau && sign & Sign::kAntiNu)
        ret += fAntiNuEToAntiNuTau.Oscillated(calc, -12, -16);

      if(flav & Flavors::kNuMuToNuE   && sign & Sign::kNu)
        ret += fNuMuToNuE.Oscillated(calc, +14, +12);
      if(flav & Flavors::kNuMuToNuE   && sign & Sign::kAntiNu)
        ret += fAntiNuMuToAntiNuE.Oscillated(calc, -14, -12);
      if(flav & Flavors::kNuMuToNuMu  && sign & Sign::kNu)
        ret += fNuMuToNuMu.Oscillated(calc, +14, +14);
      if(flav & Flavors::kNuMuToNuMu  && sign & Sign::kAntiNu)
        ret += fAntiNuMuToAntiNuMu.Oscillated(calc, -14, -14);
      if(flav & Flavors::kNuMuToNuTau && sign & Sign::kNu)
        ret += fNuMuToNuTau.Oscillated(calc, +14, +16);
      if(flav & Flavors::kNuMuToNuTau && sign & Sign::kAntiNu)
        ret += fAntiNuMuToAntiNuTau.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 += this->ComponentNC();
     if(sign & Sign::kAntiNu) ret += this->ComponentNCAnti();

    }

    return ret;
  }

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

overridevirtual

Reimplemented from ana::IPrediction.

Definition at line 381 of file PredictionExtendToPeripheral.cxx.

References ana::assert(), fAntiNuEToAntiNuE, fAntiNuEToAntiNuMu, fAntiNuEToAntiNuTau, fAntiNuMuToAntiNuE, fAntiNuMuToAntiNuMu, fAntiNuMuToAntiNuTau, fNuEToNuE, fNuEToNuMu, fNuEToNuTau, fNuMuToNuE, fNuMuToNuMu, and fNuMuToNuTau.

  {
    // beam nues and numu survival come from extending the effect on the high
    // CVN bin into the peripheral sample
    if (from == +12 && to == +12) return fNuEToNuE;
    if (from == +12 && to == +14) return fNuEToNuMu;
    if (from == +12 && to == +16) return fNuEToNuTau;

    if (from == +14 && to == +12) return fNuMuToNuE;
    if (from == +14 && to == +14) return fNuMuToNuMu;
    if (from == +14 && to == +16) return fNuMuToNuTau;

    if (from == -12 && to == -12) return fAntiNuEToAntiNuE;
    if (from == -12 && to == -14) return fAntiNuEToAntiNuMu;
    if (from == -12 && to == -16) return fAntiNuEToAntiNuTau;

    if (from == -14 && to == -12) return fAntiNuMuToAntiNuE;
    if (from == -14 && to == -14) return fAntiNuMuToAntiNuMu;
    if (from == -14 && to == -16) return fAntiNuMuToAntiNuTau;

    assert(0 && "Not reached");
  }

Spectrum ana::PredictionExtendToPeripheral::ComponentNC ( ) const

overridevirtual

Reimplemented from ana::IPrediction.

Definition at line 410 of file PredictionExtendToPeripheral.cxx.

References fNC.

Referenced by _PredictComponent().

  {
    return fNC;
  }

Spectrum ana::PredictionExtendToPeripheral::ComponentNCAnti ( ) const

overridevirtual

Reimplemented from ana::IPrediction.

Definition at line 414 of file PredictionExtendToPeripheral.cxx.

References fNCAnti.

Referenced by _PredictComponent().

  {
    return fNCAnti;
  }

Spectrum ana::PredictionExtendToPeripheral::ComponentNCTotal ( ) const

overridevirtual

Reimplemented from ana::IPrediction.

Definition at line 406 of file PredictionExtendToPeripheral.cxx.

References fNCTot.

  {
    return fNCTot;
  }

Spectrum ana::PredictionExtendToPeripheral::ExtendRecoWeight ( const Spectrum &  Core, const Spectrum &  NoExtrap  ) const

protected

Definition at line 232 of file PredictionExtendToPeripheral.cxx.

References bin, fEnergyBins, fMergePeripheral, fNBinsI, fNCore, ana::Spectrum::GetBinnings(), ana::Ratio::GetEigen(), ana::Spectrum::GetLabels(), and MergePeripheral().

Referenced by GetCore().

  {
    // NoExtrap has the nominal FD prediction in it, while Core has the
    // our decomposed prediction.  We want to copy over the decomp weights,
    // which will just be the Core / NoExtrap bin content.  That covers
    // extrapolation in energy bins.
    Eigen::ArrayXd arat = Ratio(Core,NoExtrap).GetEigen();

    const int FirstPeripheralBin = fNCore+1;
    const int LastPeripheralBin  = fNBinsI;
    const int EBinsPerSelBin = fEnergyBins;

    for (int bin = FirstPeripheralBin; bin<= LastPeripheralBin; bin++){
      // Copy decomp ratio from high CVN bin to peripheral bin
      arat[bin] = arat[bin-EBinsPerSelBin];
    }
    Ratio rat(std::move(arat), Core.GetLabels(), Core.GetBinnings());

    Spectrum extrap = NoExtrap*rat;

    // Done if we don't want to collapse energy spectrum in peripheral sample
    if (!fMergePeripheral) return extrap;

    // Now, we have to merge all energy bins in the peripheral sample
    // This is only applicable to the Nue2017 result, so assume 28 bins
    return MergePeripheral(extrap);
  }

OscillatableSpectrum ana::PredictionExtendToPeripheral::ExtendRecoWeightOscillatable ( const OscillatableSpectrum &  Core, OscillatableSpectrum  NoExtrap  ) const

protected

Definition at line 263 of file PredictionExtendToPeripheral.cxx.

References bin, fEnergyBins, fMergePeripheral, fNBinsO, fNCore, ana::Spectrum::GetBinnings(), ana::Ratio::GetEigen(), ana::Spectrum::GetLabels(), MergePeripheralOsc(), NoExtrap(), ana::ReweightableSpectrum::ReweightToRecoSpectrum(), and ana::OscillatableSpectrum::Unoscillated().

Referenced by GetCore().

  {
    // NoExtrap has the nominal FD prediction in it, while Core has the
    // our decomposed prediction.  We want to copy over the decomp weights,
    // which will just be the Core / NoExtrap bin content

    // Calling Unoscillated() will project down to the predicted reco spectra
   
    Spectrum recoCore = Core.Unoscillated();

    Spectrum recoNoExtrap = NoExtrap.Unoscillated();
    Eigen::ArrayXd arat = Ratio(recoCore,recoNoExtrap).GetEigen();

    const int FirstPeripheralBin = fNCore+1;
    const int LastPeripheralBin  = fNBinsO;

    const int EBinsPerSelBin = fEnergyBins;
    for (int bin = FirstPeripheralBin; bin<= LastPeripheralBin; bin++){
      // Copy decomp ratio from high CVN bin to peripheral bin
      arat[bin] = arat[bin-EBinsPerSelBin];
    }
    Ratio rat(std::move(arat), recoCore.GetLabels(), recoCore.GetBinnings());
    NoExtrap.ReweightToRecoSpectrum(rat*recoNoExtrap);

    // We're done if we don't want to merge peripheral sample energy dist
    if (!fMergePeripheral) return NoExtrap;

    return MergePeripheralOsc(NoExtrap);
  }

IPrediction* ana::PredictionExtendToPeripheral::GetCore ( )

inline

Definition at line 65 of file PredictionExtendToPeripheral.h.

References _PredictComponent(), ExtendRecoWeight(), ExtendRecoWeightOscillatable(), fPredCore, MergePeripheral(), MergePeripheralOsc(), NoExtrap(), ReduceHelper(), and ReduceHelperNC().

    {
      return fPredCore;
    };

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

static

Definition at line 445 of file PredictionExtendToPeripheral.cxx.

References ana::assert(), ana::bins, dir, ana::Binning::LoadFrom(), ana::LoadFrom< IPrediction >(), PredictionExtendToPeripheral(), and runNovaSAM::release.

  {
    dir = dir->GetDirectory(name.c_str()); // switch to subdir
    assert(dir);

    auto predCore = ana::LoadFrom<IPrediction>(dir, "predCore").release();
    PredictionNoExtrap* predNoExtrap = ana::LoadFrom<PredictionNoExtrap>(dir, "predNoExtrap").release();

    Binning bins= *Binning::LoadFrom(dir, "bins");

    bool mergeperipheral = true; // default

    TVectorD* vMergePeripheral = (TVectorD*)dir->Get("mergePeripheral");
    // An attempt to be backwards compatable to already generated files
    if (vMergePeripheral){
      assert(vMergePeripheral->GetNrows() == 1);
      mergeperipheral = bool((*vMergePeripheral)[0]);
    }

    delete dir;

    return std::unique_ptr<PredictionExtendToPeripheral>(new PredictionExtendToPeripheral(predCore, predNoExtrap, bins, mergeperipheral));
  }

Spectrum ana::PredictionExtendToPeripheral::MergePeripheral ( const Spectrum &  s ) const

protected

Definition at line 187 of file PredictionExtendToPeripheral.cxx.

References a, fNBinsI, fNBinsO, fNCore, fPeriBin, ana::Spectrum::GetEigen(), ana::Spectrum::GetLabels(), MECModelEnuComparisons::i, ana::Spectrum::Livetime(), ana::Spectrum::POT(), and ana::Binning::Simple().

Referenced by ExtendRecoWeight(), and GetCore().

  {
    Eigen::ArrayXd a = s.GetEigen(s.POT());
    Eigen::ArrayXd aret(fNBinsO+2);
    aret.setZero();
    for (int i = 1; i <= fNCore; i++){
      aret[i] = a[i]; // Copy Core
    }

    // Merge peripheral
    for(int i = 0; i <= 8; ++i) aret[fPeriBin] += a[fNBinsI-i];

    return Spectrum(std::move(aret),
                    HistAxis(s.GetLabels(), {Binning::Simple(fNBinsO, 0, fNBinsO)}),
                    s.POT(),
                    s.Livetime());
  }

OscillatableSpectrum ana::PredictionExtendToPeripheral::MergePeripheralOsc ( const OscillatableSpectrum &  s ) const

protected

Definition at line 206 of file PredictionExtendToPeripheral.cxx.

References fNBinsI, fNBinsO, fNCore, fPeriBin, ana::ReweightableSpectrum::GetEigen(), ana::ReweightableSpectrum::GetLabels(), ana::kTrueEnergyBins, ana::ReweightableSpectrum::Livetime(), ana::Binning::NBins(), ana::ReweightableSpectrum::POT(), and ana::Binning::Simple().

Referenced by ExtendRecoWeightOscillatable(), and GetCore().

  {
    Eigen::MatrixXd aret(kTrueEnergyBins.NBins()+2, fNBinsO+2);
    aret.setZero();
    Eigen::MatrixXd aextrap = s.GetEigen(s.POT());

    for (int binX = 1; binX <= fNCore; binX++){
      for (int binY = 1; binY <= kTrueEnergyBins.NBins(); binY++){
        aret(binY, binX) = aextrap(binY, binX);
      }
    }
    // Now, project all the bins in columns 28+ down to the 28th
    for (int binX = fNCore+1; binX <= fNBinsI; binX++){
      for (int binY = 1; binY <= kTrueEnergyBins.NBins(); binY++){
        aret(binY, fPeriBin) += aextrap(binY, binX);
      }
    }

    return OscillatableSpectrum(std::move(aret),
                                HistAxis(s.GetLabels(), {Binning::Simple(fNBinsO, 0, fNBinsO)}),
                                s.POT(),
                                s.Livetime());
  }

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

overridevirtual

Implements ana::IPrediction.

Definition at line 169 of file PredictionExtendToPeripheral.cxx.

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

Referenced by test_nue2017Prediction().

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

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

overridevirtual

Reimplemented from ana::IPrediction.

Definition at line 178 of file PredictionExtendToPeripheral.cxx.

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

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

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

overridevirtual

Implements ana::IPrediction.

Definition at line 361 of file PredictionExtendToPeripheral.cxx.

References _PredictComponent().

Referenced by Predict(), and test_nue2017Prediction().

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

Spectrum ana::PredictionExtendToPeripheral::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 371 of file PredictionExtendToPeripheral.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 79 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 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::IPrediction::PredictSyst ( osc::IOscCalc *  calc, const SystShifts &  syst  ) const

virtualinherited

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

Definition at line 49 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(), GeneratePseudoexpt(), joint_fit_2017_make_fc_slice(), joint_fit_2017_make_fc_surf(), MakeFakeData(), ana::PlotWithAreaSystErrorBand(), ana::PlotWithSystErrorBand(), ana::PlotWithSystErrorBand_Quant(), ana::SingleSampleExperiment::PredHistIncCosmics(), ana::NDPredictionNoOsc::~NDPredictionNoOsc(), 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 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(), template_basic(), and test_ana().

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

void ana::PredictionExtendToPeripheral::Reduce

(

)

Definition at line 146 of file PredictionExtendToPeripheral.cxx.

References ana::assert(), fAntiNuEToAntiNuE, fAntiNuEToAntiNuMu, fAntiNuEToAntiNuTau, fAntiNuMuToAntiNuE, fAntiNuMuToAntiNuMu, fAntiNuMuToAntiNuTau, fIsMerged, fNC, fNCAnti, fNCTot, fNuEToNuE, fNuEToNuMu, fNuEToNuTau, fNuMuToNuE, fNuMuToNuMu, fNuMuToNuTau, ReduceHelper(), and ReduceHelperNC().

  {
    assert(!fIsMerged && "Can't call Reduce() twice");

    fNuEToNuE = ReduceHelper(fNuEToNuE); 
    fAntiNuEToAntiNuE = ReduceHelper(fAntiNuEToAntiNuE); 
    fNuEToNuMu = ReduceHelper(fNuEToNuMu); 
    fAntiNuEToAntiNuMu = ReduceHelper(fAntiNuEToAntiNuMu);
    fNuEToNuTau = ReduceHelper(fNuEToNuTau); 
    fAntiNuEToAntiNuTau = ReduceHelper(fAntiNuEToAntiNuTau);
    fNuMuToNuE = ReduceHelper(fNuMuToNuE); 
    fAntiNuMuToAntiNuE = ReduceHelper(fAntiNuMuToAntiNuE); 
    fNuMuToNuMu = ReduceHelper(fNuMuToNuMu); 
    fAntiNuMuToAntiNuMu = ReduceHelper(fAntiNuMuToAntiNuMu);
    fNuMuToNuTau = ReduceHelper(fNuMuToNuTau); 
    fAntiNuMuToAntiNuTau = ReduceHelper(fAntiNuMuToAntiNuTau);
 
    fNCTot = ReduceHelperNC(fNCTot); fNC = ReduceHelperNC(fNC); fNCAnti = ReduceHelperNC(fNCAnti);

    fIsMerged = true; // dont collapse the whole universe...
  }

OscillatableSpectrum ana::PredictionExtendToPeripheral::ReduceHelper ( const OscillatableSpectrum &  s ) const

protected

Definition at line 65 of file PredictionExtendToPeripheral.cxx.

References ana::assert(), ana::bins, fNBinsO, fNCore, fNPIDBin, ana::ReweightableSpectrum::GetBinnings(), ana::ReweightableSpectrum::GetEigen(), ana::ReweightableSpectrum::GetLabels(), MECModelEnuComparisons::i, calib::j, ana::kTrueEnergyBins, PandAna.Demos.pi0_spectra::labels, livetime, ana::ReweightableSpectrum::Livetime(), ana::Binning::NBins(), pot, ana::ReweightableSpectrum::POT(), ReduceBinning, and submit_syst::x.

Referenced by GetCore(), and Reduce().

  {
    std::vector<std::string> labels = s.GetLabels();
    std::vector<Binning> bins = s.GetBinnings();
    double pot = s.POT();
    double livetime = s.Livetime();

    std::vector<int> Bins;

    for(int i = 1; i < fNBinsO+1; i++){ Bins.push_back(i); }
    int iter = fNPIDBin;

    Bins.erase(Bins.end()-2, Bins.end()); // bye peripheral
    Bins.erase(Bins.begin()+(fNCore),Bins.begin()+(fNCore+2));

    while(iter-1){ // peripheral is done
      for(int j: {9, 2, 1}){
        Bins.erase(Bins.begin() +((iter-2)* 9+j-1));
      }
      --iter;
    }  

    Eigen::MatrixXd ain = s.GetEigen(s.POT());
    assert(s.GetBinnings()[0].NBins() == fNBinsO && "Merge Peripheral Before Reduce!");
    Eigen::MatrixXd aout(kTrueEnergyBins.NBins()+2, ReduceBinning.NBins()+2);
    aout.setZero();

    int i = 0;
    for(int x: Bins){
      ++i;
      for(int j = 0; j < ain.size(); ++j){
        aout(j, i) = ain(j, x);
      }
    }

    return OscillatableSpectrum(std::move(aout),
                                HistAxis(labels, bins),
                                pot, livetime);

  }

Spectrum ana::PredictionExtendToPeripheral::ReduceHelperNC ( const Spectrum &  s ) const

protected

Definition at line 108 of file PredictionExtendToPeripheral.cxx.

References ana::assert(), ana::bins, fNBinsO, fNCore, fNPIDBin, fReduceBin, ana::Spectrum::GetBinnings(), ana::Spectrum::GetEigen(), ana::Spectrum::GetLabels(), MECModelEnuComparisons::i, calib::j, PandAna.Demos.pi0_spectra::labels, ana::Spectrum::Livetime(), ana::Binning::NBins(), pot, ana::Spectrum::POT(), ReduceBinning, and submit_syst::x.

Referenced by GetCore(), and Reduce().

  {
    std::vector<std::string> labels = s.GetLabels();
    std::vector<Binning> bins = s.GetBinnings();
    double pot = s.POT();

    std::vector<int> Bins;

    for(int i = 1; i < fNBinsO+1; i++){ Bins.push_back(i); }
    int iter = fNPIDBin;

    Bins.erase(Bins.end()-2, Bins.end()); // bye peripheral
    Bins.erase(Bins.begin()+(fNCore),Bins.begin()+fNCore+2);

    while(iter-1){ // peripheral is done
      for(int j: {9, 2, 1}){
        Bins.erase(Bins.begin() + ((iter-2)* 9+j-1));
      }
      --iter;
    }

    Eigen::ArrayXd ain = s.GetEigen(pot);
    assert( ain.size()-2 == fNBinsO && "Merge Peripheral Before Reduce!");
    Eigen::ArrayXd aout(ReduceBinning.NBins()+2);
    aout.setZero();
    
    for(int x = 0; x < fReduceBin+2; ++x){
      aout[x]= ain[Bins[x-1]];
    }

    std::vector<Binning> tmp_bin = {ReduceBinning};
    return Spectrum(std::move(aout),
                    HistAxis(s.GetLabels(), tmp_bin),
                    s.POT(), s.Livetime());
  }  

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

overridevirtual

Reimplemented from ana::IPrediction.

Definition at line 420 of file PredictionExtendToPeripheral.cxx.

References dir, fBins, fMergePeripheral, fPredCore, fPredNoExtrap, ana::Binning::SaveTo(), ana::PredictionNoExtrap::SaveTo(), ana::IPrediction::SaveTo(), and tmp.

Referenced by make_nueFDprediction_kinematics_FHC(), make_nueFDprediction_kinematics_FHC_REW(), make_nueFDprediction_kinematics_RHC(), and make_nueFDprediction_kinematics_RHC_REW().

  {
    TDirectory* tmp = gDirectory;

    dir = dir->mkdir(name.c_str()); // switch to subdir
    dir->cd();

    TObjString("PredictionExtendToPeripheral").Write("type");

    fPredCore->SaveTo(dir, "predCore");
    fPredNoExtrap->SaveTo(dir, "predNoExtrap");

    fBins.SaveTo(dir, "bins");

    TVectorD vMergePeripheral(1);
    vMergePeripheral[0] = fMergePeripheral;
    vMergePeripheral.Write("mergePeripheral");

    dir->Write();
    delete dir;

    tmp->cd();
  }

Member Data Documentation

OscillatableSpectrum ana::PredictionExtendToPeripheral::fAntiNuEToAntiNuE

protected

Definition at line 112 of file PredictionExtendToPeripheral.h.

Referenced by _PredictComponent(), ComponentCC(), and Reduce().

OscillatableSpectrum ana::PredictionExtendToPeripheral::fAntiNuEToAntiNuMu

protected

Definition at line 114 of file PredictionExtendToPeripheral.h.

Referenced by _PredictComponent(), ComponentCC(), and Reduce().

OscillatableSpectrum ana::PredictionExtendToPeripheral::fAntiNuEToAntiNuTau

protected

Definition at line 116 of file PredictionExtendToPeripheral.h.

Referenced by _PredictComponent(), ComponentCC(), and Reduce().

OscillatableSpectrum ana::PredictionExtendToPeripheral::fAntiNuMuToAntiNuE

protected

Definition at line 119 of file PredictionExtendToPeripheral.h.

Referenced by _PredictComponent(), ComponentCC(), and Reduce().

OscillatableSpectrum ana::PredictionExtendToPeripheral::fAntiNuMuToAntiNuMu

protected

Definition at line 121 of file PredictionExtendToPeripheral.h.

Referenced by _PredictComponent(), ComponentCC(), and Reduce().

OscillatableSpectrum ana::PredictionExtendToPeripheral::fAntiNuMuToAntiNuTau

protected

Definition at line 123 of file PredictionExtendToPeripheral.h.

Referenced by _PredictComponent(), ComponentCC(), and Reduce().

Binning ana::PredictionExtendToPeripheral::fBins

protected

Definition at line 90 of file PredictionExtendToPeripheral.h.

Referenced by SaveTo().

const int ana::PredictionExtendToPeripheral::fEnergyBins = fNBinsI/fNPIDBin

protected

Definition at line 96 of file PredictionExtendToPeripheral.h.

Referenced by ExtendRecoWeight(), and ExtendRecoWeightOscillatable().

bool ana::PredictionExtendToPeripheral::fIsMerged

protected

Definition at line 130 of file PredictionExtendToPeripheral.h.

Referenced by Reduce().

bool ana::PredictionExtendToPeripheral::fMergePeripheral

protected

Definition at line 91 of file PredictionExtendToPeripheral.h.

Referenced by ExtendRecoWeight(), ExtendRecoWeightOscillatable(), PredictionExtendToPeripheral(), and SaveTo().

const int ana::PredictionExtendToPeripheral::fNBinsI = fBins.NBins()

protected

Definition at line 95 of file PredictionExtendToPeripheral.h.

Referenced by ExtendRecoWeight(), MergePeripheral(), and MergePeripheralOsc().

const int ana::PredictionExtendToPeripheral::fNBinsO = fNBinsI-fPeriNum

protected

Definition at line 98 of file PredictionExtendToPeripheral.h.

Referenced by ExtendRecoWeightOscillatable(), MergePeripheral(), MergePeripheralOsc(), ReduceHelper(), and ReduceHelperNC().

Spectrum ana::PredictionExtendToPeripheral::fNC

protected

Definition at line 126 of file PredictionExtendToPeripheral.h.

Referenced by ComponentNC(), and Reduce().

Spectrum ana::PredictionExtendToPeripheral::fNCAnti

protected

Definition at line 127 of file PredictionExtendToPeripheral.h.

Referenced by ComponentNCAnti(), and Reduce().

const int ana::PredictionExtendToPeripheral::fNCore = (fNPIDBin-1)*fEnergyBins

protected

Definition at line 100 of file PredictionExtendToPeripheral.h.

Referenced by ExtendRecoWeight(), ExtendRecoWeightOscillatable(), MergePeripheral(), MergePeripheralOsc(), ReduceHelper(), and ReduceHelperNC().

Spectrum ana::PredictionExtendToPeripheral::fNCTot

protected

Definition at line 125 of file PredictionExtendToPeripheral.h.

Referenced by _PredictComponent(), ComponentNCTotal(), and Reduce().

int ana::PredictionExtendToPeripheral::fNPIDBin

protected

Definition at line 92 of file PredictionExtendToPeripheral.h.

Referenced by ReduceHelper(), and ReduceHelperNC().

OscillatableSpectrum ana::PredictionExtendToPeripheral::fNuEToNuE

protected

Definition at line 111 of file PredictionExtendToPeripheral.h.

Referenced by _PredictComponent(), ComponentCC(), and Reduce().

OscillatableSpectrum ana::PredictionExtendToPeripheral::fNuEToNuMu

protected

Definition at line 113 of file PredictionExtendToPeripheral.h.

Referenced by _PredictComponent(), ComponentCC(), and Reduce().

OscillatableSpectrum ana::PredictionExtendToPeripheral::fNuEToNuTau

protected

Definition at line 115 of file PredictionExtendToPeripheral.h.

Referenced by _PredictComponent(), ComponentCC(), and Reduce().

OscillatableSpectrum ana::PredictionExtendToPeripheral::fNuMuToNuE

protected

Definition at line 118 of file PredictionExtendToPeripheral.h.

Referenced by _PredictComponent(), ComponentCC(), and Reduce().

OscillatableSpectrum ana::PredictionExtendToPeripheral::fNuMuToNuMu

protected

Definition at line 120 of file PredictionExtendToPeripheral.h.

Referenced by _PredictComponent(), ComponentCC(), and Reduce().

OscillatableSpectrum ana::PredictionExtendToPeripheral::fNuMuToNuTau

protected

Definition at line 122 of file PredictionExtendToPeripheral.h.

Referenced by _PredictComponent(), ComponentCC(), and Reduce().

const int ana::PredictionExtendToPeripheral::fPeriBin = fNCore+int((fEnergyBins-fPeriNum)/2)+1

protected

Definition at line 101 of file PredictionExtendToPeripheral.h.

Referenced by MergePeripheral(), and MergePeripheralOsc().

const int ana::PredictionExtendToPeripheral::fPeriNum = int(fEnergyBins/2)

protected

Definition at line 97 of file PredictionExtendToPeripheral.h.

IPrediction* ana::PredictionExtendToPeripheral::fPredCore

protected

Definition at line 88 of file PredictionExtendToPeripheral.h.

Referenced by GetCore(), SaveTo(), and ana::PredictExtendOwning::~PredictExtendOwning().

PredictionNoExtrap* ana::PredictionExtendToPeripheral::fPredNoExtrap

protected

Definition at line 89 of file PredictionExtendToPeripheral.h.

Referenced by SaveTo(), and ana::PredictExtendOwning::~PredictExtendOwning().

const int ana::PredictionExtendToPeripheral::fReduceBin = fNBinsI - 2* fNPIDBin - fEnergyBins +1

protected

Definition at line 108 of file PredictionExtendToPeripheral.h.

Referenced by ReduceHelperNC().

const Binning ana::PredictionExtendToPeripheral::ReduceBinning = Binning::Simple(fReduceBin, 0, fReduceBin)

protected

Definition at line 109 of file PredictionExtendToPeripheral.h.

Referenced by ReduceHelper(), and ReduceHelperNC().

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-05-06/3FlavorAna/Prediction/PredictionExtendToPeripheral.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-05-06/3FlavorAna/Prediction/PredictionExtendToPeripheral.cxx