Dear Yang,
Q1. You are correct about discrete energy-deposition, but also about
events on spot, such as stopped particles. I don't thick you will get
anything else (other events in ENDRAW) for charged ions with low kinetic
energy. It is possible that electrons don't generated because their
energy lower that EMFCUT threshold.
But if you know W (average energy, required to produce one ion pair for
your medium gas) you can calculate average number of electron/ion pairs
on each part of the track.
Q2. I have never done such thing, but if you have got an information
about ion tracks (position, direction, kinetic energy) you can use any
other program to track ion. I thick SRIM in GarfieldPP can only
transport alpha particles.
Best regards,
Mikhail
On 13.11.2019 10:07, yangt_at_ihep.ac.cn wrote:
>>
>> Dear Mikhail,
>>
> Thanks for your patient interpretation.
>
> 1. In the endraw_new.f routine, I guess the following statements
> score the discrete energy-deposition ( energy deposition on the spot)
> events:
>
> IF (MREG .EQ. NTARGT) THEN
>
> IF ((JTRACK .EQ. -6) .OR. (JTRACK .EQ. -2)) THEN
> WRITE (82, *) 'END = ', ICODE, JTRACK, XSCO, YSCO, ZSCO, RULL
> END IF
> END IF
> RETURN
> If they are, why we get only the events of energy-deposition=0 (ie,
> RULL=0)? Does it mean the discrete deposition only happens at the
> kinetic energy~0? And the number of discrete deposition is rather
> smaller compared to the total deposition events.
>
> C Energy deposition on the track
> IF (MREG .EQ. NTARGT) THEN
> WRITE ( 82, *) 'NTRACK = ', NTRACK
> WRITE ( 82, *) 'MTRACK = ', MTRACK
> WRITE ( 82, *) ICODE, JTRACK, J0TRCK
> WRITE ( 82, *) (XTRACK(J), YTRACK(J), ZTRACK(J), J = 1, NTRACK),
> & (DTRACK(J), J = 1, MTRACK)
> WRITE ( 82, *)
> END IF
>
> the above code is about the continuous deposition, it doesn't specify
> the particle type, I don't know why we only get the results of Li and
> α? I think deposition events from secondary electrons(JTRACK=3) also
> exist, why are they not be recorded?
>
>
> 2. What I want to obtain is the ionized electrons spatial
> distribution, but I don't know how to finish it. Then I use
> the indirect method--obtain the energy deposition distribution in the
> sensitive region, so I score the ions' energy deposition using mgdraw,
> but in the discussion website,I find the energy deposition should be
> essentially attributed to
> electrons.(http://www.fluka.org/web_archive/earchive/new-fluka-discuss/7931.html),
> so is it correct if I only score the energy deposition of ions,
> whether I need to add the energy deposition contribution from
> electrons? Or, if there is a direct method to obtain the electrons
> spatial distribution in the sensitive region? I think I
> have some wrong understandings, please point out them.
>
>
> Best regards!
>
> Yang
>
>
>
>
> -----原始邮件-----
> *发件人:*"Mikhail Polkovnikov <Михаил Полковников>" <pmk_at_ihep.ru>
> *发送时间:*2019-11-13 12:00:18 (星期三)
> *收件人:* yangt_at_ihep.ac.cn, fluka-discuss_at_fluka.org
> *抄送:*
> *主题:* Re: [fluka-discuss]: Re: Fw: Fw: problems of mgdraw to score
> energy deposition position and value.
>
> 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 Thu Nov 14 2019 - 17:59:52 CET