Changeset 1007 for trunk/source/processes/hadronic/models/cascade/cascade/src/G4IntraNucleiCascader.cc
- Timestamp:
- Apr 20, 2009, 5:54:05 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/hadronic/models/cascade/cascade/src/G4IntraNucleiCascader.cc
r962 r1007 158 158 theElementaryParticleCollider); 159 159 if (verboseLevel > 2) { 160 G4cout << " New particles from G4NucleiModel::generateParticleFate " 161 << new_cascad_particles.size() << G4endl; 160 G4cout << " ew particles " << new_cascad_particles.size() << G4endl; 162 161 } 163 162 … … 219 218 } 220 219 221 } else if (model.highEnergyInteraction() ) { // high energy projectile-nucleon interaction 220 } else { // interaction 221 222 222 cascad_particles.pop_back(); 223 for (G4int i = 0; i < G4int(new_cascad_particles.size()); i++) { 224 // send secondaries straight to output - no further cascade 225 // for now, don't check for Coulomb barrier? 226 output_particles.push_back(new_cascad_particles[i].getParticle() ); 227 } 228 std::pair<G4int, G4int> holes = model.getTypesOfNucleonsInvolved(); 229 theExitonConfiguration.incrementHoles(holes.first); 230 if (holes.second > 0) theExitonConfiguration.incrementHoles(holes.second); 231 232 } else { // medium energy projectile-nucleon interaction 233 cascad_particles.pop_back(); 223 234 224 for (G4int i = 0; i < G4int(new_cascad_particles.size()); i++) 235 225 cascad_particles.push_back(new_cascad_particles[i]); 226 236 227 std::pair<G4int, G4int> holes = model.getTypesOfNucleonsInvolved(); 228 237 229 theExitonConfiguration.incrementHoles(holes.first); 230 238 231 if (holes.second > 0) theExitonConfiguration.incrementHoles(holes.second); 239 } 240 } 232 233 }; 234 }; 241 235 242 236 // Cascade is finished. Check if it's OK.
Note: See TracChangeset
for help on using the changeset viewer.