ana::Hist Class Reference

#include "/cvmfs/nova.opensciencegrid.org/externals/cafanacore/v01.11/src/CAFAna/Core/Hist.h"

Public Member Functions
	Hist (const Hist &)
	Hist (Hist &&)
Hist &	operator= (const Hist &)
Hist &	operator= (Hist &&)
	~Hist ()
bool	Initialized () const
TH1D *	ToTH1 (const Binning &bins) const
bool	HasStan () const
const Eigen::ArrayXd &	GetEigen () const
const Eigen::ArrayXstan &	GetEigenStan () const
int	GetNbinsX () const
double	GetBinError (int i) const
double	Integral () const
void	Fill (const Binning &bins, double x, double w)
void	Scale (double s)
void	Scale (const stan::math::var &s)
void	ResetErrors ()
double	GetBinContent (int i) const
void	SetBinContent (int i, double x)
void	Reset ()
void	Add (const Hist &rhs, double scale=1)
void	Multiply (const Hist &rhs)
void	Divide (const Hist &rhs)
void	Write (const Binning &bins) const

Static Public Member Functions
static Hist	Uninitialized ()
static Hist	Zero (int nbins)
static Hist	ZeroSparse (int nbins)
static Hist	AdoptSparse (Eigen::SparseVector< double > &&v)
static Hist	AdoptStan (Eigen::ArrayXstan &&v)
static Hist	Adopt (Eigen::ArrayXd &&v)
static Hist	FromDirectory (TDirectory *dir)

Protected Types
enum	EType { kUninitialized, kDense, kDenseStan, kSparse }

Protected Member Functions
	Hist ()
void	Add (const Eigen::SparseVector< double > &rhs, double scale)
void	Add (const Eigen::ArrayXstan &rhs, double scale)
void	Add (const Eigen::ArrayXd &rhs, double scale)

Protected Attributes
EType	fType
Eigen::SparseVector< double >	fDataSparse
Eigen::ArrayXstan	fDataStan
Eigen::ArrayXd	fData
Eigen::ArrayXd	fSumSq
	Accumulate errors, if enabled. More...
bool	fSqrtErrs
	Special case when filled with unweighted data. More...

Detailed Description

Definition at line 29 of file Hist.h.

Member Enumeration Documentation

enum ana::Hist::EType

protected

Enumerator
kUninitialized
kDense
kDenseStan
kSparse

Definition at line 83 of file Hist.h.

{kUninitialized, kDense, kDenseStan, kSparse};

Constructor & Destructor Documentation

ana::Hist::Hist

(

const Hist &

rhs

)

Definition at line 61 of file Hist.cxx.

References ana::assert(), fData, fDataSparse, fDataStan, fSqrtErrs, fSumSq, fType, and Initialized().

    : Hist()
  {
    assert(rhs.Initialized());

    fType = rhs.fType;
    // Only one of these will actually have contents
    fDataSparse = rhs.fDataSparse;
    fDataStan   = rhs.fDataStan;
    fData       = rhs.fData;
    fSumSq      = rhs.fSumSq;
    fSqrtErrs   = rhs.fSqrtErrs;
  }

ana::Hist::Hist

(

Hist &&

rhs

)

Definition at line 76 of file Hist.cxx.

References ana::assert(), fData, fDataSparse, fDataStan, fSqrtErrs, fSumSq, fType, and std::swap().

    : Hist()
  {
    assert(rhs.Initialized());

    fType = rhs.fType;
    std::swap(fDataSparse, rhs.fDataSparse);
    std::swap(fDataStan,   rhs.fDataStan);
    std::swap(fData,       rhs.fData);
    std::swap(fSumSq,      rhs.fSumSq);
    fSqrtErrs = rhs.fSqrtErrs;
  }

ana::Hist::~Hist

(

)

Definition at line 126 of file Hist.cxx.

  {
  }

ana::Hist::Hist ( )

protected

Definition at line 36 of file Hist.cxx.

             : fType(kUninitialized), fSqrtErrs(false)
  {
  }

Member Function Documentation

void ana::Hist::Add	(	const Hist &	rhs,
		double	scale = 1
	)

Definition at line 455 of file Hist.cxx.

