#include <Intranuke2018.h>

Inheritance diagram for genie::Intranuke2018:

Collaboration diagram for genie::Intranuke2018:

Public Member Functions
	Intranuke2018 ()
	Intranuke2018 (string name)
	Intranuke2018 (string name, string config)
	~Intranuke2018 ()
virtual void	ProcessEventRecord (GHepRecord *event_rec) const
virtual void	Configure (const Registry &config)
virtual void	Configure (string param_set)
virtual string	GetINukeMode () const
virtual string	GetGenINukeMode () const
void	SetRemnA (int A)
void	SetRemnZ (int Z)
double	GetRemnA () const
double	GetRemnZ () const
double	GetR0 () const
double	GetNR () const
double	GetDelRPion () const
double	GetDelRNucleon () const
double	GetNucRmvE () const
double	GetHadStep () const
bool	GetUseOset () const
bool	GetAltOset () const
bool	GetXsecNNCorr () const
Public Member Functions inherited from genie::EventRecordVisitorI
virtual	~EventRecordVisitorI ()
Public Member Functions inherited from genie::Algorithm
virtual	~Algorithm ()
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
virtual void	LoadConfig (void)=0
void	TransportHadrons (GHepRecord *ev) const
void	GenerateVertex (GHepRecord *ev) const
bool	NeedsRescattering (const GHepParticle *p) const
bool	CanRescatter (const GHepParticle *p) const
bool	IsInNucleus (const GHepParticle *p) const
void	SetTrackingRadius (const GHepParticle *p) const
double	GenerateStep (GHepRecord ev, GHepParticle p) const
virtual void	SimulateHadronicFinalState (GHepRecord ev, GHepParticle p) const =0
virtual int	HandleCompoundNucleus (GHepRecord ev, GHepParticle p, int mom) const =0
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
double	fTrackingRadius
	tracking radius for the nucleus in the current event
TGenPhaseSpace	fGenPhaseSpace
	a phase space generator
INukeHadroData2018 *	fHadroData2018
	a collection of h+N,h+A data & calculations
AlgFactory *	fAlgf
	algorithm factory instance
const NuclearModelI *	fNuclmodel
	nuclear model used to generate fermi momentum
int	fRemnA
	remnant nucleus A
int	fRemnZ
	remnant nucleus Z
TLorentzVector	fRemnP4
	P4 of remnant system.
GEvGenMode_t	fGMode
	event generation mode (lepton+A, hadron+A, ...)
double	fR0
	effective nuclear size param
double	fNR
	param multiplying the nuclear radius, determining how far to track hadrons beyond the "nuclear boundary"
double	fNucRmvE
	binding energy to subtract from cascade nucleons
double	fDelRPion
	factor by which Pion Compton wavelength gets multiplied to become nuclear size enhancement
double	fDelRNucleon
	factor by which Nucleon Compton wavelength gets multiplied to become nuclear size enhancement
double	fHadStep
	step size for intranuclear hadron transport
double	fNucAbsFac
	absorption xsec correction factor (hN Mode)
double	fNucCEXFac
	charge exchange xsec correction factor (hN Mode)
double	fEPreEq
	threshold for pre-equilibrium reaction
double	fFermiFac
	testing parameter to modify fermi momentum
double	fFermiMomentum
	whether or not particle collision is pauli blocked
bool	fDoFermi
	whether or not to do fermi mom.
bool	fDoMassDiff
	whether or not to do mass diff. mode
bool	fDoCompoundNucleus
	whether or not to do compound nucleus considerations
bool	fUseOset
	Oset model for low energy pion in hN.
bool	fAltOset
	NuWro's table-based implementation (not recommended)
bool	fXsecNNCorr
	use nuclear medium correction for NN cross section
double	fChPionMFPScale
	tweaking factors for tuning
double	fNeutralPionMFPScale
double	fPionFracCExScale
double	fPionFracInelScale
double	fChPionFracAbsScale
double	fNeutralPionFracAbsScale
double	fPionFracPiProdScale
double	fNucleonMFPScale
double	fNucleonFracCExScale
double	fNucleonFracInelScale
double	fNucleonFracAbsScale
double	fNucleonFracPiProdScale
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)

Friends
class	IntranukeTester

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

Definition at line 54 of file Intranuke2018.h.

Constructor & Destructor Documentation

◆ Intranuke2018() [1/3]

Intranuke2018::Intranuke2018 ( )

Definition at line 96 of file Intranuke2018.cxx.

                             :
EventRecordVisitorI()
{
 
}

References genie::EventRecordVisitorI::EventRecordVisitorI().

Referenced by genie::HAIntranuke2018::HAIntranuke2018(), genie::HAIntranuke2018::HAIntranuke2018(), genie::HNIntranuke2018::HNIntranuke2018(), and genie::HNIntranuke2018::HNIntranuke2018().

◆ Intranuke2018() [2/3]

Intranuke2018::Intranuke2018 ( string name )

Definition at line 102 of file Intranuke2018.cxx.

                                        :
EventRecordVisitorI(name)
{
 
}

References genie::EventRecordVisitorI::EventRecordVisitorI().

◆ Intranuke2018() [3/3]

Intranuke2018::Intranuke2018	(	string	name,
		string	config )

Definition at line 108 of file Intranuke2018.cxx.

                                                       :
EventRecordVisitorI(name, config)
{
 
}

References genie::EventRecordVisitorI::EventRecordVisitorI().

◆ ~Intranuke2018()

Intranuke2018::~Intranuke2018 ( )

Definition at line 114 of file Intranuke2018.cxx.

{
 
}

Member Function Documentation

◆ CanRescatter()

bool Intranuke2018::CanRescatter ( const GHepParticle * p ) const

protected

Definition at line 236 of file Intranuke2018.cxx.

