|Message: Changes in inelastic scattering.||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 simulated neutron propagation. Detector is cube 200x200x200 meters filled with oxygen.
I used physics list for neutrons: G4PreCompoundModel* thePreEquilib2; thePreEquilib2 = new G4PreCompoundModel(theHandler); thePreEquilib2->SetMinEnergy(20*MeV); thePreEquilib2->SetMaxEnergy(70*MeV);
G4NeutronInelasticProcess* ninelProc = new G4NeutronInelasticProcess(); G4NeutronInelasticCrossSection* neutron_XC = new G4NeutronInelasticCrossSection(); ninelProc->AddDataSet(neutron_XC);
G4NeutronHPInelastic* theNeutronHPInelastic = new G4NeutronHPInelastic(); theNeutronHPInelastic->SetMaxEnergy(20.*MeV ); G4NeutronHPInelasticData* theHPInelasticData = new G4NeutronHPInelasticData();
G4LENeutronInelastic* LEPnModel = new G4LENeutronInelastic(); G4double LEPUpperLimit = 25*GeV; G4double LEPpnpiLimit = 9.5*GeV; LEPnModel->SetMinEnergy(LEPpnpiLimit); LEPnModel->SetMaxEnergy(LEPUpperLimit);
ninelProc->RegisterMe(theNeutronHPInelastic); ninelProc->AddDataSet(theHPInelasticData); ninelProc->RegisterMe(thePreEquilib2); pmanager->AddDiscreteProcess(ninelProc);
I was looking for energy loss in steps, defined by Inelastic scattering. When starting neutron energy is 19 MeV, i have discrete spectrum of energy loss by Inelastic, and it seems to be equal for both version of Geant4 (4.8, 4.91). But when I run neutrons with energy of 21 MeV spectrum differs. Geant4.8 has a great peak in energy < 4 MeV. It seems very strange for me, because I didn't find anything about such changes.
regards, Maxim Gonchar
|Inline Depth:||Outline Depth:||Add message:|