36{
40 {
45 }
46
47
48
49
51
53 {
54 SlowProton( originalIncident, targetNucleus );
55 delete originalTarget;
57 }
58
59
60
61
65 modifiedOriginal = *originalIncident;
66
68 ek += tkin;
71 G4double p = std::sqrt( std::abs((et-amas)*(et+amas)) );
73 if (pp > 0.0) {
76 }
77
78
79
81 ek -= tkin;
83 et = ek + amas;
84 p = std::sqrt( std::abs((et-amas)*(et+amas)) );
86 if (pp > 0.0) {
89 }
92 SlowProton( originalIncident, targetNucleus );
93 delete originalTarget;
95 }
96
99 targetParticle = *originalTarget;
102 G4bool incidentHasChanged =
false;
103 G4bool targetHasChanged =
false;
104 G4bool quasiElastic =
false;
108
109 InitialCollision(vec, vecLen, currentParticle, targetParticle,
110 incidentHasChanged, targetHasChanged);
111
113 originalIncident, originalTarget, modifiedOriginal,
114 targetNucleus, currentParticle, targetParticle,
115 incidentHasChanged, targetHasChanged, quasiElastic);
116
118 currentParticle, targetParticle,
119 incidentHasChanged );
120
121 delete originalTarget;
123}
void Initialize(G4int items)
void SetStatusChange(G4HadFinalStateStatus aS)
void SetEnergyChange(G4double anEnergy)
void SetMomentumChange(const G4ThreeVector &aV)
const G4ParticleDefinition * GetDefinition() const
G4double GetKineticEnergy() const
const G4LorentzVector & Get4Momentum() const
G4HadFinalState theParticleChange
G4double EvaporationEffects(G4double kineticEnergy)
G4double Cinema(G4double kineticEnergy)
G4DynamicParticle * ReturnTargetParticle() const
G4double GetPDGMass() const
void CalculateMomenta(G4FastVector< G4ReactionProduct, 256 > &vec, G4int &vecLen, const G4HadProjectile *originalIncident, const G4DynamicParticle *originalTarget, G4ReactionProduct &modifiedOriginal, G4Nucleus &targetNucleus, G4ReactionProduct ¤tParticle, G4ReactionProduct &targetParticle, G4bool &incidentHasChanged, G4bool &targetHasChanged, G4bool quasiElastic)
void SetUpChange(G4FastVector< G4ReactionProduct, 256 > &vec, G4int &vecLen, G4ReactionProduct ¤tParticle, G4ReactionProduct &targetParticle, G4bool &incidentHasChanged)
void SetMomentum(const G4double x, const G4double y, const G4double z)
G4double GetKineticEnergy() const
G4ThreeVector GetMomentum() const
void SetSide(const G4int sid)
void SetKineticEnergy(const G4double en)