{
// checks whether a particle that needs to be rescattered, can in fact be
// rescattered by this cascade MC
 
  assert(p);
  return  ( p->Pdg() == kPdgPiP     ||
            p->Pdg() == kPdgPiM     ||
            p->Pdg() == kPdgPi0     ||
            p->Pdg() == kPdgProton  ||
            p->Pdg() == kPdgNeutron ||
            //      p->Pdg() == kPdgGamma   ||
            p->Pdg() == kPdgKP      //||
            //      p->Pdg() == kPdgKM
          );
}

References genie::kPdgKP, genie::kPdgNeutron, genie::kPdgPi0, genie::kPdgPiM, genie::kPdgPiP, genie::kPdgProton, and genie::GHepParticle::Pdg().

Referenced by TransportHadrons().

◆ Configure() [1/2]

void Intranuke2018::Configure ( const Registry & config )

virtual

Configure the algorithm with an external registry The registry is merged with the top level registry if it is owned, Otherwise a copy of it is added with the highest priority

Reimplemented from genie::Algorithm.

Definition at line 432 of file Intranuke2018.cxx.

{
  Algorithm::Configure(config);
  this->LoadConfig();
}

References genie::Algorithm::Configure(), and LoadConfig().

◆ Configure() [2/2]

void Intranuke2018::Configure ( string config )

virtual

Configure the algorithm from the AlgoConfigPool based on param_set string given in input An algorithm contains a vector of registries coming from different xml configuration files, which are loaded according a very precise prioriy This methods will load a number registries in order of priority: 1) "Tunable" parameter set from CommonParametes. This is loaded with the highest prioriry and it is designed to be used for tuning procedure Usage not expected from the user. 2) For every string defined in "CommonParame" the corresponding parameter set will be loaded from CommonParameter.xml 3) parameter set specified by the config string and defined in the xml file of the algorithm 4) if config is not "Default" also the Default parameter set from the same xml file will be loaded Effectively this avoids the repetion of a parameter when it is not changed in the requested configuration

Reimplemented from genie::Algorithm.

Definition at line 438 of file Intranuke2018.cxx.

{
  Algorithm::Configure(param_set);
  this->LoadConfig();
}

References genie::Algorithm::Configure(), and LoadConfig().

◆ GenerateStep()

double Intranuke2018::GenerateStep	(	GHepRecord *	ev,
		GHepParticle *	p ) const

protected

Definition at line 393 of file Intranuke2018.cxx.

{
// Generate a step (in fermis) for particle p in the input event.
// Computes the mean free path L and generate an 'interaction' distance d
// from an exp(-d/L) distribution
 
  int pdgc = p->Pdg();
 
  double scale = 1.;
  if (pdgc==kPdgPiP || pdgc==kPdgPiM) {
    scale = fChPionMFPScale;
  }
  if (pdgc==kPdgPi0) {
    scale = fNeutralPionMFPScale;
  }
  else if (pdgc==kPdgProton || pdgc==kPdgNeutron) {
    scale = fNucleonMFPScale;
  }
 
  RandomGen * rnd = RandomGen::Instance();
 
  string fINukeMode = this->GetINukeMode();
  string fINukeModeGen = this->GetGenINukeMode();
 
  double L = utils::intranuke2018::MeanFreePath(p->Pdg(), *p->X4(), *p->P4(), fRemnA,
                                                fRemnZ, fDelRPion, fDelRNucleon, fUseOset, fAltOset, fXsecNNCorr, fINukeMode);
 
  LOG("Intranuke2018", pDEBUG)    << "mode= " << fINukeModeGen;
  L *= scale;
 
  double d = -1.*L * TMath::Log(rnd->RndFsi().Rndm());
 
  /*    LOG("Intranuke2018", pDEBUG)
    << "mode= " << fINukeMode << "; Mean free path = " << L << " fm / "
                              << "Generated path length = " << d << " fm";
  */
  return d;
}

References fAltOset, fChPionMFPScale, fDelRNucleon, fDelRPion, fNeutralPionMFPScale, fNucleonMFPScale, fRemnA, fRemnZ, fUseOset, fXsecNNCorr, GetGenINukeMode(), GetINukeMode(), genie::RandomGen::Instance(), genie::kPdgNeutron, genie::kPdgPi0, genie::kPdgPiM, genie::kPdgPiP, genie::kPdgProton, LOG, genie::utils::intranuke2018::MeanFreePath(), pDEBUG, and genie::RandomGen::RndFsi().

Referenced by TransportHadrons().

◆ GenerateVertex()

void Intranuke2018::GenerateVertex ( GHepRecord * ev ) const

protected

Definition at line 156 of file Intranuke2018.cxx.

{
// Sets a vertex in the nucleus periphery
// Called onlt in hadron/photon-nucleus interactions.
 
  GHepParticle * nucltgt = evrec->TargetNucleus();
  assert(nucltgt);
 
  RandomGen * rnd = RandomGen::Instance();
  TVector3 vtx(999999.,999999.,999999.);
 
  // *** For h+A events (test mode):
  // Assume a hadron beam with uniform intensity across an area,
  // so we need to choose events uniformly within that area.
  double x=999999., y=999999., epsilon = 0.001;
  double R2  = TMath::Power(fTrackingRadius,2.);
  double rp2 = TMath::Power(x,2.) + TMath::Power(y,2.);
  while(rp2 > R2-epsilon) {
      x = (fTrackingRadius-epsilon) * rnd->RndFsi().Rndm();
      y = -fTrackingRadius + 2*fTrackingRadius * rnd->RndFsi().Rndm();
      y -= ((y>0) ? epsilon : -epsilon);
      rp2 = TMath::Power(x,2.) + TMath::Power(y,2.);
  }
  vtx.SetXYZ(x,y, -1.*TMath::Sqrt(TMath::Max(0.,R2-rp2)) + epsilon);
 
  // get the actual unit vector along the incoming hadron direction
  TVector3 direction = evrec->Probe()->P4()->Vect().Unit();
 
  // rotate the vtx position
  vtx.RotateUz(direction);
 
  LOG("Intranuke2018", pNOTICE)
     << "Generated vtx @ R = " << vtx.Mag() << " fm / "
     << print::Vec3AsString(&vtx);
 
  TObjArrayIter piter(evrec);
  GHepParticle * p = 0;
  while( (p = (GHepParticle *) piter.Next()) )
  {
    if(pdg::IsPseudoParticle(p->Pdg())) continue;
    if(pdg::IsIon           (p->Pdg())) continue;
 
    p->SetPosition(vtx.x(), vtx.y(), vtx.z(), 0.);
  }
}

