54 delete originalTarget;
64 modifiedOriginal = *originalIncident;
83 p = std::sqrt(
std::abs((et-amas)*(et+amas)) );
92 delete originalTarget;
98 targetParticle = *originalTarget;
101 G4bool incidentHasChanged =
false;
102 G4bool targetHasChanged =
false;
103 G4bool quasiElastic =
false;
109 incidentHasChanged, targetHasChanged);
112 originalIncident, originalTarget, modifiedOriginal,
113 targetNucleus, currentParticle, targetParticle,
114 incidentHasChanged, targetHasChanged, quasiElastic);
117 currentParticle, targetParticle,
118 incidentHasChanged );
120 delete originalTarget;
136 massVec[0] = targetNucleus.
AtomicMass( A+1.0, Z+1.0 );
139 massVec[1] = targetNucleus.
AtomicMass( A , Z+1.0 );
140 massVec[2] = theAtomicMass;
142 if (A > 1.0 && A-1.0 > Z)
143 massVec[3] = targetNucleus.
AtomicMass( A-1.0, Z );
145 if (A > 2.0 && A-2.0 > Z)
146 massVec[4] = targetNucleus.
AtomicMass( A-2.0, Z );
148 if (A > 3.0 && Z > 1.0 && A-3.0 > Z-1.0)
149 massVec[5] = targetNucleus.
AtomicMass( A-3.0, Z-1.0 );
151 if (A > 1.0 && A-1.0 > Z+1.0)
152 massVec[6] = targetNucleus.
AtomicMass( A-1.0, Z+1.0 );
153 massVec[7] = massVec[3];
155 if (A > 1.0 && Z > 1.0)
156 massVec[8] = targetNucleus.
AtomicMass( A-1.0, Z-1.0 );
163 targetNucleus, theAtomicMass, massVec );
169 for(
G4int i=0; i<vecLen; ++i )
191 G4bool& incidentHasChanged,
198 std::vector<G4int> fsTypes;
219 incidentHasChanged =
true;
223 targetHasChanged =
true;
225 }
else if (part2 >
neu && part2 <
xi0) {
226 targetHasChanged =
true;
230 targetHasChanged =
true;
231 incidentHasChanged =
true;
248 targetHasChanged =
true;
249 }
else if (part2 ==
neu) {
251 incidentHasChanged =
true;
252 targetHasChanged =
true;
257 targetHasChanged =
true;
261 incidentHasChanged =
true;
262 if (part2 >
neu && part2 <
xi0) targetHasChanged =
true;
272 fsTypes.erase(fsTypes.begin());
273 fsTypes.erase(fsTypes.begin());
278 for(
G4int i=0; i < mult-2; ++i ) {
279 partType = fsTypes[i];
288 CheckQnums(vec, vecLen, currentParticle, targetParticle,
289 testCharge, testBaryon, testStrange);