ana::SBNWeightSyst Class Reference

#include <SBNWeightSysts.h>

Inheritance diagram for ana::SBNWeightSyst:

Classes
struct	Univs

Public Member Functions
	SBNWeightSyst (const std::string &systName)
void	Shift (double x, caf::SRSliceProxy *sr, double &weight) const override
	Perform the systematic shift. More...
Public Member Functions inherited from ana::ISyst
	ISyst (const std::string &shortName, const std::string &latexName, bool applyPenalty=true, double min=-3, double max=+3)
	ISyst (const ISyst &)=delete
	ISyst (ISyst &&rhs)=delete
virtual	~ISyst ()
void	operator= (const ISyst &)=delete
void	operator= (ISyst &&)=delete
virtual std::string	ShortName () const final
	The name printed out to the screen. More...
virtual std::string	LatexName () const final
	The name used on plots (ROOT's TLatex syntax) More...
virtual double	Penalty (double x) const
virtual bool	ApplyPenalty () const
	Should a penalty be applied for this shift? More...
virtual double	Min () const
	Return the min/max value for this syst. More...
virtual double	Max () const
virtual int	PredInterpMaxNSigma () const

Protected Member Functions
Univs	GetUnivs (double x) const

Protected Attributes
int	fIdx
std::unordered_map< double, Univs >	fUnivs

Detailed Description

Definition at line 33 of file SBNWeightSysts.h.

Constructor & Destructor Documentation

ana::SBNWeightSyst::SBNWeightSyst

(

const std::string &

systName

)

Definition at line 47 of file SBNWeightSysts.cxx.

    : ISyst(systName, systName),
      fIdx(-1)
  {
  }

Member Function Documentation

SBNWeightSyst::Univs ana::SBNWeightSyst::GetUnivs ( double  x ) const

protected

Definition at line 73 of file SBNWeightSysts.cxx.

  {
    auto it = fUnivs.find(x);
    if(it != fUnivs.end()) return it->second;

    Univs u;
    const UniverseOracle& uo = UniverseOracle::Instance();
    // Neighbours
    double x0, x1;
    u.i0 = uo.ClosestShiftIndex(ShortName(), x, ESide::kBelow, &x0);

    if(x0 == x){
      // Found an exact match (this is OK for small integer shifts, even with
      // floating point comparisons
      u.i1 = -1;
      u.w0 = 1;
      u.w1 = 0;
      fUnivs[x] = u;
      return u;
    }

    // Otherwise we're interpolating
    u.i1 = uo.ClosestShiftIndex(ShortName(), x, ESide::kAbove, &x1);
    // Interpolation weights
    u.w0 = (x1-x)/(x1-x0);
    u.w1 = (x-x0)/(x1-x0);

    //      std::cout << ShortName() << " " << x << " sigma, found indices " << u.i0 << " and " << u.i1 << " at " << x0 << " and " << x1 << ", will use weights " << u.w0 << " and " << u.w1 << std::endl;

    // If one of the neighbours wasn't found, we fall back to just using the
    // neighbour we did find. It would probably be better to find two
    // neighbours on the same side and extrapolate.
    if(u.i0 == -1){u.i0 = u.i1; u.w0 = u.w1 = 0.5;}
    if(u.i1 == -1){u.i1 = u.i0; u.w0 = u.w1 = 0.5;}

    fUnivs[x] = u;
    return u;
  }

void ana::SBNWeightSyst::Shift ( double  sigma, caf::SRSliceProxy *  sr, double &  weight  ) const

overridevirtual

Perform the systematic shift.

Parameters

sigma	Number of sigma to shift record by
sr	The record to inspect and alter
weight	Scale this weight for reweighting systematics

Implements ana::ISyst.

Definition at line 54 of file SBNWeightSysts.cxx.

  {
    if(sr->truth.index < 0) return;

    if(fIdx == -1) fIdx = UniverseOracle::Instance().SystIndex(ShortName());

    const caf::Proxy<std::vector<caf::SRMultiverse>>& wgts = sr->truth.wgt;
    if(wgts.empty()) return;

    const Univs u = GetUnivs(x);

    double y = 0;
    if(u.w0 != 0) y += u.w0 * wgts[fIdx].univ[u.i0];
    if(u.w1 != 0) y += u.w1 * wgts[fIdx].univ[u.i1];

    weight *= y;
  }

Member Data Documentation

int ana::SBNWeightSyst::fIdx

mutableprotected

Definition at line 41 of file SBNWeightSysts.h.

std::unordered_map<double, Univs> ana::SBNWeightSyst::fUnivs

mutableprotected

Definition at line 49 of file SBNWeightSysts.h.

---

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

• SBNWeightSysts.h
• SBNWeightSysts.cxx