- Timestamp:
- Nov 5, 2010, 3:45:55 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/hadronic/models/binary_cascade/src/G4PionPlusField.cc
r962 r1340 83 83 G4double pionPlusMass = G4PionPlus::PionPlus()->GetPDGMass(); 84 84 85 G4 doubleA = theNucleus->GetMassNumber();86 G4 doubleZ = theNucleus->GetCharge();87 G4double bindingEnergy = G4NucleiProperties::GetBindingEnergy( G4lrint(A), G4lrint(Z));85 G4int A = theNucleus->GetMassNumber(); 86 G4int Z = theNucleus->GetCharge(); 87 G4double bindingEnergy = G4NucleiProperties::GetBindingEnergy(A, Z); 88 88 G4double nucleusMass = Z*proton_mass_c2+(A-Z)*neutron_mass_c2+bindingEnergy; 89 89 G4double reducedMass = pionPlusMass*nucleusMass/(pionPlusMass+nucleusMass); … … 98 98 G4double G4PionPlusField::GetBarrier() 99 99 { 100 G4 doubleA = theNucleus->GetMassNumber();101 G4 doubleZ = theNucleus->GetCharge();100 G4int A = theNucleus->GetMassNumber(); 101 G4int Z = theNucleus->GetCharge(); 102 102 G4double coulombBarrier = (1.44/1.14) * MeV * Z / (1.0 + std::pow(A,1./3.)); 103 103 return coulombBarrier;
Note: See TracChangeset
for help on using the changeset viewer.