ana::OscillatableSpectrum Class Reference

Spectrum with true energy information, allowing it to be oscillated More...

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-09-23/CAFAna/Core/OscillatableSpectrum.h"

Inheritance diagram for ana::OscillatableSpectrum:

Public Member Functions
	OscillatableSpectrum (const std::string &label, const Binning &bins, SpectrumLoaderBase &loader, const Var &var, const Cut &cut, const SystShifts &shift=kNoShift, const Var &wei=kUnweighted)
	OscillatableSpectrum (SpectrumLoaderBase &loader, const HistAxis &axis, const Cut &cut, const SystShifts &shift=kNoShift, const Var &wei=kUnweighted)
	OscillatableSpectrum (const Eigen::MatrixXd &&mat, const HistAxis &recoAxis, double pot, double livetime)
	~OscillatableSpectrum ()
	OscillatableSpectrum (const OscillatableSpectrum &rhs)
	Copy constructor. More...
	OscillatableSpectrum (OscillatableSpectrum &&rhs)
OscillatableSpectrum &	operator= (const OscillatableSpectrum &rhs)
	Assignment operator. More...
OscillatableSpectrum &	operator= (OscillatableSpectrum &&rhs)
Spectrum	Unoscillated () const
Spectrum	TrueEnergy () const
Spectrum	Oscillated (osc::IOscCalc *calc, int from, int to) const
Spectrum	Oscillated (osc::IOscCalcStan *calc, int from, int to) const
OscillatableSpectrum &	operator+= (const OscillatableSpectrum &rhs)
OscillatableSpectrum	operator+ (const OscillatableSpectrum &rhs) const
OscillatableSpectrum &	operator-= (const OscillatableSpectrum &rhs)
OscillatableSpectrum	operator- (const OscillatableSpectrum &rhs) const
void	SaveTo (TDirectory *dir, const std::string &name) const
const Var &	ReweightVar () const
	The variable that will be used to fill the y-axis. More...
void	Fill (double x, double y, double w=1)
TH2D *	ToTH2 (double pot) const
double	POT () const
double	Livetime () const
Spectrum	UnWeighted () const
Spectrum	WeightingVariable () const
Spectrum	WeightedBy (const Ratio &weights) const
	Reco spectrum with truth weights applied. More...
void	ReweightToTrueSpectrum (const Spectrum &target)
	Rescale bins so that WeightingVariable will return target. More...
void	ReweightToRecoSpectrum (const Spectrum &target)
	Recale bins so that Unweighted will return target. More...
ReweightableSpectrum &	operator+= (const ReweightableSpectrum &rhs)
ReweightableSpectrum	operator+ (const ReweightableSpectrum &rhs) const
ReweightableSpectrum &	operator-= (const ReweightableSpectrum &rhs)
ReweightableSpectrum	operator- (const ReweightableSpectrum &rhs) const
void	Clear ()
unsigned int	NDimensions () const
const std::vector< std::string > &	GetLabels () const
const std::vector< Binning > &	GetBinnings () const
Binning	GetTrueBinning () const
Eigen::MatrixXd	GetEigen (double pot) const

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

Protected Member Functions
	OscillatableSpectrum (const HistAxis &recoAxis, const Binning &binsY, const Var &rwVar)
	OscillatableSpectrum ()
	Constructor for Uninitialized() More...
template<class T >
Spectrum	_Oscillated (osc::_IOscCalc< T > *calc, int from, int to) const
ReweightableSpectrum &	PlusEqualsHelper (const ReweightableSpectrum &rhs, int sign)
void	RemoveLoader (ReweightableSpectrum **)
void	AddLoader (ReweightableSpectrum **)
void	_SaveTo (TDirectory *dir, const std::string &name, const std::string &type) const

Protected Attributes
ThreadLocal< OscCache >	fCache
Var	fRWVar
	What goes on the y axis? More...
Eigen::MatrixXd	fMat
double	fPOT
double	fLivetime
HistAxis	fAxisX
Binning	fBinsY
std::string	fTrueLabel
std::set< ReweightableSpectrum ** >	fReferences

Friends
class	SpectrumLoaderBase
class	SpectrumLoader
class	NullLoader

Detailed Description

Spectrum with true energy information, allowing it to be oscillated

Definition at line 36 of file OscillatableSpectrum.h.

Constructor & Destructor Documentation

ana::OscillatableSpectrum::OscillatableSpectrum	(	const std::string &	label,
		const Binning &	bins,
		SpectrumLoaderBase &	loader,
		const Var &	var,
		const Cut &	cut,
		const SystShifts &	shift = kNoShift,
		const Var &	wei = kUnweighted
	)

Definition at line 29 of file OscillatableSpectrum.cxx.

    : ReweightableSpectrum(loader,
                           HistAxis(label, bins, var),
                           HistAxis("True Energy (GeV)", kTrueEnergyBins, kTrueE),
                           cut, shift, wei)
  {
  }

ana::OscillatableSpectrum::OscillatableSpectrum	(	SpectrumLoaderBase &	loader,
		const HistAxis &	axis,
		const Cut &	cut,
		const SystShifts &	shift = kNoShift,
		const Var &	wei = kUnweighted
	)

