simb Namespace Reference

Base utilities and modules for event generation and detector simulation. More...

Namespaces
	gtruthaux

Classes
class	GTruth
class	MCFlux
struct	MCGeneratorInfo
class	MCNeutrino
	Event generator information. More...
class	MCParticle
class	MCTrajectory
class	MCTruth
	Event generator information. More...

Typedefs
typedef enum simb::_ev_generator	Generator_t
	generator used to produce event, if applicable More...
typedef enum simb::_ev_origin	Origin_t
	event origin types More...

Enumerations
enum	flux_code_ {   kHistPlusFocus = +1, kHistMinusFocus = -1, kGenerator = 0, kNtuple = 2,   kSimple_Flux = 3, kDk2Nu = 4 }
	Which flux was used to generate this event? More...
enum	_ev_generator {   kUnknown, _ev_generator::kGENIE, _ev_generator::kCRY, _ev_generator::kGIBUU,   _ev_generator::kNuWro, _ev_generator::kMARLEY, _ev_generator::kNEUT, _ev_generator::kCORSIKA,   _ev_generator::kGEANT, _ev_generator::kNumGenerators }
	generator used to produce event, if applicable More...
enum	curr_type_ { kCC, kNC }
	Neutrino interaction categories. More...
enum	int_type_ {   kUnknownInteraction = -1, kQE = 0, kRes = 1, kDIS = 2,   kCoh = 3, kCohElastic = 4, kElectronScattering = 5, kIMDAnnihilation = 6,   kInverseBetaDecay = 7, kGlashowResonance = 8, kAMNuGamma = 9, kMEC = 10,   kDiffractive = 11, kEM = 12, kWeakMix = 13, kNuanceOffset = 1000,   kCCQE = kNuanceOffset + 1, kNCQE = kNuanceOffset + 2, kResCCNuProtonPiPlus = kNuanceOffset + 3, kResCCNuNeutronPi0 = kNuanceOffset + 4,   kResCCNuNeutronPiPlus = kNuanceOffset + 5, kResNCNuProtonPi0 = kNuanceOffset + 6, kResNCNuProtonPiPlus = kNuanceOffset + 7, kResNCNuNeutronPi0 = kNuanceOffset + 8,   kResNCNuNeutronPiMinus = kNuanceOffset + 9, kResCCNuBarNeutronPiMinus = kNuanceOffset + 10, kResCCNuBarProtonPi0 = kNuanceOffset + 11, kResCCNuBarProtonPiMinus = kNuanceOffset + 12,   kResNCNuBarProtonPi0 = kNuanceOffset + 13, kResNCNuBarProtonPiPlus = kNuanceOffset + 14, kResNCNuBarNeutronPi0 = kNuanceOffset + 15, kResNCNuBarNeutronPiMinus = kNuanceOffset + 16,   kResCCNuDeltaPlusPiPlus = kNuanceOffset + 17, kResCCNuDelta2PlusPiMinus = kNuanceOffset + 21, kResCCNuBarDelta0PiMinus = kNuanceOffset + 28, kResCCNuBarDeltaMinusPiPlus = kNuanceOffset + 32,   kResCCNuProtonRhoPlus = kNuanceOffset + 39, kResCCNuNeutronRhoPlus = kNuanceOffset + 41, kResCCNuBarNeutronRhoMinus = kNuanceOffset + 46, kResCCNuBarNeutronRho0 = kNuanceOffset + 48,   kResCCNuSigmaPlusKaonPlus = kNuanceOffset + 53, kResCCNuSigmaPlusKaon0 = kNuanceOffset + 55, kResCCNuBarSigmaMinusKaon0 = kNuanceOffset + 60, kResCCNuBarSigma0Kaon0 = kNuanceOffset + 62,   kResCCNuProtonEta = kNuanceOffset + 67, kResCCNuBarNeutronEta = kNuanceOffset + 70, kResCCNuKaonPlusLambda0 = kNuanceOffset + 73, kResCCNuBarKaon0Lambda0 = kNuanceOffset + 76,   kResCCNuProtonPiPlusPiMinus = kNuanceOffset + 79, kResCCNuProtonPi0Pi0 = kNuanceOffset + 80, kResCCNuBarNeutronPiPlusPiMinus = kNuanceOffset + 85, kResCCNuBarNeutronPi0Pi0 = kNuanceOffset + 86,   kResCCNuBarProtonPi0Pi0 = kNuanceOffset + 90, kCCDIS = kNuanceOffset + 91, kNCDIS = kNuanceOffset + 92, kUnUsed1 = kNuanceOffset + 93,   kUnUsed2 = kNuanceOffset + 94, kCCQEHyperon = kNuanceOffset + 95, kNCCOH = kNuanceOffset + 96, kCCCOH = kNuanceOffset + 97,   kNuElectronElastic = kNuanceOffset + 98, kInverseMuDecay = kNuanceOffset + 99, kMEC2p2h = kNuanceOffset + 100 }
	Neutrino interaction categories. More...
