ana::BeamSystOrWeightBase Class Reference

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

Inheritance diagram for ana::BeamSystOrWeightBase:

Public Member Functions
	BeamSystOrWeightBase (const std::string &fname, const std::string &shortName)
	~BeamSystOrWeightBase ()

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

Definition at line 14 of file BeamSysts.h.

Member Enumeration Documentation

anonymous enum

protected

Enumerator
kNue
kAntiNue
kNumu
kAntiNumu
kNumFlavors

Definition at line 25 of file BeamSysts.h.

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

anonymous enum

protected

Enumerator
kFHC
kRHC
kNumFluxType

Definition at line 22 of file BeamSysts.h.

{kFHC, kRHC, kNumFluxType};

anonymous enum

protected

Enumerator
kND
kFD
kNumDets

Definition at line 23 of file BeamSysts.h.

{kND, kFD, kNumDets};

enum ana::BeamSystOrWeightBase::ESign

protected

Enumerator
kMinus
kPlus
kNumSigns

Definition at line 24 of file BeamSysts.h.

{kMinus, kPlus, kNumSigns};

Constructor & Destructor Documentation

ana::BeamSystOrWeightBase::BeamSystOrWeightBase	(	const std::string &	fname,
		const std::string &	shortName
	)

Definition at line 259 of file BeamSysts.cxx.

    : fFileName(fname), fHistName(histName), fHistos()
  {
  }

ana::BeamSystOrWeightBase::~BeamSystOrWeightBase

(

)

Definition at line 266 of file BeamSysts.cxx.

References fHistos, MECModelEnuComparisons::i, calib::j, kNumDets, kNumFlavors, kNumFluxType, kNumSigns, and submit_hadd::l.

  {
    for(int l = 0; l < kNumFluxType; ++l){
      for(int i = 0; i < kNumDets; ++i){
    for(int j = 0; j < kNumSigns; ++j){
      for(int k = 0; k < kNumFlavors; ++k){
        delete fHistos[l][i][j][k];
        fHistos[l][i][j][k] = 0;
      }
    }
      }
    }
  }

Member Function Documentation

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

protected

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, fHistName, fHistos, flux, make_syst_table_plots::h, InitializeHistograms(), caf::Proxy< caf::SRNeutrino >::isRHC, kAntiNue, kAntiNumu, caf::kFARDET, kFD, kFHC, kND, caf::kNEARDET, kNue, kNumu, kRHC, caf::Proxy< caf::SRNeutrino >::pdgorig, and canMan::sign().

Referenced by ana::BeamWeightFunc::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

protected

Definition at line 281 of file BeamSysts.cxx.

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

Referenced by 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();
  }

Member Data Documentation

std::string ana::BeamSystOrWeightBase::fFileName

protected

Definition at line 31 of file BeamSysts.h.

Referenced by InitializeHistograms().

std::string ana::BeamSystOrWeightBase::fHistName

protected

Definition at line 32 of file BeamSysts.h.

Referenced by GetWeight(), and InitializeHistograms().

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

mutableprotected

Definition at line 34 of file BeamSysts.h.

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

---

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

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