genie::flux Namespace Reference

GENIE flux drivers. More...

Classes
class	GAstroFlux
	A base class for the concrete astrophysical neutrino flux drivers. More...
class	GAtmoFlux
	A base class for the FLUKA, BGLRS and ATMNC atmo. nu. flux drivers. The driver depends on data files provided by the atmospheric neutrino flux simulation authors in order to determine the angular and energy dependence for each neutrino species. The position of each flux neutrino [going towards a detector centered at (0,0,0)] is generated uniformly on a plane that is perpendicular to a sphere of radius Rl at the point that is determined by the generated neutrino direction (theta,phi). The size of the area of that plane, where flux neutrinos are generated, is determined by the transverse radius Rt. You can tweak Rl, Rt to match the size of your detector. Initially, neutrino coordinates are generated in a default detector coordinate system (Topocentric Horizontal Coordinate -THZ-): +z: Points towards the local zenith. +x: On same plane as local meridian, pointing south. +y: As needed to make a right-handed coordinate system. origin: detector centre Alternative user-defined topocentric systems can be defined by specifying the appropriate rotation from THZ. The driver allows minimum and maximum energy cuts. Also it provides the options to generate wither unweighted or weighted flux neutrinos (the latter giving smoother distributions at the tails). More...
class	GBGLRSAtmoFlux
	A flux driver for the Bartol Atmospheric Neutrino Flux. More...
class	GCylindTH1Flux
	A generic GENIE flux driver. Generates a 'cylindrical' neutrino beam along the input direction, with the input transverse radius and centered at the input position. The energies are generated from the input energy spectrum (TH1D). Multiple neutrino species can be generated (you will need to supply an energy spectrum for each). More...
class	GDiffuseAstroFlux
class	GFlavorMap
	GENIE interface for flavor modification. More...
class	GFlavorMixerFactory
	A class for generating concrete GFlavorMixerI derived classes based on the factory pattern. This code supplies a CPP macro which allows the classes to self-register and thus no modification of this class is needed in order to expand the list of classes it knows about. More...
class	GFlavorMixerI
	GENIE interface for flavor modification. More...
class	GFLUKAAtmoFlux
	A flux driver for the FLUKA 3-D Atmospheric Neutrino Flux. More...
class	GFluxBlender
	GENIE GFluxI adapter to allow flavor modification. More...
class	GFluxDriverFactory
	A class for generating concrete GFluxI derived classes based on the factory pattern. This code supplies a CPP macro which allows the classes to self-register and thus no modification of this class is needed in order to expand the list of classes it knows about. More...
class	GFluxExposureI
	GENIE interface for uniform flux exposure iterface. More...