enum	_ev_origin {   kUnknown, kUnknown, kBeamNeutrino, kCosmicRay,   kSuperNovaNeutrino, kSingleParticle }
	event origin types More...

Functions
std::ostream &	operator<< (std::ostream &output, const simb::GTruth &gtruth)
std::ostream &	operator<< (std::ostream &output, const simb::MCFlux &mcflux)
std::ostream &	operator<< (std::ostream &output, const simb::MCNeutrino &mcnu)
std::ostream &	operator<< (std::ostream &output, const MCParticle &particle)
std::ostream &	operator<< (std::ostream &output, const MCTrajectory &list)
std::ostream &	operator<< (std::ostream &o, simb::MCTruth const &a)

Generated particle index
This type represents the position of a specific generated particle in the generator truth record, simb::MCTruth. The particle itself can be accessed from that record via simb::MCTruth::GetParticle().
using	GeneratedParticleIndex_t = std::size_t
	Type of particle index in the generator truth record (simb::MCTruth). More...
constexpr GeneratedParticleIndex_t	NoGeneratedParticleIndex = std::numeric_limits<GeneratedParticleIndex_t>::max()
	Constant representing the absence of generator truth information. More...
bool	isGeneratedParticleIndex (GeneratedParticleIndex_t index)
	Returns whether the specified one is an acceptable generator index. More...

Detailed Description

Base utilities and modules for event generation and detector simulation.

This class stores/retrieves the additional information needed (and not in MCTruth) to recreate a genie::EventRecord for genie based event reweighting.

This class describes a particle created in the detector Monte Carlo simulation.

This class describes the trajectory of a particle created in the Monte Carlo simulation. It generally behaves like a vector< pair<TLorentzVector,TLorentzVector> >, where the first TLorentzVector is the position and the seoond is the momentum, with the following additions: - Methods Position(int) and Momentum(int) for those who are unfamiliar with the concept of "first" and "second" as used with STL pairs: sim::Trajectory* trajectory = simb::MCParticle.Trajectory(); int numberOfPonts = trajectory->size(); for (int i=0; i<numberOfPoints; ++i) { TLorentzVector position = trajectory->Position(i); TLorentzVector momentum = trajectory->Momentum(i); } The STL equivalent to the above statements (more efficient): sim::Trajectory* trajectory = simb::MCParticle.Trajectory(); for ( sim::Trajectory::const_iterator i = trajectory->begin(); i != trajectory->end(); ++i ) { const TLorentzVector& position = (*i).first; const TLorentzVector& momentum = (*i).second; } - As above, but for each position or momentum component; e.g., trajectory->X(i). - In addition to push_back(pair< TLorentzVector, TLorentzVector>), there's also push_back(TLorentzVector,TLorentzVector) and Add(TLorentzVector,TLorentzVector). They all do the same thing: add another point to the trajectory. - Print() and operator<< methods for ROOT display and ease of debugging. There are no units defined in this class. If it's used with Geant4, the units will be (mm,ns,GeV), but this class does not enforce this.

