Hi
the issue is that you asked for track length distributed binnings
for RES-NUCL scoring, hence the code interprets this request as
a track length fluence scoring which is not what you want.
If you change WHAT(1) of the USRBIN cards from 10 to 0 you'll
find that the results from your USRBIN's and your RESNUCL's estimators
are identical.
Indeed it would be advisable to stop the code if a USRBIN with RES-NUCL
is requested with scoring distributed along the track, since it is
likely not what the users wants, we'll introduce a protection for this in
the next respins/releases.
Hope this helps!
The Fluka developers
> From: Kriti Gupta <kriti.gupta_at_demcon.com>
> To: fluka discuss <fluka-discuss_at_fluka.org>
> Subject: Gas production using new FLUKA functionality in 2011.2x.8
> Date: Mon, 20 Jan 2020 14:28:46 +0000
>
> Dear FLUKA team and users,
>
> I recently asked a question about estimating light ions
production/depositi=
> on (H,D,He etc.) using USRBIN+RES-NUCL+AUXSCORE introduced in FLUKA
2011.2x=
> .8 and comparison of the results with the traditional way of estimating
iso=
> topes (RESNUCL card). I found a huge discrepancy in the numbers obtained
fr=
> om the two ways. Please find attached the input and output files for the
sa=
> me. It would be of great help if someone could point out where this
discrep=
> ancy is coming from.
>
> Thanking you,
>
> Best regards,
> Kriti Gupta
>
>
>
>
>
> -filename="I19.inp"
> TITLE
>
> * Set the defaults for precision simulations
> DEFAULTS
> PRECISIO
> *EMF-BIAS 0.001
LPBEMF
> * Define the beam characteristics
> BEAM -0.014 0.0 0.0 0.4 0.0
-1.PROTON
> * Define the beam position
> BEAMPOS 0.0 0.0 -5.
> IONTRANS -2
> PHOTONUC 1. HYDROGEN _at_LASTMAT
> PHYSICS 3.
> EVAPORAT
> PHYSICS 1.
> COALESCE
> GEOBEGIN
> COMBNAME
> 0 0
> * Black body
> SPH blkbody 0.0 0.0 0.0 100000.
> * Void sphere
> SPH void 0.0 0.0 0.0 10000.
> * Cylindrical target
> RPP degrader -1. 1. -1. 1. 0.0 0.095
> * Cylindrical target
> RPP target3 -1. 1. -1. 1. 5.095 5.105
> END
> * Black hole
> BLKBODY 5 +blkbody -void
> * Void around
> VOID 5 +void -degrader -target3
> * Degrader
> DEGRADER 5 +degrader
> * Target
> TARGET3 5 +target3
> END
> GEOEND
> MATERIAL 6. 1.88
> Graphite
> LOW-MAT Graphite 6. -3. 296.
CARBON
> MAT-PROP 20. Graphite
> DPA-ENER
> * Molybdenum-100
> MATERIAL 42. 10.22
100.MOLY100
> LOW-MAT MOLY100 42. -2. 296.
> MOLYBDEN
> * Molybdenum-98
> MATERIAL 42. 10.22
98.MOLY98
> LOW-MAT MOLY98 42. -2. 296.
> MOLYBDEN
> * Molybdenum-97
> MATERIAL 42. 10.22
97.MOLY97
> LOW-MAT MOLY97 42. -2. 296.
> MOLYBDEN
> * Molybdenum-96
> MATERIAL 42. 10.22
96.MOLY96
> LOW-MAT MOLY96 42. -2. 296.
> MOLYBDEN
> * Molybdenum-95
> MATERIAL 42. 10.22
95.MOLY95
> LOW-MAT MOLY95 42. -2. 296.
> MOLYBDEN
> * Molybdenum-94
> MATERIAL 42. 10.22
94.MOLY94
> LOW-MAT MOLY94 42. -2. 296.
> MOLYBDEN
> * Molybdenum-92
> MATERIAL 42. 10.22
92.MOLY92
> LOW-MAT MOLY92 42. -2. 296.
> MOLYBDEN
> * Material data - Molybdenum composition
> MATERIAL 10.2 MOLY
> MAT-PROP 65. MOLY
> DPA-ENER
> COMPOUND -0.95 MOLY100 -0.0404 MOLY98 -0.0026 MOLY97MOLY
> COMPOUND -0.0023 MOLY96 -0.002 MOLY95 -0.001 MOLY94MOLY
> COMPOUND -0.0017 MOLY92 MOLY
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7..
> ASSIGNMA BLCKHOLE BLKBODY
> ASSIGNMA VACUUM VOID
> ASSIGNMA Graphite DEGRADER
> ASSIGNMA MOLY TARGET3
> RESNUCLE 3. -25. 10. TARGET3
0.04resnuct
> RESNUCLE 3. -26. 10. DEGRADER
0.38resnucd
> USRBIN 10. RES-NUCL -41. 1. 1. 5.105Hdept
> USRBIN -1. -1. 5.095 1. 1. 1. &
> AUXSCORE USRBIN PROTON Hdept
> USRBIN 10. RES-NUCL -42. 1. 1. 0.095Hdepd
> USRBIN -1. -1. 0.0 1. 1. 1. &
> AUXSCORE USRBIN PROTON Hdepd
> USRBIN 10. RES-NUCL -42. 1. 1.
0.095Hedepd
> USRBIN -1. -1. 0.0 1. 1. 1. &
> AUXSCORE USRBIN 4-HELIUM Hedepd
> USRBIN 10. RES-NUCL -42. 1. 1. 0.095Ddepd
> USRBIN -1. -1. 0.0 1. 1. 1. &
> AUXSCORE USRBIN DEUTERON Ddepd
> * Set the random number seed
> RANDOMIZ 1.
> * Set the number of primary histories to be simulated in the run
> *START 300000000.
> * Set the number of primary histories to be simulated in the run
> START 300000.
> STOP
>
>
> filename: I19_25_sum.lis
> **** I19_25.rnc ****
>
>
> Total primaries run: 3000000
>
> Total weight of the primaries run: 3000000.00
>
>
> Detector n: 1 resnuct
> (Region 4 Volume: 3.99999991E-02 cmc,
> distr. type : 3 ,
> Z_max: 10 , N-Z_max: 22 , N-Z_min: -4 )
>
> Tot. response (n/cmc/pr) 4.26583327E-02 +/- 1.39916456 %
> ( --> Nuclei/pr 1.70633325E-03 +/- 1.39916456 % )
>
> **** Isotope Yield as a function of Mass Number **** ****
> (nuclei / cmc / pr) ****
>
>
> A_min: 1 - A_max: 42
>
> A: 1 4.2658333E-02 +/- 1.3991646E+00 %
>
> **** Isotope Yield as a function of Atomic Number **** ****
> (nuclei / cmc / pr) ****
>
>
> Z_min: 1 - Z_max: 10
>
> Z: 1 4.2658333E-02 +/- 1.3991646E+00 %
>
> **** Residual nuclei distribution **** **** (nuclei / cmc / pr)
> ****
>
>
>
> A \ Z 1 2 3 4 5 6
> 7 8 9 10
> 1 4.27E-02 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00
> 0.00E+00 0.00E+00 0.00E+00 0.00E+00
> +/- 1.4 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/-
> 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 %
>
>
> :
> filename : I19_26_sum.lis
> **** I19_26.rnc ****
>
>
> Total primaries run: 3000000
>
> Total weight of the primaries run: 3000000.00
>
>
> Detector n: 1 resnucd
> (Region 3 Volume: 0.379999995 cmc,
> distr. type : 3 ,
> Z_max: 10 , N-Z_max: 6 , N-Z_min: -4 )
>
> Tot. response (n/cmc/pr) 2.82517541E-02 +/- 1.09933329 %
> ( --> Nuclei/pr 1.07356664E-02 +/- 1.09933329 % )
>
> **** Isotope Yield as a function of Mass Number **** ****
> (nuclei / cmc / pr) ****
>
>
> A_min: 1 - A_max: 26
>
> A: 13 9.8245614E-05 +/- 1.1739930E+01 %
> A: 12 9.2956140E-03 +/- 1.0847383E+00 %
> A: 10 4.2982456E-05 +/- 1.1958092E+01 %
> A: 6 1.2280702E-05 +/- 2.4281045E+01 %
> A: 4 9.1228071E-03 +/- 1.8669147E+00 %
> A: 2 1.2280702E-05 +/- 3.4006802E+01 %
> A: 1 9.6675440E-03 +/- 1.0937303E+00 %
>
> **** Isotope Yield as a function of Atomic Number **** ****
> (nuclei / cmc / pr) ****
>
>
> Z_min: 1 - Z_max: 10
>
> Z: 7 6.1403510E-06 +/- 3.0491068E+01 %
> Z: 6 9.3877192E-03 +/- 1.1333162E+00 %
> Z: 5 4.2982456E-05 +/- 1.1958092E+01 %
> Z: 3 1.2280702E-05 +/- 2.4281045E+01 %
> Z: 2 9.1228071E-03 +/- 1.8669147E+00 %
> Z: 1 9.6798247E-03 +/- 1.0875541E+00 %
>
> **** Residual nuclei distribution **** **** (nuclei / cmc / pr)
> ****
>
>
>
> A \ Z 1 2 3 4 5 6
> 7 8 9 10
> 13 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 9.21E-05
> 6.14E-06 0.00E+00 0.00E+00 0.00E+00
> +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/-12.2 %
> +/-30.5 % +/- 0.0 % +/- 0.0 % +/- 0.0 %
> 12 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 9.30E-03
> 0.00E+00 0.00E+00 0.00E+00 0.00E+00
> +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 1.1 %
> +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 %
> 10 0.00E+00 0.00E+00 0.00E+00 0.00E+00 4.30E-05 0.00E+00
> 0.00E+00 0.00E+00 0.00E+00 0.00E+00
> +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/-12.0 % +/- 0.0 %
> +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 %
> 6 0.00E+00 0.00E+00 1.23E-05 0.00E+00 0.00E+00 0.00E+00
> 0.00E+00 0.00E+00 0.00E+00 0.00E+00
> +/- 0.0 % +/- 0.0 % +/-24.3 % +/- 0.0 % +/- 0.0 % +/- 0.0 %
> +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 %
> 4 0.00E+00 9.12E-03 0.00E+00 0.00E+00 0.00E+00 0.00E+00
> 0.00E+00 0.00E+00 0.00E+00 0.00E+00
> +/- 0.0 % +/- 1.9 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 %
> +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 %
> 2 1.23E-05 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00
> 0.00E+00 0.00E+00 0.00E+00 0.00E+00
> +/-34.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 %
> +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 %
> 1 9.67E-03 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00
> 0.00E+00 0.00E+00 0.00E+00 0.00E+00
> +/- 1.1 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 %
> +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 %
>
>
> filename="I19_42.bnn.lis"; size=1896;
>
> 1
> Cartesian binning n. 1 "Hdepd " , generalized particle n. 251
> X coordinate: from -1.0000E+00 to 1.0000E+00 cm, 1 bins (
> 2.0000E+00 cm wide)
> Y coordinate: from -1.0000E+00 to 1.0000E+00 cm, 1 bins (
> 2.0000E+00 cm wide)
> Z coordinate: from 0.0000E+00 to 9.5000E-02 cm, 1 bins (
> 9.5000E-02 cm wide)
> Data follow in a matrix A(ix,iy,iz), format (1(5x,1p,10(1x,e11.4)))
>
> accurate deposition along the tracks requested
> this is a track-length binning
> 9.5404E-07
>
> Percentage errors follow in a matrix A(ix,iy,iz), format
> (1(5x,1p,10(1x,e11.4)))
>
> 1.0963E+00
> 1
> Cartesian binning n. 2 "Hedepd " , generalized particle n. 251
> X coordinate: from -1.0000E+00 to 1.0000E+00 cm, 1 bins (
> 2.0000E+00 cm wide)
> Y coordinate: from -1.0000E+00 to 1.0000E+00 cm, 1 bins (
> 2.0000E+00 cm wide)
> Z coordinate: from 0.0000E+00 to 9.5000E-02 cm, 1 bins (
> 9.5000E-02 cm wide)
> Data follow in a matrix A(ix,iy,iz), format (1(5x,1p,10(1x,e11.4)))
>
> accurate deposition along the tracks requested
> this is a track-length binning
> 8.9434E-07
>
> Percentage errors follow in a matrix A(ix,iy,iz), format
> (1(5x,1p,10(1x,e11.4)))
>
> 1.8249E+00
>
> 1
> Cartesian binning n. 3 "Ddepd " , generalized particle n. 251
> X coordinate: from -1.0000E+00 to 1.0000E+00 cm, 1 bins (
> 2.0000E+00 cm wide)
> Y coordinate: from -1.0000E+00 to 1.0000E+00 cm, 1 bins (
> 2.0000E+00 cm wide)
> Z coordinate: from 0.0000E+00 to 9.5000E-02 cm, 1 bins (
> 9.5000E-02 cm wide)
> Data follow in a matrix A(ix,iy,iz), format (1(5x,1p,10(1x,e11.4)))
>
> accurate deposition along the tracks requested
> this is a track-length binning
> 1.2211E-09
>
> Percentage errors follow in a matrix A(ix,iy,iz), format
> (1(5x,1p,10(1x,e11.4)))
>
> 3.4017E+01
>
>
> __________________________________________________________________________
> You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
>
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Tue Jan 21 2020 - 14:21:34 CET