genie::flux::GPointSourceAstroFlux Class Reference

#include <GAstroFlux.h>

Inheritance diagram for genie::flux::GPointSourceAstroFlux:

Collaboration diagram for genie::flux::GPointSourceAstroFlux:

Public Member Functions
	GPointSourceAstroFlux ()
	~GPointSourceAstroFlux ()
bool	GenerateNext (void)
	generate the next flux neutrino (return false in err)
void	AddPointSource (string name, double ra, double dec, double rel_intensity)
Public Member Functions inherited from genie::flux::GAstroFlux
virtual	~GAstroFlux ()
virtual const PDGCodeList &	FluxParticles (void)
	declare list of flux neutrinos that can be generated (for init. purposes)
virtual double	MaxEnergy (void)
	declare the max flux neutrino energy that can be generated (for init. purposes)
virtual int	PdgCode (void)
	returns the flux neutrino pdg code
virtual double	Weight (void)
	returns the flux neutrino weight (if any)
virtual const TLorentzVector &	Momentum (void)
	returns the flux neutrino 4-momentum
virtual const TLorentzVector &	Position (void)
	returns the flux neutrino 4-position (note: expect SI rather than physical units)
virtual bool	End (void)
	true if no more flux nu's can be thrown (eg reaching end of beam sim ntuples)
virtual long int	Index (void)
	returns corresponding index for current flux neutrino (e.g. for a flux ntuple returns the current entry number)
virtual void	Clear (Option_t *opt)
	reset state variables based on opt
virtual void	GenerateWeighted (bool gen_weighted)
	set whether to generate weighted or unweighted neutrinos
void	ForceMinEnergy (double emin)
void	ForceMaxEnergy (double emax)
void	SetDetectorPosition (double latitude, double longitude, double depth, double size)
void	SetRelNuPopulations (double nnue=1, double nnumu=2, double nnutau=0, double nnuebar=1, double nnumubar=2, double nnutaubar=0)
void	SetEnergyPowLawIdx (double n)
void	SetUserCoordSystem (TRotation &rotation)
	rotation Topocentric Horizontal -> User-defined Topocentric Coord System
Public Member Functions inherited from genie::GFluxI
virtual	~GFluxI ()

Private Member Functions
bool	SelectSource (void)

Private Attributes
map< int, string >	fPntSrcName
	point source name
map< int, double >	fPntSrcRA
	right ascension
map< int, double >	fPntSrcDec
	declination
map< int, double >	fPntSrcRelI
	relative intensity
double	fPntSrcTotI
	sum of all relative intensities
unsigned int	fSelSourceId

Additional Inherited Members
Protected Member Functions inherited from genie::flux::GAstroFlux
	GAstroFlux ()
void	Initialize (void)
void	CleanUp (void)
void	ResetSelection (void)
Protected Member Functions inherited from genie::GFluxI
	GFluxI ()
Protected Attributes inherited from genie::flux::GAstroFlux
PDGCodeList *	fPdgCList
	declared list of neutrino pdg-codes that can be thrown by current instance
int	fgPdgC
	(current) generated nu pdg-code
TLorentzVector	fgP4
	(current) generated nu 4-momentum
TLorentzVector	fgX4
	(current) generated nu 4-position
double	fgWeight
	(current) generated nu weight
double	fMaxEvCut
	(config) user-defined maximum energy cut
double	fMinEvCut
	(config) user-defined minimum energy cut
bool	fGenWeighted
	(config) generate a weighted or unweighted flux?
double	fDetGeoLatitude
	(config) detector: geographic latitude
double	fDetGeoLongitude
	(config) detector: geographic longitude
double	fDetGeoDepth
	(config) detector: depth from surface
double	fDetSize
	(config) detector: size (detector should be enclosed in sphere of this radius)
map< int, double >	fRelNuPopulations
	(config) relative neutrino populations
TRotation	fRotGEF2THz
	(config) coord. system rotation: GEF translated to detector centre -> THZ
TRotation	fRotTHz2User
	(config) coord. system rotation: THZ -> Topocentric user-defined
TVector3	fDetCenter
TH1D *	fEnergySpectrum
TH2D *	fSolidAngleAcceptance
NuGenerator *	fNuGen
NuPropagator *	fNuPropg

Detailed Description

Definition at line 244 of file GAstroFlux.h.

Constructor & Destructor Documentation

GPointSourceAstroFlux()

GPointSourceAstroFlux::GPointSourceAstroFlux

(

)

Definition at line 553 of file GAstroFlux.cxx.

                                             :
GAstroFlux()
{
  fPntSrcName.clear();
  fPntSrcRA.  clear();
  fPntSrcDec. clear();
  fPntSrcRelI.clear();
 
  fPntSrcTotI = 0;
}