References ana::assert(), fData, fDataSparse, fDataStan, fSqrtErrs, fSumSq, fType, Initialized(), kDense, kDenseStan, and kSparse.

Referenced by ana::Spectrum::PlusEqualsHelper().

  {
    assert(Initialized());
    assert(rhs.Initialized());

    switch(rhs.fType){
    case kSparse:    Add(rhs.fDataSparse, scale); break;
    case kDenseStan: Add(rhs.fDataStan,   scale); break;
    case kDense:     Add(rhs.fData,       scale); break;
    default: abort(); // unreachable
    }

    if(scale != 1 || !rhs.fSqrtErrs) fSqrtErrs = false;

    if(fSumSq.size() > 0){
      if(rhs.fSumSq.size() > 0) fSumSq += scale * rhs.fSumSq;
      // otherwise nothing to add in
    }
    else{
      fSumSq = scale * rhs.fSumSq;
    }
  }

void ana::Hist::Add ( const Eigen::SparseVector< double > &  rhs, double  scale  )

protected

Definition at line 396 of file Hist.cxx.

References fData, fDataSparse, fDataStan, fType, kDense, kDenseStan, kSparse, and scale.

  {
    switch(fType){
    case kSparse:    fDataSparse += rhs * scale; break;
    case kDenseStan: fDataStan   += rhs * scale; break;
    case kDense:     fData       += rhs * scale; break;
    default: abort(); // unreachable
    }
  }

void ana::Hist::Add ( const Eigen::ArrayXstan &  rhs, double  scale  )

protected

Definition at line 407 of file Hist.cxx.

References fData, fDataSparse, fDataStan, fType, kDense, kDenseStan, kSparse, and scale.

  {
    switch(fType){
    case kSparse:
      fType = kDenseStan;
      fDataStan = rhs * scale;
      fDataStan += fDataSparse;
      fDataSparse.resize(0);
      break;

    case kDenseStan:
      fDataStan += rhs * scale;
      break;

    case kDense:
      fType = kDenseStan;
      fDataStan = fData + rhs * scale;
      fData.resize(0);
      break;

    default: abort(); // unreachable
    }
  }

void ana::Hist::Add ( const Eigen::ArrayXd &  rhs, double  scale  )

protected

Definition at line 432 of file Hist.cxx.

References fData, fDataSparse, fDataStan, fType, kDense, kDenseStan, kSparse, and scale.

  {
    switch(fType){
    case kSparse:
      fType = kDense;
      fData = rhs * scale;
      fData += fDataSparse;
      fDataSparse.resize(0);
      break;

    case kDenseStan:
      fDataStan += rhs * scale;
      break;

    case kDense:
      fData += rhs * scale;
      break;

    default: abort(); // unreachable
    }
  }

Hist ana::Hist::Adopt ( Eigen::ArrayXd &&  v )

static

Definition at line 149 of file Hist.cxx.

References fData, fType, kDense, runNovaSAM::ret, and registry_explorer::v.

Referenced by ana::ReweightableSpectrum::UnWeighted(), ana::ReweightableSpectrum::WeightedBy(), and ana::ReweightableSpectrum::WeightingVariable().

  {
    Hist ret;
    ret.fType = kDense;
    ret.fData = std::move(v);
    return ret;
  }

Hist ana::Hist::AdoptSparse ( Eigen::SparseVector< double > &&  v )

static

Definition at line 131 of file Hist.cxx.

References fDataSparse, fType, kSparse, runNovaSAM::ret, and registry_explorer::v.

  {
    Hist ret;
    ret.fType = kSparse;
    ret.fDataSparse = std::move(v);
    return ret;
  }

Hist ana::Hist::AdoptStan ( Eigen::ArrayXstan &&  v )

static

Definition at line 140 of file Hist.cxx.

References fDataStan, fType, kDenseStan, runNovaSAM::ret, and registry_explorer::v.

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

  {
    Hist ret;
    ret.fType = kDenseStan;
    ret.fDataStan = std::move(v);
    return ret;
  }

void ana::Hist::Divide

(

const Hist &

rhs

)

Definition at line 531 of file Hist.cxx.

References ana::assert(), om::cout, allTimeWatchdog::endl, fData, fDataStan, fSqrtErrs, fSumSq, fType, Initialized(), kDense, kDenseStan, kSparse, and util::sqr().

