fluxr::GSimpleInterface Class Reference

#include "GSimpleInterface.h"

Inheritance diagram for fluxr::GSimpleInterface:

Public Member Functions
	GSimpleInterface ()
	~GSimpleInterface ()
Long64_t	GetEntries () const
int	GetRun () const
float	GetPOT () const
TLorentzVector	GetNuPosition () const
TLorentzVector	GetNuMomentum () const
void	SetRootFile (TFile *rootFileName)
bool	FillMCFlux (Long64_t ientry, simb::MCFlux &mcflux)

Private Attributes
TTree *	fFluxTree
TTree *	fMetaTree
genie::flux::GSimpleNtpEntry *	fGSimpleEntry
genie::flux::GSimpleNtpNuMI *	fGSimpleNuMI
genie::flux::GSimpleNtpAux *	fGSimpleAux
genie::flux::GSimpleNtpMeta *	fGSimpleMeta
Long64_t	fNEntries
int	fRun
float	fPOT
TLorentzVector	fNuPos
TLorentzVector	fNuMom

Detailed Description

Definition at line 13 of file GSimpleInterface.h.

Constructor & Destructor Documentation

fluxr::GSimpleInterface::GSimpleInterface

(

)

Definition at line 6 of file GSimpleInterface.cxx.

{}

fluxr::GSimpleInterface::~GSimpleInterface

(

)

Definition at line 8 of file GSimpleInterface.cxx.

{}

Member Function Documentation

bool fluxr::GSimpleInterface::FillMCFlux ( Long64_t  ientry, simb::MCFlux &  mcflux  )

virtual

Implements fluxr::FluxInterface.

Definition at line 30 of file GSimpleInterface.cxx.

References simb::MCFlux::fdk2gen, simb::MCFlux::fevtno, fFluxTree, fGSimpleAux, fGSimpleEntry, fGSimpleMeta, fGSimpleNuMI, simb::MCFlux::fmupare, simb::MCFlux::fmuparpx, simb::MCFlux::fmuparpy, simb::MCFlux::fmuparpz, simb::MCFlux::fndecay, simb::MCFlux::fnecm, simb::MCFlux::fnenergyf, simb::MCFlux::fnenergyn, simb::MCFlux::fnimpwt, simb::MCFlux::fntype, fNuMom, fNuPos, simb::MCFlux::fnwtfar, simb::MCFlux::fnwtnear, simb::MCFlux::fpdpx, simb::MCFlux::fpdpy, simb::MCFlux::fpdpz, simb::MCFlux::fppdxdz, simb::MCFlux::fppdydz, simb::MCFlux::fppmedium, simb::MCFlux::fpppz, simb::MCFlux::fptype, simb::MCFlux::frun, simb::MCFlux::ftgen, simb::MCFlux::ftgptype, simb::MCFlux::ftptype, simb::MCFlux::ftpx, simb::MCFlux::ftpy, simb::MCFlux::ftpz, simb::MCFlux::fvx, simb::MCFlux::fvy, and simb::MCFlux::fvz.

