|Message: Re: 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)|
On Mon, 8 May 2006, Sylvia Studeny wrote: > *** Discussion title: Hadronic Processes > Email replies to PublicHyperNews@slac.stanford.edu must include: > In-Reply-To: <"/hadronprocess/491"@geant4-hn.slac.stanford.edu> > Subject: ...change this to be about your reply. > > Hi, > > 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); > > pManager->AddDiscreteProcess(thenElasticProcess); > > 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 > Hello, Material definition via NIST does not do something magic - you can dump the material list and see what is defined - it is equvivalent what you need. So, if you have problems, they are in physics. I do not see any obvious problem in your physics list. Probably, you need to study few events and identify, which process provides artifitial reflections. VI
|Inline Depth:||Outline Depth:||Add message:|