GENIE's GHEP MC event record. More...

#include <GHepRecord.h>

Inheritance diagram for genie::GHepRecord:

Collaboration diagram for genie::GHepRecord:

Public Member Functions
	GHepRecord ()

	GHepRecord (int size)

	GHepRecord (const GHepRecord &record)

	GHepRecord (TRootIOCtor *)

virtual	~GHepRecord ()

virtual Interaction *	Summary (void) const

virtual void	AttachSummary (Interaction *interaction)

virtual void	AddParticle (const GHepParticle &p)

virtual void	AddParticle (int pdg, GHepStatus_t ist, int mom1, int mom2, int dau1, int dau2, const TLorentzVector &p, const TLorentzVector &v)

virtual void	AddParticle (int pdg, GHepStatus_t ist, int mom1, int mom2, int dau1, int dau2, double px, double py, double pz, double E, double x, double y, double z, double t)

virtual GHepParticle *	Particle (int position) const

virtual GHepParticle *	FindParticle (int pdg, GHepStatus_t ist, int start) const

virtual int	ParticlePosition (int pdg, GHepStatus_t i, int start=0) const

virtual int	ParticlePosition (GHepParticle *particle, int start=0) const

virtual vector< int > *	GetStableDescendants (int position) const

GEvGenMode_t	EventGenerationMode (void) const

virtual GHepParticle *	Probe (void) const

virtual GHepParticle *	TargetNucleus (void) const

virtual GHepParticle *	RemnantNucleus (void) const

virtual GHepParticle *	HitNucleon (void) const

virtual GHepParticle *	HitElectron (void) const

virtual GHepParticle *	FinalStatePrimaryLepton (void) const

virtual GHepParticle *	FinalStateHadronicSystem (void) const

virtual int	ProbePosition (void) const

virtual int	TargetNucleusPosition (void) const

virtual int	RemnantNucleusPosition (void) const

virtual int	HitNucleonPosition (void) const

virtual int	HitElectronPosition (void) const

virtual int	FinalStatePrimaryLeptonPosition (void) const

virtual int	FinalStateHadronicSystemPosition (void) const

virtual unsigned int	NEntries (int pdg, GHepStatus_t ist, int start=0) const

virtual unsigned int	NEntries (int pdg, int start=0) const

virtual TBits *	EventFlags (void) const

virtual TBits *	EventMask (void) const

virtual bool	IsUnphysical (void) const

virtual bool	Accept (void) const

virtual double	Weight (void) const

virtual double	Probability (void) const

virtual double	XSec (void) const

virtual double	DiffXSec (void) const

virtual KinePhaseSpace_t	DiffXSecVars (void) const

virtual void	SetWeight (double wght)

virtual void	SetProbability (double prob)

virtual void	SetXSec (double xsec)

virtual void	SetDiffXSec (double xsec, KinePhaseSpace_t ps)

virtual TLorentzVector *	Vertex (void) const

virtual void	SetVertex (double x, double y, double z, double t)

virtual void	SetVertex (const TLorentzVector &vtx)

virtual void	Copy (const GHepRecord &record)

virtual void	Clear (Option_t *opt="")

virtual void	ResetRecord (void)

virtual void	CompactifyDaughterLists (void)

virtual void	RemoveIntermediateParticles (void)

void	SetUnphysEventMask (const TBits &mask)

void	Print (ostream &stream) const

Static Public Member Functions
static void	SetPrintLevel (int print_level)

static int	GetPrintLevel ()

Protected Member Functions
void	InitRecord (void)

void	CleanRecord (void)

virtual void	UpdateDaughterLists (void)

virtual bool	HasCompactDaughterList (int pos)

virtual void	SwapParticles (int i, int j)

virtual void	FinalizeDaughterLists (void)

virtual int	FirstNonInitStateEntry (void)

Protected Attributes
Interaction *	fInteraction
	attached summary information More...

TLorentzVector *	fVtx
	vertex in the detector coordinate system More...

TBits *	fEventFlags
	event flags indicating various pathologies or an unphysical event More...

TBits *	fEventMask
	an input bit-field mask allowing one to ignore bits set in fEventFlags More...

double	fWeight
	event weight More...

double	fProb
	event probability (for given flux neutrino and density-weighted path-length for target element) More...

double	fXSec
	cross section for selected event More...

double	fDiffXSec
	differential cross section for selected event kinematics More...

KinePhaseSpace_t	fDiffXSecPhSp
	specifies which differential cross-section (dsig/dQ2, dsig/dQ2dW, dsig/dxdy,...) More...

Static Protected Attributes
static int	fPrintLevel = 3

Friends
ostream &	operator<< (ostream &stream, const GHepRecord &event)

Detailed Description

GENIE's GHEP MC event record.

Author: Costas Andreopoulos <c.andreopoulos cern.ch> University of Liverpool

Created:: October 1, 2004

License:: Copyright (c) 2003-2024, The GENIE Collaboration For the full text of the license visit http://copyright.genie-mc.org

Definition at line 45 of file GHepRecord.h.

Constructor & Destructor Documentation

GHepRecord::GHepRecord ( )

Definition at line 53 of file GHepRecord.cxx.

References InitRecord().

                        :
 TClonesArray("genie::GHepParticle")
 {
   this->InitRecord();
 }

GHepRecord::GHepRecord ( int size )

Definition at line 59 of file GHepRecord.cxx.

References InitRecord().

                                :
 TClonesArray("genie::GHepParticle", size)
 {
   this->InitRecord();
 }

GHepRecord::GHepRecord ( const GHepRecord & record )

Definition at line 65 of file GHepRecord.cxx.

References Copy(), and InitRecord().

                                                 :
 TClonesArray("genie::GHepParticle", record.GetEntries())
 {
   this->InitRecord();
   this->Copy(record);
 }

GHepRecord::GHepRecord ( TRootIOCtor * )

Definition at line 72 of file GHepRecord.cxx.

                                    :
 TClonesArray("genie::GHepParticle"),
 fInteraction(0),
 fVtx(0),
 fEventFlags(0),
 fEventMask(0),
 fWeight(0.),
 fProb(0.),
 fXSec(0.),
 fDiffXSec(0.)
 {
 
 }

GHepRecord::~GHepRecord ( )

virtual

Definition at line 86 of file GHepRecord.cxx.

References CleanRecord().

 {
   this->CleanRecord();
 }

Member Function Documentation

bool GHepRecord::Accept ( void ) const

virtual

Definition at line 947 of file GHepRecord.cxx.

References fEventFlags, and fEventMask.

Referenced by Print(), and genie::EventGenerator::ProcessEventRecord().

 {
   TBits flags = *fEventFlags;
   TBits mask  = *fEventMask;
   TBits bitwiseand = flags & mask;
   bool accept = (bitwiseand.CountBits() == 0);
   return accept;
 }

void GHepRecord::AddParticle ( const GHepParticle & p )

virtual

Definition at line 499 of file GHepRecord.cxx.

References LOG, genie::GHepParticle::Pdg(), pINFO, and UpdateDaughterLists().

Referenced by genie::HNIntranuke2018::AbsorbHN(), genie::NuETargetRemnantGenerator::AddElectronNeutrino(), genie::DMETargetRemnantGenerator::AddElectronNeutrino(), genie::HadronicSystemGenerator::AddFinalHadronicSyst(), genie::AMNuGammaGenerator::AddFinalStateNeutrino(), genie::hnl::Decayer::AddInitialState(), genie::InitialStateAppender::AddNeutrino(), genie::InitialStateAppender::AddNucleus(), genie::AMNuGammaGenerator::AddPhoton(), genie::NucDeExcitationSim::AddPhoton(), genie::HEDISGenerator::AddPrimaryLepton(), genie::QELHadronicSystemGenerator::AddRecoilBaryon(), genie::IBDHadronicSystemGenerator::AddRecoilBaryon(), genie::AMNuGammaGenerator::AddRecoilNucleon(), genie::RSPPResonanceSelector::AddResonance(), genie::RESHadronicSystemGenerator::AddResonance(), genie::RSPPHadronicSystemGenerator::AddResonanceDecayProducts(), genie::InitialStateAppender::AddStruckParticle(), genie::NuETargetRemnantGenerator::AddTargetNucleusRemnant(), genie::DMETargetRemnantGenerator::AddTargetNucleusRemnant(), genie::HadronicSystemGenerator::AddTargetNucleusRemnant(), genie::QELEventGenerator::AddTargetNucleusRemnant(), genie::FermiMover::AddTargetNucleusRemnant(), genie::DMELEventGenerator::AddTargetNucleusRemnant(), genie::QELEventGeneratorSM::AddTargetNucleusRemnant(), genie::AMNuGammaGenerator::AddTargetRemnant(), genie::PrimaryLeptonGenerator::AddToEventRecord(), genie::OutgoingDarkGenerator::AddToEventRecord(), genie::COHHadronicSystemGenerator::CalculateHadronicSystem_AlvarezRuso(), genie::SKHadronicSystemGenerator::CalculateHadronicSystem_AtharSingleKaon(), genie::COHHadronicSystemGenerator::CalculateHadronicSystem_BergerSehgalFM(), genie::COHHadronicSystemGenerator::CalculateHadronicSystem_ReinSehgal(), ConvertToGHepMock(), genie::HAIntranuke::ElasHA(), genie::HAIntranuke2018::ElasHA(), genie::HNIntranuke2018::ElasHN(), genie::SecondNucleonEmissionI::EmitSecondNucleon(), genie::HNIntranuke2018::GammaInelasticHN(), genie::HNIntranuke2018::HandleCompoundNucleus(), genie::HAIntranuke::Inelastic(), genie::HAIntranuke2018::Inelastic(), genie::HAIntranuke::InelasticHA(), genie::HAIntranuke2018::InelasticHA(), genie::HNIntranuke2018::InelasticHN(), InitializeEvent(), genie::utils::intranuke::PhaseSpaceDecay(), genie::utils::intranuke2018::PhaseSpaceDecay(), genie::utils::intranuke::PionProduction(), genie::utils::intranuke2018::PionProduction(), genie::DFRHadronicSystemGenerator::ProcessEventRecord(), genie::QELEventGenerator::ProcessEventRecord(), genie::PhotonCOHGenerator::ProcessEventRecord(), genie::PhotonRESGenerator::ProcessEventRecord(), genie::NormGenerator::ProcessEventRecord(), genie::NucBindEnergyAggregator::ProcessEventRecord(), genie::DMELEventGenerator::ProcessEventRecord(), genie::QELEventGeneratorSM::ProcessEventRecord(), genie::SPPEventGenerator::ProcessEventRecord(), genie::hnl::VertexGenerator::ProcessEventRecord(), genie::hnl::FluxCreator::ProcessEventRecord(), genie::HAIntranuke::SimulateHadronicFinalState(), genie::HAIntranuke2018::SimulateHadronicFinalState(), genie::HNIntranuke2018::SimulateHadronicFinalState(), genie::Intranuke::TransportHadrons(), genie::Intranuke2018::TransportHadrons(), and genie::HadronTransporter::TransportInTransparentNuc().

 {
 // Provides a simplified method for inserting entries in the TClonesArray
 
   unsigned int pos = this->GetEntries();
 
 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
   LOG("GHEP", pINFO)
     << "Adding particle with pdgc = " << p.Pdg() << " at slot = " << pos;
 #endif
   new ((*this)[pos]) GHepParticle(p);
 
   // Update the mother's daughter list. If the newly inserted particle broke
   // compactification, then run CompactifyDaughterLists()
   this->UpdateDaughterLists();
 }

