GFLUKAAtmoFlux.h
Go to the documentation of this file.
1 //____________________________________________________________________________
2 /*!
3 
4 \class genie::flux::GFLUKAAtmoFlux
5 
6 \brief A flux driver for the FLUKA 3-D Atmospheric Neutrino Flux
7 
8 \ref Astrop.Phys.19 (2003) p.269; hep-ph/0207035; hep-ph/9907408
9  Alfredo.Ferrari <Alfredo.Ferrari@cern.ch>
10  Paola.Sala <Paola.Sala@cern.ch>
11  Giuseppe Battistoni <Giuseppe.Battistoni@mi.infn.it>
12  Teresa Montaruli <Teresa.Montaruli@ba.infn.it>
13 
14  To be able to use this flux driver you will need to download the
15  flux data from: http://pcbat1.mi.infn.it/~battist/neutrino.html
16 
17  Please note that this class expects to read flux files formatted as
18  described in the above FLUKA flux page.
19  Each file contains 3 columns:
20  - neutrino energy (GeV) at bin centre
21  - neutrino cos(zenith angle) at bin centre
22  - neutrino flux (#neutrinos /GeV /m^2 /sec /sr)
23  The flux is given in 40 bins of cos(zenith angle) from -1.0 to 1.0
24  (bin width = 0.05) and 61 equally log-spaced energy bins (20 bins per
25  decade), with Emin = 0.100 GeV.
26 
27 \author Costas Andreopoulos <costas.andreopoulos \at stfc.ac.uk>
28  University of Liverpool & STFC Rutherford Appleton Lab
29 
30 \created July 3, 2005
31 
32 \cpright Copyright (c) 2003-2019, The GENIE Collaboration
33  For the full text of the license visit http://copyright.genie-mc.org
34  or see $GENIE/LICENSE
35 */
36 //____________________________________________________________________________
37 
38 #ifndef _GENIE_FLUKA_ATMO_FLUX_H_
39 #define _GENIE_FLUKA_ATMO_FLUX_H_
40 
41 #include "Tools/Flux/GAtmoFlux.h"
42 
43 namespace genie {
44 namespace flux {
45 
46 // Number of cos(zenith) and energy bins in flux simulation
47 const unsigned int kGFlk3DNumCosThetaBins = 40;
48 const double kGFlk3DCosThetaMin = -1.0;
49 const double kGFlk3DCosThetaMax = 1.0;
50 const unsigned int kGFlk3DNumLogEvBins = 61;
51 const unsigned int kGFlk3DNumLogEvBinsPerDecade = 20;
52 const double kGFlk3DEvMin = 0.100; // GeV
53 
54 class GFLUKAAtmoFlux: public GAtmoFlux {
55 
56 public :
59 
60  //
61  // Most implementation is derived from the base GAtmoFlux
62  // The concrete driver is only required to implement a function for
63  // loading the input data files
64  //
65 
66 private:
67 
68  void SetBinSizes (void);
69  bool FillFluxHisto (int nu_pdg, string filename);
70 };
71 
72 } // flux namespace
73 } // genie namespace
74 
75 #endif // _GFLUKA_ATMO_3D_FLUX_I_H_
THE MAIN GENIE PROJECT NAMESPACE
Definition: GeneratorBase.h:8
const double kGFlk3DCosThetaMax
string filename
Definition: shutoffs.py:106
Loaders::FluxType flux
const unsigned int kGFlk3DNumCosThetaBins
const unsigned int kGFlk3DNumLogEvBins
const unsigned int kGFlk3DNumLogEvBinsPerDecade
bool FillFluxHisto(int nu_pdg, string filename)
A flux driver for the FLUKA 3-D Atmospheric Neutrino Flux.
const double kGFlk3DCosThetaMin
const double kGFlk3DEvMin
A base class for the FLUKA, BGLRS and ATMNC atmo. nu. flux drivers. The driver depends on data files ...
Definition: GAtmoFlux.h:61