References epsilon, fTrackingRadius, genie::RandomGen::Instance(), genie::pdg::IsIon(), genie::pdg::IsPseudoParticle(), LOG, genie::GHepParticle::P4(), genie::GHepParticle::Pdg(), pNOTICE, genie::GHepRecord::Probe(), genie::RandomGen::RndFsi(), genie::GHepParticle::SetPosition(), genie::GHepRecord::TargetNucleus(), and genie::utils::print::Vec3AsString().

Referenced by ProcessEventRecord().

◆ GetAltOset()

bool genie::Intranuke2018::GetAltOset ( ) const

inline

Definition at line 92 of file Intranuke2018.h.

92{ return fAltOset; }

References fAltOset.

Referenced by genie::utils::intranuke2018::ProbSurvival().

◆ GetDelRNucleon()

double genie::Intranuke2018::GetDelRNucleon ( ) const

inline

Definition at line 86 of file Intranuke2018.h.

86{ return fDelRNucleon; }

References fDelRNucleon.

Referenced by genie::utils::intranuke2018::ProbSurvival().

◆ GetDelRPion()

double genie::Intranuke2018::GetDelRPion ( ) const

inline

Definition at line 85 of file Intranuke2018.h.

85{ return fDelRPion; }

References fDelRPion.

Referenced by genie::utils::intranuke2018::ProbSurvival().

◆ GetGenINukeMode()

virtual string genie::Intranuke2018::GetGenINukeMode ( ) const

inlinevirtual

Reimplemented in genie::HAIntranuke2018, and genie::HNIntranuke2018.

Definition at line 73 of file Intranuke2018.h.

73{return "XX";};

Referenced by GenerateStep().

◆ GetHadStep()

double genie::Intranuke2018::GetHadStep ( ) const

inline

Definition at line 89 of file Intranuke2018.h.

89{ return fHadStep; }

genie::Intranuke2018::fHadStep

double fHadStep

step size for intranuclear hadron transport

Definition Intranuke2018.h:130

References fHadStep.

Referenced by genie::utils::intranuke2018::ProbSurvival().

◆ GetINukeMode()

virtual string genie::Intranuke2018::GetINukeMode ( ) const

inlinevirtual

Reimplemented in genie::HAIntranuke2018, and genie::HNIntranuke2018.

Definition at line 72 of file Intranuke2018.h.

72{return "XX2018";};

Referenced by GenerateStep(), and genie::utils::intranuke2018::ProbSurvival().

◆ GetNR()

double genie::Intranuke2018::GetNR ( ) const

inline

Definition at line 83 of file Intranuke2018.h.

83{ return fNR; }

genie::Intranuke2018::fNR

double fNR

param multiplying the nuclear radius, determining how far to track hadrons beyond the "nuclear bounda...

Definition Intranuke2018.h:126

References fNR.

Referenced by genie::utils::intranuke2018::ProbSurvival().

◆ GetNucRmvE()

double genie::Intranuke2018::GetNucRmvE ( ) const

inline

Definition at line 88 of file Intranuke2018.h.

88{ return fNucRmvE; }

genie::Intranuke2018::fNucRmvE

double fNucRmvE

binding energy to subtract from cascade nucleons

Definition Intranuke2018.h:127

References fNucRmvE.

◆ GetR0()

double genie::Intranuke2018::GetR0 ( ) const

inline

Definition at line 82 of file Intranuke2018.h.

82{ return fR0; }

genie::Intranuke2018::fR0

double fR0

effective nuclear size param

Definition Intranuke2018.h:125

References fR0.

Referenced by genie::utils::intranuke2018::ProbSurvival().

◆ GetRemnA()

double genie::Intranuke2018::GetRemnA ( ) const

inline

Definition at line 79 of file Intranuke2018.h.

79{ return fRemnA; }

References fRemnA.

Referenced by genie::utils::intranuke2018::ProbSurvival().

◆ GetRemnZ()

double genie::Intranuke2018::GetRemnZ ( ) const

inline

Definition at line 80 of file Intranuke2018.h.

80{ return fRemnZ; }

References fRemnZ.

Referenced by genie::utils::intranuke2018::ProbSurvival().

◆ GetUseOset()

bool genie::Intranuke2018::GetUseOset ( ) const

inline

Definition at line 91 of file Intranuke2018.h.

91{ return fUseOset; }

References fUseOset.

Referenced by genie::utils::intranuke2018::ProbSurvival().

◆ GetXsecNNCorr()

bool genie::Intranuke2018::GetXsecNNCorr ( ) const

inline

Definition at line 93 of file Intranuke2018.h.

93{ return fXsecNNCorr; }

References fXsecNNCorr.

Referenced by genie::utils::intranuke2018::ProbSurvival().

◆ HandleCompoundNucleus()

virtual int genie::Intranuke2018::HandleCompoundNucleus	(	GHepRecord *	ev,
		GHepParticle *	p,
		int	mom ) const

protectedpure virtual

Implemented in genie::HAIntranuke2018, and genie::HNIntranuke2018.

◆ IsInNucleus()

bool Intranuke2018::IsInNucleus ( const GHepParticle * p ) const

protected

Definition at line 253 of file Intranuke2018.cxx.