Referenced by ana::Ratio::operator/=(), ana::Spectrum::operator/=(), and ana::Ratio::Ratio().

  {
    assert(Initialized());
    assert(rhs.Initialized());

    if(fType == kSparse || rhs.fType == kSparse){
      std::cout << "Hist::Divide() not implemented for sparse vectors" << std::endl;
      abort();
    }

    if(fType == kDenseStan){
      if(rhs.fType == kDenseStan){
        fDataStan /= rhs.fDataStan;
      }
      else{
        fDataStan /= rhs.fData;
      }
    }
    else{
      if(rhs.fType == kDenseStan){
        fType = kDenseStan;
        fDataStan = fData / rhs.fDataStan;
        fData.resize(0);
      }
      else{
        fData /= rhs.fData;
      }
    }

    if(fType == kDense && rhs.fType == kDense){
      if(fSumSq.size() > 0){
        if(rhs.fSumSq.size() > 0){
          fSumSq = fSumSq / util::sqr(rhs.fData) + rhs.fSumSq * util::sqr(fData) / util::sqr(util::sqr(rhs.fData));
        }
        else{
          fSumSq /= util::sqr(rhs.fData);
        }
      }
      else if(rhs.fSumSq.size() > 0){
        fSumSq = rhs.fSumSq * util::sqr(fData) / util::sqr(util::sqr(rhs.fData));
      }
    }
    else{
      // Didn't bother to implement error prop for dividing stan/sparse
      // hists
      fSumSq.resize(0);
    }

    fSqrtErrs = false;
  }

void ana::Hist::Fill	(	const Binning &	bins,
		double	x,
		double	w
	)

Definition at line 265 of file Hist.cxx.

References ana::assert(), bin, om::cout, allTimeWatchdog::endl, fData, fDataSparse, ana::Binning::FindBin(), fSqrtErrs, fSumSq, fType, gStatErrs, Initialized(), kDense, kDenseStan, kSparse, and w.

Referenced by ana::Spectrum::Fill().

  {
    assert(Initialized());

    if(w != 1) fSqrtErrs = false;

    if(fType == kDense && gStatErrs && fSumSq.size() == 0) fSumSq.resize(fData.size());

    switch(fType){
    case kSparse:
      fDataSparse.coeffRef(bins.FindBin(x)) += w;
      break;
    case kDenseStan:
      std::cout << "Hist::Fill() not supported for stan vars" << std::endl;
      abort();
    case kDense:
      {
        const int bin = bins.FindBin(x);
        fData[bin] += w;
        if(gStatErrs){
          if(fSumSq.size() == 0) fSumSq.resize(fData.size());
          fSumSq[bin] += w*w;
        }
      }
      break;
    default:
      abort(); // unreachable
    }
  }

Hist ana::Hist::FromDirectory ( TDirectory *  dir )

static

Definition at line 158 of file Hist.cxx.

References ana::assert(), d, fData, fDataSparse, fSumSq, fType, gStatErrs, make_syst_table_plots::h, MECModelEnuComparisons::i, compare_h5_caf::idx, kDense, kSparse, nbins, runNovaSAM::ret, util::sqr(), and submit_syst::y.

Referenced by ana::Spectrum::LoadFrom().

  {
    Hist ret;

    TH1D* h = (TH1D*)dir->Get("hist");
    THnSparseD* hSparse = (THnSparseD*)dir->Get("hist_sparse");
    assert(bool(h) != bool(hSparse));

    if(h){
      ret.fType = kDense;
      ret.fData = Eigen::Map<Eigen::ArrayXd>(h->GetArray(), h->GetNbinsX()+2);

      if(gStatErrs){
        const int N = h->GetNbinsX()+2;
        ret.fSumSq.resize(N);
        for(int i = 0; i < N; ++i) ret.fSumSq[i] = util::sqr(h->GetBinError(i));
      }
    }
    if(hSparse){
      ret.fType = kSparse;

      // Number of bins there would be in a dense version
      int nbins = 1;
      const int N = hSparse->GetNdimensions();
      for(int d = 0; d < N; ++d)
        nbins *= hSparse->GetAxis(d)->GetNbins()+2;

      ret.fDataSparse.resize(nbins);
      ret.fDataSparse.setZero();
      for(int i = 0; i < hSparse->GetNbins(); ++i){
        int idx[N];
        const double y = hSparse->GetBinContent(i, idx);
        int idx1d = idx[0];
        int mult = 1;
        for(int d = 1; d < N; ++d){
          mult *= hSparse->GetAxis(d-1)->GetNbins()+2;
          idx1d += mult*idx[d];
        }
        ret.fDataSparse.coeffRef(idx1d) = y;
      }
    }

    delete h;
    delete hSparse;

    return ret;
  }