Definition at line 43 of file OscillatableSpectrum.cxx.

    : ReweightableSpectrum(loader, axis,
                           HistAxis("True Energy (GeV)", kTrueEnergyBins, kTrueE),
                           cut, shift, wei)
  {
  }

ana::OscillatableSpectrum::OscillatableSpectrum	(	const Eigen::MatrixXd &&	mat,
		const HistAxis &	recoAxis,
		double	pot,
		double	livetime
	)

Definition at line 55 of file OscillatableSpectrum.cxx.

    : ReweightableSpectrum(std::move(mat), recoAxis,
                           HistAxis("True Energy (GeV)", kTrueEnergyBins, kTrueE),
                           pot, livetime)
  {
  }

ana::OscillatableSpectrum::~OscillatableSpectrum

(

)

Definition at line 65 of file OscillatableSpectrum.cxx.

  {
  }

ana::OscillatableSpectrum::OscillatableSpectrum

(

const OscillatableSpectrum &

rhs

)

Copy constructor.

Definition at line 70 of file OscillatableSpectrum.cxx.

References ana::assert(), fCache, and ana::ReweightableSpectrum::fReferences.

    : ReweightableSpectrum(rhs)
  {
    if(rhs.fCache->hash){
      fCache->spect = rhs.fCache->spect;
      fCache->hash = std::make_unique<TMD5>(*rhs.fCache->hash);
    }

    assert( rhs.fReferences.empty() ); // Copying with pending loads is unexpected
  }

ana::OscillatableSpectrum::OscillatableSpectrum

(

OscillatableSpectrum &&

rhs

)

Definition at line 82 of file OscillatableSpectrum.cxx.

References ana::assert(), and fCache.

    : ReweightableSpectrum(rhs)
  {
    if(rhs.fCache->hash){
      fCache->spect = std::move(rhs.fCache->spect);
      fCache->hash = std::move(rhs.fCache->hash);
    }

    assert( rhs.fReferences.empty() ); // Copying with pending loads is unexpected
  }

ana::OscillatableSpectrum::OscillatableSpectrum ( const HistAxis &  recoAxis, const Binning &  binsY, const Var &  rwVar  )

inlineprotected

Definition at line 102 of file OscillatableSpectrum.h.

      : ReweightableSpectrum(recoAxis, binsY, rwVar)
    {
    }

ana::OscillatableSpectrum::OscillatableSpectrum ( )

inlineprotected

Constructor for Uninitialized()

Definition at line 110 of file OscillatableSpectrum.h.

    {
    }

Member Function Documentation

template<class T >

Spectrum ana::OscillatableSpectrum::_Oscillated ( osc::_IOscCalc< T > *  calc, int  from, int  to  ) const

protected

Definition at line 137 of file OscillatableSpectrum.cxx.

References fCache, osc::_IOscCalc< T >::GetParamsHash(), samweb_client.utility::hash, runNovaSAM::ret, and ana::ReweightableSpectrum::WeightedBy().

Referenced by operator=(), and Oscillated().

  {
    TMD5* hash = calc->GetParamsHash();
    if(hash && fCache->hash && *hash == *fCache->hash){
      delete hash;
      return fCache->spect;
    }

    const OscCurve curve(calc, from, to);
    const Spectrum ret = WeightedBy(curve);
    if(hash){
      fCache->spect = ret;
      fCache->hash.reset(hash);
    }

    return ret;
  }

void ana::ReweightableSpectrum::_SaveTo ( TDirectory *  dir, const std::string &  name, const std::string &  type  ) const

protectedinherited

Definition at line 314 of file ReweightableSpectrum.cxx.

References dir, ana::ReweightableSpectrum::fAxisX, ana::ReweightableSpectrum::fLivetime, genie::utils::style::Format(), ana::ReweightableSpectrum::fPOT, ana::ReweightableSpectrum::fTrueLabel, ana::GenericHistAxis< T >::GetBinnings(), ana::GenericHistAxis< T >::GetLabels(), make_syst_table_plots::h, MECModelEnuComparisons::i, ana::GenericHistAxis< T >::NDimensions(), tmp, ana::ReweightableSpectrum::ToTH2(), and Write().

Referenced by ana::ReweightableSpectrum::ReweightableSpectrum(), ana::NDOscillatableSpectrum::SaveTo(), ana::ReweightableSpectrum::SaveTo(), and SaveTo().

  {
    TDirectory* tmp = gDirectory;

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

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

    TH2* h = ToTH2(fPOT);
    h->Write("hist");

    TH1D hPot("", "", 1, 0, 1);
    hPot.Fill(.5, fPOT);
    hPot.Write("pot");
    TH1D hLivetime("", "", 1, 0, 1);
    hLivetime.Fill(.5, fLivetime);
    hLivetime.Write("livetime");

    for(unsigned int i = 0; i < fAxisX.NDimensions(); ++i){
      TObjString(fAxisX.GetLabels()[i].c_str()).Write(TString::Format("label%d", i).Data());
      fAxisX.GetBinnings()[i].SaveTo(dir, TString::Format("bins%d", i).Data());
    }

    TObjString(fTrueLabel.c_str()).Write("true_label");

    dir->Write();
    delete dir;

    delete h;

    tmp->cd();
  }

