ana::NueAcceptSystBkg2018FHC Class Reference

FHC systs. More...

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/3FlavorAna/Systs/NueAcceptSysts.h"

Inheritance diagram for ana::NueAcceptSystBkg2018FHC:

Public Member Functions
	NueAcceptSystBkg2018FHC ()
void	Shift (double sigma, caf::SRProxy *sr, double &weight) const override
	FHC implementation. More...
virtual const std::string &	ShortName () const final
	The name printed out to the screen. More...
virtual const std::string &	LatexName () 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

FHC systs.

Definition at line 14 of file NueAcceptSysts.h.

Constructor & Destructor Documentation

ana::NueAcceptSystBkg2018FHC::NueAcceptSystBkg2018FHC ( )

inline

Definition at line 17 of file NueAcceptSysts.h.

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

: ISyst("accept_bkg_FHC", "Acceptance Bkg FHC") {}

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.

{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().

{return fLatexName;}

void ana::NueAcceptSystBkg2018FHC::Shift ( double  sigma, caf::SRProxy *  sr, double &  weight  ) const

overridevirtual

FHC implementation.

FHC values

Reimplemented from ana::ISyst.

Definition at line 34 of file NueAcceptSysts.cxx.

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

Referenced by NueAcceptSystBkg2018FHC(), ana::NueAcceptSystBkg2018RHC::NueAcceptSystBkg2018RHC(), ana::NueAcceptSystSignalKin2018FHC::NueAcceptSystSignalKin2018FHC(), ana::NueAcceptSystSignalKin2018RHC::NueAcceptSystSignalKin2018RHC(), ana::NueAcceptSystSignalKin2020FHC::NueAcceptSystSignalKin2020FHC(), and ana::NueAcceptSystSignalKin2020RHC::NueAcceptSystSignalKin2020RHC().

  {
    // Nue acceptance systs of bkg from docdb-27935 ND subcomponents extrapolation
    // maximal core / peripheral FHC (RHC) difference: +- 1.5% / +-1.2% (+- 4.1% / +- 2.3%)

    /// FHC values
    const double core = .015;
    const double peri = .012;

    // Check if FHC, otherwise left unaltered
    if(kIsRHC(sr)) return;

    // Check if in FD, otherwise left unaltered
    if(sr->hdr.det != caf::kFARDET) return;

    // Check if neutrino, but not numu->nue signal, otherwise left unaltered
    if(sr->mc.nnu == 0) return;
    if(abs(sr->mc.nu[0].pdg) == 12 && abs(sr->mc.nu[0].pdgorig) == 14 && sr->mc.nu[0].iscc) return;
    // Left unaltered if tau neutrino
    if(abs(sr->mc.nu[0].pdg) == 16 && sr->mc.nu[0].iscc) return;

    if(sigma == 0) return;

    double kFDBkgCorr = core;

    // Check whether in peripheral or core sample
    if(kNue2018FDPeripheral(sr))
      kFDBkgCorr = peri;

    weight *= 1+kFDBkgCorr*sigma;
  }

virtual const std::string& ana::ISyst::ShortName ( ) const

inlinefinalvirtualinherited

The name printed out to the screen.

Definition at line 27 of file ISyst.h.

References ana::ISyst::fShortName.

Referenced by ana::PredictionSystJoint2018::AddNormSyst(), ana::PredictionSystJointDemo::AddNormSyst(), ana::PredictionSyst3Flavor2020::AddNormSyst(), ana::CovMxManager::AddSystematic(), ana::NuISyst::BaseName(), ana::PredictionInterp::DebugPlot(), DrawBins(), DrawPulls(), DrawSystShifts(), GetGENIEShiftLabels(), ana::SystShifts::GetShift(), ana::Nus17SystFromHist::LoadHists(), ana::NusSystFromHist::LoadHists(), ana::Nus18SystFromHist::LoadHists(), MakeCovarSim(), MakeCovMx(), MakeGENIEMatrix(), plot_3flavor_withsysts(), ana::PredictionSystJoint2018::PredictionSystJoint2018(), ana::PredictionSystJointDemo::PredictionSystJointDemo(), ana::MCMCSamples::SampleValue(), ana::NuISyst::SaveTo(), ana::SystShifts::SetShift(), ana::SortSystsName(), SystsGENIEAna(), test_stanfit_withsysts(), ana::MCMCSample::Val(), and ana::MCMCSamples::VarOffset().

{return fShortName;}

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

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

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/3FlavorAna/Systs/NueAcceptSysts.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/3FlavorAna/Systs/NueAcceptSysts.cxx