Re: 131-Iodine

attached mail follows:

Dear Alfredo,
we compared the emission spectrum you reported with the one we used,
taken from
http://www.doseinfo-radar.com/RADARDecay.html

We found some difference in the CE brabching ratios, and only minor
differences in betas and gammas. Even if a simulation should be done to
compare the possible impact of such slight differences, we don't think
that this could be the main reason for our problem.

In Fluka, we used the .inp file and the source.f file attached.

We read the energy deposited in each region in the .out file attached
(path through Flair interface: Scoring -> Events by region).

Furthermore, the heading report for Energy the unit "GeV/cm**3", but I
think it should be "GeV". In fact, if I try to extend the target region
from 1 m^3 to 8 m^3 the Energy remains the same, as it should be if it
is expressed as energy deposited (GeV), and not energy density
(GeV/cm^3).
For this test, I intentionally enlarged the target region, in order to
have in both cases an absorbed fraction almost unitary.

Thank you,

Ernesto

On mer, 2009-09-23 at 17:02 +0200, Alfredo Ferrari wrote:
> Hi
>
> the emission spectra used by FLUKA for 131-I are reported below. At a
> first glance I don't see anything strange, but I let you check with the
> relevant databases if there is something suspicious (normally our
> data are adapted from the Brookhaven database). I have a suspicion (maybe
> just a suspicion) about a possible misunderstading, could you please send
> me the input and explain how you scored the energy?
>
> Ciao!
> Alfredo
>
> A: 131 Sym.: I Z: 53 m: 0 T 1/2 6.947E+05 (s)
>
> Decay ch. # BR into A Sym. Z m T 1/2 (s)
> B- 1 0.9891 -> 131 Xe 54 0 1.000E+38
> B-m 2 0.0109 -> 131 Xe 54 1 1.028E+06
>
> 15 CE lines recorded:
> Energy: 0.003 MeV, BR: 7.230 %
> Energy: 0.025 MeV, BR: 0.830 %
> Energy: 0.046 MeV, BR: 4.250 %
> Energy: 0.075 MeV, BR: 0.564 %
> Energy: 0.079 MeV, BR: 0.115 %
> Energy: 0.129 MeV, BR: 0.818 %
> Energy: 0.143 MeV, BR: 0.062 %
> Energy: 0.158 MeV, BR: 0.382 %
> Energy: 0.163 MeV, BR: 0.088 %
> Energy: 0.250 MeV, BR: 0.306 %
> Energy: 0.279 MeV, BR: 0.054 %
> Energy: 0.330 MeV, BR: 1.910 %
> Energy: 0.359 MeV, BR: 0.299 %
> Energy: 0.363 MeV, BR: 0.062 %
> Energy: 0.602 MeV, BR: 0.035 %
>
> 18 gamma lines recorded:
> Energy: 0.004 MeV, BR: 0.780 %
> Energy: 0.029 MeV, BR: 1.880 %
> Energy: 0.030 MeV, BR: 3.450 %
> Energy: 0.034 MeV, BR: 0.323 %
> Energy: 0.034 MeV, BR: 0.624 %
> Energy: 0.034 MeV, BR: 0.189 %
> Energy: 0.080 MeV, BR: 2.620 %
> Energy: 0.177 MeV, BR: 0.269 %
> Energy: 0.272 MeV, BR: 0.058 %
> Energy: 0.284 MeV, BR: 6.120 %
> Energy: 0.318 MeV, BR: 0.077 %
> Energy: 0.326 MeV, BR: 0.273 %
> Energy: 0.364 MeV, BR: 81.500 %
> Energy: 0.405 MeV, BR: 0.055 %
> Energy: 0.503 MeV, BR: 0.359 %
> Energy: 0.637 MeV, BR: 7.160 %
> Energy: 0.643 MeV, BR: 0.217 %
> Energy: 0.723 MeV, BR: 1.770 %
>
> 6 beta+/beta- emissions recorded:
> B- end point (<E>): 0.248 ( 0.069) MeV, BR: 2.080 %
> B- end point (<E>): 0.304 ( 0.087) MeV, BR: 0.645 %
> B- end point (<E>): 0.334 ( 0.097) MeV, BR: 7.230 %
> B- end point (<E>): 0.606 ( 0.192) MeV, BR: 89.600 %
> B- end point (<E>): 0.630 ( 0.200) MeV, BR: 0.050 %
> B- end point (<E>): 0.807 ( 0.283) MeV, BR: 0.390 %
>
>
>
> +----------------------------------------------------------------------+
> | Alfredo Ferrari || Tel.: +41.22.76.76119 |
> | CERN-EN/STI || Fax.: +41.22.76.69474 |
> | 1211 Geneva 23 || e-mail: Alfredo.Ferrari_at_cern.ch |
> | Switzerland || |
> +----------------------------------------------------------------------+

