ana::GetWtMinervaND Class Reference

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/CAFAna/Vars/Vars.h"

Public Member Functions
	GetWtMinervaND ()
double	operator() (const caf::SRProxy *sr) const

Protected Member Functions
void	InitializeMinervaND () const

Protected Attributes
TH1D *	fRatiosm [4] = {0}

Detailed Description

Definition at line 110 of file Vars.h.

Constructor & Destructor Documentation

ana::GetWtMinervaND::GetWtMinervaND ( )

inline

Definition at line 113 of file Vars.h.

References InitializeMinervaND(), operator()(), and sr.

{ }

Member Function Documentation

void ana::GetWtMinervaND::InitializeMinervaND ( ) const

protected

Definition at line 186 of file Vars.cxx.

References om::cerr, allTimeWatchdog::endl, fin, ana::FindCAFAnaDir(), plot_validation_datamc::fname, fRatiosm, and string.

Referenced by GetWtMinervaND(), and operator()().

                                                 {
    if(fRatiosm[0]) return;
    std::string fname = FindCAFAnaDir()+"/data/beam/MinervaFluxWtfromLeo.root";
    TFile fin(fname.c_str(), "read");
    if(fin.IsZombie()) {
      std::cerr << "Could not open " << fname << std::endl;
      return;
    }
    fRatiosm[0] = (TH1D*)fin.Get("NDnumu")->Clone();
    fRatiosm[0]->SetDirectory(0);
    fRatiosm[1] = (TH1D*)fin.Get("NDanumu")->Clone();
    fRatiosm[1]->SetDirectory(0);
    fRatiosm[2] = (TH1D*)fin.Get("NDnue")->Clone();
    fRatiosm[2]->SetDirectory(0);
    fRatiosm[3] = (TH1D*)fin.Get("NDanue")->Clone();
    fRatiosm[3]->SetDirectory(0);
  }

double ana::GetWtMinervaND::operator()

(

const caf::SRProxy *

sr

)

const

Definition at line 158 of file Vars.cxx.

References caf::Proxy< caf::SRHeader >::det, fRatiosm, caf::Proxy< caf::StandardRecord >::hdr, InitializeMinervaND(), caf::Proxy< caf::StandardRecord >::mc, caf::Proxy< caf::SRTruthBranch >::nnu, caf::Proxy< caf::SRTruthBranch >::nu, and ana::weight.

Referenced by GetWtMinervaND(), ana::GetWtMippNA49NOvAFD::GetWtMippNA49NOvAFD(), and ana::GetWtMippNA49NOvAND::GetWtMippNA49NOvAND().

  {
    InitializeMinervaND();
    if(sr->hdr.det!=1 || sr->mc.nnu==0) return 1;

    double te = sr->mc.nu[0].E;
    float weight = 1;
    if(sr->mc.nu[0].pdgorig == 14) {  // weights taken from Leo ratio plots
      if(te < 10){weight = fRatiosm[0]->Interpolate(te);}
      else {weight = 1;}
    }
    else if(sr->mc.nu[0].pdgorig == -14) {
      if(te < 10) {weight = fRatiosm[1]->Interpolate(te);}
      else {weight = 1;}
    }
    else if(sr->mc.nu[0].pdgorig == 12) {
      if(te < 10) {weight = fRatiosm[2]->Interpolate(te);}
      else {weight = 1;}
    }
    else if(sr->mc.nu[0].pdgorig == -12) {
      if(te < 10) {weight = fRatiosm[3]->Interpolate(te);}
      else {weight = 1;}
    }
    else weight = 1;
    //std::cout<<"energy: "<<te<<" weight "<<weight<<std::endl;
    return weight;
  }

Member Data Documentation

TH1D* ana::GetWtMinervaND::fRatiosm[4] = {0}

mutableprotected

Definition at line 117 of file Vars.h.

Referenced by InitializeMinervaND(), and operator()().

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/CAFAna/Vars/Vars.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/CAFAna/Vars/Vars.cxx