{
// check whether the input particle is still within the nucleus
//
  return (p->X4()->Vect().Mag() < fTrackingRadius + fHadStep);
}

References fHadStep, fTrackingRadius, and genie::GHepParticle::X4().

Referenced by genie::HNIntranuke2018::HandleCompoundNucleus(), and TransportHadrons().

◆ LoadConfig()

virtual void genie::Intranuke2018::LoadConfig ( void )

protectedpure virtual

Implemented in genie::HAIntranuke2018, and genie::HNIntranuke2018.

Referenced by Configure(), and Configure().

◆ NeedsRescattering()

bool Intranuke2018::NeedsRescattering ( const GHepParticle * p ) const

protected

Definition at line 217 of file Intranuke2018.cxx.

{
// checks whether the particle should be rescattered
 
  assert(p);
 
  if(fGMode == kGMdHadronNucleus ||
     fGMode == kGMdPhotonNucleus) {
    // hadron/photon-nucleus scattering propagate the incoming particle
    return (
      (p->Status() == kIStInitialState || p->Status() == kIStHadronInTheNucleus)
       && !pdg::IsIon(p->Pdg()));
  }
  else {
   // attempt to rescatter anything marked as 'hadron in the nucleus'
   return (p->Status() == kIStHadronInTheNucleus);
  }
}

References fGMode, genie::pdg::IsIon(), genie::kGMdHadronNucleus, genie::kGMdPhotonNucleus, genie::kIStHadronInTheNucleus, genie::kIStInitialState, genie::GHepParticle::Pdg(), and genie::GHepParticle::Status().

Referenced by TransportHadrons().

◆ ProcessEventRecord()

void Intranuke2018::ProcessEventRecord ( GHepRecord * event_rec ) const

virtual

Implements genie::EventRecordVisitorI.

Reimplemented in genie::HAIntranuke2018, and genie::HNIntranuke2018.

Definition at line 119 of file Intranuke2018.cxx.

{
  // Do not continue if there is no nuclear target
  GHepParticle * nucltgt = evrec->TargetNucleus();
  if (!nucltgt) {
    LOG("HNIntranuke2018", pINFO) << "No nuclear target found - INTRANUKE exits";
    return;
  }
 
  // Decide tracking radius for the current nucleus (few * R0 * A^1/3)
  this->SetTrackingRadius(nucltgt);
 
  // Understand what the event generation mode is (hadron/photon-nucleus,
  // lepton-nucleus, nucleon decay) from the input event.
  // The determined mode has an effect on INTRANUKE behaviour (how to lookup
  // the residual nucleus, whether to set an intranuclear vtx etc) but it
  // does not affect the INTRANUKE physics.
  fGMode = evrec->EventGenerationMode();
 
  // For lepton-nucleus scattering and for nucleon decay intranuclear vtx
  // position (in the target nucleus coord system) is set elsewhere.
  // This method only takes effect in hadron/photon-nucleus interactions.
  // In this special mode, an interaction vertex is set at the periphery
  // of the target nucleus.
  if(fGMode == kGMdHadronNucleus ||
     fGMode == kGMdPhotonNucleus)
  {
    this->GenerateVertex(evrec);
  }
 
  // Now transport all hadrons outside the tracking radius.
  // Stepping part is common for both HA and HN.
  // Once it has been estabished that an interaction takes place then
  // HA and HN specific code takes over in order to simulate the final state.
  this->TransportHadrons(evrec);
}

References genie::GHepRecord::EventGenerationMode(), fGMode, GenerateVertex(), genie::kGMdHadronNucleus, genie::kGMdPhotonNucleus, LOG, pINFO, SetTrackingRadius(), genie::GHepRecord::TargetNucleus(), and TransportHadrons().

Referenced by genie::HAIntranuke2018::ProcessEventRecord(), and genie::HNIntranuke2018::ProcessEventRecord().

◆ SetRemnA()

void genie::Intranuke2018::SetRemnA ( int A )

inline

Definition at line 76 of file Intranuke2018.h.

76{ fRemnA = A; }

References fRemnA.

◆ SetRemnZ()

void genie::Intranuke2018::SetRemnZ ( int Z )

inline

Definition at line 77 of file Intranuke2018.h.

77{ fRemnZ = Z; }

References fRemnZ.

◆ SetTrackingRadius()

void Intranuke2018::SetTrackingRadius ( const GHepParticle * p ) const

protected

Definition at line 202 of file Intranuke2018.cxx.

{
  assert(p && pdg::IsIon(p->Pdg()));
  double A = p->A();
  fTrackingRadius = fR0 * TMath::Power(A, 1./3.);
 
  // multiply that by some input factor so that hadrons are tracked
  // beyond the nuclear 'boundary' since the nuclear density distribution
  // is not zero there
  fTrackingRadius *= fNR;
 
  LOG("Intranuke2018", pNOTICE)
      << "Setting tracking radius to R = " << fTrackingRadius;
}

References genie::GHepParticle::A(), fNR, fR0, fTrackingRadius, genie::pdg::IsIon(), LOG, genie::GHepParticle::Pdg(), and pNOTICE.

Referenced by ProcessEventRecord().

◆ SimulateHadronicFinalState()

virtual void genie::Intranuke2018::SimulateHadronicFinalState	(	GHepRecord *	ev,
		GHepParticle *	p ) const

protectedpure virtual

Implemented in genie::HAIntranuke2018, and genie::HNIntranuke2018.

Referenced by TransportHadrons().

◆ TransportHadrons()

void Intranuke2018::TransportHadrons ( GHepRecord * ev ) const

protected

Definition at line 260 of file Intranuke2018.cxx.

