Public Member Functions | List of all members
ana::NumuRelEnergyScaleSyst Class Reference

#include "/cvmfs/"

Inheritance diagram for ana::NumuRelEnergyScaleSyst:

Public Member Functions

 NumuRelEnergyScaleSyst ()
void Shift (double sigma, caf::SRProxy *sr, double &weight) const override
 Perform the systematic shift. More...
virtual const std::stringShortName () const final
 The name printed out to the screen. More...
virtual const std::stringLatexName () const final
 The name used on plots (ROOT's TLatex syntax) More...
virtual void TruthShift (double sigma, caf::SRNeutrinoProxy *nu, double &weight) const
virtual bool IsGenieReweight () const
 GENIE reweights can only provide +/-1,2sigma. More...

Detailed Description

Relative energy scale systematic see table 2 in DocDB 13641

Definition at line 139 of file NumuSysts.h.

Constructor & Destructor Documentation

ana::NumuRelEnergyScaleSyst::NumuRelEnergyScaleSyst ( )

Definition at line 142 of file NumuSysts.h.

References ana::getAllDirectNumuSysts2017(), ana::getAllDirectNumuSysts2018(), ana::kNumukRelEnergyScaleSyst, ana::NumuGEANTNormSyst::Shift(), sigma(), sr, and ana::weight.

142 : ISyst("RelEScale", "Relative Energy Scale") {}
ISyst(const std::string &shortName, const std::string &latexName)
Definition: ISyst.cxx:10

Member Function Documentation

virtual bool ana::ISyst::IsGenieReweight ( ) const

GENIE reweights can only provide +/-1,2sigma.

Reimplemented in ana::SummedSyst.

Definition at line 56 of file ISyst.h.

56 {return false;}
virtual const std::string& ana::ISyst::LatexName ( ) const

The name used on plots (ROOT's TLatex syntax)

Definition at line 30 of file ISyst.h.

References ana::ISyst::fLatexName, ana::ISyst::Shift(), sigma(), sr, and ana::weight.

Referenced by ana::PredictionInterp::DebugPlotColz(), GetGENIEShiftLabels(), ana::NuISyst::SaveTo(), SystsGENIEAna(), and WriteSystName().

30 {return fLatexName;}
std::string fLatexName
Definition: ISyst.h:60
void ana::NumuRelEnergyScaleSyst::Shift ( double  sigma,
caf::SRProxy sr,
double &  weight 
) const

Perform the systematic shift.

Override this function if your systematic depends on non-SRNeutrino quantities. If it is SRNeutrino-only, implement the other function, and let this default forward to you when necessary.

sigmaNumber of sigma to shift record by
srThe record to inspect and alter
weightScale this weight for reweighting systematics

5% difference between detectors

Reimplemented from ana::ISyst.

Definition at line 98 of file NumuSysts.cxx.

References caf::Proxy< caf::SRSlice >::calE, caf::Proxy< caf::SRHeader >::det, caf::Proxy< caf::StandardRecord >::energy, caf::Proxy< caf::StandardRecord >::hdr, caf::kFARDET, caf::Proxy< caf::SREnergyBranch >::numu, scale, caf::Proxy< caf::StandardRecord >::slc, caf::Proxy< caf::SRNumuEnergy >::trkccE, caf::Proxy< caf::SRNumuEnergy >::trknonqeE, and caf::Proxy< caf::SRNumuEnergy >::trkqeE.

Referenced by ana::NumuEnergyScaleSyst::Shift().

99  {
100  if (sr->hdr.det != caf::kFARDET) return;
102  const double scale = 1 + .025*sigma /// 5% difference between detectors
103  *( (sr->hdr.det==caf::kFARDET) ? 1 : -1 );
105  sr->slc.calE *= scale;
106  sr->energy.numu.trknonqeE *= scale;
107  sr->energy.numu.trkqeE *= scale;
108  sr->energy.numu.trkccE *= scale;
109  }
caf::Proxy< float > trkccE
Definition: SRProxy.h:195
Far Detector at Ash River.
Definition: SREnums.h:11
caf::Proxy< float > trknonqeE
Definition: SRProxy.h:196
caf::Proxy< caf::SRHeader > hdr
Definition: SRProxy.h:2137
caf::Proxy< caf::SRNumuEnergy > numu
Definition: SRProxy.h:214
caf::Proxy< caf::SREnergyBranch > energy
Definition: SRProxy.h:2136
Double_t scale
Definition: plot.C:25
caf::Proxy< float > trkqeE
Definition: SRProxy.h:197
double sigma(TH1F *hist, double percentile)
caf::Proxy< caf::SRSlice > slc
Definition: SRProxy.h:2142
caf::Proxy< float > calE
Definition: SRProxy.h:1292
caf::Proxy< caf::Det_t > det
Definition: SRProxy.h:232
virtual const std::string& ana::ISyst::ShortName ( ) const
virtual void ana::ISyst::TruthShift ( double  sigma,
caf::SRNeutrinoProxy nu,
double &  weight 
) const

For systematics that deal only with the neutrino truth and not any reconstruction/PID details. Systematics defined this way will work on nuTree-derived spectra too (e.g. denominators of efficiencies).

Reimplemented in demo::DemoSyst1, ana::BeamSyst, ana::GenericSystComponentScale< T >, ana::GenericSystComponentScale< T >, ana::GenericSystComponentScale< T >, and ana::NOvARwgtSyst.

Definition at line 46 of file ISyst.h.

Referenced by ana::ISyst::Shift().

49  {
50  // Implement this function if your systematic depends only
51  // SRNeutrino. Left blank by default, since systematics using other
52  // information can do nothing sensible to the nuTree.
53  }

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