Changeset 962 for trunk/source/processes/hadronic/models/chiral_inv_phase_space/interface/src/G4QElectronNuclearCrossSection.cc
- Timestamp:
- Apr 6, 2009, 12:30:29 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/hadronic/models/chiral_inv_phase_space/interface/src/G4QElectronNuclearCrossSection.cc
r819 r962 25 25 // 26 26 // 27 // $Id: G4QElectronNuclearCrossSection.cc,v 1.1 1.10.1 2008/04/23 14:57:22 gcosmoExp $28 // GEANT4 tag $Name: geant4-09-0 1-patch-02 $27 // $Id: G4QElectronNuclearCrossSection.cc,v 1.13 2008/10/24 19:25:41 dennis Exp $ 28 // GEANT4 tag $Name: geant4-09-02-ref-02 $ 29 29 // 30 30 // … … 267 267 //G4double mT= G4QPDGCode(111).GetNuclMass(Z,N,0); 268 268 G4double mT= 0.; 269 if(G4NucleiProperties Table::IsInTable(Z,A))270 mT=G4NucleiProperties::GetNuclearMass(A,Z)/MeV;269 if(G4NucleiProperties::IsInStableTable(A,Z)) 270 mT = G4NucleiProperties::GetNuclearMass(A,Z)/MeV; 271 271 else return 0.; // If it is not in the Table of Stable Nuclei, then the Threshold=0 272 272 // --------- 273 273 G4double mP= infEn; 274 274 //if(Z) mP= G4QPDGCode(111).GetNuclMass(Z-1,N,0); 275 if(Z&&G4NucleiProperties Table::IsInTable(Z-1,A-1))276 mP =G4NucleiProperties::GetNuclearMass(A-1.,Z-1.)/MeV; // ResNucMass for a proton275 if(Z&&G4NucleiProperties::IsInStableTable(A-1,Z-1)) 276 mP = G4NucleiProperties::GetNuclearMass(A-1.,Z-1.)/MeV; // ResNucMass for a proton 277 277 G4double mN= infEn; 278 278 //if(N) mN= G4QPDGCode(111).GetNuclMass(Z,N-1,0); 279 if(N&&G4NucleiProperties Table::IsInTable(Z,A-1))280 mN =G4NucleiProperties::GetNuclearMass(A-1.,Z-0.)/MeV; // ResNucMass for a neutron279 if(N&&G4NucleiProperties::IsInStableTable(A-1,Z)) 280 mN = G4NucleiProperties::GetNuclearMass(A-1.,Z-0.)/MeV; // ResNucMass for a neutron 281 281 282 282 G4double mA= infEn; 283 if(N>1&&Z>1&&G4NucleiProperties Table::IsInTable(Z-2,A-4))284 mA =G4NucleiProperties::GetNuclearMass(A-4.,Z-2.)/MeV; // ResNucMass for an alpha283 if(N>1&&Z>1&&G4NucleiProperties::IsInStableTable(A-4,Z-2)) 284 mA = G4NucleiProperties::GetNuclearMass(A-4.,Z-2.)/MeV; // ResNucMass for an alpha 285 285 286 286 G4double dP= mP +mProt - mT;
Note: See TracChangeset
for help on using the changeset viewer.