genie::flux::GMonoEnergeticFlux Class Reference

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:

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 <constantinos.andreopoulos cern.ch> University of Liverpool & STFC Rutherford Appleton Laboratory

Feb 08, 2008

Copyright (c) 2003-2020, 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.

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

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

Definition at line 41 of file GMonoEnergeticFlux.cxx.

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

GMonoEnergeticFlux::~GMonoEnergeticFlux

(

)

Definition at line 48 of file GMonoEnergeticFlux.cxx.

{
  this->CleanUp();
}

Member Function Documentation

void GMonoEnergeticFlux::CleanUp ( void  )

private

Definition at line 130 of file GMonoEnergeticFlux.cxx.

{
  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.

{
// 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.

{ 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.

{ return *fPdgCList; }

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.

{
  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.

{
// 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.

{ return -1;         }

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

Definition at line 94 of file GMonoEnergeticFlux.cxx.

{
  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.

{
  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.

{ return  fMaxEv;    }

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.

{ return  fgP4;      }

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

inlinevirtual

returns the flux neutrino pdg code

Implements genie::GFluxI.

Definition at line 50 of file GMonoEnergeticFlux.h.

{ return  fgPdgC;    }

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.

{ return  fgX4;      }

void GMonoEnergeticFlux::SetBeamSpot

(

const TVector3 &

spot

)

Definition at line 160 of file GMonoEnergeticFlux.cxx.

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

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

Definition at line 137 of file GMonoEnergeticFlux.cxx.

{
  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.

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

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

Definition at line 147 of file GMonoEnergeticFlux.cxx.

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

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

inlinevirtual

returns the flux neutrino weight (if any)

Implements genie::GFluxI.

Definition at line 51 of file GMonoEnergeticFlux.h.

{ return  1.0;       }

Member Data Documentation

TLorentzVector genie::flux::GMonoEnergeticFlux::fgP4

private

running generated nu 4-momentum

Definition at line 79 of file GMonoEnergeticFlux.h.

int genie::flux::GMonoEnergeticFlux::fgPdgC

private

running generated nu pdg-code

Definition at line 78 of file GMonoEnergeticFlux.h.

TLorentzVector genie::flux::GMonoEnergeticFlux::fgX4

private

running generated nu 4-position

Definition at line 80 of file GMonoEnergeticFlux.h.

double genie::flux::GMonoEnergeticFlux::fMaxEv

private

maximum energy

Definition at line 76 of file GMonoEnergeticFlux.h.

PDGCodeList* genie::flux::GMonoEnergeticFlux::fPdgCList

private

list of neutrino pdg-codes

Definition at line 77 of file GMonoEnergeticFlux.h.

map<int, double> genie::flux::GMonoEnergeticFlux::fProb

private

Definition at line 81 of file GMonoEnergeticFlux.h.

double genie::flux::GMonoEnergeticFlux::fProbMax

private

Definition at line 82 of file GMonoEnergeticFlux.h.

---

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

• Generator/src/Tools/Flux/GMonoEnergeticFlux.h
• Generator/src/Tools/Flux/GMonoEnergeticFlux.cxx