[fluka-discuss]: [How to use TRACKR common variables]

From: Kyung-Nam Lee <leeronaldo_at_naver.com>
Date: Wed, 3 Sep 2014 16:24:08 +0900

Dear, FLUKA advanced user.

 

I would like to calculate Dose Averaged LET at each voxel.

Because given fluka scoring commands cannot give the quantity, I try to use
user routine(comscw.f).

Dose Averaged LET can be defined like this formula.

cid:image002.jpg_at_01CFC77B.9AC5E8C0

(LET_d: Dose Averaged LET, LET: particle LET, D: absorbed dose by ith
particle, D_total: Total dose at voxel)

 

For calc. of Dose averaged LET, I make a modified comscw.f file. This file
is attached on this email.

I think TRACKR common variables can be used for that work. So, the contents
is like this.

========================= (The partial contents of comscw.f, grammar: not
considered exactly) =================

If (ISCRNG .EQ. 1) THEN

       If (JSCRNG .EQ. 2) THEN

 
COMSCW=(100.0*Dtrack(1)/Ttrack(1))*(1.602176462**(-7.0)*Dtrack(1)/RHO(MEDFLK
(MREG,1)))/Rull

       END If

END If

============================================================================
===================================

As you know, this code cannot be used because we cannot know D_total not
yet.

 

So, most of all, I would like to ask whether the usage of Dtrack and Ttrack
on this code is used properly or not.

I think Dtrack(i) and Ttrack(i) means energy deposition of the ith
deposition event and length of the ith track segment, respectively.

But, I input "1" temporarily as a index number because I cannot know how to
get index number "i" from commom variables.

If you give proper usage about using Dtrack and Ttrack, it will be great
help for modifying the code.

 

Also, for Dose avg. LET, Summation value of LET_i*Dose_i at each voxel needs
to be divided by D_total(dose at voxel).

But, D_total cannot know unless all interactions are finished.

So, I would like to ask how to handle this.

 

Thank you for your reading.

I'm looking forward to your kind reply.^^

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: <mailto:leeronaldo_at_naver.com> leeronaldo_at_naver.com (with big file),
<mailto:leeronaldo_at_kirams.re.kr> leeronaldo_at_kirams.re.kr

======================================================================

 









image001.jpg
(image/jpeg attachment: image001.jpg)

Received on Wed Sep 03 2014 - 10:42:13 CEST

This archive was generated by hypermail 2.3.0 : Wed Sep 03 2014 - 10:42:14 CEST