From: Paola Sala (Paola.Sala@cern.ch)
Date: Mon Mar 20 2006 - 13:09:11 CET
Hi,
the problem is in the meaning of
the variable IJ in comscw: it
refers to the "generalized" particle type,
that means that it is set to 208 ( or 211 for e.m.) when
the quantity to be scored is energy deposition.
You should use the JTRACK variable from common trackr
to select energy deposition by protons.
Paola
-----Original Message-----
From: owner-fluka-discuss@fisica.unimi.it on behalf of Fan, Ph.D., Jiajin
Sent: Fri 3/17/2006 12:14 AM
To: fluka-discuss@mi.infn.it
Subject: Problems about weighting deposited energy
Hi All,
I want to calculate proton dose using USRBIN and weight deposited energy by
using COMSCW.f. But I can not get results as I expected.
First, according to FLUKA manual, parameter ISCRNG in COMSCW.f should be 1
when we calculate energy density bin (208.), but in my calculation it always
be 2, I can not find out why it's not consist with manual.
Second, I change COMSCW to different values but always got same output in
USRBIN.
Is there anything wrong in my input file or COMSCW.f file?
Thank you.
In my input file, define USERWEIG and USRBIN as:
USERWEIG 0. 0. 4. 0. 0. 1.0
USRBIN 10.0 208. 60.0 .5 .5 100.
ENERGY
USRBIN -.5 -.5 0.0 1. 1. 1000.
&
And the contents in COMSCW.f are:
* ======== In order to compute RBE included dose ======== *
if ( ij .eq. 1 ) then
enerProton = etrack - 0.93827231
write(*,*) iscrng
if ( iscrng .eq. 1 ) then
if ( enerProton .lt. 0.1) COMSCW = 12.0
end if
end if
*******************************************
Jiajin Fan, Ph.D.
Research Associate, Radiation Oncology
Fox Chase Cancer Center
432 Rhawn Street, Philadelphia, PA 19111
Tel: (215) 214-1711
Fax:(215) 728-1148
Email: jiajin.fan@fccc.edu
*******************************************
This archive was generated by hypermail 2.1.6 : Mon Mar 20 2006 - 13:52:23 CET