genie::GHepRecord Class Reference

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
TLorentzVector *	fVtx
	vertex in the detector coordinate system
TBits *	fEventFlags
	event flags indicating various pathologies or an unphysical event
TBits *	fEventMask
	an input bit-field mask allowing one to ignore bits set in fEventFlags
double	fWeight
	event weight
double	fProb
	event probability (for given flux neutrino and density-weighted path-length for target element)
double	fXSec
	cross section for selected event
double	fDiffXSec
	differential cross section for selected event kinematics
KinePhaseSpace_t	fDiffXSecPhSp
	specifies which differential cross-section (dsig/dQ2, dsig/dQ2dW, dsig/dxdy,...)

Static Protected Attributes
static int	fPrintLevel

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

Detailed Description

GENIE's GHEP MC event record.

Author

Costas Andreopoulos <c.andreopoulos \at cern.ch> University of Liverpool

Created:\n October 1, 2004

License:\n Copyright (c) 2003-2025, 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() [1/4]

GHepRecord::GHepRecord

(

)

Definition at line 53 of file GHepRecord.cxx.

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

References InitRecord().

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

GHepRecord() [2/4]

GHepRecord::GHepRecord

(

int

size

)

Definition at line 59 of file GHepRecord.cxx.

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

References InitRecord().

GHepRecord() [3/4]

GHepRecord::GHepRecord

(

const GHepRecord &

record

)

Definition at line 65 of file GHepRecord.cxx.

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

References Copy(), GHepRecord(), and InitRecord().

GHepRecord() [4/4]

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

References fDiffXSec, fEventFlags, fEventMask, fInteraction, fProb, fVtx, fWeight, and fXSec.

~GHepRecord()

GHepRecord::~GHepRecord ( )

virtual

Definition at line 86 of file GHepRecord.cxx.

{
  this->CleanRecord();
}

References CleanRecord().

Member Function Documentation

Accept()

bool GHepRecord::Accept ( void ) const

virtual

Definition at line 947 of file GHepRecord.cxx.

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

References fEventFlags, and fEventMask.

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

AddParticle() [1/3]

void GHepRecord::AddParticle ( const GHepParticle & p )

virtual

Definition at line 499 of file GHepRecord.cxx.

{
// 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();
}

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

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

AddParticle() [2/3]

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.

{
// 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();
}

References LOG, pINFO, and UpdateDaughterLists().

AddParticle() [3/3]

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.

{
// 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();
}

References LOG, pINFO, and UpdateDaughterLists().

AttachSummary()

void GHepRecord::AttachSummary ( Interaction * interaction )

virtual

Definition at line 99 of file GHepRecord.cxx.

{
  fInteraction = interaction;
}

References fInteraction.

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

CleanRecord()

void GHepRecord::CleanRecord ( void )

protected

Definition at line 866 of file GHepRecord.cxx.

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

References Clear(), LOG, and pDEBUG.

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

Clear()

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

virtual

Definition at line 883 of file GHepRecord.cxx.

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

References fEventFlags, fEventMask, fInteraction, and fVtx.

Referenced by CleanRecord().

CompactifyDaughterLists()

void GHepRecord::CompactifyDaughterLists ( void )

virtual

compact

Definition at line 651 of file GHepRecord.cxx.

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

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().

Copy()

void GHepRecord::Copy ( const GHepRecord & record )

virtual

Reimplemented in genie::EventRecord.

Definition at line 907 of file GHepRecord.cxx.

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

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

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

DiffXSec()

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

inlinevirtual

Definition at line 127 of file GHepRecord.h.

{ return fDiffXSec; }

References fDiffXSec.

DiffXSecVars()

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

inlinevirtual

Definition at line 128 of file GHepRecord.h.

{ return fDiffXSecPhSp; }

References fDiffXSecPhSp.

EventFlags()

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

inlinevirtual

Definition at line 117 of file GHepRecord.h.

{ return fEventFlags; }

References fEventFlags.