double ana::Hist::GetBinContent

(

int

i

)

const

Definition at line 349 of file Hist.cxx.

References ana::assert(), fData, fDataSparse, fDataStan, fType, MECModelEnuComparisons::i, Initialized(), kDense, kDenseStan, and kSparse.

Referenced by GetBinError(), ana::Spectrum::Mean(), ana::Spectrum::MockData(), and ana::Ratio::Ratio().

  {
    assert(Initialized());

    switch(fType){
    case kSparse: return fDataSparse.coeff(i);
    case kDenseStan: return fDataStan[i].val();
    case kDense: return fData[i];
    default:
      abort(); // unreachable
    }
  }

double ana::Hist::GetBinError

(

int

i

)

const

Definition at line 241 of file Hist.cxx.

References ana::assert(), fSqrtErrs, fSumSq, GetBinContent(), Initialized(), and std::sqrt().

Referenced by ana::Spectrum::Integral(), and ToTH1().

  {
    assert(Initialized());

    if(fSqrtErrs) return sqrt(GetBinContent(i));
    if(fSumSq.size() > 0) return sqrt(fSumSq[i]);

    return 0;
  }

const Eigen::ArrayXd& ana::Hist::GetEigen ( ) const

inline

Definition at line 53 of file Hist.h.

References ana::assert().

Referenced by ana::Ratio::GetEigen(), and ana::Spectrum::GetEigen().

{assert(fType == kDense); return fData;}

const Eigen::ArrayXstan& ana::Hist::GetEigenStan ( ) const

inline

Definition at line 54 of file Hist.h.

References ana::assert(), Divide(), Fill(), MECModelEnuComparisons::i, Integral(), slidt::Reset(), scale, Scale(), SetBinContent(), w, Write(), and submit_syst::x.

Referenced by ana::Ratio::GetEigenStan(), and ana::Spectrum::GetEigenStan().

{assert(fType == kDenseStan); return fDataStan;}

int ana::Hist::GetNbinsX

(

)

const

Definition at line 228 of file Hist.cxx.

References ana::assert(), fData, fDataSparse, fDataStan, fType, Initialized(), kDense, kDenseStan, and kSparse.

Referenced by ana::Spectrum::Integral(), ana::Spectrum::MockData(), and ana::Ratio::Ratio().

  {
    assert(Initialized());

    switch(fType){
    case kSparse:    return fDataSparse.size()-2;
    case kDenseStan: return fDataStan  .size()-2;
    case kDense:     return fData      .size()-2;
    default: abort(); // unreachable
    }
  }

bool ana::Hist::HasStan ( ) const

inline

Definition at line 52 of file Hist.h.

Referenced by ana::Ratio::HasStan(), and ana::Spectrum::HasStan().

{return fType == kDenseStan;}

bool ana::Hist::Initialized ( ) const

inline

Definition at line 39 of file Hist.h.

References bins, dir, nbins, registry_explorer::v, and Zero().

Referenced by Add(), Divide(), Fill(), GetBinContent(), GetBinError(), GetNbinsX(), Hist(), Integral(), Multiply(), operator=(), ana::Spectrum::PlusEqualsHelper(), Scale(), SetBinContent(), ToTH1(), and Write().

{return fType != kUninitialized;}

double ana::Hist::Integral

(

)

const

Definition at line 252 of file Hist.cxx.

References ana::assert(), fData, fDataSparse, fDataStan, fType, Initialized(), kDense, kDenseStan, and kSparse.

