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) {
414 for (i = 1; i <= 12; ++i) {
496 en = ek + (rmnve1 + rmnve2) /
G4float(2.);
497 r__1 = en * en - rmnve1 * rmnve2;
498 pcm = r__1 > 0 ? std::sqrt(r__1) : 0;
540 black = (targ *
G4float(1.25) +
552 for (i = 1; i <= nbl; ++i) {
569 ekin1 = tex - (ekin2 - ekin1);
597 for (i = 1; i <= nt; ++i) {
627 for (i = 1; i <= nbl; ++i) {
644 ekin1 = tex - (ekin2 - ekin1);
671 for (i = 3; i <= nt; ++i) {
709 cfa =
G4float(.13043478260869565);
738 evapEnergy1 *= ran1 * gfa +
G4float(1.);
739 if (evapEnergy1 <
G4float(0.)) {
748 while ((ret_val = evapEnergy1 +
evapEnergy3) >= ek1) {
G4double condition(const G4ErrorSymMatrix &m)
static G4Pow * GetInstance()
void DeRegisterExtraProcess(G4VProcess *)
G4double GetMeanLifeTime(const G4Track &, G4ForceCondition *)
G4ParticleDefinition * pdefNeutron
void GenerateSecondaries()
G4GHEKinematicsVector * gkin
void SetParticleDef(G4ParticleDefinition *c)
void SetMomentum(const G4ThreeVector &momentum)
G4ParticleDefinition * pdefGamma
G4double powN(G4double x, G4int n) const
std::vector< G4Element * > G4ElementVector
G4AntiNeutronAnnihilationAtRest(const G4AntiNeutronAnnihilationAtRest &)
G4ParticleDefinition * pdefTriton
CLHEP::Hep3Vector G4ThreeVector
static G4HadronicProcessStore * Instance()
const G4DynamicParticle * GetDynamicParticle() const
~G4AntiNeutronAnnihilationAtRest()
G4ParticleDefinition * GetParticleDef()
G4ParticleDefinition * pdefPionZero
const G4String & GetName() const
const G4ThreeVector & GetPosition() const
void Poisso(G4float, G4int *)
void DumpInfo(G4int mode=0) const
G4GHEKinematicsVector * GetSecondaryKinematics()
G4double theNumberOfInteractionLengthLeft
G4ParticleDefinition * pdefAntiNeutron
void SetTouchableHandle(const G4TouchableHandle &apValue)
G4ParticleDefinition * pdefAlpha
void BuildPhysicsTable(const G4ParticleDefinition &)
virtual void ResetNumberOfInteractionLengthLeft()
#define G4HadronicDeprecate(name)
const G4ElementVector * GetElementVector() const
void ProposeLocalEnergyDeposit(G4double anEnergyPart)
void SetKineticEnergyAndUpdate(G4double ekin)
G4GHEKinematicsVector * pv
G4VParticleChange * AtRestDoIt(const G4Track &, const G4Step &)
void RegisterParticleForExtraProcess(G4VProcess *, const G4ParticleDefinition *)
G4GLOB_DLL std::ostream G4cout
void AntiNeutronAnnihilation(G4int *)
G4int GetNumberOfSecondaries()
G4double AtRestGetPhysicalInteractionLength(const G4Track &, G4ForceCondition *)
G4double currentInteractionLength
void SetProcessSubType(G4int)
G4ParticleDefinition * pdefDeuteron
G4double GetGlobalTime() const
const G4String & GetProcessName() const
const G4double * GetAtomicNumDensityVector() const
const G4TouchableHandle & GetTouchableHandle() const
G4Material * GetMaterial() const
void SetEnergyAndUpdate(G4double e)
G4ParticleDefinition * pdefPionPlus
void RegisterExtraProcess(G4VProcess *)
G4bool IsApplicable(const G4ParticleDefinition &)
G4double G4Log(G4double x)
void PreparePhysicsTable(const G4ParticleDefinition &)
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
virtual void Initialize(const G4Track &)
void SetEnergy(G4double e)
void SetMomentumAndUpdate(G4ParticleMomentum mom)
void SetNumberOfSecondaries(G4int totSecondaries)
G4GHEKinematicsVector result
G4ParticleDefinition * pdefProton
const G4double x[NPOINTSGL]
G4ParticleChange aParticleChange
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
void SetMass(G4double mas)
void AddSecondary(G4Track *aSecondary)
size_t GetNumberOfElements() const
G4GHEKinematicsVector * eve
void ProposeTrackStatus(G4TrackStatus status)
void SetDefinition(const G4ParticleDefinition *aParticleDefinition)
G4ThreeVector G4ParticleMomentum
G4ParticleDefinition * pdefPionMinus
void PrintInfo(const G4ParticleDefinition *)