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

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-11-25/3FlavorAna/Systs/EnergySysts2019.h"

Inheritance diagram for ana::PileupMuESyst2019:
ana::ISyst

Public Member Functions

 PileupMuESyst2019 (const std::string &shortName, const std::string &latexName, const double shiftND, const double shiftMuCat)
 
 PileupMuESyst2019 (const double shiftND, const double shiftMuCat)
 
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 fShiftND
 
const double fShiftMuCat
 

Detailed Description

Definition at line 138 of file EnergySysts2019.h.

Constructor & Destructor Documentation

ana::PileupMuESyst2019::PileupMuESyst2019 ( const std::string shortName,
const std::string latexName,
const double  shiftND,
const double  shiftMuCat 
)
inline

Definition at line 141 of file EnergySysts2019.h.

144  : ISyst(shortName, latexName), fShiftND(shiftND), fShiftMuCat(shiftMuCat) {}
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::PileupMuESyst2019::PileupMuESyst2019 ( const double  shiftND,
const double  shiftMuCat 
)
inline

Definition at line 146 of file EnergySysts2019.h.

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

147  : ISyst("PileupMuESyst2019", "Neutron Pile-up 2019"),
148  fShiftND(shiftND), fShiftMuCat(shiftMuCat) {}
ISyst(const std::string &shortName, const std::string &latexName)
Definition: ISyst.cxx:10

Member Function Documentation

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

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
inlinefinalvirtualinherited

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::PileupMuESyst2019::Shift ( double  sigma,
caf::SRProxy sr,
double &  weight 
) const
overridevirtual

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.

Parameters
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 76 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::kNEARDET, caf::Proxy< caf::SRNumuEnergy >::ndtrklenact, caf::Proxy< caf::SRNumuEnergy >::ndtrklencat, caf::Proxy< caf::SREnergyBranch >::numu, ana::BolombMuESyst2019::Shift(), sigma(), caf::Proxy< caf::SRKalman >::tracks, and caf::Proxy< caf::StandardRecord >::trk.

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

77  {
78  if(sr->trk.kalman.tracks.empty()) return;
79 
80  // This is a one-sided systematic, it can never *shorten* a track.
81  // Implementing a one-sided error like this works fine, as it leaves
82  // the objective function smooth and continuous.
83  if(sigma < 0) return;
84 
85  if(sr->hdr.det == caf::kNEARDET){
86  // Note subtle difference between "ndtrklenact" and "ndtrklencat"...
87  if(sr->energy.numu.ndtrklenact > 0 && sr->energy.numu.ndtrklencat == 0){
88  // Add for both ends of the track in the main ND
90  }
91  else if(sr->energy.numu.ndtrklenact == 0 && sr->energy.numu.ndtrklencat > 0){
92  // Add for both ends of the track in the Muon Catcher
94  }
95  else if(sr->energy.numu.ndtrklenact > 0 && sr->energy.numu.ndtrklencat > 0){
96  // Add for one end in the main detector and one in the Muon Catcher
99  }
100  // and of couse if both are zero because there's no track, do nothing
101  }
102  }
Near Detector underground.
Definition: SREnums.h:10
caf::Proxy< caf::SRHeader > hdr
Definition: SRProxy.h:2136
caf::Proxy< caf::SRNumuEnergy > numu
Definition: SRProxy.h:213
caf::Proxy< caf::SREnergyBranch > energy
Definition: SRProxy.h:2135
caf::Proxy< caf::SRTrackBranch > trk
Definition: SRProxy.h:2144
double sigma(TH1F *hist, double percentile)
caf::Proxy< caf::SRKalman > kalman
Definition: SRProxy.h:1796
caf::Proxy< float > ndtrklenact
Definition: SRProxy.h:185
caf::Proxy< std::vector< caf::SRKalmanTrack > > tracks
Definition: SRProxy.h:1779
caf::Proxy< float > ndtrklencat
Definition: SRProxy.h:186
caf::Proxy< caf::Det_t > det
Definition: SRProxy.h:231
virtual const std::string& ana::ISyst::ShortName ( ) const
inlinefinalvirtualinherited
virtual void ana::ISyst::TruthShift ( double  sigma,
caf::SRNeutrinoProxy nu,
double &  weight 
) const
inlinevirtualinherited

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  }

Member Data Documentation

const double ana::PileupMuESyst2019::fShiftMuCat
private

Definition at line 152 of file EnergySysts2019.h.

const double ana::PileupMuESyst2019::fShiftND
private

Definition at line 152 of file EnergySysts2019.h.


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