ana::EnsembleRatio Class Reference

#include <EnsembleRatio.h>

Public Member Functions
	EnsembleRatio (const EnsembleSpectrum &nom, const EnsembleSpectrum &denom)
Ratio	Nominal () const
unsigned int	NUniverses () const
Ratio	Universe (unsigned int i) const
TGraphAsymmErrors *	ErrorBand () const
EnsembleRatio &	operator*= (const EnsembleRatio &rhs)
EnsembleRatio	operator* (const EnsembleRatio &rhs) const
EnsembleRatio &	operator/= (const EnsembleRatio &rhs)
EnsembleRatio	operator/ (const EnsembleRatio &rhs) const

Protected Attributes
Ratio	fNom
std::vector< Ratio >	fUnivs

Detailed Description

Definition at line 13 of file EnsembleRatio.h.

Constructor & Destructor Documentation

ana::EnsembleRatio::EnsembleRatio	(	const EnsembleSpectrum &	nom,
		const EnsembleSpectrum &	denom
	)

Definition at line 11 of file EnsembleRatio.cxx.

    : fNom(num.Nominal(), denom.Nominal())
  {
    assert(num.NUniverses() == denom.NUniverses());
    fUnivs.reserve(num.NUniverses());
    for(unsigned int i = 0; i < num.NUniverses(); ++i)
      fUnivs.emplace_back(num.Universe(i), denom.Universe(i));
  }

Member Function Documentation

TGraphAsymmErrors * ana::EnsembleRatio::ErrorBand

(

)

const

Definition at line 22 of file EnsembleRatio.cxx.

  {
    std::unique_ptr<TH1D> hnom(fNom.ToTH1());

    std::vector<std::unique_ptr<TH1D>> hunivs;
    hunivs.reserve(fUnivs.size());
    for(const Ratio& u: fUnivs) hunivs.emplace_back(u.ToTH1());

    TGraphAsymmErrors* g = new TGraphAsymmErrors;

    for(int binIdx = 0; binIdx < hnom->GetNbinsX()+2; ++binIdx){
      const double xnom = hnom->GetXaxis()->GetBinCenter(binIdx);
      const double ynom = hnom->GetBinContent(binIdx);
      g->SetPoint(binIdx, xnom, ynom);

      const double dx = hnom->GetXaxis()->GetBinWidth(binIdx);

      std::vector<double> ys;
      ys.reserve(hunivs.size());
      for(const std::unique_ptr<TH1D>& hu: hunivs){
        ys.push_back(hu->GetBinContent(binIdx));
      }

      // 1 sigma
      const double y0 = FindQuantile(.5-0.6827/2, ys);
      const double y1 = FindQuantile(.5+0.6827/2, ys);

      // It's theoretically possible for the central value to be outside the
      // error bands - clamp to zero in that case
      g->SetPointError(binIdx, dx/2, dx/2,
                       std::max(ynom-y0, 0.),
                       std::max(y1-ynom, 0.));
    } // end for binIdx

    return g;
  }

Ratio ana::EnsembleRatio::Nominal ( ) const

inline

Definition at line 18 of file EnsembleRatio.h.

{return fNom;}

unsigned int ana::EnsembleRatio::NUniverses ( ) const

inline

Definition at line 19 of file EnsembleRatio.h.

{return fUnivs.size();}

EnsembleRatio ana::EnsembleRatio::operator*

(

const EnsembleRatio &

rhs

)

const

Definition at line 70 of file EnsembleRatio.cxx.

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

EnsembleRatio & ana::EnsembleRatio::operator*=

(

const EnsembleRatio &

rhs

)

Definition at line 60 of file EnsembleRatio.cxx.

  {
    assert(rhs.NUniverses() == fUnivs.size());
    fNom *= rhs.fNom;
    for(unsigned int i = 0; i < fUnivs.size(); ++i)
      fUnivs[i] *= rhs.fUnivs[i];
    return *this;
  }

EnsembleRatio ana::EnsembleRatio::operator/

(

const EnsembleRatio &

rhs

)

const

Definition at line 88 of file EnsembleRatio.cxx.

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

EnsembleRatio & ana::EnsembleRatio::operator/=

(

const EnsembleRatio &

rhs

)

Definition at line 78 of file EnsembleRatio.cxx.

  {
    assert(rhs.NUniverses() == fUnivs.size());
    fNom /= rhs.fNom;
    for(unsigned int i = 0; i < fUnivs.size(); ++i)
      fUnivs[i] /= rhs.fUnivs[i];
    return *this;
  }

Ratio ana::EnsembleRatio::Universe ( unsigned int  i ) const

inline

Definition at line 20 of file EnsembleRatio.h.

    {
      assert(i < fUnivs.size());
      return fUnivs[i];
    }

Member Data Documentation

Ratio ana::EnsembleRatio::fNom

protected

Definition at line 35 of file EnsembleRatio.h.

std::vector<Ratio> ana::EnsembleRatio::fUnivs

protected

Definition at line 36 of file EnsembleRatio.h.

---

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

• EnsembleRatio.h
• EnsembleRatio.cxx