ana::Spectrum Class Reference

Representation of a spectrum in any variable, with associated POT. More...

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-07-07/CAFAna/Core/Spectrum.h"

Public Types
enum	ESparse { kDense, kSparse, kFloat }

Public Member Functions
	Spectrum (const std::string &label, const Binning &bins, SpectrumLoaderBase &loader, const Var &var, const Cut &cut, const SystShifts &shift=kNoShift, const Var &wei=kUnweighted)
	Spectrum (const std::string &label, const Binning &bins, SpectrumLoaderBase &loader, const MultiVar &var, const Cut &cut, const SystShifts &shift=kNoShift, const Var &wei=kUnweighted)
	Spectrum (const std::string &label, const Binning &bins, SpectrumLoaderBase &loader, const NuTruthVar &var, const NuTruthCut &cut, const SystShifts &shifts=kNoShift, const NuTruthVar &wei=kNuTruthUnweighted)
	Spectrum (const std::string &label, SpectrumLoaderBase &loader, const Binning &binsx, const NuTruthVar &varx, const Binning &binsy, const NuTruthVar &vary, const NuTruthCut &cut, const SystShifts &shifts=kNoShift, const NuTruthVar &wei=kNuTruthUnweighted)
	Spectrum (SpectrumLoaderBase &loader, const HistAxis &axis, const Cut &cut, const SystShifts &shift=kNoShift, const Var &wei=kUnweighted)
	Spectrum (SpectrumLoaderBase &loader, const MultiVarHistAxis &axis, const Cut &cut, const SystShifts &shift=kNoShift, const Var &wei=kUnweighted)
	Spectrum (SpectrumLoaderBase &loader, const NuTruthHistAxis &axis, const NuTruthCut &cut, const SystShifts &shifts=kNoShift, const NuTruthVar &weight=kNuTruthUnweighted)
	Spectrum (const std::string &label, const Binning &bins, ESparse sparse=kDense)
	Spectrum (const std::string &label, double pot, double livetime, const Binning &bins)
	Spectrum (TH1 *h, const std::vector< std::string > &labels, const std::vector< Binning > &bins, double pot, double livetime)
	Copies h. More...
	Spectrum (TH1 *h, double pot, double livetime)
	Spectrum (std::unique_ptr< TH1D > h, const std::vector< std::string > &labels, const std::vector< Binning > &bins, double pot, double livetime)
	Takes possession of h. More...
	Spectrum (std::unique_ptr< TH1F > h, const std::vector< std::string > &labels, const std::vector< Binning > &bins, double pot, double livetime)
	Spectrum (const std::string &label, SpectrumLoaderBase &loader, const Binning &binsx, const Var &varx, const Binning &binsy, const Var &vary, const Cut &cut, const SystShifts &shift=kNoShift, const Var &wei=kUnweighted)
	2D Spectrum of two Vars More...
	Spectrum (const std::string &label, SpectrumLoaderBase &loader, const Binning &binsx, const MultiVar &varx, const Binning &binsy, const MultiVar &vary, const Cut &cut, const SystShifts &shift=kNoShift, const Var &wei=kUnweighted)
	2D Spectrum of two MultiVars More...
	Spectrum (SpectrumLoaderBase &loader, const HistAxis &xAxis, const HistAxis &yAxis, const Cut &cut, const SystShifts &shift=kNoShift, const Var &wei=kUnweighted)
	2D Spectrum taking 2 HistAxis More...
	Spectrum (const std::string &xLabel, const std::string &yLabel, SpectrumLoaderBase &loader, const Binning &binsx, const Var &varx, const Binning &binsy, const Var &vary, const Cut &cut, const SystShifts &shift=kNoShift, const Var &wei=kUnweighted)
	Spectrum (const std::string &xLabel, const std::string &yLabel, SpectrumLoaderBase &loader, const Binning &binsx, const MultiVar &varx, const Binning &binsy, const MultiVar &vary, const Cut &cut, const SystShifts &shift=kNoShift, const Var &wei=kUnweighted)
	Spectrum (const std::string &label, SpectrumLoaderBase &loader, const Binning &binsx, const Var &varx, const Binning &binsy, const Var &vary, const Binning &binsz, const Var &varz, const Cut &cut, const SystShifts &shift=kNoShift, const Var &wei=kUnweighted, ESparse sparse=kDense)
	3D Spectrum of three Vars More...
	Spectrum (const std::string &xLabel, const std::string &yLabel, const std::string &zLabel, SpectrumLoaderBase &loader, const Binning &binsx, const Var &varx, const Binning &binsy, const Var &vary, const Binning &binsz, const Var &varz, const Cut &cut, const SystShifts &shift=kNoShift, const Var &wei=kUnweighted, ESparse sparse=kDense)
	Spectrum (SpectrumLoaderBase &loader, const HistAxis &xAxis, const HistAxis &yAxis, const HistAxis &zAxis, const Cut &cut, const SystShifts &shift=kNoShift, const Var &wei=kUnweighted, ESparse sparse=kDense)
	3D Spectrum taking 3 HistAxis More...
virtual	~Spectrum ()
	Spectrum (const Spectrum &rhs)
	Spectrum (Spectrum &&rhs)
Spectrum &	operator= (const Spectrum &rhs)
Spectrum &	operator= (Spectrum &&rhs)
	Spectrum (const SpectrumStan &rhs)
	Copy conversion from SpectrumStan. More...
Spectrum &	operator= (const SpectrumStan &rhs)
void	Fill (double x, double w=1)
TH1D *	ToTH1 (double exposure, Color_t col=kBlack, Style_t style=kSolid, EExposureType expotype=kPOT, EBinType bintype=kBinContent) const
	Histogram made from this Spectrum, scaled to some exposure. More...
TH1D *	ToTH1 (double exposure, EExposureType expotype, EBinType bintype=kBinContent) const
	Histogram made from this Spectrum, scaled to some exposure. More...
TH2 *	ToTH2 (double exposure, EExposureType expotype=kPOT, EBinType bintype=kBinContent) const
	Spectrum must be 2D to obtain TH2. More...
TH2 *	ToTH2NormX (double exposure, EExposureType expotype=kPOT) const
	Spectrum must be 2D to obtain TH2. Normalized to X axis. More...
TH3 *	ToTH3 (double exposure, EExposureType expotype=kPOT, EBinType bintype=kBinContent) const
	Spectrum must be 3D to obtain TH3. More...
double	Integral (double exposure, double *err=0, EExposureType expotype=kPOT) const
	Return total number of events scaled to pot. More...
double	Mean () const
	Return mean of 1D histogram. More...
Spectrum	MockData (double pot, int idx=0) const
	Mock data is FakeData with Poisson fluctuations applied. More...
Spectrum	FakeData (double pot) const
	Fake data is a MC spectrum scaled to the POT expected in the data. More...
Spectrum	FakeData (double pot, double livetime) const
double	POT () const
double	Livetime () const
	Seconds. For informational purposes only. No calculations use this. More...
void	OverridePOT (double newpot)
	DO NOT USE UNLESS YOU ARE 110% CERTAIN THERE ISN'T A BETTER WAY! More...
void	OverrideLivetime (double newlive)
	DO NOT USE UNLESS YOU ARE 110% CERTAIN THERE ISN'T A BETTER WAY! More...
void	Clear ()
void	Scale (double c)
	Multiply this spectrum by a constant c. More...
Spectrum &	operator+= (const Spectrum &rhs)
Spectrum	operator+ (const Spectrum &rhs) const
Spectrum &	operator-= (const Spectrum &rhs)
Spectrum	operator- (const Spectrum &rhs) const
Spectrum &	operator*= (const Ratio &rhs)
Spectrum	operator* (const Ratio &rhs) const
Spectrum &	operator/= (const Ratio &rhs)
Spectrum	operator/ (const Ratio &rhs) const
void	SaveTo (TDirectory *dir) const
unsigned int	NDimensions () const
std::vector< std::string >	GetLabels () const
std::vector< Binning >	GetBinnings () const

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

Protected Member Functions
	Spectrum (const std::vector< std::string > &labels, const std::vector< Binning > &bins, ESparse sparse=kDense)
Binning	GetMultiDBinning () const
void	ConstructHistogram (ESparse sparse=kDense)
void	RemoveLoader (SpectrumLoaderBase *)
void	AddLoader (SpectrumLoaderBase *)
Spectrum &	PlusEqualsHelper (const Spectrum &rhs, int sign)
	Helper for operator+= and operator-=. More...

Protected Attributes
TH1D *	fHistD
TH1F *	fHistF
THnSparseD *	fHistSparse
double	fPOT
double	fLivetime
std::set< SpectrumLoaderBase * >	fLoaderCount
	This count is maintained by SpectrumLoader, as a sanity check. More...
std::vector< std::string >	fLabels
std::vector< Binning >	fBins

Friends
class	SpectrumStan
class	SpectrumLoaderBase
class	SpectrumLoader
class	NullLoader

Detailed Description

Representation of a spectrum in any variable, with associated POT.

Definition at line 37 of file Spectrum.h.

Member Enumeration Documentation

enum ana::Spectrum::ESparse

Enumerator
kDense
kSparse
kFloat

Definition at line 45 of file Spectrum.h.

{kDense, kSparse, kFloat};

Constructor & Destructor Documentation