Referenced by GetNuMomentum().

  {
    if (!fFluxTree->GetEntry(ientry)) return false;

    fNuPos = TLorentzVector(
      fGSimpleEntry->vtxx * 100, fGSimpleEntry->vtxy * 100, fGSimpleEntry->vtxz * 100, 0);
    fNuMom =
      TLorentzVector(fGSimpleEntry->px, fGSimpleEntry->py, fGSimpleEntry->pz, fGSimpleEntry->E);

    //stealing code from GenieHelper (nutools/EventGeneratorBase/GENIE/GENIEHelper.cxx)
    flux.fntype = fGSimpleEntry->pdg;
    flux.fnimpwt = fGSimpleEntry->wgt;
    flux.fdk2gen = fGSimpleEntry->dist;
    flux.fnenergyn = flux.fnenergyf = fGSimpleEntry->E;

    if (fGSimpleNuMI) {
      flux.frun = fGSimpleNuMI->run;
      flux.fevtno = fGSimpleNuMI->evtno;
      flux.ftpx = fGSimpleNuMI->tpx;
      flux.ftpy = fGSimpleNuMI->tpy;
      flux.ftpz = fGSimpleNuMI->tpz;
      flux.ftptype = fGSimpleNuMI->tptype; // converted to PDG
      flux.fvx = fGSimpleNuMI->vx;
      flux.fvy = fGSimpleNuMI->vy;
      flux.fvz = fGSimpleNuMI->vz;

      flux.fndecay = fGSimpleNuMI->ndecay;
      flux.fppmedium = fGSimpleNuMI->ppmedium;

      flux.fpdpx = fGSimpleNuMI->pdpx;
      flux.fpdpy = fGSimpleNuMI->pdpy;
      flux.fpdpz = fGSimpleNuMI->pdpz;

      double apppz = fGSimpleNuMI->pppz;
      if (TMath::Abs(fGSimpleNuMI->pppz) < 1.0e-30) apppz = 1.0e-30;
      flux.fppdxdz = fGSimpleNuMI->pppx / apppz;
      flux.fppdydz = fGSimpleNuMI->pppy / apppz;
      flux.fpppz = fGSimpleNuMI->pppz;

      flux.fptype = fGSimpleNuMI->ptype;
    }

    // anything useful stuffed into vdbl or vint?
    // need to check the metadata  auxintname, auxdblname
    if (fGSimpleAux && fGSimpleMeta) {
      // references just for reducing complexity
      const std::vector<std::string>& auxdblname = fGSimpleMeta->auxdblname;
      const std::vector<std::string>& auxintname = fGSimpleMeta->auxintname;
      const std::vector<int>& auxint = fGSimpleAux->auxint;
      const std::vector<double>& auxdbl = fGSimpleAux->auxdbl;

      for (size_t id = 0; id < auxdblname.size(); ++id) {
        if ("muparpx" == auxdblname[id]) flux.fmuparpx = auxdbl[id];
        if ("muparpy" == auxdblname[id]) flux.fmuparpy = auxdbl[id];
        if ("muparpz" == auxdblname[id]) flux.fmuparpz = auxdbl[id];
        if ("mupare" == auxdblname[id]) flux.fmupare = auxdbl[id];
        if ("necm" == auxdblname[id]) flux.fnecm = auxdbl[id];
        if ("nimpwt" == auxdblname[id]) flux.fnimpwt = auxdbl[id];
        if ("fgXYWgt" == auxdblname[id]) { flux.fnwtnear = flux.fnwtfar = auxdbl[id]; }
      }
      for (size_t ii = 0; ii < auxintname.size(); ++ii) {
        if ("tgen" == auxintname[ii]) flux.ftgen = auxint[ii];
        if ("tgptype" == auxintname[ii]) flux.ftgptype = auxint[ii];
      }
    }
    return true;
  }

Long64_t fluxr::GSimpleInterface::GetEntries ( ) const

inlinevirtual

Implements fluxr::FluxInterface.

Definition at line 18 of file GSimpleInterface.h.

References fNEntries.

{ return fNEntries; };

TLorentzVector fluxr::GSimpleInterface::GetNuMomentum ( ) const

inlinevirtual

Implements fluxr::FluxInterface.

Definition at line 22 of file GSimpleInterface.h.

References FillMCFlux(), fNuMom, and SetRootFile().

{ return fNuMom; };

TLorentzVector fluxr::GSimpleInterface::GetNuPosition ( ) const

inlinevirtual

Implements fluxr::FluxInterface.

Definition at line 21 of file GSimpleInterface.h.

References fNuPos.

{ return fNuPos; };

float fluxr::GSimpleInterface::GetPOT ( ) const

inlinevirtual

Implements fluxr::FluxInterface.

Definition at line 20 of file GSimpleInterface.h.

References fPOT.

{ return fPOT; };

int fluxr::GSimpleInterface::GetRun ( ) const