Referenced by ana::Spectrum::Integral(), and ana::Spectrum::PlusEqualsHelper().

  {
    assert(Initialized());

    switch(fType){
    case kSparse:    return fDataSparse.sum();
    case kDenseStan: return fDataStan  .sum().val();
    case kDense:     return fData      .sum();
    default: abort(); // unreachable
    }
  }

void ana::Hist::Multiply

(

const Hist &

rhs

)

Definition at line 479 of file Hist.cxx.

References ana::assert(), om::cout, allTimeWatchdog::endl, fData, fDataStan, fSqrtErrs, fSumSq, fType, Initialized(), kDense, kDenseStan, kSparse, and util::sqr().

Referenced by ana::Ratio::operator*=(), and ana::Spectrum::operator*=().

  {
    assert(Initialized());
    assert(rhs.Initialized());

    if(fType == kSparse || rhs.fType == kSparse){
      std::cout << "Hist::Multiply() not implemented for sparse vectors" << std::endl;
      abort();
    }

    if(fType == kDenseStan){
      if(rhs.fType == kDenseStan){
        fDataStan *= rhs.fDataStan;
      }
      else{
        fDataStan *= rhs.fData;
      }
    }
    else{
      if(rhs.fType == kDenseStan){
        fType = kDenseStan;
        fDataStan = fData * rhs.fDataStan;
        fData.resize(0);
      }
      else{
        fData *= rhs.fData;
      }
    }

    if(fType == kDense && rhs.fType == kDense){
      if(fSumSq.size() > 0){
        if(rhs.fSumSq.size() > 0){
          fSumSq = fSumSq * util::sqr(rhs.fData) + rhs.fSumSq * util::sqr(fData);
        }
        else{
          fSumSq *= util::sqr(rhs.fData);
        }
      }
      else if(rhs.fSumSq.size() > 0){
        fSumSq = rhs.fSumSq * util::sqr(fData);
      }
    }
    else{
      // Didn't bother to implement error prop for multiplying stan/sparse
      // hists
      fSumSq.resize(0);
    }

    fSqrtErrs = false;
  }

Hist & ana::Hist::operator=

(

const Hist &

rhs

)

Definition at line 90 of file Hist.cxx.

References ana::assert(), fData, fDataSparse, fDataStan, fSqrtErrs, fSumSq, fType, and Initialized().

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

    DontAddDirectory guard;

    assert(rhs.Initialized());

    fType = rhs.fType;
    fDataSparse = rhs.fDataSparse;
    fDataStan   = rhs.fDataStan;
    fData       = rhs.fData;
    fSumSq      = rhs.fSumSq;
    fSqrtErrs   = rhs.fSqrtErrs;

    return *this;
  }

Hist & ana::Hist::operator=

(

Hist &&

rhs

)

Definition at line 109 of file Hist.cxx.

References ana::assert(), fData, fDataSparse, fDataStan, fSqrtErrs, fSumSq, fType, and std::swap().

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

    assert(rhs.Initialized());

    fType = rhs.fType;
    std::swap(fDataSparse, rhs.fDataSparse);
    std::swap(fDataStan,   rhs.fDataStan);
    std::swap(fData,       rhs.fData);
    std::swap(fSumSq,      rhs.fSumSq);
    fSqrtErrs = rhs.fSqrtErrs;

    return *this;
  }

void ana::Hist::Reset

(

void

)

Definition at line 382 of file Hist.cxx.

References fData, fDataSparse, fDataStan, fSqrtErrs, fSumSq, fType, kDense, kDenseStan, and kSparse.

Referenced by ana::Spectrum::Clear().

  {
    switch(fType){
    case kSparse:    fDataSparse.setZero(); break;
    case kDenseStan: fDataStan  .setZero(); break;
    case kDense:     fData      .setZero(); break;
    default: ; // OK?
    }

    fSqrtErrs = true;
    fSumSq.resize(0);
  }

void ana::Hist::ResetErrors

(

)

Definition at line 342 of file Hist.cxx.

References fSqrtErrs, and fSumSq.

Referenced by ana::Spectrum::AsimovData(), and ana::Spectrum::MockData().

  {
    fSumSq.resize(0);
    fSqrtErrs = true;
  }

void ana::Hist::Scale

(

double

s

)

Definition at line 296 of file Hist.cxx.

