ana::SystShifts Class Reference

Simple record of shifts applied to systematic parameters. More...

#include <SystShifts.h>

Public Member Functions
	SystShifts ()
	SystShifts (const ISyst *syst, double shift)
	SystShifts (const std::map< const ISyst *, double > &shifts)
bool	IsNominal () const
void	SetShift (const ISyst *syst, double shift)
	Shifts are 0=nominal, -1,+1 = 1 sigma shifts. More...
double	GetShift (const ISyst *syst) const
void	ResetToNominal ()
double	Penalty () const
	Penalty term for chi-squared fits. More...
void	Shift (caf::SRSliceProxy *slc, double &weight) const
std::string	ShortName () const
	Brief description of component shifts, for printing to screen. More...
std::string	LatexName () const
	Long description of component shifts, for plot labels. More...
int	ID () const
	SystShifts with the same set of systs should have the same ID. More...
std::vector< const ISyst * >	ActiveSysts () const
void	SaveTo (TDirectory *dir) const

Static Public Member Functions
static SystShifts	Nominal ()
static SystShifts	RandomThrow (const std::vector< const ISyst * > &systs)
static std::unique_ptr < SystShifts >	LoadFrom (TDirectory *dir)

Protected Attributes
std::unordered_map< const ISyst *, double >	fSysts
int	fID

Static Protected Attributes
static int	fgNextID = 1
	The next unused ID. More...

Detailed Description

Simple record of shifts applied to systematic parameters.

Definition at line 16 of file SystShifts.h.

Constructor & Destructor Documentation

ana::SystShifts::SystShifts

(

)

Definition at line 20 of file SystShifts.cxx.

                         : fID(0)
  {
  }

ana::SystShifts::SystShifts	(	const ISyst *	syst,
		double	shift
	)

Definition at line 25 of file SystShifts.cxx.

    : fID(fgNextID++)
  {
    if(shift != 0) fSysts.emplace(syst, shift);
  }

ana::SystShifts::SystShifts

(

const std::map< const ISyst *, double > &

shifts

)

Definition at line 32 of file SystShifts.cxx.

    : fID(fgNextID++)
  {
    for(auto it: shifts) if(it.second != 0) fSysts.emplace(it.first, it.second);
  }

Member Function Documentation

std::vector< const ISyst * > ana::SystShifts::ActiveSysts

(

)

const

Definition at line 120 of file SystShifts.cxx.

  {
    std::vector<const ISyst*> ret;
    for(auto it: fSysts) ret.push_back(it.first);
    return ret;
  }

double ana::SystShifts::GetShift

(

const ISyst *

syst

)

const

Definition at line 56 of file SystShifts.cxx.

  {
    assert(syst);

    auto it = fSysts.find(syst);
    return (it == fSysts.end()) ? 0 : it->second;
  }

int ana::SystShifts::ID ( ) const

inline

SystShifts with the same set of systs should have the same ID.

Definition at line 44 of file SystShifts.h.

{return fID;}

bool ana::SystShifts::IsNominal ( ) const

inline

Definition at line 26 of file SystShifts.h.

{return fSysts.empty();}

std::string ana::SystShifts::LatexName

(

)

const

Long description of component shifts, for plot labels.

Definition at line 106 of file SystShifts.cxx.

  {
    if(fSysts.empty()) return "Nominal";

    std::string ret;
    for(auto it: fSysts){
      if(!ret.empty()) ret += ", ";
      ret += it.first->LatexName() + TString::Format(" = %+g", it.second).Data();
    }

    return ret;
  }

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

static

Definition at line 151 of file SystShifts.cxx.

  {
    TObjString* tag = (TObjString*)dir->Get("type");
    assert(tag);
    assert(tag->GetString() == "SystShifts");

    auto ret = std::make_unique<SystShifts>();

    TH1* h = (TH1*)dir->Get("vals");
    if(h){ // no histogram means nominal
      for(int i = 1; i <= h->GetNbinsX(); ++i){
        ret->SetShift(SystRegistry::ShortNameToSyst(h->GetXaxis()->GetBinLabel(i)),
                      h->GetBinContent(i));
      }
    }

    return ret;
  }

static SystShifts ana::SystShifts::Nominal ( )

inlinestatic

Definition at line 23 of file SystShifts.h.

{return SystShifts();}

double ana::SystShifts::Penalty

(

)

const

Penalty term for chi-squared fits.

Definition at line 73 of file SystShifts.cxx.

  {
    double ret = 0;
    // Systematics are all expressed in terms of sigmas
    for(auto it: fSysts) {
      // Only apply a penalty for systematics where this has been requested
      if (it.first->ApplyPenalty())
        ret += util::sqr(it.second);
    }
    return ret;
  }

SystShifts ana::SystShifts::RandomThrow ( const std::vector< const ISyst * > &  systs )

static

Definition at line 39 of file SystShifts.cxx.

  {
    SystShifts ret;
    for(const ISyst* s: systs) ret.SetShift(s, gRandom->Gaus(0, 1));
    return ret;
  }

void ana::SystShifts::ResetToNominal

(

)

Definition at line 65 of file SystShifts.cxx.

  {
    fID = 0;

    fSysts.clear();
  }

void ana::SystShifts::SaveTo

(

TDirectory *

dir

)

const

Definition at line 128 of file SystShifts.cxx.

  {
    TDirectory* tmp = gDirectory;

    dir->cd();
    TObjString("SystShifts").Write("type");

    // Don't write any histogram for the nominal case
    if(!fSysts.empty()){
      TH1D h("", "", fSysts.size(), 0, fSysts.size());
      int ibin = 0;
      for(auto it: fSysts){
        ++ibin;
        h.GetXaxis()->SetBinLabel(ibin, it.first->ShortName().c_str());
        h.SetBinContent(ibin, it.second);
      }
      h.Write("vals");
    }

    tmp->cd();
  }

void ana::SystShifts::SetShift	(	const ISyst *	syst,
		double	shift
	)

Shifts are 0=nominal, -1,+1 = 1 sigma shifts.

Definition at line 47 of file SystShifts.cxx.

  {
    fID = fgNextID++;

    fSysts.erase(syst);
    if(shift != 0) fSysts.emplace(syst, shift);
  }

void ana::SystShifts::Shift	(	caf::SRSliceProxy *	slc,
		double &	weight
	)		const

Definition at line 86 of file SystShifts.cxx.

  {
    for(auto it: fSysts) it.first->Shift(it.second, sr, weight);
  }

std::string ana::SystShifts::ShortName

(

)

const

Brief description of component shifts, for printing to screen.

Definition at line 92 of file SystShifts.cxx.

  {
    if(fSysts.empty()) return "nominal";

    std::string ret;
    for(auto it: fSysts){
      if(!ret.empty()) ret += ",";
      ret += it.first->ShortName() + TString::Format("=%+g", it.second).Data();
    }

    return ret;
  }

Member Data Documentation

int ana::SystShifts::fgNextID = 1

staticprotected

The next unused ID.

Definition at line 58 of file SystShifts.h.

int ana::SystShifts::fID

protected

Definition at line 56 of file SystShifts.h.

std::unordered_map<const ISyst*, double> ana::SystShifts::fSysts

protected

Definition at line 54 of file SystShifts.h.

---

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

• SystShifts.h
• SystShifts.cxx