Ignore:
Timestamp:
Apr 6, 2009, 12:30:29 PM (15 years ago)
Author:
garnier
Message:

update processes

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/source/processes/hadronic/models/chiral_inv_phase_space/interface/src/G4QElectronNuclearCrossSection.cc

    r819 r962  
    2525//
    2626//
    27 // $Id: G4QElectronNuclearCrossSection.cc,v 1.11.10.1 2008/04/23 14:57:22 gcosmo Exp $
    28 // GEANT4 tag $Name: geant4-09-01-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 $
    2929//
    3030//
     
    267267  //G4double mT= G4QPDGCode(111).GetNuclMass(Z,N,0);
    268268  G4double mT= 0.;
    269   if(G4NucleiPropertiesTable::IsInTable(Z,A))
    270                                             mT=G4NucleiProperties::GetNuclearMass(A,Z)/MeV;
     269  if(G4NucleiProperties::IsInStableTable(A,Z))
     270                               mT = G4NucleiProperties::GetNuclearMass(A,Z)/MeV;
    271271  else return 0.;       // If it is not in the Table of Stable Nuclei, then the Threshold=0
    272272  // ---------
    273273  G4double mP= infEn;
    274274  //if(Z) mP= G4QPDGCode(111).GetNuclMass(Z-1,N,0);
    275   if(Z&&G4NucleiPropertiesTable::IsInTable(Z-1,A-1))
    276             mP=G4NucleiProperties::GetNuclearMass(A-1.,Z-1.)/MeV; // ResNucMass for a proton
     275  if(Z&&G4NucleiProperties::IsInStableTable(A-1,Z-1))
     276            mP = G4NucleiProperties::GetNuclearMass(A-1.,Z-1.)/MeV; // ResNucMass for a proton
    277277  G4double mN= infEn;
    278278  //if(N) mN= G4QPDGCode(111).GetNuclMass(Z,N-1,0);
    279   if(N&&G4NucleiPropertiesTable::IsInTable(Z,A-1))
    280             mN=G4NucleiProperties::GetNuclearMass(A-1.,Z-0.)/MeV;  // ResNucMass for a neutron
     279  if(N&&G4NucleiProperties::IsInStableTable(A-1,Z))
     280            mN = G4NucleiProperties::GetNuclearMass(A-1.,Z-0.)/MeV;  // ResNucMass for a neutron
    281281
    282282  G4double mA= infEn;
    283   if(N>1&&Z>1&&G4NucleiPropertiesTable::IsInTable(Z-2,A-4))
    284             mA=G4NucleiProperties::GetNuclearMass(A-4.,Z-2.)/MeV; // ResNucMass for an alpha
     283  if(N>1&&Z>1&&G4NucleiProperties::IsInStableTable(A-4,Z-2))
     284            mA = G4NucleiProperties::GetNuclearMass(A-4.,Z-2.)/MeV; // ResNucMass for an alpha
    285285
    286286  G4double dP= mP +mProt - mT;
Note: See TracChangeset for help on using the changeset viewer.