References ana::assert(), fData, fDataSparse, fDataStan, fSqrtErrs, fSumSq, fType, Initialized(), kDense, kDenseStan, and kSparse.

Referenced by ana::Spectrum::AsimovData(), ana::Ratio::Ratio(), and ana::Spectrum::Scale().

  {
    assert(Initialized());

    switch(fType){
    case kSparse:    fDataSparse *= s; break;
    case kDenseStan: fDataStan   *= s; break;
    case kDense:     fData       *= s; break;
    default:
      abort(); // unreachable
    }

    if(s != 1) fSqrtErrs = false;
    fSumSq *= s;
  }

void ana::Hist::Scale

(

const stan::math::var &

s

)

Definition at line 313 of file Hist.cxx.

References ana::assert(), fData, fDataSparse, fDataStan, fSqrtErrs, fSumSq, fType, Initialized(), kDense, kDenseStan, kSparse, and stan::math::var::val().

  {
    assert(Initialized());
    switch (fType){
    case kSparse:
      fDataStan = fDataSparse * s;
      fDataSparse.resize(0);
      break;

    case kDense:
      fDataStan = fData * s;
      fData.resize(0);
      break;

    case kDenseStan:
      fDataStan *= s;
      break;

    default:
      abort(); // unreachable
    }

    fType = kDenseStan;

    if(s != 1) fSqrtErrs = false;
    fSumSq *= s.val();
  }

void ana::Hist::SetBinContent	(	int	i,
		double	x
	)

Definition at line 363 of file Hist.cxx.

References ana::assert(), om::cout, allTimeWatchdog::endl, fData, fDataSparse, fSqrtErrs, fSumSq, fType, MECModelEnuComparisons::i, Initialized(), kDense, kDenseStan, kSparse, and submit_syst::x.

Referenced by ana::Spectrum::MockData(), and ana::Ratio::Ratio().

  {
    assert(Initialized());

    switch(fType){
    case kSparse: fDataSparse.coeffRef(i) = x; break;
    case kDenseStan:
      std::cout << "Hist::SetBinContent() not implemented for stan vars" << std::endl;
      abort();
    case kDense: fData[i] = x; break;
    default:
      abort(); // unreachable
    }

    fSqrtErrs = false;
    if(fSumSq.size() > 0) fSumSq[i] = 0;
  }

TH1D * ana::Hist::ToTH1

(

const Binning &

bins

)

const

Definition at line 207 of file Hist.cxx.

References ana::assert(), fData, fDataSparse, fDataStan, fSqrtErrs, fSumSq, fType, GetBinError(), MECModelEnuComparisons::i, Initialized(), kDense, kDenseStan, kSparse, ana::MakeTH1D(), ana::Binning::NBins(), runNovaSAM::ret, ana::UniqueName(), and febshutoff_auto::val.

Referenced by ana::Ratio::ToTH1(), ana::Spectrum::ToTH1(), and Write().

  {
    assert(Initialized());

    TH1D* ret = MakeTH1D(UniqueName().c_str(), "", bins);

    for(int i = 0; i < bins.NBins()+2; ++i){
      switch(fType){
      case kDense:     ret->SetBinContent(i, fData[i]); break;
      case kDenseStan: ret->SetBinContent(i, fDataStan[i].val()); break;
        // Interface requires returning literally a TH1D in any case
      case kSparse:    ret->SetBinContent(i, fDataSparse.coeff(i)); break;
      default:
        abort(); // unreachable
      }
      if(fSumSq.size() > 0  || fSqrtErrs) ret->SetBinError(i, GetBinError(i));
    }
    return ret;
  }

static Hist ana::Hist::Uninitialized ( )

inlinestatic

Definition at line 38 of file Hist.h.

{return Hist();}

void ana::Hist::Write

(

const Binning &

bins

)

const

Definition at line 583 of file Hist.cxx.

References ana::assert(), om::cout, allTimeWatchdog::endl, fDataSparse, fSqrtErrs, fType, make_syst_table_plots::h, compare_h5_caf::idx, Initialized(), ana::Binning::IsSimple(), it, kDense, kDenseStan, kSparse, ana::Binning::Max(), ana::Binning::Min(), getGoodRuns4SAM::n, ana::Binning::NBins(), std::sqrt(), ToTH1(), and x1.