void ana::ReweightableSpectrum::AddLoader ( ReweightableSpectrum **  ref )

protectedinherited

Definition at line 302 of file ReweightableSpectrum.cxx.

References ana::ReweightableSpectrum::fReferences.

Referenced by ana::ReweightableSpectrum::ReweightableSpectrum().

  {
    fReferences.insert(ref);
  }

void ana::ReweightableSpectrum::Clear ( )

inherited

Definition at line 290 of file ReweightableSpectrum.cxx.

References ana::ReweightableSpectrum::fMat.

Referenced by ana::PredictionExtrapSum::ComponentCC(), ana::ReweightableSpectrum::Livetime(), and Uninitialized().

  {
    fMat.setZero();
  }

void ana::ReweightableSpectrum::Fill ( double  x, double  y, double  w = 1  )

inherited

Definition at line 122 of file ReweightableSpectrum.cxx.

References ana::ReweightableSpectrum::fAxisX, ana::ReweightableSpectrum::fBinsY, ana::Binning::FindBin(), ana::ReweightableSpectrum::fMat, and ana::GenericHistAxis< T >::GetBins1D().

Referenced by ana::ReweightableSpectrum::ReweightVar(), and Uninitialized().

  {
    fMat(fBinsY.FindBin(y), fAxisX.GetBins1D().FindBin(x)) += w;
  }

const std::vector<Binning>& ana::ReweightableSpectrum::GetBinnings ( ) const

inlineinherited

Definition at line 87 of file ReweightableSpectrum.h.

References ana::ReweightableSpectrum::fAxisX, and ana::GenericHistAxis< T >::GetBinnings().

Referenced by ana::PredictionCombinePeriods::ComponentCC(), ana::MichelDecomp::MCToDCMPComp(), and ana::PredictionExtendToPeripheral::ReduceHelper().

{return fAxisX.GetBinnings();}

Eigen::MatrixXd ana::ReweightableSpectrum::GetEigen ( double  pot ) const

inlineinherited

Definition at line 90 of file ReweightableSpectrum.h.

References ana::ReweightableSpectrum::fMat, and ana::ReweightableSpectrum::fPOT.

Referenced by ana::PredictionCombinePeriods::ComponentCC(), ana::PredictionExtendToPeripheral::MergePeripheralOsc(), and ana::PredictionExtendToPeripheral::ReduceHelper().

{return fMat * pot/fPOT;}

const std::vector<std::string>& ana::ReweightableSpectrum::GetLabels ( ) const

inlineinherited

Definition at line 86 of file ReweightableSpectrum.h.

References ana::ReweightableSpectrum::fAxisX, and ana::GenericHistAxis< T >::GetLabels().

Referenced by ana::PredictionCombinePeriods::ComponentCC(), ana::MichelDecomp::MCToDCMPComp(), ana::PredictionExtendToPeripheral::MergePeripheralOsc(), and ana::PredictionExtendToPeripheral::ReduceHelper().

{return fAxisX.GetLabels();}

Binning ana::ReweightableSpectrum::GetTrueBinning ( ) const

inlineinherited

Definition at line 88 of file ReweightableSpectrum.h.

References ana::ReweightableSpectrum::fBinsY.

{return fBinsY;}

double ana::ReweightableSpectrum::Livetime ( ) const

inlineinherited

Definition at line 52 of file ReweightableSpectrum.h.

References ana::ReweightableSpectrum::Clear(), dir, ana::ReweightableSpectrum::fLivetime, ana::ReweightableSpectrum::LoadFrom(), ana::ReweightableSpectrum::operator+(), ana::ReweightableSpectrum::operator+=(), ana::ReweightableSpectrum::operator-(), ana::ReweightableSpectrum::operator-=(), ana::ReweightableSpectrum::ReweightToRecoSpectrum(), ana::ReweightableSpectrum::ReweightToTrueSpectrum(), ana::ReweightableSpectrum::SaveTo(), ana::ReweightableSpectrum::UnWeighted(), ana::ReweightableSpectrum::WeightedBy(), ana::ReweightableSpectrum::WeightingVariable(), and weights.

Referenced by ana::PredictionExtendToPeripheral::MergePeripheralOsc(), and ana::PredictionExtendToPeripheral::ReduceHelper().

{return fLivetime;}

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

static

Definition at line 214 of file OscillatableSpectrum.cxx.

References ana::assert(), ana::bins, dir, genie::utils::style::Format(), MECModelEnuComparisons::i, ana::kNoCut, ana::kNullLoader, label, PandAna.Demos.pi0_spectra::labels, ana::Binning::LoadFrom(), runNovaSAM::ret, cvnie::subdir, and getGoodRuns4SAM::tag.

