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

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-12-03/3FlavorAna/Systs/NueExtrapSysts.h"

Inheritance diagram for ana::NueExtrapSystBkg2017:
ana::ISyst

Public Member Functions

 NueExtrapSystBkg2017 ()
 
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

Definition at line 11 of file NueExtrapSysts.h.

Constructor & Destructor Documentation

ana::NueExtrapSystBkg2017::NueExtrapSystBkg2017 ( )
inline

Definition at line 14 of file NueExtrapSysts.h.

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

14 : ISyst("extrap_bkg", "Extrapolation Bkg") {}
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::NueExtrapSystBkg2017::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 22 of file NueExtrapSysts.cxx.

References abs(), caf::Proxy< caf::SRHeader >::det, caf::Proxy< caf::StandardRecord >::hdr, caf::kFARDET, ana::kNue2017FD, ana::kNue2017FDPeripheral, caf::Proxy< caf::StandardRecord >::mc, caf::Proxy< caf::SRTruthBranch >::nnu, caf::Proxy< caf::SRTruthBranch >::nu, and sigma().

Referenced by NueExtrapSystBkg2017(), and ana::NueExtrapSystSignalKin2017::NueExtrapSystSignalKin2017().

23  {
24  //Nue extrapolation systs of bkg from docdb-22478 ND subcomponents extrapolation
25  //maximal +-1.2 % in core sample and +-1.3 % in peripheral sample
26  const double core = .012;
27  const double peri = .013;
28 
29  //Check if in FD, otherwise left unaltered
30  if(sr->hdr.det != caf::kFARDET) return;
31 
32  //Check if neutrino, but not numu->nue signal, otherwise left unaltered
33  if(sr->mc.nnu == 0) return;
34  if(abs(sr->mc.nu[0].pdg) == 12 && abs(sr->mc.nu[0].pdgorig) == 14 && sr->mc.nu[0].iscc) return;
35 
36  if(sigma == 0) return;
37 
38  //Check whether in peripheral or core sample
39  if(kNue2017FDPeripheral(sr))
40  weight *= 1+peri*sigma;
41  else if(kNue2017FD(sr))
42  weight *= 1+core*sigma;
43  }
Far Detector at Ash River.
Definition: SREnums.h:11
const Var weight
caf::Proxy< caf::SRHeader > hdr
Definition: SRProxy.h:2137
caf::Proxy< std::vector< caf::SRNeutrino > > nu
Definition: SRProxy.h:618
void abs(TH1 *hist)
caf::Proxy< short int > nnu
Definition: SRProxy.h:617
const Cut kNue2017FD
Full FD nue selection for 2017 analysis.
Definition: NueCuts2017.h:130
const Cut kNue2017FDPeripheral
Definition of the FD peripheral sample around our core selected events.
Definition: NueCuts2017.h:152
double sigma(TH1F *hist, double percentile)
caf::Proxy< caf::SRTruthBranch > mc
Definition: SRProxy.h:2138
caf::Proxy< caf::Det_t > det
Definition: SRProxy.h:232
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  }

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