|Message: neutron_HP_Physics and NIST_Material||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'm doing calculations with low energy neutrons. I started with example /runandevent/Re02 and changed a few things. All was running fine so far. Now I wanted to include physics which is better adapted to my problem. So I changed the neutron physics into:
// neutron Physics pManager = G4Neutron::Neutron()->GetProcessManager();
// add process G4HadronElasticProcess* thenElasticProcess = new G4HadronElasticProcess(); thenElasticProcess->AddDataSet(new G4NeutronHPElasticData);
G4LElastic* thenElasticModel = new G4LElastic(); thenElasticModel->SetMinEnergy(20.*MeV); G4NeutronHPElastic* theNeutronHPElasticModel= new G4NeutronHPElastic(); theNeutronHPElasticModel ->SetMaxEnergy(20.*MeV); thenElasticProcess->RegisterMe(thenElasticModel); thenElasticProcess->RegisterMe(theNeutronHPElasticModel);
G4NeutronInelasticProcess* theNeutronInelasticProcess = new G4NeutronInelasticProcess(); theNeutronInelasticProcess->AddDataSet(new G4NeutronHPInelasticData);
G4NeutronHPInelastic* theNeutronHPInelasticModel= new G4NeutronHPInelastic(); theNeutronHPInelasticModel ->SetMaxEnergy(20.*MeV); G4BinaryCascade* theNeutronBCModel = new G4BinaryCascade(); theNeutronBCModel->SetMinEnergy(20.*MeV); theNeutronBCModel->SetMaxEnergy(6.*GeV); G4LENeutronInelastic* theNeutronLEPModel = new G4LENeutronInelastic(); theNeutronLEPModel->SetMinEnergy(4.*GeV); G4HENeutronInelastic* theNeutronHEPModel = new G4HENeutronInelastic(); theNeutronInelasticProcess->RegisterMe(theNeutronHPInelasticModel); theNeutronInelasticProcess->RegisterMe(theNeutronBCModel); theNeutronInelasticProcess->RegisterMe(theNeutronLEPModel); theNeutronInelasticProcess->RegisterMe(theNeutronHEPModel); pManager->AddDiscreteProcess(theNeutronInelasticProcess);
G4HadronFissionProcess* thenFission = new G4HadronFissionProcess(); thenFission->AddDataSet(new G4NeutronHPFissionData);
G4LFission* thenFissionModel = new G4LFission(); thenFissionModel->SetMinEnergy(20.*MeV); G4NeutronHPFission* theNeutronHPFissionModel= new G4NeutronHPFission(); theNeutronHPFissionModel ->SetMaxEnergy(20.*MeV); thenFission->RegisterMe(thenFissionModel); thenFission->RegisterMe(theNeutronHPFissionModel); pManager->AddDiscreteProcess(thenFission);
G4HadronCaptureProcess* thenCapture = new G4HadronCaptureProcess(); thenCapture->AddDataSet(new G4NeutronHPCaptureData);
G4LCapture* thenCaptureModel = new G4LCapture(); thenCaptureModel->SetMinEnergy(20.*MeV); G4NeutronHPCapture* theNeutronHPCaptureModel= new G4NeutronHPCapture(); theNeutronHPCaptureModel ->SetMaxEnergy(20.*MeV); thenCapture->RegisterMe(thenCaptureModel); thenCapture->RegisterMe(theNeutronHPCaptureModel); pManager->AddDiscreteProcess(thenCapture);
(Basically I added the HP-Physics). Now I have the Problem, that as long as I use G4_Water as material all looks good. But when I change the material to G4_TISSUE_SOFT_ICRP (or some other ICRP-material), the tracks start looking very strange, the neutrons get reflected at the surface and there are a lot of electrons produced.
My question is: Did I make a mistake in the neutron-physics list? Or do I have to add something for the NIST-materials? In the end I wanted to use several materials out of the NIST-list, so it would be very convenient if there is a way to use it!
Thanks a lot, Sylvia
|Inline Depth:||Outline Depth:||Add message:|