Referenced by ana::FluxReweight::LoadFrom(), ana::TrivialExtrap::LoadFrom(), ana::NoReweight::LoadFrom(), ana::TruthReweight::LoadFrom(), and ana::RecoReweight::LoadFrom().

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

    DontAddDirectory guard;

    TObjString* tag = (TObjString*)dir->Get("type");
    assert(tag);
    assert(tag->GetString() == "OscillatableSpectrum");
    delete tag;

    TH2D* spect = (TH2D*)dir->Get("hist");
    assert(spect);
    TH1* hPot = (TH1*)dir->Get("pot");
    assert(hPot);
    TH1* hLivetime = (TH1*)dir->Get("livetime");
    assert(hLivetime);

    std::vector<std::string> labels;
    std::vector<Binning> bins;

    for(int i = 0; ; ++i){
      const std::string subname = TString::Format("bins%d", i).Data();
      TDirectory* subdir = dir->GetDirectory(subname.c_str());
      if(!subdir) break;
      delete subdir;
      bins.push_back(*Binning::LoadFrom(dir, subname));
      TObjString* label = (TObjString*)dir->Get(TString::Format("label%d", i));
      labels.push_back(label ? label->GetString().Data() : "");
      delete label;
    }

    delete dir;

    auto ret = std::make_unique<OscillatableSpectrum>(kNullLoader,
                                                      HistAxis(labels, bins),
                                                      kNoCut);

    // ROOT histogram storage is row-major, but Eigen is column-major by
    // default
    typedef Eigen::Matrix<double, Eigen::Dynamic, Eigen:: Dynamic, Eigen::RowMajor> MatRowMajor;
    ret->fMat = Eigen::Map<MatRowMajor>(spect->GetArray(),
                                        ret->fMat.rows(),
                                        ret->fMat.cols());

    delete spect;

    ret->fPOT = hPot->Integral(0, -1);
    ret->fLivetime = hLivetime->Integral(0, -1);

    delete hPot;
    delete hLivetime;

    return ret;
  }

unsigned int ana::ReweightableSpectrum::NDimensions ( ) const

inlineinherited

Definition at line 85 of file ReweightableSpectrum.h.

References ana::ReweightableSpectrum::fAxisX, and ana::GenericHistAxis< T >::NDimensions().

{return fAxisX.NDimensions();}

ReweightableSpectrum ana::ReweightableSpectrum::operator+ ( const ReweightableSpectrum &  rhs ) const

inherited

Definition at line 268 of file ReweightableSpectrum.cxx.

References runNovaSAM::ret.

Referenced by ana::ReweightableSpectrum::Livetime().

  {
    ReweightableSpectrum ret = *this;
    ret += rhs;
    return ret;
  }

OscillatableSpectrum ana::OscillatableSpectrum::operator+

(

const OscillatableSpectrum &

rhs

)

const

Definition at line 181 of file OscillatableSpectrum.cxx.

References runNovaSAM::ret.

  {
    OscillatableSpectrum ret = *this;
    ret += rhs;
    return ret;
  }

ReweightableSpectrum & ana::ReweightableSpectrum::operator+= ( const ReweightableSpectrum &  rhs )

inherited

Definition at line 262 of file ReweightableSpectrum.cxx.

References ana::ReweightableSpectrum::PlusEqualsHelper().

Referenced by ana::ReweightableSpectrum::Livetime(), and operator+=().

  {
    return PlusEqualsHelper(rhs, +1);
  }

OscillatableSpectrum & ana::OscillatableSpectrum::operator+=

(

const OscillatableSpectrum &

rhs

)

Definition at line 170 of file OscillatableSpectrum.cxx.

References fCache, and ana::ReweightableSpectrum::operator+=().

  {
    ReweightableSpectrum::operator+=(rhs);

    // invalidate
    fCache->hash.reset(nullptr);

    return *this;
  }

ReweightableSpectrum ana::ReweightableSpectrum::operator- ( const ReweightableSpectrum &  rhs ) const

inherited

Definition at line 282 of file ReweightableSpectrum.cxx.

References runNovaSAM::ret.

Referenced by ana::ReweightableSpectrum::Livetime().

  {
    ReweightableSpectrum ret = *this;
    ret -= rhs;
    return ret;
  }

OscillatableSpectrum ana::OscillatableSpectrum::operator-

(

const OscillatableSpectrum &

rhs

)

const

Definition at line 200 of file OscillatableSpectrum.cxx.

References runNovaSAM::ret.

  {
    OscillatableSpectrum ret = *this;
    ret -= rhs;
    return ret;
  }

ReweightableSpectrum & ana::ReweightableSpectrum::operator-= ( const ReweightableSpectrum &  rhs )

inherited

Definition at line 276 of file ReweightableSpectrum.cxx.

References ana::ReweightableSpectrum::PlusEqualsHelper().

Referenced by ana::ReweightableSpectrum::Livetime(), and operator-=().

  {
    return PlusEqualsHelper(rhs, -1);
  }

OscillatableSpectrum & ana::OscillatableSpectrum::operator-=

(

const OscillatableSpectrum &

rhs

)

Definition at line 189 of file OscillatableSpectrum.cxx.

References fCache, and ana::ReweightableSpectrum::operator-=().

  {
    ReweightableSpectrum::operator-=(rhs);

    // invalidate
    fCache->hash.reset(nullptr);

    return *this;
  }

OscillatableSpectrum & ana::OscillatableSpectrum::operator=

(

const OscillatableSpectrum &

rhs

)

Assignment operator.

Definition at line 94 of file OscillatableSpectrum.cxx.