Typedef Documentation

using simb::GeneratedParticleIndex_t = typedef std::size_t

Type of particle index in the generator truth record (simb::MCTruth).

Definition at line 30 of file simb.h.

typedef enum simb::_ev_generator simb::Generator_t

generator used to produce event, if applicable

typedef enum simb::_ev_origin simb::Origin_t

event origin types

Enumeration Type Documentation

enum simb::_ev_generator

strong

generator used to produce event, if applicable

Enumerator
kUnknown
kGENIE
kCRY
kGIBUU
kNuWro
kMARLEY
kNEUT
kCORSIKA
kGEANT
kNumGenerators

Definition at line 18 of file MCGeneratorInfo.h.

  {
    kUnknown,
    kGENIE,
    kCRY,
    kGIBUU,
    kNuWro,
    kMARLEY,
    kNEUT,
    kCORSIKA,
    kGEANT,
    kNumGenerators, //  this should always be the last entry
  } Generator_t;

enum simb::_ev_origin

event origin types

Enumerator
kUnknown	???
kUnknown
kBeamNeutrino	Beam neutrinos.
kCosmicRay	Cosmic rays.
kSuperNovaNeutrino	Supernova neutrinos.
kSingleParticle	single particles thrown at the detector

Definition at line 21 of file MCTruth.h.

                         {
    kUnknown,           ///< ???
    kBeamNeutrino,      ///< Beam neutrinos
    kCosmicRay,         ///< Cosmic rays
    kSuperNovaNeutrino, ///< Supernova neutrinos
    kSingleParticle     ///< single particles thrown at the detector
  } Origin_t;

enum simb::curr_type_

Neutrino interaction categories.

Enumerator
kCC
kNC

Definition at line 73 of file MCNeutrino.h.

                 {
    kCC,
    kNC
  };

enum simb::flux_code_

Which flux was used to generate this event?

Enumerator
kHistPlusFocus	Flux for positive horn focus.
kHistMinusFocus	Flux for negative horn focus.
kGenerator	A bogus flux assumed by the generator.
kNtuple	Full flux simulation ntuple.
kSimple_Flux	A simplified flux ntuple for quick running.
kDk2Nu	Unified ntuple flux format (replaces 2)

Definition at line 17 of file MCFlux.h.

                 {
    kHistPlusFocus  = +1, ///< Flux for positive horn focus
    kHistMinusFocus = -1, ///< Flux for negative horn focus
    kGenerator      =  0, ///< A bogus flux assumed by the generator
    kNtuple         =  2, ///< Full flux simulation ntuple
    kSimple_Flux    =  3, ///< A simplified flux ntuple for quick running
    kDk2Nu          =  4  ///< Unified ntuple flux format (replaces 2)
  };

enum simb::int_type_

Neutrino interaction categories.