void GHepRecord::AddParticle	(	int	pdg,
		GHepStatus_t	ist,
		int	mom1,
		int	mom2,
		int	dau1,
		int	dau2,
		const TLorentzVector &	p,
		const TLorentzVector &	v
	)

virtual

Definition at line 516 of file GHepRecord.cxx.

References LOG, pINFO, and UpdateDaughterLists().

 {
 // Provides a 'simplified' method for inserting entries in the TClonesArray
 
   unsigned int pos = this->GetEntries();
 
 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
   LOG("GHEP", pINFO)
            << "Adding particle with pdgc = " << pdg << " at slot = " << pos;
 #endif
   new ((*this)[pos]) GHepParticle(pdg,status, mom1,mom2,dau1,dau2, p, v);
 
   // Update the mother's daughter list. If the newly inserted particle broke
   // compactification, then run CompactifyDaughterLists()
   this->UpdateDaughterLists();
 }

void GHepRecord::AddParticle	(	int	pdg,
		GHepStatus_t	ist,
		int	mom1,
		int	mom2,
		int	dau1,
		int	dau2,
		double	px,
		double	py,
		double	pz,
		double	E,
		double	x,
		double	y,
		double	z,
		double	t
	)

virtual

Definition at line 535 of file GHepRecord.cxx.

References LOG, pINFO, and UpdateDaughterLists().

 {
 // Provides a 'simplified' method for inserting entries in the TClonesArray
 
   unsigned int pos = this->GetEntries();
 
 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
   LOG("GHEP", pINFO)
            << "Adding particle with pdgc = " << pdg << " at slot = " << pos;
 #endif
   new ( (*this)[pos] ) GHepParticle (
             pdg, status, mom1, mom2, dau1, dau2, px, py, pz, E, x, y, z, t);
 
   // Update the mother's daughter list. If the newly inserted particle broke
   // compactification, then run CompactifyDaughterLists()
   this->UpdateDaughterLists();
 }

void GHepRecord::AttachSummary ( Interaction * interaction )

virtual

Definition at line 99 of file GHepRecord.cxx.

References fInteraction.

Referenced by InitializeEvent(), genie::PhysInteractionSelector::SelectInteraction(), and genie::ToyInteractionSelector::SelectInteraction().

 {
   fInteraction = interaction;
 }

void GHepRecord::CleanRecord ( void )

protected

Definition at line 866 of file GHepRecord.cxx.

References Clear(), LOG, and pDEBUG.

Referenced by ResetRecord(), and ~GHepRecord().

 {
 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
   LOG("GHEP", pDEBUG) << "Cleaning up GHepRecord";
 #endif
   this->Clear("C");
 }

void GHepRecord::Clear ( Option_t * opt = "" )

virtual

Definition at line 883 of file GHepRecord.cxx.

References fEventFlags, fEventMask, fInteraction, and fVtx.

Referenced by CleanRecord().

 {
   if (fInteraction) delete fInteraction;
   fInteraction=0;
 
   if (fVtx) delete fVtx;
   fVtx=0;
 
   if(fEventFlags) delete fEventFlags;
   fEventFlags=0;
 
   if(fEventMask) delete fEventMask;
   fEventMask=0;
 
   TClonesArray::Clear(opt);
 
 //  if (fInteraction) delete fInteraction;
 //  delete fVtx;
 //
 //  delete fEventFlags;
 //
 //  TClonesArray::Clear(opt);
 }

void GHepRecord::CompactifyDaughterLists ( void )

virtual

compact

Definition at line 651 of file GHepRecord.cxx.

References FinalizeDaughterLists(), genie::GHepParticle::FirstDaughter(), genie::GHepParticle::FirstMother(), HasCompactDaughterList(), genie::GHepParticle::LastDaughter(), LOG, Particle(), pINFO, pNOTICE, genie::GHepParticle::SetFirstDaughter(), genie::GHepParticle::SetLastDaughter(), and SwapParticles().

Referenced by UpdateDaughterLists().

 {
   int n = this->GetEntries();
   if(n<1) return;
 
   int i = this->Particle(n-1)->FirstMother();
   if(i<0) return;
 
   //  for(int i=0; i<n; i++) {
      bool compact = this->HasCompactDaughterList(i);
 
 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
      LOG("GHEP", pNOTICE)
         << "Particle's " << i << " daughter list is "
         << ((compact) ? "compact" : "__not__ compact");
 #endif
 
      if(!compact) {
         GHepParticle * p = this->Particle(i);
 
         int dau1 = p->FirstDaughter();
         int dau2 = p->LastDaughter();
         int ndau = dau2-dau1+1;
         int ndp  = dau2+1;
         if(dau1==-1) {ndau=0;}
 
         int curr_pos = n-1;
         while (curr_pos > ndp) {
            this->SwapParticles(curr_pos,curr_pos-1);
            curr_pos--;
         }
         if(ndau>0) {
           this->Particle(i)->SetFirstDaughter(dau1);
           this->Particle(i)->SetLastDaughter(dau2+1);
         } else {
           this->Particle(i)->SetFirstDaughter(-1);
           this->Particle(i)->SetLastDaughter(-1);
         }
         this->FinalizeDaughterLists();
 
      } //!compact
 
 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
      LOG("GHEP", pINFO)
         << "Done ompactifying daughter-list for particle at slot: " << i;
 #endif
      //  }
 }

void GHepRecord::Copy ( const GHepRecord & record )

virtual

Definition at line 907 of file GHepRecord.cxx.

References EventFlags(), EventMask(), fDiffXSec, fDiffXSecPhSp, fEventFlags, fEventMask, fInteraction, fProb, fVtx, fWeight, fXSec, ResetRecord(), and Vertex().

Referenced by genie::EventRecord::Copy(), GHepRecord(), and genie::EventGenerator::ProcessEventRecord().

 {
   // clean up
   this->ResetRecord();
 
   // copy event record entries
   unsigned int ientry = 0;
   GHepParticle * p = 0;
   TIter ghepiter(&record);
   while ( (p = (GHepParticle *) ghepiter.Next()) )
                               new ( (*this)[ientry++] ) GHepParticle(*p);
 
   // copy summary
   fInteraction = new Interaction( *record.fInteraction );
 
   // copy flags & mask
   *fEventFlags = *(record.EventFlags());
   *fEventMask  = *(record.EventMask());
 
   // copy vtx position
   TLorentzVector * v = record.Vertex();
   fVtx->SetXYZT(v->X(),v->Y(),v->Z(),v->T());
 
   // copy weights & xsecs
   fWeight       = record.fWeight;
   fProb         = record.fProb;
   fXSec         = record.fXSec;
   fDiffXSec     = record.fDiffXSec;
   fDiffXSecPhSp = record.fDiffXSecPhSp;
 }

virtual double genie::GHepRecord::DiffXSec ( void ) const

inlinevirtual

Definition at line 127 of file GHepRecord.h.

References fDiffXSec.

127 { return fDiffXSec; }

genie::GHepRecord::fDiffXSec

double fDiffXSec

differential cross section for selected event kinematics

Definition: GHepRecord.h:181

virtual KinePhaseSpace_t genie::GHepRecord::DiffXSecVars ( void ) const

