Date: Fri, 28 Feb 2014 10:14:09 -0800
Dear Francesco,
I have one more question which is on assigning volumes to region and not
specific to scoring DOSE.
I am not clear what "write the volumes (in fixed format: 10 fields per
region, 7 regions per row)" means and I can't find it in the manual.
See the images for the GEOBEGIN and REGION cards. It is true that flair
1.1-3 doesn't show the 680 cm3 as the volume for the LIC_oc region but
that is not what we were referring to. We were referring to the results
in the .out files where the LIC_oc region volume is reported as 1 under
the region summary report:
1Region # name volume ENERGY Density EM-ENRGY
Density Star Density Star Density
cm**3 GeV/cm**3
GeV/cm**3 Stars/cm**3 Stars/cm**3
/one beam
particle /one beam particle /one beam particle /one beam
particle
8 LIC_ic 1.0000D+00 1.250409D-06
1.250409D-06 0.0000D+00 0.00000D+00
I have truncated the numbers to make things fit in one line. Are we
supposed to see 680 under the volume column? If so, what did we do wrong?
Thank you and best wishes,
Mina
On 14-02-27 02:02 AM, Francesco Cerutti wrote:
>
> Hallo Mina & Martin,
>
> Coming now to your case 2b, if you follow the manual and set IVLFLG =
> 3 in the geometry title card and write the volumes (in fixed format:
> 10 fields per region, 7 regions per row) before the GEOEND card (after
> the region END card), everything works as expected: the SCORE values
> will get normalized by your volumes and the latter ones will appear in
> the volume column. You still got default volume values of 1cm3 due to
> a Flair bug, not translating correctly what was input in the region
> metacard under the Volume label. This is going to be fixed in the next
> Flair version (by the way, Flair 2 is coming!), but it gives me the
> nasty opportunity to remind people that in case of problems and
> unexpected behaviors one should look at the input file as actually fed
> to FLUKA, leaving for a little while the wonderful Flair world and
> using a trivial text editor to inspect what is underneath.
>
> Best wishes
>
> Francesco
>
> **************************************************
> Francesco Cerutti
> CERN-EN/STI
> CH-1211 Geneva 23
> Switzerland
> tel. ++41 22 7678962
> fax ++41 22 7668854
>
> On Thu, 27 Feb 2014, Mina Nozar wrote:
>
>> Hello everyone,
>>
>> We are trying to score Absorbed Dose in a given region. We have
>> tried two scoring methods:
>> 1) via a USRBIN card with REGION binning and DOSE as the quantity scored
>> and
>> 2) via a SCORE card and Energy as the quantity scored, with no volume
>> input (case 2a) and with volume input (case 2b).
>>
>> According to the manual, we should get
>> - GeV/g per primary for option 1
>> - GeV per primary for option 2a
>> - GeV/cm3 per primary for option 2b
>>
>>
>> As for the volume setting for the SCORE card, we set IVFLG to 3 in
>> the GEOBIN title card and inputted the region volume in the geometry
>> region card for the region.
>>
>>
>> We are seeing some discrepancies. The value we get from SCORE (with
>> no volume setting, case 2a) agrees with the value from the USRBIN, if
>> we divide the SCORE value by the density. This implies that the
>> SCORE value is GeV/cm3 per primary which doesn't agree with what is
>> in the manual.
>>
>> Furthermore, when we do set the volume, following the instructions
>> above, we still see a value of 1.000000000D+00 for the region's
>> volume and the same deposited energy value as in case 2a.
>>
>> This is a source of confusion for us and we are eager to find whether
>> we are missing something. Can someone shed some light on this please?
>> Is there another way to score absorbed dose in a region?
>>
>> Thank you very much,
>> Mina & Martin
>>
>>
>
(image/jpeg attachment: GEOBEGIN_card.jpg)
(image/jpeg attachment: REGION_card.jpg)