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

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

@@ -28 +28 @@
 // A prototype of the low energy neutron transport model.
 //
+// 080721 To be "ClearHistories" effective, the selection scheme of angular distribution is modified by T. Koi
+//
+//
 #include "G4NeutronHPContEnergyAngular.hh"
 
@@ -55 +58 @@
      // This is the cause of the He3 problem !!!!!!!!
      // See to it, if you can improve this.
-     G4double random = G4UniformRand();
-     G4double deltaE = theAngular[it].GetEnergy()-theAngular[it-1].GetEnergy();
-     G4double offset = theAngular[it].GetEnergy()-anEnergy;
-     if(random<offset/deltaE) it--; 
+     //080714 TK commnet Randomizing use angular distribution
+     //080714 TK Always use the upper side distribution. enabling ClearHistories method.
+     //G4double random = G4UniformRand();
+     //G4double deltaE = theAngular[it].GetEnergy()-theAngular[it-1].GetEnergy();
+     //G4double offset = theAngular[it].GetEnergy()-anEnergy;
+     //if(random<offset/deltaE) it--;
      theAngular[it].SetTarget(GetTarget());
      theAngular[it].SetTargetCode(theTargetCode);
@@ -84 +89 @@
    return result;
 }
+
+
+ 
+void G4NeutronHPContEnergyAngular::ClearHistories()
+{ 
+   if ( theAngular!= NULL )
+   { 
+      for ( G4int i = 0 ; i< nEnergy ; i++ ) 
+         theAngular[i].ClearHistories(); 
+   } 
+}