Abstract class. Is used to pass some commonly recurring methods to all concrete implementations of the EventRecordVisitorI interface generating the hadronic system for a specific processes (QEL,DIS, RES,...) More...

#include <HadronicSystemGenerator.h>

Inheritance diagram for genie::HadronicSystemGenerator:

[legend]

Collaboration diagram for genie::HadronicSystemGenerator:

[legend]

Public Member Functions
void	AddTargetNucleusRemnant (GHepRecord *event_rec) const
void	AddFinalHadronicSyst (GHepRecord *event_rec) const
void	PreHadronTransportDecays (GHepRecord *event_rec) const
TLorentzVector	Hadronic4pLAB (GHepRecord *event_rec) const
TLorentzVector	MomentumTransferLAB (GHepRecord *event_rec) const
TVector3	HCM2LAB (GHepRecord *event_rec) const
int	HadronShowerCharge (GHepRecord *event_rec) const
int	ResonanceCharge (GHepRecord *event_rec) const
Public Member Functions inherited from genie::EventRecordVisitorI
virtual	~EventRecordVisitorI ()
virtual void	ProcessEventRecord (GHepRecord *event_rec) const =0
Public Member Functions inherited from genie::Algorithm
virtual	~Algorithm ()
virtual void	Configure (const Registry &config)
virtual void	Configure (string config)
virtual void	FindConfig (void)
virtual const Registry &	GetConfig (void) const
Registry *	GetOwnedConfig (void)
virtual const AlgId &	Id (void) const
	Get algorithm ID.
virtual AlgStatus_t	GetStatus (void) const
	Get algorithm status.
virtual bool	AllowReconfig (void) const
virtual AlgCmp_t	Compare (const Algorithm *alg) const
	Compare with input algorithm.
virtual void	SetId (const AlgId &id)
	Set algorithm ID.
virtual void	SetId (string name, string config)
const Algorithm *	SubAlg (const RgKey &registry_key) const
void	AdoptConfig (void)
void	AdoptSubstructure (void)
virtual void	Print (ostream &stream) const
	Print algorithm info.

Protected Member Functions
	HadronicSystemGenerator ()
	HadronicSystemGenerator (string name)
	HadronicSystemGenerator (string name, string config)
	~HadronicSystemGenerator ()
Protected Member Functions inherited from genie::EventRecordVisitorI
	EventRecordVisitorI ()
	EventRecordVisitorI (string name)
	EventRecordVisitorI (string name, string config)
Protected Member Functions inherited from genie::Algorithm
	Algorithm ()
	Algorithm (string name)
	Algorithm (string name, string config)
void	Initialize (void)
void	DeleteConfig (void)
void	DeleteSubstructure (void)
Registry *	ExtractLocalConfig (const Registry &in) const
Registry *	ExtractLowerConfig (const Registry &in, const string &alg_key) const
	Split an incoming configuration Registry into a block valid for the sub-algo identified by alg_key.
template<class T>
bool	GetParam (const RgKey &name, T &p, bool is_top_call=true) const
template<class T>
bool	GetParamDef (const RgKey &name, T &p, const T &def) const
template<class T>
int	GetParamVect (const std::string &comm_name, std::vector< T > &v, bool is_top_call=true) const
	Handle to load vectors of parameters.
int	GetParamVectKeys (const std::string &comm_name, std::vector< RgKey > &k, bool is_top_call=true) const
template<class T>
int	GetParamMat (const std::string &comm_name, TMatrixT< T > &mat, bool is_top_call=true) const
	Handle to load matrix of parameters.
template<class T>
int	GetParamMatSym (const std::string &comm_name, TMatrixTSym< T > &mat, bool is_top_call=true) const
int	GetParamMatKeys (const std::string &comm_name, std::vector< RgKey > &k, bool is_top_call=true) const
int	AddTopRegistry (Registry *rp, bool owns=true)
	add registry with top priority, also update ownership
int	AddLowRegistry (Registry *rp, bool owns=true)
	add registry with lowest priority, also update ownership
int	MergeTopRegistry (const Registry &r)
int	AddTopRegisties (const vector< Registry * > &rs, bool owns=false)
	Add registries with top priority, also udated Ownerships.

Protected Attributes
const EventRecordVisitorI *	fPreINukeDecayer
Protected Attributes inherited from genie::Algorithm
bool	fAllowReconfig
bool	fOwnsSubstruc
	true if it owns its substructure (sub-algs,...)
AlgId	fID
	algorithm name and configuration set
vector< Registry * >	fConfVect
vector< bool >	fOwnerships
	ownership for every registry in fConfVect
