ana::MichelTaggingSyst2020 Class Reference

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

Inheritance diagram for ana::MichelTaggingSyst2020:

Public Member Functions
	MichelTaggingSyst2020 ()
void	Shift (double sigma, caf::SRProxy *sr, double &weight) const override
	Perform the systematic shift. 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...

Private Member Functions
int	Hash (caf::SRProxy *sr) const

Private Attributes
double	fProbUP = .08
double	fProbDO = .05
TRandom3	fMRnd

Detailed Description

Definition at line 10 of file MichelTaggingSyst.h.

Constructor & Destructor Documentation

ana::MichelTaggingSyst2020::MichelTaggingSyst2020 ( )

inline

Definition at line 13 of file MichelTaggingSyst.h.

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

: ISyst("michel_tagging2020", "Michel Electrons Tagging Uncertainty") {}

Member Function Documentation

int ana::MichelTaggingSyst2020::Hash ( caf::SRProxy *  sr ) const

private

Definition at line 57 of file MichelTaggingSyst.cxx.

References caf::Proxy< caf::SRHeader >::cycle, caf::Proxy< caf::SRHeader >::evt, caf::Proxy< caf::StandardRecord >::hdr, makeTrainCVSamples::int, caf::Proxy< caf::SRHeader >::run, seed, caf::Proxy< caf::SRHeader >::subevt, caf::Proxy< caf::SRHeader >::subrun, and registry_explorer::v.

Referenced by ana::MichelTaggingSyst2018::MichelTaggingSyst2018(), MichelTaggingSyst2020(), Shift(), and ana::MichelTaggingSyst2018::Shift().

  {
    std::size_t seed = 0;
    std::vector<unsigned int> evtNums {sr->hdr.run, sr->hdr.subrun, (unsigned int)(sr->hdr.cycle), sr->hdr.evt, sr->hdr.subevt};
    for (auto v : evtNums)
      seed ^= std::hash<unsigned int>{}(v) + 0x9e3779b9 + (seed<<6) + (seed>>2);
    return seed;
  }

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::MichelTaggingSyst2020::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

sigma	Number of sigma to shift record by
sr	The record to inspect and alter
weight	Scale this weight for reweighting systematics

Reimplemented from ana::ISyst.

Definition at line 15 of file MichelTaggingSyst.cxx.

References std::abs(), om::cout, caf::Proxy< caf::SRSlcME >::deltat, caf::Proxy< caf::SRTrkME >::deltat, allTimeWatchdog::endl, fMRnd, fProbDO, fProbUP, Hash(), samweb_client.utility::hash, caf::Proxy< caf::StandardRecord >::me, caf::Proxy< caf::SRSlcME >::mid, caf::Proxy< caf::SRTrkME >::mid, caf::Proxy< caf::SRMichelE >::nkalman, caf::Proxy< caf::SRMichelE >::nslc, caf::Proxy< caf::SRParticleTruth >::pdg, PandAna.Demos.tute_pid_validation::slc, caf::Proxy< caf::SRMichelE >::slc, caf::Proxy< caf::SRMichelE >::trkkalman, caf::Proxy< caf::SRSlcME >::truth, and caf::Proxy< caf::SRTrkME >::truth.

Referenced by ana::MichelTaggingSyst2018::MichelTaggingSyst2018(), and MichelTaggingSyst2020().

  {
    if(sr->me.nslc != sr->me.slc.size()){
      std::cout << "me.nslc (" << sr->me.nslc << ") does not match "
                << "me.slc.size() (" << sr->me.slc.size() << ")..."
                << std::endl;
    }

    if(sr->me.nkalman != sr->me.trkkalman.size()){
      std::cout << "me.nkalman (" << sr->me.nkalman << ") does not match "
                << "me.trkkalman.size() (" << sr->me.trkkalman.size() << ")..."
                << std::endl;
    }

    // for the random numbers
    int hash = Hash(sr);
    fMRnd.SetSeed(hash);

    // are shifting the MEs UP or DOWN?
    double fProb = (sigma >= 0) ? fProbUP : fProbDO;

    // force to pass/fail MID and deltat criteria in me.slc
    for(caf::SRSlcMEProxy& slc: sr->me.slc){
      if(std::abs(slc.truth.pdg) == 11 &&
         fMRnd.Uniform() < std::abs(sigma*fProb)){
        slc.mid = (sigma >= 0) ? 4.0 : -1.0;
        if(slc.deltat <= 1200.) slc.deltat = 1201.;
      }
    }

    // force to pass/fail MID and deltat criteria in trkkalman
    for(caf::SRTrkMEProxy& trk: sr->me.trkkalman){
      if(std::abs(trk.truth.pdg) == 11 &&
         fMRnd.Uniform() < std::abs(sigma*fProb)){
        trk.mid = (sigma >= 0) ? 4.0 : -1.0;
        if(trk.deltat <= 1200.) trk.deltat = 1201.;
      }
    }
  }

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

    {
      // 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.
    }

Member Data Documentation

TRandom3 ana::MichelTaggingSyst2020::fMRnd

mutableprivate

Definition at line 25 of file MichelTaggingSyst.h.

Referenced by Shift(), and ana::MichelTaggingSyst2018::Shift().

double ana::MichelTaggingSyst2020::fProbDO = .05

private

Definition at line 23 of file MichelTaggingSyst.h.

Referenced by Shift(), and ana::MichelTaggingSyst2018::Shift().

double ana::MichelTaggingSyst2020::fProbUP = .08

private

Definition at line 22 of file MichelTaggingSyst.h.

Referenced by Shift(), and ana::MichelTaggingSyst2018::Shift().

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-12-03/3FlavorAna/Systs/MichelTaggingSyst.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-12-03/3FlavorAna/Systs/MichelTaggingSyst.cxx