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

From: Kyung-Nam Lee <>
Date: Tue, 2 Sep 2014 11:35:27 +0900

Dear, Dr. Vittorio Boccone.


Most of all, I really appreciate your kind and quick reply.


At the first question, I can know exactly due to your advice about scoring
time. It means scoring time is each deposition step.


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.

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.

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 ?


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 ?


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)


#215-4 Gongneung-dong Nowon-gu

Seoul 139-709, Korea

Tel: 82-2-970-1737

Fax: 82-2-970-1962


E-Mail: (with big file),



From: Dr. Vittorio Boccone []
Sent: Tuesday, September 02, 2014 5:32 AM
To:; >>
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
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
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 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
Received on Tue Sep 02 2014 - 05:54:57 CEST

This archive was generated by hypermail 2.3.0 : Tue Sep 02 2014 - 05:54:59 CEST