AlgStatus_t	fStatus
	algorithm execution status
AlgMap *	fOwnedSubAlgMp
	local pool for owned sub-algs (taken out of the factory pool)

Additional Inherited Members
Static Public Member Functions inherited from genie::Algorithm
static string	BuildParamVectKey (const std::string &comm_name, unsigned int i)
static string	BuildParamVectSizeKey (const std::string &comm_name)
static string	BuildParamMatKey (const std::string &comm_name, unsigned int i, unsigned int j)
static string	BuildParamMatRowSizeKey (const std::string &comm_name)
static string	BuildParamMatColSizeKey (const std::string &comm_name)

Detailed Description

Abstract class. Is used to pass some commonly recurring methods to all concrete implementations of the EventRecordVisitorI interface generating the hadronic system for a specific processes (QEL,DIS, RES,...)

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

Created:\n July 16, 2005

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 31 of file HadronicSystemGenerator.h.

Constructor & Destructor Documentation

◆ HadronicSystemGenerator() [1/3]

HadronicSystemGenerator::HadronicSystemGenerator ( )

protected

Definition at line 26 of file HadronicSystemGenerator.cxx.

                                                 :
EventRecordVisitorI()
{
 
}

References genie::EventRecordVisitorI::EventRecordVisitorI().

◆ HadronicSystemGenerator() [2/3]

HadronicSystemGenerator::HadronicSystemGenerator ( string name )

protected

Definition at line 32 of file HadronicSystemGenerator.cxx.

                                                            :
EventRecordVisitorI(name)
{
 
}

References genie::EventRecordVisitorI::EventRecordVisitorI().

◆ HadronicSystemGenerator() [3/3]

HadronicSystemGenerator::HadronicSystemGenerator	(	string	name,
		string	config )

protected

Definition at line 38 of file HadronicSystemGenerator.cxx.

                                                                          :
EventRecordVisitorI(name, config)
{
 
}

References genie::EventRecordVisitorI::EventRecordVisitorI().

◆ ~HadronicSystemGenerator()

HadronicSystemGenerator::~HadronicSystemGenerator ( )

protected

Definition at line 44 of file HadronicSystemGenerator.cxx.

{
 
}

Member Function Documentation

◆ AddFinalHadronicSyst()

void HadronicSystemGenerator::AddFinalHadronicSyst ( GHepRecord * event_rec ) const

Definition at line 49 of file HadronicSystemGenerator.cxx.

{
// Adds a GHEP entry for the sum of the f/s hadronic system.
// Intended for DIS hadronic system generators.
 
  TLorentzVector p4 = this->Hadronic4pLAB(evrec);
  LOG("HadronicVtx", pNOTICE) << "\n HadrSyst [LAB]: " << P4AsString(&p4);
 
  TLorentzVector v4(0,0,0,0);
  int mom = evrec->HitNucleonPosition();
 
  evrec->AddParticle(
       kPdgHadronicSyst, kIStDISPreFragmHadronicState, mom,-1,-1,-1, p4, v4);
 
  // update the interaction summary
  evrec->Summary()->KinePtr()->SetHadSystP4(p4);
}

References genie::GHepRecord::AddParticle(), Hadronic4pLAB(), genie::GHepRecord::HitNucleonPosition(), genie::Interaction::KinePtr(), genie::kIStDISPreFragmHadronicState, genie::kPdgHadronicSyst, LOG, genie::utils::print::P4AsString(), pNOTICE, genie::Kinematics::SetHadSystP4(), and genie::GHepRecord::Summary().

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

◆ AddTargetNucleusRemnant()

void HadronicSystemGenerator::AddTargetNucleusRemnant ( GHepRecord * event_rec ) const

Definition at line 67 of file HadronicSystemGenerator.cxx.

