- Timestamp:
- Apr 6, 2009, 12:30:29 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/hadronic/stopping/src/G4PiMinusStopAbsorption.cc
r819 r962 38 38 #include "globals.hh" 39 39 #include "Randomize.hh" 40 #include "G4NucleiPropertiesTable.hh"41 40 #include "G4NucleiProperties.hh" 42 41 #include "G4ParticleTypes.hh" … … 94 93 } 95 94 96 G4double binding = G4NucleiProperties Table::GetBindingEnergy(static_cast<G4int>(_nucleusZ) ,static_cast<G4int>(_nucleusA)) / _nucleusA;95 G4double binding = G4NucleiProperties::GetBindingEnergy(static_cast<G4int>(_nucleusA) ,static_cast<G4int>(_nucleusZ)) / _nucleusA; 97 96 G4double mass = G4NucleiProperties::GetNuclearMass(static_cast<G4int>(newA),static_cast<G4int>(newZ)); 98 97 … … 185 184 } 186 185 187 G4double productBinding = (G4NucleiProperties Table::GetBindingEnergy(static_cast<G4int>(_nucleusZ),static_cast<G4int>(_nucleusA)) / _nucleusA) * nAbsorptionProducts;186 G4double productBinding = (G4NucleiProperties::GetBindingEnergy(static_cast<G4int>(_nucleusA),static_cast<G4int>(_nucleusZ)) / _nucleusA) * nAbsorptionProducts; 188 187 G4double mass = G4NucleiProperties::GetNuclearMass(_nucleusA - (nP + nN),_nucleusZ - nP); 189 188 G4double pNucleus = pProducts.mag(); 190 189 G4double eNucleus = std::sqrt(pNucleus*pNucleus + mass*mass); 191 190 G4double tNucleus = eNucleus - mass; 192 G4double temp = G4NucleiPropertiesTable::GetBindingEnergy(static_cast<G4int>(_nucleusZ - nP),static_cast<G4int>(_nucleusA - (nP + nN))) - 193 G4NucleiPropertiesTable::GetBindingEnergy(static_cast<G4int>(_nucleusZ),static_cast<G4int>(_nucleusA)); 191 G4double temp = 192 G4NucleiProperties::GetBindingEnergy(static_cast<G4int>(_nucleusA - (nP + nN)),static_cast<G4int>(_nucleusZ - nP)) - 193 G4NucleiProperties::GetBindingEnergy(static_cast<G4int>(_nucleusA),static_cast<G4int>(_nucleusZ)); 194 194 energy = productEnergy + productBinding + tNucleus; 195 195
Note: See TracChangeset
for help on using the changeset viewer.