Enumerator
kUnknownInteraction	n.b.: this group is similar but not quite, entirely unlike GENIE ScatteringType convention
kQE
kRes
kDIS
kCoh
kCohElastic
kElectronScattering
kIMDAnnihilation
kInverseBetaDecay
kGlashowResonance
kAMNuGamma
kMEC
kDiffractive
kEM
kWeakMix
kNuanceOffset	offset to account for adding in Nuance codes to this enum
kCCQE	charged current quasi-elastic
kNCQE	neutral current quasi-elastic
kResCCNuProtonPiPlus	resonant charged current, nu p -> l- p pi+
kResCCNuNeutronPi0	resonant charged current, nu n -> l- p pi0
kResCCNuNeutronPiPlus	resonant charged current, nu n -> l- n pi+
kResNCNuProtonPi0	resonant neutral current, nu p -> nu p pi0
kResNCNuProtonPiPlus	resonant neutral current, nu p -> nu p pi+
kResNCNuNeutronPi0	resonant neutral current, nu n -> nu n pi0
kResNCNuNeutronPiMinus	resonant neutral current, nu n -> nu p pi-
kResCCNuBarNeutronPiMinus	resonant charged current, nubar n -> l+ n pi-
kResCCNuBarProtonPi0	resonant charged current, nubar p -> l+ n pi0
kResCCNuBarProtonPiMinus	resonant charged current, nubar p -> l+ p pi-
kResNCNuBarProtonPi0	resonant charged current, nubar p -> nubar p pi0
kResNCNuBarProtonPiPlus	resonant charged current, nubar p -> nubar n pi+
kResNCNuBarNeutronPi0	resonant charged current, nubar n -> nubar n pi0
kResNCNuBarNeutronPiMinus	resonant charged current, nubar n -> nubar p pi-
kResCCNuDeltaPlusPiPlus
kResCCNuDelta2PlusPiMinus
kResCCNuBarDelta0PiMinus
kResCCNuBarDeltaMinusPiPlus
kResCCNuProtonRhoPlus
kResCCNuNeutronRhoPlus
kResCCNuBarNeutronRhoMinus
kResCCNuBarNeutronRho0
kResCCNuSigmaPlusKaonPlus
kResCCNuSigmaPlusKaon0
kResCCNuBarSigmaMinusKaon0
kResCCNuBarSigma0Kaon0
kResCCNuProtonEta
kResCCNuBarNeutronEta
kResCCNuKaonPlusLambda0
kResCCNuBarKaon0Lambda0
kResCCNuProtonPiPlusPiMinus
kResCCNuProtonPi0Pi0
kResCCNuBarNeutronPiPlusPiMinus
kResCCNuBarNeutronPi0Pi0
kResCCNuBarProtonPi0Pi0
kCCDIS	charged current deep inelastic scatter
kNCDIS	charged current deep inelastic scatter
kUnUsed1
kUnUsed2
kCCQEHyperon
kNCCOH
kCCCOH	charged current coherent pion
kNuElectronElastic	neutrino electron elastic scatter
kInverseMuDecay	inverse muon decay
kMEC2p2h	extension of nuance encoding for MEC / 2p2h