References ana::assert(), fCache, ana::ReweightableSpectrum::fReferences, and ana::ReweightableSpectrum::operator=().

  {
    if(this == &rhs) return *this;

    ReweightableSpectrum::operator=(rhs);

    if(rhs.fCache->hash){
      fCache->spect = rhs.fCache->spect;
      fCache->hash = std::make_unique<TMD5>(*rhs.fCache->hash);
    }
    else{
      fCache->hash.reset();
    }

    assert( rhs.fReferences.empty() ); // Copying with pending loads is unexpected
    assert( fReferences.empty() ); // Copying with pending loads is unexpected

    return *this;
  }

OscillatableSpectrum & ana::OscillatableSpectrum::operator=

(

OscillatableSpectrum &&

rhs

)

Definition at line 115 of file OscillatableSpectrum.cxx.

References _Oscillated(), ana::assert(), fCache, ana::ReweightableSpectrum::fReferences, and ana::ReweightableSpectrum::operator=().

  {
    if(this == &rhs) return *this;

    ReweightableSpectrum::operator=(rhs);

    if(rhs.fCache->hash){
      fCache->spect = std::move(rhs.fCache->spect);
      fCache->hash = std::move(rhs.fCache->hash);
    }
    else{
      fCache->hash.reset();
    }

    assert( rhs.fReferences.empty() ); // Copying with pending loads is unexpected
    assert( fReferences.empty() ); // Copying with pending loads is unexpected

    return *this;
  }

Spectrum ana::OscillatableSpectrum::Oscillated	(	osc::IOscCalc *	calc,
		int	from,
		int	to
	)		const

Definition at line 156 of file OscillatableSpectrum.cxx.

References _Oscillated().

Referenced by ana::PredictionExtendToPeripheral::_PredictComponent(), demo3(), make_extrap_figure(), numu_cut_flow(), ana::FDPredictionSterile::Predict(), ana::PredictionSterile::PredictComponent(), and ana::FDPredictionSterile::PredictComponent().

  {
    return _Oscillated(calc, from, to);
  }

Spectrum ana::OscillatableSpectrum::Oscillated	(	osc::IOscCalcStan *	calc,
		int	from,
		int	to
	)		const

Definition at line 163 of file OscillatableSpectrum.cxx.

References _Oscillated().

  {
    return _Oscillated(calc, from, to);
  }

ReweightableSpectrum & ana::ReweightableSpectrum::PlusEqualsHelper ( const ReweightableSpectrum &  rhs, int  sign  )

protectedinherited

Definition at line 191 of file ReweightableSpectrum.cxx.

References ana::AlmostEqual(), om::cout, allTimeWatchdog::endl, ana::ReweightableSpectrum::fLivetime, ana::ReweightableSpectrum::fMat, and ana::ReweightableSpectrum::fPOT.

Referenced by ana::ReweightableSpectrum::operator+=(), ana::ReweightableSpectrum::operator-=(), and ana::ReweightableSpectrum::ReweightableSpectrum().

  {
    // In this case it would be OK to have no POT/livetime
    if(rhs.fMat.sum() == 0) return *this;


    if((!fPOT && !fLivetime) || (!rhs.fPOT && !rhs.fLivetime)){
      std::cout << "Error: can't sum ReweightableSpectrum with no POT or livetime."
                << fPOT << " " << rhs.fPOT
                << std::endl;
      //      abort();
      return *this;
    }

    if(!fLivetime && !rhs.fPOT){
      std::cout << "Error: can't sum ReweightableSpectrum with POT ("
                << fPOT << ") but no livetime and ReweightableSpectrum with livetime ("
                << rhs.fLivetime << " sec) but no POT." << std::endl;
      abort();
    }

    if(!fPOT && !rhs.fLivetime){
      std::cout << "Error: can't sum ReweightableSpectrum with livetime ("
                << fLivetime << " sec) but no POT and ReweightableSpectrum with POT ("
                << rhs.fPOT << ") but no livetime." << std::endl;
      abort();
    }

    // And now there are still a bunch of good cases to consider

    if(fPOT && rhs.fPOT){
      // Scale by POT when possible
      fMat += rhs.fMat * sign*fPOT/rhs.fPOT;

      if(fLivetime && rhs.fLivetime){
        // If POT/livetime ratios match, keep regular lifetime, otherwise zero
        // it out.
        if(AlmostEqual(fLivetime*rhs.fPOT, rhs.fLivetime*fPOT))
          fLivetime = 0;
      }
      if(!fLivetime && rhs.fLivetime){
        // If the RHS has a livetime and we don't, copy it in (suitably scaled)
        fLivetime = rhs.fLivetime * fPOT/rhs.fPOT;
      }
      // Otherwise, keep our own livetime (if any)

      return *this;
    }

    if(fLivetime && rhs.fLivetime){
      // Scale by livetime, the only thing in common
      fMat += rhs.fMat * sign*fLivetime/rhs.fLivetime;

      if(!fPOT && rhs.fPOT){
        // If the RHS has a POT and we don't, copy it in (suitably scaled)
        fPOT = rhs.fPOT * fLivetime/rhs.fLivetime;
      }
      // Otherwise, keep our own POT (if any)

      return *this;
    }

    // That should have been all the cases. I definitely want to know what
    // happened if it wasn't.
    std::cout << "ReweightableSpectrum::operator+=(). How did we get here? "
              << fPOT << " " << fLivetime << " "
              << rhs.fPOT << " " << rhs.fLivetime << std::endl;
    abort();
  }