ana::Spectrum::Spectrum	(	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 47 of file Spectrum.cxx.

References ana::SpectrumLoaderBase::AddSpectrum().

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

    : Spectrum(label, bins)
  {
    loader.AddSpectrum(*this, var, cut, shift, wei);
  }

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

Definition at line 59 of file Spectrum.cxx.

References ana::SpectrumLoaderBase::AddSpectrum().

    : Spectrum(label, bins)
  {
    loader.AddSpectrum(*this, var, cut, shift, wei);
  }

ana::Spectrum::Spectrum	(	const std::string &	label,
		const Binning &	bins,
		SpectrumLoaderBase &	loader,
		const NuTruthVar &	var,
		const NuTruthCut &	cut,
		const SystShifts &	shifts = kNoShift,
		const NuTruthVar &	wei = kNuTruthUnweighted
	)

Definition at line 71 of file Spectrum.cxx.

References ana::SpectrumLoaderBase::AddSpectrum().

    : Spectrum(label, bins)
  {
    loader.AddSpectrum(*this, var, cut, shifts, wei);
  }

ana::Spectrum::Spectrum	(	const std::string &	label,
		SpectrumLoaderBase &	loader,
		const Binning &	binsx,
		const NuTruthVar &	varx,
		const Binning &	binsy,
		const NuTruthVar &	vary,
		const NuTruthCut &	cut,
		const SystShifts &	shifts = kNoShift,
		const NuTruthVar &	wei = kNuTruthUnweighted
	)

Definition at line 83 of file Spectrum.cxx.

References cut, label, demo0::loader, and ana::Var2D().

    : Spectrum({label, ""}, {binsx, binsy}) // TODO why no labelY?
  {
    loader.AddSpectrum(*this, Var2D(varx, binsx, vary, binsy), cut, shifts, wei);
  }

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

Definition at line 95 of file Spectrum.cxx.

References ana::SpectrumLoaderBase::AddSpectrum(), cut, ana::GenericHistAxis< T >::GetMultiDVar(), and plotSysts::shift.

    : Spectrum(axis.GetLabels(), axis.GetBinnings())
  {
    loader.AddSpectrum(*this, axis.GetMultiDVar(), cut, shift, wei);
  }

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

Definition at line 106 of file Spectrum.cxx.

References ana::SpectrumLoaderBase::AddSpectrum(), cut, ana::GenericHistAxis< T >::GetMultiDVar(), and plotSysts::shift.

    : Spectrum(axis.GetLabels(), axis.GetBinnings())
  {
    loader.AddSpectrum(*this, axis.GetMultiDVar(), cut, shift, wei);
  }

ana::Spectrum::Spectrum	(	SpectrumLoaderBase &	loader,
		const NuTruthHistAxis &	axis,
		const NuTruthCut &	cut,
		const SystShifts &	shifts = kNoShift,
		const NuTruthVar &	weight = kNuTruthUnweighted
	)

Definition at line 117 of file Spectrum.cxx.

References ana::SpectrumLoaderBase::AddSpectrum(), cut, ana::GenericHistAxis< T >::GetMultiDVar(), and ana::weight.

    : Spectrum(axis.GetLabels(), axis.GetBinnings())
  {
    loader.AddSpectrum(*this, axis.GetMultiDVar(), cut, shifts, weight);
  }

ana::Spectrum::Spectrum	(	const std::string &	label,
		const Binning &	bins,
		ESparse	sparse = kDense
	)

Definition at line 37 of file Spectrum.cxx.

References ConstructHistogram().

    : fHistD(0), fHistF(0), fHistSparse(0), fPOT(0), fLivetime(0),
      fLabels(1, label), fBins(1, bins)
  {
    ConstructHistogram(sparse);
  }

ana::Spectrum::Spectrum	(	const std::string &	label,
		double	pot,
		double	livetime,
		const Binning &	bins
	)

Definition at line 128 of file Spectrum.cxx.

References fLivetime, fPOT, livetime, and pot.

    : Spectrum(label, bins)
  {
    fPOT = pot;
    fLivetime = livetime;
  }

ana::Spectrum::Spectrum	(	TH1 *	h,
		const std::vector< std::string > &	labels,
		const std::vector< Binning > &	bins,
		double	pot,
		double	livetime
	)

Copies h.

Definition at line 137 of file Spectrum.cxx.

References ana::HistCache::Copy(), fHistD, fHistF, GetMultiDBinning(), and ana::HistCache::New().

    : fHistD(0), fHistF(0), fHistSparse(0), fPOT(pot), fLivetime(livetime), fLabels(labels), fBins(bins)
  {
    if(!h){
      fHistD = 0;
      fHistF = 0;
      return;
    }

    DontAddDirectory guard;

    const TString className = h->ClassName();

    if(className == "TH1D"){
      // Shortcut if types match
      fHistD = HistCache::Copy((TH1D*)h, GetMultiDBinning());
    }
    else if(className == "TH1F"){
      fHistF = HistCache::Copy((TH1F*)h, GetMultiDBinning());
    }
    else{
      fHistD = HistCache::New("", GetMultiDBinning());
      fHistD->Add(h);
    }
  }

ana::Spectrum::Spectrum	(	TH1 *	h,
		double	pot,
		double	livetime
	)

Definition at line 167 of file Spectrum.cxx.

    : Spectrum(h,
               std::vector<std::string>(1, h->GetXaxis()->GetTitle()),
               std::vector<Binning>(1, Binning::FromTAxis(h->GetXaxis())),
               pot, livetime)
  {
  }

ana::Spectrum::Spectrum	(	std::unique_ptr< TH1D >	h,
		const std::vector< std::string > &	labels,
		const std::vector< Binning > &	bins,
		double	pot,
		double	livetime
	)

Takes possession of h.

Definition at line 176 of file Spectrum.cxx.

References fHistD.

    : fHistD(h.release()), fHistF(0), fHistSparse(0), fPOT(pot), fLivetime(livetime), fLabels(labels), fBins(bins)
  {
    fHistD->SetDirectory(0);
  }

ana::Spectrum::Spectrum	(	std::unique_ptr< TH1F >	h,
		const std::vector< std::string > &	labels,
		const std::vector< Binning > &	bins,
		double	pot,
		double	livetime
	)

Definition at line 186 of file Spectrum.cxx.

References fHistF.

    : fHistD(0), fHistF(h.release()), fHistSparse(0), fPOT(pot), fLivetime(livetime), fLabels(labels), fBins(bins)
  {
    fHistF->SetDirectory(0);
  }

ana::Spectrum::Spectrum	(	const std::string &	label,
		SpectrumLoaderBase &	loader,
		const Binning &	binsx,
		const Var &	varx,
		const Binning &	binsy,
		const Var &	vary,
		const Cut &	cut,
		const SystShifts &	shift = kNoShift,
		const Var &	wei = kUnweighted
	)

2D Spectrum of two Vars

Definition at line 196 of file Spectrum.cxx.

    : Spectrum(label, "", loader, binsx, varx, binsy, vary, cut, shift, wei)
  {
    // TODO do we want this variant when there's one with a labelY just below?
  }

ana::Spectrum::Spectrum	(	const std::string &	label,
		SpectrumLoaderBase &	loader,
		const Binning &	binsx,
		const MultiVar &	varx,
		const Binning &	binsy,
		const MultiVar &	vary,
		const Cut &	cut,
		const SystShifts &	shift = kNoShift,
		const Var &	wei = kUnweighted
	)

2D Spectrum of two MultiVars

Definition at line 208 of file Spectrum.cxx.

    : Spectrum(label, "", loader, binsx, varx, binsy, vary, cut, shift, wei)
  {
  }

ana::Spectrum::Spectrum	(	SpectrumLoaderBase &	loader,
		const HistAxis &	xAxis,
		const HistAxis &	yAxis,
		const Cut &	cut,
		const SystShifts &	shift = kNoShift,
		const Var &	wei = kUnweighted
	)

2D Spectrum taking 2 HistAxis

Definition at line 219 of file Spectrum.cxx.

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

    : Spectrum(xAxis.GetLabels()[0], yAxis.GetLabels()[0],
               loader,
               xAxis.GetBinnings()[0], xAxis.GetVars()[0],
               yAxis.GetBinnings()[0], yAxis.GetVars()[0],
               cut, shift, wei)
  {
    // TODO - do we want to keep this variant around?
    assert(xAxis.NDimensions() == 1);
    assert(yAxis.NDimensions() == 1);
  }

ana::Spectrum::Spectrum	(	const std::string &	xLabel,
		const std::string &	yLabel,
		SpectrumLoaderBase &	loader,
		const Binning &	binsx,
		const Var &	varx,
		const Binning &	binsy,
		const Var &	vary,
		const Cut &	cut,
		const SystShifts &	shift = kNoShift,
		const Var &	wei = kUnweighted
	)

Definition at line 237 of file Spectrum.cxx.

References cut, demo0::loader, plotSysts::shift, and ana::Var2D().

    : Spectrum({xLabel, yLabel}, {binsx, binsy})
  {
    Var multiDVar = Var2D(varx, binsx, vary, binsy);

    loader.AddSpectrum(*this, multiDVar, cut, shift, wei);
  }

ana::Spectrum::Spectrum	(	const std::string &	xLabel,
		const std::string &	yLabel,
		SpectrumLoaderBase &	loader,
		const Binning &	binsx,
		const MultiVar &	varx,
		const Binning &	binsy,
		const MultiVar &	vary,
		const Cut &	cut,
		const SystShifts &	shift = kNoShift,
		const Var &	wei = kUnweighted
	)

Definition at line 253 of file Spectrum.cxx.

References cut, demo0::loader, ana::MultiVar2D(), and plotSysts::shift.

    : Spectrum({xLabel, yLabel}, {binsx, binsy})
  {
    MultiVar multiDVar = MultiVar2D(varx, binsx, vary, binsy);
    
    loader.AddSpectrum(*this, multiDVar, cut, shift, wei);
  }

ana::Spectrum::Spectrum	(	const std::string &	label,
		SpectrumLoaderBase &	loader,
		const Binning &	binsx,
		const Var &	varx,
		const Binning &	binsy,
		const Var &	vary,
		const Binning &	binsz,
		const Var &	varz,
		const Cut &	cut,
		const SystShifts &	shift = kNoShift,
		const Var &	wei = kUnweighted,
		ESparse	sparse = kDense
	)

3D Spectrum of three Vars

Definition at line 269 of file Spectrum.cxx.

    : Spectrum(label, "", "", loader, binsx, varx, binsy, vary, binsz, varz, cut, shift, wei, sparse)
  {
    // TODO do we want this variant when there's one with a labelY and labelZ
    // just below?
  }

ana::Spectrum::Spectrum	(	const std::string &	xLabel,
		const std::string &	yLabel,
		const std::string &	zLabel,
		SpectrumLoaderBase &	loader,
		const Binning &	binsx,
		const Var &	varx,
		const Binning &	binsy,
		const Var &	vary,
		const Binning &	binsz,
		const Var &	varz,
		const Cut &	cut,
		const SystShifts &	shift = kNoShift,
		const Var &	wei = kUnweighted,
		ESparse	sparse = kDense
	)

Definition at line 284 of file Spectrum.cxx.

References cut, demo0::loader, plotSysts::shift, and ana::Var3D().

    : Spectrum({xLabel, yLabel, zLabel}, {binsx, binsy, binsz}, sparse)
  {
    Var multiDVar = Var3D(varx, binsx, vary, binsy, varz, binsz);

    loader.AddSpectrum(*this, multiDVar, cut, shift, wei);
  }

ana::Spectrum::Spectrum	(	SpectrumLoaderBase &	loader,
		const HistAxis &	xAxis,
		const HistAxis &	yAxis,
		const HistAxis &	zAxis,
		const Cut &	cut,
		const SystShifts &	shift = kNoShift,
		const Var &	wei = kUnweighted,
		ESparse	sparse = kDense
	)

3D Spectrum taking 3 HistAxis

Definition at line 303 of file Spectrum.cxx.

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

    : Spectrum(xAxis.GetLabels()[0], loader,
               xAxis.GetBinnings()[0], xAxis.GetVars()[0],
               yAxis.GetBinnings()[0], yAxis.GetVars()[0],
               zAxis.GetBinnings()[0], zAxis.GetVars()[0],
               cut, shift, wei, sparse)
  {
    // TODO - do we want to keep this variant around?
    assert(xAxis.NDimensions() == 1);
    assert(yAxis.NDimensions() == 1);
    assert(zAxis.NDimensions() == 1);
  }

ana::Spectrum::~Spectrum ( )

virtual

Definition at line 324 of file Spectrum.cxx.

References om::cerr, ana::HistCache::Delete(), allTimeWatchdog::endl, fHistD, fHistF, fHistSparse, fLoaderCount, GetMultiDBinning(), and demo0::loader.

Referenced by operator=().

  {
    if(fHistD && fHistD->GetDirectory()){
      static bool once = true;
      if(once){
        once = false;
        std::cerr << "Spectrum's fHist (" << fHistD << ") is associated with a directory (" << fHistD->GetDirectory() << " " << fHistD->GetDirectory()->GetName() << " " << fHistD->GetDirectory()->GetTitle() << "). How did that happen?" << std::endl;
      }
    }

    for (SpectrumLoaderBase* loader : fLoaderCount)
    { loader->RemoveSpectrum(this); }

    if(fHistD) HistCache::Delete(fHistD, GetMultiDBinning());
    if(fHistF) HistCache::Delete(fHistF, GetMultiDBinning());

    delete fHistSparse;
  }

ana::Spectrum::Spectrum

(

const Spectrum &

rhs

)

Definition at line 344 of file Spectrum.cxx.

References ana::assert(), ana::HistCache::Copy(), fHistD, fHistF, fHistSparse, fLoaderCount, and GetMultiDBinning().

                                       :
    fHistD(0),
    fHistF(0),
    fHistSparse(0),
    fPOT(rhs.fPOT),
    fLivetime(rhs.fLivetime),
    fLabels(rhs.fLabels),
    fBins(rhs.fBins)
  {
    DontAddDirectory guard;

    assert(rhs.fHistD || rhs.fHistF || rhs.fHistSparse);
    if(rhs.fHistD)
      fHistD = HistCache::Copy(rhs.fHistD, rhs.GetMultiDBinning());
    if(rhs.fHistF)
      fHistF = HistCache::Copy(rhs.fHistF, rhs.GetMultiDBinning());
    if(rhs.fHistSparse){
      // Doesn't exist?
      // fHistSparse = new THnSparseD(*rhs.fHistSparse);
      fHistSparse = (THnSparseD*)rhs.fHistSparse->Clone();
    }

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

ana::Spectrum::Spectrum

(

Spectrum &&

rhs

)

Definition at line 370 of file Spectrum.cxx.

References ana::assert(), fHistD, fHistF, fHistSparse, fLoaderCount, and demo0::loader.

                                  :
    fHistD(0),
    fHistF(0),
    fHistSparse(0),
    fPOT(rhs.fPOT),
    fLivetime(rhs.fLivetime),
    fLabels(rhs.fLabels),
    fBins(rhs.fBins)
  {
    assert(rhs.fHistD || rhs.fHistF || rhs.fHistSparse);

    if(rhs.fHistD){
      fHistD = rhs.fHistD;
      rhs.fHistD = 0;
    }
    if(rhs.fHistF){
      fHistF = rhs.fHistF;
      rhs.fHistF = 0;
    }
    if(rhs.fHistSparse){
      fHistSparse = rhs.fHistSparse;
      rhs.fHistSparse = 0;
    }

    for (auto & loader : rhs.fLoaderCount)
    {
      loader->ReplaceSpectrum(&rhs, this);
      fLoaderCount.insert(loader);
    }
    rhs.fLoaderCount.clear();

  }

ana::Spectrum::Spectrum

(

const SpectrumStan &

rhs

)

Copy conversion from SpectrumStan.

Definition at line 404 of file Spectrum.cxx.

References ana::assert(), ConstructHistogram(), fHistD, fHistF, ana::SpectrumStan::GetBinnings(), ana::gFloatMode, ana::SpectrumStan::POT(), ana::SpectrumStan::ToBins(), and febshutoff_auto::val.

    : fHistD(nullptr), fHistF(nullptr), fHistSparse(nullptr),
      fPOT(rhs.POT()), fLivetime(rhs.Livetime()),
      fLabels(rhs.GetLabels()), fBins(rhs.GetBinnings())
  {
    // for now we don't have any other use cases besides 1D...
    assert(rhs.GetBinnings().size() == 1);

    ConstructHistogram();

    // now fill in the bin content...
    auto binC = rhs.ToBins(rhs.POT());
    for (int binIdx = 0; binIdx < rhs.GetBinnings()[0].NBins()+2; binIdx++) // be sure to get under- and overflow
    {
      if (gFloatMode)
        fHistF->SetBinContent(binIdx, binC[binIdx].val());
      else
        fHistD->SetBinContent(binIdx, binC[binIdx].val());
    }
  }

ana::Spectrum::Spectrum ( const std::vector< std::string > &  labels, const std::vector< Binning > &  bins, ESparse  sparse = kDense  )

protected

Definition at line 26 of file Spectrum.cxx.

References ConstructHistogram().

    : fHistD(0), fHistF(0), fHistSparse(0), fPOT(0), fLivetime(0),
      fLabels(labels), fBins(bins)
  {
    ConstructHistogram(sparse);
  }

Member Function Documentation

void ana::Spectrum::AddLoader ( SpectrumLoaderBase *  p )

protected

Definition at line 901 of file Spectrum.cxx.

References fLoaderCount.

Referenced by ana::SpectrumLoaderBase::AddSpectrum(), and GetBinnings().

  { fLoaderCount.insert(p); }

void ana::Spectrum::Clear

(

)

Definition at line 889 of file Spectrum.cxx.

References fHistD, fHistF, and fHistSparse.

Referenced by ana::PredictionInterp::_PredictComponentSyst(), ana::PredictionExtrapSum::ComponentNC(), ana::PredictionExtrapSum::ComponentNCAnti(), ana::PredictionExtrapSum::ComponentNCTotal(), GetNueCosmics2018(), GetNueCosmics2019(), GetNueCosmicsFuture(), ana::PredictionInterp::InitFits(), make_plots(), OverrideLivetime(), plot_nueFD_Signal_REWvsNOM_pTExtrap_FHC(), plot_nueFD_Signal_REWvsNOM_pTExtrap_RHC(), ana::PlotSpectra(), PlotStack(), ana::PredictionTemplateShift::PredictComponent(), ana::PredictionSterile::PredictComponent(), ana::PredictionNueRebinSA::PredictComponent(), ana::FDPredictionSterile::PredictComponent(), ana::PredictionExtrapSum::PredictComponent(), ana::NDPredictionSterile::PredictComponent(), ana::PredictionNoOsc::PredictComponent(), and ana::PredictionExtendToPeripheral::PredictComponent().

  {
    if(fHistD) fHistD->Reset();
    if(fHistF) fHistF->Reset();
    if(fHistSparse) fHistSparse->Reset();
  }

void ana::Spectrum::ConstructHistogram ( ESparse  sparse = kDense )

protected

Definition at line 524 of file Spectrum.cxx.

References ana::assert(), fHistD, fHistF, fHistSparse, GetMultiDBinning(), ana::gFloatMode, ana::Binning::IsSimple(), kDense, kFloat, kSparse, ana::Binning::Max(), ana::Binning::Min(), nbins, ana::Binning::NBins(), ana::HistCache::New(), ana::HistCache::NewF(), ana::UniqueName(), make_syst_table_plots::xmax, and make_mec_shifts_plots::xmin.

Referenced by GetBinnings(), and Spectrum().

  {
    DontAddDirectory guard;

    assert(!fHistD && !fHistF && !fHistSparse);

    const Binning bins1D = GetMultiDBinning();

    if(sparse == kSparse){
      assert(bins1D.IsSimple());
      const int nbins = bins1D.NBins();
      const double xmin = bins1D.Min();
      const double xmax = bins1D.Max();
      fHistSparse = new THnSparseD(UniqueName().c_str(), UniqueName().c_str(),
                                   1, &nbins, &xmin, &xmax);

      // Ensure errors get accumulated properly
      fHistSparse->Sumw2();
    }
    else if(sparse == kDense && !gFloatMode){
      fHistD = HistCache::New("", bins1D);
      fHistF = 0;

      // first clear the sum-of-weights.
      // can be pre-filled with 0s if you happen to have made several Spectra
      // out of the same Binning and a preceding one is still in the cache.
      fHistD->Sumw2(false);
      fHistD->Sumw2();      // Ensure errors get accumulated properly
    }
    else if(sparse == kFloat || gFloatMode){
      fHistF = HistCache::NewF("", bins1D);
      fHistD = 0;

      // see above
      fHistF->Sumw2(false);
      fHistF->Sumw2();      // Ensure errors get accumulated properly
    }
  }

Spectrum ana::Spectrum::FakeData

(

double

pot

)

const

Fake data is a MC spectrum scaled to the POT expected in the data.

Use for sensitivity plots and testing fit convergence

Definition at line 830 of file Spectrum.cxx.

References fHistD, fHistF, fHistSparse, fPOT, pot, and runNovaSAM::ret.

Referenced by BlessedPlotsLoad(), caf_numu_sensitivity(), caf_numu_sensitivity_no_tau(), cc(), compare_pid_cut_sens(), contours(), cpv(), delta_precision(), delta_slices(), demo_CPT(), extrap_syst_fa(), fa_numuContours(), ana::FDPredictionSterile::FakeData(), ana::NDPredictionSterile::FakeData(), futureSig_reach_singlePOTcombo_syst(), GenerateFutureData(), ana2019::fakedata::Get2019Prediction(), get_sa_syst_mags(), get_sa_systs_mags_v2(), getContProf_Sensitivity(), GetFakeData(), getSensitivity(), hierarchy(), joint_fit_make_fc_slice(), joint_fit_make_fc_surf(), make_fa_pred(), make_plots(), ana::NumuAnalysis::MakeExperiment(), MakeExtrapSurface(), MakeFakeData(), MakeNusPrediction(), makeSystTable(), makeSystTable_reduced(), MichelDecompTest(), MockData(), NDDataMCSystBandAna(), nue_ana_basic(), numu_demo_split_plot(), numu_sig_nonmax(), octant(), plot_extrap_syst_fa(), plot_fa_results(), plot_sa_results_nodata(), PlotNus17Prediction(), plots(), saTestContours(), saveContours_addExpt(), saveContours_complete(), saveContours_oscpar(), saveContours_simple(), saveContours_systs(), sensitivity(), sensitivity_slices(), starPlot(), syst_test(), template_basic(), template_GENIE_systs(), template_nonGENIE_systs(), test_ana(), test_stanfit_statsonly(), test_stanfit_systpulls(), test_stanfit_withsysts(), th23_delta_contour(), ToTH1(), and ana::UnfoldIterative::Truth().

  {
    Spectrum ret = *this;
    if(fPOT > 0){
      if(ret.fHistD)
        ret.fHistD->Scale(pot/fPOT);
      if(ret.fHistF)
        ret.fHistF->Scale(pot/fPOT);
      if(ret.fHistSparse)
        ret.fHistSparse->Scale(pot/fPOT);
    }
    ret.fPOT = pot;

    // Drop old errors, which are based on the MC statistics, and create new
    // ones that are based on the prediction for the data
    if(ret.fHistD){
      ret.fHistD->Sumw2(false);
      ret.fHistD->Sumw2();
    }
    if(ret.fHistF){
      ret.fHistF->Sumw2(false);
      ret.fHistF->Sumw2();
    }

    return ret;
  }

Spectrum ana::Spectrum::FakeData	(	double	pot,
		double	livetime
	)		const

Use for sensitivity plots when fake cosmic data is needed. Fake cosmic spectra can be added to FakeData by desired livetime.

Definition at line 858 of file Spectrum.cxx.

References fHistD, fHistF, fHistSparse, fLivetime, fPOT, livetime, pot, and runNovaSAM::ret.

  {
    Spectrum ret = *this;

    if(fPOT > 0){
      if(ret.fHistD)
        ret.fHistD->Scale(pot/fPOT);
      if(ret.fHistF)
        ret.fHistF->Scale(pot/fPOT);
      if(ret.fHistSparse)
        ret.fHistSparse->Scale(pot/fPOT);
    }
    ret.fPOT = pot;

    // overwrite livetime for fake data
    // PlusEqualsHelper will take over the rest
    ret.fLivetime = livetime;

    if(ret.fHistD){
      ret.fHistD->Sumw2(false);
      ret.fHistD->Sumw2();
    }
    if(ret.fHistF){
      ret.fHistF->Sumw2(false);
      ret.fHistF->Sumw2();
    }

    return ret;
  }

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

Definition at line 785 of file Spectrum.cxx.

References ana::assert(), fHistD, fHistF, and fHistSparse.

Referenced by first_ana_proj(), and ana::SpectrumLoader::HandleRecord().

  {
    assert( (fHistD || fHistF || fHistSparse) && "Somehow both fHist and fHistSparse are null in Spectrum::Fill" );

    if(fHistD)
      fHistD->Fill(x, w);
    else if(fHistF)
      fHistF->Fill(x, w);
    else if(fHistSparse)
      fHistSparse->Fill(&x, w);
  }

std::vector<Binning> ana::Spectrum::GetBinnings ( ) const

inline

Definition at line 299 of file Spectrum.h.

References AddLoader(), ConstructHistogram(), fBins, GetMultiDBinning(), kDense, PlusEqualsHelper(), RemoveLoader(), canMan::sign(), and Spectrum().

Referenced by ana::FitInAnaBinsBkgdEstimator::Background(), ana::NumuCC2p2hBkgdEstimator::Background(), ana::CovMxManager::CovMxManager(), ana::TwoSampleDecomp::Decomp(), ana::FakeDecomp::Decompose(), ana::PredictionInterp::FitSumComponent(), ana::SpectrumComponents::GetBinnings(), ana::Multiverse::GetNSigmaShift(), ana::nueccinc::NueCCIncCrossSectionAnalysis::getTemplateShapeOnly1D(), ana::nueccinc::NueCCIncCrossSectionAnalysis::getTemplateShapeOnly3D(), ana::PredictionNueRebinSA::Inflate(), ana::NumuCCIncAnalysis::LoadFrom(), ana::FluxDecomp::NueEstimate(), ana::BENDecomp::NueEstimate(), ana::FluxDecomp::NueEstimateFromKa(), ana::BENDecomp::NueEstimateFromKa(), ana::FluxDecomp::NueEstimateFromPi(), ana::BENDecomp::NueEstimateFromPi(), ana::PredictionCombinePeriods::PredictComponentSyst(), ana::PredictionNumuFAHadE::PredictComponentSyst(), ana::PredictionNueRebinSA::Reduce(), ana::PredictionExtendToPeripheral::ReduceHelperNC(), ana::PredictionInterp::ShiftSpectrum(), ana::SpectrumStan::SpectrumStan(), ana::UnfoldSVD::Truth(), ana::UnfoldTikhonov::Truth(), and ana::PredictionTemplateShift::Weight().

{return fBins;}

std::vector<std::string> ana::Spectrum::GetLabels ( ) const

inline

Definition at line 298 of file Spectrum.h.

References fLabels.

Referenced by ana::FitInAnaBinsBkgdEstimator::Background(), ana::NumuCC2p2hBkgdEstimator::Background(), ana::TwoSampleDecomp::Decomp(), ana::FakeDecomp::Decompose(), ana::SpectrumComponents::GetLabels(), ana::Multiverse::GetNSigmaShift(), ana::nueccinc::NueCCIncCrossSectionAnalysis::getTemplateShapeOnly1D(), ana::nueccinc::NueCCIncCrossSectionAnalysis::getTemplateShapeOnly3D(), ana::PredictionNueRebinSA::Inflate(), ana::PredictionExtendToPeripheral::MergePeripheral(), ana::FluxDecomp::NueEstimate(), ana::BENDecomp::NueEstimate(), ana::FluxDecomp::NueEstimateFromKa(), ana::BENDecomp::NueEstimateFromKa(), ana::FluxDecomp::NueEstimateFromPi(), ana::BENDecomp::NueEstimateFromPi(), ana::CutOptimization::OptimizedSigmaOverSigma(), ana::PredictionCombinePeriods::PredictComponentSyst(), ana::PredictionNumuFAHadE::PredictComponentSyst(), ana::PredictionNueRebinSA::Reduce(), ana::PredictionExtendToPeripheral::ReduceHelperNC(), ana::PredictionInterp::ShiftSpectrum(), ana::UnfoldSVD::Truth(), ana::UnfoldTikhonov::Truth(), and ana::PredictionTemplateShift::Weight().

{return fLabels;}

Binning ana::Spectrum::GetMultiDBinning ( ) const

protected

Definition at line 510 of file Spectrum.cxx.

References ana::assert(), b, fBins, getGoodRuns4SAM::n, and ana::Binning::Simple().

Referenced by ConstructHistogram(), GetBinnings(), operator=(), Spectrum(), ToTH1(), and ~Spectrum().

  {
    // TODO same function exists in HistAxis

    assert(!fBins.empty());

    if(fBins.size() == 1) return fBins[0];

    int n = 1;
    for(const Binning& b: fBins) n *= b.NBins();
    return Binning::Simple(n, 0, n);
  }

double ana::Spectrum::Integral	(	double	exposure,
		double *	err = 0,
		EExposureType	expotype = kPOT
	)		const

Return total number of events scaled to pot.

Parameters

exposure	POT/livetime to scale to
err	The statistical error on this total (optional)
expotype	What the first parameter represents

Definition at line 753 of file Spectrum.cxx.

References fHistD, fHistF, fLivetime, fPOT, make_syst_table_plots::h, MECModelEnuComparisons::i, ana::kPOT, plotSysts::ratio, util::sqr(), and std::sqrt().

Referenced by BicountEllipse_dCP(), ana::CountingExperiment::ChiSq(), compare_pid_cut_sens(), ComputeEfficiency(), contours(), ana::DataMCAreaNormalizedRatio(), ana::SpectrumComponents::DrawLegend(), ana::DataMCPair::DrawMCNormSyst(), DataMCPair::DrawMCNormSyst(), Efficiency(), FCContour(), ana::FluxMultiverseSyst::FindSigmaBoundaries(), FitParamEffectsAna(), FitSystEffectsAna(), getHists_FNEX(), GetNueCosmics2017(), GetNueCosmics2018(), GetNueCosmics2019(), ana::GetNueCosmics2020(), GetNueCosmicsFuture(), GetNueData2017(), GetNueData2018(), GetNueData2019(), ana::GetNueData2020(), GetNuePrediction2017(), GetNuePrediction2018(), GetNuePrediction2019(), ana::GetNuePrediction2020(), GetNuePredictionFuture(), ana::GetNumuCosmics2020(), ana::GetNumuData2020(), GetSpectra(), ana::nueccinc::NueCCIncCrossSectionAnalysis::getTemplateShapeOnly1D(), ana::nueccinc::NueCCIncCrossSectionAnalysis::getTemplateShapeOnly3D(), Integral(), joint_fit_2017_make_fc_slice(), joint_fit_2017_make_fc_surf(), joint_fit_make_experiments(), joint_fit_make_fc_slice(), joint_fit_make_fc_surf(), make_fc_mass_and_oct_nersc_2018(), make_fc_mass_and_oct_nersc_2019(), make_fc_mh_nersc_2018(), make_fc_mh_nersc_2019(), make_fc_oct_nersc_2018(), make_fc_oct_nersc_2019(), make_fc_slice(), make_fc_slices_nersc_2018(), make_fc_slices_nersc_2018_stats(), make_fc_slices_nersc_2019(), make_fc_surf(), make_fc_surfaces_nersc_2018(), make_fc_surfaces_nersc_2018_stats(), make_fc_surfaces_nersc_2019(), make_nominal_xs(), make_xs(), MakeHBar(), MakeTable(), MCCompPredictionTable(), monoprob_plot(), ana::GenieMultiverseNormalizedSpectra::NormalizeSpectra(), PeripheralCuts(), ana::PlotWithAreaSystErrorBand(), predEventCountWithSystError(), PrintCount(), PrintEventCounts(), sensitivity(), sensitivity_slices(), slices(), slices_paper(), SpectrumParamEffectsAna(), ToTH1(), and ana::UnfoldSVD::Truth().

  {
    const double ratio = (expotype == kPOT) ? exposure/fPOT : exposure/fLivetime;

    TH1* h = fHistD ? (TH1*)fHistD : (TH1*)fHistF;

    if(err){
      *err = 0;

      for(int i = 0; i < h->GetNbinsX()+2; ++i){
        *err += util::sqr(h->GetBinError(i));
      }
      *err = sqrt(*err) * ratio;
    }

    // TODO how to integrate fHistSparse?

    return h->Integral(0, -1) * ratio;
  }

double ana::Spectrum::Livetime ( ) const

inline

Seconds. For informational purposes only. No calculations use this.

Definition at line 266 of file Spectrum.h.

References fLivetime.

Referenced by ana::FitInAnaBinsBkgdEstimator::Background(), ana::NumuCC2p2hBkgdEstimator::Background(), BlessedPlotsLoad(), CalcRWithSystsNus17(), contours(), Cumulative(), ana::TwoSampleDecomp::Decomp(), ana::FakeDecomp::Decompose(), FCContour(), fd_plot(), fd_plot_bless(), ana::FluxMultiverseSyst::FindSigmaBoundaries(), FitParamEffectsAna(), FitSystEffectsAna(), getBNBPlots(), GetFluxError(), getHists_FNEX(), ana::Multiverse::GetNSigmaShift(), GetNueCosmics2017(), GetNueCosmics2018(), GetNueCosmics2019(), ana::GetNueCosmics2020(), GetNueCosmicsFuture(), GetNueData2017(), GetNueData2018(), GetNueData2019(), ana::GetNueData2020(), ana::GetNumuCosmics2020(), ana::GetNumuData2020(), ana::PredictionInterp::InitFits(), joint_fit_make_experiments(), joint_fit_make_fc_slice(), joint_fit_make_fc_surf(), make_fc_nus_surfs_nersc_2018(), make_fc_nus_surfs_nersc_2019(), make_fc_points(), make_fc_slice(), make_fc_surf(), make_nus17_fc_surfs(), make_nus_fc_surfs(), make_xs(), make_xs_1D(), MakeExtrapSurface(), MakeHBar(), makeMatrixElementSurface(), MakeNusPrediction(), ana::MichelDecomp::MCToDCMPComp(), ana::PredictionExtendToPeripheral::MergePeripheral(), monoprob_plot(), mrbrem_get_reweighted_spectra(), mrbrem_plots(), plot_datamcpred(), PlotNus17Prediction(), PlotNus17PredSystsData(), PlotNus18Sideband(), PlotNusSensAna01(), ana::PlotSpectra(), PlotSysts(), ana::PredictionNumuFAHadE::PredictComponentSyst(), ana::PredictionExtendToPeripheral::ReduceHelperNC(), selection_story_plots(), sensitivity(), sensitivity_slices(), ana::PredictionInterp::ShiftSpectrum(), ana::SingleSampleExperiment::SingleSampleExperiment(), slices(), slices_paper(), timingPeak(), ana::UnfoldSVD::Truth(), ana::UnfoldTikhonov::Truth(), and ana::PredictionTemplateShift::Weight().

{return fLivetime;}

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

static

Definition at line 1066 of file Spectrum.cxx.

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

Referenced by Analyse_GetEfficiency(), Analyse_GetEfficiency_UseNEntries(), bdtstudyplotter(), BlessedPlotsAnaByPeriod(), CalcRWithSysts(), calculateComponentsNumu(), calculateWrongSignNue(), calculateWrongSignNumuQ1(), calculateWrongSignNumuQ2(), calculateWrongSignNumuQ3(), calculateWrongSignNumuQ4(), CalculateXSec(), ana::CompareNDDataOneMC(), ana::CompareNDDataTwoMC(), ana::ComparisonTable(), ana::CovMxManager::CovMxManager(), CutTableAna(), CVNphotonHist(), CVNphotonSplitHist(), demo2p5b(), DoThePlots(), doUnfolding(), EHadVisMECpairs(), FDDataMC(), FDDataMCSystBandAna(), FidOptHist(), FillFlavorContainers(), FitParamEffectsAna(), FitSystEffectsAna(), FOMCalc(), GeniePredictionToRoot(), get_data_histogram(), get_numu_data_histogram(), GetCosmics(), Plotter::GetDataPlots(), GetHist(), getHists_FNEX(), GetHistVectors(), Plotter::GetMCPlots(), ana::GetNDComponents(), ana::GetNumuCosmics2020(), GetNumuCosmicsFuture(), GetNumuData2017(), GetNumuData2018(), GetNumuData2019(), GetSpectToHist(), ana::InteractionSpectra::InteractionSpectra(), joint_fit_make_experiments(), Load1DHistFromSpec(), ana::CountingExperiment::LoadFrom(), jw::TrivialPrediction::LoadFrom(), ana::TrivialPrediction::LoadFrom(), ana::SingleSampleExperiment::LoadFrom(), ana::CrossSectionSpectra::LoadFrom(), ana::CutOptimization::LoadFrom(), ana::TrivialCrossSectionAnalysis::LoadFrom(), ana::Multiverse::LoadFrom(), ana::SpectrumComponents::LoadFrom(), ana::DataMCPair::LoadFrom(), ana::RecoReweight::LoadFrom(), ana::FluxReweight::LoadFrom(), ana::LoadFromUpDownSpectra(), LoadMaps(), make_dataMC(), make_eff_plots_areaNorm(), make_estimate_energy(), make_fc_points(), make_muonid_opt(), make_nominal_xs(), Make_NuMuCC_Inc_XS(), Make_NuMuCC_Inc_XS_v2(), make_vertex_optimiz(), MakeCutFlow(), makeEnergyEstimator(), makeFlatWeight(), MakePlots(), makeResolutionPlots(), MakeSurface(), MakeThePlots(), makeXSecPlots1D(), makeXSecPlots2D(), makeXSecPlots_TemplateFit(), mec_tuning_fitter_2020(), MRDiFStudy_FHC_Step2(), MRDiFStudy_RHC_Step2(), MuonCatcherComp_ProdPlots(), nc_bkgd_by_interaction_mode(), NDDataMC(), NDDataMCSystBandAna(), nue_bg_stat_err_calc(), nue_data_mc_validation(), nue_decomp_scales(), nue_fd_mc_validation(), numu_data_mc_validation(), OverrideLivetime(), plot_DataMCComp_numu(), plot_ND_DataMC(), plot_ND_DataMC_energybinning(), plot_nd_spectra_2018(), plot_paper_nd_datamc(), plot_recoE_numu(), plot_sa_results(), plot_spectra_2dplots(), plot_uncertainty(), plotHist_SpectrumCVNID(), PlotResolution(), Plotsidebandfittest(), Plotting_Data2DataComp(), Plotting_Data2DataComp_SingleCanvas(), Plotting_DataAndPrediction(), Plotting_DataSpectra_MakeTGraph(), Plotting_OverlayStudies(), ppfx_smooth_weights_make(), preselection_cutflow(), ProduceCompPlots(), ReMIdHist(), resolutionplotter(), ana::ResolutionScan::ResolutionScan(), SaveHistogramsFD(), SaveHistogramsND(), ShwZOptHist(), Spec2DToHist(), Spec2DtoHist(), SpecToHist(), SpectrumParamEffectsAna(), twodvtxcontplotter(), Unfold1D(), Unfold3D(), vertexstudyploter(), xsec_tot_uncert_optimization(), and xsec_uncertainty_per_bin().

  {
    DontAddDirectory guard;

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

    TH1D* spectD = 0;
    TH1F* spectF = 0;
    dir->GetObject("hist", spectD);
    dir->GetObject("hist", spectF);
    THnSparseD* spectSparse = (THnSparseD*)dir->Get("hist_sparse");
    assert(spectD || spectF || spectSparse);
    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){
      TDirectory* subdir = dir->GetDirectory(TString::Format("bins%d", i));
      if(!subdir) break;
      bins.push_back(*Binning::LoadFrom(subdir));
      TObjString* label = (TObjString*)dir->Get(TString::Format("label%d", i));
      labels.push_back(label ? label->GetString().Data() : "");
      delete subdir;
      delete label;
    }

    if(bins.empty() && labels.empty()){
      // Must be an old file. Make an attempt at backwards compatibility.
      if(spectD){
        bins.push_back(Binning::FromTAxis(spectD->GetXaxis()));
        labels.push_back(spectD->GetXaxis()->GetTitle());
      }
      else{
        bins.push_back(Binning::FromTAxis(spectSparse->GetAxis(0)));
        labels.push_back(spectSparse->GetAxis(0)->GetTitle());
      }
    }

    if(spectSparse){
      std::unique_ptr<Spectrum> ret = std::make_unique<Spectrum>((TH1*)0, labels, bins, hPot->GetBinContent(1), hLivetime->GetBinContent(1));
      ret->fHistSparse = spectSparse;

      delete spectD;
      delete spectF;
      delete hPot;
      delete hLivetime;
      return ret;
    }
    else{
      std::unique_ptr<Spectrum> ret;
      if(spectD)
        ret = std::make_unique<Spectrum>(std::unique_ptr<TH1D>(spectD), labels, bins, hPot->GetBinContent(1), hLivetime->GetBinContent(1));
      else
        ret = std::make_unique<Spectrum>(std::unique_ptr<TH1F>(spectF), labels, bins, hPot->GetBinContent(1), hLivetime->GetBinContent(1));
      delete spectSparse;
      delete hPot;
      delete hLivetime;
      return ret;
    }
  }

double ana::Spectrum::Mean

(

)

const

Return mean of 1D histogram.

Definition at line 775 of file Spectrum.cxx.

References fHistD, fHistF, and make_syst_table_plots::h.

Referenced by ana::DataMCPair::CreateSystTable(), and ToTH1().

  {
    // Allow GetMean() to work even if this histogram never had any explicit
    // Fill() calls made.
    TH1* h = fHistD ? (TH1*)fHistD : (TH1*)fHistF;
    if(h->GetEntries() == 0) h->SetEntries(1);
    return h->GetMean();
  }

Spectrum ana::Spectrum::MockData	(	double	pot,
		int	idx = 0
	)		const

Mock data is FakeData with Poisson fluctuations applied.

Use for low-budget MDCs, or just getting a sense of the expected scale of statistical variation

Definition at line 798 of file Spectrum.cxx.

References FakeData(), fHistD, fHistF, fHistSparse, make_syst_table_plots::h, MECModelEnuComparisons::i, runNovaSAM::ret, and generate_hists::rnd.

Referenced by demo4(), demo5(), FCTutorial2020(), fd_plot(), fill_col(), ana2019::fakedata::Get2019Prediction(), GetMockData(), goodness_of_fit(), joint_fit_2017_make_fc_slice(), joint_fit_2017_make_fc_surf(), joint_fit_make_fc_slice(), joint_fit_make_fc_surf(), make_fc_mass_and_oct_nersc_2018(), make_fc_mass_and_oct_nersc_2019(), make_fc_mh_nersc_2018(), make_fc_mh_nersc_2019(), make_fc_nus_surfs_nersc_2018(), make_fc_nus_surfs_nersc_2019(), make_fc_oct_nersc_2018(), make_fc_oct_nersc_2019(), make_fc_points(), make_fc_slice(), make_fc_slices_nersc_2018(), make_fc_slices_nersc_2018_stats(), make_fc_slices_nersc_2019(), make_fc_surf(), make_fc_surfaces_2020_validation(), make_fc_surfaces_nersc_2018(), make_fc_surfaces_nersc_2018_stats(), make_fc_surfaces_nersc_2019(), make_mockdata_syst_contours(), make_nus17_fc_surfs(), make_nus_fc_surfs(), selection_story_plots(), syst_plot_test(), test_saloaders(), and ToTH1().

  {
    Spectrum ret = FakeData(pot);

    TRandom3 rnd(idx); // zero seeds randomly

    if(ret.fHistD || ret.fHistF){
      TH1* h = ret.fHistD ? (TH1*)ret.fHistD : (TH1*)ret.fHistF;
      for(int i = 0; i < h->GetNbinsX()+2; ++i){
        h->SetBinContent(i, rnd.Poisson(h->GetBinContent(i)));
      }
    }
    if(ret.fHistSparse){
      for(int i = 0; i < ret.fHistSparse->GetNbins(); ++i)
        ret.fHistSparse->SetBinContent(i, rnd.Poisson(ret.fHistSparse->GetBinContent(i)));
    }

    // Drop old errors, which are based on the MC statistics, and create new
    // ones that are based on the prediction for the data
    if(ret.fHistD){
      ret.fHistD->Sumw2(false);
      ret.fHistD->Sumw2();
    }
    if(ret.fHistF){
      ret.fHistF->Sumw2(false);
      ret.fHistF->Sumw2();
    }

    return ret;
  }

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

inline

Definition at line 297 of file Spectrum.h.

References fLabels.

Referenced by ana::SpectrumComponents::NDimensions(), ana::CutOptimization::OptimizedSigmaOverSigma(), ana::Multiverse::ToHist(), and ana::CutOptimization::ToHist().

{return fLabels.size();}

Spectrum ana::Spectrum::operator*

(

const Ratio &

rhs

)

const

Definition at line 1016 of file Spectrum.cxx.

References runNovaSAM::ret.

Referenced by OverrideLivetime().

  {
    Spectrum ret = *this;
    ret *= rhs;
    return ret;
  }

Spectrum & ana::Spectrum::operator*=

(

const Ratio &

rhs

)

Definition at line 1008 of file Spectrum.cxx.

References ana::Ratio::fHist, fHistD, fHistF, and make_syst_table_plots::h.

Referenced by OverrideLivetime().

  {
    TH1* h = fHistD ? (TH1*)fHistD : (TH1*)fHistF;
    h->Multiply(rhs.fHist);
    return *this;
  }

Spectrum ana::Spectrum::operator+

(

const Spectrum &

rhs

)

const

Definition at line 986 of file Spectrum.cxx.

References runNovaSAM::ret.

Referenced by OverrideLivetime().

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

Spectrum & ana::Spectrum::operator+=

(

const Spectrum &

rhs

)

Definition at line 980 of file Spectrum.cxx.

References PlusEqualsHelper().

Referenced by OverrideLivetime().

  {
    return PlusEqualsHelper(rhs, +1);
  }

Spectrum ana::Spectrum::operator-

(

const Spectrum &

rhs

)

const

Definition at line 1000 of file Spectrum.cxx.

References runNovaSAM::ret.

Referenced by OverrideLivetime().

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

Spectrum & ana::Spectrum::operator-=

(

const Spectrum &

rhs

)

Definition at line 994 of file Spectrum.cxx.

References PlusEqualsHelper().

Referenced by OverrideLivetime().

  {
    return PlusEqualsHelper(rhs, -1);
  }

Spectrum ana::Spectrum::operator/

(

const Ratio &

rhs

)

const

Definition at line 1032 of file Spectrum.cxx.

References runNovaSAM::ret.

Referenced by OverrideLivetime().

  {
    Spectrum ret = *this;
    ret /= rhs;
    return ret;
  }

Spectrum & ana::Spectrum::operator/=

(

const Ratio &

rhs

)

Definition at line 1024 of file Spectrum.cxx.

References ana::Ratio::fHist, fHistD, fHistF, and make_syst_table_plots::h.

Referenced by OverrideLivetime().

  {
    TH1* h = fHistD ? (TH1*)fHistD : (TH1*)fHistF;
    h->Divide(rhs.fHist);
    return *this;
  }

Spectrum & ana::Spectrum::operator=

(

const Spectrum &

rhs

)

Definition at line 426 of file Spectrum.cxx.

References ana::assert(), ana::HistCache::Copy(), ana::HistCache::Delete(), fBins, fHistD, fHistF, fHistSparse, fLabels, fLivetime, fLoaderCount, fPOT, and GetMultiDBinning().

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

    DontAddDirectory guard;

    if(fHistD) HistCache::Delete(fHistD, GetMultiDBinning());
    if(fHistF) HistCache::Delete(fHistF, GetMultiDBinning());
    delete fHistSparse;

    assert(rhs.fHistD || rhs.fHistF || rhs.fHistSparse);

    if(rhs.fHistD){
      fHistD = HistCache::Copy(rhs.fHistD, rhs.GetMultiDBinning());
      fHistF = 0;
      fHistSparse = 0;
    }

    if(rhs.fHistF){
      fHistF = HistCache::Copy(rhs.fHistF, rhs.GetMultiDBinning());
      fHistD = 0;
      fHistSparse = 0;
    }

    if(rhs.fHistSparse){
      fHistSparse = (THnSparseD*)rhs.fHistSparse->Clone();
      fHistD = 0;
      fHistF = 0;
    }

    fPOT = rhs.fPOT;
    fLivetime = rhs.fLivetime;
    fLabels = rhs.fLabels;
    fBins = rhs.fBins;

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

    return *this;
  }