References clear, fPntSrcDec, fPntSrcName, fPntSrcRA, fPntSrcRelI, fPntSrcTotI, and genie::flux::GAstroFlux::GAstroFlux().

~GPointSourceAstroFlux()

GPointSourceAstroFlux::~GPointSourceAstroFlux

(

)

Definition at line 564 of file GAstroFlux.cxx.

{
 
}

Member Function Documentation

AddPointSource()

void GPointSourceAstroFlux::AddPointSource	(	string	name,
		double	ra,
		double	dec,
		double	rel_intensity )

Definition at line 574 of file GAstroFlux.cxx.

{
  bool accept =
         (ra  >= 0.      && ra  < 2.*kPi)  &&
         (dec >= -kPi/2. && dec <= kPi/2.) &&
         (rel_intensity > 0) &&
         (name.size() > 0);
 
  if(accept) {
     int id = fPntSrcName.size();
 
     fPntSrcName.insert( map<int, string>::value_type(id, name         ) );
     fPntSrcRA.  insert( map<int, double>::value_type(id, ra           ) );
     fPntSrcDec. insert( map<int, double>::value_type(id, dec          ) );
     fPntSrcRelI.insert( map<int, double>::value_type(id, rel_intensity) );
 
     fPntSrcTotI += rel_intensity;
  }
}

References fPntSrcDec, fPntSrcName, fPntSrcRA, fPntSrcRelI, fPntSrcTotI, and genie::constants::kPi.

GenerateNext()

bool GPointSourceAstroFlux::GenerateNext ( void )

virtual

generate the next flux neutrino (return false in err)

Reimplemented from genie::flux::GAstroFlux.

Definition at line 569 of file GAstroFlux.cxx.

{
  return true;
}

SelectSource()

bool GPointSourceAstroFlux::SelectSource ( void )

private

Definition at line 595 of file GAstroFlux.cxx.

{
  if(fPntSrcRelI.size() == 0) {
    return false;
  }
 
  if(fPntSrcTotI <= 0.) {
    return false;
  }
 
  unsigned int srcid = 0;
  double wght  = 0;
 
  // Generate weighted flux:
  //
  if(fGenWeighted) {
     RandomGen * rnd = RandomGen::Instance();
     unsigned int nsrc = fPntSrcName.size();
     srcid = rnd->RndFlux().Integer(nsrc);
     wght  = fPntSrcRelI[srcid] / fPntSrcTotI;
  }
  // Generate un-weighted flux:
  //
  else {
     RandomGen * rnd = RandomGen::Instance();
     double xsum = 0.;
     double xrnd = fPntSrcTotI * rnd->RndFlux().Uniform();
     map<int,double>::const_iterator piter = fPntSrcRelI.begin();
     for( ; piter != fPntSrcRelI.end(); ++piter) {
       xsum += piter->second;
       if(xrnd < xsum) {
         srcid = piter->first;
         break;
       }
     }
     wght = 1.;
  }
 
  //
  fSelSourceId  = srcid;
 
  //
  fgWeight *= wght;
 
  return true;
}

References genie::flux::GAstroFlux::fGenWeighted, genie::flux::GAstroFlux::fgWeight, fPntSrcName, fPntSrcRelI, fPntSrcTotI, fSelSourceId, genie::RandomGen::Instance(), and genie::RandomGen::RndFlux().

Member Data Documentation

fPntSrcDec

map<int, double> genie::flux::GPointSourceAstroFlux::fPntSrcDec

private

declination

Definition at line 261 of file GAstroFlux.h.

Referenced by AddPointSource(), and GPointSourceAstroFlux().

fPntSrcName

map<int, string> genie::flux::GPointSourceAstroFlux::fPntSrcName

private

point source name

Definition at line 259 of file GAstroFlux.h.

Referenced by AddPointSource(), GPointSourceAstroFlux(), and SelectSource().

fPntSrcRA

map<int, double> genie::flux::GPointSourceAstroFlux::fPntSrcRA

private

right ascension

Definition at line 260 of file GAstroFlux.h.

Referenced by AddPointSource(), and GPointSourceAstroFlux().

fPntSrcRelI

map<int, double> genie::flux::GPointSourceAstroFlux::fPntSrcRelI

private

relative intensity

Definition at line 262 of file GAstroFlux.h.

Referenced by AddPointSource(), GPointSourceAstroFlux(), and SelectSource().

fPntSrcTotI

double genie::flux::GPointSourceAstroFlux::fPntSrcTotI

private

sum of all relative intensities

Definition at line 263 of file GAstroFlux.h.

Referenced by AddPointSource(), GPointSourceAstroFlux(), and SelectSource().

fSelSourceId

unsigned int genie::flux::GPointSourceAstroFlux::fSelSourceId

private

Definition at line 265 of file GAstroFlux.h.

Referenced by SelectSource().

---

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

• GAstroFlux.h
• GAstroFlux.cxx