genie::flux::GAstroFlux::NuGenerator Class Reference

#include <GAstroFlux.h>

Public Member Functions
	NuGenerator ()
	~NuGenerator ()
bool	SelectNuPdg (bool weighted, const map< int, double > &nupdgpdf, int &nupdg, double &wght)
bool	SelectEnergy (bool weighted, TH1D &log10epdf, double log10emin, double log10emax, double &log10e, double &wght)
bool	SelectOrigin (bool weighted, TH2D &opdf, double &phi, double &costheta, double &wght)

Detailed Description

Definition at line 199 of file GAstroFlux.h.

Constructor & Destructor Documentation

genie::flux::GAstroFlux::NuGenerator::NuGenerator ( )

inline

Definition at line 201 of file GAstroFlux.h.

{}

genie::flux::GAstroFlux::NuGenerator::~NuGenerator ( )

inline

Definition at line 202 of file GAstroFlux.h.

{}

Member Function Documentation

bool GAstroFlux::NuGenerator::SelectEnergy	(	bool	weighted,
		TH1D &	log10epdf,
		double	log10emin,
		double	log10emax,
		double &	log10e,
		double &	wght
	)

Definition at line 411 of file GAstroFlux.cxx.

References genie::RandomGen::Instance(), and genie::RandomGen::RndFlux().

{
// select neutrino energy
//

  log10E   = -9999999;
  wght     = 0;

  if(log10Emax <= log10Emin) {
    return false;
  }

  // Generate weighted flux:
  //
  if(weighted) {
     RandomGen * rnd = RandomGen::Instance();
     log10E  = log10Emin + (log10Emax-log10Emin) * rnd->RndFlux().Rndm();
     wght    = log10Epdf.GetBinContent(log10Epdf.FindBin(log10E));
  }

  // Generate un-weighted flux:
  //
  else {
     do {
       log10E = log10Epdf.GetRandom();
     }
     while(log10E < log10Emin || log10E > log10Emax);
     wght = 1.;
  }

  return true;
}

bool GAstroFlux::NuGenerator::SelectNuPdg	(	bool	weighted,
		const map< int, double > &	nupdgpdf,
		int &	nupdg,
		double &	wght
	)

Definition at line 363 of file GAstroFlux.cxx.

References genie::RandomGen::Instance(), and genie::RandomGen::RndFlux().

{
// select neutrino species based on relative neutrino species populations
//
  nupdg = 0;
  wght  = 0;

  if(nupdgpdf.size() == 0) {
    return false;
  }

  RandomGen * rnd = RandomGen::Instance();

  // Generate weighted flux:
  //
  if(weighted) {
     unsigned int nnu = nupdgpdf.size();
     unsigned int inu = rnd->RndFlux().Integer(nnu);
     map<int,double>::const_iterator iter = nupdgpdf.begin();
     advance(iter,inu);
     nupdg = iter->first;
     wght  = iter->second;
  }
  // Generate un-weighted flux:
  //
  else {
     double xsum = 0.;
     double xrnd = rnd->RndFlux().Uniform();
     map<int,double>::const_iterator iter = nupdgpdf.begin();
     for( ; iter != nupdgpdf.end(); ++iter) {
       xsum += iter->second;
       if(xrnd < xsum) {
         nupdg = iter->first;
         break;
       }
     }
     wght = 1.;
  }

  if(nupdg==0) {
    return false;
  }

  return true;
}

bool GAstroFlux::NuGenerator::SelectOrigin	(	bool	weighted,
		TH2D &	opdf,
		double &	phi,
		double &	costheta,
		double &	wght
	)

Definition at line 446 of file GAstroFlux.cxx.

References genie::RandomGen::Instance(), genie::constants::kPi, and genie::RandomGen::RndFlux().

{
  wght     = 0;
  costheta = -999999;
  phi      = -999999;

  // Generate weighted flux:
  //
  if(weighted) {
     RandomGen * rnd = RandomGen::Instance();
     phi      = 2.*kPi * rnd->RndFlux().Rndm();
     costheta = -1. + 2.*rnd->RndFlux().Rndm();
     wght     = opdf.GetBinContent(opdf.FindBin(phi,costheta));
  }

  // Generate un-weighted flux:
  //
  else {
     opdf.GetRandom2(phi,costheta);
     wght = 1.;
  }

  return true;
}

---

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

• GAstroFlux.h
• GAstroFlux.cxx