{
// transport all hadrons outside the nucleus
 
  int inucl = -1;
  fRemnA = -1;
  fRemnZ = -1;
 
  //  Get 'nuclear environment' at the beginning of hadron transport
  //  and keep track of the remnant nucleus A,Z
 
  if(fGMode == kGMdHadronNucleus ||
     fGMode == kGMdPhotonNucleus)
  {
     inucl = evrec->TargetNucleusPosition();
  }
  else if(fGMode == kGMdLeptonNucleus ||
          fGMode == kGMdDarkMatterNucleus ||
          fGMode == kGMdNucleonDecay) {
    inucl = evrec->RemnantNucleusPosition();
  }
 
  LOG("Intranuke2018", pNOTICE)
     << "Propagating hadrons within nucleus found in position = " << inucl;
  GHepParticle * nucl = evrec->Particle(inucl);
  if(!nucl) {
    LOG("Intranuke2018", pERROR)
       << "No nucleus found in position = " << inucl;
    LOG("Intranuke2018", pERROR)
       << *evrec;
    return;
  }
 
  fRemnA = nucl->A();
  fRemnZ = nucl->Z();
 
  LOG("Intranuke2018", pNOTICE)
      << "Nucleus (A,Z) = (" << fRemnA << ", " << fRemnZ << ")";
 
  const TLorentzVector & p4nucl = *(nucl->P4());
  fRemnP4 = p4nucl;
 
  // Loop over GHEP and run intranuclear rescattering on handled particles
  TObjArrayIter piter(evrec);
  GHepParticle * p = 0;
  int icurr = -1;
 
  while( (p = (GHepParticle *) piter.Next()) )
  {
    icurr++;
 
    // Check whether the particle needs rescattering, otherwise skip it
    if( ! this->NeedsRescattering(p) ) continue;
 
    LOG("Intranuke2018", pNOTICE)
      << " >> Stepping a " << p->Name()
                        << " with kinetic E = " << p->KinE() << " GeV";
 
    // Rescatter a clone, not the original particle
    GHepParticle * sp = new GHepParticle(*p);
 
    // Set clone's mom to be the hadron that was cloned
    sp->SetFirstMother(icurr);
 
    // Check whether the particle can be rescattered
    if(!this->CanRescatter(sp)) {
 
       // if I can't rescatter it, I will just take it out of the nucleus
       LOG("Intranuke2018", pNOTICE)
              << "... Current version can't rescatter a " << sp->Name();
       sp->SetFirstMother(icurr);
       sp->SetStatus(kIStStableFinalState);
       evrec->AddParticle(*sp);
       delete sp;
       continue; // <-- skip to next GHEP entry
    }
 
    // Start stepping particle out of the nucleus
    bool has_interacted = false;
    while ( this-> IsInNucleus(sp) )
    {
      // advance the hadron by a step
      utils::intranuke2018::StepParticle(sp, fHadStep);
 
      // check whether it interacts
      double d = this->GenerateStep(evrec,sp);
      has_interacted = (d<fHadStep);
      if(has_interacted) break;
    }//stepping
 
    if(has_interacted && fRemnA>0)  {
        // the particle interacts - simulate the hadronic interaction
      LOG("Intranuke2018", pNOTICE)
          << "Particle has interacted at location:  "
          << sp->X4()->Vect().Mag() << " / nucl rad= " << fTrackingRadius;
        this->SimulateHadronicFinalState(evrec,sp);
    } else if(has_interacted && fRemnA<=0) {
        // nothing left to interact with!
      LOG("Intranuke2018", pNOTICE)
          << "*** Nothing left to interact with, escaping.";
        sp->SetStatus(kIStStableFinalState);
        evrec->AddParticle(*sp);
        evrec->Particle(sp->FirstMother())->SetRescatterCode(1);
    } else {
        // the exits the nucleus without interacting - Done with it!
        LOG("Intranuke2018", pNOTICE)
          << "*** Hadron escaped the nucleus! Done with it.";
        sp->SetStatus(kIStStableFinalState);
        evrec->AddParticle(*sp);
        evrec->Particle(sp->FirstMother())->SetRescatterCode(1);
    }
    delete sp;
 
    // Current snapshot
    //LOG("Intranuke2018", pINFO) << "Current event record snapshot: " << *evrec;
 
  }// GHEP entries
 
  // Add remnant nucleus - that 'hadronic blob' has all the remaining hadronic
  // 4p not  put explicitly into the simulated particles
  TLorentzVector v4(0.,0.,0.,0.);
  GHepParticle remnant_nucleus(
    kPdgHadronicBlob, kIStFinalStateNuclearRemnant, inucl,-1,-1,-1, fRemnP4, v4);
  evrec->AddParticle(remnant_nucleus);
  // Mark the initial remnant nucleus as an intermediate state
  // Don't do that in the hadron/photon-nucleus scatterig mode since the initial
  // remnant nucleus and the target nucleus coincide.
  if(fGMode != kGMdHadronNucleus &&
     fGMode != kGMdPhotonNucleus) {
     evrec->Particle(inucl)->SetStatus(kIStIntermediateState);
  }
}

Referenced by ProcessEventRecord().

◆ IntranukeTester

friend class IntranukeTester

friend

Definition at line 56 of file Intranuke2018.h.

References IntranukeTester.

Referenced by IntranukeTester.

Member Data Documentation

◆ fAlgf

AlgFactory* genie::Intranuke2018::fAlgf

protected

algorithm factory instance

Definition at line 117 of file Intranuke2018.h.

◆ fAltOset

bool genie::Intranuke2018::fAltOset

protected

NuWro's table-based implementation (not recommended)

Definition at line 140 of file Intranuke2018.h.

Referenced by GenerateStep(), GetAltOset(), genie::HAIntranuke2018::LoadConfig(), and genie::HNIntranuke2018::LoadConfig().

◆ fChPionFracAbsScale

double genie::Intranuke2018::fChPionFracAbsScale

protected

Definition at line 147 of file Intranuke2018.h.

Referenced by genie::HAIntranuke2018::HadronFateHA(), and genie::HAIntranuke2018::LoadConfig().

◆ fChPionMFPScale

