<div dir="ltr"><div dir="ltr"><br></div><div class="gmail_quote"><div dir="ltr" class="gmail_attr"><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><br clear="all"><div><div dir="ltr" class="gmail-m_-6309180306086249390gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><span style="color:rgb(51,51,51);font-family:arial,helvetica;font-size:medium">I am trying to create 14 MeV neutrons from deuterium beam on deuterium target. For this, I used the physics list QGSP_BIC_AllHP. with dueteron energy of 11 MeV</span> . </div><div dir="ltr"><br></div><div><span style="color:rgb(51,51,51);font-family:arial,helvetica;font-size:medium">I expect to see a 14.0 MeV neutrons because of the reaction's positive Q-value 3.27 MeV. All I see now is breakup neutrons without the 14MeV peak</span><br></div><div><span style="color:rgb(51,51,51);font-family:arial,helvetica;font-size:medium"><br></span></div><div><div><img src="cid:ii_jx6jr5nd0" alt="image.png" width="542" height="269"><br></div></div><div><span style="color:rgb(51,51,51);font-family:arial,helvetica;font-size:medium"><br></span></div><div><span style="color:rgb(51,51,51);font-family:arial,helvetica;font-size:medium">Please what am i doing wrong? can some help ?? thanks </span></div><div dir="ltr"><br></div><div>Physics list </div><div><br></div><div><br></div><div>/gate/physics/addPhysicsList QGSP_BERT_HP <br>/gate/physics/addPhysicsList QGSP_BIC_AllHP </div><div>/gate/physics/addPhysicsList FTFP_BERT_HP <br></div><div><br></div><div>## Neutron Capture<br>/gate/physics/addProcess NeutronCapture<br>/gate/physics/processes/NeutronCapture/setModel G4LCapture<br>#/gate/physics/processes/NeutronCapture/setModel G4NeutronHPorLCapture<br>/gate/physics/processes/NeutronCapture/setModel G4NeutronHPCapture<br>/gate/physics/processes/NeutronCapture/setDataSet G4NeutronHPCaptureData<br><br># inelastic<br>/gate/physics/addProcess NeutronInelastic<br>/gate/physics/processes/NeutronInelastic/setModel PreCompound<br>#/gate/physics/processes/NeutronInelastic/setModel G4NutronInelastic<br>/gate/physics/processes/NeutronInelastic/setDataSet G4NeutronHPInelasticData<br><br># Fission<br>/gate/physics/addProcess Fission<br>/gate/physics/processes/Fission/setModel G4LFission<br>/gate/physics/processes/Fission/setDataSet G4NeutronHPFissionData<br>#proton.deuteron,triton and alpha<br>/gate/physics/addProcess IonInelastic Default<br>/gate/physics/processes/IonInelastic/setModel G4BinaryLightIonReaction Default<br>/gate/physics/processes/IonInelastic/setModel G4LEDeuteronInelastic deuteron<br>/gate/physics/processes/IonInelastic/setModel G4LETritonInelastic triton<br>/gate/physics/processes/IonInelastic/setModel G4LEAlphaInelastic alpha<br>/gate/physics/processes/IonInelastic/G4BinaryLightIonReaction/setEmin 0 MeV deuteron<br>/gate/physics/processes/IonInelastic/G4BinaryLightIonReaction/setEmax 20 GeV deuteron<br>/gate/physics/processes/IonInelastic/G4BinaryLightIonReaction/setEmin 0 MeV triton<br>/gate/physics/processes/IonInelastic/G4BinaryLightIonReaction/setEmax 20 GeV triton<br>/gate/physics/processes/IonInelastic/G4BinaryLightIonReaction/setEmin 0 MeV alpha<br>/gate/physics/processes/IonInelastic/G4BinaryLightIonReaction/setEmax 20 GeV alpha<br>/gate/physics/processes/IonInelastic/G4LEDeuteronInelastic/setEmin 0 MeV deuteron<br>/gate/physics/processes/IonInelastic/G4LEDeuteronInelastic/setEmax 80 MeV deuteron<br>/gate/physics/processes/IonInelastic/G4LETritonInelastic/setEmin 0 MeV triton<br>/gate/physics/processes/IonInelastic/G4LETritonInelastic/setEmax 80 MeV triton<br>/gate/physics/processes/IonInelastic/G4LEAlphaInelastic/setEmin 0 MeV alpha<br>/gate/physics/processes/IonInelastic/G4LEAlphaInelastic/setEmax 80 MeV alpha<br>/gate/physics/processes/IonInelastic/setDataSet G4IonsShenCrossSection GenericIon<br>/gate/physics/processes/IonInelastic/setDataSet G4TripathiLightCrossSection deuteron<br>/gate/physics/processes/IonInelastic/setDataSet G4TripathiLightCrossSection triton<br>/gate/physics/processes/IonInelastic/setDataSet G4TripathiLightCrossSection alpha<br></div><div><br></div><div><br></div><div dir="ltr"><div><div>Sincerely, </div></div><div>Inno<br></div></div></div></div></div></div></div></div>
</blockquote></div></div>