Definition at line 79 of file MCNeutrino.h.

                {
    kUnknownInteraction        =   -1, ///< n.b.: this group is similar but not quite, entirely unlike GENIE ScatteringType convention
    kQE                        =    0,
    kRes                       =    1,
    kDIS                       =    2,
    kCoh                       =    3,
    kCohElastic                =    4,
    kElectronScattering        =    5,
    kIMDAnnihilation           =    6,
    kInverseBetaDecay          =    7,
    kGlashowResonance          =    8,
    kAMNuGamma                 =    9,
    kMEC                       =   10,
    kDiffractive               =   11,
    kEM                        =   12,
    kWeakMix                   =   13,
    kNuanceOffset              = 1000,                ///< offset to account for adding in Nuance codes to this enum
    kCCQE                      = kNuanceOffset +  1,  ///< charged current quasi-elastic
    kNCQE                      = kNuanceOffset +  2,  ///< neutral current quasi-elastic
    kResCCNuProtonPiPlus       = kNuanceOffset +  3,  ///< resonant charged current, nu p -> l- p pi+
    kResCCNuNeutronPi0         = kNuanceOffset +  4,  ///< resonant charged current, nu n -> l- p pi0
    kResCCNuNeutronPiPlus      = kNuanceOffset +  5,  ///< resonant charged current, nu n -> l- n pi+
    kResNCNuProtonPi0          = kNuanceOffset +  6,  ///< resonant neutral current, nu p -> nu p pi0
    kResNCNuProtonPiPlus       = kNuanceOffset +  7,  ///< resonant neutral current, nu p -> nu p pi+
    kResNCNuNeutronPi0         = kNuanceOffset +  8,  ///< resonant neutral current, nu n -> nu n pi0
    kResNCNuNeutronPiMinus     = kNuanceOffset +  9,  ///< resonant neutral current, nu n -> nu p pi-
    kResCCNuBarNeutronPiMinus  = kNuanceOffset + 10,  ///< resonant charged current, nubar n -> l+ n pi-
    kResCCNuBarProtonPi0       = kNuanceOffset + 11,  ///< resonant charged current, nubar p -> l+ n pi0
    kResCCNuBarProtonPiMinus   = kNuanceOffset + 12,  ///< resonant charged current, nubar p -> l+ p pi-
    kResNCNuBarProtonPi0       = kNuanceOffset + 13,  ///< resonant charged current, nubar p -> nubar p pi0
    kResNCNuBarProtonPiPlus    = kNuanceOffset + 14,  ///< resonant charged current, nubar p -> nubar n pi+
    kResNCNuBarNeutronPi0      = kNuanceOffset + 15,  ///< resonant charged current, nubar n -> nubar n pi0
    kResNCNuBarNeutronPiMinus  = kNuanceOffset + 16,  ///< resonant charged current, nubar n -> nubar p pi-
    kResCCNuDeltaPlusPiPlus    = kNuanceOffset + 17,
    kResCCNuDelta2PlusPiMinus  = kNuanceOffset + 21,
    kResCCNuBarDelta0PiMinus   = kNuanceOffset + 28,
    kResCCNuBarDeltaMinusPiPlus= kNuanceOffset + 32,
    kResCCNuProtonRhoPlus      = kNuanceOffset + 39,
    kResCCNuNeutronRhoPlus     = kNuanceOffset + 41,
    kResCCNuBarNeutronRhoMinus = kNuanceOffset + 46,
    kResCCNuBarNeutronRho0     = kNuanceOffset + 48,
    kResCCNuSigmaPlusKaonPlus  = kNuanceOffset + 53,
    kResCCNuSigmaPlusKaon0     = kNuanceOffset + 55,
    kResCCNuBarSigmaMinusKaon0 = kNuanceOffset + 60,
    kResCCNuBarSigma0Kaon0     = kNuanceOffset + 62,
    kResCCNuProtonEta          = kNuanceOffset + 67,
    kResCCNuBarNeutronEta      = kNuanceOffset + 70,
    kResCCNuKaonPlusLambda0    = kNuanceOffset + 73,
    kResCCNuBarKaon0Lambda0    = kNuanceOffset + 76,
    kResCCNuProtonPiPlusPiMinus= kNuanceOffset + 79,
    kResCCNuProtonPi0Pi0       = kNuanceOffset + 80,
    kResCCNuBarNeutronPiPlusPiMinus = kNuanceOffset + 85,
    kResCCNuBarNeutronPi0Pi0   = kNuanceOffset + 86,
    kResCCNuBarProtonPi0Pi0    = kNuanceOffset + 90,
    kCCDIS                     = kNuanceOffset + 91,  ///< charged current deep inelastic scatter
    kNCDIS                     = kNuanceOffset + 92,  ///< charged current deep inelastic scatter
    kUnUsed1                   = kNuanceOffset + 93,
    kUnUsed2                   = kNuanceOffset + 94,
    kCCQEHyperon               = kNuanceOffset + 95,
    kNCCOH                     = kNuanceOffset + 96,
    kCCCOH                     = kNuanceOffset + 97,  ///< charged current coherent pion
    kNuElectronElastic         = kNuanceOffset + 98,  ///< neutrino electron elastic scatter
    kInverseMuDecay            = kNuanceOffset + 99,  ///< inverse muon decay
    kMEC2p2h                   = kNuanceOffset + 100  ///< extension of nuance encoding for MEC / 2p2h
  };

Function Documentation

bool simb::isGeneratedParticleIndex ( GeneratedParticleIndex_t  index )

inline

Returns whether the specified one is an acceptable generator index.

Definition at line 37 of file simb.h.

    { return index != NoGeneratedParticleIndex; }

std::ostream& simb::operator<<	(	std::ostream &	output,
		const MCTrajectory &	list
	)

