Message: Problems with neutron physics list in an energy range Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Problems with neutron physics list in an energy range 

Forum: Hadronic Processes
Date: 20 Dec, 2010
From: Manuel Castro Avila <Manuel Castro Avila>

Hi everyone,

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???

Thanks,

Manuel

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

1 Idea: Re: Problems with neutron physics list in an energy range   (Vladimir Ivanchenko - 04 Jan, 2011)
 Add Message Add Message
to: "Problems with neutron physics list in an energy range"

 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 ]