Referenced by genie::SKHadronicSystemGenerator::CalculateHadronicSystem_AtharSingleKaon(), genie::SKKinematicsGenerator::CalculateKin_AtharSingleKaon(), Copy(), genie::FermiMover::KickHitNucleon(), genie::KineGeneratorWithCache::MaxXSec(), genie::DFRHadronicSystemGenerator::ProcessEventRecord(), genie::DFRKinematicsGenerator::ProcessEventRecord(), genie::DFRPrimaryLeptonGenerator::ProcessEventRecord(), genie::DISKinematicsGenerator::ProcessEventRecord(), genie::DISPrimaryLeptonGenerator::ProcessEventRecord(), genie::DMDISKinematicsGenerator::ProcessEventRecord(), genie::DMDISOutgoingDarkGenerator::ProcessEventRecord(), genie::DMEKinematicsGenerator::ProcessEventRecord(), genie::DMELEventGenerator::ProcessEventRecord(), genie::DMELKinematicsGenerator::ProcessEventRecord(), genie::HEDISKinematicsGenerator::ProcessEventRecord(), genie::HELeptonKinematicsGenerator::ProcessEventRecord(), genie::IBDKinematicsGenerator::ProcessEventRecord(), genie::NuEKinematicsGenerator::ProcessEventRecord(), genie::PauliBlocker::ProcessEventRecord(), genie::QELEventGenerator::ProcessEventRecord(), genie::QELEventGeneratorSM::ProcessEventRecord(), genie::QELKinematicsGenerator::ProcessEventRecord(), genie::RESKinematicsGenerator::ProcessEventRecord(), genie::SKPrimaryLeptonGenerator::ProcessEventRecord(), genie::SPPEventGenerator::ProcessEventRecord(), genie::DMELKinematicsGenerator::SpectralFuncExperimentalCode(), genie::QELKinematicsGenerator::SpectralFuncExperimentalCode(), and genie::COHKinematicsGenerator::throwOnTooManyIterations().

EventGenerationMode()

GEvGenMode_t GHepRecord::EventGenerationMode

(

void

)

const

Definition at line 209 of file GHepRecord.cxx.

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

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::Intranuke2025::ProcessEventRecord(), genie::Intranuke::ProcessEventRecord(), and TargetNucleusPosition().

EventMask()

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

inlinevirtual

Definition at line 118 of file GHepRecord.h.

{ return fEventMask; }

References fEventMask.

Referenced by Copy().

FinalizeDaughterLists()

void GHepRecord::FinalizeDaughterLists ( void )

protectedvirtual

Definition at line 799 of file GHepRecord.cxx.

{
// 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);
  }
}

References genie::GHepParticle::FirstMother().

Referenced by CompactifyDaughterLists().

FinalStateHadronicSystem()

GHepParticle * GHepRecord::FinalStateHadronicSystem ( void ) const

virtual

Definition at line 342 of file GHepRecord.cxx.

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

References FinalStateHadronicSystemPosition(), and Particle().

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

FinalStateHadronicSystemPosition()

int GHepRecord::FinalStateHadronicSystemPosition ( void ) const

virtual

Definition at line 471 of file GHepRecord.cxx.

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

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

Referenced by FinalStateHadronicSystem().

FinalStatePrimaryLepton()

GHepParticle * GHepRecord::FinalStatePrimaryLepton ( void ) const

virtual

Definition at line 333 of file GHepRecord.cxx.

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

References FinalStatePrimaryLeptonPosition(), and Particle().

Referenced by genie::DMETargetRemnantGenerator::AddElectronNeutrino(), genie::NuETargetRemnantGenerator::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::DMDISOutgoingDarkGenerator::ProcessEventRecord(), genie::SKPrimaryLeptonGenerator::ProcessEventRecord(), genie::OutgoingDarkGenerator::SetPolarization(), genie::utils::SetPrimaryLeptonPolarization(), and genie::hnl::Decayer::UpdateEventRecord().

FinalStatePrimaryLeptonPosition()

int GHepRecord::FinalStatePrimaryLeptonPosition ( void ) const

virtual

Definition at line 459 of file GHepRecord.cxx.

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

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

Referenced by FinalStatePrimaryLepton().

FindParticle()

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

virtual

Definition at line 118 of file GHepRecord.cxx.

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

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

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

FirstNonInitStateEntry()

int GHepRecord::FirstNonInitStateEntry ( void )

protectedvirtual

Definition at line 739 of file GHepRecord.cxx.

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

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

GetPrintLevel()

int GHepRecord::GetPrintLevel ( )

static

Definition at line 960 of file GHepRecord.cxx.

{
  return fPrintLevel;
}

References fPrintLevel.

GetStableDescendants()

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

virtual

Definition at line 174 of file GHepRecord.cxx.

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

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

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

HasCompactDaughterList()

bool GHepRecord::HasCompactDaughterList ( int pos )

protectedvirtual

Definition at line 700 of file GHepRecord.cxx.

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

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

Referenced by CompactifyDaughterLists().

HitElectron()

GHepParticle * GHepRecord::HitElectron ( void ) const

virtual

Definition at line 323 of file GHepRecord.cxx.

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