Definition at line 67 of file MCTrajectory.cxx.

  {
    // Determine a field width for the voxel number.
    MCTrajectory::size_type numberOfTrajectories = list.size();
    int numberOfDigits = (int) std::log10( (double) numberOfTrajectories ) + 1;

    // A simple header.
    output.width( numberOfDigits );
    output << "#" << ": < position (x,y,z,t), momentum (Px,Py,Pz,E) >" << std::endl; 

    // Write each trajectory point on a separate line.
    MCTrajectory::size_type nTrajectory = 0;
    for ( MCTrajectory::const_iterator trajectory = list.begin(); trajectory != list.end(); ++trajectory, ++nTrajectory ){
        output.width( numberOfDigits );
        output << nTrajectory << ": "
               << "< (" << (*trajectory).first.X() 
               << "," << (*trajectory).first.Y() 
               << "," << (*trajectory).first.Z() 
               << "," << (*trajectory).first.T() 
               << ") , (" << (*trajectory).second.Px() 
               << "," << (*trajectory).second.Py() 
               << "," << (*trajectory).second.Pz() 
               << "," << (*trajectory).second.E() 
               << ") >" << std::endl;
      }

    return output;
  }

std::ostream& simb::operator<<	(	std::ostream &	o,
		simb::MCTruth const &	a
	)

Definition at line 70 of file MCTruth.cxx.

  {
    if(a.Origin() == kCosmicRay) 
      o << "This is a cosmic ray event" << std::endl;
    else if(a.Origin() == kBeamNeutrino){
      o << "This is a beam neutrino event" << std::endl;
      o << a.GetNeutrino();
    }
    else if(a.Origin() == kSuperNovaNeutrino){ 
      o << "This is a supernova neutrino event" << std::endl;
      o << a.GetNeutrino();
    }  

    for (int i = 0; i < a.NParticles(); ++i)
      o << i << " " << a.GetParticle(i) << std::endl;

    return o;
  }

std::ostream& simb::operator<<	(	std::ostream &	output,
		const simb::MCNeutrino &	mcnu
	)

Definition at line 80 of file MCNeutrino.cxx.

  {
    output << " neutrino =         " << mcnu.Nu().PdgCode()    << std::endl
           << " neutrino energy =  " << mcnu.Nu().E()          << std::endl
           << " CCNC =             " << mcnu.CCNC()            << std::endl
           << " mode =             " << mcnu.Mode()            << std::endl
           << " interaction type = " << mcnu.InteractionType() << std::endl
           << " target =           " << mcnu.Target()          << std::endl 
           << " nucleon =          " << mcnu.HitNuc()          << std::endl 
           << " quark =            " << mcnu.HitQuark()        << std::endl 
           << " W =                " << mcnu.W()               << std::endl 
           << " X =                " << mcnu.X()               << std::endl 
           << " Y =                " << mcnu.Y()               << std::endl 
           << " Q^2 =              " << mcnu.QSqr()            << std::endl;

    return output;
  }

std::ostream& simb::operator<<	(	std::ostream &	output,
		const simb::GTruth &	gtruth
	)

