Changeset 1347 for trunk/source/processes/hadronic/models/neutron_hp/src/G4NeutronHPInelasticBaseFS.cc

Timestamp:

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

Author:

garnier

Message:

geant4 tag 9.4

File:

• trunk/source/processes/hadronic/models/neutron_hp/src/G4NeutronHPInelasticBaseFS.cc (modified) (4 diffs)

trunk/source/processes/hadronic/models/neutron_hp/src/G4NeutronHPInelasticBaseFS.cc

@@ -31 +31 @@
 //        Introduce theNDLDataA,Z which has A and Z of NDL data by T. Koi
 // 081024 G4NucleiPropertiesTable:: to G4NucleiProperties::
+// 101111 Add Special treatment for Be9(n,2n)Be8(2a) case by T. Koi
 //
 #include "G4NeutronHPInelasticBaseFS.hh"
@@ -360 +361 @@
   else if(theEnergyAngData!=0)
   {
+
     G4double theGammaEnergy = theEnergyAngData->GetTotalMeanEnergy();
     G4double anEnergy = boosted.GetKineticEnergy();
@@ -371 +373 @@
     G4double eBindHe3 = G4NucleiProperties::GetBindingEnergy(3,2);
     G4double eBindA = G4NucleiProperties::GetBindingEnergy(4,2);
+    G4int ia=0;
     for(i=0; i<tmpHadrons->size(); i++)
     {
@@ -396 +399 @@
       {
         eBindProducts+=eBindA;
-      }
-    }
+        ia++; 
+      }
+    }
+
     theGammaEnergy += eBindProducts;
+
+//101111 
+//Special treatment for Be9 + n -> 2n + Be8 -> 2n + a + a
+if ( (G4int)(theBaseZ+eps) == 4 && (G4int)(theBaseA+eps) == 9 )
+{
+   // This only valid for G4NDL3.13,,,
+   if ( std::abs( theNuclearMassDifference -   
+        ( G4NucleiProperties::GetBindingEnergy( 8 , 4 ) -
+        G4NucleiProperties::GetBindingEnergy( 9 , 4 ) ) ) < 1*keV 
+      && ia == 2 )
+   {
+      theGammaEnergy -= (2*eBindA);
+   }
+}
     
     G4ReactionProductVector * theOtherPhotons = 0;