References HitElectronPosition(), and Particle().

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

HitElectronPosition()

int GHepRecord::HitElectronPosition ( void ) const

virtual

Definition at line 441 of file GHepRecord.cxx.

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

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

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

HitNucleon()

GHepParticle * GHepRecord::HitNucleon ( void ) const

virtual

Definition at line 313 of file GHepRecord.cxx.

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

References HitNucleonPosition(), and Particle().

Referenced by genie::AMNuGammaGenerator::AddPhoton(), genie::AMNuGammaGenerator::AddRecoilNucleon(), genie::HadronicSystemGenerator::AddTargetNucleusRemnant(), genie::AMNuGammaGenerator::AddTargetRemnant(), genie::NucDeExcitationSim::ArgonTargetSim(), genie::SKHadronicSystemGenerator::CalculateHadronicSystem_AtharSingleKaon(), genie::NucDeExcitationSim::CarbonTargetSim(), genie::SecondNucleonEmissionI::EmitSecondNucleon(), genie::HadronicSystemGenerator::Hadronic4pLAB(), genie::FermiMover::KickHitNucleon(), genie::NucDeExcitationSim::OxygenTargetSim(), genie::DFRHadronicSystemGenerator::ProcessEventRecord(), genie::DMELEventGenerator::ProcessEventRecord(), genie::DMELKinematicsGenerator::ProcessEventRecord(), genie::FermiMover::ProcessEventRecord(), genie::PauliBlocker::ProcessEventRecord(), genie::QELEventGenerator::ProcessEventRecord(), genie::QELEventGeneratorSM::ProcessEventRecord(), genie::QELKinematicsGenerator::ProcessEventRecord(), genie::SpectralFunction2p2h::ProcessEventRecord(), genie::SPPEventGenerator::ProcessEventRecord(), genie::SRCNuclearRecoil::ProcessEventRecord(), genie::DMELKinematicsGenerator::SpectralFuncExperimentalCode(), and genie::QELKinematicsGenerator::SpectralFuncExperimentalCode().

HitNucleonPosition()

int GHepRecord::HitNucleonPosition ( void ) const

virtual

Definition at line 418 of file GHepRecord.cxx.

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

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::IBDHadronicSystemGenerator::AddRecoilBaryon(), genie::QELHadronicSystemGenerator::AddRecoilBaryon(), genie::AMNuGammaGenerator::AddRecoilNucleon(), genie::RESHadronicSystemGenerator::AddResonance(), genie::RSPPResonanceSelector::AddResonance(), genie::SKHadronicSystemGenerator::CalculateHadronicSystem_AtharSingleKaon(), genie::QELEventGeneratorSuSA::GenerateNucleon(), HitNucleon(), genie::DFRHadronicSystemGenerator::ProcessEventRecord(), genie::DMELEventGenerator::ProcessEventRecord(), genie::QELEventGenerator::ProcessEventRecord(), genie::QELEventGeneratorSM::ProcessEventRecord(), and genie::SPPEventGenerator::ProcessEventRecord().

InitRecord()

void GHepRecord::InitRecord ( void )

protected

Definition at line 837 of file GHepRecord.cxx.

{
#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);
}

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

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

IsUnphysical()

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

inlinevirtual

Definition at line 119 of file GHepRecord.h.

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

References fEventFlags.

Referenced by GenerateEvent(), and Print().

NEntries() [1/2]

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

virtual

Definition at line 477 of file GHepRecord.cxx.

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

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

NEntries() [2/2]

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

virtual

Definition at line 488 of file GHepRecord.cxx.

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

References genie::GHepParticle::Pdg().

Particle()

GHepParticle * GHepRecord::Particle ( int position ) const

virtual

Definition at line 104 of file GHepRecord.cxx.

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

References LOG, and pINFO.

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

ParticlePosition() [1/2]

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

virtual

Definition at line 156 of file GHepRecord.cxx.

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

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

ParticlePosition() [2/2]

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

virtual

Definition at line 137 of file GHepRecord.cxx.

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

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

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

Print()

void GHepRecord::Print

(

ostream &

stream

)

const

Definition at line 965 of file GHepRecord.cxx.

{
  // 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";
}

