|Message: QGSP_BIC_HP HPNeutronInelastic pb with kinematics||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)|
Here is my problem. I am using geant4 version 4.9.1.p03. In my simulation,14 MeV neutrons are impinging 28Si or 16O targets. I was first interested in (n,alpha) reactions but now i am also looking at elastic as well as inelastic scattering. As i had horrible doubt about the physics list i was using, i have decided yesterday to export the PhysiscsList of hadronic/Hadr01 in my code using QGSP_BIC_HP via the command /testhadr/addPhysics QGSP_BIC_HP. Unfortunetly, the results i have are the same as the ones i had before. You will find in attachment a postcript file in which i have plotted the lab energy of the scattered 28Si versus its lab emission angle as an example. The elastic scattering is well treated and the kinetic energy ranges from 0 up to 1.8 MeV as expected. At this stage, all is ok!!!! However, as you can see on the right figure (zoom of the left one), i can distinguish curves that may result from inelastic scattering on excited states of 28Si (i guess it just because i do not see any other processes that will give me 28Si in the final state as there is no fusion evaporation in HP. right?). To obtain such a bidim, i have inserted a condition in SteppingAction.cc asking for recording events if particleName == 28Si[0.0]. I really don t understand the shape of these curves. In case of inelastic scattering ( i have not checked this but it must be) the lab emission angle has to be less than 90 deg which means that these curves must be under the elastic one. Furthemore, their shape is quite strange. Normally, if inelastic ones, for a choosen excited state, the curve must begin at 90 deg with no kinetic energy up to a maximum theta_max and then goes back to 0 Mev energy at 0 deg lab emission . I must do something wrong for sure but i am a little bit lost at this stage. I really need help. Can you tell me where am i wrong? > Thanks in advance. > jean-etienne
|Inline Depth:||Outline Depth:||Add message:|