GENIEGenerator
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Public Member Functions | Private Member Functions | List of all members
genie::QELHadronicSystemGenerator Class Reference

Generates the final state hadronic system in v QEL interactions. Is a concrete implementation of the EventRecordVisitorI interface. More...

#include <QELHadronicSystemGenerator.h>

Inheritance diagram for genie::QELHadronicSystemGenerator:
Inheritance graph
[legend]
Collaboration diagram for genie::QELHadronicSystemGenerator:
Collaboration graph
[legend]

Public Member Functions

 QELHadronicSystemGenerator ()
 
 QELHadronicSystemGenerator (string config)
 
 ~QELHadronicSystemGenerator ()
 
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 RegistryGetConfig (void) const
 
RegistryGetOwnedConfig (void)
 
virtual const AlgIdId (void) const
 Get algorithm ID. More...
 
virtual AlgStatus_t GetStatus (void) const
 Get algorithm status. More...
 
virtual bool AllowReconfig (void) const
 
virtual AlgCmp_t Compare (const Algorithm *alg) const
 Compare with input algorithm. More...
 
virtual void SetId (const AlgId &id)
 Set algorithm ID. More...
 
virtual void SetId (string name, string config)
 
const AlgorithmSubAlg (const RgKey &registry_key) const
 
void AdoptConfig (void)
 
void AdoptSubstructure (void)
 
virtual void Print (ostream &stream) const
 Print algorithm info. More...
 

Private Member Functions

void AddRecoilBaryon (GHepRecord *event_rec) const
 

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)
 
RegistryExtractLocalConfig (const Registry &in) const
 
