Message: Trouble with G4FissLib Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Trouble with G4FissLib 

Forum: Hadronic Processes
Date: 19 Dec, 2012
From: Matthew Blackston <Matthew Blackston>

I'm using version 4.9.6 and I'm getting warnings and a fatal error related to induced 
fission when I use G4FissLib.

My goal was to do a comparison of G4NeutronHPFission and G4FissLib for induced
fission. To do the test, I modeled the physics list on QGSP_BIC_HP but created a 
copy of G4NeutronHPBuilder, and modified the following method:

void myNeutronHPBuilder::
Build(G4HadronFissionProcess * aP)
{
//  if(theHPFission == 0) theHPFission = new G4NeutronHPFission;
//  theHPFission->SetMinEnergy(theMin);
//  theHPFission->SetMaxEnergy(theMax);
//  if(theHPFissionData==0) theHPFissionData=new G4NeutronHPFissionData;
//  aP->AddDataSet(theHPFissionData);
//  aP->RegisterMe(theHPFission);

    //LLNL model
    if(theHPFission == 0) myHPFission = new G4FissLib;
    myHPFission->SetMinEnergy(theMin);
    myHPFission->SetMaxEnergy(theMax);
    if(theHPFissionData==0) theHPFissionData=new G4NeutronHPFissionData;
    aP->AddDataSet(theHPFissionData);
    aP->RegisterMe(myHPFission);
}
 
When I do this, I get warnings like: 


Warning: try to access G4NeutronHPReactionWhiteBoard before opening

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : had012
      issued by : G4HadronicProcess:CheckResult()
Warning: Bad energy non-conservation detected, will re-sample the interaction
 Process / Model: nFission / HadronicModel
 Primary: neutron (2112), E= 948.439, target nucleus (92,235)
 E(initial - final) = 217010 MeV.

*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------


And then finally a fatal exception:


-------- EEEE ------- G4Exception-START -------- EEEE -------
*** G4Exception : had006
      issued by : G4HadronicProcess::PostStepDoIt
Call for HadronicModel
Target element Uranium  Z= 92  A= 235
Unrecoverable error in the method ApplyYourself of nFission
TrackID= 1  ParentID= 0  neutron
Ekin(GeV)= 0.00887385;  direction= (-1.62156e-05,-1.39439e-06,1)
Position(mm)= (-0.0146393,-0.00125885,902.795);  material U
PhysicalVolume  <test_physV>
 ApplyYourself does not completed after 100 attempts

*** Fatal Exception *** core dump ***
-------- EEEE -------- G4Exception-END --------- EEEE -------



G4NeutronHPReactionWhiteBoard and G4NeutronHPManager have only been 
added to version 9.6. There are several changes to the G4NeutronHPFission 
class related to this, but I don't see similar changes in the G4FissLib class. I'm 
wondering if the warnings and errors are related to other ways in which the 
G4FissLib class have not been updated with newer releases.

Am I introducing G4FissLib incorrectly or is there something that needs to
be fixed?

Thanks,
Matthew

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

1 None: Re: Trouble with G4FissLib   (Matthew Blackston - 19 Dec, 2012)
 Add Message Add Message
to: "Trouble with G4FissLib"

 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 ]