42 #define MAX_SECONDARIES 100
140 G4cout <<
"G4AntiNeutronAnnihilationAtRestProcess::AtRestGetPhysicalInteractionLength ";
177 for (
G4int i1=0; i1 < numberOfElements; i1++ )
179 normalization += theAtomicNumberDensity[i1] ;
184 for (
G4int i2=0; i2 < numberOfElements; i2++ )
186 runningSum += theAtomicNumberDensity[i2];
188 if (random<=runningSum)
194 if (random>runningSum)
201 G4cout <<
"G4AntiNeutronAnnihilationAtRest::AtRestDoIt is invoked " <<
G4endl;
220 G4Track* aNewTrack =
new G4Track( aNewParticle, localtime*
s, position );
281 for (l = 1; l <=
ntot; ++l) {
322 ran1 = xav + ran1 * std::sqrt(xav);
336 for (i = 1; i <= fivex; ++i) {
397 for (i = 2; i <= j; ++i) {
429 G4float ran1, ran2, ekin, tkin;
486 en = ek + (rmnve1 + rmnve2) /
G4float(2.);
487 r__1 = en * en - rmnve1 * rmnve2;
488 pcm = r__1 > 0 ? std::sqrt(r__1) : 0;
528 black = (targ *
G4float(1.25) +
540 for (i = 1; i <= nbl; ++i) {
557 ekin1 = tex - (ekin2 - ekin1);
585 for (i = 1; i <= nt; ++i) {
615 for (i = 1; i <= nbl; ++i) {
632 ekin1 = tex - (ekin2 - ekin1);
659 for (i = 3; i <= nt; ++i) {
673 G4float cfa, gfa, ran1, ran2, ekin1, atno3;
697 cfa =
G4float(.13043478260869565);
726 evapEnergy1 *= ran1 * gfa +
G4float(1.);
727 if (evapEnergy1 < 0.0
f) {
736 while ((ret_val = evapEnergy1 +
evapEnergy3) >= ek1) {
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
void Poisso(G4float, G4int *)
void SetEnergy(G4double e)
G4ParticleDefinition * pdefTriton
void SetMass(G4double mas)
virtual void ResetNumberOfInteractionLengthLeft()
G4ParticleDefinition * pdefNeutron
void RegisterExtraProcess(G4VProcess *)
G4ParticleDefinition * pdefAntiNeutron
G4GHEKinematicsVector * gkin
G4AntiNeutronAnnihilationAtRest(const G4AntiNeutronAnnihilationAtRest &)
void SetKineticEnergyAndUpdate(G4double ekin)
void SetEnergyAndUpdate(G4double e)
const G4ElementVector * GetElementVector() const
const G4TouchableHandle & GetTouchableHandle() const
G4GHEKinematicsVector result
G4double condition(const G4ErrorSymMatrix &m)
G4double GetPDGMass() const
G4GHEKinematicsVector * pv
G4double G4Log(G4double x)
G4double powN(G4double x, G4int n) const
void AddSecondary(G4Track *aSecondary)
G4ParticleChange aParticleChange
G4double currentInteractionLength
const G4String & GetName() const
~G4AntiNeutronAnnihilationAtRest()
static G4Pow * GetInstance()
G4double GetMeanLifeTime(const G4Track &, G4ForceCondition *)
void SetParticleDef(G4ParticleDefinition *c)
virtual void Initialize(const G4Track &)
G4int GetNumberOfSecondaries()
G4ParticleDefinition * pdefAlpha
G4GHEKinematicsVector * GetSecondaryKinematics()
void BuildPhysicsTable(const G4ParticleDefinition &)
G4double GetGlobalTime() const
void SetTouchableHandle(const G4TouchableHandle &apValue)
G4ParticleDefinition * pdefGamma
const G4String & GetProcessName() const
const G4ThreeVector & GetPosition() const
void PreparePhysicsTable(const G4ParticleDefinition &)
G4double AtRestGetPhysicalInteractionLength(const G4Track &, G4ForceCondition *)
G4Material * GetMaterial() const
G4GHEKinematicsVector * eve
void PrintInfo(const G4ParticleDefinition *)
static G4HadronicProcessStore * Instance()
void SetMomentumAndUpdate(G4ParticleMomentum mom)
void DeRegisterExtraProcess(G4VProcess *)
G4ParticleDefinition * pdefPionZero
std::vector< G4Element * > G4ElementVector
void GenerateSecondaries()
void SetProcessSubType(G4int)
void SetMomentum(const G4ThreeVector &momentum)
void SetDefinition(const G4ParticleDefinition *aParticleDefinition)
G4GLOB_DLL std::ostream G4cout
G4ParticleDefinition * pdefPionMinus
#define G4HadronicDeprecate(name)
const G4double * GetAtomicNumDensityVector() const
G4ParticleDefinition * GetParticleDef()
G4ParticleDefinition * pdefPionPlus
G4double theNumberOfInteractionLengthLeft
void ProposeLocalEnergyDeposit(G4double anEnergyPart)
G4ParticleDefinition * pdefDeuteron
G4ParticleDefinition * pdefProton
void DumpInfo(G4int mode=0) const
void ProposeTrackStatus(G4TrackStatus status)
static constexpr double GeV
const G4DynamicParticle * GetDynamicParticle() const
G4VParticleChange * AtRestDoIt(const G4Track &, const G4Step &)
void RegisterParticleForExtraProcess(G4VProcess *, const G4ParticleDefinition *)
size_t GetNumberOfElements() const
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
void AntiNeutronAnnihilation(G4int *)
G4bool IsApplicable(const G4ParticleDefinition &)
void SetNumberOfSecondaries(G4int totSecondaries)