inlinevirtual

Definition at line 128 of file GHepRecord.h.

References fDiffXSecPhSp.

128 { return fDiffXSecPhSp; }

genie::GHepRecord::fDiffXSecPhSp

KinePhaseSpace_t fDiffXSecPhSp

specifies which differential cross-section (dsig/dQ2, dsig/dQ2dW, dsig/dxdy,...)

Definition: GHepRecord.h:182

virtual TBits* genie::GHepRecord::EventFlags ( void ) const

inlinevirtual

Definition at line 117 of file GHepRecord.h.

References fEventFlags.

117 { return fEventFlags; }

genie::GHepRecord::fEventFlags

TBits * fEventFlags

event flags indicating various pathologies or an unphysical event

Definition: GHepRecord.h:174

GEvGenMode_t GHepRecord::EventGenerationMode ( void ) const

Definition at line 209 of file GHepRecord.cxx.

References genie::pdg::IsAntiDarkMatter(), genie::pdg::IsDarkMatter(), genie::pdg::IsHadron(), genie::pdg::IsHNL(), genie::pdg::IsIon(), genie::pdg::IsLepton(), genie::pdg::IsNucleon(), genie::kGMdDarkMatterNucleus, genie::kGMdHadronNucleus, genie::kGMdHNLDecay, genie::kGMdLeptonNucleus, genie::kGMdNucleonDecay, genie::kGMdPhotonNucleus, genie::kGMdUnknown, genie::kIStDecayedState, genie::kIStInitialState, genie::kPdgGamma, Particle(), genie::GHepParticle::Pdg(), and genie::GHepParticle::Status().

Referenced by ProbePosition(), genie::Intranuke2018::ProcessEventRecord(), genie::Intranuke::ProcessEventRecord(), and TargetNucleusPosition().

 {
   GHepParticle * p0 = this->Particle(0);
   if(!p0) return kGMdUnknown;
   GHepParticle * p1 = this->Particle(1);
   if(!p1) return kGMdUnknown;
 
   int p0pdg = p0->Pdg();
   GHepStatus_t p0st = p0->Status();
   int p1pdg = p1->Pdg();
   GHepStatus_t p1st = p1->Status();
 
   // In lepton+nucleon/nucleus mode, the 1st entry in the event record
   // is a charged or neutral lepton with status code = kIStInitialState
   if( pdg::IsLepton(p0pdg) && p0st == kIStInitialState )
   {
     return kGMdLeptonNucleus;
   }
 
   // In dark matter mode, the 1st entry in the event record is a dark
   // matter particle
   if( (pdg::IsDarkMatter(p0pdg) || pdg::IsAntiDarkMatter(p0pdg)) && p0st == kIStInitialState )
   {
     return kGMdDarkMatterNucleus;
   }
 
   // In hadron+nucleon/nucleus mode, the 1st entry in the event record
   // is a hadron with status code = kIStInitialState and the 2nd entry
   // is a nucleon or nucleus with status code = kIStInitialState
   if( pdg::IsHadron(p0pdg) && p0st == kIStInitialState )
   {
     if( (pdg::IsIon(p1pdg) || pdg::IsNucleon(p1pdg)) && p1st == kIStInitialState)
     {
        return kGMdHadronNucleus;
     }
   }
 
   // As above, with a photon as a probe
   if( p0pdg == kPdgGamma && p0st == kIStInitialState )
   {
     if( (pdg::IsIon(p1pdg) || pdg::IsNucleon(p1pdg)) && p1st == kIStInitialState)
     {
        return kGMdPhotonNucleus;
     }
   }
 
   // In nucleon decay mode,
   // - [if the decayed nucleon was a bound one] the 1st entry in the event
   //   record is a nucleus with status code = kIStInitialState and the
   //   2nd entry is a nucleon with code = kIStDecayedState
   // - [if the decayed nucleon was a free one] the first entry in the event
   //   record is a nucleon with status code = kIStInitialState and it has a
   //   single daughter which is a nucleon with status code = kIStDecayedState.
 
   if( pdg::IsIon(p0pdg)     && p0st == kIStInitialState &&
       pdg::IsNucleon(p1pdg) && p1st == kIStDecayedState)
   {
      return kGMdNucleonDecay;
   }
   if( pdg::IsNucleon(p0pdg) && p0st == kIStInitialState &&
       pdg::IsNucleon(p1pdg) && p1st == kIStDecayedState)
   {
      return kGMdNucleonDecay;
   }
 
   // In HNL decay mode, the first entry is a HNL.
   
   if( pdg::IsHNL(p0pdg) && p0st == kIStInitialState )
   {
      return kGMdHNLDecay;
   }
 
   return kGMdUnknown;
 }

virtual TBits* genie::GHepRecord::EventMask ( void ) const

inlinevirtual

Definition at line 118 of file GHepRecord.h.

References fEventMask.

Referenced by Copy().

118 { return fEventMask; }

genie::GHepRecord::fEventMask

TBits * fEventMask

an input bit-field mask allowing one to ignore bits set in fEventFlags

Definition: GHepRecord.h:175

void GHepRecord::FinalizeDaughterLists ( void )

protectedvirtual

Definition at line 799 of file GHepRecord.cxx.

References genie::GHepParticle::FirstMother().

Referenced by CompactifyDaughterLists().

 {
 // Update all daughter-lists based on particle 'first mother' field.
 // To work correctly, the daughter-lists must have been compactified first.
 
   GHepParticle * p1 = 0;
   TIter iter1(this);
   int i1=0;
   while( (p1 = (GHepParticle *)iter1.Next()) ) {
     int dau1 = -1;
     int dau2 = -1;
     GHepParticle * p2 = 0;
     TIter iter2(this);
     int i2=0;
     while( (p2 = (GHepParticle *)iter2.Next()) ) {
 
        if(p2->FirstMother() == i1) {
           dau1 = (dau1<0) ? i2 : TMath::Min(dau1,i2);
           dau2 = (dau2<0) ? i2 : TMath::Max(dau2,i2);
        }
        i2++;
     }
     i1++;
     p1 -> SetFirstDaughter (dau1);
     p1 -> SetLastDaughter  (dau2);
   }
 }

GHepParticle * GHepRecord::FinalStateHadronicSystem ( void ) const

virtual

Definition at line 342 of file GHepRecord.cxx.

References FinalStateHadronicSystemPosition(), and Particle().

Referenced by genie::DISHadronicSystemGenerator::SimulateFormationZone().

 {
 // Returns the GHepParticle representing the sum of the DIS pre-fragm f/s
 // hadronic system, or 0 if it does not exist.
 
   int ipos = this->FinalStateHadronicSystemPosition();
   if(ipos>-1) return this->Particle(ipos);
   return 0;
 }

int GHepRecord::FinalStateHadronicSystemPosition ( void ) const

virtual

Definition at line 471 of file GHepRecord.cxx.

References genie::kIStDISPreFragmHadronicState, genie::kPdgHadronicSyst, and ParticlePosition().

Referenced by FinalStateHadronicSystem().

 {
   return this->ParticlePosition(
                         kPdgHadronicSyst,kIStDISPreFragmHadronicState,0);
 }

GHepParticle * GHepRecord::FinalStatePrimaryLepton ( void ) const

virtual

Definition at line 333 of file GHepRecord.cxx.

References FinalStatePrimaryLeptonPosition(), and Particle().

Referenced by genie::NuETargetRemnantGenerator::AddElectronNeutrino(), genie::DMETargetRemnantGenerator::AddElectronNeutrino(), genie::COHHadronicSystemGenerator::CalculateHadronicSystem_AlvarezRuso(), genie::SKHadronicSystemGenerator::CalculateHadronicSystem_AtharSingleKaon(), genie::COHHadronicSystemGenerator::CalculateHadronicSystem_BergerSehgalFM(), genie::COHHadronicSystemGenerator::CalculateHadronicSystem_ReinSehgal(), genie::hnl::Decayer::GenerateDecayProducts(), genie::HadronicSystemGenerator::Hadronic4pLAB(), genie::HadronicSystemGenerator::MomentumTransferLAB(), genie::DFRHadronicSystemGenerator::ProcessEventRecord(), genie::DFRPrimaryLeptonGenerator::ProcessEventRecord(), genie::DISPrimaryLeptonGenerator::ProcessEventRecord(), genie::SKPrimaryLeptonGenerator::ProcessEventRecord(), genie::DMDISOutgoingDarkGenerator::ProcessEventRecord(), genie::OutgoingDarkGenerator::SetPolarization(), genie::utils::SetPrimaryLeptonPolarization(), and genie::hnl::Decayer::UpdateEventRecord().

 {
 // Returns the GHepParticle representing the final state primary lepton.
 
   int ipos = this->FinalStatePrimaryLeptonPosition();
   if(ipos>-1) return this->Particle(ipos);
   return 0;
 }

int GHepRecord::FinalStatePrimaryLeptonPosition ( void ) const

virtual

Definition at line 459 of file GHepRecord.cxx.

References genie::GHepParticle::FirstDaughter(), and Probe().

Referenced by FinalStatePrimaryLepton().

 {
 // Returns the GHEP position GHepParticle representing the final state
 // primary lepton.
 
   GHepParticle * probe = this->Probe();
   if(!probe) return -1;
 
   int ifsl = probe->FirstDaughter();
   return ifsl;
 }

