Generates all the non-primary lepton final state particles in neutrino-electron events. Is a concrete implementation of the EventRecordVisitorI interface. More...

#include <NuETargetRemnantGenerator.h>

Inheritance diagram for genie::NuETargetRemnantGenerator:

[legend]

Collaboration diagram for genie::NuETargetRemnantGenerator:

[legend]

Public Member Functions
	NuETargetRemnantGenerator ()
	NuETargetRemnantGenerator (string config)
	~NuETargetRemnantGenerator ()
void	ProcessEventRecord (GHepRecord *evrec) 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	AddElectronNeutrino (GHepRecord *evrec) const
void	AddTargetNucleusRemnant (GHepRecord *evrec) 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

Generates all the non-primary lepton final state particles in neutrino-electron events. Is a concrete implementation of the EventRecordVisitorI interface.

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

Created:\n July 17, 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 27 of file NuETargetRemnantGenerator.h.

Constructor & Destructor Documentation

◆ NuETargetRemnantGenerator() [1/2]

NuETargetRemnantGenerator::NuETargetRemnantGenerator ( )

Definition at line 23 of file NuETargetRemnantGenerator.cxx.

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

References genie::EventRecordVisitorI::EventRecordVisitorI().

◆ NuETargetRemnantGenerator() [2/2]

NuETargetRemnantGenerator::NuETargetRemnantGenerator ( string config )

Definition at line 29 of file NuETargetRemnantGenerator.cxx.

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

References genie::EventRecordVisitorI::EventRecordVisitorI().

◆ ~NuETargetRemnantGenerator()

NuETargetRemnantGenerator::~NuETargetRemnantGenerator ( )

Definition at line 35 of file NuETargetRemnantGenerator.cxx.

{
 
}

Member Function Documentation

◆ AddElectronNeutrino()

void NuETargetRemnantGenerator::AddElectronNeutrino ( GHepRecord * evrec ) const

private

Definition at line 46 of file NuETargetRemnantGenerator.cxx.

{
  //-- Get all initial & final state particles 4-momenta (in the LAB frame)
 
  //incoming v:
  GHepParticle * nu = evrec->Probe();
 
  //struck particle
  GHepParticle * el = evrec->HitElectron();
 
  //final state primary lepton:
  GHepParticle * l = evrec->FinalStatePrimaryLepton();
 
  assert(nu);
  assert(el);
  assert(l);
 
  //-- Force energy conservation
  //   Pv(Ev,pxv,pyv,pzv) + Pe(En,pxn,pyn,pzn) = Pfsl(El,pxl,pyl,pzl) + Px
 
  const TLorentzVector & p4v = *(nu->P4());
  const TLorentzVector & p4e = *(el->P4());
  const TLorentzVector & p4l = *(l->P4());
  const TLorentzVector & p4  = p4v + p4e - p4l;
 
  //-- Vtx position
  const TLorentzVector & vtx = *(nu->X4());
 
  LOG("NuETargetRemnant", pINFO) << "Adding final state lepton from e- vtx";
 
  const ProcessInfo & proc_info = evrec->Summary()->ProcInfo();
  int mom  = evrec->HitElectronPosition();
  int pdgc = 0;
  if      (proc_info.IsNuElectronElastic()) pdgc = evrec->Summary()->InitState().ProbePdg();
  else if (proc_info.IsInverseMuDecay()) pdgc = kPdgNuE;
  else if (proc_info.IsIMDAnnihilation()) pdgc = kPdgAntiNuMu;
  assert(pdgc!=0);
  evrec->AddParticle(
           pdgc,kIStStableFinalState, mom,-1,-1,-1, p4, vtx);
}

References genie::GHepRecord::AddParticle(), genie::GHepRecord::FinalStatePrimaryLepton(), genie::GHepRecord::HitElectron(), genie::GHepRecord::HitElectronPosition(), genie::Interaction::InitState(), genie::ProcessInfo::IsIMDAnnihilation(), genie::ProcessInfo::IsInverseMuDecay(), genie::ProcessInfo::IsNuElectronElastic(), genie::kIStStableFinalState, genie::kPdgAntiNuMu, genie::kPdgNuE, LOG, genie::GHepParticle::P4(), pINFO, genie::GHepRecord::Probe(), genie::InitialState::ProbePdg(), genie::Interaction::ProcInfo(), genie::GHepRecord::Summary(), and genie::GHepParticle::X4().

Referenced by ProcessEventRecord().

◆ AddTargetNucleusRemnant()

void NuETargetRemnantGenerator::AddTargetNucleusRemnant ( GHepRecord * evrec ) const

private

Definition at line 87 of file NuETargetRemnantGenerator.cxx.

{
// add the remnant nuclear target at the GHEP record
 
  LOG("NuETargetRemnant", pDEBUG) << "Adding final state nucleus";
 
  //-- get A,Z for initial state nucleus
  Interaction * interaction = evrec->Summary();
  const InitialState & init_state = interaction->InitState();
 
  bool is_nucleus = init_state.Tgt().IsNucleus();
  if (!is_nucleus) {
    LOG("NuETargetRemnant", pDEBUG)
               << "Initial state not a nucleus - no remnant nucleus to add";
    return;
  }
  int A = init_state.Tgt().A();
  int Z = init_state.Tgt().Z();
 
  int    ipdgc = pdg::IonPdgCode(A, Z);
  double mass  = PDGLibrary::Instance()->Find(ipdgc)->Mass();
 
  //-- Add the nucleus to the event record
  LOG("NuETargetRemnant", pINFO)
          << "Adding nucleus [A = " << A << ", Z = " << Z
                                            << ", pdgc = " << ipdgc << "]";
 
  int mom = evrec->TargetNucleusPosition();
  evrec->AddParticle(
             ipdgc,kIStStableFinalState, mom,-1,-1,-1, 0,0,0,mass, 0,0,0,0);
}

References genie::Target::A(), genie::GHepRecord::AddParticle(), genie::PDGLibrary::Find(), genie::Interaction::InitState(), genie::PDGLibrary::Instance(), genie::pdg::IonPdgCode(), genie::Target::IsNucleus(), genie::kIStStableFinalState, LOG, pDEBUG, pINFO, genie::GHepRecord::Summary(), genie::GHepRecord::TargetNucleusPosition(), genie::InitialState::Tgt(), and genie::Target::Z().

Referenced by ProcessEventRecord().

◆ ProcessEventRecord()

void NuETargetRemnantGenerator::ProcessEventRecord ( GHepRecord * evrec ) const

virtual

Implements genie::EventRecordVisitorI.

Definition at line 40 of file NuETargetRemnantGenerator.cxx.

{
  this -> AddElectronNeutrino     (evrec);
  this -> AddTargetNucleusRemnant (evrec);
}

References AddElectronNeutrino(), and AddTargetNucleusRemnant().

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

Public Member Functions

Private Member Functions

Additional Inherited Members