Changeset 1347 for trunk/source/processes/hadronic/cross_sections/src/G4EMDissociationCrossSection.cc

Timestamp:

Dec 22, 2010, 3:52:27 PM (13 years ago)

Author:

garnier

Message:

geant4 tag 9.4

File:

• trunk/source/processes/hadronic/cross_sections/src/G4EMDissociationCrossSection.cc (modified) (5 diffs)

trunk/source/processes/hadronic/cross_sections/src/G4EMDissociationCrossSection.cc

@@ -56 +56 @@
 // Beta release
 //
-// 30. May 2005, J.P. Wellisch removed a compilation warning on gcc 3.4 for 
+// 30 May 2005, J.P. Wellisch removed a compilation warning on gcc 3.4 for 
 //               geant4 7.1.
+// 09 November 2010, V.Ivanchenko make class applicable for Hydrogen but 
+//                   set cross section for Hydrogen to zero  
 //
 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
@@ -97 +99 @@
 G4bool
 G4EMDissociationCrossSection::IsIsoApplicable(const G4DynamicParticle* theDynamicParticle,
-                                              G4int /*ZZ*/, G4int AA)
+                                              G4int /*ZZ*/, G4int/* AA*/)
 {
 //
@@ -107 +109 @@
 //
   if (G4ParticleTable::GetParticleTable()->GetIonTable()->
-    IsIon(theDynamicParticle->GetDefinition()) && AA > 1)
+      IsIon(theDynamicParticle->GetDefinition()) /*&& AA > 1*/)
     return true;
   else
@@ -127 +129 @@
   G4int nIso = theElement->GetNumberOfIsotopes();
   G4double crossSection = 0;
+
+  // VI protection for Hydrogen
+  if(theElement->GetZ() < 1.5) { return crossSection; }
      
   if (nIso) {
@@ -156 +161 @@
                                                    G4int ZZ, G4int AA, G4double /*temperature*/)
 {
+  // VI protection for Hydrogen
+  if(ZZ <= 1) { return 0.0; }
+
 //
 // Get relevant information about the projectile and target (A, Z) and