References Accept(), genie::utils::print::BoolAsYNString(), genie::units::cm2, genie::GHepFlags::Describe(), genie::GHepParticle::E(), fDiffXSec, fDiffXSecPhSp, fEventFlags, fEventMask, fInteraction, genie::GHepParticle::FirstDaughter(), genie::GHepParticle::FirstMother(), fPrintLevel, fWeight, fXSec, genie::GHepParticle::GetPolarization(), genie::GHepParticle::IsOnMassShell(), IsUnphysical(), genie::kIStFinalStateNuclearRemnant, genie::kIStInitialState, genie::kIStStableFinalState, genie::kPSQ2fE, genie::kPSQ2vpfE, genie::kPSWQ2ctpphipfE, genie::kPSWQ2fE, genie::kPSxyfE, genie::kPSxytfE, genie::kPSyfE, genie::GHepParticle::LastDaughter(), genie::GHepParticle::LastMother(), genie::GHepParticle::Mass(), genie::GHepParticle::Name(), genie::GHepParticle::P4(), genie::GHepParticle::Pdg(), genie::GHepParticle::PolzIsSet(), Probe(), genie::GHepParticle::Px(), genie::GHepParticle::Py(), genie::GHepParticle::Pz(), genie::GHepParticle::RescatterCode(), genie::GHepParticle::Status(), Vertex(), genie::GHepParticle::Vt(), genie::GHepParticle::Vx(), genie::GHepParticle::Vy(), and genie::GHepParticle::Vz().

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

Probability()

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

inlinevirtual

Definition at line 125 of file GHepRecord.h.

{ return fProb;     }

References fProb.

Probe()

GHepParticle * GHepRecord::Probe ( void ) const

virtual

Definition at line 284 of file GHepRecord.cxx.

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

References Particle(), and ProbePosition().

Referenced by genie::DMETargetRemnantGenerator::AddElectronNeutrino(), genie::NuETargetRemnantGenerator::AddElectronNeutrino(), genie::MECGenerator::AddFinalStateLepton(), genie::AMNuGammaGenerator::AddFinalStateNeutrino(), genie::AMNuGammaGenerator::AddPhoton(), genie::HEDISGenerator::AddPrimaryLepton(), genie::IBDHadronicSystemGenerator::AddRecoilBaryon(), genie::QELHadronicSystemGenerator::AddRecoilBaryon(), genie::AMNuGammaGenerator::AddRecoilNucleon(), genie::RESHadronicSystemGenerator::AddResonance(), genie::OutgoingDarkGenerator::AddToEventRecord(), genie::PrimaryLeptonGenerator::AddToEventRecord(), genie::COHHadronicSystemGenerator::CalculateHadronicSystem_AlvarezRuso(), genie::SKHadronicSystemGenerator::CalculateHadronicSystem_AtharSingleKaon(), genie::COHHadronicSystemGenerator::CalculateHadronicSystem_BergerSehgalFM(), genie::COHHadronicSystemGenerator::CalculateHadronicSystem_ReinSehgal(), FinalStatePrimaryLeptonPosition(), FindhAFate(), genie::hnl::Decayer::GenerateDecayProducts(), genie::Intranuke2018::GenerateVertex(), genie::Intranuke2025::GenerateVertex(), genie::Intranuke::GenerateVertex(), genie::HadronicSystemGenerator::Hadronic4pLAB(), genie::HAIntranuke2018::InelasticHA(), genie::HAIntranuke2025::InelasticHA(), genie::HAIntranuke::InelasticHA(), genie::HadronicSystemGenerator::MomentumTransferLAB(), Print(), genie::DFRHadronicSystemGenerator::ProcessEventRecord(), genie::DMELEventGenerator::ProcessEventRecord(), genie::DMEOutgoingDarkGenerator::ProcessEventRecord(), genie::NormGenerator::ProcessEventRecord(), genie::NuEPrimaryLeptonGenerator::ProcessEventRecord(), genie::OutgoingDarkGenerator::ProcessEventRecord(), genie::PrimaryLeptonGenerator::ProcessEventRecord(), genie::QELEventGenerator::ProcessEventRecord(), genie::QELEventGeneratorSM::ProcessEventRecord(), genie::SPPEventGenerator::ProcessEventRecord(), genie::QELEventGeneratorSuSA::SelectLeptonKinematics(), genie::MECGenerator::SelectNSVLeptonKinematics(), and genie::MECGenerator::SelectSuSALeptonKinematics().

ProbePosition()

int GHepRecord::ProbePosition ( void ) const

virtual

Definition at line 352 of file GHepRecord.cxx.

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

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

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

RemnantNucleus()

GHepParticle * GHepRecord::RemnantNucleus ( void ) const

virtual

Definition at line 303 of file GHepRecord.cxx.

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

References Particle(), and RemnantNucleusPosition().

RemnantNucleusPosition()

int GHepRecord::RemnantNucleusPosition ( void ) const

virtual

