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	GDk2NuFlux
	An implementation of the GENIE GFluxI interface ("flux driver") encapsulating reading/processing the "dk2nu" tree structure. More...

class	GDk2NuFluxXMLHelper

class	GFlavorMap
	GENIE interface for flavor modification. More...

class	GFlavorMixerFactory

class	GFlavorMixerI
	GENIE interface for flavor modification. More...

class	GFlavorOsc
	GENIE interface for flavor modification. More...

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

class	GFluxBlender

class	GFluxDriverFactory

class	GFluxExposureI
	GENIE interface for uniform flux exposure iterface. More...

class	GFluxFileConfigI
	GENIE interface for uniform flux exposure iterface. More...

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 33 of file GFlavorMixerFactory.h.

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

Definition at line 36 of file GFluxDriverFactory.h.

Enumeration Type Documentation

enum genie::flux::EExposure

Enumerator
kUnknown
kPOTs
kSeconds
kNExposureTypes

Definition at line 25 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 975 of file GSimpleNtpFlux.cxx.

References genie::flux::GSimpleNtpEntry::dist, genie::flux::GSimpleNtpEntry::E, genie::flux::GSimpleNtpEntry::metakey, genie::flux::GSimpleNtpEntry::pdg, genie::flux::GSimpleNtpEntry::px, genie::flux::GSimpleNtpEntry::py, genie::flux::GSimpleNtpEntry::pz, runNovaSAM::stream, genie::flux::GSimpleNtpEntry::vtxx, genie::flux::GSimpleNtpEntry::vtxy, genie::flux::GSimpleNtpEntry::vtxz, and genie::flux::GSimpleNtpEntry::wgt.

     {
       stream << "\nGSimpleNtpEntry " 
              << " PDG " << entry.pdg 
              << " wgt " << entry.wgt
              << " ( metakey " << entry.metakey << " )"
              << "\n   vtx [" << entry.vtxx << "," << entry.vtxy << ","
              << entry.vtxz << "] 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 990 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 1007 of file GSimpleNtpFlux.cxx.

References genie::flux::GSimpleNtpAux::auxdbl, and genie::flux::GSimpleNtpAux::auxint.

     {
       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 1023 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 1183 of file GJPARCNuFlux.cxx.

Referenced by genie::flux::GJPARCNuFluxPassThroughInfo::~GJPARCNuFluxPassThroughInfo().

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