Spectrum & ana::Spectrum::operator=

(

Spectrum &&

rhs

)

Definition at line 467 of file Spectrum.cxx.

References ana::assert(), ana::HistCache::Delete(), fBins, fHistD, fHistF, fHistSparse, fLabels, fLivetime, fLoaderCount, fPOT, GetMultiDBinning(), and demo0::loader.

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

    if(fHistD) HistCache::Delete(fHistD, GetMultiDBinning());
    if(fHistF) HistCache::Delete(fHistF, GetMultiDBinning());
    delete fHistSparse;

    assert(rhs.fHistD || rhs.fHistF || rhs.fHistSparse);

    fHistD = rhs.fHistD;
    fHistF = rhs.fHistF;
    fHistSparse = rhs.fHistSparse;

    fPOT = rhs.fPOT;
    fLivetime = rhs.fLivetime;
    fLabels = rhs.fLabels;
    fBins = rhs.fBins;

    rhs.fHistD = 0;
    rhs.fHistF = 0;
    rhs.fHistSparse = 0;

    for (auto & loader : rhs.fLoaderCount)
    {
      loader->ReplaceSpectrum(&rhs, this);
      fLoaderCount.insert(loader);
    }
    rhs.fLoaderCount.clear();

    return *this;
  }

Spectrum & ana::Spectrum::operator=