class	GFluxFileConfigI
class	GHAKKMAtmoFlux
	A driver for the HAKKM 3-D atmospheric neutrino flux (commonly known as the `Honda flux') More...
class	GJPARCNuFlux
	A GENIE flux driver encapsulating the JPARC neutrino flux. It reads-in the official JPARC neutrino flux ntuples. More...
class	GJPARCNuFluxPassThroughInfo
class	GMonoEnergeticFlux
	A simple GENIE flux driver for monoenergetic neutrinos along the z direction. Can handle a mix of neutrinos with their corresponding weight. This trivial flux driver was added in so that single energy neutrinos can be easily used with the event generation driver that can handle a target mix or detailed geometries. More...
class	GNuMIFlux
	A GENIE flux driver encapsulating the NuMI neutrino flux. It reads-in the official GNUMI neutrino flux ntuples. Supports both geant3 and geant4 formats. More...
class	GNuMIFluxPassThroughInfo
class	GNuMIFluxXMLHelper
class	GPointSourceAstroFlux
class	GSimpleNtpAux
class	GSimpleNtpEntry
class	GSimpleNtpFlux
	A GENIE flux driver using a simple ntuple format. More...
class	GSimpleNtpMeta
class	GSimpleNtpNuMI

Typedefs
typedef genie::flux::GFlavorMixerI *(*	GFlavorMixerICtorFuncPtr_t) ()
typedef genie::GFluxI *(*	GFluxICtorFuncPtr_t) ()
typedef enum genie::flux::EExposure	Exposure_t

Enumerations
enum	EExposure { kUnknown = 0, kPOTs = 1, kSeconds = 2, kNExposureTypes }

Functions
ostream &	operator<< (ostream &stream, const genie::flux::GJPARCNuFluxPassThroughInfo &info)
ostream &	operator<< (ostream &stream, const genie::flux::GNuMIFluxPassThroughInfo &info)
ostream &	operator<< (ostream &stream, const genie::flux::GSimpleNtpEntry &entry)
ostream &	operator<< (ostream &stream, const genie::flux::GSimpleNtpNuMI &numi)
ostream &	operator<< (ostream &stream, const genie::flux::GSimpleNtpAux &aux)
ostream &	operator<< (ostream &stream, const genie::flux::GSimpleNtpMeta &meta)

Variables
const double	kAstroDefMaxEv = 1E+20 * units::GeV
const double	kAstroDefMinEv = 1E-3 * units::GeV
const int	kAstroNlog10EvBins = 1000
const int	kAstroNCosThetaBins = 500
const int	kAstroNPhiBins = 500
const unsigned int	kBGLRS3DNumCosThetaBins = 20
const double	kBGLRS3DCosThetaMin = -1.0
const double	kBGLRS3DCosThetaMax = 1.0
const unsigned int	kBGLRS3DNumLogEvBinsLow = 40
const unsigned int	kBGLRS3DNumLogEvBinsPerDecadeLow = 20
const unsigned int	kBGLRS3DNumLogEvBinsHigh = 30
const unsigned int	kBGLRS3DNumLogEvBinsPerDecadeHigh = 10
const double	kBGLRS3DEvMin = 0.1
const unsigned int	kGFlk3DNumCosThetaBins = 40
const double	kGFlk3DCosThetaMin = -1.0
const double	kGFlk3DCosThetaMax = 1.0
const unsigned int	kGFlk3DNumLogEvBins = 61
const unsigned int	kGFlk3DNumLogEvBinsPerDecade = 20
const double	kGFlk3DEvMin = 0.100
const unsigned int	kGHnd3DNumCosThetaBins = 20
const double	kGHnd3DCosThetaMin = -1.0
const double	kGHnd3DCosThetaMax = 1.0
const unsigned int	kGHnd3DNumPhiBins = 12
const double	kGHnd3DPhiMin = 0.0
const double	kGHnd3DPhiMax = 360.0
const unsigned int	kGHnd3DNumLogEvBins = 101
const unsigned int	kGHnd3DNumLogEvBinsPerDecade = 20
const double	kGHnd3DEvMin = 0.1
const int	fNgmax = 12

Detailed Description

GENIE flux drivers.

Typedef Documentation

typedef enum genie::flux::EExposure genie::flux::Exposure_t

typedef genie::flux::GFlavorMixerI*(* genie::flux::GFlavorMixerICtorFuncPtr_t) ()

Definition at line 41 of file GFlavorMixerFactory.h.

typedef genie::GFluxI*(* genie::flux::GFluxICtorFuncPtr_t) ()

Definition at line 44 of file GFluxDriverFactory.h.

Enumeration Type Documentation

enum genie::flux::EExposure

Enumerator
kUnknown
kPOTs
kSeconds
kNExposureTypes

Definition at line 30 of file GFluxExposureI.h.

                         {
    kUnknown   = 0,
    kPOTs      = 1,  // particles (protons) on target
    kSeconds   = 2,  // exposure duration
    kNExposureTypes  //
  } Exposure_t;

Function Documentation

ostream & genie::flux::operator<<	(	ostream &	stream,
		const genie::flux::GSimpleNtpEntry &	entry
	)

Definition at line 967 of file GSimpleNtpFlux.cxx.

    {
      stream << "\nGSimpleNtpEntry "
             << " PDG " << entry.pdg
             << " wgt " << entry.wgt
             << " ( metakey " << entry.metakey << " )"
             << "\n   vtx [" << entry.vtxx << "," << entry.vtxy << ","
             << entry.vtxz << ", t=" << entry.vtxt << "] dist " << entry.dist
             << "\n   p4  [" << entry.px << "," << entry.py << ","
             << entry.pz << "," << entry.E << "]";
     return stream;
  }

ostream & genie::flux::operator<<	(	ostream &	stream,
		const genie::flux::GSimpleNtpNuMI &	numi
	)

Definition at line 982 of file GSimpleNtpFlux.cxx.

{
  stream << "\nGSimpleNtpNuMI "
         << "run " << numi.run
         << " evtno " << numi.evtno
         << " entryno " << numi.entryno
         << "\n   ndecay " << numi.ndecay  << " ptype " << numi.ptype
         << "\n   tptype " << numi.tptype << " ppmedium " << numi.ppmedium
         << "\n  tp[xyz] [" << numi.tpx  << "," << numi.tpy  << "," << numi.tpz  << "]"
         << "\n   v[xyz] [" << numi.vx   << "," << numi.vy   << "," << numi.vz   << "]"
         << "\n  pd[xyz] [" << numi.pdpx << "," << numi.pdpy << "," << numi.pdpz << "]"
         << "\n  pp[xyz] [" << numi.pppx << "," << numi.pppy << "," << numi.pppz << "]"
    ;
  return stream;
}

ostream & genie::flux::operator<<	(	ostream &	stream,
		const genie::flux::GSimpleNtpAux &	aux
	)

Definition at line 999 of file GSimpleNtpFlux.cxx.

    {
      stream << "\nGSimpleNtpAux ";
      size_t nInt = aux.auxint.size();
      stream << "\n   ints: ";
      for (size_t ijInt=0; ijInt < nInt; ++ijInt)
        stream << " " << aux.auxint[ijInt];
      size_t nDbl = aux.auxdbl.size();
      stream << "\n   doubles: ";
      for (size_t ijDbl=0; ijDbl < nDbl; ++ijDbl)
        stream << " " << aux.auxdbl[ijDbl];

     return stream;
  }

ostream & genie::flux::operator<<	(	ostream &	stream,
		const genie::flux::GSimpleNtpMeta &	meta
	)

Definition at line 1015 of file GSimpleNtpFlux.cxx.

    {
      size_t nf = meta.pdglist.size();
      stream << "\nGSimpleNtpMeta " << nf << " flavors: ";
      for (size_t i=0; i<nf; ++i) stream << " " << meta.pdglist[i];

      //stream << "\nGSimpleNtpMeta " << meta.nflavors
      //       << " flavors: ";
      //for (int i=0; i< meta.nflavors; ++i) stream << " " << meta.flavor[i];

      stream << "\n maxEnergy " << meta.maxEnergy
             << " min/maxWgt " << meta.minWgt << "/" << meta.maxWgt
             << " protons " << meta.protons
             << " metakey " << meta.metakey
             << "\n windowBase [" << meta.windowBase[0] << ","
             << meta.windowBase[1] << "," << meta.windowBase[2] << "]"
             << "\n windowDir1 [" << meta.windowDir1[0] << ","
             << meta.windowDir1[1] << "," << meta.windowDir1[2] << "]"
             << "\n windowDir2 [" << meta.windowDir2[0] << ","
             << meta.windowDir2[1] << "," << meta.windowDir2[2] << "]";

      size_t nInt = meta.auxintname.size();
      if ( nInt > 0 ) stream << "\n aux ints:    ";
      for (size_t ijInt=0; ijInt < nInt; ++ijInt)
        stream << " " << meta.auxintname[ijInt];

      size_t nDbl = meta.auxdblname.size();
      if ( nDbl > 0 ) stream << "\n aux doubles: ";
      for (size_t ijDbl=0; ijDbl < nDbl; ++ijDbl)
        stream << " " << meta.auxdblname[ijDbl];

      size_t nfiles = meta.infiles.size();
      stream << "\n " << nfiles << " input files: ";
      UInt_t nprint = TMath::Min(UInt_t(nfiles),
                                 genie::flux::GSimpleNtpMeta::mxfileprint);
      for (UInt_t ifiles=0; ifiles < nprint; ++ifiles)
        stream << "\n    " << meta.infiles[ifiles];

      stream << "\n input seed: " << meta.seed;

     return stream;
  }

ostream & genie::flux::operator<<	(	ostream &	stream,
		const genie::flux::GJPARCNuFluxPassThroughInfo &	info
	)

Definition at line 1078 of file GJPARCNuFlux.cxx.

    {
      stream << "\n idfd       = " << info.idfd
             << "\n norm       = " << info.norm
             << "\n flux entry = " << info.fluxentry
             << "\n flux file  = " << info.fluxfilename
             << "\n Enu        = " << info.Enu
             << "\n geant code = " << info.ppid
             << "\n (pdg code) = " << pdg::GeantToPdg(info.ppid)
             << "\n decay mode = " << info.mode
             << "\n nvtx0      = " << info.nvtx0
             << "\n |momentum| @ decay       = " << info.ppi
             << "\n position_vector @ decay  = ("
                 << info.xpi[0] << ", "
                 << info.xpi[1] << ", "
                 << info.xpi[2]  << ")"
             << "\n direction_vector @ decay = ("
                 << info.npi[0] << ", "
                 << info.npi[1] << ", "
                 << info.npi[2]  << ")"
             << "\n |momentum| @ prod.       = " << info.ppi0
             << "\n position_vector @ prod.  = ("
                 << info.xpi0[0] << ", "
                 << info.xpi0[1] << ", "
                 << info.xpi0[2]  << ")"
             << "\n direction_vector @ prod. = ("
                 << info.npi0[0] << ", "
                 << info.npi0[1] << ", "
                 << info.npi0[2]  << ")"
             << "\n cospibm = " << info.cospibm
             << "\n cospi0bm = " << info.cospi0bm
             << "\n Plus additional info if flux version is later than 07a"
             << endl;

    return stream;
  }

ostream & genie::flux::operator<<	(	ostream &	stream,
		const genie::flux::GNuMIFluxPassThroughInfo &	info
	)

Definition at line 2053 of file GNuMIFlux.cxx.

    {
      // stream << "\n ndecay   = " << info.ndecay << std::endl;
      stream << "\nGNuMIFlux run " << info.run << " evtno " << info.evtno
             << " (pcodes " << info.pcodes << " units " << info.units << ")"
             << "\n random dk: dx/dz " << info.ndxdz
             << " dy/dz " << info.ndydz
             << " pz " <<  info.npz << " E " << info.nenergy
             << "\n near00 dk: dx/dz " << info.ndxdznea
             << " dy/dz " << info.ndydznea
             << " E " <<  info.nenergyn << " wgt " << info.nwtnear
             << "\n far00  dk: dx/dz " << info.ndxdzfar
             << " dy/dz " << info.ndydzfar
             << " E " <<  info.nenergyf << " wgt " << info.nwtfar
             << "\n norig " << info.norig << " ndecay " << info.ndecay
             << " ntype " << info.ntype
             << "\n had vtx " << info.vx << " " << info.vy << " " << info.vz
             << "\n parent p3 @ dk " << info.pdpx << " " << info.pdpy << " " << info.pdpz
             << "\n parent prod: dx/dz " << info.ppdxdz
             << " dy/dz " << info.ppdydz
             << " pz " << info.pppz << " E " << info.ppenergy
             << "\n ppmedium " << info.ppmedium << " ptype " << info.ptype
             << " ppvtx " << info.ppvx << " " << info.ppvy << " " << info.ppvz
             << "\n mu parent p4 " << info.muparpx << " " << info.muparpy
             << " " << info.muparpz << " " << info.mupare
             << "\n necm " << info.necm << " nimpwt " << info.nimpwt
             << "\n point x,y,z " << info.xpoint << " " << info.ypoint
             << " " << info.zpoint
             << "\n tv x,y,z " << info.tvx << " " << info.tvy << " " << info.tvz
             << "\n tptype " << info.tptype << " tgen " << info.tgen
             << " tgptype " << info.tgptype
             << "\n tgp px,py,pz " << info.tgppx << " " << info.tgppy
             << " " << info.tgppz
             << "\n tpriv x,y,z " << info.tprivx << " " << info.tprivy
             << " " << info.tprivz
             << "\n beam x,y,z " << info.beamx << " " << info.beamy
             << " " << info.beamz
             << "\n beam px,py,pz " << info.beampx << " " << info.beampy
             << " " << info.beampz
        ;

#ifndef SKIP_MINERVA_MODS
      //=========================================
      // The following was inserted by MINERvA
      //=========================================
      stream << "\nNeutrino History : ntrajectories " << info.ntrajectory
             << "\n (trkID, pdg) of nu parents: [";
      Int_t ntraj = info.ntrajectory;
      if ( info.overflow ) ntraj = GNuMIFluxPassThroughInfo::MAX_N_TRAJ;

      for ( Int_t itt = 0; itt < ntraj; ++itt )
        stream << "(" << info.trackId[itt-1] << ", " << info.pdgcode[itt] << ") ";
      stream << "]\n";
      //END of minerva additions
#endif

      stream << "\nCurrent: pdg " << info.fgPdgC
             << " xywgt " << info.fgXYWgt
             << "\n p4 (beam): " << utils::print::P4AsShortString(&info.fgP4)
             << "\n x4 (beam): " << utils::print::X4AsString(&info.fgX4)
             << "\n p4 (user): " << utils::print::P4AsShortString(&info.fgP4User)
             << "\n x4 (user): " << utils::print::X4AsString(&info.fgX4User);
        ;
#ifdef GNUMI_TEST_XY_WGT
      stream << "\n" << xypartials::GetStaticInstance();
#endif


  /*
  //std::cout << "GNuMIFlux::PrintCurrent ....." << std::endl;
  //LOG("Flux", pINFO)
  LOG("Flux", pNOTICE)
    << "Current Leaf Values: "
    << " run " << fLf_run << " evtno " << fLf_evtno << "\n"
    << " NenergyN " << fLf_NenergyN << " NWtNear " << fLf_NWtNear
    << " NenergyF " << fLf_NenergyF << " NWtFar  " << fLf_NWtFar << "\n"
    << " Norig " << fLf_Norig << " Ndecay " << fLf_Ndecay << " Ntype " << fLf_Ntype << "\n"
    << " Vxyz " << fLf_Vx << " " << fLf_Vy << " " << fLf_Vz
    << " pdPxyz " << fLf_pdPx << " " << fLf_pdPy << " " << fLf_pdPz << "\n"
    << " pp dxdz " << fLf_ppdxdz << " dydz " << fLf_ppdydz << " pz " << fLf_pppz << "\n"
    << " pp energy " << fLf_ppenergy << " medium " << fLf_ppmedium
    << " ptype " << fLf_ptype
    << " ppvxyz " << fLf_ppvx << " " << fLf_ppvy << " " << fLf_ppvz << "\n"
    << " muparpxyz " << fLf_muparpx << " " << fLf_muparpy << " " << fLf_muparpz
    << " mupare " << fLf_mupare << "\n"
    << ;
  */

     return stream;
  }

Variable Documentation

const int genie::flux::fNgmax = 12

Definition at line 151 of file GJPARCNuFlux.h.

const double genie::flux::kAstroDefMaxEv = 1E+20 * units::GeV

Definition at line 116 of file GAstroFlux.h.

const double genie::flux::kAstroDefMinEv = 1E-3 * units::GeV

Definition at line 117 of file GAstroFlux.h.

const int genie::flux::kAstroNCosThetaBins = 500

Definition at line 119 of file GAstroFlux.h.

const int genie::flux::kAstroNlog10EvBins = 1000

Definition at line 118 of file GAstroFlux.h.

const int genie::flux::kAstroNPhiBins = 500

Definition at line 120 of file GAstroFlux.h.

const double genie::flux::kBGLRS3DCosThetaMax = 1.0

Definition at line 52 of file GBGLRSAtmoFlux.h.

const double genie::flux::kBGLRS3DCosThetaMin = -1.0

Definition at line 51 of file GBGLRSAtmoFlux.h.

const double genie::flux::kBGLRS3DEvMin = 0.1

Definition at line 57 of file GBGLRSAtmoFlux.h.

const unsigned int genie::flux::kBGLRS3DNumCosThetaBins = 20

Definition at line 50 of file GBGLRSAtmoFlux.h.

const unsigned int genie::flux::kBGLRS3DNumLogEvBinsHigh = 30

Definition at line 55 of file GBGLRSAtmoFlux.h.

const unsigned int genie::flux::kBGLRS3DNumLogEvBinsLow = 40

Definition at line 53 of file GBGLRSAtmoFlux.h.

const unsigned int genie::flux::kBGLRS3DNumLogEvBinsPerDecadeHigh = 10

Definition at line 56 of file GBGLRSAtmoFlux.h.

const unsigned int genie::flux::kBGLRS3DNumLogEvBinsPerDecadeLow = 20

Definition at line 54 of file GBGLRSAtmoFlux.h.

const double genie::flux::kGFlk3DCosThetaMax = 1.0

Definition at line 49 of file GFLUKAAtmoFlux.h.

const double genie::flux::kGFlk3DCosThetaMin = -1.0

Definition at line 48 of file GFLUKAAtmoFlux.h.

const double genie::flux::kGFlk3DEvMin = 0.100

Definition at line 52 of file GFLUKAAtmoFlux.h.

const unsigned int genie::flux::kGFlk3DNumCosThetaBins = 40

Definition at line 47 of file GFLUKAAtmoFlux.h.

const unsigned int genie::flux::kGFlk3DNumLogEvBins = 61

Definition at line 50 of file GFLUKAAtmoFlux.h.

const unsigned int genie::flux::kGFlk3DNumLogEvBinsPerDecade = 20

Definition at line 51 of file GFLUKAAtmoFlux.h.

const double genie::flux::kGHnd3DCosThetaMax = 1.0

Definition at line 37 of file GHAKKMAtmoFlux.h.

const double genie::flux::kGHnd3DCosThetaMin = -1.0

Definition at line 36 of file GHAKKMAtmoFlux.h.

const double genie::flux::kGHnd3DEvMin = 0.1

Definition at line 43 of file GHAKKMAtmoFlux.h.

const unsigned int genie::flux::kGHnd3DNumCosThetaBins = 20

Definition at line 35 of file GHAKKMAtmoFlux.h.

const unsigned int genie::flux::kGHnd3DNumLogEvBins = 101

Definition at line 41 of file GHAKKMAtmoFlux.h.

const unsigned int genie::flux::kGHnd3DNumLogEvBinsPerDecade = 20

Definition at line 42 of file GHAKKMAtmoFlux.h.

const unsigned int genie::flux::kGHnd3DNumPhiBins = 12

Definition at line 38 of file GHAKKMAtmoFlux.h.

const double genie::flux::kGHnd3DPhiMax = 360.0

Definition at line 40 of file GHAKKMAtmoFlux.h.

const double genie::flux::kGHnd3DPhiMin = 0.0

Definition at line 39 of file GHAKKMAtmoFlux.h.