Referenced by ana::Spectrum::SaveTo().

  {
    assert(Initialized());

    if(fType == kDenseStan){
      std::cout << "Hist::Write() not implemented (impossible?) for stan vars" << std::endl;
      abort();
    }
    if(fType == kDense){
      TH1D* h = ToTH1(bins);
      h->Write("hist");
      delete h;
    }
    if(fType == kSparse){
      const int n = bins.NBins();
      const double x0 = bins.IsSimple() ? bins.Min() : 0;
      const double x1 = bins.IsSimple() ? bins.Max() : bins.NBins();
      THnSparseD* h = new THnSparseD("", "", 1, &n, &x0, &x1);

      for(Eigen::SparseVector<double>::InnerIterator it(fDataSparse); it; ++it){
        const int idx = it.index();
        h->SetBinContent(&idx, it.value());
        if(fSqrtErrs) h->SetBinError(&idx, sqrt(it.value()));
      }

      h->Write("hist_sparse");
      delete h;
    }
  }

Hist ana::Hist::Zero ( int  nbins )

static

Definition at line 41 of file Hist.cxx.

References fData, fSqrtErrs, fType, kDense, runNovaSAM::ret, and Zero().

Referenced by ana::Spectrum::Spectrum().

  {
    Hist ret;
    ret.fType = kDense;
    ret.fData = Eigen::ArrayXd::Zero(nbins+2);
    ret.fSqrtErrs = true;
    return ret;
  }

Hist ana::Hist::ZeroSparse ( int  nbins )

static

Definition at line 51 of file Hist.cxx.

References fDataSparse, fSqrtErrs, fType, kSparse, and runNovaSAM::ret.

Referenced by ana::Spectrum::Spectrum().

  {
    Hist ret;
    ret.fType = kSparse;
    ret.fDataSparse = Eigen::SparseVector<double>(nbins+2);
    ret.fSqrtErrs = true;
    return ret;
  }

Member Data Documentation

Eigen::ArrayXd ana::Hist::fData

protected

Definition at line 88 of file Hist.h.

Referenced by Add(), Adopt(), Divide(), Fill(), FromDirectory(), GetBinContent(), GetNbinsX(), Hist(), Integral(), Multiply(), operator=(), Reset(), Scale(), SetBinContent(), ToTH1(), and Zero().

Eigen::SparseVector<double> ana::Hist::fDataSparse

protected

Definition at line 86 of file Hist.h.

Referenced by Add(), AdoptSparse(), Fill(), FromDirectory(), GetBinContent(), GetNbinsX(), Hist(), Integral(), operator=(), Reset(), Scale(), SetBinContent(), ToTH1(), Write(), and ZeroSparse().

Eigen::ArrayXstan ana::Hist::fDataStan

protected

Definition at line 87 of file Hist.h.

Referenced by Add(), AdoptStan(), Divide(), GetBinContent(), GetNbinsX(), Hist(), Integral(), Multiply(), operator=(), Reset(), Scale(), and ToTH1().

bool ana::Hist::fSqrtErrs

protected

Special case when filled with unweighted data.

Definition at line 90 of file Hist.h.

Referenced by Add(), Divide(), Fill(), GetBinError(), Hist(), Multiply(), operator=(), Reset(), ResetErrors(), Scale(), SetBinContent(), ToTH1(), Write(), Zero(), and ZeroSparse().

Eigen::ArrayXd ana::Hist::fSumSq

protected

Accumulate errors, if enabled.

Definition at line 89 of file Hist.h.

Referenced by Add(), Divide(), Fill(), FromDirectory(), GetBinError(), Hist(), Multiply(), operator=(), Reset(), ResetErrors(), Scale(), SetBinContent(), and ToTH1().

EType ana::Hist::fType

protected

Definition at line 84 of file Hist.h.

Referenced by Add(), Adopt(), AdoptSparse(), AdoptStan(), Divide(), Fill(), FromDirectory(), GetBinContent(), GetNbinsX(), Hist(), Integral(), Multiply(), operator=(), Reset(), Scale(), SetBinContent(), ToTH1(), Write(), Zero(), and ZeroSparse().

---

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

• CAFAna/Core/Hist.h
• CAFAna/Core/Hist.cxx