Definition at line 397 of file GHepRecord.cxx.

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

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::Intranuke2018::TransportHadrons(), genie::Intranuke2025::TransportHadrons(), and genie::Intranuke::TransportHadrons().

RemoveIntermediateParticles()

void GHepRecord::RemoveIntermediateParticles ( void )

virtual

Definition at line 617 of file GHepRecord.cxx.

{
  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();
}

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().

ResetRecord()

void GHepRecord::ResetRecord ( void )

virtual

Definition at line 874 of file GHepRecord.cxx.

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

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

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

SetDiffXSec()

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

inlinevirtual

Definition at line 133 of file GHepRecord.h.

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

References fDiffXSec, and fDiffXSecPhSp.

Referenced by genie::COHKinematicsGenerator::CalculateKin_AlvarezRuso(), genie::SKKinematicsGenerator::CalculateKin_AtharSingleKaon(), genie::COHKinematicsGenerator::CalculateKin_BergerSehgal(), genie::COHKinematicsGenerator::CalculateKin_BergerSehgalFM(), genie::COHKinematicsGenerator::CalculateKin_ReinSehgal(), genie::KineGeneratorWithCache::MaxXSec(), genie::DFRKinematicsGenerator::ProcessEventRecord(), genie::DISKinematicsGenerator::ProcessEventRecord(), genie::DMDISKinematicsGenerator::ProcessEventRecord(), genie::DMEKinematicsGenerator::ProcessEventRecord(), genie::DMELEventGenerator::ProcessEventRecord(), genie::DMELKinematicsGenerator::ProcessEventRecord(), genie::HEDISKinematicsGenerator::ProcessEventRecord(), genie::HELeptonKinematicsGenerator::ProcessEventRecord(), genie::IBDKinematicsGenerator::ProcessEventRecord(), genie::NormGenerator::ProcessEventRecord(), genie::NuEKinematicsGenerator::ProcessEventRecord(), genie::QELEventGenerator::ProcessEventRecord(), genie::QELEventGeneratorSM::ProcessEventRecord(), genie::QELKinematicsGenerator::ProcessEventRecord(), genie::RESKinematicsGenerator::ProcessEventRecord(), genie::SPPEventGenerator::ProcessEventRecord(), genie::DMELKinematicsGenerator::SpectralFuncExperimentalCode(), and genie::QELKinematicsGenerator::SpectralFuncExperimentalCode().

SetPrintLevel()

void GHepRecord::SetPrintLevel ( int print_level )

static

Definition at line 956 of file GHepRecord.cxx.

{
  fPrintLevel = print_level;
}

References fPrintLevel.

Referenced by Initialize(), main(), and TestDecay().

SetProbability()

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

inlinevirtual

Definition at line 131 of file GHepRecord.h.

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

References fProb.

Referenced by main().

SetUnphysEventMask()

void GHepRecord::SetUnphysEventMask

(

const TBits &

mask

)

Definition at line 938 of file GHepRecord.cxx.

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

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

SetVertex() [1/2]

void GHepRecord::SetVertex ( const TLorentzVector & vtx )

virtual

Definition at line 832 of file GHepRecord.cxx.

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

References fVtx.

SetVertex() [2/2]

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

virtual

Definition at line 827 of file GHepRecord.cxx.

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

References fVtx.

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

SetWeight()

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

inlinevirtual

Definition at line 130 of file GHepRecord.h.

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

References fWeight.

Referenced by genie::COHKinematicsGenerator::CalculateKin_AlvarezRuso(), genie::SKKinematicsGenerator::CalculateKin_AtharSingleKaon(), genie::COHKinematicsGenerator::CalculateKin_ReinSehgal(), main(), genie::CascadeReweight::ProcessEventRecord(), genie::DFRKinematicsGenerator::ProcessEventRecord(), genie::DISKinematicsGenerator::ProcessEventRecord(), genie::DMDISKinematicsGenerator::ProcessEventRecord(), genie::DMEKinematicsGenerator::ProcessEventRecord(), genie::DMELKinematicsGenerator::ProcessEventRecord(), genie::hnl::FluxCreator::ProcessEventRecord(), genie::hnl::VertexGenerator::ProcessEventRecord(), genie::IBDKinematicsGenerator::ProcessEventRecord(), genie::NuEKinematicsGenerator::ProcessEventRecord(), genie::QELEventGeneratorSM::ProcessEventRecord(), genie::QELKinematicsGenerator::ProcessEventRecord(), genie::RESKinematicsGenerator::ProcessEventRecord(), and genie::SPPEventGenerator::ProcessEventRecord().

SetXSec()

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

