From: <yangt_at_ihep.ac.cn>
Date: Tue, 12 Nov 2019 08:27:59 +0800 (GMT+08:00)
Dear Mikhail,
Thanks for your reply, I am very grateful to you. I read your modified subroutine file, and the source.f file remains unchanged, is it? For the endraw_new.f file, you add code at the front of any entry to sorce the continuous energy deposition, so does my former endraw.f file only scores the discrete energy deposition and that leads to a unrealistic small energy-deposition number? I run your modified routine and it works, and the event number is reasonable. I look the output, and I have some questions needed to confirm.
1. I list the point energy-deposition (attachment:"end.png"), and I find the last column corresponding to energy deposition is 0, which is unreasonable, and I see the corresponding statement is "WRITE (82, *) 'END = ', ICODE, JTRACK, XSCO, YSCO, ZSCO, RULL", so the variable "RULL" should be changed? and how to do?
2. The continuous energy-deposition shows that the z coordinate increases to some value and then restarts at z~0 (attachment:"track.png"), does it mean a ion trajectory is end or leave the sensitive region and a new track is open? indeed, I always don't know how to get the particle trajectory in FLUKA, not like as SRIM or Geant4, the particle track can be seen. so if I need to plot the particle trajectory in postprocessing, I can use mgdraw routine like this? or is there a build-in card or method in fluka?
Best regards!
Yang
-----原始邮件-----
发件人:"Mikhail Polkovnikov <Михаил Полковников>" <pmk_at_ihep.ru>
发送时间:2019-11-10 02:55:00 (星期日)
收件人: yangt_at_ihep.ac.cn, fluka-discuss_at_fluka.org
抄送:
主题: [fluka-discuss]: Re: Fw: Fw: problems of mgdraw to score energy deposition position and value.
On 04.11.2019 9:57, yangt_at_ihep.ac.cn wrote:
Dear Mikhail,
I am sorry to bother you, however, what i querry is an old problem listed on the discussion website：http://www.fluka.org/web_archive/earchive/new-fluka-discuss/15973.html
I submited the problems several days ago, but I don't receive any reply and I cannot solve it yet. (attachmets: detector2.inp source.f endraw.f...)
Recently I have made some simulations using mgdraw.f subroutine to score the energy deposition.But there exists some unsolved problems/errors.
What I want to score is the continuous energy deposition and discrete/point energy deposition, I don't know if both the two types of deposition can be accessed only with mgdraw.f. in my simulation, "source.f" samples the primary from file "detector001_ions.dat" which has 169 charged particles(71#α+ 98# Li+).
1) But I only obtain 97 energy-deposition events using "endraw.f" ("detector2001_endraw.dat", I find I define a wrong momentum in source.f in previous version submitted on the website and I correct it in the attachment),as we known, a charged particle will undergoes many energy-deposition or ionization events along its trajectory, so the total number of energy-deposition events should be much greater than the incident ions' number, but I only get 97 events which is even significantly lower than the ion's number (which is 169 seen in detector001_ions.dat).
2) And another strange thing is it obtains much more Li+ deposition events(jtrack=-2) than α seen in file "energy depositon in dat", it is also strange, what's wrong?
If you have some suggestions, please tell me, and publish the answer on the discussion website,I think many people have the same problem.
any suggestion is appreciable!
Best regards!
Yang.
Dear Yang,
1) Continuous energy loss goes along the track (SUBROUTINE MGDRAW), energy deposition on the spot (ENTRY ENDRAW). You will a lot of energy deposition events on the track, see attachment.
2) Lithium is heavier than alpha particle, and its energy is less.
I have changes some elements in attachment.
Best regards,
Mikhail
__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
Received on Tue Nov 12 2019 - 02:51:24 CET
