genie::InitialStateAppender Class Reference

Appends the initial state information to the event record. Is a concerete implementation of the EventRecordVisitorI interface. More...

#include <InitialStateAppender.h>

Inheritance diagram for genie::InitialStateAppender:

Collaboration diagram for genie::InitialStateAppender:

Public Member Functions
	InitialStateAppender ()
	InitialStateAppender (string config)
	~InitialStateAppender ()
void	ProcessEventRecord (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	AddNeutrino (GHepRecord *event_rec) const
void	AddNucleus (GHepRecord *event_rec) const
void	AddStruckParticle (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::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::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

Appends the initial state information to the event record. Is a concerete implementation of the EventRecordVisitorI interface.

Author

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

Created:\n October 04, 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 26 of file InitialStateAppender.h.

Constructor & Destructor Documentation

InitialStateAppender() [1/2]

InitialStateAppender::InitialStateAppender

(

)

Definition at line 26 of file InitialStateAppender.cxx.

                                           :
EventRecordVisitorI("genie::InitialStateAppender")
{
 
}

References genie::EventRecordVisitorI::EventRecordVisitorI().

InitialStateAppender() [2/2]

InitialStateAppender::InitialStateAppender

(

string

config

)

Definition at line 32 of file InitialStateAppender.cxx.

                                                        :
EventRecordVisitorI("genie::InitialStateAppender", config)
{
 
}

References genie::EventRecordVisitorI::EventRecordVisitorI().

~InitialStateAppender()

InitialStateAppender::~InitialStateAppender

(

)

Definition at line 38 of file InitialStateAppender.cxx.

{
 
}

Member Function Documentation

AddNeutrino()

void InitialStateAppender::AddNeutrino ( GHepRecord * event_rec ) const

private

Definition at line 66 of file InitialStateAppender.cxx.

{
  Interaction * interaction = evrec->Summary();
  const InitialState & init_state = interaction->InitState();
 
  TLorentzVector * p4 = init_state.GetProbeP4(kRfLab);
  const TLorentzVector v4(0.,0.,0.,0.);
 
  int pdgc = init_state.ProbePdg();
 
  LOG("ISApp", pINFO) << "Adding neutrino [pdgc = " << pdgc << "]";
 
  evrec->AddParticle(pdgc,kIStInitialState, -1,-1,-1,-1, *p4, v4);
 
  delete p4;
}

References genie::GHepRecord::AddParticle(), genie::InitialState::GetProbeP4(), genie::Interaction::InitState(), genie::kIStInitialState, genie::kRfLab, LOG, pINFO, genie::InitialState::ProbePdg(), and genie::GHepRecord::Summary().

Referenced by ProcessEventRecord().

AddNucleus()

void InitialStateAppender::AddNucleus ( GHepRecord * event_rec ) const

private

Definition at line 83 of file InitialStateAppender.cxx.

{
  Interaction * interaction = evrec->Summary();
  const InitialState & init_state = interaction->InitState();
  const ProcessInfo & proc_info   = interaction->ProcInfo();
 
  bool is_nucleus = init_state.Tgt().IsNucleus() || proc_info.IsNorm();
  if(!is_nucleus && !proc_info.IsGlashowResonance() && !proc_info.IsPhotonCoherent()) {
    LOG("ISApp", pINFO)
         << "Not an interaction with a nuclear target - no nucleus to add";
    return;
  }
  int    A    = init_state.Tgt().A();
  int    Z    = init_state.Tgt().Z();
  int    pdgc = pdg::IonPdgCode(A, Z);
  double M    = PDGLibrary::Instance()->Find(pdgc)->Mass();
 
  LOG("ISApp", pINFO)
          << "Adding nucleus [A = " << A << ", Z = " << Z
                                              << ", pdg = " << pdgc << "]";
 
  evrec->AddParticle(pdgc,kIStInitialState,-1,-1,-1,-1, 0,0,0,M, 0,0,0,0);
}

References genie::Target::A(), genie::GHepRecord::AddParticle(), genie::PDGLibrary::Find(), genie::Interaction::InitState(), genie::PDGLibrary::Instance(), genie::pdg::IonPdgCode(), genie::ProcessInfo::IsGlashowResonance(), genie::ProcessInfo::IsNorm(), genie::Target::IsNucleus(), genie::ProcessInfo::IsPhotonCoherent(), genie::kIStInitialState, LOG, pINFO, genie::Interaction::ProcInfo(), genie::GHepRecord::Summary(), genie::InitialState::Tgt(), and genie::Target::Z().

Referenced by ProcessEventRecord().

AddStruckParticle()

void InitialStateAppender::AddStruckParticle ( GHepRecord * event_rec ) const

private

Definition at line 107 of file InitialStateAppender.cxx.

{
  Interaction * interaction = evrec->Summary();
  const InitialState & init_state = interaction->InitState();
  const ProcessInfo & proc_info   = interaction->ProcInfo();
 
  // EDIT: Add the dark matter scattering off electron here
  bool hit_e = proc_info.IsInverseMuDecay()    ||
               proc_info.IsIMDAnnihilation()   ||
               proc_info.IsNuElectronElastic() ||
               proc_info.IsDarkMatterElectronElastic() ||
               proc_info.IsGlashowResonance();
 
  if(hit_e) {
    int    pdgc = kPdgElectron;
    double mass = PDGLibrary::Instance()->Find(pdgc)->Mass();
    const TLorentzVector p4(0,0,0, mass);
    const TLorentzVector v4(0.,0.,0.,0.);
 
    LOG("ISApp", pINFO) << "Adding struck electron";
    evrec->AddParticle(pdgc, kIStInitialState, 1, -1, -1, -1, p4, v4);
    return;
  }
 
  int pdgc = init_state.Tgt().HitNucPdg();
 
  if(pdgc != 0) {
 
    bool is_nucleus = init_state.Tgt().IsNucleus();
 
    GHepStatus_t ist   = (is_nucleus) ? kIStNucleonTarget : kIStInitialState;
    int          imom1 = (is_nucleus) ? 1 : -1;
    int          imom2 = -1;
 
    const TLorentzVector p4(init_state.Tgt().HitNucP4());
    const TLorentzVector v4(0.,0.,0.,0.);
 
    LOG("ISApp", pINFO)<< "Adding struck nucleon [pdgc = " << pdgc << "]";
 
    evrec->AddParticle(pdgc, ist, imom1, imom2, -1, -1, p4, v4);
 
  }//if struck nucleon was set
}

References genie::GHepRecord::AddParticle(), genie::PDGLibrary::Find(), genie::Target::HitNucP4(), genie::Target::HitNucPdg(), genie::Interaction::InitState(), genie::PDGLibrary::Instance(), genie::ProcessInfo::IsDarkMatterElectronElastic(), genie::ProcessInfo::IsGlashowResonance(), genie::ProcessInfo::IsIMDAnnihilation(), genie::ProcessInfo::IsInverseMuDecay(), genie::Target::IsNucleus(), genie::ProcessInfo::IsNuElectronElastic(), genie::kIStInitialState, genie::kIStNucleonTarget, genie::kPdgElectron, LOG, pINFO, genie::Interaction::ProcInfo(), genie::GHepRecord::Summary(), and genie::InitialState::Tgt().

Referenced by ProcessEventRecord().

ProcessEventRecord()

void InitialStateAppender::ProcessEventRecord ( GHepRecord * event_rec ) const

virtual

Implements genie::EventRecordVisitorI.

Definition at line 43 of file InitialStateAppender.cxx.

{
// Adds the initial state particles at the event record (the order is
// significant)
 
  LOG("ISApp", pINFO) << "Adding the initial state to the event record";
 
  //-- add the incoming neutrino to the event record
  this->AddNeutrino(evrec);
 
  //-- add the nuclear target at the event record (if any)
  this->AddNucleus(evrec);
 
  //-- add the struck nucleon to the event record (if any)
  //   It is added with status-code = 0 (init state) if the target was a
  //   free nucleon, or with a status-code = 11 (nucleon target) if the
  //   target was a nucleus.
  //   If the interaction was ve- elastic, inverse muon decay or Glashow
  //   resonance then it will add the target e- instead.
  this->AddStruckParticle(evrec);
  LOG("ISApp", pINFO) << *evrec;
}

References AddNeutrino(), AddNucleus(), AddStruckParticle(), LOG, and pINFO.

---

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

• InitialStateAppender.h
• InitialStateAppender.cxx