inlinevirtual

Definition at line 132 of file GHepRecord.h.

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

References fXSec.

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

Summary()

Interaction * GHepRecord::Summary ( void ) const

virtual

Definition at line 91 of file GHepRecord.cxx.

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

References fInteraction, LOG, and pWARN.

Referenced by genie::DMETargetRemnantGenerator::AddElectronNeutrino(), genie::NuETargetRemnantGenerator::AddElectronNeutrino(), genie::HadronicSystemGenerator::AddFinalHadronicSyst(), genie::COHDNuEventGenerator::AddFinalStateDarkNeutrino(), genie::InitialStateAppender::AddNeutrino(), genie::InitialStateAppender::AddNucleus(), genie::HEDISGenerator::AddPrimaryLepton(), genie::IBDHadronicSystemGenerator::AddRecoilBaryon(), genie::QELHadronicSystemGenerator::AddRecoilBaryon(), genie::RSPPResonanceSelector::AddResonance(), genie::RSPPHadronicSystemGenerator::AddResonanceDecayProducts(), genie::InitialStateAppender::AddStruckParticle(), genie::DMETargetRemnantGenerator::AddTargetNucleusRemnant(), genie::NuETargetRemnantGenerator::AddTargetNucleusRemnant(), genie::OutgoingDarkGenerator::AddToEventRecord(), genie::PrimaryLeptonGenerator::AddToEventRecord(), genie::COHHadronicSystemGenerator::CalculateHadronicSystem_AlvarezRuso(), genie::SKHadronicSystemGenerator::CalculateHadronicSystem_AtharSingleKaon(), genie::COHHadronicSystemGenerator::CalculateHadronicSystem_BergerSehgalFM(), genie::COHHadronicSystemGenerator::CalculateHadronicSystem_ReinSehgal(), genie::COHKinematicsGenerator::CalculateKin_AlvarezRuso(), genie::SKKinematicsGenerator::CalculateKin_AtharSingleKaon(), genie::COHKinematicsGenerator::CalculateKin_BergerSehgal(), genie::COHKinematicsGenerator::CalculateKin_BergerSehgalFM(), genie::COHKinematicsGenerator::CalculateKin_ReinSehgal(), genie::COHPrimaryLeptonGenerator::CalculatePrimaryLepton_AlvarezRuso(), genie::RESHadronicSystemGenerator::GetResonancePdgCode(), genie::HadronicSystemGenerator::HadronShowerCharge(), main(), genie::KineGeneratorWithCache::MaxXSec(), genie::OutgoingDarkGenerator::NucRestFrame2Lab(), genie::PrimaryLeptonGenerator::NucRestFrame2Lab(), genie::DFRHadronicSystemGenerator::ProcessEventRecord(), genie::DFRKinematicsGenerator::ProcessEventRecord(), genie::DISHadronicSystemGenerator::ProcessEventRecord(), genie::DISKinematicsGenerator::ProcessEventRecord(), genie::DMDISKinematicsGenerator::ProcessEventRecord(), genie::DMEKinematicsGenerator::ProcessEventRecord(), genie::DMELEventGenerator::ProcessEventRecord(), genie::DMELKinematicsGenerator::ProcessEventRecord(), genie::DMEOutgoingDarkGenerator::ProcessEventRecord(), genie::FermiMover::ProcessEventRecord(), genie::HEDISKinematicsGenerator::ProcessEventRecord(), genie::HELeptonKinematicsGenerator::ProcessEventRecord(), genie::IBDKinematicsGenerator::ProcessEventRecord(), genie::NormGenerator::ProcessEventRecord(), genie::NucDeExcitationSim::ProcessEventRecord(), genie::NuEKinematicsGenerator::ProcessEventRecord(), genie::NuEPrimaryLeptonGenerator::ProcessEventRecord(), genie::OutgoingDarkGenerator::ProcessEventRecord(), genie::PauliBlocker::ProcessEventRecord(), genie::PrimaryLeptonGenerator::ProcessEventRecord(), genie::QELEventGenerator::ProcessEventRecord(), genie::QELEventGeneratorSM::ProcessEventRecord(), genie::QELEventGeneratorSuSA::ProcessEventRecord(), genie::QELKinematicsGenerator::ProcessEventRecord(), genie::RESKinematicsGenerator::ProcessEventRecord(), genie::RSPPResonanceSelector::ProcessEventRecord(), genie::SPPEventGenerator::ProcessEventRecord(), genie::VertexGenerator::ProcessEventRecord(), genie::RSPPResonanceSelector::SelectResonance(), genie::DMELKinematicsGenerator::SpectralFuncExperimentalCode(), and genie::QELKinematicsGenerator::SpectralFuncExperimentalCode().

