Changeset 1315 for trunk/source/processes/hadronic/models/chiral_inv_phase_space/cross_sections/src/G4QPhotonNuclearCrossSection.cc

Timestamp:

Jun 18, 2010, 11:42:07 AM (14 years ago)

Author:

garnier

Message:

update geant4-09-04-beta-cand-01 interfaces-V09-03-09 vis-V09-03-08

File:

• trunk/source/processes/hadronic/models/chiral_inv_phase_space/cross_sections/src/G4QPhotonNuclearCrossSection.cc (modified) (2 diffs)

trunk/source/processes/hadronic/models/chiral_inv_phase_space/cross_sections/src/G4QPhotonNuclearCrossSection.cc

@@ -26 +26 @@
 //
 // The lust update: M.V. Kossov, CERN/ITEP(Moscow) 17-June-02
-// GEANT4 tag $Name: geant4-09-03 $
+// GEANT4 tag $Name: geant4-09-04-beta-cand-01 $
 //
 //
@@ -276 +276 @@
   //G4double mT= G4QPDGCode(111).GetNuclMass(Z,N,0);
   G4double mT= 0.;
-  if(G4NucleiProperties::IsInStableTable(A,Z))
+  if(Z&&G4NucleiProperties::IsInStableTable(A,Z))
                            mT = G4NucleiProperties::GetNuclearMass(A,Z)/MeV;
+  else return 0.;       // If it is not in the Table of Stable Nuclei, then the Threshold=0
   G4double mP= infEn;
-  if(Z&&G4NucleiProperties::IsInStableTable(A-1,Z-1))
+  if(A>1&&Z&&G4NucleiProperties::IsInStableTable(A-1,Z-1))
           mP = G4NucleiProperties::GetNuclearMass(A-1,Z-1)/MeV;// ResNucMass for a proton
 
   G4double mN= infEn;
-  if(N&&G4NucleiProperties::IsInStableTable(A-1,Z))
+  if(A>1&&G4NucleiProperties::IsInStableTable(A-1,Z))
           mN = G4NucleiProperties::GetNuclearMass(A-1,Z)/MeV;  // ResNucMass for a neutron
 
   G4double mA= infEn;
-  if(N>1&&Z>1&&G4NucleiProperties::IsInStableTable(A-4,Z-2))
+  if(A>4&&Z>1&&G4NucleiProperties::IsInStableTable(A-4,Z-2))
           mA=G4NucleiProperties::GetNuclearMass(A-4,Z-2)/MeV;  // ResNucMass for an alpha