GHepParticle * GHepRecord::FindParticle	(	int	pdg,
		GHepStatus_t	ist,
		int	start
	)		const

virtual

Definition at line 118 of file GHepRecord.cxx.

References LOG, genie::GHepParticle::Pdg(), pINFO, and genie::GHepParticle::Status().

Referenced by genie::NucBindEnergyAggregator::ProcessEventRecord().

 {
 // Returns the first GHepParticle with the input pdg-code and status
 // starting from the specified position of the event record.
 
   int nentries = this->GetEntries();
   for(int i = start; i < nentries; i++) {
      GHepParticle * p = (GHepParticle *) (*this)[i];
      if(p->Status() == status && p->Pdg() == pdg) return p;
   }
 
   LOG("GHEP", pINFO)
     << "No particle found with: (pos >= " << start
     << ", pdg = " << pdg << ", ist = " << status << ")";
 
   return 0;
 }

int GHepRecord::FirstNonInitStateEntry ( void )

protectedvirtual

Definition at line 739 of file GHepRecord.cxx.

References genie::kIStInitialState, genie::kIStNucleonTarget, and genie::GHepParticle::Status().

 {
   GHepParticle * p = 0;
   TIter iter(this);
   int pos = 0;
   while( (p = (GHepParticle *)iter.Next()) ) {
     int ist = p->Status();
     if(ist != kIStInitialState && ist != kIStNucleonTarget) return pos;
     pos++;
   }
   return pos;
 }

int GHepRecord::GetPrintLevel ( )

static

Definition at line 960 of file GHepRecord.cxx.

References fPrintLevel.

 {
   return fPrintLevel;
 }

vector< int > * GHepRecord::GetStableDescendants ( int position ) const

virtual

Definition at line 174 of file GHepRecord.cxx.

References genie::GHepParticle::FirstMother(), genie::kIStStableFinalState, Particle(), and genie::GHepParticle::Status().

Referenced by genie::utils::intranuke::PreEquilibrium(), and genie::utils::intranuke2018::PreEquilibrium().

 {
 // Returns a list of all stable descendants of the GHEP entry in the input
 // slot. The user adopts the output vector.
 
   // return null if particle index is out of range
   int nentries = this->GetEntries();
   if(position<0 || position>=nentries) return 0;
 
   vector<int> * descendants = new vector<int>;
 
   // return itself if it is a stable final state particle
   if(this->Particle(position)->Status() == kIStStableFinalState) {
     descendants->push_back(position);
     return descendants;
   }
 
   for(int i = 0; i < nentries; i++) {
     if(i==position) continue;
     GHepParticle * p = (GHepParticle *) (*this)[i];
     if(p->Status() != kIStStableFinalState) continue;
     bool is_descendant=false;
     int mom = p->FirstMother();
     while(mom>-1) {
       if(mom==position) is_descendant=true;
       if(is_descendant) {
         descendants->push_back(i);
         break;
       }
       mom = this->Particle(mom)->FirstMother();
     }
   }
   return descendants;
 }

bool GHepRecord::HasCompactDaughterList ( int pos )

protectedvirtual

Definition at line 700 of file GHepRecord.cxx.

References genie::GHepParticle::FirstMother(), LOG, pDEBUG, and pINFO.

Referenced by CompactifyDaughterLists().

 {
 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
   LOG("GHEP", pDEBUG) << "Examining daughter-list of particle at: " << pos;
 #endif
   vector<int> daughters;
   GHepParticle * p = 0;
   TIter iter(this);
   int i=0;
   while( (p = (GHepParticle *)iter.Next()) ) {
     if(p->FirstMother() == pos) {
 
 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
        LOG("GHEP", pDEBUG) << "Particle at: " << i << " is a daughter";
 #endif
        daughters.push_back(i);
     }
     i++;
   }
 
   bool is_compact = true;
   if(daughters.size()>1) {
     sort(daughters.begin(), daughters.end());
     vector<int>::iterator diter = daughters.begin();
     int prev = *diter;
     for(; diter != daughters.end(); ++diter) {
        int curr = *diter;
        is_compact = is_compact && (TMath::Abs(prev-curr)<=1);
        prev = curr;
     }
   }
 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
   LOG("GHEP", pINFO)
       << "Daughter-list of particle at: " << pos << " is "
                             << (is_compact ? "" : "not") << " compact";
 #endif
   return is_compact;
 }

GHepParticle * GHepRecord::HitElectron ( void ) const

virtual

Definition at line 323 of file GHepRecord.cxx.

References HitElectronPosition(), and Particle().

Referenced by genie::NuETargetRemnantGenerator::AddElectronNeutrino(), and genie::DMETargetRemnantGenerator::AddElectronNeutrino().

 {
 // Returns the GHepParticle representing the struck electron, or 0 if it does
 // not exist.
 
   int ipos = this->HitElectronPosition();
   if(ipos>-1) return this->Particle(ipos);
   return 0;
 }

int GHepRecord::HitElectronPosition ( void ) const

virtual

Definition at line 441 of file GHepRecord.cxx.

References genie::pdg::IsElectron(), genie::kIStInitialState, Particle(), genie::GHepParticle::Pdg(), genie::GHepParticle::Status(), and TargetNucleus().

Referenced by genie::NuETargetRemnantGenerator::AddElectronNeutrino(), genie::DMETargetRemnantGenerator::AddElectronNeutrino(), and HitElectron().

 {
 // Returns the GHEP position of the GHepParticle representing a hit electron.
 // Same as above..
 
   GHepParticle * nucleus = this->TargetNucleus();
 
   int ipos = (nucleus) ? 2 : 1;
 
   GHepParticle * p = this->Particle(ipos);
   if(!p) return -1;
 
   bool ise = pdg::IsElectron(p->Pdg());
   if(ise && p->Status()==kIStInitialState) return ipos;
 
   return -1;
 }

GHepParticle * GHepRecord::HitNucleon ( void ) const

virtual

Definition at line 313 of file GHepRecord.cxx.

References HitNucleonPosition(), and Particle().

 {
 // Returns the GHepParticle representing the struck nucleon, or 0 if it does
 // not exist.
 
   int ipos = this->HitNucleonPosition();
   if(ipos>-1) return this->Particle(ipos);
   return 0;
 }

int GHepRecord::HitNucleonPosition ( void ) const

virtual

Definition at line 418 of file GHepRecord.cxx.

References genie::pdg::Is2NucleonCluster(), genie::pdg::IsNucleon(), genie::kIStInitialState, genie::kIStNucleonTarget, Particle(), genie::GHepParticle::Pdg(), genie::GHepParticle::Status(), and TargetNucleus().

Referenced by genie::HadronicSystemGenerator::AddFinalHadronicSyst(), genie::QELHadronicSystemGenerator::AddRecoilBaryon(), genie::IBDHadronicSystemGenerator::AddRecoilBaryon(), genie::AMNuGammaGenerator::AddRecoilNucleon(), genie::RSPPResonanceSelector::AddResonance(), genie::RESHadronicSystemGenerator::AddResonance(), genie::SKHadronicSystemGenerator::CalculateHadronicSystem_AtharSingleKaon(), HitNucleon(), genie::DFRHadronicSystemGenerator::ProcessEventRecord(), genie::QELEventGenerator::ProcessEventRecord(), genie::DMELEventGenerator::ProcessEventRecord(), genie::QELEventGeneratorSM::ProcessEventRecord(), and genie::SPPEventGenerator::ProcessEventRecord().

 {
 // Returns the GHEP position of the GHepParticle representing the hit nucleon.
 // If a struck nucleon is set it will be at slot 2 (for scattering off nuclear
 // targets) or at slot 1 (for free nucleon scattering).
 // If the struck nucleon is not set (eg coherent scattering, ve- scattering)
 // it returns 0.
 
   GHepParticle * nucleus = this->TargetNucleus();
 
   int          ipos = (nucleus) ? 2 : 1;
   GHepStatus_t ist  = (nucleus) ? kIStNucleonTarget : kIStInitialState;
 
   GHepParticle * p = this->Particle(ipos);
   if(!p) return -1;
 
 //  bool isN = pdg::IsNeutronOrProton(p->Pdg());
   bool isN = pdg::IsNucleon(p->Pdg()) || pdg::Is2NucleonCluster(p->Pdg());
   if(isN && p->Status()==ist) return ipos;
 
   return -1;
 }

void GHepRecord::InitRecord ( void )

protected

Definition at line 837 of file GHepRecord.cxx.

References fDiffXSec, fDiffXSecPhSp, fEventFlags, fEventMask, fInteraction, fProb, fVtx, fWeight, fXSec, genie::kPSNull, LOG, genie::GHepFlags::NFlags(), pDEBUG, and pINFO.

Referenced by GHepRecord(), and ResetRecord().

 {
 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
   LOG("GHEP", pDEBUG) << "Initializing GHepRecord";
 #endif
   fInteraction  = 0;
   fWeight       = 1.;
   fProb         = 1.;
   fXSec         = 0.;
   fDiffXSec     = 0.;
   fDiffXSecPhSp = kPSNull;
   fVtx          = new TLorentzVector(0,0,0,0);
 
   fEventFlags  = new TBits(GHepFlags::NFlags());
   fEventFlags -> ResetAllBits(false);
 
   fEventMask   = new TBits(GHepFlags::NFlags());
 //fEventMask  -> ResetAllBits(true);
   for(unsigned int i = 0; i < GHepFlags::NFlags(); i++) {
    fEventMask->SetBitNumber(i, true);
   }
 
   LOG("GHEP", pINFO)
     << "Initialised unphysical event mask (bits: " << GHepFlags::NFlags() - 1
     << " -> 0) : " << *fEventMask;
 
   this->SetOwner(true);
 }