double genie::Intranuke2018::fChPionMFPScale

protected

tweaking factors for tuning

Definition at line 143 of file Intranuke2018.h.

Referenced by GenerateStep(), genie::HAIntranuke2018::LoadConfig(), and genie::HNIntranuke2018::LoadConfig().

◆ fDelRNucleon

double genie::Intranuke2018::fDelRNucleon

protected

factor by which Nucleon Compton wavelength gets multiplied to become nuclear size enhancement

Definition at line 129 of file Intranuke2018.h.

Referenced by GenerateStep(), GetDelRNucleon(), genie::HAIntranuke2018::LoadConfig(), and genie::HNIntranuke2018::LoadConfig().

◆ fDelRPion

double genie::Intranuke2018::fDelRPion

protected

factor by which Pion Compton wavelength gets multiplied to become nuclear size enhancement

Definition at line 128 of file Intranuke2018.h.

Referenced by GenerateStep(), GetDelRPion(), genie::HAIntranuke2018::LoadConfig(), and genie::HNIntranuke2018::LoadConfig().

◆ fDoCompoundNucleus

bool genie::Intranuke2018::fDoCompoundNucleus

protected

whether or not to do compound nucleus considerations

Definition at line 138 of file Intranuke2018.h.

Referenced by genie::HNIntranuke2018::FateWeight(), genie::HNIntranuke2018::HandleCompoundNucleus(), genie::HAIntranuke2018::LoadConfig(), and genie::HNIntranuke2018::LoadConfig().

◆ fDoFermi

bool genie::Intranuke2018::fDoFermi

protected

whether or not to do fermi mom.

Definition at line 136 of file Intranuke2018.h.

Referenced by genie::HNIntranuke2018::AbsorbHN(), genie::HNIntranuke2018::ElasHN(), genie::HNIntranuke2018::GammaInelasticHN(), genie::HNIntranuke2018::HandleCompoundNucleus(), genie::HAIntranuke2018::Inelastic(), genie::HAIntranuke2018::InelasticHA(), genie::HNIntranuke2018::InelasticHN(), genie::HAIntranuke2018::LoadConfig(), genie::HNIntranuke2018::LoadConfig(), genie::HNIntranuke2018::SimulateHadronicFinalState(), and genie::HAIntranuke2018::SimulateHadronicFinalStateKinematics().

◆ fDoMassDiff

bool genie::Intranuke2018::fDoMassDiff

protected

whether or not to do mass diff. mode

Definition at line 137 of file Intranuke2018.h.

◆ fEPreEq

double genie::Intranuke2018::fEPreEq

protected

threshold for pre-equilibrium reaction

Definition at line 133 of file Intranuke2018.h.

Referenced by genie::HNIntranuke2018::HandleCompoundNucleus(), genie::HAIntranuke2018::LoadConfig(), and genie::HNIntranuke2018::LoadConfig().

◆ fFermiFac

double genie::Intranuke2018::fFermiFac

protected

testing parameter to modify fermi momentum

Definition at line 134 of file Intranuke2018.h.

Referenced by genie::HNIntranuke2018::AbsorbHN(), genie::HNIntranuke2018::ElasHN(), genie::HNIntranuke2018::GammaInelasticHN(), genie::HNIntranuke2018::HandleCompoundNucleus(), genie::HAIntranuke2018::Inelastic(), genie::HAIntranuke2018::InelasticHA(), genie::HNIntranuke2018::InelasticHN(), genie::HAIntranuke2018::LoadConfig(), genie::HNIntranuke2018::LoadConfig(), genie::HNIntranuke2018::SimulateHadronicFinalState(), and genie::HAIntranuke2018::SimulateHadronicFinalStateKinematics().

◆ fFermiMomentum

double genie::Intranuke2018::fFermiMomentum

protected

whether or not particle collision is pauli blocked

Definition at line 135 of file Intranuke2018.h.

Referenced by genie::HNIntranuke2018::AbsorbHN(), genie::HAIntranuke2018::Inelastic(), genie::HNIntranuke2018::InelasticHN(), genie::HAIntranuke2018::LoadConfig(), and genie::HNIntranuke2018::LoadConfig().

◆ fGenPhaseSpace

TGenPhaseSpace genie::Intranuke2018::fGenPhaseSpace

mutableprotected

a phase space generator

Definition at line 115 of file Intranuke2018.h.

◆ fGMode

GEvGenMode_t genie::Intranuke2018::fGMode

mutableprotected

event generation mode (lepton+A, hadron+A, ...)

Definition at line 122 of file Intranuke2018.h.

Referenced by NeedsRescattering(), ProcessEventRecord(), and TransportHadrons().

◆ fHadroData2018

INukeHadroData2018* genie::Intranuke2018::fHadroData2018

protected

a collection of h+N,h+A data & calculations

Definition at line 116 of file Intranuke2018.h.

Referenced by genie::HNIntranuke2018::AbsorbHN(), genie::HNIntranuke2018::ElasHN(), genie::HNIntranuke2018::GammaInelasticHN(), genie::HAIntranuke2018::HadronFateHA(), genie::HNIntranuke2018::HadronFateHN(), genie::HAIntranuke2018::Inelastic(), genie::HAIntranuke2018::InelasticHA(), genie::HAIntranuke2018::LoadConfig(), and genie::HNIntranuke2018::LoadConfig().

◆ fHadStep

double genie::Intranuke2018::fHadStep

protected

step size for intranuclear hadron transport

Definition at line 130 of file Intranuke2018.h.

Referenced by GetHadStep(), IsInNucleus(), genie::HAIntranuke2018::LoadConfig(), genie::HNIntranuke2018::LoadConfig(), and TransportHadrons().

◆ fNeutralPionFracAbsScale

double genie::Intranuke2018::fNeutralPionFracAbsScale

protected

Definition at line 148 of file Intranuke2018.h.