inlinevirtual

Implements fluxr::FluxInterface.

Definition at line 19 of file GSimpleInterface.h.

References fRun.

{ return fRun; };

void fluxr::GSimpleInterface::SetRootFile

(

TFile *

rootFileName

)

Definition at line 10 of file GSimpleInterface.cxx.

References fFluxTree, fGSimpleAux, fGSimpleEntry, fGSimpleMeta, fGSimpleNuMI, fMetaTree, fNEntries, fPOT, and fRun.

Referenced by GetNuMomentum().

  {
    fFluxTree = dynamic_cast<TTree*>(fluxInputFile->Get("flux"));
    fMetaTree = dynamic_cast<TTree*>(fluxInputFile->Get("meta"));
    fGSimpleEntry = new genie::flux::GSimpleNtpEntry;
    fGSimpleNuMI = new genie::flux::GSimpleNtpNuMI;
    fGSimpleAux = new genie::flux::GSimpleNtpAux;
    fGSimpleMeta = new genie::flux::GSimpleNtpMeta;
    fFluxTree->SetBranchAddress("entry", &fGSimpleEntry);
    fFluxTree->SetBranchAddress("numi", &fGSimpleNuMI);
    fFluxTree->SetBranchAddress("aux", &fGSimpleAux);
    fMetaTree->SetBranchAddress("meta", &fGSimpleMeta);

    fNEntries = fFluxTree->GetEntries();
    fFluxTree->GetEntry(0);
    fRun = fGSimpleNuMI->run;
    fMetaTree->GetEntry(0);
    fPOT = fGSimpleMeta->protons;
  }

Member Data Documentation

TTree* fluxr::GSimpleInterface::fFluxTree

private

Definition at line 28 of file GSimpleInterface.h.

Referenced by FillMCFlux(), and SetRootFile().

genie::flux::GSimpleNtpAux* fluxr::GSimpleInterface::fGSimpleAux

private

Definition at line 32 of file GSimpleInterface.h.

Referenced by FillMCFlux(), and SetRootFile().

genie::flux::GSimpleNtpEntry* fluxr::GSimpleInterface::fGSimpleEntry

private

Definition at line 30 of file GSimpleInterface.h.

Referenced by FillMCFlux(), and SetRootFile().

genie::flux::GSimpleNtpMeta* fluxr::GSimpleInterface::fGSimpleMeta

private

Definition at line 33 of file GSimpleInterface.h.

Referenced by FillMCFlux(), and SetRootFile().

genie::flux::GSimpleNtpNuMI* fluxr::GSimpleInterface::fGSimpleNuMI

private

Definition at line 31 of file GSimpleInterface.h.

Referenced by FillMCFlux(), and SetRootFile().

TTree* fluxr::GSimpleInterface::fMetaTree

private

Definition at line 29 of file GSimpleInterface.h.

Referenced by SetRootFile().

Long64_t fluxr::GSimpleInterface::fNEntries

private

Definition at line 34 of file GSimpleInterface.h.

Referenced by GetEntries(), and SetRootFile().

TLorentzVector fluxr::GSimpleInterface::fNuMom

private

Definition at line 38 of file GSimpleInterface.h.

Referenced by FillMCFlux(), and GetNuMomentum().

TLorentzVector fluxr::GSimpleInterface::fNuPos

private

Definition at line 37 of file GSimpleInterface.h.

Referenced by FillMCFlux(), and GetNuPosition().

float fluxr::GSimpleInterface::fPOT

private

Definition at line 36 of file GSimpleInterface.h.

Referenced by GetPOT(), and SetRootFile().

int fluxr::GSimpleInterface::fRun

private

Definition at line 35 of file GSimpleInterface.h.

Referenced by GetRun(), and SetRootFile().

---

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

• larsim/v09_43_00/source/larsim/PPFXFluxReader/GSimpleInterface.h
• larsim/v09_43_00/source/larsim/PPFXFluxReader/GSimpleInterface.cxx