RegistryExtractLowerConfig (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. More...
 
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. More...
 
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. More...
 
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 More...
 
int AddLowRegistry (Registry *rp, bool owns=true)
 add registry with lowest priority, also update ownership More...
 
int MergeTopRegistry (const Registry &r)
 
int AddTopRegisties (const vector< Registry * > &rs, bool owns=false)
 Add registries with top priority, also udated Ownerships. More...
 
- Protected Attributes inherited from genie::HadronicSystemGenerator
const EventRecordVisitorIfPreINukeDecayer
 
- Protected Attributes inherited from genie::Algorithm
bool fAllowReconfig
 
bool fOwnsSubstruc
 true if it owns its substructure (sub-algs,...) More...
 
AlgId fID
 algorithm name and configuration set More...
 
vector< Registry * > fConfVect
 
vector< bool > fOwnerships
 ownership for every registry in fConfVect More...
 
AlgStatus_t fStatus
 algorithm execution status More...
 
AlgMapfOwnedSubAlgMp
 local pool for owned sub-algs (taken out of the factory pool) More...
 

Detailed Description

Generates the final state hadronic system in v QEL interactions. Is a concrete implementation of the EventRecordVisitorI interface.

Author
Costas Andreopoulos <c.andreopoulos cern.ch> University of Liverpool
Created:
October 03, 2004
License:
Copyright (c) 2003-2024, The GENIE Collaboration For the full text of the license visit http://copyright.genie-mc.org

Definition at line 26 of file QELHadronicSystemGenerator.h.

Constructor & Destructor Documentation

QELHadronicSystemGenerator::QELHadronicSystemGenerator ( )

Definition at line 23 of file QELHadronicSystemGenerator.cxx.

23  :
24 HadronicSystemGenerator("genie::QELHadronicSystemGenerator")
25 {
26 
27 }
QELHadronicSystemGenerator::QELHadronicSystemGenerator ( string  config)

Definition at line 29 of file QELHadronicSystemGenerator.cxx.

29  :
30 HadronicSystemGenerator("genie::QELHadronicSystemGenerator", config)
31 {
32 
33 }
QELHadronicSystemGenerator::~QELHadronicSystemGenerator ( )

Definition at line 35 of file QELHadronicSystemGenerator.cxx.

36 {
37 
38 }

Member Function Documentation

void QELHadronicSystemGenerator::AddRecoilBaryon ( GHepRecord event_rec) const
private

Definition at line 51 of file QELHadronicSystemGenerator.cxx.

References genie::GHepRecord::AddParticle(), genie::XclsTag::CharmHadronPdg(), genie::Interaction::ExclTag(), genie::HadronicSystemGenerator::Hadronic4pLAB(), genie::GHepRecord::HitNucleonPosition(), genie::Interaction::InitState(), genie::XclsTag::IsCharmEvent(), genie::Target::IsNucleus(), genie::XclsTag::IsStrangeEvent(), genie::kIStHadronInTheNucleus, genie::kIStStableFinalState, LOG, genie::GHepRecord::Particle(), pINFO, genie::GHepRecord::Probe(), genie::Interaction::RecoilNucleonPdg(), genie::GHepParticle::RemovalEnergy(), genie::GHepParticle::SetRemovalEnergy(), genie::XclsTag::StrangeHadronPdg(), genie::GHepRecord::Summary(), genie::InitialState::Tgt(), and genie::GHepParticle::X4().

Referenced by ProcessEventRecord().

52 {
53  //-- Determine the pdg & status code of the recoil baryon
54  Interaction * interaction = evrec->Summary();
55  const XclsTag & xcls = interaction->ExclTag();
56  int pdgc = 0;
57  if (xcls.IsCharmEvent()) { pdgc = xcls.CharmHadronPdg(); }
58  else if(xcls.IsStrangeEvent()) { pdgc = xcls.StrangeHadronPdg(); }
59  else { pdgc = interaction->RecoilNucleonPdg(); }
60  assert(pdgc!=0);
61 
62  //-- Determine the status code
63  const Target & tgt = interaction->InitState().Tgt();
64  GHepStatus_t ist = (tgt.IsNucleus()) ?
66 
67  //-- Get the vtx position
68  GHepParticle * neutrino = evrec->Probe();
69  const TLorentzVector & vtx = *(neutrino->X4());
70 
71  //-- Get nucleon 4-momentum (in the LAB frame) & position
72  TLorentzVector p4 = this->Hadronic4pLAB(evrec);
73 
74  //-- Get mother position
75  int mom = evrec->HitNucleonPosition();
76 
77  //-- Add the final state recoil baryon at the EventRecord
78  LOG("QELHadronicVtx", pINFO)
79  << "Adding recoil baryon [pdgc = " << pdgc << "]";
80 
81  GHepParticle p(pdgc, ist, mom,-1,-1,-1, p4, vtx);
82  double w = ( xcls.IsCharmEvent() || xcls.IsStrangeEvent()) ?
83  0. : evrec->Particle(mom)->RemovalEnergy();
84 
85  p.SetRemovalEnergy(w);
86  evrec->AddParticle(p);
87 }
TLorentzVector Hadronic4pLAB(GHepRecord *event_rec) const
int RecoilNucleonPdg(void) const
recoil nucleon pdg
int CharmHadronPdg(void) const
Definition: XclsTag.h:52
bool IsStrangeEvent(void) const
Definition: XclsTag.h:53
bool IsNucleus(void) const
Definition: Target.cxx:272
Contains minimal information for tagging exclusive processes.
Definition: XclsTag.h:39
bool IsCharmEvent(void) const
Definition: XclsTag.h:50
Summary information for an interaction.
Definition: Interaction.h:56
int StrangeHadronPdg(void) const
Definition: XclsTag.h:55
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
Definition: Messenger.h:96
A Neutrino Interaction Target. Is a transparent encapsulation of quite different physical systems suc...
Definition: Target.h:40
#define pINFO
Definition: Messenger.h:62
const TLorentzVector * X4(void) const
Definition: GHepParticle.h:79
const XclsTag & ExclTag(void) const
Definition: Interaction.h:72
const InitialState & InitState(void) const
Definition: Interaction.h:69
const Target & Tgt(void) const
Definition: InitialState.h:66
STDHEP-like event record entry that can fit a particle or a nucleus.
Definition: GHepParticle.h:39
enum genie::EGHepStatus GHepStatus_t
void QELHadronicSystemGenerator::ProcessEventRecord ( GHepRecord event_rec) const
virtual

Implements genie::EventRecordVisitorI.

Definition at line 40 of file QELHadronicSystemGenerator.cxx.

References AddRecoilBaryon().

41 {
42 // This method generates the final state hadronic system
43 
44  // Add the recoil baryon
45  // (p or n - Lambda_c+,Sigma_c+,Sigma_c++ in charm/QEL)
46  // Its 4-momentum is computed by requiring the energy + momentum to be
47  // conserved.
48  this->AddRecoilBaryon(evrec);
49 }
void AddRecoilBaryon(GHepRecord *event_rec) const

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