62 currentParticle.SetKineticEnergy( ek );
64 G4double p = std::sqrt( std::abs((et-amas)*(et+amas)) );
65 G4double pp = currentParticle.GetMomentum().mag();
69 currentParticle.SetMomentum( momentum * (p/pp) );
76 currentParticle.SetKineticEnergy( ek );
78 p = std::sqrt( std::abs((et-amas)*(et+amas)) );
79 pp = currentParticle.GetMomentum().mag();
83 currentParticle.SetMomentum( momentum * (p/pp) );
89 targetParticle.SetSide( -1 );
90 G4bool incidentHasChanged =
false;
91 G4bool targetHasChanged =
false;
92 G4bool quasiElastic =
false;
98 if( currentParticle.GetKineticEnergy() > cutOff )
100 incidentHasChanged, targetHasChanged);
103 originalIncident, originalTarget, modifiedOriginal,
104 targetNucleus, currentParticle, targetParticle,
105 incidentHasChanged, targetHasChanged, quasiElastic );
108 currentParticle, targetParticle,
109 incidentHasChanged );
111 delete originalTarget;
127 G4bool& incidentHasChanged,
134 std::vector<G4int> fsTypes;
145 partType = fsTypes[0];
146 if (partType !=
pro) {
147 targetHasChanged =
true;
158 partType = fsTypes[0];
159 if (partType !=
neu) {
160 targetHasChanged =
true;
171 fsTypes.erase(fsTypes.begin());
176 for(
G4int i=0; i < mult-1; ++i ) {
177 partType = fsTypes[i];
178 if (partType ==
pip) {
184 incidentHasChanged =
true;
186 partType = fsTypes[choose];
189 fsTypes.erase(fsTypes.begin()+choose);
196 for(
G4int i=0; i < mult-2; ++i ) {
197 partType = fsTypes[i];
210 CheckQnums(vec, vecLen, currentParticle, targetParticle,
211 testCharge, testBaryon, testStrange);
std::vector< G4int > GetFSPartTypesForPipN(G4int mult, G4double KE) const
void Initialize(G4int items)
void SetUpChange(G4FastVector< G4ReactionProduct, 256 > &vec, G4int &vecLen, G4ReactionProduct ¤tParticle, G4ReactionProduct &targetParticle, G4bool &incidentHasChanged)
G4double Cinema(G4double kineticEnergy)
void SetMomentumChange(const G4ThreeVector &aV)
void SetDefinition(const G4ParticleDefinition *aParticleDefinition)
void SetMomentum(const G4double x, const G4double y, const G4double z)
void InitialCollision(G4FastVector< G4ReactionProduct, 256 > &vec, G4int &vecLen, G4ReactionProduct ¤tParticle, G4ReactionProduct &targetParticle, G4bool &incidentHasChanged, G4bool &targetHasChanged)
const G4ParticleDefinition * GetDefinition() const
std::vector< G4int > GetFSPartTypesForPipP(G4int mult, G4double KE) 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)
G4double GetPDGMass() const
G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
void CheckQnums(G4FastVector< G4ReactionProduct, 256 > &vec, G4int &vecLen, G4ReactionProduct ¤tParticle, G4ReactionProduct &targetParticle, G4double Q, G4double B, G4double S)
G4DynamicParticle * ReturnTargetParticle() const
void SetEnergyChange(G4double anEnergy)
G4double GetKineticEnergy() const
G4ParticleDefinition * GetDefinition() const
G4int GetMultiplicityT12(G4double KE) const
G4int GetMultiplicityT32(G4double KE) const
void SetSide(const G4int sid)
const G4LorentzVector & Get4Momentum() const
G4double GetKineticEnergy() const
G4ParticleDefinition * particleDef[18]
void SetElement(G4int anIndex, Type *anElement)
const G4ParticleDefinition * GetDefinition() const
G4double EvaporationEffects(G4double kineticEnergy)
static constexpr double GeV
G4HadFinalState theParticleChange
void SetMayBeKilled(const G4bool f)
void SetStatusChange(G4HadFinalStateStatus aS)