Changeset 1347 for trunk/source/processes/hadronic/models/neutron_hp/src/G4NeutronHPInelasticCompFS.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/G4NeutronHPInelasticCompFS.cc (modified) (2 diffs)

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

@@ -41 +41 @@
 //        add two_body_reaction
 // 100909 add safty 
+// 101111 add safty for _nat_ data case in Binary reaction, but break conservation  
 //
 #include "G4NeutronHPInelasticCompFS.hh"
@@ -681 +682 @@
    G4double AA = hadron->GetDefinition()->GetPDGMass() / proj->GetDefinition()->GetPDGMass(); 
    G4double E1 = proj->GetKineticEnergy();
-   G4double beta = std::sqrt ( A*(A+1-AA)/AA*(1+(1+A)/A*Q/E1) );
+
+// 101111
+// In _nat_ data (Q+E1) could become negative value, following line is safty for this case.
+   //if ( (Q+E1) < 0 ) 
+   if ( ( 1 + (1+A)/A*Q/E1 ) < 0 ) 
+   {
+// 1.0e-6 eV is additional safty for numeric precision 
+      Q = -( A/(1+A)*E1 ) + 1.0e-6*eV;
+   }
+
+   G4double beta = std::sqrt ( A*(A+1-AA)/AA*( 1 + (1+A)/A*Q/E1 ) );
    G4double gamma = AA/(A+1-AA)*beta;
    G4double E3 = AA/std::pow((1+A),2)*(beta*beta+1+2*beta*mu)*E1;