ana::UniverseOracle Class Reference

#include <UniverseOracle.h>

Public Member Functions
bool	SystExists (const std::string &name) const
std::vector< std::string >	Systs () const
	List of all known syst names. More...
const std::vector< float > &	ShiftsForSyst (const std::string &name) const
	List of shifts for this syst in each universe. More...
unsigned int	ParameterSetIndex (const std::string &name) const
	Which index in the weights array corresponds to this parameter set? More...
unsigned int	SystIndex (const std::string &name) const
	Which index in the weights array corresponds to the shifting of just this syst (in any parameter set)? More...
unsigned int	ClosestShiftIndex (const std::string &name, double shift, ESide side=ESide::kEither, double *trueShift=0) const
	Within that entry, which index corresponds most closely to 'shift'? More...

Static Public Member Functions
static UniverseOracle &	Instance ()

Protected Member Functions
	UniverseOracle ()

Protected Attributes
std::map< std::string, unsigned int >	fPSetIdxs
std::map< std::string, unsigned int >	fSystIdxs
std::map< std::string, std::vector< float > >	fShiftVals

Detailed Description

Definition at line 16 of file UniverseOracle.h.

Constructor & Destructor Documentation

ana::UniverseOracle::UniverseOracle ( )

protected

Definition at line 75 of file UniverseOracle.cxx.

  {
    const caf::SRGlobal global = GetSRGlobal();
    std::cout << "\nSystematic weights in file:" << std::endl;
    PrintSRGlobal(global);

    for(unsigned int i = 0; i < global.wgts.size(); ++i){
      const caf::SRWeightPSet& pset = global.wgts[i];

      // Save the pset index in all cases
      fPSetIdxs[pset.name] = i;

      // Only save the remaining fields in parameter sets where only a single
      // knob is shifted
      if(pset.map.size() != 1) continue;

      // Save which position in the vector this was
      fSystIdxs[pset.map[0].param.name] = i;
      // Save all the knob values
      fShiftVals[pset.map[0].param.name] = pset.map[0].vals;
    }
  }

Member Function Documentation

unsigned int ana::UniverseOracle::ClosestShiftIndex	(	const std::string &	name,
		double	shift,
		ESide	side = ESide::kEither,
		double *	trueShift = 0
	)		const

Within that entry, which index corresponds most closely to 'shift'?

Definition at line 143 of file UniverseOracle.cxx.

  {
    const std::vector<float>& v = ShiftsForSyst(name);
    int bestIdx = -1;
    double bestDist;
    for(unsigned int i = 0; i < v.size(); ++i){
      if(side == ESide::kBelow && v[i] > shift) continue;
      if(side == ESide::kAbove && v[i] < shift) continue;
      const double dv = fabs(v[i]-shift);
      if(bestIdx == -1 || dv < bestDist){
        bestIdx = i;
        bestDist = dv;
        if(trueShift) *trueShift = v[i];
      }
    }
    return unsigned(bestIdx);
  }

UniverseOracle & ana::UniverseOracle::Instance ( )

static

Definition at line 68 of file UniverseOracle.cxx.

  {
    static UniverseOracle uo;
    return uo;
  }

unsigned int ana::UniverseOracle::ParameterSetIndex

(

const std::string &

name

)

const

Which index in the weights array corresponds to this parameter set?

Definition at line 121 of file UniverseOracle.cxx.

  {
    auto it = fPSetIdxs.find(name);
    if(it == fPSetIdxs.end()){
      std::cout << "UniverseOracle: pset '" << name << "' not known" << std::endl;
      abort();
    }
    return it->second;
  }

const std::vector< float > & ana::UniverseOracle::ShiftsForSyst

(

const std::string &

name

)

const

List of shifts for this syst in each universe.

Definition at line 114 of file UniverseOracle.cxx.

  {
    assert(SystExists(name));
    return fShiftVals.find(name)->second;
  }

bool ana::UniverseOracle::SystExists

(

const std::string &

name

)

const

Definition at line 99 of file UniverseOracle.cxx.

  {
    return fShiftVals.find(name) != fShiftVals.end();
  }

unsigned int ana::UniverseOracle::SystIndex

(

const std::string &

name

)

const

Which index in the weights array corresponds to the shifting of just this syst (in any parameter set)?

Definition at line 132 of file UniverseOracle.cxx.

  {
    auto it = fSystIdxs.find(name);
    if(it == fSystIdxs.end()){
      std::cout << "UniverseOracle: syst '" << name << "' not known" << std::endl;
      abort();
    }
    return it->second;
  }

std::vector< std::string > ana::UniverseOracle::Systs

(

)

const

List of all known syst names.

Definition at line 105 of file UniverseOracle.cxx.

  {
    std::vector<std::string> ret;
    ret.reserve(fShiftVals.size());
    for(auto it: fShiftVals) ret.push_back(it.first);
    return ret;
  }

Member Data Documentation

std::map<std::string, unsigned int> ana::UniverseOracle::fPSetIdxs

protected

Definition at line 42 of file UniverseOracle.h.

std::map<std::string, std::vector<float> > ana::UniverseOracle::fShiftVals

protected

Definition at line 44 of file UniverseOracle.h.

std::map<std::string, unsigned int> ana::UniverseOracle::fSystIdxs

protected

Definition at line 43 of file UniverseOracle.h.

---

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

• UniverseOracle.h
• UniverseOracle.cxx