ana::NuMIFluxSyst Class Reference

#include <NuMIFluxSysts.h>

Inheritance diagram for ana::NuMIFluxSyst:

Public Member Functions
virtual	~NuMIFluxSyst ()
void	Shift (double sigma, caf::SRSliceProxy *slc, double &weight) const override
	Perform the systematic shift. More...
Public Member Functions inherited from ana::ISyst
	ISyst (const std::string &shortName, const std::string &latexName, bool applyPenalty=true, double min=-3, double max=+3)
	ISyst (const ISyst &)=delete
	ISyst (ISyst &&rhs)=delete
virtual	~ISyst ()
void	operator= (const ISyst &)=delete
void	operator= (ISyst &&)=delete
virtual std::string	ShortName () const final
	The name printed out to the screen. More...
virtual std::string	LatexName () const final
	The name used on plots (ROOT's TLatex syntax) More...
virtual double	Penalty (double x) const
virtual bool	ApplyPenalty () const
	Should a penalty be applied for this shift? More...
virtual double	Min () const
	Return the min/max value for this syst. More...
virtual double	Max () const
virtual int	PredInterpMaxNSigma () const

Protected Member Functions
	NuMIFluxSyst (const std::string &dir, const std::string &prefix, const std::string &name)

Protected Attributes
std::string	fHistName
std::string	fName
TH1 *	fScale [2][2][2]

Friends
const NuMIFluxSyst *	GetNuMIFluxSyst (const std::string &, const std::string &, const std::string &)

Detailed Description

Definition at line 12 of file NuMIFluxSysts.h.

Constructor & Destructor Documentation

ana::NuMIFluxSyst::~NuMIFluxSyst ( )

virtual

Definition at line 18 of file NuMIFluxSysts.cxx.

  {
    for(int i = 0; i < 2; ++i)
      for(int j = 0; j < 2; ++j)
        for(int k = 0; k < 2; ++k)
          delete fScale[i][j][k];
  }

ana::NuMIFluxSyst::NuMIFluxSyst ( const std::string &  dir, const std::string &  prefix, const std::string &  name  )

inlineprotected

Definition at line 24 of file NuMIFluxSysts.h.

      : ISyst("numi_"+name, "NuMI flux: "+name),
        fHistName(dir+"/"+prefix+name), fName(name), fScale()
    {
    }

Member Function Documentation

void ana::NuMIFluxSyst::Shift ( double  sigma, caf::SRSliceProxy *  sr, double &  weight  ) const

overridevirtual

Perform the systematic shift.

Parameters

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

Implements ana::ISyst.

Definition at line 27 of file NuMIFluxSysts.cxx.

  {
    if(slc->truth.index < 0) return; // No neutrino

    // TODO switch to regular detector flag as soon as it's filled reliably
    if(slc->truth.baseline < 500){
      std::cout << "Using NuMIFluxSyst on what appears not to be an icarus file (baseline = " << slc->truth.baseline << ")!" << std::endl;
      abort();
    }

    if(!fScale[0][0][0]){
      const char* sbndata = getenv("SBNDATA_DIR");
      if(!sbndata){
        std::cout << "NuMIFluxSysts: $SBNDATA_DIR environment variable not set. Please setup the sbndata product." << std::endl;
        abort();
      }

      const std::string fname = std::string(sbndata)+"/beamData/NuMIdata/icarus_numi_flux_syst_ana.root";
      TFile f(fname.c_str());
      if(f.IsZombie()){
        std::cout << "NuMIFluxSysts: Failed to open " << fname << std::endl;
        abort();
      }

      for(int hcIdx: {0, 1}){
        for(int flavIdx: {0, 1}){
          for(int signIdx: {0, 1}){
            std::string hName = fHistName;
            if(hcIdx == 0) hName += "_fhc"; else hName += "_rhc";
            if(flavIdx == 0) hName += "_nue"; else hName += "_numu";
            if(signIdx == 1) hName += "bar";

            TH1* h = (TH1*)f.Get(hName.c_str());
            if(!h){
              std::cout << "NuMIFluxSysts: failed to find " << hName << " in " << f.GetName() << std::endl;
              abort();
            }
            h = (TH1*)h->Clone(UniqueName().c_str());
            h->SetDirectory(0);

            fScale[hcIdx][flavIdx][signIdx] = h;
          }
        }
      }
    } // end if

    const int flav = abs(slc->truth.initpdg);
    if(flav != 12 && flav != 14) return;

    const int hcIdx = 0; // TODO TODO always FHC
    const int flavIdx = (flav == 12) ? 0 : 1;
    const int signIdx = (slc->truth.initpdg > 0) ? 0 : 1;

    TH1* h = fScale[hcIdx][flavIdx][signIdx];
    assert(h);

    const int bin = h->FindBin(slc->truth.E);

    if(bin == 0 || bin == h->GetNbinsX()+1) return;

    weight *= 1 + sigma*h->GetBinContent(bin);
  }

Friends And Related Function Documentation

const NuMIFluxSyst* GetNuMIFluxSyst ( const std::string &  , const std::string &  , const std::string &    )

friend

Definition at line 91 of file NuMIFluxSysts.cxx.

  {
    // Make sure we always give the same one back
    static std::map<std::string, const NuMIFluxSyst*> cache;

    const std::string key = dir+"/"+prefix+name;
    if(cache.count(key) == 0) cache[key] = new NuMIFluxSyst(dir, prefix, name);

    return cache[key];
  }

Member Data Documentation

std::string ana::NuMIFluxSyst::fHistName

protected

Definition at line 32 of file NuMIFluxSysts.h.

std::string ana::NuMIFluxSyst::fName

protected

Definition at line 32 of file NuMIFluxSysts.h.

TH1* ana::NuMIFluxSyst::fScale[2][2][2]

mutableprotected

Definition at line 34 of file NuMIFluxSysts.h.

---

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

• NuMIFluxSysts.h
• NuMIFluxSysts.cxx