double ana::ReweightableSpectrum::POT ( ) const

inlineinherited

Definition at line 51 of file ReweightableSpectrum.h.

References ana::ReweightableSpectrum::fPOT.

Referenced by ana::PredictionExtendToPeripheral::MergePeripheralOsc(), ana::PlotAllRecoTrue(), and ana::PredictionExtendToPeripheral::ReduceHelper().

{return fPOT;}

void ana::ReweightableSpectrum::RemoveLoader ( ReweightableSpectrum **  ref )

protectedinherited

Definition at line 296 of file ReweightableSpectrum.cxx.

References ana::ReweightableSpectrum::fReferences.

Referenced by ana::ReweightableSpectrum::ReweightableSpectrum().

  {
    fReferences.erase(ref);
  }

void ana::ReweightableSpectrum::ReweightToRecoSpectrum ( const Spectrum &  target )

inherited

Recale bins so that Unweighted will return target.

Definition at line 180 of file ReweightableSpectrum.cxx.

References ana::ReweightableSpectrum::fMat, ana::Ratio::GetEigen(), ratio(), and ana::ReweightableSpectrum::UnWeighted().

Referenced by ana::FluxReweight::Eval(), ana::TruthReweight::Eval(), ana::RecoReweight::Eval(), ana::PredictionExtendToPeripheral::ExtendRecoWeightOscillatable(), ana::ReweightableSpectrum::Livetime(), ana::TruthReweight::SavePlots(), ana::UnfoldIterative::Truth(), and Uninitialized().

  {
    // This is a big component of what extrapolations do, so it should be fast

    const Ratio ratio(target, UnWeighted());
    // We want to multiply all the columns by this ratio

    fMat *= ratio.GetEigen().matrix().asDiagonal();
  }

void ana::ReweightableSpectrum::ReweightToTrueSpectrum ( const Spectrum &  target )

inherited

Rescale bins so that WeightingVariable will return target.

Definition at line 169 of file ReweightableSpectrum.cxx.

References ana::ReweightableSpectrum::fMat, ana::Ratio::GetEigen(), ratio(), and ana::ReweightableSpectrum::WeightingVariable().

Referenced by ana::TruthReweight::Eval(), ana::ReweightableSpectrum::Livetime(), ana::UnfoldIterative::Truth(), and Uninitialized().

  {
    // This is a big component of what extrapolations do, so it should be fast

    const Ratio ratio(target, WeightingVariable());
    // We want to multiply all the rows by this ratio, so left-multiply

    fMat = ratio.GetEigen().matrix().asDiagonal() * fMat;
  }

const Var& ana::ReweightableSpectrum::ReweightVar ( ) const

inlineinherited

The variable that will be used to fill the y-axis.

By convention, return zero if the information can't be obtained, and this event will be skipped.

Definition at line 45 of file ReweightableSpectrum.h.

References ana::ReweightableSpectrum::Fill(), ana::ReweightableSpectrum::fRWVar, pot, ana::ReweightableSpectrum::ToTH2(), w, submit_syst::x, and submit_syst::y.

Referenced by ana::SpectrumLoaderBase::AddReweightableSpectrum().

{return fRWVar;}

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

Definition at line 208 of file OscillatableSpectrum.cxx.

References ana::ReweightableSpectrum::_SaveTo().

Referenced by ana::TruthReweight::SavePlots(), ana::TrivialExtrap::SaveTo(), ana::NoReweight::SaveTo(), ana::TruthReweight::SaveTo(), and ana::RecoReweight::SaveTo().

  {
    _SaveTo(dir, name, "OscillatableSpectrum");
  }

TH2D * ana::ReweightableSpectrum::ToTH2 ( double  pot ) const

inherited

Definition at line 94 of file ReweightableSpectrum.cxx.

References ana::assert(), ana::ReweightableSpectrum::fAxisX, ana::ReweightableSpectrum::fBinsY, ana::ReweightableSpectrum::fMat, ana::ReweightableSpectrum::fPOT, ana::ReweightableSpectrum::fTrueLabel, ana::GenericHistAxis< T >::GetBins1D(), ana::GenericHistAxis< T >::GetLabel1D(), MECModelEnuComparisons::i, calib::j, ana::MakeTH2D(), runNovaSAM::ret, and ana::UniqueName().