virtual bool genie::GHepRecord::IsUnphysical ( void ) const

inlinevirtual

Definition at line 119 of file GHepRecord.h.

References fEventFlags.

Referenced by GenerateEvent(), and Print().

119 { return (fEventFlags->CountBits()>0); }

genie::GHepRecord::fEventFlags

TBits * fEventFlags

event flags indicating various pathologies or an unphysical event

Definition: GHepRecord.h:174

unsigned int GHepRecord::NEntries	(	int	pdg,
		GHepStatus_t	ist,
		int	start = `0`
	)		const

virtual

Definition at line 477 of file GHepRecord.cxx.

References genie::GHepParticle::Pdg(), and genie::GHepParticle::Status().

 {
   unsigned int nentries = 0;
 
   for(int i = start; i < this->GetEntries(); i++) {
      GHepParticle * p = (GHepParticle *) (*this)[i];
      if(p->Pdg()==pdg && p->Status()==ist) nentries++;
   }
   return nentries;
 }

unsigned int GHepRecord::NEntries	(	int	pdg,
		int	start = `0`
	)		const

virtual

Definition at line 488 of file GHepRecord.cxx.

References genie::GHepParticle::Pdg().

 {
   unsigned int nentries = 0;
 
   for(int i = start; i < this->GetEntries(); i++) {
      GHepParticle * p = (GHepParticle *) (*this)[i];
      if(p->Pdg()==pdg) nentries++;
   }
   return nentries;
 }

GHepParticle * GHepRecord::Particle ( int position ) const

virtual

Definition at line 104 of file GHepRecord.cxx.

References LOG, and pINFO.

Referenced by AcceptEvent(), genie::AMNuGammaGenerator::AddFinalStateNeutrino(), genie::hnl::Decayer::AddInitialState(), genie::NucDeExcitationSim::AddPhoton(), genie::QELHadronicSystemGenerator::AddRecoilBaryon(), genie::IBDHadronicSystemGenerator::AddRecoilBaryon(), genie::RSPPHadronicSystemGenerator::AddResonanceDecayProducts(), genie::QELEventGenerator::AddTargetNucleusRemnant(), genie::FermiMover::AddTargetNucleusRemnant(), genie::DMELEventGenerator::AddTargetNucleusRemnant(), genie::QELEventGeneratorSM::AddTargetNucleusRemnant(), genie::AMNuGammaGenerator::AddTargetRemnant(), CompactifyDaughterLists(), EventGenerationMode(), FinalStateHadronicSystem(), FinalStatePrimaryLepton(), genie::hnl::Decayer::GenerateDecayProducts(), GetStableDescendants(), HitElectron(), HitElectronPosition(), HitNucleon(), HitNucleonPosition(), main(), genie::utils::ghep::NeutReactionCode(), genie::utils::intranuke::PreEquilibrium(), genie::utils::intranuke2018::PreEquilibrium(), Probe(), genie::HEDISGenerator::ProcessEventRecord(), genie::NucBindEnergyAggregator::ProcessEventRecord(), genie::PauliBlocker::ProcessEventRecord(), genie::hnl::VertexGenerator::ProcessEventRecord(), genie::hnl::FluxCreator::ProcessEventRecord(), RemnantNucleus(), RemnantNucleusPosition(), genie::HAIntranuke::SimulateHadronicFinalState(), genie::HAIntranuke2018::SimulateHadronicFinalState(), genie::HNIntranuke2018::SimulateHadronicFinalState(), genie::HAIntranuke::SimulateHadronicFinalStateKinematics(), genie::HAIntranuke2018::SimulateHadronicFinalStateKinematics(), SwapParticles(), TargetNucleus(), TargetNucleusPosition(), genie::Intranuke::TransportHadrons(), genie::Intranuke2018::TransportHadrons(), UpdateDaughterLists(), and genie::hnl::Decayer::UpdateEventRecord().

 {
 // Returns the GHepParticle from the specified position of the event record.
 
   if( position >=0 && position < this->GetEntries() ) {
      GHepParticle * particle = (GHepParticle *) (*this)[position];
      if(particle) return particle;
   }
   LOG("GHEP", pINFO)
     << "No particle found with: (pos = " << position << ")";
 
   return 0;
 }

int GHepRecord::ParticlePosition	(	int	pdg,
		GHepStatus_t	i,
		int	start = `0`
	)		const

virtual

Definition at line 137 of file GHepRecord.cxx.

References LOG, genie::GHepParticle::Pdg(), pINFO, and genie::GHepParticle::Status().

Referenced by FinalStateHadronicSystemPosition(), genie::utils::intranuke::PhaseSpaceDecay(), and genie::utils::intranuke2018::PhaseSpaceDecay().

 {
 // Returns the position of the first GHepParticle with the input pdg-code
 // and status starting from the specified position of the event record.
 
   int nentries = this->GetEntries();
   for(int i = start; i < nentries; i++) {
      GHepParticle * p = (GHepParticle *) (*this)[i];
      if(p->Status() == status && p->Pdg() == pdg) return i;
   }
 
   LOG("GHEP", pINFO)
     << "No particle found with: (pos >= " << start
     << ", pdg = " << pdg << ", ist = " << status << ")";
 
   return -1;
 }

int GHepRecord::ParticlePosition	(	GHepParticle *	particle,
		int	start = `0`
	)		const

virtual

Definition at line 156 of file GHepRecord.cxx.

References genie::GHepParticle::Compare(), LOG, and pINFO.

 {
 // Returns the position of the first match with the specified GHepParticle
 // starting from the specified position of the event record.
 
   int nentries = this->GetEntries();
   for(int i = start; i < nentries; i++) {
      GHepParticle * p = (GHepParticle *) (*this)[i];
      if( p->Compare(particle) ) return i;
   }
 
   LOG("GHEP", pINFO)
     << "No particle found with pos >= " << start
     << " matching particle: " << *particle;
 
   return -1;
 }

void GHepRecord::Print ( ostream & stream ) const

Definition at line 965 of file GHepRecord.cxx.