Definition at line 89 of file GTruth.cxx.

  {
    output
      << "GTruth:" << std::endl
      //  123456789012
      << " Vertex      " << gtruthaux::stringifyTLorentzVector(gtruth.fVertex)  << std::endl
      << " weight      " << std::setw(11) << gtruth.fweight       << " "
      << "  prob       " << std::setw(11) << gtruth.fprobability  << std::endl
      << " Xsec        " << std::setw(11) << gtruth.fXsec         << " "
      << "  DiffXsec   " << std::setw(11) << gtruth.fDiffXsec     << " "
      << "  GPhaseSpace " << std::setw(9) << gtruth.fGPhaseSpace << std::endl

      << " probe       " << std::setw(11) << gtruth.fProbePDG     << std::endl
      << " ProbeP4     " << gtruthaux::stringifyTLorentzVector(gtruth.fProbeP4) << std::endl
      << " TgtP4       " << gtruthaux::stringifyTLorentzVector(gtruth.fTgtP4) << std::endl

      << " Z A PDG (nuc,qrk) " << std::setw(4)  << gtruth.ftgtZ         << " "
                         << std::setw(4)  << gtruth.ftgtA         << " "
                         << std::setw(11) << gtruth.ftgtPDG       << " ("
                         << std::setw(11) << gtruth.fHitNucPDG    << " "
                         << std::setw(11) << gtruth.fHitQrkPDG    << ") "
      << " IsSeaQuark  " << std::setw(4)  << ((gtruth.fIsSeaQuark)?"yes":"no") << std::endl
      << " HitNucP4    " << gtruthaux::stringifyTLorentzVector(gtruth.fHitNucP4) << " "
      << " HitNucPos   " << std::setw(11) << gtruth.fHitNucPos    << std::endl

      << " Gscatter    " << std::setw(11) << gtruth.fGscatter     << " "
      << "  Gint       " << std::setw(11) << gtruth.fGint         << std::endl

      << "  Q2 q2      " << std::setw(11) << gtruth.fgQ2          << " "
                         << std::setw(11) << gtruth.fgq2          << std::endl
      << "  W  T       " << std::setw(11) << gtruth.fgW           << " "
                         << std::setw(11) << gtruth.fgT           << std::endl
      << "  X  Y       " << std::setw(11) << gtruth.fgX           << " "
                         << std::setw(11) << gtruth.fgY           << std::endl
      << "  Wrun       " << std::setw(11) << gtruth.fgWrun        << std::endl
      << " FSlepton    " << gtruthaux::stringifyTLorentzVector(gtruth.fFSleptonP4) << std::endl
      << " FShadSyst   " << gtruthaux::stringifyTLorentzVector(gtruth.fFShadSystP4) << std::endl

      << "  IsCharm    " << std::setw(4)  << ((gtruth.fIsCharm)?"yes":"no")   << " "
      << "  CharmPDG   " << std::setw(6)  << gtruth.fCharmHadronPdg         << " "
      << "  IsStrange  " << std::setw(4)  << ((gtruth.fIsStrange)?"yes":"no") << " "
      << "  StrangePDG " << std::setw(6)  << gtruth.fStrangeHadronPdg         << std::endl
      << "  Np Nn      " << std::setw(4)  << gtruth.fNumProton    << " "
                         << std::setw(4)  << gtruth.fNumNeutron   << " "
      << "  Npi(0,+,-) " << std::setw(4)  << gtruth.fNumPi0       << " "
                         << std::setw(4)  << gtruth.fNumPiPlus    << " "
                         << std::setw(4)  << gtruth.fNumPiMinus   << std::endl
      << "  NSingleGammas " << std::setw(4)  << gtruth.fNumSingleGammas
      << "  NRho(0,+,-) " << std::setw(4)  << gtruth.fNumRho0       << " "
                          << std::setw(4)  << gtruth.fNumRhoPlus    << " "
                          << std::setw(4)  << gtruth.fNumRhoMinus   << std::endl
      << "  ResNum     " << std::setw(4)  << gtruth.fResNum       << " "
      << "  DecayMode  " << std::setw(4)  << gtruth.fDecayMode    << std::endl
      << "  FinalQuarkPdg " << std::setw(4) << gtruth.fFinalQuarkPdg << " "
      << "  FinalLeptonPdg " << std::setw(4) << gtruth.fFinalLeptonPdg

      << std::endl ;

    return output;
  }

std::ostream& simb::operator<<	(	std::ostream &	output,
		const MCParticle &	particle
	)

