ana::NueSyst Class Reference

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

Inheritance diagram for ana::NueSyst:

Public Member Functions
	NueSyst (std::string name, std::string label, std::string sample_name="", std::map< int, std::pair< TH1 *, TH1 * > > shifts={})
	Implementation of standard ISyst::Shift function for NueSyst. More...
virtual void	Shift (double sigma, caf::SRProxy *sr, double &weight) const override
	Implementation of standard ISyst::Shift function for NueSyst. More...
virtual void	SaveTo (TDirectory *dir) const override
	SaveTo implementation for NueSyst. More...
std::string	SampleName () const
std::string	BaseName () const
	Return base name of systematic. 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...

Static Public Member Functions
static std::unique_ptr< NueSyst >	LoadFrom (TDirectory *dir)
	LoadFrom implementation for NueSyst. More...

Protected Member Functions
virtual double	WeightFor (SystChannel channel, double sigma, double calE) const
	Get weight for a given energy, channel and sigma shift. More...

Protected Attributes
std::string	fSampleName
std::map< int, std::pair< TH1 *, TH1 * > >	fShifts

Detailed Description

Definition at line 18 of file CovMxSysts.h.

Constructor & Destructor Documentation

ana::NueSyst::NueSyst	(	std::string	name,
		std::string	label,
		std::string	sample_name = "",
		std::map< int, std::pair< TH1 *, TH1 * > >	shifts = {}
	)

Implementation of standard ISyst::Shift function for NueSyst.

Definition at line 117 of file CovMxSysts.cxx.

    : NuISyst(name, label, sampleName, shifts)
  {}

Member Function Documentation

string ana::NuISyst::BaseName ( ) const

inherited

Return base name of systematic.

Definition at line 44 of file CovMxSysts.cxx.

References ana::NuISyst::fSampleName, and ana::ISyst::ShortName().

                                 {
    if (fSampleName.empty()) return ShortName(); // If there's no sample prefix to remove
    return ShortName().substr(fSampleName.size()+1);
  }

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;}

unique_ptr< NueSyst > ana::NueSyst::LoadFrom ( TDirectory *  dir )

static

LoadFrom implementation for NueSyst.

Definition at line 171 of file CovMxSysts.cxx.

References ana::assert(), findDuplicateFiles::key, label, make_pair(), next(), and getGoodRuns4SAM::tag.

                                                       {

    DontAddDirectory guard;

    // Check this is a NueSyst
    TObjString* tag = (TObjString*)dir->Get("type");
    assert(tag);
    string sTag = tag->GetString().Data();
    assert(sTag == "NueSyst");

    // Load names
    TObjString* name = (TObjString*)dir->Get("name");
    assert(name);
    string sName = name->GetString().Data();
    TObjString* label = (TObjString*)dir->Get("label");
    assert(label);
    string sLabel = label->GetString().Data();
    TObjString* sampleName = (TObjString*)dir->Get("samplename");
    assert(sampleName);
    string sSampleName = sampleName->GetString().Data();

    map<int, pair<TH1*, TH1*>> shifts;

    // Load all shifts
    TIter next(dir->GetDirectory("Sigmas")->GetListOfKeys());
    TKey* key;
    while ((key = (TKey*)next())) {
      TDirectory* sigmaDir = (TDirectory*)key->ReadObj();
      const TVectorD vSigma = *(TVectorD*)sigmaDir->Get("sigma");
      TH1* hSignal = (TH1*)sigmaDir->Get("sig");
      TH1* hBackground = (TH1*)sigmaDir->Get("bkg");
      shifts[vSigma[0]] = make_pair(hSignal, hBackground);
    }

    return make_unique<NueSyst>(sName, sLabel, sSampleName, shifts);

  } // function NumuSyst::SaveTo

std::string ana::NuISyst::SampleName ( ) const

inlineinherited

Definition at line 25 of file CovMxSysts.h.

References dir, ana::SaveTo(), sigma(), and string.

{ return fSampleName; };

void ana::NueSyst::SaveTo ( TDirectory *  dir ) const

overridevirtual

SaveTo implementation for NueSyst.

Reimplemented from ana::NuISyst.

Definition at line 159 of file CovMxSysts.cxx.

References ana::NuISyst::SaveTo(), and tmp.

                                            {

    TDirectory* tmp = gDirectory;
    dir->cd();
    TObjString("NueSyst").Write("type");
    tmp->cd();
    NuISyst::SaveTo(dir);

  } // function NueSyst::SaveTo

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

overridevirtual

Implementation of standard ISyst::Shift function for NueSyst.

Reimplemented from ana::ISyst.

Definition at line 124 of file CovMxSysts.cxx.

