Message: possible bug fix in NeutronHPFission model Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None possible bug fix in NeutronHPFission model 

Forum: Hadronic Processes
Date: 19 Nov, 2007
From: Victor_Makarov <Victor_Makarov>

As I've already posted in Bugzilla, NeutronHPFission model produces too many fission neutrons with energies around 20 MeV. After some investigations I found out that error hides in G4NeutronHPVector which returns energy value to G4NeutronHPArbitaryTab. I'm not very good in Geant 4 programming, so I could not locate the error itself in G4NeutronHPVector. But I tried this: I changed next lines of G4NeutronHPArbitaryTab.cc

theBuffer.Merge(theManager.GetScheme(low), anEnergy, 
               theDistFunc+low, theDistFunc+high);
return theBuffer.Sample();

to these:

theBuffer.Merge(theManager.GetScheme(low), anEnergy, 
               theDistFunc+low, theDistFunc+high);
G4double ee;
ee = theBuffer.Sample();
while (ee>19.0*MeV)  ee = theBuffer.Sample();
return ee;

The result was quite good - energy distribution became identical to LLNL fission model's and to Geant 3's. Of course, I'll be very thankful if one can correct G4NeutronHPVector.cc itself and share the solution. Comments are welcome!

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: possible bug fix in NeutronHPFission model   (Koi, Tatsumi - 06 Dec, 2007)
 Add Message Add Message
to: "possible bug fix in NeutronHPFission model"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]