A flux driver for the FLUKA 3-D Atmospheric Neutrino Flux. More...

#include <GFLUKAAtmoFlux.h>

Inheritance diagram for genie::flux::GFLUKAAtmoFlux:

Collaboration diagram for genie::flux::GFLUKAAtmoFlux:

Public Member Functions
	GFLUKAAtmoFlux ()

	~GFLUKAAtmoFlux ()

Public Member Functions inherited from genie::flux::GAtmoFlux
virtual	~GAtmoFlux ()

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 bool	GenerateNext (void)
	generate the next flux neutrino (return false in err) 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...

double	Enu (void)

double	Energy (void)

double	CosTheta (void)

double	CosZenith (void)

long int	NFluxNeutrinos (void) const
	Number of flux nu's generated. Not the same as the number of nu's thrown towards the geometry (if there are cuts). More...

void	ForceMinEnergy (double emin)

void	ForceMaxEnergy (double emax)

void	SetSpectralIndex (double index)

void	SetRadii (double Rlongitudinal, double Rtransverse)

double	GetFluxSurfaceArea (void)

double	GetLongitudinalRadius (void)

double	GetTransverseRadius (void)

void	SetUserCoordSystem (TRotation &rotation)
	Rotation: Topocentric Horizontal -> User-defined Topocentric Coord System. More...

void	AddFluxFile (int neutrino_pdg, string filename)

void	AddFluxFile (string filename)

bool	LoadFluxData (void)

TH3D *	GetFluxHistogram (int flavour)

double	GetTotalFlux (void)

double	GetTotalFluxInEnergyRange (void)

double	GetFlux (int flavour)

double	GetFlux (int flavour, double energy)

double	GetFlux (int flavour, double energy, double costh)

double	GetFlux (int flavour, double energy, double costh, double phi)

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

Private Member Functions
void	SetBinSizes (void)

bool	FillFluxHisto (int nu_pdg, string filename)

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

bool	GenerateNext_1try (void)

void	Initialize (void)

void	CleanUp (void)

void	ResetSelection (void)

double	MinEnergy (void)

TH3D *	CreateFluxHisto (string name, string title)

void	ZeroFluxHisto (TH3D *hist)

void	AddAllFluxes (void)

int	SelectNeutrino (double Ev, double costheta, double phi)

TH3D *	CreateNormalisedFluxHisto (TH3D *hist)

Protected Member Functions inherited from genie::GFluxI
	GFluxI ()

Protected Attributes inherited from genie::flux::GAtmoFlux
double	fMaxEv
	maximum energy (in input flux files) More...

PDGCodeList *	fPdgCList
	input list of neutrino pdg-codes 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	fWeight
	current generated nu weight More...

long int	fNNeutrinos
	number of flux neutrinos thrown so far More...

double	fMaxEvCut
	user-defined cut: maximum energy More...

double	fMinEvCut
	user-defined cut: minimum energy More...

double	fRl
	defining flux neutrino generation surface: longitudinal radius More...

double	fRt
	defining flux neutrino generation surface: transverse radius More...

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

unsigned int	fNumPhiBins
	number of phi bins in input flux data files More...

unsigned int	fNumCosThetaBins
	number of cos(theta) bins in input flux data files More...

unsigned int	fNumEnergyBins
	number of energy bins in input flux data files More...

double *	fPhiBins
	phi bins in input flux data files More...

double *	fCosThetaBins
	cos(theta) bins in input flux data files More...

double *	fEnergyBins
	energy bins in input flux data files More...

bool	fGenWeighted
	generate a weighted or unweighted flux? More...

double	fSpectralIndex
	power law function used for weighted flux More...

bool	fInitialized
	flag to check that initialization is run More...

TH3D *	fTotalFluxHisto
	flux = f(Ev,cos8,phi) summed over neutrino species More...

double	fTotalFluxHistoIntg
	fFluxSum2D integral More...

map< int, TH3D * >	fFluxHistoMap
	flux = f(Ev,cos8,phi) for each neutrino species More...

map< int, TH3D * >	fRawFluxHistoMap
	flux = f(Ev,cos8,phi) for each neutrino species More...

