A simple GENIE flux driver for monoenergetic neutrinos along the z direction. Can handle a mix of neutrinos with their corresponding weight. This trivial flux driver was added in so that single energy neutrinos can be easily used with the event generation driver that can handle a target mix or detailed geometries. More...

#include <GMonoEnergeticFlux.h>

Inheritance diagram for genie::flux::GMonoEnergeticFlux:

[legend]

Collaboration diagram for genie::flux::GMonoEnergeticFlux:

[legend]

Public Member Functions
	GMonoEnergeticFlux ()

	GMonoEnergeticFlux (double Ev, int pdg)

	GMonoEnergeticFlux (double Ev, const map< int, double > &numap)

	~GMonoEnergeticFlux ()

const PDGCodeList &	FluxParticles (void)
	declare list of flux neutrinos that can be generated (for init. purposes) More...

double	MaxEnergy (void)
	declare the max flux neutrino energy that can be generated (for init. purposes) More...

bool	GenerateNext (void)
	generate the next flux neutrino (return false in err) More...

int	PdgCode (void)
	returns the flux neutrino pdg code More...

double	Weight (void)
	returns the flux neutrino weight (if any) More...

const TLorentzVector &	Momentum (void)
	returns the flux neutrino 4-momentum More...

const TLorentzVector &	Position (void)
	returns the flux neutrino 4-position (note: expect SI rather than physical units) More...

bool	End (void)
	true if no more flux nu's can be thrown (eg reaching end of beam sim ntuples) More...

long int	Index (void)
	returns corresponding index for current flux neutrino (e.g. for a flux ntuple returns the current entry number) More...

void	Clear (Option_t *opt)
	reset state variables based on opt More...

void	GenerateWeighted (bool gen_weighted)
	set whether to generate weighted or unweighted neutrinos More...

void	SetDirectionCos (double dx, double dy, double dz)

void	SetRayOrigin (double x, double y, double z)

void	SetNuDirection (const TVector3 &direction)

void	SetBeamSpot (const TVector3 &spot)

void	Initialize (double Ev, int pdg)

void	Initialize (double Ev, const map< int, double > &numap)

Public Member Functions inherited from genie::GFluxI
virtual	~GFluxI ()

Private Member Functions
void	CleanUp (void)

Private Attributes
double	fMaxEv
	maximum energy More...

PDGCodeList *	fPdgCList
	list of neutrino pdg-codes More...

int	fgPdgC
	running generated nu pdg-code More...

TLorentzVector	fgP4
	running generated nu 4-momentum More...

TLorentzVector	fgX4
	running generated nu 4-position More...

map< int, double >	fProb

double	fProbMax

Additional Inherited Members
Protected Member Functions inherited from genie::GFluxI
	GFluxI ()

Detailed Description

A simple GENIE flux driver for monoenergetic neutrinos along the z direction. Can handle a mix of neutrinos with their corresponding weight. This trivial flux driver was added in so that single energy neutrinos can be easily used with the event generation driver that can handle a target mix or detailed geometries.

Author: Costas Andreopoulos <c.andreopoulos cern.ch> University of Liverpool

Created:: Feb 08, 2008

License:: Copyright (c) 2003-2024, The GENIE Collaboration For the full text of the license visit http://copyright.genie-mc.org

Definition at line 38 of file GMonoEnergeticFlux.h.

Constructor & Destructor Documentation

GMonoEnergeticFlux::GMonoEnergeticFlux ( )

Definition at line 28 of file GMonoEnergeticFlux.cxx.

                                        :
 GFluxI()
 {
   // default ctor for consistency with GFluxDriverFactory needs
   // up to user to call Initialize() to set energy and flavor(s)
 }

GMonoEnergeticFlux::GMonoEnergeticFlux	(	double	Ev,
		int	pdg
	)

Definition at line 35 of file GMonoEnergeticFlux.cxx.

References Initialize().

                                                          :
 GFluxI()
 {
   this->Initialize(Ev,pdg);
 }

GMonoEnergeticFlux::GMonoEnergeticFlux	(	double	Ev,
		const map< int, double > &	numap
	)

Definition at line 41 of file GMonoEnergeticFlux.cxx.

References Initialize().

                                                  :
 GFluxI()
 {
   this->Initialize(Ev,numap);
 }

GMonoEnergeticFlux::~GMonoEnergeticFlux ( )

Definition at line 48 of file GMonoEnergeticFlux.cxx.

References CleanUp().

 {
   this->CleanUp();
 }

Member Function Documentation

void GMonoEnergeticFlux::CleanUp ( void )

private

Definition at line 130 of file GMonoEnergeticFlux.cxx.

References fPdgCList, LOG, and pNOTICE.

Referenced by ~GMonoEnergeticFlux().

 {
   LOG("Flux", pNOTICE) << "Cleaning up...";
 
   if (fPdgCList) delete fPdgCList;
 }

void GMonoEnergeticFlux::Clear ( Option_t * opt )

virtual

reset state variables based on opt

Implements genie::GFluxI.

Definition at line 77 of file GMonoEnergeticFlux.cxx.

References LOG, and pERROR.

 {
 // Dummy clear method needed to conform to GFluxI interface
 //
   LOG("Flux", pERROR) <<
       "No Clear(Option_t * opt) method implemented for opt: "<< opt;
 }

bool genie::flux::GMonoEnergeticFlux::End ( void )

inlinevirtual

true if no more flux nu's can be thrown (eg reaching end of beam sim ntuples)

Implements genie::GFluxI.

Definition at line 54 of file GMonoEnergeticFlux.h.

54 { return false; }

const PDGCodeList& genie::flux::GMonoEnergeticFlux::FluxParticles ( void )

inlinevirtual

declare list of flux neutrinos that can be generated (for init. purposes)

Implements genie::GFluxI.

Definition at line 47 of file GMonoEnergeticFlux.h.

References fPdgCList.

47 { return *fPdgCList; }

genie::flux::GMonoEnergeticFlux::fPdgCList

PDGCodeList * fPdgCList

list of neutrino pdg-codes

Definition: GMonoEnergeticFlux.h:77

bool GMonoEnergeticFlux::GenerateNext ( void )

virtual

generate the next flux neutrino (return false in err)

Implements genie::GFluxI.

Definition at line 53 of file GMonoEnergeticFlux.cxx.

References fgP4, fgPdgC, fgX4, fProb, fProbMax, genie::RandomGen::Instance(), LOG, genie::utils::print::P4AsShortString(), pINFO, genie::RandomGen::RndFlux(), and genie::utils::print::X4AsString().

 {
   RandomGen * rnd = RandomGen::Instance();
   double p = fProbMax * rnd->RndFlux().Rndm();
 
   map<int,double>::const_iterator iter;
   for(iter = fProb.begin(); iter != fProb.end(); ++iter) {
      int    nupdgc = iter->first;
      double prob   = iter->second;
      if (p<prob) {
        fgPdgC = nupdgc;
        break;
      }
   }
 
   LOG("Flux", pINFO)
         << "Generated neutrino: "
         << "\n pdg-code: " << fgPdgC
         << "\n p4: " << utils::print::P4AsShortString(&fgP4)
         << "\n x4: " << utils::print::X4AsString(&fgX4);
 
   return true;
 }

void GMonoEnergeticFlux::GenerateWeighted ( bool gen_weighted )

virtual

set whether to generate weighted or unweighted neutrinos

Implements genie::GFluxI.

Definition at line 85 of file GMonoEnergeticFlux.cxx.

References LOG, and pERROR.

 {
 // Dummy implementation needed to conform to GFluxI interface
 //
   LOG("Flux", pERROR) <<
       "No GenerateWeighted(bool gen_weighted) method implemented for " <<
       "gen_weighted: " << gen_weighted;
 }

long int genie::flux::GMonoEnergeticFlux::Index ( void )

inlinevirtual

returns corresponding index for current flux neutrino (e.g. for a flux ntuple returns the current entry number)

Implements genie::GFluxI.

Definition at line 55 of file GMonoEnergeticFlux.h.

55 { return -1; }

void GMonoEnergeticFlux::Initialize	(	double	Ev,
		int	pdg
	)

Definition at line 94 of file GMonoEnergeticFlux.cxx.

Referenced by GMonoEnergeticFlux().

 {
   map<int,double> numap;
   numap.insert( map<int, double>::value_type(pdg, 1.) );
 
   this->Initialize(Ev,numap);
 }

void GMonoEnergeticFlux::Initialize	(	double	Ev,
		const map< int, double > &	numap
	)

Definition at line 102 of file GMonoEnergeticFlux.cxx.

