simb Namespace Reference

ART objects. 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

ART objects.

Base utilities and modules for event generation and detector simulation.

Version

Id

PrimaryParticleInformation.h,v 1.7 2012-09-20 21:47:05 greenc Exp

Author

selig.nosp@m.man@.nosp@m.nevis.nosp@m..col.nosp@m.umbia.nosp@m..edu, brebe.nosp@m.l@fn.nosp@m.al.go.nosp@m.v

PrimaryParticleInformation 10-Sep-2007 Bill Seligman

11-Feb-2009 selig.nosp@m.man@.nosp@m.nevis.nosp@m..col.nosp@m.umbia.nosp@m..edu Revised for LArSoft

Purpose: This class is "attached" to the G4PrimaryParticle. It's used to save the MCTruth object associated with the event.

Background: Read this carefully, because this class probably doesn't do what you think it does.

Geant4 has various "truth" classes: G4Event, G4Track, G4PrimaryVertex, G4PrimaryParticle, etc. For all of these classes, Geant4 provides a facility for the user to include additional information that's "attached" to the class in question.

In this case, this class defines additional information to included with the G4PrimaryParticle class. In particular, it stores the pointer to the simb::MCTruth object that was the source of the G4PrimaryParticle information.

The reason why this class is necessary for the G4Base application is that it allows the ParticleListAction class access to the MCTruth pointer during Geant4's tracking.

Version

Id

ReweightAna.h,v 1.2 2012-09-19 18:12:09 nsmayer Exp

Author

nsmay.nosp@m.er@f.nosp@m.nal.g.nosp@m.ov

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,      
    kCosmicRay,         
    kSuperNovaNeutrino, 
    kSingleParticle     
  } 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, 
    kHistMinusFocus = -1, 
    kGenerator      =  0, 
    kNtuple         =  2, 
    kSimple_Flux    =  3, 
    kDk2Nu          =  4  
  };

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, 
    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  
  };

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.

References NoGeneratedParticleIndex.

Referenced by sim::GeneratedParticleInfo::hasGeneratedParticleIndex(), larg4::ParticleListAction::ParticleInfo_t::isPrimary(), and larg4::ParticleListActionService::ParticleInfo_t::isPrimary().

    { return index != NoGeneratedParticleIndex; }

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

Definition at line 67 of file MCTrajectory.cxx.

References simb::MCTrajectory::begin(), simb::MCTrajectory::end(), and simb::MCTrajectory::size().

  {
    // 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.

References simb::MCTruth::GetNeutrino(), simb::MCTruth::GetParticle(), kBeamNeutrino, kCosmicRay, kSuperNovaNeutrino, simb::MCTruth::NParticles(), and simb::MCTruth::Origin().

  {
    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.

References simb::MCNeutrino::CCNC(), simb::MCParticle::E(), simb::MCNeutrino::HitNuc(), simb::MCNeutrino::HitQuark(), simb::MCNeutrino::InteractionType(), simb::MCNeutrino::Mode(), simb::MCNeutrino::Nu(), simb::MCParticle::PdgCode(), simb::MCNeutrino::QSqr(), simb::MCNeutrino::Target(), simb::MCNeutrino::W(), simb::MCNeutrino::X(), and simb::MCNeutrino::Y().

  {
    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.

References simb::GTruth::fCharmHadronPdg, simb::GTruth::fDecayMode, simb::GTruth::fDiffXsec, simb::GTruth::fFinalLeptonPdg, simb::GTruth::fFinalQuarkPdg, simb::GTruth::fFShadSystP4, simb::GTruth::fFSleptonP4, simb::GTruth::fGint, simb::GTruth::fGPhaseSpace, simb::GTruth::fgQ2, simb::GTruth::fgq2, simb::GTruth::fGscatter, simb::GTruth::fgT, simb::GTruth::fgW, simb::GTruth::fgWrun, simb::GTruth::fgX, simb::GTruth::fgY, simb::GTruth::fHitNucP4, simb::GTruth::fHitNucPDG, simb::GTruth::fHitNucPos, simb::GTruth::fHitQrkPDG, simb::GTruth::fIsCharm, simb::GTruth::fIsSeaQuark, simb::GTruth::fIsStrange, simb::GTruth::fNumNeutron, simb::GTruth::fNumPi0, simb::GTruth::fNumPiMinus, simb::GTruth::fNumPiPlus, simb::GTruth::fNumProton, simb::GTruth::fNumRho0, simb::GTruth::fNumRhoMinus, simb::GTruth::fNumRhoPlus, simb::GTruth::fNumSingleGammas, simb::GTruth::fprobability, simb::GTruth::fProbeP4, simb::GTruth::fProbePDG, simb::GTruth::fResNum, simb::GTruth::fStrangeHadronPdg, simb::GTruth::ftgtA, simb::GTruth::fTgtP4, simb::GTruth::ftgtPDG, simb::GTruth::ftgtZ, simb::GTruth::fVertex, simb::GTruth::fweight, simb::GTruth::fXsec, and simb::gtruthaux::stringifyTLorentzVector().

  {
    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 157 of file MCParticle.cxx.

References simb::MCParticle::E(), simb::MCParticle::Mass(), simb::MCParticle::Mother(), simb::MCParticle::NumberTrajectoryPoints(), simb::MCParticle::PdgCode(), simb::MCParticle::Process(), simb::MCParticle::Px(), simb::MCParticle::Py(), simb::MCParticle::Pz(), simb::MCParticle::StatusCode(), simb::MCParticle::T(), simb::MCParticle::TrackId(), simb::MCParticle::Vx(), simb::MCParticle::Vy(), and simb::MCParticle::Vz().

  {
    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.

References simb::MCFlux::fdk2gen, simb::MCFlux::fevtno, simb::MCFlux::fmupare, simb::MCFlux::fmuparpx, simb::MCFlux::fmuparpy, simb::MCFlux::fmuparpz, simb::MCFlux::fndecay, simb::MCFlux::fnecm, simb::MCFlux::fnenergyf, simb::MCFlux::fnenergyn, simb::MCFlux::fnimpwt, simb::MCFlux::fntype, simb::MCFlux::fnwtfar, simb::MCFlux::fnwtnear, simb::MCFlux::fpdpx, simb::MCFlux::fpdpy, simb::MCFlux::fpdpz, simb::MCFlux::fppdxdz, simb::MCFlux::fppdydz, simb::MCFlux::fppmedium, simb::MCFlux::fpppz, simb::MCFlux::fptype, simb::MCFlux::frun, simb::MCFlux::ftptype, simb::MCFlux::ftpx, simb::MCFlux::ftpy, simb::MCFlux::ftpz, simb::MCFlux::fvx, simb::MCFlux::fvy, and simb::MCFlux::fvz.

  {
    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.

Referenced by larg4::ParticleListAction::ParticleInfo_t::clear(), larg4::ParticleListActionService::ParticleInfo_t::clear(), larg4::ParticleListActionService::GetPrimaryTruthIndex(), larg4::ParticleListAction::GetPrimaryTruthIndex(), isGeneratedParticleIndex(), g4b::PrimaryParticleInformation::IsInMCTruth(), g4b::PrimaryParticleInformation::MCParticleIndex(), larg4::ParticleListAction::PreTrackingAction(), and larg4::ParticleListActionService::preUserTrackingAction().