*$ CREATE SOURCE.FOR
*COPY SOURCE
*
*=== source ===========================================================*
*
      SUBROUTINE SOURCE ( NOMORE )

      INCLUDE '(DBLPRC)'
      INCLUDE '(DIMPAR)'
      INCLUDE '(IOUNIT)'
*
*----------------------------------------------------------------------*
* *
* Copyright (C) 1990-2009 by Alfredo Ferrari & Paola Sala *
* All Rights Reserved. *
* *
* *
* New source for FLUKA9x-FLUKA20xy: *
* *
* Created on 07 january 1990 by Alfredo Ferrari & Paola Sala *
* Infn - Milan *
* *
* Last change on 08-feb-09 by Alfredo Ferrari *
* *
* This is just an example of a possible user written source routine. *
* note that the beam card still has some meaning - in the scoring the *
* maximum momentum used in deciding the binning is taken from the *
* beam momentum. Other beam card parameters are obsolete. *
* *
* Output variables: *
* *
* Nomore = if > 0 the run will be terminated *
* *
*----------------------------------------------------------------------*
*
      INCLUDE '(BEAMCM)'
      INCLUDE '(FHEAVY)'
      INCLUDE '(FLKSTK)'
      INCLUDE '(IOIOCM)'
      INCLUDE '(LTCLCM)'
      INCLUDE '(PAPROP)'
      INCLUDE '(SOURCM)'
      INCLUDE '(SUMCOU)'
*
      LOGICAL LFIRST
*
      SAVE LFIRST
      DATA LFIRST / .TRUE. /
*======================================================================*
* *
* BASIC VERSION *
* *
*======================================================================*
      NOMORE = 0
* +-------------------------------------------------------------------*
* | First call initializations:
      IF ( LFIRST ) THEN
* | *** The following 3 cards are mandatory ***
         TKESUM = ZERZER
         LFIRST = .FALSE.
         LUSSRC = .TRUE.
* | *** User initialization ***
      END IF
* |
* +-------------------------------------------------------------------*
* Push one source particle to the stack. Note that you could as well
* push many but this way we reserve a maximum amount of space in the
* stack for the secondaries to be generated
* Npflka is the stack counter: of course any time source is called it
* must be =0
      NPFLKA = NPFLKA + 1
* Wt is the weight of the particle
      WTFLK (NPFLKA) = ONEONE
      WEIPRI = WEIPRI + WTFLK (NPFLKA)
* Particle type (1=proton.....). Ijbeam is the type set by the BEAM
* card
* +-------------------------------------------------------------------*
* | (Radioactive) isotope:
      IF ( IJBEAM .EQ. -2 .AND. LRDBEA ) THEN
         IARES = IPROA
         IZRES = IPROZ
         IISRES = IPROM
         CALL STISBM ( IARES, IZRES, IISRES )
         IJHION = IPROZ * 1000 + IPROA
         IJHION = IJHION * 100 + KXHEAV
         IONID = IJHION
         CALL DCDION ( IONID )
         CALL SETION ( IONID )
* |
* +-------------------------------------------------------------------*
* | Heavy ion:
      ELSE IF ( IJBEAM .EQ. -2 ) THEN
         IJHION = IPROZ * 1000 + IPROA
         IJHION = IJHION * 100 + KXHEAV
         IONID = IJHION
         CALL DCDION ( IONID )
         CALL SETION ( IONID )
         ILOFLK (NPFLKA) = IJHION
* | Flag this is prompt radiation
         LRADDC (NPFLKA) = .FALSE.
* | Group number for "low" energy neutrons, set to 0 anyway
         IGROUP (NPFLKA) = 0