Referenced by ana::ReweightableSpectrum::_SaveTo(), demo2(), demo3(), ana::MichelDecomp::GetSum(), ana::MichelDecomp::GetTemplateContent(), ana::MichelDecomp::MDCMPLogLikelihood(), ana::PlotAllRecoTrue(), ana::CrossSectionAnalysis::PlotRecoToTrueMatrix(), ana::SingleNucAnalysis::PlotRecoToTrueMatrix2D(), ana::NumuCC2p2hAnalysis::PlotRecoToTrueMatrix2D(), ana::ReweightableSpectrum::ReweightVar(), ana::MichelDecomp::SaveTempPlots(), ana::UnfoldSVD::Truth(), Unfold1D(), Unfold3D(), ana::UnfoldTikhonov::UnfoldTikhonov(), and Uninitialized().

  {
    // Could have a file temporarily open
    DontAddDirectory guard;

    TH2D* ret = MakeTH2D(UniqueName().c_str(), "", fAxisX.GetBins1D(), fBinsY);

    for(int i = 0; i < fMat.rows(); ++i){
      for(int j = 0; j < fMat.cols(); ++j){
        ret->SetBinContent(j, i, fMat(i, j));
      }
    }

    if(fPOT){
      ret->Scale(pot/fPOT);
    }
    else{
      // How did it get events with no POT?
      assert(ret->Integral() == 0);
    }

    ret->GetXaxis()->SetTitle(fAxisX.GetLabel1D().c_str());
    ret->GetYaxis()->SetTitle(fTrueLabel.c_str());

    return ret;
  }

Spectrum ana::OscillatableSpectrum::TrueEnergy ( ) const

inline

Definition at line 86 of file OscillatableSpectrum.h.

References calc, dir, ana::LoadFrom(), novaddt::operator+(), novaddt::operator+=(), novaddt::operator-(), and ana::SaveTo().

Referenced by ana::TruthReweight::Eval(), and ana::TruthReweight::SavePlots().

{return WeightingVariable();}

static OscillatableSpectrum ana::OscillatableSpectrum::Uninitialized ( )

inlinestatic

The only valid thing to do with such a spectrum is to assign something else into it.

Definition at line 63 of file OscillatableSpectrum.h.

References ana::ReweightableSpectrum::Clear(), ana::ReweightableSpectrum::Fill(), ana::ReweightableSpectrum::ReweightToRecoSpectrum(), ana::ReweightableSpectrum::ReweightToTrueSpectrum(), and ana::ReweightableSpectrum::ToTH2().

{return OscillatableSpectrum();}

Spectrum ana::OscillatableSpectrum::Unoscillated ( ) const

inline

Definition at line 85 of file OscillatableSpectrum.h.

Referenced by ana::TruthReweight::Eval(), ana::RecoReweight::Eval(), ana::PredictionExtendToPeripheral::ExtendRecoWeightOscillatable(), ana::ModularExtrapSterile::NCTotalComponent(), ND_DataMC(), nue_decomp_scales(), and ana::TruthReweight::SavePlots().

{return UnWeighted();}

Spectrum ana::ReweightableSpectrum::UnWeighted ( ) const

inherited

Definition at line 140 of file ReweightableSpectrum.cxx.

References ana::Hist::Adopt(), ana::ReweightableSpectrum::fAxisX, ana::ReweightableSpectrum::fLivetime, ana::ReweightableSpectrum::fMat, ana::ReweightableSpectrum::fPOT, and ana::ProjectionX().

Referenced by ana::ReweightableSpectrum::Livetime(), ana::ReweightableSpectrum::ReweightToRecoSpectrum(), and ana::UnfoldSVD::Truth().

  {
    return Spectrum(Hist::Adopt(ProjectionX(fMat)), fAxisX, fPOT, fLivetime);
  }

Spectrum ana::ReweightableSpectrum::WeightedBy ( const Ratio &  weights ) const

inherited

Reco spectrum with truth weights applied.

Definition at line 152 of file ReweightableSpectrum.cxx.

References ana::Hist::Adopt(), ana::Hist::AdoptStan(), ana::ReweightableSpectrum::fAxisX, ana::ReweightableSpectrum::fLivetime, ana::ReweightableSpectrum::fMat, ana::ReweightableSpectrum::fPOT, ana::Ratio::GetEigen(), ana::Ratio::GetEigenStan(), and ana::Ratio::HasStan().

Referenced by _Oscillated(), ana::ReweightableSpectrum::Livetime(), and ana::NDOscillatableSpectrum::Oscillated().

  {
    if(!ws.HasStan()){
      const Eigen::VectorXd& vec = ws.GetEigen();

      return Spectrum(Hist::Adopt(Eigen::ArrayXd(vec.transpose() * fMat)),
                      fAxisX, fPOT, fLivetime);
    }
    else{
      const Eigen::VectorXstan& vec = ws.GetEigenStan();

      return Spectrum(Hist::AdoptStan(vec.transpose() * fMat),
                      fAxisX, fPOT, fLivetime);
    }
  }

Spectrum ana::ReweightableSpectrum::WeightingVariable ( ) const

inherited

Definition at line 146 of file ReweightableSpectrum.cxx.

References ana::Hist::Adopt(), ana::ReweightableSpectrum::fAxisX, ana::ReweightableSpectrum::fLivetime, ana::ReweightableSpectrum::fMat, ana::ReweightableSpectrum::fPOT, and ana::ProjectionY().

Referenced by ana::ReweightableSpectrum::Livetime(), ana::ReweightableSpectrum::ReweightToTrueSpectrum(), ana::UnfoldIterative::Truth(), and ana::UnfoldSVD::Truth().

  {
    return Spectrum(Hist::Adopt(ProjectionY(fMat)), fAxisX, fPOT, fLivetime);
  }

Friends And Related Function Documentation

friend class NullLoader

friend

Definition at line 41 of file OscillatableSpectrum.h.

friend class SpectrumLoader

friend

