ana::BoosterFluxHadronSyst Class Reference

#include <BoosterFluxSysts.h>

Inheritance diagram for ana::BoosterFluxHadronSyst:

Public Member Functions
virtual	~BoosterFluxHadronSyst ()
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
	BoosterFluxHadronSyst (int i)

Protected Attributes
int	fIdx
TH1 *	fScale [3]

Friends
const BoosterFluxHadronSyst *	GetBoosterFluxHadronSyst (unsigned int)

Detailed Description

Definition at line 11 of file BoosterFluxSysts.h.

Constructor & Destructor Documentation

ana::BoosterFluxHadronSyst::~BoosterFluxHadronSyst ( )

virtual

Definition at line 15 of file BoosterFluxSysts.cxx.

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

ana::BoosterFluxHadronSyst::BoosterFluxHadronSyst ( int  i )

inlineprotected

Definition at line 21 of file BoosterFluxSysts.h.

      : ISyst(TString::Format("booster_%i", i).Data(),
              TString::Format("Booster flux hadron syst #%i", i).Data()),
        fIdx(i), fScale() {}

Member Function Documentation

void ana::BoosterFluxHadronSyst::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 21 of file BoosterFluxSysts.cxx.

{
  // Only implemented for numus so far
  if(/*sr->mc.nnu == 0 ||*/ slc->truth.initpdg != 14) return;

  if(!fScale[0]){
    TFile f((FindCAFAnaDir() + "/Systs/booster_flux_shifts.root").c_str());
    assert(!f.IsZombie());

    for(int det = 0; det < 3; ++det){
      std::string detStr;
      if(det == 0) detStr = "nd";
      if(det == 1) detStr = "ub";
      if(det == 2) detStr = "fd";

      TH1*& h = fScale[det];

      h = (TH1*)f.Get(TString::Format("syst%d/%s_numu", fIdx, detStr.c_str()).Data());
      h = (TH1*)h->Clone(UniqueName().c_str());
      assert(h);
      h->SetDirectory(0);
    }
  } // end if

  // TODO use the regular detector flag as soon as it's filled reliably
  int det;
  const double L = slc->truth.baseline;
  if(L < 150) det = 0;
  else if(L < 500) det = 1;
  else det = 2;

  TH1* h = fScale[det];
  assert(h);

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

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

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

Friends And Related Function Documentation

const BoosterFluxHadronSyst* GetBoosterFluxHadronSyst ( unsigned  int )

friend

Definition at line 63 of file BoosterFluxSysts.cxx.

{
  // Make sure we always give the same one back
  static std::vector<const BoosterFluxHadronSyst*> cache;

  if(i >= cache.size()) cache.resize(i+1);

  if(!cache[i]) cache[i] = new BoosterFluxHadronSyst(i);

  return cache[i];
}

Member Data Documentation

int ana::BoosterFluxHadronSyst::fIdx

protected

Definition at line 26 of file BoosterFluxSysts.h.

TH1* ana::BoosterFluxHadronSyst::fScale[3]

mutableprotected

Definition at line 28 of file BoosterFluxSysts.h.

---

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

• BoosterFluxSysts.h
• BoosterFluxSysts.cxx