vector< int >	fFluxFlavour
	input flux file for each neutrino species More...

vector< string >	fFluxFile
	input flux file for each neutrino species More...

Detailed Description

A flux driver for the FLUKA 3-D Atmospheric Neutrino Flux.

References:: Astrop.Phys.19 (2003) p.269; hep-ph/0207035; hep-ph/9907408 Alfredo.Ferrari Alfre.nosp@m.do.F.nosp@m.errar.nosp@m.i@ce.nosp@m.rn.ch Paola.Sala Paola.nosp@m..Sal.nosp@m.a@cer.nosp@m.n.ch Giuseppe Battistoni Giuse.nosp@m.ppe..nosp@m.Batti.nosp@m.ston.nosp@m.i@mi..nosp@m.infn.nosp@m..it Teresa Montaruli Teres.nosp@m.a.Mo.nosp@m.ntaru.nosp@m.li@b.nosp@m.a.inf.nosp@m.n.it

To be able to use this flux driver you will need to download the flux data from: http://pcbat1.mi.infn.it/~battist/neutrino.html

Please note that this class expects to read flux files formatted as described in the above FLUKA flux page. Each file contains 3 columns:

neutrino energy (GeV) at bin centre
neutrino cos(zenith angle) at bin centre
neutrino flux (#neutrinos /GeV /m^2 /sec /sr) The flux is given in 40 bins of cos(zenith angle) from -1.0 to 1.0 (bin width = 0.05) and 61 equally log-spaced energy bins (20 bins per decade), with Emin = 0.100 GeV.

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

Created:: July 3, 2005

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

Definition at line 54 of file GFLUKAAtmoFlux.h.

Constructor & Destructor Documentation

GFLUKAAtmoFlux::GFLUKAAtmoFlux ( )

Definition at line 31 of file GFLUKAAtmoFlux.cxx.

References Initialize(), LOG, and pNOTICE.

                                :
 GAtmoFlux()
 {
   LOG("Flux", pNOTICE)
     << "Instantiating the GENIE FLUKA atmospheric neutrino flux driver";
 
   this->Initialize();
   this->SetBinSizes();
 }

GFLUKAAtmoFlux::~GFLUKAAtmoFlux ( )

Definition at line 41 of file GFLUKAAtmoFlux.cxx.

 {
 
 }

Member Function Documentation

bool GFLUKAAtmoFlux::FillFluxHisto	(	int	nu_pdg,
		string	filename
	)

privatevirtual

Implements genie::flux::GAtmoFlux.

Definition at line 109 of file GFLUKAAtmoFlux.cxx.

References genie::flux::GAtmoFlux::fRawFluxHistoMap, genie::constants::kPi, LOG, pERROR, pINFO, and pNOTICE.

 {
   LOG("Flux", pNOTICE)
     << "Loading FLUKA atmospheric flux for neutrino: " << nu_pdg
     << " from file: " << filename;
 
   TH3D* histo = 0;
   std::map<int,TH3D*>::iterator myMapEntry = fRawFluxHistoMap.find(nu_pdg);
   if( myMapEntry != fRawFluxHistoMap.end() ){
       histo = myMapEntry->second;
   }
   if(!histo) {
      LOG("Flux", pERROR) << "Null flux histogram!";
      return false;
   }
   ifstream flux_stream(filename.c_str(), ios::in);
   if(!flux_stream) {
      LOG("Flux", pERROR) << "Could not open file: " << filename;
      return false;
   }
 
   int    ibin;
   double energy, costheta, flux;
   char   j1, j2;
 
   double scale = 1.0; // 1.0 [m^2], OR 1.0e-4 [cm^2]
 
   while ( !flux_stream.eof() ) {
     flux = 0.0;
     flux_stream >> energy >> j1 >> costheta >> j2 >> flux;
     if( flux>0.0 ){
       LOG("Flux", pINFO)
         << "Flux[Ev = " << energy
         << ", cos8 = " << costheta << "] = " << flux;
       // note: reversing the Fluka sign convention for zenith angle
       //       1 phi bin
       ibin = histo->FindBin( (Axis_t)energy, (Axis_t)(-costheta), (Axis_t)kPi );
       histo->SetBinContent( ibin, (Stat_t)(scale*flux) );
     }
   }
   return true;
 }

void GFLUKAAtmoFlux::SetBinSizes ( void )

private

Definition at line 46 of file GFLUKAAtmoFlux.cxx.

References genie::flux::GAtmoFlux::fCosThetaBins, genie::flux::GAtmoFlux::fEnergyBins, genie::flux::GAtmoFlux::fMaxEv, genie::flux::GAtmoFlux::fNumCosThetaBins, genie::flux::GAtmoFlux::fNumEnergyBins, genie::flux::GAtmoFlux::fNumPhiBins, genie::flux::GAtmoFlux::fPhiBins, genie::flux::kGFlk3DCosThetaMax, genie::flux::kGFlk3DCosThetaMin, genie::flux::kGFlk3DEvMin, genie::flux::kGFlk3DNumCosThetaBins, genie::flux::kGFlk3DNumLogEvBins, genie::flux::kGFlk3DNumLogEvBinsPerDecade, genie::constants::kPi, LOG, and pDEBUG.

 {
 // Generate the correct cos(theta) and energy bin sizes
 // The flux is given in 40 bins of cos(zenith angle) from -1.0 to 1.0
 // (bin width = 0.05) and 61 equally log-spaced energy bins (20 bins
 // per decade), with Emin = 0.100 GeV.
 //
 
   fPhiBins       = new double [2];
   fCosThetaBins  = new double [kGFlk3DNumCosThetaBins  + 1];
   fEnergyBins    = new double [kGFlk3DNumLogEvBins     + 1];
 
   fPhiBins[0] = 0;
   fPhiBins[1] = 2.*kPi;
 
   double dcostheta =
       (kGFlk3DCosThetaMax - kGFlk3DCosThetaMin) /
       (double) kGFlk3DNumCosThetaBins;
 
   double logEmax = TMath::Log10(1.);
   double logEmin = TMath::Log10(kGFlk3DEvMin);
   double dlogE =
       (logEmax - logEmin) /
       (double) kGFlk3DNumLogEvBinsPerDecade;
 
   for(unsigned int i=0; i<= kGFlk3DNumCosThetaBins; i++) {
      fCosThetaBins[i] = kGFlk3DCosThetaMin + i * dcostheta;
      if(i != kGFlk3DNumCosThetaBins) {
        LOG("Flux", pDEBUG)
          << "FLUKA 3d flux: CosTheta bin " << i+1
          << ": lower edge = " << fCosThetaBins[i];
      } else {
        LOG("Flux", pDEBUG)
          << "FLUKA 3d flux: CosTheta bin " << kGFlk3DNumCosThetaBins
          << ": upper edge = " << fCosThetaBins[kGFlk3DNumCosThetaBins];
      }
   }
 
   for(unsigned int i=0; i<= kGFlk3DNumLogEvBins; i++) {
      fEnergyBins[i] = TMath::Power(10., logEmin + i*dlogE);
      if(i != kGFlk3DNumLogEvBins) {
        LOG("Flux", pDEBUG)
          << "FLUKA 3d flux: Energy bin " << i+1
          << ": lower edge = " << fEnergyBins[i];
      } else {
        LOG("Flux", pDEBUG)
          << "FLUKA 3d flux: Energy bin " << kGFlk3DNumLogEvBins
          << ": upper edge = " << fEnergyBins[kGFlk3DNumLogEvBins];
      }
   }
 
   for(unsigned int i=0; i< kGFlk3DNumLogEvBins; i++) {
        LOG("Flux", pDEBUG)
          << "FLUKA 3d flux: Energy bin " << i+1
          << ": bin centre = " << (fEnergyBins[i] + fEnergyBins[i+1])/2.;
   }
 
   fNumPhiBins      = 1;
   fNumCosThetaBins = kGFlk3DNumCosThetaBins;
   fNumEnergyBins   = kGFlk3DNumLogEvBins;
   fMaxEv           = fEnergyBins[fNumEnergyBins];
 }

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

Public Member Functions

Private Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation