#include <G4INCLEventInfo.hh>

Public Member Functions
	EventInfo ()

void	reset ()
	Reset the EventInfo members. More...

void	remnantToParticle (const G4int remnantIndex)
	Move a remnant to the particle array. More...

Data Fields
Short_t	nParticles
	Number of particles in the final state. More...

Short_t	A [maxSizeParticles]
	Particle mass number. More...

Short_t	Z [maxSizeParticles]
	Particle charge number. More...

Float_t	EKin [maxSizeParticles]
	Particle kinetic energy [MeV]. More...

Float_t	px [maxSizeParticles]
	Particle momentum, x component [MeV/c]. More...

Float_t	py [maxSizeParticles]
	Particle momentum, y component [MeV/c]. More...

Float_t	pz [maxSizeParticles]
	Particle momentum, z component [MeV/c]. More...

Float_t	theta [maxSizeParticles]
	Particle momentum polar angle [radians]. More...

Float_t	phi [maxSizeParticles]
	Particle momentum azimuthal angle [radians]. More...

Short_t	origin [maxSizeParticles]
	Origin of the particle. More...

Float_t	emissionTime [maxSizeParticles]
	Emission time [fm/c]. More...

std::vector< std::string >	history
	History of the particle. More...

Int_t	nRemnants
	Number of remnants. More...

Short_t	ARem [maxSizeRemnants]
	Remnant mass number. More...

Short_t	ZRem [maxSizeRemnants]
	Remnant charge number. More...

Float_t	EStarRem [maxSizeRemnants]
	Remnant excitation energy [MeV]. More...

Float_t	JRem [maxSizeRemnants]
	Remnant spin [ $\hbar$ ]. More...

Float_t	EKinRem [maxSizeRemnants]
	Remnant kinetic energy [MeV]. More...

Float_t	pxRem [maxSizeRemnants]
	Remnant momentum, x component [MeV/c]. More...

Float_t	pyRem [maxSizeRemnants]
	Remnant momentum, y component [MeV/c]. More...

Float_t	pzRem [maxSizeRemnants]
	Remnant momentum, z component [MeV/c]. More...

Float_t	thetaRem [maxSizeRemnants]
	Remnant momentum polar angle [radians]. More...

Float_t	phiRem [maxSizeRemnants]
	Remnant momentum azimuthal angle [radians]. More...

Float_t	jxRem [maxSizeRemnants]
	Remnant angular momentum, x component [ $\hbar$ ]. More...

Float_t	jyRem [maxSizeRemnants]
	Remnant angular momentum, y component [ $\hbar$ ]. More...

Float_t	jzRem [maxSizeRemnants]
	Remnant angular momentum, z component [ $\hbar$ ]. More...

Int_t	projectileType
	Projectile particle type. More...

Short_t	At
	Mass number of the target nucleus. More...

Short_t	Zt
	Charge number of the target nucleus. More...

Short_t	Ap
	Mass number of the projectile nucleus. More...

Short_t	Zp
	Charge number of the projectile nucleus. More...

Float_t	Ep
	Projectile kinetic energy given as input. More...

Float_t	impactParameter
	Impact parameter [fm]. More...

Int_t	nCollisions
	Number of accepted two-body collisions. More...

Float_t	stoppingTime
	Cascade stopping time [fm/c]. More...

Float_t	EBalance
	Energy-conservation balance [MeV]. More...

Float_t	pLongBalance
	Longitudinal momentum-conservation balance [MeV/c]. More...

Float_t	pTransBalance
	Transverse momentum-conservation balance [MeV/c]. More...

Short_t	nCascadeParticles
	Number of cascade particles. More...

Bool_t	transparent
	True if the event is transparent. More...

Bool_t	forcedCompoundNucleus
	True if the event is a forced CN. More...

Bool_t	nucleonAbsorption
	True if the event is a nucleon absorption. More...

Bool_t	pionAbsorption
	True if the event is a pion absorption. More...

Int_t	nDecays
	Number of accepted Delta decays. More...

Int_t	nBlockedCollisions
	Number of two-body collisions blocked by Pauli or CDPP. More...

Int_t	nBlockedDecays
	Number of decays blocked by Pauli or CDPP. More...

Float_t	effectiveImpactParameter
	Effective (Coulomb-distorted) impact parameter [fm]. More...

Bool_t	deltasInside
	Event involved deltas in the nucleus at the end of the cascade. More...

Bool_t	forcedDeltasInside
	Event involved forced delta decays inside the nucleus. More...

Bool_t	forcedDeltasOutside
	Event involved forced delta decays outside the nucleus. More...

Bool_t	clusterDecay
	Event involved cluster decay. More...

Float_t	firstCollisionTime
	Time of the first collision [fm/c]. More...

Float_t	firstCollisionXSec
	Cross section of the first collision (mb) More...

Float_t	firstCollisionSpectatorPosition
	Position of the spectator on the first collision (fm) More...

Float_t	firstCollisionSpectatorMomentum
	Momentum of the spectator on the first collision (fm) More...

Bool_t	firstCollisionIsElastic
	True if the first collision was elastic. More...

Int_t	nReflectionAvatars
	Number of reflection avatars. More...

Int_t	nCollisionAvatars
	Number of collision avatars. More...

Int_t	nDecayAvatars
	Number of decay avatars. More...

Int_t	nUnmergedSpectators
	Number of dynamical spectators that were merged back into the projectile remnant. More...

Int_t	nEnergyViolationInteraction
	Number of attempted collisions/decays for which the energy-conservation algorithm failed to find a solution. More...

Static Public Attributes
static G4ThreadLocal Int_t	eventNumber = 0
	Number of the event. More...

static const Short_t	maxSizeRemnants = 10
	Maximum array size for remnants. More...

static const Short_t	maxSizeParticles = 1000
	Maximum array size for emitted particles. More...

Detailed Description

Definition at line 66 of file G4INCLEventInfo.hh.

Constructor & Destructor Documentation

G4INCL::EventInfo::EventInfo ( )

inline

Definition at line 67 of file G4INCLEventInfo.hh.

References A, ARem, EKin, EKinRem, emissionTime, EStarRem, JRem, jxRem, jyRem, jzRem, maxSizeParticles, maxSizeRemnants, origin, phi, phiRem, px, pxRem, py, pyRem, pz, pzRem, theta, thetaRem, Z, and ZRem.

                   :
         nParticles(0),
         nRemnants(0),
         projectileType(0),
         At(0),
         Zt(0),
         Ap(0),
         Zp(0),
         Ep((Float_t)0.0),
         impactParameter((Float_t)0.0),
         nCollisions(0),
         stoppingTime((Float_t)0.0),
         EBalance((Float_t)0.0),
         pLongBalance((Float_t)0.0),
         pTransBalance((Float_t)0.0),
         nCascadeParticles(0),
         transparent(false),
         forcedCompoundNucleus(false),
         nucleonAbsorption(false),
         pionAbsorption(false),
         nDecays(0),
         nBlockedCollisions(0),
         nBlockedDecays(0),
         effectiveImpactParameter((Float_t)0.0),
         deltasInside(false),
         forcedDeltasInside(false),
         forcedDeltasOutside(false),
         clusterDecay(false),
         firstCollisionTime((Float_t)0.0),
         firstCollisionXSec((Float_t)0.0),
         firstCollisionSpectatorPosition((Float_t)0.0),
         firstCollisionSpectatorMomentum((Float_t)0.0),
         firstCollisionIsElastic(false),
         nReflectionAvatars(0),
         nCollisionAvatars(0),
         nDecayAvatars(0),
         nUnmergedSpectators(0),
         nEnergyViolationInteraction(0)
 #ifdef INCL_INVERSE_KINEMATICS
 
 #endif
       {
         std::fill_n(A, maxSizeParticles, 0);
         std::fill_n(Z, maxSizeParticles, 0);
         std::fill_n(EKin, maxSizeParticles, (Float_t)0.0);
         std::fill_n(px, maxSizeParticles, (Float_t)0.0);
         std::fill_n(py, maxSizeParticles, (Float_t)0.0);
         std::fill_n(pz, maxSizeParticles, (Float_t)0.0);
         std::fill_n(theta, maxSizeParticles, (Float_t)0.0);
         std::fill_n(phi, maxSizeParticles, (Float_t)0.0);
         std::fill_n(origin, maxSizeParticles, 0);
         std::fill_n(emissionTime, maxSizeParticles, (Float_t)0.0);
         std::fill_n(ARem, maxSizeRemnants, 0);
         std::fill_n(ZRem, maxSizeRemnants, 0);
         std::fill_n(EStarRem, maxSizeRemnants, (Float_t)0.0);
         std::fill_n(JRem, maxSizeRemnants, (Float_t)0.0);
         std::fill_n(EKinRem, maxSizeRemnants, (Float_t)0.0);
         std::fill_n(pxRem, maxSizeRemnants, (Float_t)0.0);
         std::fill_n(pyRem, maxSizeRemnants, (Float_t)0.0);
         std::fill_n(pzRem, maxSizeRemnants, (Float_t)0.0);
         std::fill_n(thetaRem, maxSizeRemnants, (Float_t)0.0);
         std::fill_n(phiRem, maxSizeRemnants, (Float_t)0.0);
         std::fill_n(jxRem, maxSizeRemnants, (Float_t)0.0);
         std::fill_n(jyRem, maxSizeRemnants, (Float_t)0.0);
         std::fill_n(jzRem, maxSizeRemnants, (Float_t)0.0);
 #ifdef INCL_INVERSE_KINEMATICS
         std::fill_n(EKinPrime, maxSizeParticles, (Float_t)0.0);
         std::fill_n(pzPrime, maxSizeParticles, (Float_t)0.0);
         std::fill_n(thetaPrime, maxSizeParticles, (Float_t)0.0);
 #endif
       }

Member Function Documentation

void G4INCL::EventInfo::remnantToParticle ( const G4int remnantIndex )

Move a remnant to the particle array.

Definition at line 81 of file G4INCLEventInfo.cc.

References A, ARem, EKin, EKinRem, emissionTime, history, nParticles, origin, phi, G4INCL::Math::pi, px, pxRem, py, pyRem, pz, pzRem, stoppingTime, theta, Z, and ZRem.

                                                             {
     A[nParticles] = ARem[remnantIndex];
     Z[nParticles] = ZRem[remnantIndex];
     emissionTime[nParticles] = stoppingTime;
 
     px[nParticles] = pxRem[remnantIndex];
     py[nParticles] = pyRem[remnantIndex];
     pz[nParticles] = pzRem[remnantIndex];
 
     const G4double plab = std::sqrt(pxRem[remnantIndex]*pxRem[remnantIndex]
                                   +pyRem[remnantIndex]*pyRem[remnantIndex]
                                   +pzRem[remnantIndex]*pzRem[remnantIndex]);
     G4double pznorm = pzRem[remnantIndex]/plab;
     if(pznorm>1.)
       pznorm = 1.;
     else if(pznorm<-1.)
       pznorm = -1.;
     theta[nParticles] = 180.*std::acos(pznorm)/G4INCL::Math::pi;
     phi[nParticles] = 180.*std::atan2(pyRem[remnantIndex],pxRem[remnantIndex])/G4INCL::Math::pi;
 
     EKin[nParticles] = EKinRem[remnantIndex];
     origin[nParticles] = -1; // Origin: cascade
     history.push_back(""); // history
     nParticles++;
 // assert(history.size()==(unsigned int)nParticles);
   }

void G4INCL::EventInfo::reset ( )

inline

Reset the EventInfo members.

Definition at line 301 of file G4INCLEventInfo.hh.

References Ap, At, clusterDecay, deltasInside, EBalance, effectiveImpactParameter, Ep, firstCollisionIsElastic, firstCollisionSpectatorMomentum, firstCollisionSpectatorPosition, firstCollisionTime, firstCollisionXSec, forcedCompoundNucleus, forcedDeltasInside, forcedDeltasOutside, history, impactParameter, nBlockedCollisions, nBlockedDecays, nCascadeParticles, nCollisionAvatars, nCollisions, nDecayAvatars, nDecays, nEnergyViolationInteraction, nParticles, nReflectionAvatars, nRemnants, nucleonAbsorption, nUnmergedSpectators, pionAbsorption, pLongBalance, projectileType, pTransBalance, stoppingTime, transparent, Zp, and Zt.

                    {
         nParticles = 0;
         history.clear();
         nRemnants = 0;
         projectileType = 0;
         At = 0;
         Zt = 0;
         Ap = 0;
         Zp = 0;
         Ep = (Float_t)0.0;
         impactParameter = (Float_t)0.0;
         nCollisions = 0;
         stoppingTime = (Float_t)0.0;
         EBalance = (Float_t)0.0;
         pLongBalance = (Float_t)0.0;
         pTransBalance = (Float_t)0.0;
         nCascadeParticles = 0;
         transparent = false;
         forcedCompoundNucleus = false;
         nucleonAbsorption = false;
         pionAbsorption = false;
         nDecays = 0;
         nBlockedCollisions = 0;
         nBlockedDecays = 0;
         effectiveImpactParameter = (Float_t)0.0;
         deltasInside = false;
         forcedDeltasInside = false;
         forcedDeltasOutside = false;
         clusterDecay = false;
         firstCollisionTime = (Float_t)0.0;
         firstCollisionXSec = (Float_t)0.0;
         firstCollisionSpectatorPosition = (Float_t)0.0;
         firstCollisionSpectatorMomentum = (Float_t)0.0;
         firstCollisionIsElastic = false;
         nReflectionAvatars = 0;
         nCollisionAvatars = 0;
         nDecayAvatars = 0;
         nUnmergedSpectators = 0;
         nEnergyViolationInteraction = 0;
 #ifdef INCL_INVERSE_KINEMATICS
 
 #endif
       }

Field Documentation

Short_t G4INCL::EventInfo::A[maxSizeParticles]

Particle mass number.

Definition at line 151 of file G4INCLEventInfo.hh.

Referenced by G4INCLXXInterface::ApplyYourself(), EventInfo(), G4INCL::Nucleus::fillEventInfo(), and remnantToParticle().

Short_t G4INCL::EventInfo::Ap

Mass number of the projectile nucleus.

Definition at line 227 of file G4INCLEventInfo.hh.

Referenced by G4INCL::Nucleus::getConservationBalance(), and reset().

Short_t G4INCL::EventInfo::ARem[maxSizeRemnants]

Remnant mass number.

Definition at line 195 of file G4INCLEventInfo.hh.

Referenced by G4INCLXXInterface::ApplyYourself(), EventInfo(), G4INCL::Nucleus::fillEventInfo(), and remnantToParticle().

Short_t G4INCL::EventInfo::At

Mass number of the target nucleus.

Definition at line 223 of file G4INCLEventInfo.hh.

Referenced by G4INCL::Nucleus::getConservationBalance(), and reset().

Bool_t G4INCL::EventInfo::clusterDecay

Event involved cluster decay.

Definition at line 269 of file G4INCLEventInfo.hh.

Referenced by reset().

Bool_t G4INCL::EventInfo::deltasInside

Event involved deltas in the nucleus at the end of the cascade.

Definition at line 263 of file G4INCLEventInfo.hh.

Referenced by reset().

Float_t G4INCL::EventInfo::EBalance

Energy-conservation balance [MeV].

Definition at line 239 of file G4INCLEventInfo.hh.

Referenced by reset().

Float_t G4INCL::EventInfo::effectiveImpactParameter

Effective (Coulomb-distorted) impact parameter [fm].

Definition at line 261 of file G4INCLEventInfo.hh.

Referenced by reset().

Float_t G4INCL::EventInfo::EKin[maxSizeParticles]

Particle kinetic energy [MeV].

Definition at line 155 of file G4INCLEventInfo.hh.

Referenced by G4INCLXXInterface::ApplyYourself(), EventInfo(), G4INCL::Nucleus::fillEventInfo(), and remnantToParticle().

Float_t G4INCL::EventInfo::EKinRem[maxSizeRemnants]

Remnant kinetic energy [MeV].

Definition at line 203 of file G4INCLEventInfo.hh.

Referenced by G4INCLXXInterface::ApplyYourself(), EventInfo(), G4INCL::Nucleus::fillEventInfo(), and remnantToParticle().

