#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) More...

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) More...

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

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

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

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

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

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

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

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

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

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

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

Private Member Functions
bool	SelectSource (void)

Private Attributes
map< int, string >	fPntSrcName
	point source name More...

map< int, double >	fPntSrcRA
	right ascension More...

map< int, double >	fPntSrcDec
	declination More...

map< int, double >	fPntSrcRelI
	relative intensity More...

double	fPntSrcTotI
	sum of all relative intensities More...

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

int	fgPdgC
	(current) generated nu pdg-code More...

TLorentzVector	fgP4
	(current) generated nu 4-momentum More...

TLorentzVector	fgX4
	(current) generated nu 4-position More...

double	fgWeight
	(current) generated nu weight More...

double	fMaxEvCut
	(config) user-defined maximum energy cut More...

double	fMinEvCut
	(config) user-defined minimum energy cut More...

bool	fGenWeighted
	(config) generate a weighted or unweighted flux? More...

double	fDetGeoLatitude
	(config) detector: geographic latitude More...

double	fDetGeoLongitude
	(config) detector: geographic longitude More...

double	fDetGeoDepth
	(config) detector: depth from surface More...

double	fDetSize
	(config) detector: size (detector should be enclosed in sphere of this radius) More...

map< int, double >	fRelNuPopulations
	(config) relative neutrino populations More...

TRotation	fRotGEF2THz
	(config) coord. system rotation: GEF translated to detector centre -> THZ More...

TRotation	fRotTHz2User
	(config) coord. system rotation: THZ -> Topocentric user-defined More...

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 ( )

Definition at line 553 of file GAstroFlux.cxx.

References clear, fPntSrcDec, fPntSrcName, fPntSrcRA, fPntSrcRelI, and fPntSrcTotI.

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

GPointSourceAstroFlux::~GPointSourceAstroFlux ( )

Definition at line 564 of file GAstroFlux.cxx.

 {
 
 }

Member Function Documentation

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

Definition at line 574 of file GAstroFlux.cxx.

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

 {
   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;
   }
 }

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;
 }

bool GPointSourceAstroFlux::SelectSource ( void )

private

Definition at line 595 of file GAstroFlux.cxx.

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

 {
   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;
 }