SwapParticles()

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

protectedvirtual

Definition at line 752 of file GHepRecord.cxx.

{
#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);
      }
    }
  }
}

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

Referenced by CompactifyDaughterLists().

TargetNucleus()

GHepParticle * GHepRecord::TargetNucleus ( void ) const

virtual

Definition at line 293 of file GHepRecord.cxx.

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

References Particle(), and TargetNucleusPosition().

Referenced by genie::HNIntranuke2018::AbsorbHN(), genie::HNIntranuke2025::AbsorbHN(), genie::AMNuGammaGenerator::AddRecoilNucleon(), genie::CEvNSEventGenerator::AddRecoilNucleus(), genie::COHDNuEventGenerator::AddRecoilNucleus(), genie::DMELEventGenerator::AddTargetNucleusRemnant(), genie::FermiMover::AddTargetNucleusRemnant(), genie::HadronicSystemGenerator::AddTargetNucleusRemnant(), genie::QELEventGenerator::AddTargetNucleusRemnant(), genie::QELEventGeneratorSM::AddTargetNucleusRemnant(), genie::AMNuGammaGenerator::AddTargetRemnant(), genie::OutgoingDarkGenerator::AddToEventRecord(), genie::PrimaryLeptonGenerator::AddToEventRecord(), genie::COHHadronicSystemGenerator::CalculateHadronicSystem_AlvarezRuso(), genie::COHHadronicSystemGenerator::CalculateHadronicSystem_BergerSehgalFM(), genie::COHHadronicSystemGenerator::CalculateHadronicSystem_ReinSehgal(), genie::HAIntranuke2018::ElasHA(), genie::HAIntranuke2025::ElasHA(), genie::HAIntranuke::ElasHA(), genie::HNIntranuke2018::ElasHN(), genie::HNIntranuke2025::ElasHN(), genie::utils::intranuke::Equilibrium(), genie::HNIntranuke2018::GammaInelasticHN(), genie::HNIntranuke2025::GammaInelasticHN(), genie::Intranuke2018::GenerateVertex(), genie::Intranuke2025::GenerateVertex(), genie::Intranuke::GenerateVertex(), HitElectronPosition(), HitNucleonPosition(), genie::HAIntranuke2018::Inelastic(), genie::HAIntranuke2025::Inelastic(), genie::HAIntranuke::Inelastic(), genie::HAIntranuke2018::InelasticHA(), genie::HAIntranuke2025::InelasticHA(), genie::HAIntranuke::InelasticHA(), genie::FermiMover::KickHitNucleon(), genie::utils::intranuke::PreEquilibrium(), genie::DMELEventGenerator::ProcessEventRecord(), genie::HadronTransporter::ProcessEventRecord(), genie::Intranuke2018::ProcessEventRecord(), genie::Intranuke2025::ProcessEventRecord(), genie::Intranuke::ProcessEventRecord(), genie::NormGenerator::ProcessEventRecord(), genie::NucBindEnergyAggregator::ProcessEventRecord(), genie::NucDeExcitationSim::ProcessEventRecord(), genie::PauliBlocker::ProcessEventRecord(), genie::QELEventGenerator::ProcessEventRecord(), genie::SRCNuclearRecoil::ProcessEventRecord(), genie::VertexGenerator::ProcessEventRecord(), RemnantNucleusPosition(), genie::DISHadronicSystemGenerator::SimulateFormationZone(), genie::utils::intranuke2018::ThreeBodyKinematics(), genie::utils::intranuke2025::ThreeBodyKinematics(), genie::utils::intranuke::ThreeBodyKinematics(), genie::utils::intranuke2018::TwoBodyCollision(), genie::utils::intranuke2025::TwoBodyCollision(), and genie::utils::intranuke::TwoBodyCollision().

TargetNucleusPosition()

int GHepRecord::TargetNucleusPosition ( void ) const

virtual

Definition at line 370 of file GHepRecord.cxx.

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

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::COHDNuEventGenerator::AddFinalStateDarkNeutrino(), genie::CEvNSEventGenerator::AddRecoilNucleus(), genie::COHDNuEventGenerator::AddRecoilNucleus(), genie::DMELEventGenerator::AddTargetNucleusRemnant(), genie::DMETargetRemnantGenerator::AddTargetNucleusRemnant(), genie::FermiMover::AddTargetNucleusRemnant(), genie::HadronicSystemGenerator::AddTargetNucleusRemnant(), genie::NuETargetRemnantGenerator::AddTargetNucleusRemnant(), genie::QELEventGenerator::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::Intranuke2018::TransportHadrons(), genie::Intranuke2025::TransportHadrons(), and genie::Intranuke::TransportHadrons().

