ana::BeamWeightFunc Class Reference

Incorporate into a weight Var like Var wei(BeamWeight(... More...

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-03-08/CAFAna/Systs/BeamSysts.h"

Inheritance diagram for ana::BeamWeightFunc:

Public Member Functions
	BeamWeightFunc (const std::string &fname, const std::string &histname)
double	operator() (const caf::SRNeutrinoProxy *sr) const

Protected Types
enum	{ kFHC, kRHC, kNumFluxType }
enum	{ kND, kFD, kNumDets }
enum	ESign { kMinus, kPlus, kNumSigns }
enum	{   kNue, kAntiNue, kNumu, kAntiNumu,   kNumFlavors }

Protected Member Functions
void	InitializeHistograms () const
double	GetWeight (const caf::SRNeutrinoProxy *nu, ESign sign) const

Protected Attributes
std::string	fFileName
std::string	fHistName
TH1D *	fHistos [kNumFluxType][kNumDets][kNumSigns][kNumFlavors]

Detailed Description

Incorporate into a weight Var like Var wei(BeamWeight(...

Definition at line 55 of file BeamSysts.h.

Member Enumeration Documentation

anonymous enum

protectedinherited

Enumerator
kNue
kAntiNue
kNumu
kAntiNumu
kNumFlavors

Definition at line 25 of file BeamSysts.h.

{kNue, kAntiNue, kNumu, kAntiNumu, kNumFlavors};

anonymous enum

protectedinherited

Enumerator
kFHC
kRHC
kNumFluxType

Definition at line 22 of file BeamSysts.h.

{kFHC, kRHC, kNumFluxType};

anonymous enum

protectedinherited

Enumerator
kND
kFD
kNumDets

Definition at line 23 of file BeamSysts.h.

{kND, kFD, kNumDets};

enum ana::BeamSystOrWeightBase::ESign

protectedinherited

Enumerator
kMinus
kPlus
kNumSigns

Definition at line 24 of file BeamSysts.h.

{kMinus, kPlus, kNumSigns};

Constructor & Destructor Documentation

ana::BeamWeightFunc::BeamWeightFunc ( const std::string &  fname, const std::string &  histname  )

inline

Definition at line 58 of file BeamSysts.h.

References ana::GetPPFXFluxUnivWgtSmooth_NT(), ana::GetPPFXUnivSyst(), and sr.

Referenced by ana::GetPPFXFluxUnivWgtSmooth_NT().

      : BeamSystOrWeightBase(fname, histname)
    {
    }

Member Function Documentation

double ana::BeamSystOrWeightBase::GetWeight ( const caf::SRNeutrinoProxy *  nu, ESign  sign  ) const

protectedinherited

Definition at line 350 of file BeamSysts.cxx.

References om::cerr, fillBadChanDBTables::det, caf::Proxy< caf::SRNeutrino >::det, caf::Proxy< caf::SRNeutrino >::E, allTimeWatchdog::endl, energy, ana::BeamSystOrWeightBase::fHistName, ana::BeamSystOrWeightBase::fHistos, flux, make_syst_table_plots::h, ana::BeamSystOrWeightBase::InitializeHistograms(), caf::Proxy< caf::SRNeutrino >::isRHC, ana::BeamSystOrWeightBase::kAntiNue, ana::BeamSystOrWeightBase::kAntiNumu, caf::kFARDET, ana::BeamSystOrWeightBase::kFD, ana::BeamSystOrWeightBase::kFHC, ana::BeamSystOrWeightBase::kND, caf::kNEARDET, ana::BeamSystOrWeightBase::kNue, ana::BeamSystOrWeightBase::kNumu, ana::BeamSystOrWeightBase::kRHC, caf::Proxy< caf::SRNeutrino >::pdgorig, and canMan::sign().

Referenced by operator()(), and ana::BeamSyst::TruthShift().

  {
    InitializeHistograms();

    int flux = kFHC;
    if(nu->isRHC) flux = kRHC;

    // Need different histogram per detector
    int det = 0;
    switch (nu->det)
    {
      case caf::kNEARDET : det = kND; break;
      case caf::kFARDET : det = kFD; break;
      default: std::cerr <<"Wrong detector; ignore" << std::endl; return 1;
    }

    //Set the different weights for p/m sigma and neutrino flavor/energy

    double energy = nu->E; // True neutrino energy

    //Find the neutrino flavor
    int flav = 0;
    switch(nu->pdgorig)
    {
      case  14: flav = kNumu;     break;
      case  12: flav = kNue;      break;
      case -14: flav = kAntiNumu; break;
      case -12: flav = kAntiNue;  break;
      default: std::cerr << "Wrong nu flavor; ignore" << std::endl; return 1;
    }

    //Find the right histogram for the detector and sign
    TH1D* h = fHistos[flux][det][sign][flav];
    if (h == 0){
      std::cerr << fHistName+": Can't find desired histogram; ignore"<< std:: endl;
      return 1;
    }
    //Only use weights if the energy is in the range of the histogram
    if (energy > h->GetXaxis()->GetXmin() &&
        energy < h->GetXaxis()->GetXmax() ){
      return h->Interpolate(energy);
    }

    //Main use case is for PPFX smoothened weights. 
    //Return the weight in the last bin if the energy is larger than the range of the histogram 
    //Generalized for all BeamWeightFuncs and BeamSysts
    if (energy > h->GetXaxis()->GetXmax())
      return h->GetBinContent(h->GetNbinsX());
    
    return 1;
  }

void ana::BeamSystOrWeightBase::InitializeHistograms ( ) const

protectedinherited

Definition at line 281 of file BeamSysts.cxx.

References om::cerr, fillBadChanDBTables::det, allTimeWatchdog::endl, ana::BeamSystOrWeightBase::fFileName, ana::BeamSystOrWeightBase::fHistName, ana::BeamSystOrWeightBase::fHistos, fin, flux, compareCafs::histName, ana::BeamSystOrWeightBase::kAntiNue, ana::BeamSystOrWeightBase::kAntiNumu, ana::BeamSystOrWeightBase::kFD, ana::BeamSystOrWeightBase::kFHC, ana::BeamSystOrWeightBase::kMinus, ana::BeamSystOrWeightBase::kND, ana::BeamSystOrWeightBase::kNue, ana::BeamSystOrWeightBase::kNumu, ana::BeamSystOrWeightBase::kPlus, ana::BeamSystOrWeightBase::kRHC, and canMan::sign().

Referenced by ana::BeamSystOrWeightBase::GetWeight().

  {
    // Already been called
    if(fHistos[0][0][0][0]) return;

    //open file
    TFile fin (fFileName.c_str(),"read");
    if(fin.IsZombie()){
      std::cerr << "Warning: couldn't open " << fFileName << std::endl;
      return;
    }

    //Check that desired systematic is in file
    // Loop through the detector folder, over all histograms

    TIter iterHist(gDirectory->GetListOfKeys());
    TKey* keyHist;
    int foundHisto=0;
    bool isSeparatedByFlavor = false;

    while((keyHist = (TKey*)iterHist())) {
      TString histName = keyHist->GetName(); // Get a histogram name
      if(histName.Contains(fHistName)){   //is it the desired syst
    foundHisto++;

    int flux=0;
    if(histName.Contains("FHC")) flux = kFHC;
    if(histName.Contains("RHC")) flux = kRHC;

        int det = 0;
        if(histName.Contains("ND")) det = kND;
        if(histName.Contains("FD")) det = kFD;

        int sign = 0;
        if(histName.Contains("minus")) sign = kMinus;
        if(histName.Contains("plus" )) sign = kPlus;

        int flav = 0;
        if(histName.Contains("numu"))  flav = kNumu;
        if(histName.Contains("nue"))   flav = kNue;
        if(histName.Contains("anumu")) flav = kAntiNumu;
        if(histName.Contains("anue"))  flav = kAntiNue;

        //store relevant histograms
    fHistos[flux][det][sign][flav] = (TH1D*) fin.Get(histName)->Clone();
        // disassociate it from the file it came from
    fHistos[flux][det][sign][flav]->SetDirectory(0);
      }
      if(histName.Contains("numu")){
    isSeparatedByFlavor=true;
      }
    }

    if (foundHisto==0) {
      std::cerr << "Beam systematic histogram " << fHistName
        << " is not in file " << fFileName <<"; aborting" << std::endl;
      abort();
    }

    if (!isSeparatedByFlavor) {
      std::cerr << "Beam systematic "<< fHistName
        << " doesn't have required flavor information; aborting"
        << std::endl;
      abort();
    }
    fin.Close();
  }

double ana::BeamWeightFunc::operator()

(

const caf::SRNeutrinoProxy *

sr

)

const

Definition at line 23 of file BeamSysts.cxx.

References ana::BeamSystOrWeightBase::GetWeight(), and ana::BeamSystOrWeightBase::kPlus.

  {
    // if (sr->mc.nnu != 1)
    //   return 1.;

    return GetWeight(sr, kPlus);
  }

Member Data Documentation

std::string ana::BeamSystOrWeightBase::fFileName

protectedinherited

Definition at line 31 of file BeamSysts.h.

Referenced by ana::BeamSystOrWeightBase::InitializeHistograms().

std::string ana::BeamSystOrWeightBase::fHistName

protectedinherited

Definition at line 32 of file BeamSysts.h.

Referenced by ana::BeamSystOrWeightBase::GetWeight(), and ana::BeamSystOrWeightBase::InitializeHistograms().

TH1D* ana::BeamSystOrWeightBase::fHistos[kNumFluxType][kNumDets][kNumSigns][kNumFlavors]

mutableprotectedinherited

Definition at line 34 of file BeamSysts.h.

Referenced by ana::BeamSystOrWeightBase::GetWeight(), ana::BeamSystOrWeightBase::InitializeHistograms(), and ana::BeamSystOrWeightBase::~BeamSystOrWeightBase().

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-03-08/CAFAna/Systs/BeamSysts.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-03-08/CAFAna/Systs/BeamSysts.cxx