Public Member Functions | Private Attributes | List of all members
ana::BolombMuESyst2019 Class Reference

#include "/cvmfs/"

Inheritance diagram for ana::BolombMuESyst2019:

Public Member Functions

 BolombMuESyst2019 (const std::string &shortName, const std::string &latexName, const double shift)
 BolombMuESyst2019 (const double shift)
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...

Private Attributes

const double fShift

Detailed Description

Definition at line 161 of file EnergySysts2019.h.

Constructor & Destructor Documentation

ana::BolombMuESyst2019::BolombMuESyst2019 ( const std::string shortName,
const std::string latexName,
const double  shift 

Definition at line 164 of file EnergySysts2019.h.

167  : ISyst(shortName, latexName), fShift(shift) {}
string shortName
THUMBNAIL BLOCK: We need to make a thumbnail for each.
ISyst(const std::string &shortName, const std::string &latexName)
Definition: ISyst.cxx:10
ana::BolombMuESyst2019::BolombMuESyst2019 ( const double  shift)

Definition at line 169 of file EnergySysts2019.h.

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

170  : ISyst("BolombMuESyst2019", "Bodek and Coulomb 2019"),
171  fShift(shift) {}
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::BolombMuESyst2019::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

Reimplemented from ana::ISyst.

Definition at line 106 of file EnergySysts2019.cxx.

References caf::Proxy< caf::SRHeader >::det, caf::Proxy< caf::StandardRecord >::energy, caf::Proxy< caf::StandardRecord >::hdr, caf::Proxy< caf::SRTrackBranch >::kalman, caf::kFARDET, caf::kNEARDET, std::max(), caf::Proxy< caf::SRNumuEnergy >::ndtrklenact, caf::Proxy< caf::SREnergyBranch >::numu, caf::Proxy< caf::SRKalman >::tracks, and caf::Proxy< caf::StandardRecord >::trk.

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

107  {
108  if(sr->trk.kalman.tracks.empty()) return;
110  // This is a one-sided systematic, it can never *lengthen* a track.
111  if(sigma < 0) return;
113  // The systematic has to be applied as a length, not as an energy, so
114  // here's the conversion constant, which is precise enough for these
115  // purposes. It is as per plots in doc-26649, discarding subtlety. This
116  // is the full physical length in the detector, averaged over air and
117  // whatnot. Since this correction pertains to the interaction vertex, we
118  // don't need to worry about the Muon Catcher.
119  const double cm_per_mev = 0.44;
121  const double cmshift = fShift * cm_per_mev;
123  // In principle, this error applies equally to nu_e and nu_tau, but no
124  // attempt is made here to adjust anything but the muon.
125  if(sr->hdr.det == caf::kNEARDET){
126  sr->energy.numu.ndtrklenact =
127  std::max(sr->energy.numu.ndtrklenact + cmshift*sigma, 0.);
128  }
129  else if(sr->hdr.det == caf::kFARDET){
130  sr->trk.kalman.tracks[0].len =
131  std::max(sr->trk.kalman.tracks[0].len + cmshift*sigma, 0.);
132  }
133  }
Near Detector underground.
Definition: SREnums.h:10
T max(const caf::Proxy< T > &a, T b)
Far Detector at Ash River.
Definition: SREnums.h:11
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
caf::Proxy< caf::SRTrackBranch > trk
Definition: SRProxy.h:2145
double sigma(TH1F *hist, double percentile)
caf::Proxy< caf::SRKalman > kalman
Definition: SRProxy.h:1797
caf::Proxy< float > ndtrklenact
Definition: SRProxy.h:186
caf::Proxy< std::vector< caf::SRKalmanTrack > > tracks
Definition: SRProxy.h:1780
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 ana::BeamSyst, demo::DemoSyst1, ana::GenericSystComponentScale< T >, ana::GenericSystComponentScale< T >, ana::GenericSystComponentScale< T >, ana::ReinteractionSyst, 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  }

Member Data Documentation

const double ana::BolombMuESyst2019::fShift

Definition at line 175 of file EnergySysts2019.h.

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