|Message: Problems with neutron physics list in an energy range||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
I am using the following neutron physic list into my simulation:
// Inelastic G4NeutronInelasticProcess* theInelasticProcess = new 4NeutronInelasticProcess(); G4LENeutronInelastic* model1 = new G4LENeutronInelastic(); theInelasticProcess->RegisterMe(model1); G4NeutronHPInelasticData* theNeutronData1 = new G4NeutronHPInelasticData; theInelasticProcess->AddDataSet(theNeutronData1); pmanager->AddDiscreteProcess(theInelasticProcess);
// Capture G4HadronCaptureProcess* theCaptureProcess = new G4HadronCaptureProcess; G4LCapture* captureModel = new G4LCapture(); theCaptureProcess->RegisterMe(captureModel); G4NeutronHPCaptureData* theNeutronData3 = new G4NeutronHPCaptureData; theCaptureProcess->AddDataSet(theNeutronData3); pmanager->AddDiscreteProcess(theCaptureProcess);
// elastic scattering G4HadronElasticProcess* theNeutronElasticProcess = new G4HadronElasticProcess; G4LElastic* theElasticModel1 = new G4LElastic; theNeutronElasticProcess->RegisterMe(theElasticModel1); G4NeutronHPElasticData * theNeutronData = new G4NeutronHPElasticData; theNeutronElasticProcess->AddDataSet(theNeutronData); pmanager->AddDiscreteProcess(theNeutronElasticProcess);
//Fission G4HadronFissionProcess* fissionProc = new G4HadronFissionProcess(); G4LFission* fissionModel = new G4LFission(); fissionProc->RegisterMe(fissionModel); pmanager->AddDiscreteProcess(fissionProc);
So, I ran some simulations using neutrons with energy between 0.8-100 Mev and all worked without problems. Now I try to use neutrons between 100-4000 MeV and when I run the simulation, appears an error:
### Run 1 start. Start Run processing. 0 Segmentation fault
if I comment all the Inelastic physics, the simulation run without problems. I checked about the Inelastic model G4LENeutronInelastic, and it works in the energy range which I use the neutrons (100-4000 MeV).
I would like your help, because I do not know what is wrong???
|Inline Depth:||Outline Depth:||Add message:|