Referenced by genie::operator<<(), and genie::EventRecord::Print().

 {
   // Print levels:
   //  0 -> prints particle list
   //  1 -> prints particle list + event flags
   //  2 -> prints particle list + event flags + wght/xsec
   //  3 -> prints particle list + event flags + wght/xsec + summary
   // 10 -> as in level 0 but showing particle positions too
   // 11 -> as in level 1 but showing particle positions too
   // 12 -> as in level 2 but showing particle positions too
   // 13 -> as in level 3 but showing particle positions too
 
    bool accept_input_print_level =
        (fPrintLevel >= 0 && fPrintLevel <= 3) ||
        (fPrintLevel >=10 && fPrintLevel <=13);
 
   int printlevel = (accept_input_print_level) ? fPrintLevel : 3;
   int printlevel_orig = printlevel;
 
   bool showpos = false;
   if(printlevel >= 10) {
      printlevel-=10;
      showpos=true;
   }
 
   stream << "\n\n|";
   stream << setfill('-') << setw(115) << "|";
 
   stream << "\n|GENIE GHEP Event Record [print level: "
          << setfill(' ') << setw(3) << printlevel_orig << "]"
          << setfill(' ') << setw(73) << "|";
 
   stream << "\n|";
   stream << setfill('-') << setw(115) << "|";
 
   stream << "\n| ";
   stream << setfill(' ') << setw(6)  << "Idx | "
          << setfill(' ') << setw(16) << "Name | "
          << setfill(' ') << setw(6)  << "Ist | "
          << setfill(' ') << setw(13) << "PDG | "
          << setfill(' ') << setw(12) << "Mother  | "
          << setfill(' ') << setw(12) << "Daughter  | "
          << setfill(' ') << setw(10) << ((showpos) ? "Px(x) |" : "Px | ")
          << setfill(' ') << setw(10) << ((showpos) ? "Py(y) |" : "Py | ")
          << setfill(' ') << setw(10) << ((showpos) ? "Pz(z) |" : "Pz | ")
          << setfill(' ') << setw(10) << ((showpos) ?  "E(t) |" :  "E | ")
          << setfill(' ') << setw(10) << "m  | ";
 
   stream << "\n|";
   stream << setfill('-') << setw(115) << "|";
 
   GHepParticle * p = 0;
   TObjArrayIter piter(this);
   TVector3 polarization(0,0,0);
 
   unsigned int idx = 0;
 
   double sum_E  = 0;
   double sum_px = 0;
   double sum_py = 0;
   double sum_pz = 0;
 
   while( (p = (GHepParticle *) piter.Next()) ) {
 
      stream << "\n| ";
      stream << setfill(' ') << setw(3)  << idx++               << " | ";
      stream << setfill(' ') << setw(13) << p->Name()           << " | ";
      stream << setfill(' ') << setw(3)  << p->Status()         << " | ";
      stream << setfill(' ') << setw(10) << p->Pdg()            << " | ";
      stream << setfill(' ') << setw(3)  << p->FirstMother()    << " | ";
      stream << setfill(' ') << setw(3)  << p->LastMother()     << " | ";
      stream << setfill(' ') << setw(3)  << p->FirstDaughter()  << " | ";
      stream << setfill(' ') << setw(3)  << p->LastDaughter()   << " | ";
      stream << std::fixed << setprecision(3);
      stream << setfill(' ') << setw(7)  << p->Px()             << " | ";
      stream << setfill(' ') << setw(7)  << p->Py()             << " | ";
      stream << setfill(' ') << setw(7)  << p->Pz()             << " | ";
      stream << setfill(' ') << setw(7)  << p->E()              << " | ";
 
      if (p->IsOnMassShell())
         stream << setfill(' ') << setw(7)  << p->Mass()        << " | ";
      else
         stream << setfill('*') << setw(7)  << p->Mass()        << " | M = "
                << p->P4()->M() << " ";
 
      if (p->PolzIsSet()) {
        p->GetPolarization(polarization);
        stream << "P = (" << polarization.x() << "," << polarization.y()
               << "," << polarization.z() << ")";
      }
 
      if (p->RescatterCode() != -1) {
        stream << "FSI = " << p->RescatterCode();
      }
 
      // plot particle position if requested
      if(showpos) {
        stream << "\n| ";
        stream << setfill(' ') << setw(6)  << " | ";
        stream << setfill(' ') << setw(16) << " | ";
        stream << setfill(' ') << setw(6)  << " | ";
        stream << setfill(' ') << setw(13) << " | ";
        stream << setfill(' ') << setw(6)  << " | ";
        stream << setfill(' ') << setw(6)  << " | ";
        stream << setfill(' ') << setw(6)  << " | ";
        stream << setfill(' ') << setw(6)  << " | ";
        stream << std::fixed  << setprecision(3);
        stream << setfill(' ') << setw(7)  << p->Vx()  << " | ";
        stream << setfill(' ') << setw(7)  << p->Vy()  << " | ";
        stream << setfill(' ') << setw(7)  << p->Vz()  << " | ";
        stream << setfill(' ') << setw(7)  << p->Vt()  << " | ";
        stream << setfill(' ') << setw(10) << " | ";
      }
 
      // compute P4_{final} - P4_{nitial}
 
      if(p->Status() == kIStStableFinalState ||
         p->Status() == kIStFinalStateNuclearRemnant) {
 
           sum_E  += p->E();
           sum_px += p->Px();
           sum_py += p->Py();
           sum_pz += p->Pz();
      } else
      if(p->Status() == kIStInitialState) {
        /*
      if(p->Status() == kIStInitialState || p->Status() == kIStNucleonTarget) {
        */
           sum_E  -= p->E();
           sum_px -= p->Px();
           sum_py -= p->Py();
           sum_pz -= p->Pz();
      }
 
   } // loop over particles
 
   stream << "\n|";
   stream << setfill('-') << setw(115) << "|";
 
   // Print SUMS
   stream << "\n| ";
   stream << setfill(' ') << setw(17) << "Fin-Init:  "
          << setfill(' ') << setw(6)  << "      "
          << setfill(' ') << setw(18) << "      "
          << setfill(' ') << setw(12) << "          "
          << setfill(' ') << setw(12) << "          | ";
   stream << std::fixed  << setprecision(3);
   stream << setfill(' ') << setw(7)  << sum_px  << " | ";
   stream << setfill(' ') << setw(7)  << sum_py  << " | ";
   stream << setfill(' ') << setw(7)  << sum_pz  << " | ";
   stream << setfill(' ') << setw(7)  << sum_E   << " | ";
   stream << setfill(' ') << setw(10)  << "   | ";
 
   stream << "\n|";
   stream << setfill('-') << setw(115) << "|";
 
   // Print vertex
 
   GHepParticle * probe = this->Probe();
   if(probe){
     stream << "\n| ";
     stream << setfill(' ') << setw(17) << "Vertex:    ";
     stream << setfill(' ') << setw(11)
                        << ((probe) ? probe->Name() : "unknown probe") << " @ (";
 
     stream << std::fixed  << setprecision(5);
     stream << "x = " << setfill(' ') << setw(11) << this->Vertex()->X() << " m, ";
     stream << "y = " << setfill(' ') << setw(11) << this->Vertex()->Y() << " m, ";
     stream << "z = " << setfill(' ') << setw(11) << this->Vertex()->Z() << " m, ";
     stream << std::scientific << setprecision(6);
     stream << "t = " << setfill(' ') << setw(15) << this->Vertex()->T() << " s) ";
     stream << std::fixed  << setprecision(3);
     stream << setfill(' ') << setw(2)  << "|";
 
     stream << "\n|";
     stream << setfill('-') << setw(115) << "|";
   }
 
   // Print FLAGS
 
   if(printlevel>=1) {
     stream << "\n| ";
     stream << "Err flag [bits:" << fEventFlags->GetNbits()-1 << "->0] : "
            << *fEventFlags << "    |  "
            << "1st set: " << setfill(' ') << setw(56)
            << ( this->IsUnphysical() ?
                  GHepFlags::Describe(GHepFlag_t(fEventFlags->FirstSetBit())) :
                  "none") << " | ";
     stream << "\n| ";
     stream << "Err mask [bits:" << fEventMask->GetNbits()-1 << "->0] : "
            << *fEventMask << "    |  "
            << "Is unphysical: " << setfill(' ') << setw(5)
            << utils::print::BoolAsYNString(this->IsUnphysical()) << " |   "
            << "Accepted: " << setfill(' ') << setw(5)
            << utils::print::BoolAsYNString(this->Accept())
            << "                          |";
     stream << "\n|";
     stream << setfill('-') << setw(115) << "|";
   }
 
   if(printlevel>=2) {
     stream << "\n| ";
     stream << std::scientific << setprecision(5);
 
     stream << "sig(Ev) = "
            << setfill(' ') << setw(17) << fXSec/units::cm2
            << " cm^2  |";
 
     switch(fDiffXSecPhSp) {
       case ( kPSyfE   ) :
         stream << " dsig(y;E)/dy =          " << setfill(' ') << setw(13) << fDiffXSec/units::cm2 << " cm^2       |";
         break;
       case ( kPSxyfE  ) :
         stream << " d2sig(x,y;E)/dxdy =     " << setfill(' ') << setw(13) << fDiffXSec/units::cm2 << " cm^2       |";
         break;
       case ( kPSxytfE ) :
         stream << " d3sig(x,y,t;E)/dxdydt = " << setfill(' ') << setw(13) << fDiffXSec/units::cm2 << " cm^2/GeV^2 |";
         break;
       case ( kPSQ2fE  ) :
         stream << " dsig(Q2;E)/dQ2 =        " << setfill(' ') << setw(13) << fDiffXSec/units::cm2 << " cm^2/GeV^2 |";
         break;
       case ( kPSQ2vpfE  ) :
         stream << " dsig(Q2,v,p;E)/dQ2dvdp =" << setfill(' ') << setw(13) << fDiffXSec/units::cm2 << " cm^2/GeV^4 |";
         break;
       case ( kPSWQ2fE ) :
         stream << " d2sig(W,Q2;E)/dWdQ2 =   " << setfill(' ') << setw(13) << fDiffXSec/units::cm2 << " cm^2/GeV^3 |";
         break;
       case ( kPSWQ2ctpphipfE ) :
         stream << " d4sig(W,Q2,cos(theta),phi;E)/dWdQ2dOmega =   " << setfill(' ') << setw(13) << fDiffXSec/units::cm2 << " cm^2/GeV^3 |";
         break;
       default :
         stream << " dsig(Ev;{K_s})/dK   =   " << setfill(' ') << setw(13) << fDiffXSec/units::cm2 << " cm^2/{K}   |";
     }
     stream << " Weight = "
            << setfill(' ') << setw(16)
            << std::fixed << setprecision(5)
            << fWeight
            << " |";
 
     stream << "\n|";
     stream << setfill('-') << setw(115) << "|";
   }
 
   stream << "\n";
   stream << setfill(' ');
 
   if(printlevel==3) {
      if(fInteraction) stream << *fInteraction;
      else             stream << "NULL Interaction!" << endl;
   }
   stream << "\n";
 }

virtual double genie::GHepRecord::Probability ( void ) const

inlinevirtual

Definition at line 125 of file GHepRecord.h.

References fProb.

125 { return fProb; }

genie::GHepRecord::fProb

double fProb

event probability (for given flux neutrino and density-weighted path-length for target element) ...

Definition: GHepRecord.h:179

GHepParticle * GHepRecord::Probe ( void ) const

virtual

Definition at line 284 of file GHepRecord.cxx.

References Particle(), and ProbePosition().

 {
 // Returns the GHepParticle representing the probe (neutrino, e,...).
 
   int ipos = this->ProbePosition();
   if(ipos>-1) return this->Particle(ipos);
   return 0;
 }

int GHepRecord::ProbePosition ( void ) const

virtual

Definition at line 352 of file GHepRecord.cxx.

References EventGenerationMode(), genie::kGMdDarkMatterNucleus, genie::kGMdHadronNucleus, genie::kGMdHNLDecay, genie::kGMdLeptonNucleus, and genie::kGMdPhotonNucleus.

Referenced by genie::COHDNuEventGenerator::AddFinalStateDarkNeutrino(), genie::HEDISGenerator::AddPrimaryLepton(), genie::PrimaryLeptonGenerator::AddToEventRecord(), genie::OutgoingDarkGenerator::AddToEventRecord(), Probe(), genie::QELEventGenerator::ProcessEventRecord(), genie::NormGenerator::ProcessEventRecord(), genie::DMELEventGenerator::ProcessEventRecord(), genie::QELEventGeneratorSM::ProcessEventRecord(), and genie::SPPEventGenerator::ProcessEventRecord().

 {
 // Returns the GHEP position of the GHepParticle representing the probe
 // (neutrino, e,...).
 
   // The probe is *always* at slot 0.
   GEvGenMode_t mode = this->EventGenerationMode();
   if(mode == kGMdLeptonNucleus ||
      mode == kGMdDarkMatterNucleus ||
      mode == kGMdHadronNucleus ||
      mode == kGMdPhotonNucleus ||
      mode == kGMdHNLDecay)
   {
     return 0;
   }
   return -1;
 }

GHepParticle * GHepRecord::RemnantNucleus ( void ) const

virtual

Definition at line 303 of file GHepRecord.cxx.

References Particle(), and RemnantNucleusPosition().

 {
 // Returns the GHepParticle representing the remnant nucleus,
 // or 0 if it does not exist.
 
   int ipos = this->RemnantNucleusPosition();
   if(ipos>-1) return this->Particle(ipos);
   return 0;
 }

int GHepRecord::RemnantNucleusPosition ( void ) const

virtual

Definition at line 397 of file GHepRecord.cxx.

References genie::GHepParticle::FirstDaughter(), genie::pdg::IsIon(), genie::kIStStableFinalState, genie::GHepParticle::LastDaughter(), Particle(), genie::GHepParticle::Pdg(), genie::GHepParticle::Status(), and TargetNucleus().

Referenced by genie::HEDISGenerator::ProcessEventRecord(), genie::NucBindEnergyAggregator::ProcessEventRecord(), RemnantNucleus(), genie::Intranuke::TransportHadrons(), and genie::Intranuke2018::TransportHadrons().

 {
 // Returns the GHEP position of the GHepParticle representing the remnant
 // nucleus - or -1 if the interaction takes place at a free nucleon.
 
   GHepParticle * p = this->TargetNucleus();
   if(!p) return -1;
 
   int dau1 = p->FirstDaughter();
   int dau2 = p->LastDaughter();
 
   if(dau1==-1 && dau2==-1) return -1;
 
   for(int i=dau1; i<=dau2; i++) {
     GHepParticle * dp = this->Particle(i);
     int dpdgc = dp->Pdg();
     if(pdg::IsIon(dpdgc) && dp->Status()==kIStStableFinalState) return i;
   }
   return -1;
 }

void GHepRecord::RemoveIntermediateParticles ( void )

virtual

Definition at line 617 of file GHepRecord.cxx.

References genie::kIStInitialState, genie::kIStNucleonTarget, genie::kIStStableFinalState, LOG, genie::GHepParticle::Name(), pDEBUG, pNOTICE, genie::GHepParticle::SetFirstDaughter(), genie::GHepParticle::SetFirstMother(), genie::GHepParticle::SetLastDaughter(), genie::GHepParticle::SetLastMother(), and genie::GHepParticle::Status().

 {
   LOG("GHEP", pNOTICE) << "Removing all intermediate particles from GHEP";
   this->Compress();
 
   int i=0;
   GHepParticle * p = 0;
 
   TIter iter(this);
   while( (p = (GHepParticle *)iter.Next()) ) {
 
     if(!p) continue;
     GHepStatus_t ist = p->Status();
 
     bool keep = (ist==kIStInitialState) ||
                 (ist==kIStStableFinalState) || (ist==kIStNucleonTarget);
     if(keep) {
        p->SetFirstDaughter(-1);
        p->SetLastDaughter(-1);
        p->SetFirstMother(-1);
        p->SetLastMother(-1);
     } else {
        LOG("GHEP", pNOTICE)
            << "Removing: " << p->Name() << " from slot: " << i;
        this->RemoveAt(i);
     }
     i++;
   }
 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
   LOG("GHEP", pDEBUG) << "Compressing GHEP record to remove empty slots";
 #endif
   this->Compress();
 }

void GHepRecord::ResetRecord ( void )

virtual

Definition at line 874 of file GHepRecord.cxx.

References CleanRecord(), InitRecord(), LOG, and pDEBUG.

Referenced by Copy(), and genie::EventGenerator::ProcessEventRecord().

 {
 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
   LOG("GHEP", pDEBUG) << "Reseting GHepRecord";
 #endif
   this->CleanRecord();
   this->InitRecord();
 }

virtual void genie::GHepRecord::SetDiffXSec	(	double	xsec,
		KinePhaseSpace_t	ps
	)

inlinevirtual

Definition at line 133 of file GHepRecord.h.

References fDiffXSec, fDiffXSecPhSp, and genie::units::ps.

   { fDiffXSecPhSp = ps;
     fDiffXSec = (xsec>0) ? xsec : 0.;
   }

void GHepRecord::SetPrintLevel ( int print_level )

static

Definition at line 956 of file GHepRecord.cxx.

References fPrintLevel.

Referenced by main(), and TestDecay().

 {
   fPrintLevel = print_level;
 }

virtual void genie::GHepRecord::SetProbability ( double prob )

inlinevirtual

Definition at line 131 of file GHepRecord.h.

References fProb.

Referenced by main().

131 { fProb = (prob>0) ? prob : 0.; }

genie::GHepRecord::fProb

double fProb

event probability (for given flux neutrino and density-weighted path-length for target element) ...

Definition: GHepRecord.h:179

void GHepRecord::SetUnphysEventMask ( const TBits & mask )

Definition at line 938 of file GHepRecord.cxx.

References fEventMask, LOG, genie::GHepFlags::NFlags(), and pINFO.

 {
  *fEventMask = mask;
 
   LOG("GHEP", pINFO)
     << "Setting unphysical event mask (bits: " << GHepFlags::NFlags() - 1
     << " -> 0) : " << *fEventMask;
 }

void GHepRecord::SetVertex	(	double	x,
		double	y,
		double	z,
		double	t
	)

virtual

Definition at line 827 of file GHepRecord.cxx.

References fVtx.

Referenced by genie::hnl::VertexGenerator::ProcessEventRecord(), and genie::hnl::FluxCreator::ProcessEventRecord().

 {
   fVtx->SetXYZT(x,y,z,t);
 }

void GHepRecord::SetVertex ( const TLorentzVector & vtx )

virtual

Definition at line 832 of file GHepRecord.cxx.

References fVtx.

 {
   fVtx->SetXYZT(vtx.X(),vtx.Y(),vtx.Z(),vtx.T());
 }

virtual void genie::GHepRecord::SetWeight ( double wght )

inlinevirtual

Definition at line 130 of file GHepRecord.h.

References fWeight.

130 { fWeight = (wght>0) ? wght : 0.; }

genie::GHepRecord::fWeight

double fWeight

event weight

Definition: GHepRecord.h:178

virtual void genie::GHepRecord::SetXSec ( double xsec )

inlinevirtual

Definition at line 132 of file GHepRecord.h.

References fXSec.

Referenced by genie::hnl::FluxCreator::ProcessEventRecord(), and genie::PhysInteractionSelector::SelectInteraction().

132 { fXSec = (xsec>0) ? xsec : 0.; }

genie::GHepRecord::fXSec

double fXSec

cross section for selected event

Definition: GHepRecord.h:180

Interaction * GHepRecord::Summary ( void ) const

virtual

Definition at line 91 of file GHepRecord.cxx.

References fInteraction, LOG, and pWARN.

 {
   if(!fInteraction) {
     LOG("GHEP", pWARN) << "Returning NULL interaction";
   }
   return fInteraction;
 }

void GHepRecord::SwapParticles	(	int	i,
		int	j
	)

protectedvirtual

Definition at line 752 of file GHepRecord.cxx.

References genie::GHepParticle::Copy(), genie::GHepParticle::HasDaughters(), LOG, genie::GHepParticle::Name(), Particle(), pINFO, and genie::GHepParticle::SetFirstMother().

Referenced by CompactifyDaughterLists().

 {
 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
   LOG("GHEP", pINFO) << "Swapping GHepParticles : " << i << " <--> " << j;
 #endif
   int n = this->GetEntries();
   assert(i>=0 && j>=0 && i<n && j<n);
 
   if(i==j) return;
 
   GHepParticle * pi  = this->Particle(i);
   GHepParticle * pj  = this->Particle(j);
   GHepParticle * tmp = new GHepParticle(*pi);
 
   pi->Copy(*pj);
   pj->Copy(*tmp);
 
   delete tmp;
 
   // tell their daughters
   if(pi->HasDaughters()) {
 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
     LOG("GHEP", pINFO)
       << pi->Name() << "(previously at pos: " << j
       << ") is now at pos: " << i << " -> Notify daughters";
 #endif
     for(int k=0; k<n; k++) {
       if(this->Particle(k)->FirstMother()==j) {
         this->Particle(k)->SetFirstMother(i);
       }
     }
   }
 
   if(pj->HasDaughters()) {
 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
     LOG("GHEP", pINFO)
       << pj->Name() << "(previously at pos: " << i
       << ") is now at pos: " << j << " -> Notify daughters";
 #endif
     for(int k=0; k<n; k++) {
       if(this->Particle(k)->FirstMother()==i) {
         this->Particle(k)->SetFirstMother(j);
       }
     }
   }
 }

GHepParticle * GHepRecord::TargetNucleus ( void ) const

virtual

Definition at line 293 of file GHepRecord.cxx.

References Particle(), and TargetNucleusPosition().

 {
 // Returns the GHepParticle representing the target / initial state nucleus,
 // or 0 if it does not exist.
 
   int ipos = this->TargetNucleusPosition();
   if(ipos>-1) return this->Particle(ipos);
   return 0;
 }

int GHepRecord::TargetNucleusPosition ( void ) const

virtual

Definition at line 370 of file GHepRecord.cxx.

References EventGenerationMode(), genie::pdg::IsIon(), genie::kGMdDarkMatterNucleus, genie::kGMdHadronNucleus, genie::kGMdLeptonNucleus, genie::kGMdNucleonDecay, genie::kGMdPhotonNucleus, genie::kIStInitialState, Particle(), genie::GHepParticle::Pdg(), and genie::GHepParticle::Status().

Referenced by genie::NuETargetRemnantGenerator::AddTargetNucleusRemnant(), genie::DMETargetRemnantGenerator::AddTargetNucleusRemnant(), genie::HadronicSystemGenerator::AddTargetNucleusRemnant(), genie::QELEventGenerator::AddTargetNucleusRemnant(), genie::FermiMover::AddTargetNucleusRemnant(), genie::DMELEventGenerator::AddTargetNucleusRemnant(), genie::QELEventGeneratorSM::AddTargetNucleusRemnant(), genie::AMNuGammaGenerator::AddTargetRemnant(), genie::COHHadronicSystemGenerator::CalculateHadronicSystem_AlvarezRuso(), genie::COHHadronicSystemGenerator::CalculateHadronicSystem_BergerSehgalFM(), genie::COHHadronicSystemGenerator::CalculateHadronicSystem_ReinSehgal(), genie::SecondNucleonEmissionI::EmitSecondNucleon(), genie::NormGenerator::ProcessEventRecord(), TargetNucleus(), genie::Intranuke::TransportHadrons(), and genie::Intranuke2018::TransportHadrons().

 {
 // Returns the GHEP position of the GHepParticle representing the target
 // nucleus - or -1 if the interaction takes place at a free nucleon.
 
   GEvGenMode_t mode = this->EventGenerationMode();
 
   if(mode == kGMdLeptonNucleus ||
      mode == kGMdDarkMatterNucleus ||
      mode == kGMdHadronNucleus ||
      mode == kGMdPhotonNucleus)
   {
      GHepParticle * p = this->Particle(1); // If exists, it will be at slot 1
      if(!p) return -1;
      int pdgc = p->Pdg();
      if(pdg::IsIon(pdgc) && p->Status()==kIStInitialState) return 1;
   }
   if(mode == kGMdNucleonDecay) {
      GHepParticle * p = this->Particle(0); // If exists, it will be at slot 0
      if(!p) return -1;
      int pdgc = p->Pdg();
      if(pdg::IsIon(pdgc) && p->Status()==kIStInitialState) return 0;
   }
 
   return -1;
 }

void GHepRecord::UpdateDaughterLists ( void )

protectedvirtual

Definition at line 556 of file GHepRecord.cxx.

References CompactifyDaughterLists(), genie::GHepParticle::FirstDaughter(), genie::GHepParticle::FirstMother(), genie::GHepParticle::LastDaughter(), LOG, Particle(), pINFO, pNOTICE, genie::GHepParticle::SetFirstDaughter(), and genie::GHepParticle::SetLastDaughter().

Referenced by AddParticle().

 {
   int pos = this->GetEntries() - 1; // position of last entry
 
 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
   LOG("GHEP", pINFO)
      << "Updating the daughter-list for the mother of particle at: " << pos;
 #endif
 
   GHepParticle * p = this->Particle(pos);
   assert(p);
 
   int mom_pos = p->FirstMother();
 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
   LOG("GHEP", pINFO) << "Mother particle is at slot: " << mom_pos;
 #endif
   if(mom_pos==-1) return; // may not have mom (eg init state)
   GHepParticle * mom = this->Particle(mom_pos);
   if(!mom) return; // may not have mom (eg init state)
 
   int dau1 = mom->FirstDaughter();
   int dau2 = mom->LastDaughter();
 
   // handles the case where the daughter list was initially empty
   if(dau1 == -1) {
      mom->SetFirstDaughter(pos);
      mom->SetLastDaughter(pos);
  #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
      LOG("GHEP", pINFO)
         << "Done! Daughter-list is compact: [" << pos << ", " << pos << "]";
  #endif
      return;
   }
   // handles the case where the new daughter is added at the slot just before
   // an already compact daughter list
   if(pos == dau1-1) {
      mom->SetFirstDaughter(pos);
  #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
      LOG("GHEP", pINFO)
        << "Done! Daughter-list is compact: [" << pos << ", " << dau2 << "]";
  #endif
      return;
   }
   // handles the case where the new daughter is added at the slot just after
   // an already compact daughter list
   if(pos == dau2+1) {
      mom->SetLastDaughter(pos);
  #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
      LOG("GHEP", pINFO)
        << "Done! Daughter-list is compact: [" << dau1 << ", " << pos << "]";
  #endif
      return;
   }
 
   // If you are here, then the last particle insertion broke the daughter
   // list compactification - Run the compactifier
   LOG("GHEP", pNOTICE)
       << "Daughter-list is not compact - Running compactifier";
   this->CompactifyDaughterLists();
 }

virtual TLorentzVector* genie::GHepRecord::Vertex ( void ) const

inlinevirtual

Definition at line 140 of file GHepRecord.h.

References fVtx.

Referenced by genie::hnl::Decayer::AddInitialState(), Copy(), Print(), and genie::hnl::VertexGenerator::ProcessEventRecord().

140 { return fVtx; }

genie::GHepRecord::fVtx

TLorentzVector * fVtx

vertex in the detector coordinate system

Definition: GHepRecord.h:171

virtual double genie::GHepRecord::Weight ( void ) const

inlinevirtual

Definition at line 124 of file GHepRecord.h.

References fWeight.

124 { return fWeight; }

genie::GHepRecord::fWeight

double fWeight

event weight

Definition: GHepRecord.h:178

virtual double genie::GHepRecord::XSec ( void ) const

inlinevirtual

Definition at line 126 of file GHepRecord.h.

References fXSec.

126 { return fXSec; }

genie::GHepRecord::fXSec

double fXSec

cross section for selected event

Definition: GHepRecord.h:180

Friends And Related Function Documentation

ostream& operator<<	(	ostream &	stream,
		const GHepRecord &	event
	)

friend

Definition at line 46 of file GHepRecord.cxx.

  {
    rec.Print(stream);
    return stream;
  }

Member Data Documentation

double genie::GHepRecord::fDiffXSec

protected

differential cross section for selected event kinematics

Definition at line 181 of file GHepRecord.h.

Referenced by Copy(), DiffXSec(), InitRecord(), Print(), and SetDiffXSec().

KinePhaseSpace_t genie::GHepRecord::fDiffXSecPhSp

protected

specifies which differential cross-section (dsig/dQ2, dsig/dQ2dW, dsig/dxdy,...)

Definition at line 182 of file GHepRecord.h.

Referenced by Copy(), DiffXSecVars(), InitRecord(), Print(), and SetDiffXSec().

TBits* genie::GHepRecord::fEventFlags

protected

event flags indicating various pathologies or an unphysical event

Definition at line 174 of file GHepRecord.h.

Referenced by Accept(), Clear(), Copy(), EventFlags(), InitRecord(), IsUnphysical(), and Print().

TBits* genie::GHepRecord::fEventMask

protected

an input bit-field mask allowing one to ignore bits set in fEventFlags

Definition at line 175 of file GHepRecord.h.

Referenced by Accept(), Clear(), Copy(), EventMask(), InitRecord(), Print(), and SetUnphysEventMask().

Interaction* genie::GHepRecord::fInteraction

protected

attached summary information

Definition at line 168 of file GHepRecord.h.

Referenced by AttachSummary(), Clear(), Copy(), InitRecord(), Print(), and Summary().

int GHepRecord::fPrintLevel = 3

staticprotected

Definition at line 196 of file GHepRecord.h.

Referenced by GetPrintLevel(), Print(), and SetPrintLevel().

double genie::GHepRecord::fProb

protected

event probability (for given flux neutrino and density-weighted path-length for target element)

Definition at line 179 of file GHepRecord.h.

Referenced by Copy(), InitRecord(), Probability(), and SetProbability().

TLorentzVector* genie::GHepRecord::fVtx

protected

vertex in the detector coordinate system

Definition at line 171 of file GHepRecord.h.

Referenced by Clear(), Copy(), InitRecord(), SetVertex(), and Vertex().

double genie::GHepRecord::fWeight

protected

event weight

Definition at line 178 of file GHepRecord.h.

Referenced by Copy(), InitRecord(), Print(), SetWeight(), and Weight().

double genie::GHepRecord::fXSec

protected

cross section for selected event

Definition at line 180 of file GHepRecord.h.

Referenced by Copy(), InitRecord(), Print(), SetXSec(), and XSec().

The documentation for this class was generated from the following files:

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Static Protected Attributes

Friends

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Friends And Related Function Documentation

Member Data Documentation