(

const SpectrumStan &

rhs

)

Definition at line 501 of file Spectrum.cxx.

References Spectrum(), and ~Spectrum().

  {
    // use placement new to just overwrite this object using the relevant constructor
    this->~Spectrum();
    new (this) Spectrum(rhs);
    return *this;
  }

void ana::Spectrum::OverrideLivetime ( double  newlive )

inline

DO NOT USE UNLESS YOU ARE 110% CERTAIN THERE ISN'T A BETTER WAY!

Definition at line 272 of file Spectrum.h.

References make_syst_table_plots::c, Clear(), dir, fLivetime, LoadFrom(), operator*(), operator*=(), operator+(), operator+=(), operator-(), operator-=(), operator/(), operator/=(), SaveTo(), and Scale().

Referenced by getBNBPlots(), make_fc_nus_surfs_nersc_2019(), make_fc_points(), sensitivity(), sensitivity_slices(), and timingPeak().

{fLivetime = newlive;}

void ana::Spectrum::OverridePOT ( double  newpot )

inline

DO NOT USE UNLESS YOU ARE 110% CERTAIN THERE ISN'T A BETTER WAY!

Definition at line 269 of file Spectrum.h.

References fPOT.

Referenced by BlessedPlotsLoad(), ana::DataMCAreaNormalizedRatio(), ana::DataMCPair::DrawMCNormSyst(), DataMCPair::DrawMCNormSyst(), ana::FluxMultiverseSyst::FindSigmaBoundaries(), ana2019::fakedata::Get2019Prediction(), ana::nueccinc::NueCCIncCrossSectionAnalysis::getTemplateShapeOnly1D(), ana::nueccinc::NueCCIncCrossSectionAnalysis::getTemplateShapeOnly3D(), make_fc_nus_surfs_nersc_2018(), make_fc_nus_surfs_nersc_2019(), make_fc_points(), MakeExtrapSurface(), MakeHBar(), MakeNusPrediction(), PlotNus17Prediction(), ana::PlotSpectra(), PlotStack(), PlotSysts(), and SpectrumParamEffectsAna().

