Generates the 'final state' hadronic system in v SPP interactions. It adds the remnant nucleus (if any) and the baryon resonance decay products at the GHEP record. The resonance decay products are pre- determined since in this thread we generate exclusive SPP reactions. The module uses a simple phase space decay. Is a concrete implementation of the EventRecordVisitorI interface. More...

#include <RSPPHadronicSystemGenerator.h>

Inheritance diagram for genie::RSPPHadronicSystemGenerator:

[legend]

Collaboration diagram for genie::RSPPHadronicSystemGenerator:

[legend]

Public Member Functions
	RSPPHadronicSystemGenerator ()
	RSPPHadronicSystemGenerator (string config)
	~RSPPHadronicSystemGenerator ()
void	ProcessEventRecord (GHepRecord *event_rec) const
Public Member Functions inherited from genie::HadronicSystemGenerator
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 ()
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.

Private Member Functions
void	AddResonanceDecayProducts (GHepRecord *event_rec) const

Private Attributes
TGenPhaseSpace	fPhaseSpaceGenerator

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)
Protected Member Functions inherited from genie::HadronicSystemGenerator
	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 inherited from genie::HadronicSystemGenerator
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)

Detailed Description

Generates the 'final state' hadronic system in v SPP interactions. It adds the remnant nucleus (if any) and the baryon resonance decay products at the GHEP record. The resonance decay products are pre- determined since in this thread we generate exclusive SPP reactions. The module uses a simple phase space decay. Is a concrete implementation of the EventRecordVisitorI interface.

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

Created:\n November 23, 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 32 of file RSPPHadronicSystemGenerator.h.

Constructor & Destructor Documentation

◆ RSPPHadronicSystemGenerator() [1/2]

RSPPHadronicSystemGenerator::RSPPHadronicSystemGenerator ( )

Definition at line 30 of file RSPPHadronicSystemGenerator.cxx.

                                                         :
HadronicSystemGenerator("genie::RSPPHadronicSystemGenerator")
{
 
}

References genie::HadronicSystemGenerator::HadronicSystemGenerator().

◆ RSPPHadronicSystemGenerator() [2/2]

RSPPHadronicSystemGenerator::RSPPHadronicSystemGenerator ( string config )

Definition at line 36 of file RSPPHadronicSystemGenerator.cxx.

                                                                     :
HadronicSystemGenerator("genie::RSPPHadronicSystemGenerator", config)
{
 
}

References genie::HadronicSystemGenerator::HadronicSystemGenerator().

◆ ~RSPPHadronicSystemGenerator()

RSPPHadronicSystemGenerator::~RSPPHadronicSystemGenerator ( )

Definition at line 42 of file RSPPHadronicSystemGenerator.cxx.

{
 
}

Member Function Documentation

◆ AddResonanceDecayProducts()

void RSPPHadronicSystemGenerator::AddResonanceDecayProducts ( GHepRecord * event_rec ) const

private

Definition at line 55 of file RSPPHadronicSystemGenerator.cxx.

{
// generate momenta for the baryon resonance decay products and add them at
// the event record
 
  //-- find out which SPP channel we are generating
  Interaction * interaction = evrec->Summary();
  SppChannel_t spp_channel = SppChannel::FromInteraction(interaction);
 
  //-- get the final state nucleon and pion
  int nuc_pdgc = SppChannel::FinStateNucleon (spp_channel);
  int pi_pdgc  = SppChannel::FinStatePion    (spp_channel);
 
  //-- get the total 4-p for the two-hadron system (= parent resonance 4-p)
 
  const InitialState & init_state = interaction->InitState();
  bool is_nucleus = init_state.Tgt().IsNucleus();
 
  //-- access the resonance entry at the GHEP record
  int res_pos = 0;
  if(is_nucleus) res_pos = 4;
  else           res_pos = 3;
 
  GHepParticle * res  = evrec->Particle(res_pos);
  const TLorentzVector & x4 = *(res->X4());
 
  //-- mark the resonance as decayed
  res->SetStatus(kIStDecayedState);
 
  //-- generate 4-p for the two-hadron system
  double mnuc = PDGLibrary::Instance() -> Find(nuc_pdgc) -> Mass();
  double mpi  = PDGLibrary::Instance() -> Find(pi_pdgc)  -> Mass();
 
  double mass[2] = { mnuc, mpi };
 
  TLorentzVector * p4 = res->GetP4();
 
  LOG("RESHadronicVtx", pINFO)
                 << "\n RES 4-P = " << utils::print::P4AsString(p4);
 
  bool is_permitted = fPhaseSpaceGenerator.SetDecay(*p4, 2, mass);
  assert(is_permitted);
 
  fPhaseSpaceGenerator.Generate();
 
  //-- add the two hadrons at the event record
  TLorentzVector & p4_nuc = *fPhaseSpaceGenerator.GetDecay(0);
  TLorentzVector & p4_pi  = *fPhaseSpaceGenerator.GetDecay(1);
  TLorentzVector vdummy(0,0,0,0); // dummy 'vertex'
 
  // decide the particle status
  GHepStatus_t ist = (is_nucleus) ?
                              kIStHadronInTheNucleus : kIStStableFinalState;
  int mom = res_pos;
  evrec->AddParticle(nuc_pdgc, ist, mom,-1,-1,-1, p4_nuc, x4);
  evrec->AddParticle(pi_pdgc,  ist, mom,-1,-1,-1, p4_pi,  x4);
  delete p4;
}

References genie::GHepRecord::AddParticle(), genie::SppChannel::FinStateNucleon(), genie::SppChannel::FinStatePion(), fPhaseSpaceGenerator, genie::SppChannel::FromInteraction(), genie::GHepParticle::GetP4(), genie::Interaction::InitState(), genie::PDGLibrary::Instance(), genie::Target::IsNucleus(), genie::kIStDecayedState, genie::kIStHadronInTheNucleus, genie::kIStStableFinalState, LOG, genie::utils::print::P4AsString(), genie::GHepRecord::Particle(), pINFO, genie::GHepParticle::SetStatus(), genie::GHepRecord::Summary(), genie::InitialState::Tgt(), and genie::GHepParticle::X4().

Referenced by ProcessEventRecord().

◆ ProcessEventRecord()

void RSPPHadronicSystemGenerator::ProcessEventRecord ( GHepRecord * event_rec ) const

virtual

Implements genie::EventRecordVisitorI.

Definition at line 47 of file RSPPHadronicSystemGenerator.cxx.

{
// This method generates the final state hadronic system
 
  //-- Add the baryon resonance decay products at the event record
  this->AddResonanceDecayProducts(evrec);
}

References AddResonanceDecayProducts().

Member Data Documentation

◆ fPhaseSpaceGenerator

TGenPhaseSpace genie::RSPPHadronicSystemGenerator::fPhaseSpaceGenerator

mutableprivate

Definition at line 45 of file RSPPHadronicSystemGenerator.h.

Referenced by AddResonanceDecayProducts().

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

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members