evgen::ldm::IMeVPrtlFlux Class Reference

IMeVPrtlFlux interface class definiton. More...

#include <IMeVPrtlFlux.h>

Inheritance diagram for evgen::ldm::IMeVPrtlFlux:

Public Member Functions
virtual	~IMeVPrtlFlux () noexcept=default
	Virtual Destructor. More...
virtual bool	MakeFlux (const simb::MCFlux &mcflux, MeVPrtlFlux &flux, double &weight)=0
	IMeVPrtlFlux (const fhicl::ParameterSet &pset)
Public Member Functions inherited from evgen::ldm::IMeVPrtlStage
virtual	~IMeVPrtlStage () noexcept
	Virtual Destructor. More...
	IMeVPrtlStage (const char *name)
virtual void	configure (const fhicl::ParameterSet &)=0
	Interface for configuring the particular algorithm tool. More...
virtual double	MaxWeight ()=0
TVector3	RandomUnitVector ()
double	GetRandom ()
const char *	Name ()

Protected Member Functions
TLorentzVector	BeamOrigin ()
double	EnuLab (double enucm, TLorentzVector meson_mom, TLorentzVector meson_pos)

Protected Attributes
evgb::EvtTimeShiftI *	fTimeShiftMethod
TRotation	fBeam2Det
TVector3	fBeamOrigin
std::string	fSpillTimeConfig
Protected Attributes inherited from evgen::ldm::IMeVPrtlStage
CLHEP::HepRandomEngine *	fEngine
const char *	fName

Detailed Description

IMeVPrtlFlux interface class definiton.

Definition at line 39 of file IMeVPrtlFlux.h.

Constructor & Destructor Documentation

virtual evgen::ldm::IMeVPrtlFlux::~IMeVPrtlFlux ( )

virtualdefaultnoexcept

Virtual Destructor.

evgen::ldm::IMeVPrtlFlux::IMeVPrtlFlux ( const fhicl::ParameterSet &  pset )

inline

Definition at line 49 of file IMeVPrtlFlux.h.

    {
      // rotation matrix
      std::vector<double> rotation = pset.get<std::vector<double>>("Beam2DetectorRotation");
      assert(rotation.size() == 9);
      // get the three axes
      TVector3 beam2DetX(rotation[0], rotation[3], rotation[6]);
      TVector3 beam2DetY(rotation[1], rotation[4], rotation[7]);
      TVector3 beam2DetZ(rotation[2], rotation[5], rotation[8]);

      fBeam2Det.RotateAxes(beam2DetX, beam2DetY, beam2DetZ);

      // beam origin
      std::vector<double> origin = pset.get<std::vector<double>>("BeamOrigin");
      assert(origin.size() == 3 || origin.size() == 6);
      // beam origin specified in detector rotation-frame
      if (origin.size() == 3) {
        fBeamOrigin.SetXYZ(origin[0], origin[1], origin[2]);
      } 
      // beam origin specified in beam rotation-frame
      else if (origin.size() == 6) {
        TVector3 userpos = TVector3(origin[0], origin[1], origin[2]);
        TVector3 beampos = TVector3(origin[3], origin[4], origin[5]);
        fBeamOrigin = userpos - fBeam2Det * beampos;
      }

      // get random seed for stuff
      unsigned seed = art::ServiceHandle<rndm::NuRandomService>()->getSeed();

      // use the time-shifting tools from GENIE
      fTimeShiftMethod = NULL;
      if (fSpillTimeConfig != "") {
        fTimeShiftMethod = evgb::EvtTimeShiftFactory::Instance().GetEvtTimeShift(fSpillTimeConfig);
        if ( fTimeShiftMethod ) {
          if ( ! fTimeShiftMethod->IsRandomGeneratorSeeded() ) {
            fTimeShiftMethod->GetRandomGenerator()->SetSeed(seed);
          }
          fTimeShiftMethod->PrintConfig();
        } 
        else {
          evgb::EvtTimeShiftFactory::Instance().Print();
        }
      }
      if (fTimeShiftMethod) {
        std::cout << "Timing Config:\n";
        fTimeShiftMethod->PrintConfig();
        std::cout << std::endl;
      }
      std::cout << "Neutrino TIF: " << (fBeamOrigin.Mag()/Constants::Instance().c_cm_per_ns) << std::endl;
    }

Member Function Documentation

TLorentzVector evgen::ldm::IMeVPrtlFlux::BeamOrigin ( )

inlineprotected

Definition at line 107 of file IMeVPrtlFlux.h.

                              {
    double toff = fTimeShiftMethod ? fTimeShiftMethod->TimeOffset() : 0.;

    // TODO: what to do here? For now -- don't shift time at all
    //
    // subtract out the delay of neutrinos reaching the beam
    // double neutrino_tif = fBeamOrigin.Mag()/Constants::Instance().c_cm_per_ns;
    // toff -= neutrino_tif;
    return TLorentzVector(fBeamOrigin, toff);
  }

double evgen::ldm::IMeVPrtlFlux::EnuLab ( double  enucm, TLorentzVector  meson_mom, TLorentzVector  meson_pos  )

inlineprotected

Definition at line 119 of file IMeVPrtlFlux.h.

                                                                                  { // all in detector coordinates
    // Assume neutrino travels to center of detector
    double costh = meson_mom.Vect().Unit().Dot(-meson_pos.Vect().Unit());

    // Scale factor
    double M = 1. / (meson_mom.Gamma() * (1 - meson_mom.Beta() * costh));

    return M * enucm;
  }

virtual bool evgen::ldm::IMeVPrtlFlux::MakeFlux ( const simb::MCFlux &  mcflux, MeVPrtlFlux &  flux, double &  weight  )

pure virtual

Implemented in evgen::ldm::Kaon2HiggsFlux, evgen::ldm::Kaon2HNLFlux, and evgen::ldm::MonoEnergyHiggsFlux.

Member Data Documentation

TRotation evgen::ldm::IMeVPrtlFlux::fBeam2Det

protected

Definition at line 103 of file IMeVPrtlFlux.h.

TVector3 evgen::ldm::IMeVPrtlFlux::fBeamOrigin

protected

Definition at line 104 of file IMeVPrtlFlux.h.

std::string evgen::ldm::IMeVPrtlFlux::fSpillTimeConfig

protected

Definition at line 105 of file IMeVPrtlFlux.h.

evgb::EvtTimeShiftI* evgen::ldm::IMeVPrtlFlux::fTimeShiftMethod

protected

Definition at line 102 of file IMeVPrtlFlux.h.

---

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

• IMeVPrtlFlux.h