{
// add the remnant nuclear target at the GHEP record
 
  LOG("HadronicVtx", pDEBUG) << "Adding final state nucleus";
 
  //-- skip for non nuclear targets
  GHepParticle * nucleus = evrec->TargetNucleus();
  if (!nucleus) {
    LOG("HadronicVtx", pDEBUG)
               << "Initial state not a nucleus - no remnant nucleus to add";
    return;
  }
 
  //-- compute A,Z for final state nucleus & get its PDG code and its mass
  GHepParticle * nucleon = evrec->HitNucleon();
  assert(nucleon);
  int  npdgc = nucleon->Pdg();
  bool is_p  = pdg::IsProton(npdgc);
  int A = nucleus->A();
  int Z = nucleus->Z();
  if (is_p) Z--;
  A--;
  TParticlePDG * particle = 0;
  int ipdgc = pdg::IonPdgCode(A, Z);
  particle = PDGLibrary::Instance()->Find(ipdgc);
  if(!particle) {
      LOG("HadronicVtx", pFATAL)
          << "No particle with [A = " << A << ", Z = " << Z
                            << ", pdgc = " << ipdgc << "] in PDGLibrary!";
      assert(particle);
  }
  double Mf  = particle->Mass();   // remnant nucleus rest mass
  double Mf2 = TMath::Power(Mf,2);
 
  //-- Has opposite momentum from the struck nucleon
  double px = -1.* nucleon->Px();
  double py = -1.* nucleon->Py();
  double pz = -1.* nucleon->Pz();
  double E  = TMath::Sqrt(Mf2 + nucleon->P4()->Vect().Mag2());
 
  //-- Add the nucleus to the event record
  LOG("HadronicVtx", pINFO)
       << "Adding nucleus [A = " << A << ", Z = " << Z
                                           << ", pdgc = " << ipdgc << "]";
 
  int mom = evrec->TargetNucleusPosition();
  evrec->AddParticle(
           ipdgc,kIStStableFinalState, mom,-1,-1,-1, px,py,pz,E, 0,0,0,0);
}

References genie::GHepParticle::A(), genie::GHepRecord::AddParticle(), genie::PDGLibrary::Find(), genie::GHepRecord::HitNucleon(), genie::PDGLibrary::Instance(), genie::pdg::IonPdgCode(), genie::pdg::IsProton(), genie::kIStStableFinalState, LOG, genie::GHepParticle::P4(), pDEBUG, genie::GHepParticle::Pdg(), pFATAL, pINFO, genie::GHepParticle::Px(), genie::GHepParticle::Py(), genie::GHepParticle::Pz(), genie::GHepRecord::TargetNucleus(), genie::GHepRecord::TargetNucleusPosition(), and genie::GHepParticle::Z().

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

◆ Hadronic4pLAB()

TLorentzVector HadronicSystemGenerator::Hadronic4pLAB ( GHepRecord * event_rec ) const

Definition at line 127 of file HadronicSystemGenerator.cxx.

{
// Returns the final state hadronic system 4-p in LAB
 
  GHepParticle * nu = evrec->Probe(); // incoming v
  GHepParticle * N = evrec->HitNucleon();  // struck nucleon
  GHepParticle * l = evrec->FinalStatePrimaryLepton();  // f/s primary lepton
 
  assert(nu);
  assert(N);
  assert(l);
 
  LOG("HadronicVtx", pINFO)
      << "\n v [LAB]: " << P4AsString( nu->P4() )
      << "\n N [LAB]: " << P4AsString( N->P4()  )
      << "\n l [LAB]: " << P4AsString( l->P4()  );
 
  //-- Compute the Final State Hadronic System 4p (PX = Pv + PN - Pl)
 
  const TLorentzVector & p4nu = *(nu->P4());
  const TLorentzVector & p4N  = *(N ->P4());
  const TLorentzVector & p4l  = *(l ->P4());
 
  TLorentzVector pX4 = p4nu + p4N - p4l;
 
  LOG("HadronicVtx", pINFO) << "\n HadrSyst [LAB]: " << P4AsString(&pX4);
 
  return pX4;
}

References genie::GHepRecord::FinalStatePrimaryLepton(), genie::GHepRecord::HitNucleon(), LOG, genie::GHepParticle::P4(), genie::utils::print::P4AsString(), pINFO, and genie::GHepRecord::Probe().

Referenced by AddFinalHadronicSyst(), genie::IBDHadronicSystemGenerator::AddRecoilBaryon(), genie::QELHadronicSystemGenerator::AddRecoilBaryon(), genie::RESHadronicSystemGenerator::AddResonance(), genie::RSPPResonanceSelector::AddResonance(), HCM2LAB(), and genie::DISHadronicSystemGenerator::ProcessEventRecord().

◆ HadronShowerCharge()

int HadronicSystemGenerator::HadronShowerCharge ( GHepRecord * event_rec ) const

Definition at line 193 of file HadronicSystemGenerator.cxx.

{
// Returns the hadron shower charge in units of +e
// eg in v n -> l- X the hadron shower charge is +1
 
  int hadronShowerCharge = 0;
 
  Interaction * interaction = evrec->Summary();
  const InitialState & init_state = interaction->InitState();
 
  int hit_nucleon = init_state.Tgt().HitNucPdg();
 
  assert( pdg::IsProton(hit_nucleon) || pdg::IsNeutron(hit_nucleon) );
 
  double qfsl  = interaction->FSPrimLepton()->Charge() / 3.;
  double qp    = interaction->InitState().Probe()->Charge() / 3.;
  double qnuc  = PDGLibrary::Instance()->Find(hit_nucleon)->Charge() / 3.;
 
  // probe + nucleon - primary final state lepton
  hadronShowerCharge = (int) (qp + qnuc - qfsl);
 
  return hadronShowerCharge;
}

References genie::PDGLibrary::Find(), genie::Interaction::FSPrimLepton(), genie::Target::HitNucPdg(), genie::Interaction::InitState(), genie::PDGLibrary::Instance(), genie::pdg::IsNeutron(), genie::pdg::IsProton(), genie::InitialState::Probe(), genie::GHepRecord::Summary(), and genie::InitialState::Tgt().

Referenced by ResonanceCharge().

◆ HCM2LAB()

TVector3 HadronicSystemGenerator::HCM2LAB ( GHepRecord * event_rec ) const

Definition at line 177 of file HadronicSystemGenerator.cxx.

{
// Velocity for the Hadronic CM -> LAB active Lorentz transform
 
  TLorentzVector pH = this->Hadronic4pLAB(evrec);
 
  //-- Compute the velocity of the LAB frame in the Final State Hadronic
  //   CM Frame (PxH/EH, PyH/EH, PzH/EH)
 
  TVector3 beta = pH.BoostVector();
 
  LOG("HadronicVtx", pINFO) << "beta (HCM->LAB): " << Vec3AsString(&beta);
 
  return beta;
}

References Hadronic4pLAB(), LOG, pINFO, and genie::utils::print::Vec3AsString().

◆ MomentumTransferLAB()

TLorentzVector HadronicSystemGenerator::MomentumTransferLAB ( GHepRecord * event_rec ) const

Definition at line 158 of file HadronicSystemGenerator.cxx.

{
  GHepParticle * nu = evrec->Probe();  // incoming v
  GHepParticle * l = evrec->FinalStatePrimaryLepton();  // f/s primary lepton
 
  assert(nu);
  assert(l);
 
  const TLorentzVector & p4nu = *(nu->P4());
  const TLorentzVector & p4l  = *(l ->P4());
 
  TLorentzVector pq4 = p4nu - p4l; // q
 
  LOG("HadronicVtx", pNOTICE)
                      << "\n Momentum Transfer [LAB]: " << P4AsString(&pq4);
  return pq4;
}

References genie::GHepRecord::FinalStatePrimaryLepton(), LOG, genie::GHepParticle::P4(), genie::utils::print::P4AsString(), pNOTICE, and genie::GHepRecord::Probe().

◆ PreHadronTransportDecays()

void HadronicSystemGenerator::PreHadronTransportDecays ( GHepRecord * event_rec ) const

Definition at line 119 of file HadronicSystemGenerator.cxx.

{
  if(fPreINukeDecayer) {
      fPreINukeDecayer->ProcessEventRecord(evrec);
  }
}

References fPreINukeDecayer.

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

◆ ResonanceCharge()

int HadronicSystemGenerator::ResonanceCharge ( GHepRecord * event_rec ) const

Definition at line 217 of file HadronicSystemGenerator.cxx.

{
  return this->HadronShowerCharge(evrec);
}

References HadronShowerCharge().

Referenced by genie::RSPPResonanceSelector::AddResonance(), genie::RESHadronicSystemGenerator::GetResonancePdgCode(), and genie::RSPPResonanceSelector::SelectResonance().

Member Data Documentation

◆ fPreINukeDecayer

const EventRecordVisitorI* genie::HadronicSystemGenerator::fPreINukeDecayer

protected

Definition at line 59 of file HadronicSystemGenerator.h.

Referenced by genie::DISHadronicSystemGenerator::LoadConfig(), and PreHadronTransportDecays().

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

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HadronicSystemGenerator() [1/3]

◆ HadronicSystemGenerator() [2/3]

◆ HadronicSystemGenerator() [3/3]

◆ ~HadronicSystemGenerator()

Member Function Documentation

◆ AddFinalHadronicSyst()

◆ AddTargetNucleusRemnant()

◆ Hadronic4pLAB()

◆ HadronShowerCharge()

◆ HCM2LAB()

◆ MomentumTransferLAB()

◆ PreHadronTransportDecays()

◆ ResonanceCharge()

Member Data Documentation

◆ fPreINukeDecayer