ana::EnergyScaleSyst Class Reference

#include <EnergySysts.h>

Inheritance diagram for ana::EnergyScaleSyst:

Public Member Functions
	EnergyScaleSyst (EnergyScaleSystTerm _term, EnergyScaleSystParticle _part, EnergyScaleSystDetector _detector, double _uncertainty, const std::string &name, const std::string &latexName)
void	Shift (double sigma, 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

Private Attributes
EnergyScaleSystTerm	term
EnergyScaleSystParticle	part
EnergyScaleSystDetector	detector
double	uncertainty

Detailed Description

Definition at line 33 of file EnergySysts.h.

Constructor & Destructor Documentation

ana::EnergyScaleSyst::EnergyScaleSyst ( EnergyScaleSystTerm  _term, EnergyScaleSystParticle  _part, EnergyScaleSystDetector  _detector, double  _uncertainty, const std::string &  name, const std::string &  latexName  )

inline

Definition at line 36 of file EnergySysts.h.

                                                                                                                                                                                         : 
      ISyst(name, latexName), term(_term), part(_part), detector(_detector), uncertainty(_uncertainty) {}

Member Function Documentation

void ana::EnergyScaleSyst::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 6 of file EnergySysts.cxx.

  {
    double scale = uncertainty * sigma;
    auto& det = sr->truth.det;
    bool all = (detector == EnergyScaleSystDetector::kAll);
    bool nd = (detector == EnergyScaleSystDetector::kSBND && det == caf::kSBND);
    bool ub = (detector == EnergyScaleSystDetector::kMicroBooNE && det != caf::kSBND && det != caf::kICARUS);
    bool fd = (detector == EnergyScaleSystDetector::kICARUS && det == caf::kICARUS);
    bool detector_cut = all || nd || ub || fd;
    if(!sr->truth.iscc || abs(sr->truth.pdg) != 14 || isnan(sr->fake_reco.nuE) || !detector_cut) 
      return ;
    double particle_energy = 0.0;
    switch(part) {
    case EnergyScaleSystParticle::kMuon:
      particle_energy += sr->fake_reco.lepton.ke;
      break;
    case EnergyScaleSystParticle::kHadron:
      for(const auto& hadron: sr->fake_reco.hadrons) {
        particle_energy += hadron.ke;
      }
      break;
    case EnergyScaleSystParticle::kNeutron:
      for(const auto& hadron: sr->fake_reco.hadrons) {
        if(hadron.pid == 2112) {
          particle_energy += hadron.ke;
        }
      }
      break;
    case EnergyScaleSystParticle::kEM:
      for(const auto& hadron: sr->fake_reco.hadrons) {
        if(hadron.pid == 111) {
          particle_energy += hadron.ke;
        }
      }
      break;
   case EnergyScaleSystParticle::kChargedHadron:
      for(const auto& hadron: sr->fake_reco.hadrons) {
        auto pid = hadron.pid;
        if(pid == 2212 || abs(pid) == 211) {
          particle_energy += hadron.ke;
        }
      }
      break;
    }
    switch(term) {
    case EnergyScaleSystTerm::kConstant:
      break;
    case EnergyScaleSystTerm::kSqrt:
      scale *= std::sqrt(particle_energy);
      break;
    case EnergyScaleSystTerm::kInverseSqrt:
      scale /= std::sqrt(particle_energy + 0.1);
      break;
    }
    sr->fake_reco.nuE += particle_energy * scale;
   
  }

Member Data Documentation

EnergyScaleSystDetector ana::EnergyScaleSyst::detector

private

Definition at line 44 of file EnergySysts.h.

EnergyScaleSystParticle ana::EnergyScaleSyst::part

private

Definition at line 43 of file EnergySysts.h.

EnergyScaleSystTerm ana::EnergyScaleSyst::term

private

Definition at line 42 of file EnergySysts.h.

double ana::EnergyScaleSyst::uncertainty

private

Definition at line 45 of file EnergySysts.h.

---

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

• EnergySysts.h
• EnergySysts.cxx