* |
* +-------------------------------------------------------------------*
* | Normal hadron:
      ELSE
         IONID = IJBEAM
         ILOFLK (NPFLKA) = IJBEAM
* | Flag this is prompt radiation
         LRADDC (NPFLKA) = .FALSE.
* | Group number for "low" energy neutrons, set to 0 anyway
         IGROUP (NPFLKA) = 0
      END IF
* |
* +-------------------------------------------------------------------*
* From this point .....
* Particle generation (1 for primaries)
      LOFLK (NPFLKA) = 1
* User dependent flag:
      LOUSE (NPFLKA) = 0
* User dependent spare variables:
      DO 100 ISPR = 1, MKBMX1
         SPAREK (ISPR,NPFLKA) = ZERZER
 100 CONTINUE
* User dependent spare flags:
      DO 200 ISPR = 1, MKBMX2
         ISPARK (ISPR,NPFLKA) = 0
 200 CONTINUE
* Save the track number of the stack particle:
      ISPARK (MKBMX2,NPFLKA) = NPFLKA
      NPARMA = NPARMA + 1
      NUMPAR (NPFLKA) = NPARMA
      NEVENT (NPFLKA) = 0
      DFNEAR (NPFLKA) = +ZERZER
* ... to this point: don't change anything
* Particle age (s)
      AGESTK (NPFLKA) = +ZERZER
      AKNSHR (NPFLKA) = -TWOTWO
* Kinetic energy of the particle (GeV)
      TKEFLK (NPFLKA) = SQRT ( PBEAM**2 + AM (IONID)**2 ) - AM (IONID)
* Particle momentum
      PMOFLK (NPFLKA) = PBEAM
* PMOFLK (NPFLKA) = SQRT ( TKEFLK (NPFLKA) * ( TKEFLK (NPFLKA)
* & + TWOTWO * AM (IONID) ) )
* Cosines (tx,ty,tz)
      TXFLK (NPFLKA) = UBEAM
      TYFLK (NPFLKA) = VBEAM
      TZFLK (NPFLKA) = WBEAM
* TZFLK (NPFLKA) = SQRT ( ONEONE - TXFLK (NPFLKA)**2
* & - TYFLK (NPFLKA)**2 )
* Polarization cosines:
      TXPOL (NPFLKA) = -TWOTWO
      TYPOL (NPFLKA) = +ZERZER
      TZPOL (NPFLKA) = +ZERZER
* Particle coordinates
      XFLK (NPFLKA) = 0.6*FLRNDM(XD)
      YFLK (NPFLKA) = 0.6*FLRNDM(YD)
      ZFLK (NPFLKA) = 0.6*FLRNDM(ZD)
* Calculate the total kinetic energy of the primaries: don't change
      IF ( ILOFLK (NPFLKA) .EQ. -2 .OR. ILOFLK (NPFLKA) .GT. 100000 )
     & THEN
         TKESUM = TKESUM + TKEFLK (NPFLKA) * WTFLK (NPFLKA)
      ELSE IF ( ILOFLK (NPFLKA) .NE. 0 ) THEN
         TKESUM = TKESUM + ( TKEFLK (NPFLKA) + AMDISC (ILOFLK(NPFLKA)) )
     & * WTFLK (NPFLKA)
      ELSE
         TKESUM = TKESUM + TKEFLK (NPFLKA) * WTFLK (NPFLKA)
      END IF
      RADDLY (NPFLKA) = ZERZER
* Here we ask for the region number of the hitting point.
* NREG (NPFLKA) = ...
* The following line makes the starting region search much more
* robust if particles are starting very close to a boundary:
      CALL GEOCRS ( TXFLK (NPFLKA), TYFLK (NPFLKA), TZFLK (NPFLKA) )
      CALL GEOREG ( XFLK (NPFLKA), YFLK (NPFLKA), ZFLK (NPFLKA),
     & NRGFLK(NPFLKA), IDISC )
* Do not change these cards:
      CALL GEOHSM ( NHSPNT (NPFLKA), 1, -11, MLATTC )
      NLATTC (NPFLKA) = MLATTC
      CMPATH (NPFLKA) = ZERZER
      CALL SOEVSV
      RETURN
*=== End of subroutine Source =========================================*
      END