Float_t G4INCL::EventInfo::emissionTime[maxSizeParticles]

Emission time [fm/c].

Definition at line 172 of file G4INCLEventInfo.hh.

Referenced by EventInfo(), G4INCL::Nucleus::fillEventInfo(), and remnantToParticle().

Float_t G4INCL::EventInfo::Ep

Projectile kinetic energy given as input.

Definition at line 231 of file G4INCLEventInfo.hh.

Referenced by reset().

Float_t G4INCL::EventInfo::EStarRem[maxSizeRemnants]

Remnant excitation energy [MeV].

Definition at line 199 of file G4INCLEventInfo.hh.

Referenced by G4INCLXXInterface::ApplyYourself(), EventInfo(), and G4INCL::Nucleus::fillEventInfo().

G4ThreadLocal Int_t G4INCL::EventInfo::eventNumber = 0

static

Number of the event.

Definition at line 140 of file G4INCLEventInfo.hh.

Bool_t G4INCL::EventInfo::firstCollisionIsElastic

True if the first collision was elastic.

Definition at line 279 of file G4INCLEventInfo.hh.

Referenced by G4INCL::Nucleus::fillEventInfo(), and reset().

Float_t G4INCL::EventInfo::firstCollisionSpectatorMomentum

Momentum of the spectator on the first collision (fm)

Definition at line 277 of file G4INCLEventInfo.hh.

Referenced by G4INCL::Nucleus::fillEventInfo(), and reset().

Float_t G4INCL::EventInfo::firstCollisionSpectatorPosition

Position of the spectator on the first collision (fm)

Definition at line 275 of file G4INCLEventInfo.hh.

Referenced by G4INCL::Nucleus::fillEventInfo(), and reset().

Float_t G4INCL::EventInfo::firstCollisionTime

Time of the first collision [fm/c].

Definition at line 271 of file G4INCLEventInfo.hh.

Referenced by G4INCL::Nucleus::fillEventInfo(), and reset().

Float_t G4INCL::EventInfo::firstCollisionXSec

Cross section of the first collision (mb)

Definition at line 273 of file G4INCLEventInfo.hh.

Referenced by G4INCL::Nucleus::fillEventInfo(), and reset().

Bool_t G4INCL::EventInfo::forcedCompoundNucleus

True if the event is a forced CN.

Definition at line 249 of file G4INCLEventInfo.hh.

Referenced by G4INCL::Nucleus::fillEventInfo(), and reset().

Bool_t G4INCL::EventInfo::forcedDeltasInside

Event involved forced delta decays inside the nucleus.

Definition at line 265 of file G4INCLEventInfo.hh.

Referenced by reset().

Bool_t G4INCL::EventInfo::forcedDeltasOutside

Event involved forced delta decays outside the nucleus.

Definition at line 267 of file G4INCLEventInfo.hh.

Referenced by reset().

std::vector<std::string> G4INCL::EventInfo::history

History of the particle.

Condensed information about the de-excitation chain of a particle. For cascade particles, it is just an empty string. For particles arising from the de-excitation of a cascade remnant, it is a string of characters. Each character represents one or more identical steps in the de-excitation process. The currently defined possible character values and their meanings are the following:

e: evaporation product E: evaporation residue m: multifragmentation a: light partner in asymmetric fission or IMF emission A: heavy partner in asymmetric fission or IMF emission f: light partner in fission F: heavy partner in fission s: saddle-to-scission emission n: non-statistical emission (decay)

Definition at line 191 of file G4INCLEventInfo.hh.

Referenced by G4INCL::Nucleus::fillEventInfo(), remnantToParticle(), and reset().

Float_t G4INCL::EventInfo::impactParameter

Impact parameter [fm].

Definition at line 233 of file G4INCLEventInfo.hh.

Referenced by reset().

Float_t G4INCL::EventInfo::JRem[maxSizeRemnants]

Remnant spin [ $\hbar$ ].

Definition at line 201 of file G4INCLEventInfo.hh.

Referenced by EventInfo(), and G4INCL::Nucleus::fillEventInfo().

Float_t G4INCL::EventInfo::jxRem[maxSizeRemnants]

Remnant angular momentum, x component [ $\hbar$ ].