References febshutoff_auto::chan, om::cout, caf::Proxy< caf::SRHeader >::det, allTimeWatchdog::endl, caf::Proxy< caf::StandardRecord >::hdr, caf::Proxy< caf::SRSpill >::isFHC, caf::Proxy< caf::SRSpill >::isRHC, ana::kBkg, caf::kFARDET, caf::kNEARDET, ana::kNue2020AnaBin, ana::kNue2020FDAllSamples_ML, ana::kNue2020ND, ana::kSig, ana::kUnknown, caf::Proxy< caf::StandardRecord >::mc, caf::Proxy< caf::SRTruthBranch >::nnu, caf::Proxy< caf::SRTruthBranch >::nu, canMan::sign(), caf::Proxy< caf::StandardRecord >::spill, and ana::NuISyst::WeightFor().

                                                                        {

    if (sr->hdr.det == caf::kNEARDET && !kNue2020ND(sr)) return;
    if (sr->hdr.det == caf::kFARDET && !kNue2020FDAllSamples_ML(sr)) return;

    // Find the energy estimate for this slice
    double calE = kNue2020AnaBin(sr);

    // Find the signal or background oscillation channel
    SystChannel chan = SystChannel::kUnknown;
    int sign;
    if (sr->spill.isFHC) sign = 1;
    else if (sr->spill.isRHC) sign = -1;
    else throw runtime_error("Event is neither FHC nor RHC!");

    if (sr->mc.nnu == 0) {
      chan = SystChannel::kBkg;
    }
    else if (sr->mc.nu[0].iscc && sr->mc.nu[0].pdg == sign*12 && sr->mc.nu[0].pdgorig == sign*14) {
      chan = SystChannel::kSig;
    }
    else if (!sr->mc.nu[0].iscc || sr->mc.nu[0].pdg != sign*12 || sr->mc.nu[0].pdgorig != sign*14) {
      chan = SystChannel::kBkg;
    }
    if (chan == SystChannel::kUnknown) {
      cout << "Error: Unknown oscillation channel." << endl;
      abort();
    }

    weight *= WeightFor(chan, sigma, calE);

  } // function NueSyst::Shift

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

double ana::NuISyst::WeightFor ( SystChannel  channel, double  sigma, double  calE  ) const

protectedvirtualinherited

Get weight for a given energy, channel and sigma shift.

Definition at line 80 of file CovMxSysts.cxx.

References ana::assert(), bin, febshutoff_auto::end, ana::NuISyst::fShifts, it, ana::kSig, std::max(), and runNovaSAM::ret.

Referenced by ana::NumuSyst::Shift(), ana::NCSyst::Shift(), and Shift().

                                                                                {

    // If it's a one-sided shift (ie. no -1 sigma weights) then we want to return
    // 1 for a negative sigma
    if (sigma < 0 and not fShifts.count(-1)) return 1;
    const int bin = fShifts.at(0).first->FindBin(calE);

    auto low = fShifts.begin();
    auto end = fShifts.end();
    std::advance(end, -2);
    if (sigma < fShifts.begin()->first)
      low = fShifts.begin();
    else if (sigma >= end->first)
      low = end;
    else {
      for (auto it = fShifts.begin(); it != end; ++it) {
        if (sigma >= it->first) {
          low = it;
          break;
        }
      }
    }

    auto high = low;
    ++high;

    // // Why would we have templates differing by more than 1 sigma?
    // // fracpart below assumes this
    assert(high->first - low->first == 1);

    TH1* hLow = nullptr;
    TH1* hHigh = nullptr;

    if (channel == SystChannel::kSig) {
      hLow = low->second.first;
      hHigh = high->second.first;
    }
    else {
      hLow = low->second.second;
      hHigh = high->second.second;
    }

    const double fracpart = sigma - low->first;
    const double ret = (fracpart*hHigh->GetBinContent(bin)) +
      ((1-fracpart)*hLow->GetBinContent(bin));

    return std::max(0., ret); // Keep the LL from blowing up

  } // function NuISyst::WeightFor

Member Data Documentation

std::string ana::NuISyst::fSampleName

protectedinherited

Definition at line 31 of file CovMxSysts.h.

Referenced by ana::NuISyst::BaseName(), and ana::NuISyst::SaveTo().

std::map<int, std::pair<TH1*, TH1*> > ana::NuISyst::fShifts

protectedinherited

Definition at line 32 of file CovMxSysts.h.

Referenced by ana::NuISyst::SaveTo(), ana::NuISyst::WeightFor(), and ana::NuISyst::~NuISyst().

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-11-30/3FlavorAna/Systs/CovMxSysts.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-11-30/3FlavorAna/Systs/CovMxSysts.cxx