ana::OscillatableSpectrum Class Reference

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

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-11-28/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
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
const std::vector< Binning > &	GetTrueBinnings () 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 ()
	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
Eigen::MatrixXd	fMat
double	fPOT
double	fLivetime
LabelsAndBins	fAxisX
LabelsAndBins	fAxisY
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 39 of file OscillatableSpectrum.cxx.

    : ReweightableSpectrum(loader,
                           HistAxis(label, bins, var),
                           HistAxis("True Energy (GeV)", kTrueEnergyBins, kTrueE),
                           kHasNu && 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 53 of file OscillatableSpectrum.cxx.

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

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

Definition at line 65 of file OscillatableSpectrum.cxx.

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

ana::OscillatableSpectrum::~OscillatableSpectrum

(

)

Definition at line 75 of file OscillatableSpectrum.cxx.

  {
  }

ana::OscillatableSpectrum::OscillatableSpectrum

(

const OscillatableSpectrum &

rhs

)

Copy constructor.

Definition at line 80 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 92 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 ( )

inlineprotected

Constructor for Uninitialized()

Definition at line 103 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 147 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 285 of file ReweightableSpectrum.cxx.

References dir, ana::ReweightableSpectrum::fAxisX, ana::ReweightableSpectrum::fAxisY, ana::ReweightableSpectrum::fLivetime, genie::utils::style::Format(), ana::ReweightableSpectrum::fPOT, ana::LabelsAndBins::GetBinnings(), ana::LabelsAndBins::GetLabels(), make_syst_table_plots::h, MECModelEnuComparisons::i, ana::LabelsAndBins::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());
    }

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

    dir->Write();
    delete dir;

    delete h;

    tmp->cd();
  }

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

protectedinherited

Definition at line 273 of file ReweightableSpectrum.cxx.

References ana::ReweightableSpectrum::fReferences.

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

  {
    fReferences.insert(ref);
  }

void ana::ReweightableSpectrum::Clear ( void  )

inherited

Definition at line 261 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 93 of file ReweightableSpectrum.cxx.

References ana::ReweightableSpectrum::fAxisX, ana::ReweightableSpectrum::fAxisY, ana::Binning::FindBin(), ana::ReweightableSpectrum::fMat, and ana::LabelsAndBins::GetBins1D().

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

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

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

inlineinherited

Definition at line 83 of file ReweightableSpectrum.h.

References ana::ReweightableSpectrum::fAxisX, and ana::LabelsAndBins::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 86 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 82 of file ReweightableSpectrum.h.

References ana::ReweightableSpectrum::fAxisX, and ana::LabelsAndBins::GetLabels().

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

{return fAxisX.GetLabels();}

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

inlineinherited

Definition at line 84 of file ReweightableSpectrum.h.

References ana::ReweightableSpectrum::fAxisY, and ana::LabelsAndBins::GetBinnings().

{return fAxisY.GetBinnings();}

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(), string, 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 224 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, string, 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 81 of file ReweightableSpectrum.h.

References ana::ReweightableSpectrum::fAxisX, and ana::LabelsAndBins::NDimensions().

{return fAxisX.NDimensions();}

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

inherited

Definition at line 239 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 191 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 233 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 180 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 253 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 210 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 247 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 199 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 104 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 125 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 166 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 173 of file OscillatableSpectrum.cxx.

References _Oscillated().

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

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

protectedinherited

Definition at line 162 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 267 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 151 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 140 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;
  }

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

Definition at line 218 of file OscillatableSpectrum.cxx.

References ana::ReweightableSpectrum::_SaveTo().

Referenced by ana::TruthReweight::SavePlots(), ana::FluxReweight::SaveTo(), 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 65 of file ReweightableSpectrum.cxx.

References ana::assert(), ana::ReweightableSpectrum::fAxisX, ana::ReweightableSpectrum::fAxisY, ana::ReweightableSpectrum::fMat, ana::ReweightableSpectrum::fPOT, ana::LabelsAndBins::GetBins1D(), ana::LabelsAndBins::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::MichelDecomp::SaveTempPlots(), ana::UnfoldSVD::Truth(), Unfold1D(), Unfold3D(), ana::UnfoldTikhonov::UnfoldTikhonov(), ana::ReweightableSpectrum::Uninitialized(), and Uninitialized().

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

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

    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(fAxisY.GetLabel1D().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-(), ana::SaveTo(), and string.

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::FluxReweight::Eval(), ana::TruthReweight::Eval(), ana::RecoReweight::Eval(), ana::PredictionExtendToPeripheral::ExtendRecoWeightOscillatable(), ana::ModularExtrapSterile::NCTotalComponent(), ND_DataMC(), nue_decomp_scales(), ana::FluxReweight::SavePlots(), and ana::TruthReweight::SavePlots().

{return UnWeighted();}

Spectrum ana::ReweightableSpectrum::UnWeighted ( ) const

inherited

Definition at line 111 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 123 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 117 of file ReweightableSpectrum.cxx.

References ana::Hist::Adopt(), ana::ReweightableSpectrum::fAxisY, 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)), fAxisY, 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

LabelsAndBins ana::ReweightableSpectrum::fAxisX

protectedinherited

Definition at line 117 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(), and ana::ReweightableSpectrum::WeightedBy().

LabelsAndBins ana::ReweightableSpectrum::fAxisY

protectedinherited

Definition at line 117 of file ReweightableSpectrum.h.

Referenced by ana::ReweightableSpectrum::_SaveTo(), ana::ReweightableSpectrum::Fill(), ana::ReweightableSpectrum::GetTrueBinnings(), ana::ReweightableSpectrum::operator=(), ana::ReweightableSpectrum::ToTH2(), and ana::ReweightableSpectrum::WeightingVariable().

ThreadLocal<OscCache> ana::OscillatableSpectrum::fCache

mutableprotected

Definition at line 109 of file OscillatableSpectrum.h.

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

double ana::ReweightableSpectrum::fLivetime

protectedinherited

Definition at line 115 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 113 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 114 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 121 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().

---

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

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