{fPOT = newpot;}

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

protected

Helper for operator+= and operator-=.

Definition at line 905 of file Spectrum.cxx.

References ana::AlmostEqual(), om::cout, allTimeWatchdog::endl, fHistD, fHistF, fHistSparse, fLivetime, fPOT, make_syst_table_plots::h, and rh.

Referenced by GetBinnings(), operator+=(), and operator-=().

  {
    // In this case it would be OK to have no POT/livetime
    if(rhs.fHistD && rhs.fHistD->Integral(0, -1) == 0) return *this;
    if(rhs.fHistF && rhs.fHistF->Integral(0, -1) == 0) return *this;

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

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

    if(!fPOT && !rhs.fLivetime){
      std::cout << "Error: can't sum Spectrum with livetime ("
                << fLivetime << " sec) but no POT and Spectrum 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
      TH1* h = fHistD ? (TH1*)fHistD : (TH1*)fHistF;
      TH1* rh = rhs.fHistD ? (TH1*)rhs.fHistD : (TH1*)rhs.fHistF;
      if(rh) h->Add(rh, sign*fPOT/rhs.fPOT);
      if(rhs.fHistSparse) fHistSparse->Add(rhs.fHistSparse, 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
      TH1* h = fHistD ? (TH1*)fHistD : (TH1*)fHistF;
      TH1* rh = rhs.fHistD ? (TH1*)rhs.fHistD : (TH1*)rhs.fHistF;
      if(rh) h->Add(rh, sign*fLivetime/rhs.fLivetime);
      if(rhs.fHistSparse) fHistSparse->Add(rhs.fHistSparse, 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 << "Spectrum::operator+=(). How did we get here? "
              << fPOT << " " << fLivetime << " "
              << rhs.fPOT << " " << rhs.fLivetime << std::endl;
    abort();
  }

double ana::Spectrum::POT ( ) const

inline

Definition at line 263 of file Spectrum.h.

References fPOT.

Referenced by ana::nueccinc::NueCCIncCrossSectionAnalysis::AddEnhancedSample(), ana::DataMCPair::AddExposure(), ana::SpectrumLoaderMockData::AddSpectrum(), alternate_cosrej_table(), caf_numu_fd_validation_data(), CalcChi2(), CalcR(), CalcRWithSystsNus17(), calculate_avg_nslcs(), calculateComponentsNumu(), calculateWrongSignNue(), calculateWrongSignNumuQ1(), calculateWrongSignNumuQ2(), calculateWrongSignNumuQ3(), calculateWrongSignNumuQ4(), ana::CountingExperiment::ChiSq(), ana::SingleSampleExperiment::ChiSq(), ana::CovMxExperiment::ChiSq(), compare_pid_cut_sens(), CompareDecompDataMC(), ComputeEfficiency(), Cumulative(), ana::DataMCAreaNormalizedRatio(), ana::DataMCComparison(), Plotter::DataMCComparison(), ana::DataMCComparisonAreaNormalized(), ana::DataMCComparisonComponents(), ana::TwoSampleDecomp::Decomp(), ana::FakeDecomp::Decompose(), demo6(), ana::SingleSampleExperiment::Derivative(), draw_decomp_plots(), ana::DataMCPair::DrawData(), ana::SpectrumComponents::DrawLegend(), ana::DataMCPair::DrawMCComponents(), ana::DataMCPair::DrawMCNormSyst(), DataMCPair::DrawMCNormSyst(), ana::DataMCPair::DrawMCSyst(), drawPlot(), ana::TwoSampleDecomp::DrawSigBkgOverlay(), ana::TwoSampleDecomp::DrawTwoSamplesWithRatios(), EHadVisMECpairs(), fa_nd_mc_prediction_full_syst(), FCContour(), fd_plot(), fd_plot_bless(), fill_col(), ana::FluxMultiverseSyst::FindSigmaBoundaries(), FitParamEffectsAna(), ana::PredictionInterp::FitSumComponent(), FitSystEffectsAna(), GetExtrap(), getHists_FNEX(), ana::Multiverse::GetNSigmaShift(), GetNueData2017(), GetNueData2018(), GetNueData2019(), ana::GetNueData2020(), ana::GetNumuData2020(), GetSpectra(), ana::MichelDecomp::GetSum(), ana::MichelDecomp::GetTemplateContent(), ana::nueccinc::NueCCIncCrossSectionAnalysis::getTemplateShapeOnly1D(), ana::nueccinc::NueCCIncCrossSectionAnalysis::getTemplateShapeOnly3D(), ana::PredictionNueRebinSA::Inflate(), ana::PredictionInterp::InitFits(), joint_fit_make_fc_slice(), joint_fit_make_fc_surf(), ana::SingleSampleExperiment::LogLikelihood(), Make2DPlot(), make_fc_points(), make_fc_slice(), make_fc_surf(), make_muonid_opt(), make_nus17_fc_surfs(), make_nus_fc_surfs(), make_plots(), make_vertex_optimiz(), make_xs(), make_xs_1D(), makeEnergyEstimator(), makeFlatWeight(), MakeHBar(), makeMatrixElementSurface(), MakePlot(), MakePlots(), ana::MakeQuantileHistogram(), ana::FluxDecomp::MakeWeightsNumuFromKaon(), ana::BENDecomp::MakeWeightsNumuFromKaon(), ana::FluxDecomp::MakeWeightsNumuFromPion(), ana::BENDecomp::MakeWeightsNumuFromPion(), ana::MichelDecomp::MCToDCMPComp(), mec_tuning_fitter_2020(), ana::PredictionExtendToPeripheral::MergePeripheral(), MichelDecompTest(), mrbrem_get_reweighted_spectra(), mrbrem_plots(), mre_blessed(), mre_comp_split(), mre_example(), nd_plot_bless(), ana::GenieMultiverseNormalizedSpectra::NormalizeSpectra(), nue_bg_stat_err_calc(), nue_cosrej_study(), nue_decomp_scales(), nue_decomp_scales_forMakeMichel(), ana::FluxDecomp::NueEstimate(), ana::BENDecomp::NueEstimate(), ana::FluxDecomp::NueEstimateFromKa(), ana::BENDecomp::NueEstimateFromKa(), ana::FluxDecomp::NueEstimateFromPi(), ana::BENDecomp::NueEstimateFromPi(), ana::CutOptimization::OptimizedSigmaOverSigma(), plot(), plot_2d(), plot_2d_vars(), plot_3NDvsFD(), plot_3NDvsFD_FHC(), plot_3NDvsFD_RHC(), plot_datamc_ND_numu(), plot_datamc_ND_numu_REW(), plot_datamcpred(), plot_diff(), plot_fa_results(), plot_nd_data_mc(), plot_ND_numu_NOMvsREW(), plot_NDvsFD_REW(), plot_NDvsFD_weights(), plot_NDvsFD_weights_FHC(), plot_NDvsFD_weights_RHC(), plot_predictions(), plot_recoE_numu(), plot_time(), ana::PlotAllSelectionDecomposition(), ana::PlotAllSignalEstimates(), ana::CrossSectionAnalysis::PlotBackgroundEstimate(), ana::NumuCCIncAnalysis::PlotBackgroundEstimate(), ana::SingleNucAnalysis::PlotBackgroundEstimate2D(), ana::NumuCC2p2hAnalysis::PlotBackgroundEstimate2D(), ana::CrossSectionAnalysis::PlotData(), ana::NumuCCIncAnalysis::PlotData(), ana::SingleNucAnalysis::PlotData2D(), ana::NumuCC2p2hAnalysis::PlotData2D(), ana::CutOptimization::PlotDebug(), PlotNueDecompFourBins(), PlotNus17PredSystsData(), PlotNus18Sideband(), PlotNusSensAna01(), ana::CrossSectionAnalysis::PlotRecoToTrueMatrix(), ana::SingleNucAnalysis::PlotRecoToTrueMatrix2D(), ana::NumuCC2p2hAnalysis::PlotRecoToTrueMatrix2D(), ana::CrossSectionAnalysis::PlotSignalEstimate(), ana::SingleNucAnalysis::PlotSignalEstimate2D(), ana::NumuCC2p2hAnalysis::PlotSignalEstimate2D(), ana::PlotSpectra(), PlotStack(), ana::PlotStack(), PlotSyst(), PlotSysts(), ana::CrossSectionAnalysis::PlotUnfoldedSignal(), ana::SingleNucAnalysis::PlotUnfoldedSignal2D(), ana::NumuCC2p2hAnalysis::PlotUnfoldedSignal2D(), PlotVertices(), ana::SingleSampleExperiment::PredHistIncCosmics(), ana::NDPredictionSterile::Predict(), ana::NDPredictionSterile::PredictComponent(), ana::PredictionNumuFAHadE::PredictComponentSyst(), PrintEventCounts(), ana::TwoSampleDecomp::RatioCalc(), reco_minus_true_panels(), ana::PredictionNueRebinSA::Reduce(), ana::PredictionExtendToPeripheral::ReduceHelperNC(), ana::CrossSectionAnalysis::Result(), ana::NumuCCIncAnalysis::Result(), ana::NumuCC2p2hAnalysis::Result1DFluxInt(), ana::SingleNucAnalysis::Result2D(), ana::NumuCC2p2hAnalysis::Result2D(), ana::NumuCCIncAnalysis::Result2D(), ana::NumuCCIncAnalysis::ReturnHists(), ana::MichelDecomp::SaveCompPlots(), ana::FluxDecomp::SavePlots(), ana::BENDecomp::SavePlots(), ana::TruthReweight::SavePlots(), ana::RecoReweight::SavePlots(), ana::FluxDecomp::SavePlotsKa(), ana::BENDecomp::SavePlotsKa(), ana::FluxDecomp::SavePlotsPi(), ana::BENDecomp::SavePlotsPi(), saveS(), saveS1(), ana::MichelDecomp::SaveTempPlots(), selection_story_plots(), sensitivity(), ana::PredictionInterp::ShiftSpectrum(), ana::SimpleFOM(), SpectrumParamEffectsAna(), ana::SpectrumStan::SpectrumStan(), TableNueNDComponents(), TableNueNDDataMC(), test_beam_errorband(), test_micheldecomp(), test_sam(), test_sam_project(), test_stanfit_systpulls(), ana::Multiverse::ToHist(), ana::CutOptimization::ToHist(), ana::ICrossSectionAnalysis::ToHist(), ana::UnfoldIterative::Truth(), ana::UnfoldSVD::Truth(), ana::UnfoldTikhonov::Truth(), and ana::PredictionTemplateShift::Weight().

{return fPOT;}

void ana::Spectrum::RemoveLoader ( SpectrumLoaderBase *  p )

protected

Definition at line 897 of file Spectrum.cxx.

References fLoaderCount.

Referenced by GetBinnings().

  { fLoaderCount.erase(p); }

void ana::Spectrum::SaveTo

(

TDirectory *

dir

)

const

Definition at line 1040 of file Spectrum.cxx.

References fBins, fHistD, fHistF, fHistSparse, fLabels, fLivetime, genie::utils::style::Format(), fPOT, MECModelEnuComparisons::i, and tmp.

Referenced by Ana2017_box_opening_macro(), Ana2017_sb_opening_macro(), Ana2018_box_opening_macro(), angle(), bdtstudyspectrums(), BlessedPlotsLoad(), caf_nue_data_mc(), ccpiinc_mc_studies(), containmentstudy(), CosmicPred(), ana::CovMxManager::CovMxManager(), CutTableLoad(), CVNphoton(), CVNphotonSplit(), datamc_ND_numu_kinematics(), datamc_ND_numu_kinematics_FHC(), datamc_ND_numu_kinematics_FHC_REW(), datamc_ND_numu_kinematics_REW(), datamc_ND_numu_kinematics_RHC(), datamc_ND_numu_kinematics_RHC_REW(), dataprocess_numuccinc(), estimate_energy(), fd_data_mc(), FDDataMCSystBandLoad(), FidOpt(), FidWShwCuts(), FillSpectra(), FitSystEffectsLoad(), get_cosmic_spectra(), get_data_and_cosmic(), get_data_histogram(), get_numi_data_histogram(), getData(), getFlatWeightSpectra(), GetNDDataMC(), getSpectra_ForFitting(), make_fd_data_spec(), make_mockdata_syst_contours(), make_nue_fd_datamc(), make_nue_firstana_exposure(), make_pi0_xcheck(), make_pid(), MakeCosmics(), MakeFakeData(), makeFakeDataFluxes(), MakeNus17Prediction(), MakeNus18SidebandPred(), MakeNusPrediction(), makeRealDataFluxes(), makeXSecPlots_TemplateFit(), mrbrem_get_initial_spectra(), mrbrem_get_reweighted_spectra(), MRDiFStudy_FHC_Step1(), MRDiFStudy_RHC_Step1(), mre_blessed(), mre_comp_split(), mre_sel_efficiency(), ncpi0HistoGrid2(), nd_data_mc(), ND_DataMC(), NDDataMCSystBandLoad(), nue_bg_stat_errs_load(), nue_decomp_scales(), nus17_box_opening(), nus17_fd_cut_tables(), nus17_fd_cut_tables2D(), nus18_box_opening(), nus_ana01_box_opening(), nus_ana01_sideband_box_opening(), OverrideLivetime(), pion_multiverse(), plot_fa_results(), preselection_cutflow(), ReMId(), resolutionscript(), resolutionstudy(), sa_box_opening_macro(), sa_fd_cut_tables(), sa_hesideband_box_opening_macro(), saveFDMCHists(), saveSpectraForUnf(), ana::xsec::UnfoldingVariable::SaveSpectrums(), ana::PredictionXSecTuning::SaveTo(), ana::CountingExperiment::SaveTo(), ana::PredictionNoOsc::SaveTo(), ana::PredictionNueRebinSA::SaveTo(), jw::TrivialPrediction::SaveTo(), ana::CheatDecomp::SaveTo(), ana::FitInAnaBinsBkgdEstimator::SaveTo(), ana::NumuCC2p2hBkgdEstimator::SaveTo(), ana::nueccinc::NueCCIncMRECorrection::SaveTo(), ana::NCDecomp::SaveTo(), ana::NueDecomp::SaveTo(), ana::TrivialPrediction::SaveTo(), ana::SingleSampleExperiment::SaveTo(), ana::CrossSectionSpectra::SaveTo(), ana::NumuDecomp::SaveTo(), ana::PredictionScaleComp::SaveTo(), ana::ProportionalDecomp::SaveTo(), ana::CutOptimization::SaveTo(), ana::TrivialCrossSectionAnalysis::SaveTo(), ana::FakeDecomp::SaveTo(), ana::TrivialExtrap::SaveTo(), ana::CrossSectionAnalysis::SaveTo(), ana::TwoSampleDecomp::SaveTo(), ana::SpectrumComponents::SaveTo(), ana::NumuCCIncAnalysis::SaveTo(), ana::FluxDecomp::SaveTo(), ana::nueccinc::NueCCIncCrossSectionAnalysis::SaveTo(), ana::MichelDecomp::SaveTo(), ana::BENDecomp::SaveTo(), ana::DataMCPair::SaveTo(), ana::RecoReweight::SaveTo(), ana::nueccinc::NueCCIncEnhancedSamples::SaveTo(), ana::FluxReweight::SaveTo(), ShwZOpt(), sidebandfittest(), test_beam_errorband(), test_fluxhadr_prod_weights_Flux(), test_genieweights(), test_predictionscalecomp(), twodstudyvtxcont(), validationscript(), and vertexstudy().

  {
    TDirectory* tmp = gDirectory;
    dir->cd();

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

    if(fHistD) fHistD->Write("hist");
    if(fHistF) fHistF->Write("hist");
    if(fHistSparse) fHistSparse->Write("hist_sparse");
    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 < fBins.size(); ++i){
      TObjString(fLabels[i].c_str()).Write(TString::Format("label%d", i).Data());
      fBins[i].SaveTo(dir->mkdir(TString::Format("bins%d", i)));
    }

    tmp->cd();
  }

void ana::Spectrum::Scale

(

double

c

)

Multiply this spectrum by a constant c.

Definition at line 745 of file Spectrum.cxx.

References fHistD, fHistF, and fHistSparse.

Referenced by ComputeEfficiency(), make_fc_nus_surfs_nersc_2018(), make_fc_nus_surfs_nersc_2019(), make_xs(), make_xs_1D(), MakeExtrapSurface(), ana::FluxDecomp::MakeWeightsNumuFromKaon(), ana::BENDecomp::MakeWeightsNumuFromKaon(), ana::GenieMultiverseNormalizedSpectra::NormalizeSpectra(), OverrideLivetime(), PlotNus17Prediction(), ana::PredictionAddRock::PredictComponentSyst(), ana::NumuCCIncAnalysis::Result(), ana::NumuCCIncAnalysis::Result2D(), and ana::PredictionExtendToPeripheral::Scale().

  {
    if(fHistD) fHistD->Scale(c);
    if(fHistF) fHistF->Scale(c);
    if(fHistSparse) fHistSparse->Scale(c);
  }

TH1D * ana::Spectrum::ToTH1	(	double	exposure,
		Color_t	col = kBlack,
		Style_t	style = kSolid,
		EExposureType	expotype = kPOT,
		EBinType	bintype = kBinContent
	)		const

Histogram made from this Spectrum, scaled to some exposure.

Parameters

exposure	POT or livetime (seconds)
col	Histogram color (default black)
style	Histogram line style (default solid)
expotype	How to interpret exposure (kPOT (default) or kLivetime)

Definition at line 564 of file Spectrum.cxx.

References ana::HistCache::Copy(), om::cout, allTimeWatchdog::endl, fBins, fHistD, fHistF, fHistSparse, fLabels, fLivetime, fPOT, GetMultiDBinning(), MECModelEnuComparisons::i, ana::kBinDensity, ana::kLivetime, ana::kPOT, submit_hadd::l, label, livetime, ana::HistCache::New(), pot, and runNovaSAM::ret.

Referenced by ana::nueccinc::NueCCIncCrossSectionAnalysis::AddEnhancedSample(), ana::AddErrorInQuadrature(), ana::PredictionSystJoint2018::AddNormSyst(), ana::PredictionSystJointDemo::AddNormSyst(), ana::PredictionSyst3Flavor2020::AddNormSyst(), ana::PredictionSystNueSecondAna::AddNormSysts(), ana::PredictionSystNue2017::AddNormSysts(), ana::PredictionSystNumu2017::AddNormSysts(), alternate_cosrej_table(), AnalyzeNus18Pred(), ana::FitInAnaBinsBkgdEstimator::Background(), ana::NumuCC2p2hBkgdEstimator::Background(), bin_composition_pie_chart(), caf_numu_fd_validation_data(), caf_numu_fd_validation_MC(), caf_numu_fd_validation_MC_no_tau(), calcAlphaBetaEachBin(), CalcChi2(), ana::nueccinc::NueCCIncMRECorrection::CalcEfficiencyCorrection(), CalcR(), CalcRWithSystsNus17(), calculate_avg_nslcs(), calculateComponentsNumu(), calculateWrongSignNue(), calculateWrongSignNumuQ1(), calculateWrongSignNumuQ2(), calculateWrongSignNumuQ3(), calculateWrongSignNumuQ4(), cc(), check_predinterp(), check_predinterp_numu(), ana::OscCovMxExperiment::ChiSq(), ana::SingleSampleExperiment::ChiSq(), ana::CovMxExperiment::ChiSq(), Compare(), compare_pid_cut_sens(), CompareDecompDataMC(), CompareMCCompPrediction(), ComparePredictions(), ana::ComparePredictions(), ana::ModularExtrapComponent::ComparisonPlot(), ana::ComparisonTable(), ana::PredictionInterp::ComponentDerivative(), Cumulative(), CVNCuts(), ana::DataMCComparison(), Plotter::DataMCComparison(), ana::DataMCComparisonAreaNormalized(), ana::DataMCComparisonComponents(), ana::TwoSampleDecomp::Decomp(), ana::FakeDecomp::Decompose(), ana::MichelDecomp::Decompose(), demo1(), demo3(), demo5(), demo6(), demo_flat(), demo_nueNumuSysts(), ana::PredictionAddRock::Derivative(), ana::SingleSampleExperiment::Derivative(), ana::DataMCPair::DrawData(), drawPlot(), ana::TwoSampleDecomp::DrawSigBkgOverlay(), DrawSurface(), drawSystsShiftingNDdata(), drawSystsShiftingNDdata_updatedAna(), ana::TwoSampleDecomp::DrawTwoSamplesWithRatios(), efficiency(), efficiencySA(), EHadVisMECpairs(), fa_fd_data_mc_systs(), fa_nd_mc_prediction_full_syst(), FCContour(), fd_plot(), fd_plot_bless(), fill_col(), ana::FluxMultiverseSyst::FindSigmaBoundaries(), first_ana_proj(), FitParamEffectsAna(), ana::PredictionInterp::FitSumComponent(), FitSystEffectsAna(), ana::ModularExtrapComponent::FormSmartRatio(), GenerateFutureData(), GeniePredictionToRoot(), GetBackgroundStatisticalUncertainty(), GetBackgroundSystematicUncertainty(), ana::GetBFSystBands(), GetBG(), ana::GetBG(), getBNBPlots(), getContProf(), getContProf_Sensitivity(), GetDenominator(), ana::nueccinc::NueCCIncMRECorrection::getEfficiency1D(), GetEfficiencyDenominator(), GetEfficiencySystematicUncertainty(), GetExtrap(), GetFakeData(), ana::nueccinc::NueCCIncCrossSectionAnalysis::getFlux1D(), GetFluxError(), ana::nueccinc::NueCCIncCrossSectionAnalysis::getFullTemplate(), getHists_FNEX(), GetNC(), ana::GetNC(), ana::GetNDComponents(), ana::GetNDCompsFromDecomp(), ana::GetNDDecompsFromDecomp(), ana::Multiverse::GetNSigmaShift(), GetNueCosmics2017(), GetNueCosmics2018(), GetNueCosmics2019(), GetNueCosmicsFuture(), GetNumuCosmicsFuture(), ana::nueccinc::NueCCIncMRECorrection::getPreselected1D(), GetQuantilePredictionHist(), ana::nueccinc::NueCCIncMRECorrection::getSelected1D(), GetSelectedStatisticalUncertainty(), GetSig(), GetSpectra(), ana::GetSpectrum(), ana::GetSystBands(), ana::nueccinc::NueCCIncCrossSectionAnalysis::getTemplateShapeOnly1D(), ana::nueccinc::NueCCIncCrossSectionAnalysis::getTemplateShapeOnly3D(), getTimePeakPlots(), ana::PredictionNueRebinSA::Inflate(), Integral(), joint_fit_make_fc_slice(), joint_fit_make_fc_surf(), ana::SingleSampleExperiment::LogLikelihood(), make_extrap_figure(), make_fc_nus_surfs_nersc_2018(), make_fc_points(), make_fc_slice(), make_fc_surf(), make_muonid_opt(), make_nus17_fc_surfs(), make_nus_fc_surfs(), make_plots(), make_vertex_optimiz(), make_xs_1D(), MakeCovMx(), makeEnergyEstimator(), makeFlatWeight(), makeMatrixElementSurface(), MakePeriCutPlot(), MakePlot(), MakePlots(), MakeResultPlot(), MakeSurface(), MakeTable(), ana::FluxDecomp::MakeWeightsNumuFromKaon(), ana::BENDecomp::MakeWeightsNumuFromKaon(), makeXSecPlots1D(), makeXSecPlots2D(), ana::MichelDecomp::MCToDCMPComp(), ana::PredictionExtendToPeripheral::MergePeripheral(), MichelDecompTest(), mrbrem_get_reweighted_spectra(), mrbrem_plots(), mre_blessed(), mre_comp_split(), myPlotWithSystErrorBand(), nd_plot_bless(), NDHists(), nhits(), nue_ana_basic(), nue_bg_stat_err_calc(), nue_cosrej_study(), nue_decomp_scales(), nue_decomp_scales_forMakeMichel(), nue_pid_effs(), nue_pid_effs_miniprod(), nue_pid_effs_paper_numu_neweff(), ana::FluxDecomp::NueEstimate(), ana::BENDecomp::NueEstimate(), ana::BENDecomp::NueEstimateFromKa(), numu_cut_flow(), numu_demo_split_plot(), PeripheralCuts(), pi0_xcheck(), plot(), plot_3NDvsFD(), plot_3NDvsFD_FHC(), plot_3NDvsFD_RHC(), plot_datamc_ND_numu(), plot_datamc_ND_numu_REW(), plot_datamcpred(), plot_diff(), plot_fa_results(), plot_fd_datamc(), plot_kinematics_cafana(), plot_nd_data_mc(), plot_ND_numu_NOMvsREW(), plot_NDvsFD_REW(), plot_NDvsFD_weights(), plot_NDvsFD_weights_FHC(), plot_NDvsFD_weights_RHC(), plot_nue_xsec_pred(), plot_nueFD_Signal_REWvsNOM_FHC(), plot_nueFD_Signal_REWvsNOM_pTExtrap_FHC(), plot_nueFD_Signal_REWvsNOM_pTExtrap_RHC(), plot_nueFD_Signal_REWvsNOM_RHC(), plot_pi0_xcheck(), plot_pid(), plot_predictions(), plot_recoE_numu(), plot_sa_results(), plot_time(), ana::PlotAllSelectionDecomposition(), ana::PlotAllSignalEstimates(), ana::CrossSectionAnalysis::PlotBackgroundEstimate(), ana::NumuCCIncAnalysis::PlotBackgroundEstimate(), PlotComp(), ana::CrossSectionAnalysis::PlotData(), ana::NumuCCIncAnalysis::PlotData(), ana::CutOptimization::PlotDebug(), ana::CrossSectionAnalysis::PlotFluxEstimate(), ana::SingleNucAnalysis::PlotFluxEstimate2D(), ana::NumuCC2p2hAnalysis::PlotFluxEstimate2D(), PlotNuePredictionFourBins(), PlotNumuPredData(), PlotNus17PredSystsData(), PlotNus18Sideband(), PlotNusSensAna01(), PlotPionPlots(), PlotPurEff(), ana::PlotPurEff(), plots(), ana::PlotSpectra(), PlotStack(), ana::PlotStack(), PlotSyst(), PlotSysts(), Plotting_DataAndPrediction(), ana::CrossSectionAnalysis::PlotUnfoldedSignal(), ana::PlotWithSystErrorBand(), ana::PlotWithSystErrorBand_Quant(), ana::PlotWithSystErrorBandTwoPreds(), ana::SingleSampleExperiment::PredHistIncCosmics(), ana::PredictionCombinePeriods::PredictComponentSyst(), ana::PredictionNumuFAHadE::PredictComponentSyst(), PredictCosmic(), preselection_cutflow(), printline(), ana::SpectrumComponents::Purity(), ana::Ratio::Ratio(), ratio(), ana::TwoSampleDecomp::RatioCalc(), reco_minus_true_panels(), ana::PredictionNueRebinSA::Reduce(), ana::PredictionExtendToPeripheral::ReduceHelperNC(), resolution2018(), ana::CrossSectionAnalysis::Result(), ana::NumuCCIncAnalysis::Result(), ana::NumuCC2p2hAnalysis::Result1DFluxInt(), ana::NumuCC2p2hAnalysis::Result2D(), ana::SingleNucAnalysis::Result2D(), ana::NumuCCIncAnalysis::Result2D(), ana::NumuCCIncAnalysis::ReturnHists(), rhc_sensitivity(), ana::MichelDecomp::SaveCompPlots(), SaveDCMPPlots(), SaveOscPrediction(), ana::FluxDecomp::SavePlots(), ana::BENDecomp::SavePlots(), ana::TruthReweight::SavePlots(), ana::RecoReweight::SavePlots(), ana::FluxReweight::SavePlots(), ana::FluxDecomp::SavePlotsKa(), ana::BENDecomp::SavePlotsKa(), ana::BENDecomp::SavePlotsPi(), saveS(), selection_story_plots(), ana::PredictionInterp::ShiftSpectrum(), ana::SimpleFOM(), ana::SingleSampleExperiment::SingleSampleExperiment(), SpectrumParamEffectsAna(), ana::SpectrumStan::SpectrumStan(), spects(), sterile_demo(), syst_plot_test(), TableNuePredictionFourBins(), template_basic(), test_ana(), test_beam_errorband(), test_fluxhadr_prod_weights_Flux(), test_genie_systs(), test_genieweights(), test_micheldecomp(), test_nue2017Prediction(), test_nue2018_fitter(), test_nueextrapsyst(), test_nuwro(), test_saloaders(), test_sam(), test_sam_project(), test_stanfit_systpulls(), TestPred(), timingPeak(), ana::FluxMultiverseSyst::ToAreaNormalizedTH1(), ana::GenieMultiverseSpectra::ToAreaNormalizedTH1(), ana::Multiverse::ToHist(), ana::CutOptimization::ToHist(), ana::ICrossSectionAnalysis::ToHist(), ana::FluxMultiverseSyst::ToTH1(), ana::GenieMultiverseSpectra::ToTH1(), ToTH1(), ana::ToTH2(), ana::ToTH3(), Toy_analyses(), ana::UnfoldSVD::Truth(), ana::UnfoldTikhonov::Truth(), Unfold1D(), Unfold3D(), UnfoldInOut(), and ana::PredictionTemplateShift::Weight().

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

    TH1D* ret = 0;
    if(fHistD){
      ret = HistCache::Copy(fHistD, GetMultiDBinning());
    }
    else if(fHistF){
      ret = HistCache::New(fHistF->GetTitle(), GetMultiDBinning());
      ret->Add(fHistF);
    }
    else{
      ret = fHistSparse->Projection(0);
    }

    if(expotype == kPOT){
      const double pot = exposure;
      if(fPOT){
        ret->Scale(pot/fPOT);
      }
      else{
        // Allow zero POT if there are also zero events
        if(ret->Integral() > 0){
          std::cout << "Error: Spectrum with " << ret->Integral()
                    << " entries has zero POT, no way to scale to "
                    << exposure << " POT.";
          if(fLivetime > 0){
            std::cout << " Spectrum has " << fLivetime << " seconds livetime. "
                      << "Did you mean to pass kLivetime to ToTH1()?";
          }
          std::cout << std::endl;
          abort();
        }
      }
    }
    if(expotype == kLivetime){
      const double livetime = exposure;
      if(fLivetime){
        ret->Scale(livetime/fLivetime);
      }
      else{
        // Allow zero exposure if there are also zero events
        if(ret->Integral() > 0){
          std::cout << "Error: Spectrum with " << ret->Integral()
                    << " entries has zero livetime, no way to scale to "
                    << livetime << " seconds.";
          if(fPOT > 0){
            std::cout << " Spectrum has " << fPOT << " POT. "
                      << "Did you mean to pass kPOT to ToTH1()?";
          }
          std::cout << std::endl;
          abort();
        }
      }
    }

    std::string label;
    for(const std::string& l: fLabels) label += l + " and ";
    label.resize(label.size()-5); // drop the last "and"
    ret->GetXaxis()->SetTitle(label.c_str());

    if(fBins.size() == 1){
      for(unsigned int i = 0; i < fBins[0].Labels().size(); ++i){
        ret->GetXaxis()->SetBinLabel(i+1, fBins[0].Labels()[i].c_str());
      }
    }

    ret->GetYaxis()->SetTitle("Events");


    ret->SetLineColor(col);
    ret->SetMarkerColor(col);
    ret->SetLineStyle(style);

    if(bintype == kBinDensity) ret->Scale(1, "width");

    // Allow GetMean() and friends to work even if this histogram never had any
    // explicit Fill() calls made.
    if(ret->GetEntries() == 0) ret->SetEntries(1);

    return ret;
  }

TH1D* ana::Spectrum::ToTH1 ( double  exposure, EExposureType  expotype, EBinType  bintype = kBinContent  ) const

inline

Histogram made from this Spectrum, scaled to some exposure.

Parameters

exposure	POT or livetime (seconds)
expotype	How to interpret exposure (kPOT (default) or kLivetime)

Definition at line 221 of file Spectrum.h.

References FakeData(), compare_h5_caf::idx, Integral(), ana::kBinContent, ana::kPOT, Mean(), MockData(), ToTH1(), ToTH2(), ToTH2NormX(), and ToTH3().

    {
      return ToTH1(exposure, kBlack, kSolid, expotype, bintype);
    }

TH2 * ana::Spectrum::ToTH2	(	double	exposure,
		EExposureType	expotype = kPOT,
		EBinType	bintype = kBinContent
	)		const

Spectrum must be 2D to obtain TH2.

Definition at line 652 of file Spectrum.cxx.

References visualisationForPaperMasterPlot::ax, om::cout, allTimeWatchdog::endl, fBins, fLabels, MECModelEnuComparisons::i, ana::kBinDensity, runNovaSAM::ret, and ana::ToTH2().

Referenced by ana::nueccinc::NueCCIncMRECorrection::CalcEfficiencyCorrection(), demo6(), ana::nueccinc::NueCCIncCrossSectionAnalysis::doUnfolding1D(), ana::nueccinc::NueCCIncCrossSectionAnalysis::doUnfolding2D(), drawLongTerm(), drawTimePlots(), drawVsPOT(), FD_Data_PosComp(), FD_plots(), GetHistsFD(), GetHistsND(), GetSpectra2D(), getTimePeakPlots(), ana::nueccinc::NueCCIncCrossSectionAnalysis::getUnfoldingMatrix(), HadEFit(), Make2DPlot(), make_quantiles_histogram_2020(), makeEnergyEstimator(), MakePlots(), ana::MakeQuantileHistogram(), mec_tuning(), mec_tuning_fitter_2020(), MichelDecompTest(), mre_example(), MuonFit(), plot_2d(), plot_2d_vars(), plot_pi0_xcheck(), ana::SingleNucAnalysis::PlotBackgroundEstimate2D(), ana::NumuCC2p2hAnalysis::PlotBackgroundEstimate2D(), ana::SingleNucAnalysis::PlotData2D(), ana::NumuCC2p2hAnalysis::PlotData2D(), ana::SingleNucAnalysis::PlotUnfoldedSignal2D(), ana::NumuCC2p2hAnalysis::PlotUnfoldedSignal2D(), PlotVertices(), PositionComparison(), print_tables(), PrintPlot(), saveS1(), ana::Multiverse::ToHist(), ana::CutOptimization::ToHist(), ana::ICrossSectionAnalysis::ToHist(), ToTH1(), ana::GenieMultiverseSpectra::ToTH2(), and ToTH2NormX().

  {
    if(fBins.size() != 2){
      std::cout << "Error: This Spectrum does not appear to be 2D." << std::endl;
      abort();
    }

    TH2* ret = ana::ToTH2(*this, exposure, expotype, fBins[0], fBins[1]);

    for(int iax = 0; iax < 2; ++iax){
      TAxis* ax = ret->GetXaxis();
      if(iax == 1) ax = ret->GetYaxis();

      ax->SetTitle(fLabels[iax].c_str());

      for(unsigned int i = 0; i < fBins[iax].Labels().size(); ++i){
        ax->SetBinLabel(i+1, fBins[iax].Labels()[i].c_str());
      }
    }

    if(bintype == kBinDensity) ret->Scale(1, "width");

    // Allow GetMean() and friends to work even if this histogram never had any
    // explicit Fill() calls made.
    if(ret->GetEntries() == 0) ret->SetEntries(1);

    return ret;
  }

TH2 * ana::Spectrum::ToTH2NormX	(	double	exposure,
		EExposureType	expotype = kPOT
	)		const

Spectrum must be 2D to obtain TH2. Normalized to X axis.

If no entries in the column, skip normalization

Definition at line 682 of file Spectrum.cxx.

References fBins, MECModelEnuComparisons::i, calib::j, nbinsx, norm, and ToTH2().

Referenced by ToTH1().

  {
    TH2* xyhist = ToTH2(exposure, expotype);
    if(!xyhist) return nullptr;

    const int nbinsx = fBins[0].NBins();
    const int nbinsy = fBins[1].NBins();

    // Normalize 2D histogram to X-axis spectrum
    for(int i=1; i<=nbinsx; ++i){
      double norm = 0.0;
      for(int j=1; j<=nbinsy; ++j){
        norm += xyhist->GetBinContent(i, j);
      }
      /// If no entries in the column, skip normalization
      if(norm < 0.0000001) continue;

      norm = 1.0 / norm;
      for(int j=1; j<=nbinsy; ++j){
        xyhist->SetBinContent(i,j, xyhist->GetBinContent(i, j) * norm);
      }
    }
    
    // Allow GetMean() and friends to work even if this histogram never had any
    // explicit Fill() calls made.
    if(xyhist->GetEntries() == 0) xyhist->SetEntries(1);

    return xyhist;
  }

TH3 * ana::Spectrum::ToTH3	(	double	exposure,
		EExposureType	expotype = kPOT,
		EBinType	bintype = kBinContent
	)		const

Spectrum must be 3D to obtain TH3.

Definition at line 713 of file Spectrum.cxx.

References visualisationForPaperMasterPlot::ax, om::cout, allTimeWatchdog::endl, fBins, fLabels, MECModelEnuComparisons::i, ana::kBinDensity, runNovaSAM::ret, and ana::ToTH3().

Referenced by ana::FluxDecomp::MakeWeightsNumuFromPion(), ana::BENDecomp::MakeWeightsNumuFromPion(), ana::FluxDecomp::NueEstimateFromPi(), ana::BENDecomp::NueEstimateFromPi(), ana::FluxDecomp::SavePlotsPi(), ana::BENDecomp::SavePlotsPi(), ana::Multiverse::ToHist(), ana::CutOptimization::ToHist(), ana::ICrossSectionAnalysis::ToHist(), and ToTH1().

  {
    if(fBins.size() != 3){
      std::cout << "Error: This Spectrum does not appear to be 3D." << std::endl;
      abort();
    }

    TH3* ret = ana::ToTH3(*this, exposure, expotype,
                          fBins[0], fBins[1], fBins[2]);

    for(int iax = 0; iax < 3; ++iax){
      TAxis* ax = ret->GetXaxis();
      if(iax == 1) ax = ret->GetYaxis();
      if(iax == 2) ax = ret->GetZaxis();

      ax->SetTitle(fLabels[iax].c_str());

      for(unsigned int i = 0; i < fBins[iax].Labels().size(); ++i){
        ax->SetBinLabel(i+1, fBins[iax].Labels()[i].c_str());
      }
    }

    if(bintype == kBinDensity) ret->Scale(1, "width");

    // Allow GetMean() and friends to work even if this histogram never had any
    // explicit Fill() calls made.
    if(ret->GetEntries() == 0) ret->SetEntries(1);

    return ret;
  }

Friends And Related Function Documentation

friend class NullLoader

friend

Definition at line 43 of file Spectrum.h.

friend class SpectrumLoader

friend

Definition at line 42 of file Spectrum.h.

friend class SpectrumLoaderBase

friend

Definition at line 41 of file Spectrum.h.

friend class SpectrumStan

friend

Definition at line 40 of file Spectrum.h.

Member Data Documentation

std::vector<Binning> ana::Spectrum::fBins

protected

Definition at line 326 of file Spectrum.h.

Referenced by GetBinnings(), GetMultiDBinning(), operator=(), SaveTo(), ToTH1(), ToTH2(), ToTH2NormX(), and ToTH3().

TH1D* ana::Spectrum::fHistD

protected

Definition at line 316 of file Spectrum.h.

Referenced by Clear(), ConstructHistogram(), FakeData(), Fill(), Integral(), Mean(), MockData(), operator*=(), operator/=(), operator=(), PlusEqualsHelper(), SaveTo(), Scale(), Spectrum(), ToTH1(), and ~Spectrum().

TH1F* ana::Spectrum::fHistF

protected

Definition at line 317 of file Spectrum.h.

Referenced by Clear(), ConstructHistogram(), FakeData(), Fill(), Integral(), Mean(), MockData(), operator*=(), operator/=(), operator=(), PlusEqualsHelper(), SaveTo(), Scale(), Spectrum(), ToTH1(), and ~Spectrum().

THnSparseD* ana::Spectrum::fHistSparse

protected

Definition at line 318 of file Spectrum.h.

Referenced by Clear(), ConstructHistogram(), FakeData(), Fill(), MockData(), operator=(), PlusEqualsHelper(), SaveTo(), Scale(), Spectrum(), ToTH1(), and ~Spectrum().

std::vector<std::string> ana::Spectrum::fLabels

protected

Definition at line 325 of file Spectrum.h.

Referenced by GetLabels(), NDimensions(), operator=(), SaveTo(), ToTH1(), ToTH2(), and ToTH3().

double ana::Spectrum::fLivetime

protected

Definition at line 320 of file Spectrum.h.

Referenced by FakeData(), Integral(), Livetime(), operator=(), OverrideLivetime(), PlusEqualsHelper(), SaveTo(), Spectrum(), ana::SpectrumLoader::StoreExposures(), and ToTH1().

std::set<SpectrumLoaderBase*> ana::Spectrum::fLoaderCount

protected

This count is maintained by SpectrumLoader, as a sanity check.

Definition at line 323 of file Spectrum.h.

Referenced by AddLoader(), operator=(), RemoveLoader(), Spectrum(), ana::SpectrumStan::SpectrumStan(), and ~Spectrum().

double ana::Spectrum::fPOT

protected

Definition at line 319 of file Spectrum.h.

Referenced by FakeData(), Integral(), operator=(), OverridePOT(), PlusEqualsHelper(), POT(), SaveTo(), Spectrum(), ana::SpectrumLoader::StoreExposures(), and ToTH1().

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-07-07/CAFAna/Core/Spectrum.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-07-07/CAFAna/Core/Spectrum.cxx