Fwd: RE: [fluka-discuss]: How to calculate LET(linear energy transfer) with voxel-by-voxel basis

From: Vittorio Boccone <dr.vittorio.boccone_at_ieee.org>
Date: Wed, 3 Sep 2014 21:31:50 +0200

---------- Forwarded message ----------
From: Vittorio Boccone <dr.vittorio.boccone_at_ieee.org>
Date: Wed, Sep 3, 2014 at 12:39 AM
Subject: Re: RE: [fluka-discuss]: How to calculate LET(linear energy
transfer) with voxel-by-voxel basis
To: Kyung-Nam Lee <leeronaldo_at_naver.com>


Dear Kyung-Nam Lee,



On Tue, Sep 2, 2014 at 4:35 AM, Kyung-Nam Lee <leeronaldo_at_naver.com> wrote:
>
> At the second question, The reason why I would like to have the deq99c.f
> file is not for calculation of effective dose.
>
> What I really want to know through deq99c.f file is proper usage of
> "COMMOM variables" about calc. of particle LET and its usage for equivalent
> dose or effective dose.
>
This external function is really deprecated and it's not used since very
long time. This mean that trying to use it might give you other problems.
This is why it is not given anymore with the code.


> I mean, for calc. of equivalent dose, as you know, we need to have
> radiation weighting factor.
>
> This weighting factor can vary for heavy ions due to varied LETs at each
> voxel(3-D bin with Cartesian coordinate).
>
> Therefore, I think the user routine file is more likely to include the
> whole process from calc. of particle LET to loading or calc. of the
> weighting factor at each voxel.
>
I really don't see the point on why you want to do it manually. Voxel
structure is only used for collecting up the scoring and defining the
"regions" of where the material is. The steps on which the interactions and
the scoring values are calculate are much smaller than the Voxel
segmentation. I suggest you to read the introductory part of the manual to
start to understand how FLUKA works.

In addition you can profit from the pdf of the FLUKA schools if you go on
the course page.

> So, if I can learn this process with useful "COMMOM variables"(e.g.
> JTRACK, XTRACK, YTRACK, ZTRACK, MTRACK. et al.), it will be great help and
> practice for me to apply to the other case.
>
> So, although scoring with user routine is rather risky way, is there good
> way for download this file ?
>

I strongly discourage you to use this old routine. It's not even good for
education purposes.
There are other routines where you can practice the usage of the COMMON
variables, look at the advanced course lectures, there's also an exercise
(with solution) that might help you on the way.


> At your comment for scoring LET by using USRYIELD,
>
> As you know, USRYIELD can give many boundary crossing quantity like USRBDX
> command. But, I would like to scoring LET at each voxel.
>
> So, I think USRYIELD is not more likely to be used for this purpose. I
> think this work can be performed by using USRBIN command and user
> routine(comscw.f or fluscw.f).
>
> Is it correct ?
>

If you want to build a custom LET map I'm afraid you need to do it by hand
http://www.fluka.org/web_archive/earchive/new-fluka-discuss/4807.html

Thank you for your reading.
>
> I'm looking forward to your advice.^^
>
> Have a nice day!!~~^^
>
>
>
> Best Regards,
>
> Kyung-Nam Lee
>
>
>
> ======================================================================
>
> Kyung-Nam Lee, M.S.
>
> Research Physicist
>
> Medical Physics Research Team
>
> Division of Heavy Ion Clinical Researches
>
> Korea Heavy Ion Medical Accelerator Project (KHIMA Project)
>
> KOREA INSTITUTE of RADIOLOGICAL & MEDICAL SCIENCES (KIRAMS)
>
> #215-4 Gongneung-dong Nowon-gu
>
> Seoul 139-709, Korea
>
> Tel: 82-2-970-1737
>
> Fax: 82-2-970-1962
>
> Mobile:82-10-6366-9161
>
> E-Mail: leeronaldo_at_naver.com (with big file), leeronaldo_at_kirams.re.kr
>
> ======================================================================
>
>
>
> *From:* Dr. Vittorio Boccone [mailto:dr.vittorio.boccone_at_ieee.org]
> *Sent:* Tuesday, September 02, 2014 5:32 AM
> *To:* leeronaldo_at_naver.com; fluka-discuss_at_fluka.org >> "
> fluka-discuss_at_fluka.org"
> *Subject:* Re: [fluka-discuss]: How to calculate LET(linear energy
> transfer) with voxel-by-voxel basis
>
>
>
> Dear Kyung-Nam Lee,
> I assume you already went through the manual Section 13.2.2
>
> First of all, I hope to ask when the routine(comscw.f) is called
>
> specifically.
>
> I mean scoring time, namely, whether the routine is called after end of
>
> "event" or end of "each interaction" or not.
>
> Is it depend on used scoring command(e.g. USRBIN, USRYIELD, USRBDX) ??
>
>
> You enable the comsw routine with from the WHAT(6) of USERWEIG.
> It's as a matter of fact a weighting factor which is calculated at scoring
> time per each deposition step and modifies directly the value scored with
> USRBIN etc...
>
>
> Also, I would like to ask how can I get the deq99.f(user routine for calc.
>
> of effective dose) file.
>
> I could not download the file at this
>
> link(http://info-fluka-discussion.web.cern.ch/info-fluka-discussion/download
>
> /deq99c.f)
>
> So, could you attach this file when you reply ?
>
> Effective dose is now calculated by default using the AUXSCORE card. We
> will correct the manual to removed the part mentioning of the usage of the
> deq99.f routine which is at the moment strongly deprecated.
>
>
> Because I am a beginner for modifying user routines, It will be great help
>
> for me if I understand the contents of the file.
>
> Then, I can know how to calculate LET by using useful "COMMON variable". In addition, I would like to actively listen to any kind of your advice about calculating Dose Averaged LET(my final goal in this process) with
>
> FLUKA user routine.
>
> I think you can just score the LET with USRYIELD and you don't really need
> to recalculate it (unless you need to do it on purpose).
> Please check
> http://www.fluka.org/web_archive/earchive/new-fluka-discuss/6829.html and
> the other LET related question on the forum, they might give you more
> answers and clarify a bit the LET scoring option you have.
>
> I am looking forward to your generous reply.
>
> It might not have been a very generous reply... but I hope it clarifies a
> bit the scoring of the LETs
> Vittorio
>
Received on Wed Sep 03 2014 - 23:07:04 CEST

This archive was generated by hypermail 2.3.0 : Wed Sep 03 2014 - 23:07:05 CEST