On 12.11.2019 3:27, yangt_at_ihep.ac.cn wrote:
>
> 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
>
Dear Yang,
Correction about number of energy deposition events between Li and
alpha. Energy deposition depends upon Z and beta value, deltaE ~
Z^2/beta^2, there are more events for Li rather than for alpha.
1. In your case, it shows that charged particles with Z=+2 and Z=+3 are
stopped ( ICODE=12 no more kinetic energy). That is why the RULL value
is 0. When your score photons and electrons the RULL value will greater
than 0 and ICODE will be different.
2. Correct. One track leaves the sensitive region and a new track
begins, you can see different JTRACK values: -2, -6. I have never drawn
track trajectory in mgdraw user routine, maybe you can use USRBIN +
AUXSCORE, but it isn't post processing.
Best regards,
Mikhail
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Wed Nov 13 2019 - 06:23:55 CET