Definition at line 151 of file MCParticle.cxx.

  {
    output << "ID=" << particle.TrackId() << ", ";
    int pdg =  particle.PdgCode();

    // Try to translate the PDG code into text.
    const TDatabasePDG* databasePDG = TDatabasePDG::Instance();
    const TParticlePDG* definition = databasePDG->GetParticle( pdg );
    // Check that the particle is known to ROOT.  If not, this is
    // not a major error; Geant4 has an internal particle coding
    // scheme for nuclei that ROOT doesn't recognize.
    if ( definition != 0 ) output << definition->GetName();
    else output << "PDG=" << pdg;
    
    output << ", mass="      << particle.Mass()
           << ", Mother ID=" << particle.Mother()
           << ", Process="   << particle.Process()
           << ", Status="    << particle.StatusCode()
           << "\nthere are " << particle.NumberTrajectoryPoints() << " trajectory points";

    if(particle.NumberTrajectoryPoints() > 0 )
      output << "\nInitial vtx (x,y,z,t)=(" << particle.Vx()
             << "," << particle.Vy()
             << "," << particle.Vz()
             << "," << particle.T()
             << "),\n Initial mom (Px,Py,Pz,E)=(" << particle.Px()
             << "," << particle.Py()
             << "," << particle.Pz()
             << "," << particle.E()
             << ")" << std::endl;
    else
      output << std::endl;

    return output;
  }

std::ostream& simb::operator<<	(	std::ostream &	output,
		const simb::MCFlux &	mcflux
	)

Definition at line 384 of file MCFlux.cxx.

  {
    output
      << "MCFlux:" << std::endl
      //  123456789012
      << "  flux job  " << std::setw(11) << mcflux.frun      << " "
      << "  pot #     " << std::setw(11) << mcflux.fevtno    << std::endl
      << "  ntype     " << std::setw(11) << mcflux.fntype    << " "
      << "  ptype     " << std::setw(11) << mcflux.fptype    << " "
      << "  impwt     " << std::setw(11) << mcflux.fnimpwt   << std::endl
      << "  ndecay    " << std::setw(11) << mcflux.fndecay   << " "
      << "  ppmedium  " << std::setw(11) << mcflux.fppmedium << " "
      << "  tptype    " << std::setw(11) << mcflux.ftptype   << std::endl
      << "  vxyz      " << std::setw(11) << mcflux.fvx       << " "
      <<                   std::setw(11) << mcflux.fvy       << " "
      <<                   std::setw(11) << mcflux.fvz       << std::endl
      << "  pdpxyz    " << std::setw(11) << mcflux.fpdpx     << " "
      <<                   std::setw(11) << mcflux.fpdpy     << " "
      <<                   std::setw(11) << mcflux.fpdpz     << std::endl
      << "  tpxyz     " << std::setw(11) << mcflux.ftpx      << " "
      <<                   std::setw(11) << mcflux.ftpy      << " "
      <<                   std::setw(11) << mcflux.ftpz      << std::endl
      << "  ppd[xy]dz " << std::setw(11) << mcflux.fppdxdz   << " "
      <<                   std::setw(11) << mcflux.fppdydz   << std::endl
      << "  pppxyz    " << std::setw(11) << mcflux.fppdxdz*mcflux.fpppz << " "
      <<                   std::setw(11) << mcflux.fppdydz*mcflux.fpppz << " "
      <<                   std::setw(11) << mcflux.fpppz     << std::endl
      << "  muparpxyz " << std::setw(11) << mcflux.fmuparpx  << " "
      <<                   std::setw(11) << mcflux.fmuparpy  << " "
      <<                   std::setw(11) << mcflux.fmuparpz  << " "
      << "  mupare    " << std::setw(11) << mcflux.fmupare   << std::endl
      << "  necm      " << std::setw(11) << mcflux.fnecm     << " "
      << "  dk2gen    " << std::setw(11) << mcflux.fdk2gen   << std::endl
      << "  near E    " << std::setw(11) << mcflux.fnenergyn << " "
      << "  wgt       " << std::setw(11) << mcflux.fnwtnear  << std::endl
      << "  far E     " << std::setw(11) << mcflux.fnenergyf << " "
      << "  wgt       " << std::setw(11) << mcflux.fnwtfar   << std::endl;
    return output;
  }

Variable Documentation

constexpr GeneratedParticleIndex_t simb::NoGeneratedParticleIndex = std::numeric_limits<GeneratedParticleIndex_t>::max()

Constant representing the absence of generator truth information.

Definition at line 34 of file simb.h.