- Contemporary messages sorted: [ by date ] [ by thread ] [ by subject ] [ by author ] [ by messages with attachments ]

From: Kyung-Nam Lee <leeronaldo_at_naver.com>

Date: Mon, 15 Sep 2014 23:44:28 +0900

Dear, Dr. Vittorio Boccone.

Most of all, I would like to say thank you for your specific comments.^^

It is great help for me to understand hooking action through fluka.

Based on your advices, I think there is no hooking action after end of each

"event" and each "cycle" on fluka.

I mean there is only hooking action after end of each interaction.

Ok!. as you mentioned, it is likely that off-line calculation is only way

for calc. of dose avg. LET.

I will calculate dose avg. LET by using possible ways you mentioned already.

Thank you for your help!!~~^^

Have a nice day!!~~^^

Best Regards,

Kyung-Nam Lee

From: owner-fluka-discuss_at_mi.infn.it [mailto:owner-fluka-discuss_at_mi.infn.it]

On Behalf Of Dr. Vittorio Boccone

Sent: Monday, September 15, 2014 6:47 AM

To: Kyung-Nam Lee; fluka-discuss

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

On 9/3/14 9:24 AM, Kyung-Nam Lee wrote:

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.

TRACKR common variables contains the history of the currently transported

particle.

========================= (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.

Do two USRBIN one with LET*D and the second with DOSE. Do the Dtot

normalization offline.

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.

I think the calculation is fine. If you are using the embedded scoring you

could just calculate the Energy deposited and

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).

You can't control that. The Voxels are just a way to create the geometry as

I mentioned in my first e-mail.

The COMSCW is just a weight which will be used to weight the value you need

for the scoring.

FLUKA will take care of propagating to the scoring structure which in your

case should be a USRBIN with Cartesian mesh.

But, D_total cannot know unless all interactions are finished.

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

At this stage I would suggest you to do a USERDUMP, save all the

information of the tracks you need and do the calculation offline.

Received on Mon Sep 15 2014 - 18:12:19 CEST

Date: Mon, 15 Sep 2014 23:44:28 +0900

Dear, Dr. Vittorio Boccone.

Most of all, I would like to say thank you for your specific comments.^^

It is great help for me to understand hooking action through fluka.

Based on your advices, I think there is no hooking action after end of each

"event" and each "cycle" on fluka.

I mean there is only hooking action after end of each interaction.

Ok!. as you mentioned, it is likely that off-line calculation is only way

for calc. of dose avg. LET.

I will calculate dose avg. LET by using possible ways you mentioned already.

Thank you for your help!!~~^^

Have a nice day!!~~^^

Best Regards,

Kyung-Nam Lee

From: owner-fluka-discuss_at_mi.infn.it [mailto:owner-fluka-discuss_at_mi.infn.it]

On Behalf Of Dr. Vittorio Boccone

Sent: Monday, September 15, 2014 6:47 AM

To: Kyung-Nam Lee; fluka-discuss

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

On 9/3/14 9:24 AM, Kyung-Nam Lee wrote:

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.

TRACKR common variables contains the history of the currently transported

particle.

========================= (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.

Do two USRBIN one with LET*D and the second with DOSE. Do the Dtot

normalization offline.

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.

I think the calculation is fine. If you are using the embedded scoring you

could just calculate the Energy deposited and

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).

You can't control that. The Voxels are just a way to create the geometry as

I mentioned in my first e-mail.

The COMSCW is just a weight which will be used to weight the value you need

for the scoring.

FLUKA will take care of propagating to the scoring structure which in your

case should be a USRBIN with Cartesian mesh.

But, D_total cannot know unless all interactions are finished.

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

At this stage I would suggest you to do a USERDUMP, save all the

information of the tracks you need and do the calculation offline.

(image/jpeg attachment: image001.jpg)

*
This archive was generated by hypermail 2.3.0
: Mon Sep 15 2014 - 18:12:21 CEST
*