|Message: Re: confused on the results about neutron yield and spectra produced from reaction (d,n)||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)|
Dear Vladimir Ivantchenko,
Thank for your quick respond very much!
>>You may want to try out $G4INSTALL/examples/extended/hadronic/Hadr01 >>where different PhysLists are available. You do not need write your code but >>try to see if the reaction of interest exist.
I have tried out example Hadr01. Here is macro file i used:
/gun/energy 200 MeV /gun/particle deuteron
/control/verbose 2 /run/verbose 2
/vis/viewer/set/viewpointThetaPhi 90 0 deg /vis/viewer/zoom 1.0 /vis/viewer/flush
/tracking/storeTrajectory 1 /vis/scene/endOfEventAction accumulate
/testhadr/TargetMat G4_Be /testhadr/TargetRadius 1 cm /testhadr/TargetLength 10 cm /testhadr/NumberDivZ 100 /testhadr/DrawTracks all /testhadr/CutsAll 1 mm /testhadr/ListPhysics QGSP_BIC
# /testhadr/HistoName test /testhadr/HistoType aida #/testhadr/HistoOption "uncompress"
/run/initialize /run/beamOn 10000
However, I could not see any neutron. Here are results:
Run terminated. Run Summary Number of events processed : 10000 User=63.547s Real=99.515s Sys=24.953s RunAction: End of run actions are started /vis/viewer/update HistoManager: End of run actions are started ======================================================== Beam particle deuteron Beam Energy(MeV) 200 Number of events 10000 Average energy deposit (MeV) 177.2 RMS(MeV) 4.95 Average number of steps 100 Average number of gamma 0 Average number of e- 0 Average number of e+ 0 Average number of neutrons 0 Average number of protons 0 Average number of antiprotons 0 Average number of pi+ & pi- 0 Average number of pi0 0 Average number of kaons 0 Average number of muons 0 Average number of deuterons+tritons 0 Average number of He3+alpha 0 Average number of ions 0 Average number of forward neutrons 0 Average number of reflected neutrons 0 Average number of leaked neutrons 0 Average number of proton leak 0 Average number of pion leak 0 ========================================================
Closing the tree... Histograms and Ntuples are saved WARNING: G4VisManager::EndOfRun: Event keeping was suspended. The number of events in the run exceeded the maximum to be kept, 100. The number of events to be kept can be changed with "/vis/scene/endOfEventAction accumulate <N>", where N is the maximum number you wish to allow. N < 0 means "unlimited". 100 events have been kept for refreshing and/or reviewing.
On the other hand, Although there is the report:
### Run 1 start ### Histo books 22 histogramsHisto: tree store : test.aida
but I could not see the histogram file test.aida.
Could you please let me know what i should do to solve this problem.
Thank you very much in advance!