UpdateDaughterLists()

void GHepRecord::UpdateDaughterLists ( void )

protectedvirtual

Definition at line 556 of file GHepRecord.cxx.

{
  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();
}

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(), AddParticle(), and AddParticle().

Vertex()

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

inlinevirtual

Definition at line 140 of file GHepRecord.h.

{ return fVtx; }

References fVtx.

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

Weight()

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

inlinevirtual

Definition at line 124 of file GHepRecord.h.

{ return fWeight;   }

References fWeight.

Referenced by genie::COHKinematicsGenerator::CalculateKin_AlvarezRuso(), genie::SKKinematicsGenerator::CalculateKin_AtharSingleKaon(), genie::COHKinematicsGenerator::CalculateKin_ReinSehgal(), genie::AGCharmPythiaBaseHadro2023::ProcessEventRecord(), genie::AGKYLowW2019::ProcessEventRecord(), genie::DFRKinematicsGenerator::ProcessEventRecord(), genie::DISKinematicsGenerator::ProcessEventRecord(), genie::DMDISKinematicsGenerator::ProcessEventRecord(), genie::DMEKinematicsGenerator::ProcessEventRecord(), genie::DMELKinematicsGenerator::ProcessEventRecord(), genie::hnl::FluxCreator::ProcessEventRecord(), genie::hnl::VertexGenerator::ProcessEventRecord(), genie::IBDKinematicsGenerator::ProcessEventRecord(), genie::NuEKinematicsGenerator::ProcessEventRecord(), genie::QELEventGeneratorSM::ProcessEventRecord(), genie::QELKinematicsGenerator::ProcessEventRecord(), genie::RESKinematicsGenerator::ProcessEventRecord(), and genie::SPPEventGenerator::ProcessEventRecord().

XSec()

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

inlinevirtual

Definition at line 126 of file GHepRecord.h.

{ return fXSec;     }

References fXSec.

Referenced by genie::COHKinematicsGenerator::CalculateKin_AlvarezRuso(), genie::COHKinematicsGenerator::CalculateKin_ReinSehgal(), genie::DFRKinematicsGenerator::ProcessEventRecord(), genie::DISKinematicsGenerator::ProcessEventRecord(), genie::DMDISKinematicsGenerator::ProcessEventRecord(), genie::DMEKinematicsGenerator::ProcessEventRecord(), genie::DMELKinematicsGenerator::ProcessEventRecord(), genie::hnl::FluxCreator::ProcessEventRecord(), genie::IBDKinematicsGenerator::ProcessEventRecord(), genie::NuEKinematicsGenerator::ProcessEventRecord(), genie::QELEventGeneratorSM::ProcessEventRecord(), genie::QELKinematicsGenerator::ProcessEventRecord(), genie::RESKinematicsGenerator::ProcessEventRecord(), and genie::SPPEventGenerator::ProcessEventRecord().

Friends And Related Symbol Documentation

operator<<

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

friend

Definition at line 46 of file GHepRecord.cxx.

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

References GHepRecord(), and Print().

Member Data Documentation

fDiffXSec

double genie::GHepRecord::fDiffXSec

protected

differential cross section for selected event kinematics

Definition at line 181 of file GHepRecord.h.

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

fDiffXSecPhSp

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().

fEventFlags

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(), GHepRecord(), InitRecord(), IsUnphysical(), and Print().

fEventMask

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(), GHepRecord(), InitRecord(), Print(), and SetUnphysEventMask().

fInteraction

Interaction* genie::GHepRecord::fInteraction

protected

attached summary information

Definition at line 168 of file GHepRecord.h.

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

fPrintLevel

int genie::GHepRecord::fPrintLevel

staticprotected

Definition at line 196 of file GHepRecord.h.

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

fProb

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(), GHepRecord(), InitRecord(), Probability(), and SetProbability().

fVtx

TLorentzVector* genie::GHepRecord::fVtx

protected

vertex in the detector coordinate system

Definition at line 171 of file GHepRecord.h.

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

fWeight

double genie::GHepRecord::fWeight

protected

event weight

Definition at line 178 of file GHepRecord.h.

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

fXSec

double genie::GHepRecord::fXSec

protected

cross section for selected event

Definition at line 180 of file GHepRecord.h.

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

---

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

• GHepRecord.h
• GHepRecord.cxx