References fgP4, fgPdgC, fgX4, fMaxEv, fPdgCList, fProb, fProbMax, LOG, pNOTICE, and genie::PDGCodeList::push_back().

 {
   LOG("Flux", pNOTICE) << "Initializing GMonoEnergeticFlux driver";
 
   fMaxEv = Ev + 0.05;
 
   fPdgCList = new PDGCodeList;
   fPdgCList->clear();
 
   fProbMax = 0;
   fProb.clear();
 
   map<int,double>::const_iterator iter;
   for(iter = numap.begin(); iter != numap.end(); ++iter) {
         int    nupdgc = iter->first;
         double nuwgt  = iter->second;
 
         fPdgCList->push_back(nupdgc);
 
         fProbMax+=nuwgt;
         fProb.insert(map<int, double>::value_type(nupdgc,fProbMax));
   }
 
   fgPdgC = 0;
   fgP4.SetPxPyPzE (0.,0.,Ev,Ev);
   fgX4.SetXYZT    (0.,0.,0.,0.);
 }

double genie::flux::GMonoEnergeticFlux::MaxEnergy ( void )

inlinevirtual

declare the max flux neutrino energy that can be generated (for init. purposes)

Implements genie::GFluxI.

Definition at line 48 of file GMonoEnergeticFlux.h.

References fMaxEv.

48 { return fMaxEv; }

genie::flux::GMonoEnergeticFlux::fMaxEv

double fMaxEv

maximum energy

Definition: GMonoEnergeticFlux.h:76

const TLorentzVector& genie::flux::GMonoEnergeticFlux::Momentum ( void )

inlinevirtual

returns the flux neutrino 4-momentum

Implements genie::GFluxI.

Definition at line 52 of file GMonoEnergeticFlux.h.

References fgP4.

52 { return fgP4; }

genie::flux::GMonoEnergeticFlux::fgP4

TLorentzVector fgP4

running generated nu 4-momentum

Definition: GMonoEnergeticFlux.h:79

int genie::flux::GMonoEnergeticFlux::PdgCode ( void )

inlinevirtual

returns the flux neutrino pdg code

Implements genie::GFluxI.

Definition at line 50 of file GMonoEnergeticFlux.h.

References fgPdgC.

50 { return fgPdgC; }

genie::flux::GMonoEnergeticFlux::fgPdgC

int fgPdgC

running generated nu pdg-code

Definition: GMonoEnergeticFlux.h:78

const TLorentzVector& genie::flux::GMonoEnergeticFlux::Position ( void )

inlinevirtual

returns the flux neutrino 4-position (note: expect SI rather than physical units)

Implements genie::GFluxI.

Definition at line 53 of file GMonoEnergeticFlux.h.

References fgX4.

53 { return fgX4; }

genie::flux::GMonoEnergeticFlux::fgX4

TLorentzVector fgX4

running generated nu 4-position

Definition: GMonoEnergeticFlux.h:80

void GMonoEnergeticFlux::SetBeamSpot ( const TVector3 & spot )

Definition at line 160 of file GMonoEnergeticFlux.cxx.

References SetRayOrigin().

 {
   SetRayOrigin(spot.x(), spot.y(), spot.z());
 }

void GMonoEnergeticFlux::SetDirectionCos	(	double	dx,
		double	dy,
		double	dz
	)

Definition at line 137 of file GMonoEnergeticFlux.cxx.

References fgP4, LOG, genie::utils::print::P3AsString(), and pNOTICE.

Referenced by SetNuDirection().

 {
   TVector3 dircos1 = TVector3(dx,dy,dz).Unit();
   LOG("Flux", pNOTICE) << "SetDirectionCos "
                        << utils::print::P3AsString(&dircos1);
   double E = fgP4.E();
   fgP4.SetVect(E*dircos1);
 
 }

void GMonoEnergeticFlux::SetNuDirection ( const TVector3 & direction )

Definition at line 155 of file GMonoEnergeticFlux.cxx.

References SetDirectionCos().

 {
   SetDirectionCos(direction.x(), direction.y(), direction.z());
 }

void GMonoEnergeticFlux::SetRayOrigin	(	double	x,
		double	y,
		double	z
	)

Definition at line 147 of file GMonoEnergeticFlux.cxx.

References fgX4, LOG, pNOTICE, and genie::utils::print::Vec3AsString().

Referenced by SetBeamSpot().

 {
   TVector3 xyz(x,y,z);
   LOG("Flux", pNOTICE) << "SetRayOrigin "
                        << utils::print::Vec3AsString(&xyz);
   fgX4.SetVect(xyz);
 }