Referenced by genie::HAIntranuke2018::HadronFateHA(), and genie::HAIntranuke2018::LoadConfig().

◆ fNeutralPionMFPScale

double genie::Intranuke2018::fNeutralPionMFPScale

protected

Definition at line 144 of file Intranuke2018.h.

Referenced by GenerateStep(), genie::HAIntranuke2018::LoadConfig(), and genie::HNIntranuke2018::LoadConfig().

◆ fNR

double genie::Intranuke2018::fNR

protected

param multiplying the nuclear radius, determining how far to track hadrons beyond the "nuclear boundary"

Definition at line 126 of file Intranuke2018.h.

Referenced by GetNR(), genie::HAIntranuke2018::LoadConfig(), genie::HNIntranuke2018::LoadConfig(), and SetTrackingRadius().

◆ fNucAbsFac

double genie::Intranuke2018::fNucAbsFac

protected

absorption xsec correction factor (hN Mode)

Definition at line 131 of file Intranuke2018.h.

Referenced by genie::HNIntranuke2018::HadronFateHN(), genie::HNIntranuke2018::HadronFateOset(), genie::HAIntranuke2018::LoadConfig(), and genie::HNIntranuke2018::LoadConfig().

◆ fNucCEXFac

double genie::Intranuke2018::fNucCEXFac

protected

charge exchange xsec correction factor (hN Mode)

Definition at line 132 of file Intranuke2018.h.

Referenced by genie::HNIntranuke2018::HadronFateHN(), genie::HNIntranuke2018::HadronFateOset(), genie::HAIntranuke2018::LoadConfig(), and genie::HNIntranuke2018::LoadConfig().

◆ fNucleonFracAbsScale

double genie::Intranuke2018::fNucleonFracAbsScale

protected

Definition at line 153 of file Intranuke2018.h.

Referenced by genie::HAIntranuke2018::HadronFateHA(), and genie::HAIntranuke2018::LoadConfig().

◆ fNucleonFracCExScale

double genie::Intranuke2018::fNucleonFracCExScale

protected

Definition at line 151 of file Intranuke2018.h.

Referenced by genie::HAIntranuke2018::HadronFateHA(), and genie::HAIntranuke2018::LoadConfig().

◆ fNucleonFracInelScale

double genie::Intranuke2018::fNucleonFracInelScale

protected

Definition at line 152 of file Intranuke2018.h.

Referenced by genie::HAIntranuke2018::HadronFateHA(), and genie::HAIntranuke2018::LoadConfig().

◆ fNucleonFracPiProdScale

double genie::Intranuke2018::fNucleonFracPiProdScale

protected

Definition at line 154 of file Intranuke2018.h.

Referenced by genie::HAIntranuke2018::HadronFateHA(), and genie::HAIntranuke2018::LoadConfig().

◆ fNucleonMFPScale

double genie::Intranuke2018::fNucleonMFPScale

protected

Definition at line 150 of file Intranuke2018.h.

Referenced by GenerateStep(), genie::HAIntranuke2018::LoadConfig(), and genie::HNIntranuke2018::LoadConfig().

◆ fNuclmodel

const NuclearModelI* genie::Intranuke2018::fNuclmodel

protected

nuclear model used to generate fermi momentum

Definition at line 118 of file Intranuke2018.h.

Referenced by genie::HNIntranuke2018::AbsorbHN(), genie::HNIntranuke2018::ElasHN(), genie::HNIntranuke2018::GammaInelasticHN(), genie::HNIntranuke2018::HandleCompoundNucleus(), genie::HAIntranuke2018::Inelastic(), genie::HAIntranuke2018::InelasticHA(), genie::HNIntranuke2018::InelasticHN(), genie::HAIntranuke2018::LoadConfig(), genie::HNIntranuke2018::LoadConfig(), genie::HNIntranuke2018::SimulateHadronicFinalState(), and genie::HAIntranuke2018::SimulateHadronicFinalStateKinematics().

◆ fNucRmvE

double genie::Intranuke2018::fNucRmvE

protected

binding energy to subtract from cascade nucleons

Definition at line 127 of file Intranuke2018.h.

Referenced by GetNucRmvE(), genie::HNIntranuke2018::HandleCompoundNucleus(), genie::HAIntranuke2018::Inelastic(), genie::HAIntranuke2018::LoadConfig(), genie::HNIntranuke2018::LoadConfig(), genie::HNIntranuke2018::SimulateHadronicFinalState(), and genie::HAIntranuke2018::SimulateHadronicFinalStateKinematics().

◆ fPionFracCExScale

double genie::Intranuke2018::fPionFracCExScale

protected

Definition at line 145 of file Intranuke2018.h.

Referenced by genie::HAIntranuke2018::HadronFateHA(), and genie::HAIntranuke2018::LoadConfig().

◆ fPionFracInelScale

double genie::Intranuke2018::fPionFracInelScale

protected

Definition at line 146 of file Intranuke2018.h.

Referenced by genie::HAIntranuke2018::HadronFateHA(), and genie::HAIntranuke2018::LoadConfig().

◆ fPionFracPiProdScale

double genie::Intranuke2018::fPionFracPiProdScale

protected

Definition at line 149 of file Intranuke2018.h.

Referenced by genie::HAIntranuke2018::HadronFateHA(), and genie::HAIntranuke2018::LoadConfig().

◆ fR0

double genie::Intranuke2018::fR0

protected

effective nuclear size param

Definition at line 125 of file Intranuke2018.h.

Referenced by GetR0(), genie::HAIntranuke2018::LoadConfig(), genie::HNIntranuke2018::LoadConfig(), and SetTrackingRadius().

◆ fRemnA

int genie::Intranuke2018::fRemnA

mutableprotected

remnant nucleus A

Definition at line 119 of file Intranuke2018.h.

