Changeset 962 for trunk/source/processes/hadronic/models/neutron_hp/src/G4NeutronHPProduct.cc

Timestamp:

Apr 6, 2009, 12:30:29 PM (15 years ago)

Author:

garnier

Message:

update processes

File:

• trunk/source/processes/hadronic/models/neutron_hp/src/G4NeutronHPProduct.cc (modified) (5 diffs)

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

@@ -27 +27 @@
 // J.P. Wellisch, Nov-1996
 // A prototype of the low energy neutron transport model.
+//
+// 080718 As for secondary photons, if its mean value has a value of integer, 
+//        then a sampling of multiplicity that based on Poisson Distribution 
+//        is not carried out and the mean is used as a multiplicity.
+//        modified by T. Koi.
+// 080721 Using ClearHistories() methodl for limiting the sum of secondary energies
+//        modified by T. Koi.
+// 080901 bug fix of too many secnodaries production in nd reactinos by T. Koi
+//
 #include "G4NeutronHPProduct.hh" 
 #include "G4Poisson.hh"
@@ -38 +47 @@
   G4int multi;
   multi = G4int(mean+0.0001);
-  if(theMassCode==0) multi = G4Poisson(mean); // @@@@gammas. please X-check this
+  //if(theMassCode==0) multi = G4Poisson(mean); // @@@@gammas. please X-check this
+  //080718
+  if ( theMassCode == 0 )
+  { 
+     if ( G4int ( mean ) == mean )
+     {
+        multi = (G4int) mean;
+     }
+     else
+     {
+        multi = G4Poisson ( mean ); 
+     }
+  }
   theDist->SetTarget(theTarget);
   theDist->SetNeutron(theNeutron);
@@ -46 +67 @@
   theCurrentMultiplicity = static_cast<G4int>(mean);
   G4ReactionProduct * tmp;
+  theDist->ClearHistories();
   for(i=0;i<multi;i++)
   {
@@ -56 +78 @@
     delete  tmp;
   }
+/*
+//080901 TK Comment out, too many secondaries are produced in deuteron reactions
   if(theTarget->GetMass()<2*GeV) // @@@ take care of residuals in all cases
   {
@@ -62 +86 @@
     if(tmp != 0) { result->push_back(tmp); }
   }
+*/
   return result;
 }