Definition at line 40 of file OscillatableSpectrum.h.

friend class SpectrumLoaderBase

friend

Definition at line 39 of file OscillatableSpectrum.h.

Member Data Documentation

HistAxis ana::ReweightableSpectrum::fAxisX

protectedinherited

Definition at line 123 of file ReweightableSpectrum.h.

Referenced by ana::ReweightableSpectrum::_SaveTo(), ana::ReweightableSpectrum::Fill(), ana::ReweightableSpectrum::GetBinnings(), ana::ReweightableSpectrum::GetLabels(), ana::ReweightableSpectrum::NDimensions(), ana::ReweightableSpectrum::operator=(), ana::ReweightableSpectrum::ToTH2(), ana::ReweightableSpectrum::UnWeighted(), ana::ReweightableSpectrum::WeightedBy(), and ana::ReweightableSpectrum::WeightingVariable().

Binning ana::ReweightableSpectrum::fBinsY

protectedinherited

Definition at line 125 of file ReweightableSpectrum.h.

Referenced by ana::ReweightableSpectrum::Fill(), ana::ReweightableSpectrum::GetTrueBinning(), ana::ReweightableSpectrum::operator=(), ana::ReweightableSpectrum::ReweightableSpectrum(), and ana::ReweightableSpectrum::ToTH2().

ThreadLocal<OscCache> ana::OscillatableSpectrum::fCache

mutableprotected

Definition at line 116 of file OscillatableSpectrum.h.

Referenced by _Oscillated(), operator+=(), operator-=(), operator=(), and OscillatableSpectrum().

double ana::ReweightableSpectrum::fLivetime

protectedinherited

Definition at line 121 of file ReweightableSpectrum.h.

Referenced by ana::ReweightableSpectrum::_SaveTo(), ana::ReweightableSpectrum::Livetime(), ana::ReweightableSpectrum::operator=(), ana::ReweightableSpectrum::PlusEqualsHelper(), ana::ReweightableSpectrum::ReweightableSpectrum(), ana::ReweightableSpectrum::UnWeighted(), ana::ReweightableSpectrum::WeightedBy(), and ana::ReweightableSpectrum::WeightingVariable().

Eigen::MatrixXd ana::ReweightableSpectrum::fMat

protectedinherited

Definition at line 119 of file ReweightableSpectrum.h.

Referenced by ana::ReweightableSpectrum::Clear(), ana::ReweightableSpectrum::Fill(), ana::ReweightableSpectrum::GetEigen(), ana::ReweightableSpectrum::operator=(), ana::ReweightableSpectrum::PlusEqualsHelper(), ana::ReweightableSpectrum::ReweightableSpectrum(), ana::ReweightableSpectrum::ReweightToRecoSpectrum(), ana::ReweightableSpectrum::ReweightToTrueSpectrum(), ana::NDOscillatableSpectrum::Scale(), ana::ReweightableSpectrum::ToTH2(), ana::ReweightableSpectrum::UnWeighted(), ana::ReweightableSpectrum::WeightedBy(), and ana::ReweightableSpectrum::WeightingVariable().

double ana::ReweightableSpectrum::fPOT

protectedinherited

Definition at line 120 of file ReweightableSpectrum.h.

Referenced by ana::ReweightableSpectrum::_SaveTo(), ana::ReweightableSpectrum::GetEigen(), ana::ReweightableSpectrum::operator=(), ana::ReweightableSpectrum::PlusEqualsHelper(), ana::ReweightableSpectrum::POT(), ana::ReweightableSpectrum::ReweightableSpectrum(), ana::ReweightableSpectrum::ToTH2(), ana::ReweightableSpectrum::UnWeighted(), ana::ReweightableSpectrum::WeightedBy(), and ana::ReweightableSpectrum::WeightingVariable().

std::set<ReweightableSpectrum**> ana::ReweightableSpectrum::fReferences

protectedinherited

Things that point at this ReweightableSpectrum. Maintained by SpectrumLoader

Definition at line 131 of file ReweightableSpectrum.h.

Referenced by ana::ReweightableSpectrum::AddLoader(), ana::NDOscillatableSpectrum::NDOscillatableSpectrum(), ana::ReweightableSpectrum::operator=(), ana::NDOscillatableSpectrum::operator=(), operator=(), OscillatableSpectrum(), ana::ReweightableSpectrum::RemoveLoader(), ana::ReweightableSpectrum::ReweightableSpectrum(), and ana::ReweightableSpectrum::~ReweightableSpectrum().

Var ana::ReweightableSpectrum::fRWVar

protectedinherited

What goes on the y axis?

Definition at line 117 of file ReweightableSpectrum.h.

Referenced by ana::ReweightableSpectrum::operator=(), and ana::ReweightableSpectrum::ReweightVar().

std::string ana::ReweightableSpectrum::fTrueLabel

protectedinherited

Definition at line 127 of file ReweightableSpectrum.h.

Referenced by ana::ReweightableSpectrum::_SaveTo(), ana::ReweightableSpectrum::ReweightableSpectrum(), and ana::ReweightableSpectrum::ToTH2().

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-09-23/CAFAna/Core/OscillatableSpectrum.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-09-23/CAFAna/Core/OscillatableSpectrum.cxx