Referenced by genie::HNIntranuke2018::AbsorbHN(), genie::HAIntranuke2018::ElasHA(), genie::HNIntranuke2018::ElasHN(), genie::HNIntranuke2018::FateWeight(), genie::HNIntranuke2018::GammaInelasticHN(), GenerateStep(), GetRemnA(), genie::HNIntranuke2018::HadronFateHN(), genie::HNIntranuke2018::HandleCompoundNucleus(), genie::HAIntranuke2018::Inelastic(), genie::HAIntranuke2018::InelasticHA(), genie::HNIntranuke2018::InelasticHN(), SetRemnA(), genie::HNIntranuke2018::SimulateHadronicFinalState(), genie::HAIntranuke2018::SimulateHadronicFinalStateKinematics(), and TransportHadrons().

◆ fRemnP4

TLorentzVector genie::Intranuke2018::fRemnP4

mutableprotected

P4 of remnant system.

Definition at line 121 of file Intranuke2018.h.

Referenced by genie::HNIntranuke2018::AbsorbHN(), genie::HAIntranuke2018::ElasHA(), genie::HNIntranuke2018::ElasHN(), genie::HNIntranuke2018::GammaInelasticHN(), genie::HNIntranuke2018::HandleCompoundNucleus(), genie::HAIntranuke2018::Inelastic(), genie::HAIntranuke2018::InelasticHA(), genie::HNIntranuke2018::InelasticHN(), genie::HNIntranuke2018::SimulateHadronicFinalState(), genie::HAIntranuke2018::SimulateHadronicFinalStateKinematics(), and TransportHadrons().

◆ fRemnZ

int genie::Intranuke2018::fRemnZ

mutableprotected

remnant nucleus Z

Definition at line 120 of file Intranuke2018.h.

Referenced by genie::HNIntranuke2018::AbsorbHN(), genie::HAIntranuke2018::ElasHA(), genie::HNIntranuke2018::ElasHN(), genie::HNIntranuke2018::FateWeight(), genie::HNIntranuke2018::GammaInelasticHN(), GenerateStep(), GetRemnZ(), genie::HNIntranuke2018::HadronFateHN(), genie::HNIntranuke2018::HandleCompoundNucleus(), genie::HAIntranuke2018::Inelastic(), genie::HAIntranuke2018::InelasticHA(), genie::HNIntranuke2018::InelasticHN(), SetRemnZ(), genie::HNIntranuke2018::SimulateHadronicFinalState(), genie::HAIntranuke2018::SimulateHadronicFinalStateKinematics(), and TransportHadrons().

◆ fTrackingRadius

double genie::Intranuke2018::fTrackingRadius

mutableprotected

tracking radius for the nucleus in the current event

Definition at line 114 of file Intranuke2018.h.

Referenced by GenerateVertex(), IsInNucleus(), SetTrackingRadius(), and TransportHadrons().

◆ fUseOset

bool genie::Intranuke2018::fUseOset

protected

Oset model for low energy pion in hN.

Definition at line 139 of file Intranuke2018.h.

Referenced by genie::HNIntranuke2018::AbsorbHN(), GenerateStep(), GetUseOset(), genie::HNIntranuke2018::HadronFateHN(), genie::HAIntranuke2018::LoadConfig(), and genie::HNIntranuke2018::LoadConfig().

◆ fXsecNNCorr

bool genie::Intranuke2018::fXsecNNCorr

protected

use nuclear medium correction for NN cross section

Definition at line 141 of file Intranuke2018.h.

Referenced by GenerateStep(), GetXsecNNCorr(), genie::HAIntranuke2018::LoadConfig(), and genie::HNIntranuke2018::LoadConfig().

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

Public Member Functions

Protected Member Functions

Protected Attributes

Friends

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ Intranuke2018() [1/3]

◆ Intranuke2018() [2/3]

◆ Intranuke2018() [3/3]

◆ ~Intranuke2018()

Member Function Documentation

◆ CanRescatter()

◆ Configure() [1/2]

◆ Configure() [2/2]

◆ GenerateStep()

◆ GenerateVertex()

◆ GetAltOset()

◆ GetDelRNucleon()

◆ GetDelRPion()

◆ GetGenINukeMode()

◆ GetHadStep()

◆ GetINukeMode()

◆ GetNR()

◆ GetNucRmvE()

◆ GetR0()

◆ GetRemnA()

◆ GetRemnZ()

◆ GetUseOset()

◆ GetXsecNNCorr()

◆ HandleCompoundNucleus()

◆ IsInNucleus()

◆ LoadConfig()

◆ NeedsRescattering()

◆ ProcessEventRecord()

◆ SetRemnA()

◆ SetRemnZ()

◆ SetTrackingRadius()

◆ SimulateHadronicFinalState()

◆ TransportHadrons()

Friends And Related Symbol Documentation

◆ IntranukeTester

Member Data Documentation

◆ fAlgf

◆ fAltOset

◆ fChPionFracAbsScale

◆ fChPionMFPScale

◆ fDelRNucleon

◆ fDelRPion

◆ fDoCompoundNucleus

◆ fDoFermi

◆ fDoMassDiff

◆ fEPreEq

◆ fFermiFac

◆ fFermiMomentum

◆ fGenPhaseSpace

◆ fGMode

◆ fHadroData2018

◆ fHadStep

◆ fNeutralPionFracAbsScale

◆ fNeutralPionMFPScale

◆ fNR

◆ fNucAbsFac

◆ fNucCEXFac

◆ fNucleonFracAbsScale

◆ fNucleonFracCExScale

◆ fNucleonFracInelScale

◆ fNucleonFracPiProdScale

◆ fNucleonMFPScale

◆ fNuclmodel

◆ fNucRmvE

◆ fPionFracCExScale

◆ fPionFracInelScale

◆ fPionFracPiProdScale

◆ fR0

◆ fRemnA

◆ fRemnP4

◆ fRemnZ

◆ fTrackingRadius

◆ fUseOset