From: Markus Kirk (m.kirk@gsi.de)
Date: Wed Jun 13 2007 - 16:07:16 CEST
Dear Fluka Users,
Why does my SOURCE routine cause an abort during a run when the
primaries are ions (e.g. here A=6,Z=3) ?
There is however no problem when the primaries are protons.
Please see the attached files...
The SOURCE describes a beam with an energy distribution
sampled from the probability distribution read from the file BETABE.
All other attributes (X,Y,Z and direction cosines) are copied from those
variables set by the BEAM and BEAMPOS cards.
Any ideas?
I thank you in anticipation.
Yours sincerely,
Markus
100
0.00000d0 8.28639d-06
0.0151515d0 1.05402d-05
0.0303030d0 1.33580d-05
0.0454545d0 1.68670d-05
0.0606061d0 2.12196d-05
0.0757576d0 2.65976d-05
0.0909091d0 3.32164d-05
0.106061d0 4.13302d-05
0.121212d0 5.12374d-05
0.136364d0 6.32865d-05
0.151515d0 7.78825d-05
0.166667d0 9.54935d-05
0.181818d0 0.000116657d0
0.196970d0 0.000141989d0
0.212121d0 0.000172189d0
0.227273d0 0.000208045d0
0.242424d0 0.000250447d0
0.257576d0 0.000300385d0
0.272727d0 0.000358959d0
0.287879d0 0.000427383d0
0.303030d0 0.000506984d0
0.318182d0 0.000599206d0
0.333333d0 0.000705607d0
0.348485d0 0.000827855d0
0.363636d0 0.000967722d0
0.378788d0 0.00112707d0
0.393939d0 0.00130785d0
0.409091d0 0.00151206d0
0.424242d0 0.00174174d0
0.439394d0 0.00199896d0
0.454545d0 0.00228575d0
0.469697d0 0.00260410d0
0.484848d0 0.00295592d0
0.500000d0 0.00334297d0
0.515152d0 0.00376683d0
0.530303d0 0.00422888d0
0.545455d0 0.00473019d0
0.560606d0 0.00527154d0
0.575758d0 0.00585330d0
0.590909d0 0.00647543d0
0.606061d0 0.00713743d0
0.621212d0 0.00783826d0
0.636364d0 0.00857634d0
0.651515d0 0.00934952d0
0.666667d0 0.0101550d0
0.681818d0 0.0109895d0
0.696970d0 0.0118490d0
0.712121d0 0.0127288d0
0.727273d0 0.0136238d0
0.742424d0 0.0145283d0
0.757576d0 0.0154360d0
0.772727d0 0.0163403d0
0.787879d0 0.0172342d0
0.803030d0 0.0181103d0
0.818182d0 0.0189612d0
0.833333d0 0.0197793d0
0.848485d0 0.0205571d0
0.863636d0 0.0212870d0
0.878788d0 0.0219621d0
0.893939d0 0.0225756d0
0.909091d0 0.0231210d0
0.924242d0 0.0235929d0
0.939394d0 0.0239861d0
0.954545d0 0.0242964d0
0.969697d0 0.0245206d0
0.984848d0 0.0246560d0
1.00000d0 0.0247014d0
1.01515d0 0.0246560d0
1.03030d0 0.0245206d0
1.04545d0 0.0242964d0
1.06061d0 0.0239861d0
1.07576d0 0.0235929d0
1.09091d0 0.0231210d0
1.10606d0 0.0225756d0
1.12121d0 0.0219621d0
1.13636d0 0.0212870d0
1.15152d0 0.0205571d0
1.16667d0 0.0197793d0
1.18182d0 0.0189612d0
1.19697d0 0.0181103d0
1.21212d0 0.0172342d0
1.22727d0 0.0163403d0
1.24242d0 0.0154360d0
1.25758d0 0.0145283d0
1.27273d0 0.0136238d0
1.28788d0 0.0127288d0
1.30303d0 0.0118490d0
1.31818d0 0.0109895d0
1.33333d0 0.0101550d0
1.34848d0 0.00934952d0
1.36364d0 0.00857634d0
1.37879d0 0.00783826d0
1.39394d0 0.00713743d0
1.40909d0 0.00647544d0
1.42424d0 0.00585330d0
1.43939d0 0.00527154d0
1.45455d0 0.00473019d0
1.46970d0 0.00422888d0
1.48485d0 0.00376683d0
1.50000d0 0.00334297d0
*
* M. Kirk
*
* 24 April 2007
*
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
TITLE
EURISOL DS Beta Beams. Hypothetical geometry. Worser case than expected.
**DEFAULTS NEW-DEFA
**PHYSICS 2.0 EVAPORAT
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+...
BEAM -1.5 0.0 0.0 10.0 10.0 1.PROTON
**BEAM -1.5 10.0 10.0 1.HEAVYION
* Primary beam is X(A=WHAT2,Z=WHAT1)
**HI-PROPE 3.0 6.0 0.0 0.0 0.0
*
BEAMPOS 0.0 0.0 -100.0 0.0 0.0
*
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+...
* Open the file containing the formatted SOURCE data which is read in the
* SOURCE card below. WHAT1 is the unit number of the file.
* A positive unit means formatted file is to be expected.
* The user has to make sure that the given unit number is not alread in use!
OPEN 29. READONLY
/u/kirk/codes/fluka/fluka_work/BETABE
*
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
* The user written source routine source_cpt.f is called here.
* NOTE: the unit number of the file to read must be the same as that in OPEN.
SOURCE 29. 1.0
*
**EVENTYPE 0. 0. 2. 0. 0. 0.DPMJET
*
GEOBEGIN COMBINAT
CERN PS dipole with collimator & beampipe
RPP 1-5000000.0+5000000.0-5000000.0+5000000.0-5000000.0+5000000.0
RPP 2-1000000.0+1000000.0-1000000.0+1000000.0-1000000.0+1000000.0
* Collimator
RPP 3 -30.00 +30.00 -30.00 +30.00 -50.00 +0.00
* Coil insulation
RPP 4 -30.00 +30.00 -30.00 +30.00 +0.00 +1.00
END
*
* REGIONS
*
*234x6789x123456x89 123x56789 x234567x9 1234x6789 1x345678x 12345x789 12x stop at 73
* Black hole
001 5 +1 -2
* Space around targets
002 5 +2 -3 -4
* Collimator
003 5 +3
* Insulation
004 5 +4
*
END
*
GEOEND
*
* #############################################################################
* ############################## MAKE MATERIALS ###############################
* #############################################################################
*
*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
*
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
MATERIAL 1.0 1.0079 0.000089 3.0 0.0 0. HYDROGEN
*
* The following is just an override of Oxygen, keeping the default settings...
MATERIAL 8.0 15.9994 0.00133 8.0 0.0 0. OXYGEN
*
MATERIAL 14.0 28.086 2.33 26.0 0.0 0. SILICON
*
MATERIAL 15.0 30.97376 2.2 27.0 0.0 0. PHOSPHO
*
MATERIAL 16.0 32.066 2.07 28.0 0.0 0. SULFUR
*
MATERIAL 22.0 47.867 4.50 29.0 0.0 0. TITANIUM
*
MATERIAL 24.0 51.9961 7.19 30.0 0.0 0. CHROMIUM
*
MATERIAL 28.0 58.6934 7.19 31.0 0.0 0. NICKEL
*
* Argon material n. 32, normal temperature...
MATERIAL 18.0 39.948 1.400 32.0 0.0 0. ARGON
*
* Boron material n. 33.
MATERIAL 5.0 10.8117 2.460 33.0 0.0 0. BORON
*
* ---------- MAT#34 austenitischer Stahl X6CrNiTi1810 Dichte 7.62 ----------
MATERIAL 0.0 0.0 7.62 34.0 0.0 0. AUSTEN
COMPOUND -1.48E-7 6.0 -5.20E-3 26.0 -2.58E-5 27. AUSTEN
COMPOUND -1.78E-4 28.0 -3.55E-3 29.0 -1.73E-1 30. AUSTEN
COMPOUND -7.03E-1 11.0 -1.14E-1 31.0 0.0 0. AUSTEN
*
* ---------------------------- Water MAT#35 --------------------------------
*
MATERIAL 0.0 0.0 1.0 35.0 0.0 0. WATER
COMPOUND 2.0 3.0 1.0 8.0 0.0 0. WATER
*
* ----------------------------- Air MAT#36 ---------------------------------
*
* Air (mat. 36): 3 components, N (mat. 7), O (mat. 8), gaseous Ar (mat.32)
*
*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
MATERIAL 0.0 0.0 0.001225 36.0 0.0 0. AIR
COMPOUND -0.9256E-3 7.0 -0.2837E-3 8.0 -1.572E-5 32. AIR
*
* ------------------------- Type E Glass MAT#37 ----------------------------
*
* The glass used in the insulation for the coils (epoxy fibre glass) is of type E.
* Glass is the main composition of the insulation (~80%)...so assume 100%
* The chemical compositon of type E glass is:
*
* Table 1.
* Substance Mass%
*----------------------------------------------------
* SiO2 52-56% ...take mid value 54% !!
* Alkaline Oxides 0-2% ... 1%
* CaO 16-25% ... 20.5% !!
* MgO 0-5% ... 2.5% !!
* B2O3 5-10% ... 7.5% !!
* Al2O3 12-16% ... 14% !!
* TiO2 0-0.8% ... 0.4%
* Fe2O3 0.05-0.4% ... 0.225%
* F2 0-1% ... 0.5%
*----------------------------------------------------
*
* Break table 1 into the individual effective mass fractions of the most
* significant elements marked in table 1 with a "!!". We thus obtain...
*
* Si 0.2974
* O 0.4839
* Ca 0.1732
* Mg 0.0180
* B 0.0276
*
* Araldite MY 740 Liquid, unmodified, solvent-free Bisphenol A epoxy resin
* Density at 25 degC 1.15-1.20 g/cm3 (Take mid value 1.175 g/cm3)
*
*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
MATERIAL 0.0 0.0 1.175 37.0 0.0 0. GLASS
COMPOUND -0.2974 14.0 -0.4839 8.0 -0.1732 21. GLASS
COMPOUND -0.0180 9.0 -0.0276 33.0 0.0 0. GLASS
*
*
* #############################################################################
* ####################### ASSIGN MATERIALS TO REGIONS #########################
* #############################################################################
*
* NOTE!!!
* The regions INDEXES specified in WHAT2- correspond to the strings which label
* with 3 characters, those regions defined in the geometry above. Thus take
* GREAT care when commenting out or reintroducing a region that the INDEXES are
* reordered back to a CONTINUOUS numerical order! It is also best advised to
* keep the labels in the geometry above CONSISTENT with the corresponding index.
* E.g. 001 corresponds ro region index 1.0 and so on.
*
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
*
* Black hole to region 001
ASSIGNMA 1.0 1.0
* Air to region 002. Space around targets.
ASSIGNMA 36.0 2.0
*
* Assign C to region 003. Collimator
ASSIGNMA 6.0 3.0
*
* Assign Insulation #37 to region 004.
ASSIGNMA 37.0 4.0
*
*
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
*
* MESH BINNING. PHYSICAL DOSE
* ===========================
*
* Coil insulation Unit Xmax Ymax Zmax
USRBIN 10.0 208.0 42.0 30.00 30.00 1.00 DI40
USRBIN -30.00 -30.00 0.00 1.0 1.0 1.0 &
*
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
RANDOMIZE 1.0
*
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
START 1.0E0
STOP
*
* M. Kirk
*
* 24 April 2007
*
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
TITLE
EURISOL DS Beta Beams. Hypothetical geometry. Worser case than expected.
DEFAULTS NEW-DEFA
PHYSICS 2.0 EVAPORAT
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+...
**BEAM -1.5 0.0 0.0 10.0 10.0 1.PROTON
BEAM -1.5 10.0 10.0 1.HEAVYION
* Primary beam is X(A=WHAT2,Z=WHAT1)
HI-PROPE 3.0 6.0 0.0 0.0 0.0
*
BEAMPOS 0.0 0.0 -100.0 0.0 0.0
*
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+...
* Open the file containing the formatted SOURCE data which is read in the
* SOURCE card below. WHAT1 is the unit number of the file.
* A positive unit means formatted file is to be expected.
* The user has to make sure that the given unit number is not alread in use!
OPEN 29. READONLY
/u/kirk/codes/fluka/fluka_work/BETABE
*
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
* The user written source routine source_cpt.f is called here.
* NOTE: the unit number of the file to read must be the same as that in OPEN.
SOURCE 29. 1.0
*
EVENTYPE 0. 0. 2. 0. 0. 0.DPMJET
*
GEOBEGIN COMBINAT
CERN PS dipole with collimator & beampipe
RPP 1-5000000.0+5000000.0-5000000.0+5000000.0-5000000.0+5000000.0
RPP 2-1000000.0+1000000.0-1000000.0+1000000.0-1000000.0+1000000.0
* Collimator
RPP 3 -30.00 +30.00 -30.00 +30.00 -50.00 +0.00
* Coil insulation
RPP 4 -30.00 +30.00 -30.00 +30.00 +0.00 +1.00
END
*
* REGIONS
*
*234x6789x123456x89 123x56789 x234567x9 1234x6789 1x345678x 12345x789 12x stop at 73
* Black hole
001 5 +1 -2
* Space around targets
002 5 +2 -3 -4
* Collimator
003 5 +3
* Insulation
004 5 +4
*
END
*
GEOEND
*
* #############################################################################
* ############################## MAKE MATERIALS ###############################
* #############################################################################
*
*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
*
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
MATERIAL 1.0 1.0079 0.000089 3.0 0.0 0. HYDROGEN
*
* The following is just an override of Oxygen, keeping the default settings...
MATERIAL 8.0 15.9994 0.00133 8.0 0.0 0. OXYGEN
*
MATERIAL 14.0 28.086 2.33 26.0 0.0 0. SILICON
*
MATERIAL 15.0 30.97376 2.2 27.0 0.0 0. PHOSPHO
*
MATERIAL 16.0 32.066 2.07 28.0 0.0 0. SULFUR
*
MATERIAL 22.0 47.867 4.50 29.0 0.0 0. TITANIUM
*
MATERIAL 24.0 51.9961 7.19 30.0 0.0 0. CHROMIUM
*
MATERIAL 28.0 58.6934 7.19 31.0 0.0 0. NICKEL
*
* Argon material n. 32, normal temperature...
MATERIAL 18.0 39.948 1.400 32.0 0.0 0. ARGON
*
* Boron material n. 33.
MATERIAL 5.0 10.8117 2.460 33.0 0.0 0. BORON
*
* ---------- MAT#34 austenitischer Stahl X6CrNiTi1810 Dichte 7.62 ----------
MATERIAL 0.0 0.0 7.62 34.0 0.0 0. AUSTEN
COMPOUND -1.48E-7 6.0 -5.20E-3 26.0 -2.58E-5 27. AUSTEN
COMPOUND -1.78E-4 28.0 -3.55E-3 29.0 -1.73E-1 30. AUSTEN
COMPOUND -7.03E-1 11.0 -1.14E-1 31.0 0.0 0. AUSTEN
*
* ---------------------------- Water MAT#35 --------------------------------
*
MATERIAL 0.0 0.0 1.0 35.0 0.0 0. WATER
COMPOUND 2.0 3.0 1.0 8.0 0.0 0. WATER
*
* ----------------------------- Air MAT#36 ---------------------------------
*
* Air (mat. 36): 3 components, N (mat. 7), O (mat. 8), gaseous Ar (mat.32)
*
*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
MATERIAL 0.0 0.0 0.001225 36.0 0.0 0. AIR
COMPOUND -0.9256E-3 7.0 -0.2837E-3 8.0 -1.572E-5 32. AIR
*
* ------------------------- Type E Glass MAT#37 ----------------------------
*
* The glass used in the insulation for the coils (epoxy fibre glass) is of type E.
* Glass is the main composition of the insulation (~80%)...so assume 100%
* The chemical compositon of type E glass is:
*
* Table 1.
* Substance Mass%
*----------------------------------------------------
* SiO2 52-56% ...take mid value 54% !!
* Alkaline Oxides 0-2% ... 1%
* CaO 16-25% ... 20.5% !!
* MgO 0-5% ... 2.5% !!
* B2O3 5-10% ... 7.5% !!
* Al2O3 12-16% ... 14% !!
* TiO2 0-0.8% ... 0.4%
* Fe2O3 0.05-0.4% ... 0.225%
* F2 0-1% ... 0.5%
*----------------------------------------------------
*
* Break table 1 into the individual effective mass fractions of the most
* significant elements marked in table 1 with a "!!". We thus obtain...
*
* Si 0.2974
* O 0.4839
* Ca 0.1732
* Mg 0.0180
* B 0.0276
*
* Araldite MY 740 Liquid, unmodified, solvent-free Bisphenol A epoxy resin
* Density at 25 degC 1.15-1.20 g/cm3 (Take mid value 1.175 g/cm3)
*
*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
MATERIAL 0.0 0.0 1.175 37.0 0.0 0. GLASS
COMPOUND -0.2974 14.0 -0.4839 8.0 -0.1732 21. GLASS
COMPOUND -0.0180 9.0 -0.0276 33.0 0.0 0. GLASS
*
*
* #############################################################################
* ####################### ASSIGN MATERIALS TO REGIONS #########################
* #############################################################################
*
* NOTE!!!
* The regions INDEXES specified in WHAT2- correspond to the strings which label
* with 3 characters, those regions defined in the geometry above. Thus take
* GREAT care when commenting out or reintroducing a region that the INDEXES are
* reordered back to a CONTINUOUS numerical order! It is also best advised to
* keep the labels in the geometry above CONSISTENT with the corresponding index.
* E.g. 001 corresponds ro region index 1.0 and so on.
*
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
*
* Black hole to region 001
ASSIGNMA 1.0 1.0
* Air to region 002. Space around targets.
ASSIGNMA 36.0 2.0
*
* Assign C to region 003. Collimator
ASSIGNMA 6.0 3.0
*
* Assign Insulation #37 to region 004.
ASSIGNMA 37.0 4.0
*
*
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
*
* MESH BINNING. PHYSICAL DOSE
* ===========================
*
* Coil insulation Unit Xmax Ymax Zmax
USRBIN 10.0 208.0 42.0 30.00 30.00 1.00 DI40
USRBIN -30.00 -30.00 0.00 1.0 1.0 1.0 &
*
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
RANDOMIZE 1.0
*
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
START 1.0E0
STOP
1 *====================================================================*
* *
* *
* *
* *
* FFFFF L U U K K AAA 222 000 000 666 *
* F L U U K K A A 2 2 0 0 0 0 6 6 *
* F L U U K K A A 2 0 0 0 0 6 *
* FFFF L U U KK AAAAA == 2 0 0 0 0 6666 *
* F L U U K K A A 2 0 0 0 0 6 6 *
* F L U U K K A A 2 0 0 0 0 6 6 *
* F LLLLL UUU K K A A 22222 000 000 666 *
* *
* *
* *
* Version of Fluka2006 for GNU/Linux operating system *
* *
* By : *
* *
* Alfredo Ferrari & Paola Sala INFN Milan & CERN/AB *
* *
* Alfredo.Ferrari@cern.ch Paola.Sala@mi.infn.it *
* *
* *
* Alberto Fasso` SLAC, Stanford *
* *
* Fasso@slac.stanford.edu *
* *
* *
* Johannes Ranft Siegen University *
* *
* Johannes.Ranft@cern.ch *
* *
* *
* *
* This version includes all the features of the Fluka2006 *
* *
* package *
* *
* *
* *
*====================================================================*
1
FLUKA User license:
By downloading and/or using FLUKA, you agree with the
following license/conditions/requests:
Copyright statement and license conditions
Copyright Italian National Institute for Nuclear Physics
(INFN) and European Organization for Nuclear Research (CERN)
("the FLUKA copyright holders"), 1989-2006.
All rights not expressly granted under this license are
reserved.
This software results from work performed by Alberto Fasso`,
Alfredo Ferrari, Johannes Ranft and Paola Sala.
INFN and CERN are the exclusive source of distribution of
the code, bug fixes and documentation of the FLUKA
software. Each official version of FLUKA is identified by a
numbering scheme specifying major and minor releases.
The FLUKA Coordination Committee or its delegates are able
to grant any of the permissions noted in this License
Agreement as requiring a specific consent. Any such consent
may only be granted in writing.
Installation, use, reproduction, display of the FLUKA
software ("FLUKA"), in source and binary forms, are
permitted free of charge on a non-exclusive basis for
internal scientific, non-commercial and non-weapon-related
use by non-profit organizations only. Any exercise of these
rights is subject to the following conditions:
1 Insertion of the FLUKA code, in whole or in part, into
other codes, or its translation into any other computer
language are possible only after obtaining prior written
permission. Modifications of the FLUKA code are permitted
for use by the licensee only, unless authorized in
written.
2 FLUKA is non-transferable, non-sub-licensable and may not
be distributed in any way, without express written
consent, whether in original or modified form. Site-wise
or collaboration-wise conditions can be agreed with the
FLUKA Coordination Committee.
3 Notwithstanding the above, the licensee may modify and
sub-license FLUKA User Routines to third parties in so far
as their purpose is limited to the adaptation of input and
output interfaces of FLUKA and their modification does not
circumvent, replace, add to or modify any of the functions
of FLUKA, or extract specific isolated results from any of
the individual internal physics models embedded within
FLUKA.
4 The licensee shall forthwith license all its modifications
of FLUKA to the FLUKA copyright holders, at no cost and
with no limitation of use. The licensee acknowledges that
the FLUKA copyright holders may insert such modifications
into future releases of FLUKA, subject to appropriate
acknowledgment of the licensee's contribution.
5 Any publication by the licensee with respect to FLUKA or
results obtained from it shall explicitly acknowledge
FLUKA by quoting its set of references and the FLUKA
copyright holders. The licensee shall not without prior
written permission publish documents or results based on a
modified form of FLUKA, unless the modification
exclusively concerns User Routines for the adaptation of
its input and output interfaces which comply with the same
restrictions, as defined in section 3) as those which
apply to sub-licensing. Any publication of documents or
results shall be based only on official FLUKA versions as
obtained from the FLUKA website (http://www.fluka.org) or
from any authorized mirror. Publication here implies any
legal publication to any third party, whether verbal,
electronic, visual, in writing or otherwise.
6 The licensee shall ensure that the FLUKA references,
copyright statement and license conditions are not altered
or removed from FLUKA. Any integration of any portion of
FLUKA, in modified or in unmodified form, into any other
software package must preserve the internal copyright
notices in those portions of FLUKA that have been
employed, and must reproduce such notices within any
additional global notices included along or embedded
within the software into which FLUKA has been
integrated. Any portion of FLUKA so integrated, whether
modified or unmodified shall continue to be subject to
these license conditions.
7 Nothing in this license shall be construed as to grant any
rights in any of the FLUKA versions since 1989. In
addition, users are not permitted to circumvent any
protection in prior distributions of FLUKA that provided
for a preset expiration date of the code
8 Versions or parts of the FLUKA source code, entrusted to
individuals or groups prior to the enactment of the
CERN-INFN Collaboration Agreement, which are listed in
Chapter 5 of Annex 1 of the EP-AB-INFN Scientific
Agreement (19-02-2003), together with the agreed
conditions of use, are subject to this License Agreement
in addition to any other restrictions on the scope of use
that may have been part of the initial use grant.
9 Commercial use of FLUKA, outside the scope of this
license, must be negotiated with the copyright holders
10 DISCLAIMER
THIS SOFTWARE IS PROVIDED BY THE FLUKA COPYRIGHT HOLDERS "AS
IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT
NOT LIMITED TO, IMPLIED WARRANTIES OF MERCHANTABILITY, OF
SATISFACTORY QUALITY, AND FITNESS FOR A PARTICULAR PURPOSE
OR USE ARE DISCLAIMED. THE FLUKA COPYRIGHT HOLDERS AND THE
AUTHORS MAKE NO REPRESENTATION THAT THE SOFTWARE AND
MODIFICATIONS THEREOF, WILL NOT INFRINGE ANY PATENT,
COPYRIGHT, TRADE SECRET OR OTHER PROPRIETARY RIGHT.
11 LIMITATION OF LIABILITY
THE FLUKA COPYRIGHT HOLDERS AND THE AUTHORS SHALL HAVE NO
LIABILITY FOR DIRECT, INDIRECT, SPECIAL, INCIDENTAL,
CONSEQUENTIAL, EXEMPLARY, OR PUNITIVE DAMAGES OF ANY
CHARACTER INCLUDING, WITHOUT LIMITATION, PROCUREMENT OF
SUBSTITUTE GOODS OR SERVICES, LOSS OF USE, DATA OR PROFITS,
OR BUSINESS INTERRUPTION, HOWEVER CAUSED AND ON ANY THEORY
OF CONTRACT, WARRANTY, TORT (INCLUDING NEGLIGENCE), PRODUCT
LIABILITY OR OTHERWISE, ARISING IN ANY WAY OUT OF THE USE OF
THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH
DAMAGES.
REQUESTS MADE BY THE FLUKA AUTHORS TO ALL USERS
All licensees are requested to report as soon as practical
to the Authors any errors or bugs found in any portion of
FLUKA and its accompanying documentation.
All licensees are requested to forward copies of all
comparisons that they may make between FLUKA results and
data or other codes as soon as practical. The Authors agree
to keep any such communications confidential unless
otherwise notified by the contributing user.
The Authors reserve the right of publishing any benchmarking
and/or comparisons of the distinct separate performance of
the individual internal models that are embedded into FLUKA,
whether the comparisons are with data or with other
codes. The Authors would also like to convey a general
willingness to conduct any such benchmarking efforts either
upon request or in collaboration with interested parties. In
case of doubt please contact the Authors.
Users should exclusively download FLUKA from the official
FLUKA website (http://www.fluka.org) or one of the
authorized mirror sites. Users are invited to regularly
update their FLUKA version to profit for improvements and
bug fixes.
Users are invited to use reasonably updated versions of the
code for publications. Publications of results based on
those FLUKA versions that are declared unsupported and
obsolete on the official FLUKA website shall be avoided.
Users should address any request of consent to one member of
the FLUKA Coordinating Committee, which at present is
composed as follows:
Giuseppe Battistoni Giuseppe.Battistoni@mi.infn.it
(chairman)
John Ferguson John.Ferguson@cern.ch
John Harvey John.Harvey@cern.ch
Johannes Ranft Johannes.Ranft@cern.ch
Paola Sala Paola.Sala@mi.infn.it
In accordance with the User License, the use of the FLUKA code
shall be acknowledged explicitly by quoting the following and
only the following set of references:
- A. Ferrari, P.R. Sala, A. Fasso', and J. Ranft,
"FLUKA: a multi-particle transport code",
CERN 2005-10 (2005), INFN/TC_05/11, SLAC-R-773
- A. Fasso', A. Ferrari, S. Roesler, P.R. Sala, G. Battistoni,
F. Cerutti, E. Gadioli, M.V. Garzelli, F. Ballarini,
A. Ottolenghi, A. Empl and J. Ranft,
"The physics models of FLUKA: status and recent developments",
Computing in High Energy and Nuclear Physics 2003 Conference
(CHEP2003), La Jolla, CA, USA, March 24-28, 2003,
(paper MOMT005) eConf C0303241 (2003), arXiv:hep-ph/0306267
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!! DOWNLOAD/USE OF THE FLUKA SOFTWARE IMPLIES FULL !!!!
!!!! ACCEPTANCE OF THE LICENSE AND ASSOCIATED CONDITIONS !!!!
!!!! INCLUDING THE AUTHOR'S REQUESTS !!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1 ****************************************************************************************************
FLUKA2006 Version 3.2 Oct-06 by A. Ferrari DATE: 6/13/ 7 TIME: 15:30:14
****************************************************************************************************
*---------------- ----------------*
*---------------- M. Kirk ----------------*
*---------------- ----------------*
*---------------- 24 April 2007 ----------------*
*---------------- ----------------*
*---------------- ...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8 ----------------*
***** Next control card ***** TITLE 0.000 0.000 0.000 0.000 0.000 0.000
EURISOL DS Beta Beams. Hypothetical geometry. Worser case than expected.
***** Next control card ***** DEFAULTS 0.000 0.000 0.000 0.000 0.000 0.000 NEW-DEFA
*---------------- ...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+... ----------------*
*---------------- *BEAM -1.5 0.0 0.0 10.0 10.0 1.PROTON ----------------*
*---------------- Primary beam is X(A=WHAT2,Z=WHAT1) ----------------*
*---------------- ----------------*
*---------------- ----------------*
*---------------- ...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+... ----------------*
*---------------- Open the file containing the formatted SOURCE data which is read in the ----------------*
*---------------- SOURCE card below. WHAT1 is the unit number of the file. ----------------*
*---------------- A positive unit means formatted file is to be expected. ----------------*
*---------------- The user has to make sure that the given unit number is not alread in use! ----------------*
*---------------- ----------------*
*---------------- ...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8 ----------------*
*---------------- The user written source routine source_cpt.f is called here. ----------------*
*---------------- NOTE: the unit number of the file to read must be the same as that in OPEN. ----------------*
*---------------- ----------------*
*---------------- ----------------*
*---------------- ----------------*
*---------------- ############################################################################# ----------------*
*---------------- ############################## MAKE MATERIALS ############################### ----------------*
*---------------- ############################################################################# ----------------*
*---------------- ----------------*
*---------------- 23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 ----------------*
*---------------- ----------------*
*---------------- ...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8 ----------------*
*---------------- ----------------*
*---------------- The following is just an override of Oxygen, keeping the default settings... ----------------*
*---------------- ----------------*
*---------------- ----------------*
*---------------- ----------------*
*---------------- ----------------*
*---------------- ----------------*
*---------------- ----------------*
*---------------- ----------------*
*---------------- Argon material n. 32, normal temperature... ----------------*
*---------------- ----------------*
*---------------- Boron material n. 33. ----------------*
*---------------- ----------------*
*---------------- ---------- MAT#34 austenitischer Stahl X6CrNiTi1810 Dichte 7.62 ---------- ----------------*
*---------------- ----------------*
*---------------- ---------------------------- Water MAT#35 -------------------------------- ----------------*
*---------------- ----------------*
*---------------- ----------------*
*---------------- ----------------------------- Air MAT#36 --------------------------------- ----------------*
*---------------- ----------------*
*---------------- Air (mat. 36): 3 components, N (mat. 7), O (mat. 8), gaseous Ar (mat.32) ----------------*
*---------------- ----------------*
*---------------- 23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 ----------------*
*---------------- ----------------*
*---------------- ------------------------- Type E Glass MAT#37 ---------------------------- ----------------*
*---------------- ----------------*
*---------------- The glass used in the insulation for the coils (epoxy fibre glass) is of type ----------------*
*---------------- Glass is the main composition of the insulation (~80%)...so assume 100% ----------------*
*---------------- The chemical compositon of type E glass is: ----------------*
*---------------- ----------------*
*---------------- Table 1. ----------------*
*---------------- Substance Mass% ----------------*
*---------------- ---------------------------------------------------- ----------------*
*---------------- SiO2 52-56% ...take mid value 54% !! ----------------*
*---------------- Alkaline Oxides 0-2% ... 1% ----------------*
*---------------- CaO 16-25% ... 20.5% !! ----------------*
*---------------- MgO 0-5% ... 2.5% !! ----------------*
*---------------- B2O3 5-10% ... 7.5% !! ----------------*
*---------------- Al2O3 12-16% ... 14% !! ----------------*
*---------------- TiO2 0-0.8% ... 0.4% ----------------*
*---------------- Fe2O3 0.05-0.4% ... 0.225% ----------------*
*---------------- F2 0-1% ... 0.5% ----------------*
*---------------- ---------------------------------------------------- ----------------*
*---------------- ----------------*
*---------------- Break table 1 into the individual effective mass fractions of the most ----------------*
*---------------- significant elements marked in table 1 with a "!!". We thus obtain... ----------------*
*---------------- ----------------*
*---------------- Si 0.2974 ----------------*
*---------------- O 0.4839 ----------------*
*---------------- Ca 0.1732 ----------------*
*---------------- Mg 0.0180 ----------------*
*---------------- B 0.0276 ----------------*
*---------------- ----------------*
*---------------- Araldite MY 740 Liquid, unmodified, solvent-free Bisphenol A epoxy resin ----------------*
*---------------- Density at 25 degC 1.15-1.20 g/cm3 (Take mid value 1.175 g/cm3) ----------------*
*---------------- ----------------*
*---------------- 23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 ----------------*
*---------------- ----------------*
*---------------- ----------------*
*---------------- ############################################################################# ----------------*
*---------------- ####################### ASSIGN MATERIALS TO REGIONS ######################### ----------------*
*---------------- ############################################################################# ----------------*
*---------------- ----------------*
*---------------- NOTE!!! ----------------*
*---------------- The regions INDEXES specified in WHAT2- correspond to the strings which label ----------------*
*---------------- with 3 characters, those regions defined in the geometry above. Thus take ----------------*
*---------------- GREAT care when commenting out or reintroducing a region that the INDEXES are ----------------*
*---------------- reordered back to a CONTINUOUS numerical order! It is also best advised to ----------------*
*---------------- keep the labels in the geometry above CONSISTENT with the corresponding index. ----------------*
*---------------- E.g. 001 corresponds ro region index 1.0 and so on. ----------------*
*---------------- ----------------*
*---------------- ...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8 ----------------*
*---------------- ----------------*
*---------------- Black hole to region 001 ----------------*
*---------------- Air to region 002. Space around targets. ----------------*
*---------------- ----------------*
*---------------- Assign C to region 003. Collimator ----------------*
*---------------- ----------------*
*---------------- Assign Insulation #37 to region 004. ----------------*
*---------------- ----------------*
*---------------- ----------------*
*---------------- ...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8 ----------------*
*---------------- ----------------*
*---------------- MESH BINNING. PHYSICAL DOSE ----------------*
*---------------- =========================== ----------------*
*---------------- ----------------*
*---------------- Coil insulation Unit Xmax Ymax Zmax ----------------*
*---------------- ----------------*
*---------------- ...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8 ----------------*
*---------------- ----------------*
*---------------- ...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8 ----------------*
***** Next control card ***** PHYSICS 2.000 0.000 0.000 0.000 0.000 0.000 EVAPORAT
**** Density scaling factors start at location 1 and end at 4000 (I*4 addr.) ****
***** Next control card ***** BEAM -1.500 0.000 0.000 10.00 10.00 1.000 HEAVYION
***** Next control card ***** HI-PROPE 3.000 6.000 0.000 0.000 0.000 0.000
***** Next control card ***** BEAMPOS 0.000 0.000 -100.0 0.000 0.000 0.000
***** Next control card ***** OPEN 29.00 0.000 0.000 0.000 0.000 0.000 READONLY
File to be opened: /u/kirk/codes/fluka/fluka_work/BETABE
***** Next control card ***** SOURCE 29.00 1.000 0.000 0.000 0.000 0.000
***** Next control card ***** EVENTYPE 0.000 0.000 2.000 0.000 0.000 0.000 DPMJET
***** Next control card ***** GEOBEGIN 0.000 0.000 0.000 0.000 0.000 0.000 COMBINAT
1
CERN PS dipole with collimator & beampipe
IVOPT = 0 IDBG = 0
Body data
RPP 1 -0.5000000E+07 0.5000000E+07 -0.5000000E+07 0.5000000E+07 -0.5000000E+07 0.5000000E+07 5
RPP 2 -0.1000000E+07 0.1000000E+07 -0.1000000E+07 0.1000000E+07 -0.1000000E+07 0.1000000E+07 15
*
*--------------- Collimator ---------------*
*
RPP 3 -0.3000000E+02 0.3000000E+02 -0.3000000E+02 0.3000000E+02 -0.5000000E+02 0.0000000E+00 25
*
*--------------- Coil insulation ---------------*
*
RPP 4 -0.3000000E+02 0.3000000E+02 -0.3000000E+02 0.3000000E+02 0.0000000E+00 0.1000000E+01 35
END 5 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 45
Number of bodies 4
Length of FPD-Array 50
Input region data
*
*--------------- * ---------------*
*
*
*--------------- * REGIONS ---------------*
*
*
*--------------- * ---------------*
*
*
*--------------- *234x6789x123456x89 123x56789 x234567x9 1234x6789 1x345678x 12345x789 12x stop ---------------*
*
*
*--------------- * Black hole ---------------*
*
001 5 1 -2 0 0 0 0 0 0 0 Z 1
*
*--------------- * Space around targets ---------------*
*
002 5 2 -3 -4 0 0 0 0 0 0 Z 2
*
*--------------- * Collimator ---------------*
*
003 5 3 0 0 0 0 0 0 0 0 Z 3
*
*--------------- * Insulation ---------------*
*
004 5 4 0 0 0 0 0 0 0 0 Z 4
*
*--------------- * ---------------*
*
END 0 0 0 0 0 0 0 0 0 0 Z 5
Number of input regions 4
Number of code zones 4
Length of integer array 101
CODE ZONE INPUT REGION ZONE DATA LOC. NO. OF BODIES REGION NO.
1 1 29 2 1
2 2 38 3 2
3 3 51 1 3
4 4 56 1 4
I IR1(I) IR2(I)
1 1 1
2 2 2
3 3 3
4 4 4
Interpreted body echo
Body n. 1 RPP Body.#1
-5000000. 5000000. -5000000. 5000000. -5000000. 5000000.
Body n. 2 RPP Body.#2
-1000000. 1000000. -1000000. 1000000. -1000000. 1000000.
Body n. 3 RPP Body.#3
-30.00000 30.00000 -30.00000 30.00000 -50.00000 0.000000
Body n. 4 RPP Body.#4
-30.00000 30.00000 -30.00000 30.00000 0.000000 1.000000
Interpreted region echo
Region n. 1 Reg.#1
1 -2
Region n. 2 Reg.#2
2 -3 -4
Region n. 3 Reg.#3
3
Region n. 4 Reg.#4
4
1 OPTION 0 WAS USED IN CALCULATING VOLUMES, FOR 4 REGIONS
3: INPUT VOLUMES, ANYTHING ELSE: VOLUMES = 1.0
VOLUMES (CM**3)
1 REG 1 2 3 4
VOLUME 1.000E+00 1.000E+00 1.000E+00 1.000E+00
*---------------- ----------------*
0NGEOM= 4002, NGLAST= 4276
**** Stars/Energy accumulation arrays start at location 4279 and end at 4310 (I*4 addr.) ****
GEOEND
***** Next control card ***** MATERIAL 1.000 1.008 8.9000E-05 3.000 0.000 0.000 HYDROGEN
***** Next control card ***** MATERIAL 8.000 16.00 1.3300E-03 8.000 0.000 0.000 OXYGEN
***** Next control card ***** MATERIAL 14.00 28.09 2.330 26.00 0.000 0.000 SILICON
***** Next control card ***** MATERIAL 15.00 30.97 2.200 27.00 0.000 0.000 PHOSPHO
***** Next control card ***** MATERIAL 16.00 32.07 2.070 28.00 0.000 0.000 SULFUR
***** Next control card ***** MATERIAL 22.00 47.87 4.500 29.00 0.000 0.000 TITANIUM
***** Next control card ***** MATERIAL 24.00 52.00 7.190 30.00 0.000 0.000 CHROMIUM
***** Next control card ***** MATERIAL 28.00 58.69 7.190 31.00 0.000 0.000 NICKEL
***** Next control card ***** MATERIAL 18.00 39.95 1.400 32.00 0.000 0.000 ARGON
***** Next control card ***** MATERIAL 5.000 10.81 2.460 33.00 0.000 0.000 BORON
***** Next control card ***** MATERIAL 0.000 0.000 7.620 34.00 0.000 0.000 AUSTEN
***** Next control card ***** MATERIAL 0.000 0.000 1.000 35.00 0.000 0.000 WATER
***** Next control card ***** MATERIAL 0.000 0.000 1.2250E-03 36.00 0.000 0.000 AIR
***** Next control card ***** MATERIAL 0.000 0.000 1.175 37.00 0.000 0.000 GLASS
***** Next control card ***** COMPOUND -1.4800E-07 6.000 -5.2000E-03 26.00 -2.5800E-05 27.00 AUSTEN
***** Next control card ***** COMPOUND -1.7800E-04 28.00 -3.5500E-03 29.00 -0.1730 30.00 AUSTEN
***** Next control card ***** COMPOUND -0.7030 11.00 -0.1140 31.00 0.000 0.000 AUSTEN
***** Next control card ***** COMPOUND 2.000 3.000 1.000 8.000 0.000 0.000 WATER
***** Next control card ***** COMPOUND -9.2560E-04 7.000 -2.8370E-04 8.000 -1.5720E-05 32.00 AIR
***** Next control card ***** COMPOUND -0.2974 14.00 -0.4839 8.000 -0.1732 21.00 GLASS
***** Next control card ***** COMPOUND -1.8000E-02 9.000 -2.7600E-02 33.00 0.000 0.000 GLASS
***** Next control card ***** RANDOMIZ 1.000 0.000 0.000 0.000 0.000 0.000
RM64 INITIALIZED: 98765 12345 0 0
***** Next control card ***** ASSIGNMA 1.000 1.000 0.000 0.000 0.000 0.000
***** Next control card ***** ASSIGNMA 36.00 2.000 0.000 0.000 0.000 0.000
***** Next control card ***** ASSIGNMA 6.000 3.000 0.000 0.000 0.000 0.000
***** Next control card ***** ASSIGNMA 37.00 4.000 0.000 0.000 0.000 0.000
***** Next control card ***** USRBIN 10.00 208.0 42.00 30.00 30.00 1.000 DI40
***** Next control card ***** USRBIN -30.00 -30.00 0.000 1.000 1.000 1.000 &
***** Next control card ***** START 1.000 0.000 0.000 0.000 0.000 0.000
Total time used for input reading: 0.126 s
*** Reading evaporation and nuclear data from unit: 14
*** Evaporation: using NNDC (1996) data ***
Starting location in blank common of LVL data: 4312
Last location in blank common of LVL data: 586781
Starting location in blank common of gamma data: 586782
Last location in blank common of gamma data: 689209
Starting location in blank common of beta data: 689210
Last location in blank common of beta data: 722788
Starting location in blank common of GDR data: 722789
Last location in blank common of GDR data: 777843
Starting location in blank common of (g,x) data: 777844
Last location in blank common of (g,x) data: 1080506
**** RIPL2/Ign. self-cons. T=0 N,Z-dep. level density used ****
**** RIPL-2 / Ignyatuk level density en. dep. used ****
**** with Moller, Nix self-cons set of parameters for T=oo ****
**** Original Gilbert/Cameron pairing energy used ****
**** Power exponent for grey particle correlations set at -13.****
**** Maximum Fermi momentum : 0.268371314 GeV/c ****
**** Maximum Fermi energy : 0.0376013778 GeV ****
**** Average Fermi energy : 0.022676846 GeV ****
**** Average binding energy : 0.00768006314 GeV ****
**** Nuclear well depth : 0.04528144 GeV ****
**** Excess mass for 11-B : 0.00866803993 GeV ****
**** Cameron E. m. for 11-B : 0.00792484451 GeV ****
**** Cam.El. E. m. for 11-B : 0.00851833168 GeV ****
**** My.&Sw. E. m. for 11-B : 0.00887639914 GeV ****
**** Atomic mass for 11-B : 10.255106 GeV ****
**** Nuclear mass for 11-B : 10.2525511 GeV ****
**** Excess mass for 12-C : 0. GeV ****
**** Cameron E. m. for 12-C : 0.00151353085 GeV ****
**** Cam.El. E. m. for 12-C : 0.00163579884 GeV ****
**** My.&Sw. E. m. for 12-C : -0.000383192004 GeV ****
**** Atomic mass for 12-C : 11.1779318 GeV ****
**** Nuclear mass for 12-C : 11.1748667 GeV ****
**** Excess mass for 13-C : 0.00312499981 GeV ****
**** Cameron E. m. for 13-C : 0.00287424121 GeV ****
**** Cam.El. E. m. for 13-C : 0.00348056527 GeV ****
**** My.&Sw. E. m. for 13-C : 0.00355792139 GeV ****
**** Atomic mass for 13-C : 12.1125507 GeV ****
**** Nuclear mass for 13-C : 12.1094866 GeV ****
**** Excess mass for 14-N : 0.00286339992 GeV ****
**** Cameron E. m. for 14-N : 0.00337313744 GeV ****
**** Cam.El. E. m. for 14-N : 0.0036975285 GeV ****
**** My.&Sw. E. m. for 14-N : 0.00286558713 GeV ****
**** Atomic mass for 14-N : 13.0437841 GeV ****
**** Nuclear mass for 14-N : 13.0402088 GeV ****
**** Excess mass for 15-O : 0.00285555003 GeV ****
**** Cameron E. m. for 15-O : 0.00240658456 GeV ****
**** Cam.El. E. m. for 15-O : 0.00259008165 GeV ****
**** My.&Sw. E. m. for 15-O : 0.00332844956 GeV ****
**** Atomic mass for 15-O : 13.9752703 GeV ****
**** Nuclear mass for 15-O : 13.9711847 GeV ****
**** Excess mass for 16-O : -0.00473699998 GeV ****
**** Cameron E. m. for 16-O : -0.00331055629 GeV ****
**** Cam.El. E. m. for 16-O : -0.00319385715 GeV ****
**** My.&Sw. E. m. for 16-O : -0.00483116647 GeV ****
**** Atomic mass for 16-O : 14.8991718 GeV ****
**** Nuclear mass for 16-O : 14.8950863 GeV ****
**** Excess mass for 27-Al : -0.0171969198 GeV ****
**** Cameron E. m. for 27-Al : -0.0186259765 GeV ****
**** Cam.El. E. m. for 27-Al : -0.0171393938 GeV ****
**** My.&Sw. E. m. for 27-Al : -0.017145671 GeV ****
**** Atomic mass for 27-Al : 25.1331501 GeV ****
**** Nuclear mass for 27-Al : 25.1265125 GeV ****
**** Excess mass for 28-Si : -0.0214928202 GeV ****
**** Cameron E. m. for 28-Si : -0.0220552329 GeV ****
**** Cam.El. E. m. for 28-Si : -0.0202140864 GeV ****
**** My.&Sw. E. m. for 28-Si : -0.021535717 GeV ****
**** Atomic mass for 28-Si : 26.0603485 GeV ****
**** Nuclear mass for 28-Si : 26.0532017 GeV ****
**** Excess mass for 40-Ca : -0.0348461308 GeV ****
**** Cameron E. m. for 40-Ca : -0.0357677415 GeV ****
**** Cam.El. E. m. for 40-Ca : -0.0336877368 GeV ****
**** My.&Sw. E. m. for 40-Ca : -0.0349160209 GeV ****
**** Atomic mass for 40-Ca : 37.224926 GeV ****
**** Nuclear mass for 40-Ca : 37.2147255 GeV ****
**** Excess mass for 55-Fe : -0.0574751087 GeV ****
**** Cameron E. m. for 55-Fe : -0.0595041849 GeV ****
**** Cam.El. E. m. for 55-Fe : -0.0580860823 GeV ****
**** My.&Sw. E. m. for 55-Fe : -0.0575032495 GeV ****
**** Atomic mass for 55-Fe : 51.1747131 GeV ****
**** Nuclear mass for 55-Fe : 51.1614609 GeV ****
**** Excess mass for 56-Fe : -0.0606013089 GeV ****
**** Cameron E. m. for 56-Fe : -0.0623576604 GeV ****
**** Cam.El. E. m. for 56-Fe : -0.0608849637 GeV ****
**** My.&Sw. E. m. for 56-Fe : -0.0604862086 GeV ****
**** Atomic mass for 56-Fe : 52.1030807 GeV ****
**** Nuclear mass for 56-Fe : 52.0898285 GeV ****
**** Excess mass for 107-Ag: -0.088405259 GeV ****
**** Cameron E. m. for 107-Ag: -0.0891378522 GeV ****
**** Cam.El. E. m. for 107-Ag: -0.0886852369 GeV ****
**** My.&Sw. E. m. for 107-Ag: -0.0882571116 GeV ****
**** Atomic mass for 107-Ag: 99.5814896 GeV ****
**** Nuclear mass for 107-Ag: 99.5576096 GeV ****
**** Excess mass for 132-Xe: -0.0892794058 GeV ****
**** Cameron E. m. for 132-Xe: -0.0898088515 GeV ****
**** Cam.El. E. m. for 132-Xe: -0.0892864987 GeV ****
**** My.&Sw. E. m. for 132-Xe: -0.0894251093 GeV ****
**** Atomic mass for 132-Xe: 122.867973 GeV ****
**** Nuclear mass for 132-Xe: 122.840576 GeV ****
**** Excess mass for 181-Ta: -0.0484412275 GeV ****
**** Cameron E. m. for 181-Ta: -0.0481105074 GeV ****
**** Cam.El. E. m. for 181-Ta: -0.0478131436 GeV ****
**** My.&Sw. E. m. for 181-Ta: -0.0482191741 GeV ****
**** Atomic mass for 181-Ta: 168.552032 GeV ****
**** Nuclear mass for 181-Ta: 168.515137 GeV ****
**** Excess mass for 208-Pb: -0.0217638295 GeV ****
**** Cameron E. m. for 208-Pb: -0.0184126478 GeV ****
**** Cam.El. E. m. for 208-Pb: -0.020333156 GeV ****
**** My.&Sw. E. m. for 208-Pb: -0.0209293775 GeV ****
**** Atomic mass for 208-Pb: 193.72905 GeV ****
**** Nuclear mass for 208-Pb: 193.687683 GeV ****
**** Excess mass for 209-Bi: -0.0182725303 GeV ****
**** Cameron E. m. for 209-Bi: -0.0147107309 GeV ****
**** Cam.El. E. m. for 209-Bi: -0.0171704441 GeV ****
**** My.&Sw. E. m. for 209-Bi: -0.0177380107 GeV ****
**** Atomic mass for 209-Bi: 194.664047 GeV ****
**** Nuclear mass for 209-Bi: 193.687195 GeV ****
**** Excess mass for 235-U : 0.0409132205 GeV ****
**** Cameron E. m. for 235-U : 0.0464000776 GeV ****
**** Cam.El. E. m. for 235-U : 0.0420622788 GeV ****
**** My.&Sw. E. m. for 235-U : 0.0413222089 GeV ****
**** Atomic mass for 235-U : 218.942078 GeV ****
**** Nuclear mass for 235-U : 218.895767 GeV ****
**** Excess mass for 238-U : 0.0473045185 GeV ****
**** Cameron E. m. for 238-U : 0.0524553321 GeV ****
**** Cam.El. E. m. for 238-U : 0.0481762439 GeV ****
**** My.&Sw. E. m. for 238-U : 0.0473943055 GeV ****
**** Atomic mass for 238-U : 221.74295 GeV ****
**** Nuclear mass for 238-U : 221.696655 GeV ****
**** Evaporation from residual nucleus activated ****
**** Deexcitation gamma production activated ****
**** Evaporated "heavies" transport activated ****
**** High Energy fission requested & activated ****
**** Fermi Break Up requested & activated ****
Heavy Ion option initialized inside Evvini 300630
---> requested particle production model DPMJET
implementation under work, apr-04
(evdini) 20000. 5.
0. 0. 0. 0. 0. 0.
---> DPMJET default Glauber parameter mode
0: new (default), 1: old way
change with WHAT(1) 0.
---> DPMJET output on unit=19 by default
change unit with WHAT(2) 0.
---> DPMJET output verbosity level
change with WHAT(3) 0.
---> initializing DPMJET
LOUDPM,IFDPM 19 0
****************************************************************************************************
DPMJET VERSION II.5 (Sept. 1999)
DUAL PARTON MODEL FOR HADRON NUCLEUS COLLISIONS
AND NUCLEUS NUCLEUS COLLISIONS
INCLUDING A FORMATION TIME INTRANUCLEAR CASCADE
Minijets and DTUJET like multiple soft jets
Nuclear evaporation and residual target and
projectile nuclei
****************************************************************************************************
- DPMJET-II.53 interfaced to run within FLUKA 2002 -
> no DPMJET input cards were read
---> DPMJET event generator version: II.53
---> requested particle production model RQMD
orphaned/modified version 2.4
implementation under work, oct-02
---> RQMD flag to enable input default =0
change with WHAT(4) 0.
---> RQMD default standard output unit=none
change with WHAT(5) 0.
---> RQMD default extra output unit=none
change with WHAT(6) 0.
---> initializing RQMD -rqmd-
LOU,IOU,IOE,LIN F -1 -1 F
RQMD init successfully finished
Minimum kinetic energy for DPMJET: 5.000 (GeV/n)
Minimum kinetic energy for rQMD : 0.1001 (GeV/n)
Overall minimum kinetic energy for ion nuclear interactions: 0.1001 (GeV/n)
**** Fluorescence data successfully retrieved from unit 13 ****
**** Subroutine Mulmix: medium n. 36 ****
Number of elements = 3, Density= 1.225000E-03 (g/cm**3)
0 I Z Pa F_i Rho_i
Index Atomic Atomic Proportion Proportion
Number Weight by Number by weight
1 7.00000 14.0067 0.784754 9.255849E-04
2 8.00000 15.9994 0.210573 2.836954E-04
3 18.0000 39.9480 4.673085E-03 1.571974E-05
ZTILDE,AE1O3,BLCCRA= 7.56380E+00 2.51981E+00 9.97355E-03
**** Warning!!! Least square fit for blccre failed to keep max. rel. Blcce err. below 1% ****
**** Max. error is 1.1 %, for beta2 = 0.00358 ****
ZTILDE,AE1O3,BLCCRE= 6.53935E+00 2.51981E+00 1.04506E-02
BLCC,XCC,TFFLU0,XR0FLU= 7.83319E+00 2.65738E-05 8.54719E-01 4.25526E-05
BLCCE,XCCE,TFEMF0,XR0EMF= 8.91162E+00 2.83218E-02 2.24469E+00 9.00128E-02
Particle n.: -6 Ecutm (prim. & sec.) = 3.747 GeV 3.747 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: -5 Ecutm (prim. & sec.) = 2.828 GeV 2.828 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: -4 Ecutm (prim. & sec.) = 2.829 GeV 2.829 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: -3 Ecutm (prim. & sec.) = 1.896 GeV 1.896 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: -2 Ecutm (prim. & sec.) = 5.622 GeV 5.622 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 1 Ecutm (prim. & sec.) = 0.9583 GeV 0.9583 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 2 Ecutm (prim. & sec.) = 0.9583 GeV 0.9583 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 3 Ecutm (prim. & sec.) = 2.0511E-02 GeV 2.0511E-02 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 4 Ecutm (prim. & sec.) = 2.0511E-02 GeV 2.0511E-02 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 10 Ecutm (prim. & sec.) = 0.1257 GeV 0.1257 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 11 Ecutm (prim. & sec.) = 0.1257 GeV 0.1257 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 13 Ecutm (prim. & sec.) = 0.1596 GeV 0.1596 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 14 Ecutm (prim. & sec.) = 0.1596 GeV 0.1596 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 15 Ecutm (prim. & sec.) = 0.5136 GeV 0.5136 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 16 Ecutm (prim. & sec.) = 0.5136 GeV 0.5136 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 20 Ecutm (prim. & sec.) = 1.217 GeV 1.217 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 21 Ecutm (prim. & sec.) = 1.209 GeV 1.209 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 31 Ecutm (prim. & sec.) = 1.209 GeV 1.209 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 33 Ecutm (prim. & sec.) = 1.217 GeV 1.217 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 36 Ecutm (prim. & sec.) = 1.341 GeV 1.341 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 37 Ecutm (prim. & sec.) = 1.341 GeV 1.341 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 38 Ecutm (prim. & sec.) = 1.692 GeV 1.692 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 39 Ecutm (prim. & sec.) = 1.692 GeV 1.692 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 41 Ecutm (prim. & sec.) = 1.797 GeV 1.797 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 42 Ecutm (prim. & sec.) = 1.797 GeV 1.797 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 45 Ecutm (prim. & sec.) = 1.889 GeV 1.889 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 46 Ecutm (prim. & sec.) = 1.889 GeV 1.889 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 49 Ecutm (prim. & sec.) = 1.988 GeV 1.988 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 50 Ecutm (prim. & sec.) = 1.988 GeV 1.988 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 51 Ecutm (prim. & sec.) = 2.305 GeV 2.305 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 52 Ecutm (prim. & sec.) = 2.486 GeV 2.486 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 54 Ecutm (prim. & sec.) = 2.594 GeV 2.594 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 57 Ecutm (prim. & sec.) = 2.305 GeV 2.305 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 58 Ecutm (prim. & sec.) = 2.486 GeV 2.486 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 60 Ecutm (prim. & sec.) = 2.594 GeV 2.594 GeV, Hthnsz = 1.0000E+30 GeV
**** Atomic electron contribution accounted via the Fano correction ****
**** f_Fano = 9.481E+00, Z_Fano = 7.359E+00 ****
**** Subroutine Mulmix: medium n. 6 ****
Number of elements = 1, Density= 2.00000 (g/cm**3)
0 I Z Pa F_i Rho_i
Index Atomic Atomic Proportion Proportion
Number Weight by Number by weight
1 6.00000 12.0107 1.00000 2.00000
ZTILDE,AE1O3,BLCCRA= 6.00000E+00 2.29011E+00 6.40294E-03
ZTILDE,AE1O3,BLCCRE= 5.00000E+00 2.29011E+00 6.40294E-03
BLCC,XCC,TFFLU0,XR0FLU= 1.20908E+04 9.69872E-04 5.53741E-04 3.95301E-05
BLCCE,XCCE,TFEMF0,XR0EMF= 1.41059E+04 1.04758E+00 1.41811E-03 8.36853E-02
Particle n.: -6 Ecutm (prim. & sec.) = 3.747 GeV 3.747 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: -5 Ecutm (prim. & sec.) = 2.828 GeV 2.828 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: -4 Ecutm (prim. & sec.) = 2.829 GeV 2.829 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: -3 Ecutm (prim. & sec.) = 1.896 GeV 1.896 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: -2 Ecutm (prim. & sec.) = 5.622 GeV 5.622 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 1 Ecutm (prim. & sec.) = 0.9583 GeV 0.9583 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 2 Ecutm (prim. & sec.) = 0.9583 GeV 0.9583 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 3 Ecutm (prim. & sec.) = 2.0511E-02 GeV 2.0511E-02 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 4 Ecutm (prim. & sec.) = 2.0511E-02 GeV 2.0511E-02 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 10 Ecutm (prim. & sec.) = 0.1257 GeV 0.1257 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 11 Ecutm (prim. & sec.) = 0.1257 GeV 0.1257 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 13 Ecutm (prim. & sec.) = 0.1596 GeV 0.1596 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 14 Ecutm (prim. & sec.) = 0.1596 GeV 0.1596 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 15 Ecutm (prim. & sec.) = 0.5136 GeV 0.5136 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 16 Ecutm (prim. & sec.) = 0.5136 GeV 0.5136 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 20 Ecutm (prim. & sec.) = 1.217 GeV 1.217 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 21 Ecutm (prim. & sec.) = 1.209 GeV 1.209 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 31 Ecutm (prim. & sec.) = 1.209 GeV 1.209 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 33 Ecutm (prim. & sec.) = 1.217 GeV 1.217 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 36 Ecutm (prim. & sec.) = 1.341 GeV 1.341 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 37 Ecutm (prim. & sec.) = 1.341 GeV 1.341 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 38 Ecutm (prim. & sec.) = 1.692 GeV 1.692 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 39 Ecutm (prim. & sec.) = 1.692 GeV 1.692 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 41 Ecutm (prim. & sec.) = 1.797 GeV 1.797 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 42 Ecutm (prim. & sec.) = 1.797 GeV 1.797 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 45 Ecutm (prim. & sec.) = 1.889 GeV 1.889 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 46 Ecutm (prim. & sec.) = 1.889 GeV 1.889 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 49 Ecutm (prim. & sec.) = 1.988 GeV 1.988 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 50 Ecutm (prim. & sec.) = 1.988 GeV 1.988 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 51 Ecutm (prim. & sec.) = 2.305 GeV 2.305 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 52 Ecutm (prim. & sec.) = 2.486 GeV 2.486 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 54 Ecutm (prim. & sec.) = 2.594 GeV 2.594 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 57 Ecutm (prim. & sec.) = 2.305 GeV 2.305 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 58 Ecutm (prim. & sec.) = 2.486 GeV 2.486 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 60 Ecutm (prim. & sec.) = 2.594 GeV 2.594 GeV, Hthnsz = 1.0000E+30 GeV
**** Atomic electron contribution accounted via the Fano correction ****
**** f_Fano = 9.639E+00, Z_Fano = 6.000E+00 ****
**** Subroutine Mulmix: medium n. 37 ****
Number of elements = 5, Density= 1.17500 (g/cm**3)
0 I Z Pa F_i Rho_i
Index Atomic Atomic Proportion Proportion
Number Weight by Number by weight
1 14.0000 28.0855 0.218562 0.349410
2 8.00000 15.9994 0.624263 0.568526
3 20.0000 40.0780 8.919852E-02 0.203490
4 12.0000 24.3050 1.528596E-02 2.114789E-02
5 5.00000 10.8117 5.269033E-02 3.242676E-02
**** Warning!!! Least square fit for blccra failed to keep max. rel. Blcc err. below 1% ****
**** Max. error is 4.4 %, for beta2 = 0.04066 ****
ZTILDE,AE1O3,BLCCRA= 1.36565E+01 3.42217E+00 3.01316E-02
**** Warning!!! Least square fit for blccre failed to keep max. rel. Blcce err. below 1% ****
**** Max. error is 3.7 %, for beta2 = 0.00358 ****
ZTILDE,AE1O3,BLCCRE= 1.25169E+01 3.42217E+00 3.45718E-02
BLCC,XCC,TFFLU0,XR0FLU= 8.21949E+03 1.04383E-03 8.14549E-04 5.16001E-05
BLCCE,XCCE,TFEMF0,XR0EMF= 8.98595E+03 1.08695E+00 2.22612E-03 1.08790E-01
Particle n.: -6 Ecutm (prim. & sec.) = 3.747 GeV 3.747 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: -5 Ecutm (prim. & sec.) = 2.828 GeV 2.828 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: -4 Ecutm (prim. & sec.) = 2.829 GeV 2.829 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: -3 Ecutm (prim. & sec.) = 1.896 GeV 1.896 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: -2 Ecutm (prim. & sec.) = 5.622 GeV 5.622 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 1 Ecutm (prim. & sec.) = 0.9583 GeV 0.9583 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 2 Ecutm (prim. & sec.) = 0.9583 GeV 0.9583 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 3 Ecutm (prim. & sec.) = 2.0511E-02 GeV 2.0511E-02 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 4 Ecutm (prim. & sec.) = 2.0511E-02 GeV 2.0511E-02 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 10 Ecutm (prim. & sec.) = 0.1257 GeV 0.1257 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 11 Ecutm (prim. & sec.) = 0.1257 GeV 0.1257 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 13 Ecutm (prim. & sec.) = 0.1596 GeV 0.1596 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 14 Ecutm (prim. & sec.) = 0.1596 GeV 0.1596 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 15 Ecutm (prim. & sec.) = 0.5136 GeV 0.5136 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 16 Ecutm (prim. & sec.) = 0.5136 GeV 0.5136 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 20 Ecutm (prim. & sec.) = 1.217 GeV 1.217 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 21 Ecutm (prim. & sec.) = 1.209 GeV 1.209 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 31 Ecutm (prim. & sec.) = 1.209 GeV 1.209 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 33 Ecutm (prim. & sec.) = 1.217 GeV 1.217 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 36 Ecutm (prim. & sec.) = 1.341 GeV 1.341 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 37 Ecutm (prim. & sec.) = 1.341 GeV 1.341 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 38 Ecutm (prim. & sec.) = 1.692 GeV 1.692 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 39 Ecutm (prim. & sec.) = 1.692 GeV 1.692 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 41 Ecutm (prim. & sec.) = 1.797 GeV 1.797 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 42 Ecutm (prim. & sec.) = 1.797 GeV 1.797 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 45 Ecutm (prim. & sec.) = 1.889 GeV 1.889 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 46 Ecutm (prim. & sec.) = 1.889 GeV 1.889 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 49 Ecutm (prim. & sec.) = 1.988 GeV 1.988 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 50 Ecutm (prim. & sec.) = 1.988 GeV 1.988 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 51 Ecutm (prim. & sec.) = 2.305 GeV 2.305 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 52 Ecutm (prim. & sec.) = 2.486 GeV 2.486 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 54 Ecutm (prim. & sec.) = 2.594 GeV 2.594 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 57 Ecutm (prim. & sec.) = 2.305 GeV 2.305 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 58 Ecutm (prim. & sec.) = 2.486 GeV 2.486 GeV, Hthnsz = 1.0000E+30 GeV
Particle n.: 60 Ecutm (prim. & sec.) = 2.594 GeV 2.594 GeV, Hthnsz = 1.0000E+30 GeV
**** Atomic electron contribution accounted via the Fano correction ****
**** f_Fano = 9.003E+00, Z_Fano = 1.186E+01 ****
**** Isotope tabulation data start at location 1080507 and end at 1105352 (I*4 addr.) ****
No radioactive products/decays requested
Flags for applying biasing to prompt and/or decay radiation:
Hadr/muon EM Low en. Neut.
Prompt/Decay Prompt/Decay Prompt/Decay
Inter./decay length: T F T F T F
Leading Particle : T F T F T F
Importance and WW : T F T F T F
EM transport threshold multipliers: prompt decay
1.00E+00 1.00E+00
**** Non analog absorption factors start at location 1105355 and end at 1105362 (I*4 addr.) ****
**** Biased downscattering factors start at location 1105365 and end at 1105372 (I*4 addr.) ****
**** Non analog absorption group limits start at location 1105373 and end at 1105376 (I*4 addr.) ****
**** Biased downscattering group limits start at location 1105377 and end at 1105380 (I*4 addr.) ****
**** Cut-off group limits start at location 1105381 and end at 1105384 (I*4 addr.) ****
**** Lower Weight Window limits start at location 1105387 and end at 1105394 (I*4 addr.) ****
1
***** Neutron Xsec: group limits, average energies, velocities and momenta *****
***** start at location 1105395, end at location 1106018 (I*4 addr.) *****
***** Neutron Xsec: ****
**** Group limits, average energies, velocities and momenta start at location 1105395, end at location 1106018 (I*4 addr.) *****
*** (n,p) proton production activated for Xsec mat. # 3 ***
Group cross sections storage starts at 1110051
Last location used for group xsecs 1301101
**** THE FOLLOWING VALUES ARE FROM TAPE. ****
Panini Xsec
0 Number of primary groups (Ngrpn) 72
Number of primary downscatters (Nds) 72
Number of secondary groups (Ngrpg) 22
Number of secondary downscatters (Ndsg) 22
Number of prim+sec groups (Ingp) 94
Table length (Itbl) 97
Loc of within group (sig gg) (Isgg) 4
Number of media (Nxsmed) 129
Number of coefficients (Ncoef) 6
Number of angles (Nansct) 3
1
*** Fluka to low en. xsec material correspondence: printed atomic densities are meaningless when used in a compound ***
Fluka medium Name Xsec medium atomic density Id. 1 Id. 2 Id. 3
number number ( at/(cm barn) )
1 BLCKHOLE 0 0.0000E+00 0 0 0
2 VACUUM 1000 0.0000E+00 0 0 0
6 CARBON 2 1.0028E-01 6 -2 293
7 NITROGEN 3 0.0000E+00 7 -2 293
8 OXYGEN 4 0.0000E+00 8 16 293
9 MAGNESIU 5 0.0000E+00 12 -2 293
14 SILICON 6 0.0000E+00 14 -2 293
21 CALCIUM 8 0.0000E+00 20 -2 293
32 ARGON 7 0.0000E+00 18 -2 293
33 BORON 1 0.0000E+00 5 -2 293
*** dp/dx tab. generated up to 11.09 GeV/c/n ***
**** Nuclear form factor 'a la Kelner' selected ****
**** Standard Coulomb correction selected ****
**** for charged hadron and muon bremmstrahlung ****
***** dp/dx : material number 36 "AIR " *****
***** Gas: actual (Fluka) pressure : 1.0000E+00 atm. *****
***** Average excitation energy : 8.5667E+01 eV, weighted Z/A : 4.9919E-01 *****
***** Sternheimer density effect parameters: *****
***** X0 = 1.8000, X1 = 4.0000, C = -10.5787, A = 0.2150 m = 3.0000 D0 = 0.0000 *****
***** Restricted energy loss tabulated in 50 intervals *****
***** Delta ray production activated above 1.0000E-03 GeV *****
***** dE/dx fluctuations activated for this medium, level 1 *****
***** (up to 2I discrete levels, up to 2 K-edges) *****
***** Restricted pair production energy loss added *****
***** Exp. pair production activated above 0.0000E+00 GeV *****
***** Restricted bremsstrahlung energy loss added *****
***** Exp. bremsstrahlung activated above 1.0000E-03 GeV *****
***** dp/dx : material number 6 "CARBON " *****
***** Average excitation energy : 7.8000E+01 eV, weighted Z/A : 4.9955E-01 *****
***** Sternheimer density effect parameters: *****
***** X0 = -0.0351, X1 = 2.4860, C = -2.9930, A = 0.2024 m = 3.0040 D0 = 0.1000 *****
***** Restricted energy loss tabulated in 50 intervals *****
***** Delta ray production activated above 1.0000E-03 GeV *****
***** dE/dx fluctuations activated for this medium, level 1 *****
***** (up to 2I discrete levels, up to 2 K-edges) *****
***** Restricted pair production energy loss added *****
***** Exp. pair production activated above 0.0000E+00 GeV *****
***** Restricted bremsstrahlung energy loss added *****
***** Exp. bremsstrahlung activated above 1.0000E-03 GeV *****
***** dp/dx : material number 37 "GLASS " *****
***** Average excitation energy : 1.3560E+02 eV, weighted Z/A : 4.9824E-01 *****
***** Sternheimer density effect parameters: *****
***** X0 = 0.2000, X1 = 3.0000, C = -4.6329, A = 0.1691 m = 3.0000 D0 = 0.0000 *****
***** Restricted energy loss tabulated in 50 intervals *****
***** Delta ray production activated above 1.0000E-03 GeV *****
***** dE/dx fluctuations activated for this medium, level 1 *****
***** (up to 2I discrete levels, up to 2 K-edges) *****
***** Restricted pair production energy loss added *****
***** Exp. pair production activated above 0.0000E+00 GeV *****
***** Restricted bremsstrahlung energy loss added *****
***** Exp. bremsstrahlung activated above 1.0000E-03 GeV *****
***** dp/dx tabulations in blank common start at location 1301105 *****
***** end at location 1306504 (I*4 addresses) *****
*** Range: P_max: 1.77827895 P_min 0.00398107152 ***
***** Range tabulations in blank common start at location 1306507 *****
***** end at location 1306830 (I*4 addresses) *****
***** Sigtab : material number 7 "NITROGEN" *****
***** Muon photonuclear interac. requested *****
***** Xsec tabulations in blank common start at location 1306833 *****
***** end at location 1315004 (I*4 addresses) for this medium *****
***** Sigtab : material number 8 "OXYGEN " *****
***** Muon photonuclear interac. requested *****
***** Xsec tabulations in blank common start at location 1315005 *****
***** end at location 1323176 (I*4 addresses) for this medium *****
***** Sigtab : material number 32 "ARGON " *****
***** Muon photonuclear interac. requested *****
***** Xsec tabulations in blank common start at location 1323177 *****
***** end at location 1331348 (I*4 addresses) for this medium *****
***** Sigtab : material number 36 "AIR " *****
***** Muon photonuclear interac. requested *****
***** Xsec tabulations in blank common start at location 1331349 *****
***** end at location 1339520 (I*4 addresses) for this medium *****
***** Sigtab : material number 6 "CARBON " *****
***** Muon photonuclear interac. requested *****
***** Xsec tabulations in blank common start at location 1339521 *****
***** end at location 1347692 (I*4 addresses) for this medium *****
***** Sigtab : material number 14 "SILICON " *****
***** Muon photonuclear interac. requested *****
***** Xsec tabulations in blank common start at location 1347693 *****
***** end at location 1355864 (I*4 addresses) for this medium *****
***** Sigtab : material number 21 "CALCIUM " *****
***** Muon photonuclear interac. requested *****
***** Xsec tabulations in blank common start at location 1355865 *****
***** end at location 1364036 (I*4 addresses) for this medium *****
***** Sigtab : material number 9 "MAGNESIU" *****
***** Muon photonuclear interac. requested *****
***** Xsec tabulations in blank common start at location 1364037 *****
***** end at location 1372208 (I*4 addresses) for this medium *****
***** Sigtab : material number 33 "BORON " *****
***** Muon photonuclear interac. requested *****
***** Xsec tabulations in blank common start at location 1372209 *****
***** end at location 1380380 (I*4 addresses) for this medium *****
***** Sigtab : material number 37 "GLASS " *****
***** Muon photonuclear interac. requested *****
***** Xsec tabulations in blank common start at location 1380381 *****
***** end at location 1388552 (I*4 addresses) for this medium *****
***** Total,elastic,inelastic,pair,bremss. macroscopic cross sections tabulated in 72 intervals *****
***** Xsec tabulations in blank common start at location 1306833 *****
***** end at location 1388552 (I*4 addresses) *****
***** Form Fact. tabulations in blank common start at location 1388553 *****
***** end at location 1449924 (I*4 addresses) *****
Chem. sym. : N, Z = 7., A = 14.007, PZ = 7.84754E-01 , RHOZ = 9.25585E-04
Chem. sym. : O, Z = 8., A = 15.999, PZ = 2.10573E-01 , RHOZ = 2.83695E-04
Chem. sym. : Ar, Z = 18., A = 39.948, PZ = 4.67309E-03 , RHOZ = 1.57197E-05
Chem. sym. : C, Z = 6., A = 12.011, PZ = 1.00000E+00 , RHOZ = 2.00000E+00
Chem. sym. : Ca, Z = 20., A = 40.078, PZ = 8.91985E-02 , RHOZ = 2.03490E-01
Chem. sym. : Si, Z = 14., A = 28.085, PZ = 2.18562E-01 , RHOZ = 3.49410E-01
Chem. sym. : O, Z = 8., A = 15.999, PZ = 6.24263E-01 , RHOZ = 5.68526E-01
Chem. sym. : Mg, Z = 12., A = 24.305, PZ = 1.52860E-02 , RHOZ = 2.11479E-02
Chem. sym. : B, Z = 5., A = 10.812, PZ = 5.26903E-02 , RHOZ = 3.24268E-02
*** Blank common cells from 1449925 to 1450155 allocated for an EMF R*4 tab. array
*** Blank common cells from 1450156 to 1450455 allocated for an EMF R*4 tab. array
*** Blank common cells from 1450456 to 1450771 allocated for an EMF R*4 tab. array
*** Blank common cells from 1450772 to 1451087 allocated for an EMF R*4 tab. array
*** Blank common cells from 1451088 to 1451308 allocated for an EMF R*4 tab. array
*** Blank common cells from 1451309 to 1451600 allocated for an EMF R*4 tab. array
*** Blank common cells from 1451601 to 1451909 allocated for an EMF R*4 tab. array
*** Blank common cells from 1451910 to 1452218 allocated for an EMF R*4 tab. array
*** Blank common cells from 1452219 to 1452431 allocated for an EMF R*4 tab. array
*** Blank common cells from 1452432 to 1452716 allocated for an EMF R*4 tab. array
*** Blank common cells from 1452717 to 1453018 allocated for an EMF R*4 tab. array
*** Blank common cells from 1453019 to 1453321 allocated for an EMF R*4 tab. array
*** Blank common cells from 1453322 to 1453527 allocated for an EMF R*4 tab. array
*** Blank common cells from 1453528 to 1453806 allocated for an EMF R*4 tab. array
*** Blank common cells from 1453807 to 1454103 allocated for an EMF R*4 tab. array
*** Blank common cells from 1454104 to 1454400 allocated for an EMF R*4 tab. array
*** Blank common cells from 1454401 to 1456160 allocated for an EMF R*4 tab. array
*** Blank common cells from 1456161 to 1457988 allocated for an EMF R*4 tab. array
*** Blank common cells from 1457989 to 1459827 allocated for an EMF R*4 tab. array
*** Blank common cells from 1459828 to 1461666 allocated for an EMF R*4 tab. array
*** Blank common cells from 1461667 to 1463558 allocated for an EMF R*4 tab. array
*** Blank common cells from 1463559 to 1463733 allocated for an EMF R*4 tab. array
*** Blank common cells from 1463734 to 1463971 allocated for an EMF R*4 tab. array
*** Blank common cells from 1463972 to 1464217 allocated for an EMF R*4 tab. array
*** Blank common cells from 1464218 to 1464463 allocated for an EMF R*4 tab. array
*** Blank common cells from 1464464 to 1464761 allocated for an EMF R*4 tab. array
*** Blank common cells from 1464762 to 1465078 allocated for an EMF R*4 tab. array
*** Blank common cells from 1465079 to 1465395 allocated for an EMF R*4 tab. array
*** Blank common cells from 1465396 to 1465556 allocated for an EMF R*4 tab. array
*** Blank common cells from 1465557 to 1465775 allocated for an EMF R*4 tab. array
*** Blank common cells from 1465776 to 1466000 allocated for an EMF R*4 tab. array
*** Blank common cells from 1466001 to 1466225 allocated for an EMF R*4 tab. array
*** Blank common cells from 1466226 to 1466504 allocated for an EMF R*4 tab. array
*** Blank common cells from 1466505 to 1466800 allocated for an EMF R*4 tab. array
*** Blank common cells from 1466801 to 1467097 allocated for an EMF R*4 tab. array
*** Blank common cells from 1467098 to 1467816 allocated for an EMF R*4 tab. array
*** Blank common cells from 1467817 to 1468591 allocated for an EMF R*4 tab. array
*** Blank common cells from 1468592 to 1469371 allocated for an EMF R*4 tab. array
*** Blank common cells from 1469372 to 1470151 allocated for an EMF R*4 tab. array
*** Blank common cells from 1470152 to 1470979 allocated for an EMF R*4 tab. array
*** Blank common cells from 1470980 to 1471818 allocated for an EMF R*4 tab. array
*** Blank common cells from 1471819 to 1472657 allocated for an EMF R*4 tab. array
*** Blank common cells from 1472658 to 1473542 allocated for an EMF R*4 tab. array
Electron Xsecs for EMF medium # 1(AIR ) tabulated in 322 intervals
Min kin. energy 1.000E-03 (GeV)Max kin. energy 11.1 (GeV)
*** Blank common cells from 1478781 to 1479426 allocated for an EMF tab. array
*** Blank common cells from 1479427 to 1480072 allocated for an EMF tab. array
*** Blank common cells from 1480073 to 1480718 allocated for an EMF tab. array
*** Blank common cells from 1480719 to 1481364 allocated for an EMF tab. array
*** Blank common cells from 1481365 to 1482010 allocated for an EMF tab. array
*** Blank common cells from 1482011 to 1482656 allocated for an EMF tab. array
*** Blank common cells from 1482657 to 1483302 allocated for an EMF tab. array
*** Blank common cells from 1483303 to 1483948 allocated for an EMF tab. array
*** Blank common cells from 1483949 to 1484594 allocated for an EMF tab. array
*** Blank common cells from 1484595 to 1485240 allocated for an EMF tab. array
*** Blank common cells from 1485241 to 1485886 allocated for an EMF tab. array
*** Blank common cells from 1485887 to 1486532 allocated for an EMF tab. array
*** Blank common cells from 1486533 to 1487178 allocated for an EMF tab. array
*** Blank common cells from 1487179 to 1488240 allocated for an EMF tab. array
*** Blank common cells from 1488241 to 1489302 allocated for an EMF tab. array
*** Blank common cells from 1489303 to 1490364 allocated for an EMF tab. array
*** Blank common cells from 1490365 to 1491426 allocated for an EMF tab. array
*** Blank common cells from 1491427 to 1492488 allocated for an EMF tab. array
*** Blank common cells from 1492489 to 1493550 allocated for an EMF tab. array
*** Blank common cells from 1493551 to 1494612 allocated for an EMF tab. array
Electron Xsecs for EMF medium # 2(CARBON ) tabulated in 322 intervals
Min kin. energy 1.000E-03 (GeV)Max kin. energy 11.1 (GeV)
*** Blank common cells from 1499849 to 1500494 allocated for an EMF tab. array
*** Blank common cells from 1500495 to 1501140 allocated for an EMF tab. array
*** Blank common cells from 1501141 to 1501786 allocated for an EMF tab. array
*** Blank common cells from 1501787 to 1502432 allocated for an EMF tab. array
*** Blank common cells from 1502433 to 1503078 allocated for an EMF tab. array
*** Blank common cells from 1503079 to 1503724 allocated for an EMF tab. array
*** Blank common cells from 1503725 to 1504370 allocated for an EMF tab. array
*** Blank common cells from 1504371 to 1505016 allocated for an EMF tab. array
*** Blank common cells from 1505017 to 1505662 allocated for an EMF tab. array
*** Blank common cells from 1505663 to 1506308 allocated for an EMF tab. array
*** Blank common cells from 1506309 to 1506954 allocated for an EMF tab. array
*** Blank common cells from 1506955 to 1507600 allocated for an EMF tab. array
*** Blank common cells from 1507601 to 1508246 allocated for an EMF tab. array
*** Blank common cells from 1508247 to 1509308 allocated for an EMF tab. array
*** Blank common cells from 1509309 to 1510370 allocated for an EMF tab. array
*** Blank common cells from 1510371 to 1511432 allocated for an EMF tab. array
*** Blank common cells from 1511433 to 1512494 allocated for an EMF tab. array
*** Blank common cells from 1512495 to 1513556 allocated for an EMF tab. array
*** Blank common cells from 1513557 to 1514618 allocated for an EMF tab. array
*** Blank common cells from 1514619 to 1515680 allocated for an EMF tab. array
Electron Xsecs for EMF medium # 3(GLASS ) tabulated in 322 intervals
Min kin. energy 1.000E-03 (GeV)Max kin. energy 11.1 (GeV)
*** Blank common cells from 1520917 to 1521562 allocated for an EMF tab. array
*** Blank common cells from 1521563 to 1522208 allocated for an EMF tab. array
*** Blank common cells from 1522209 to 1522854 allocated for an EMF tab. array
*** Blank common cells from 1522855 to 1523500 allocated for an EMF tab. array
*** Blank common cells from 1523501 to 1524146 allocated for an EMF tab. array
*** Blank common cells from 1524147 to 1524792 allocated for an EMF tab. array
*** Blank common cells from 1524793 to 1525438 allocated for an EMF tab. array
*** Blank common cells from 1525439 to 1526084 allocated for an EMF tab. array
*** Blank common cells from 1526085 to 1526730 allocated for an EMF tab. array
*** Blank common cells from 1526731 to 1527376 allocated for an EMF tab. array
*** Blank common cells from 1527377 to 1528022 allocated for an EMF tab. array
*** Blank common cells from 1528023 to 1528668 allocated for an EMF tab. array
*** Blank common cells from 1528669 to 1529314 allocated for an EMF tab. array
*** Blank common cells from 1529315 to 1530376 allocated for an EMF tab. array
*** Blank common cells from 1530377 to 1531438 allocated for an EMF tab. array
*** Blank common cells from 1531439 to 1532500 allocated for an EMF tab. array
*** Blank common cells from 1532501 to 1533562 allocated for an EMF tab. array
*** Blank common cells from 1533563 to 1534624 allocated for an EMF tab. array
*** Blank common cells from 1534625 to 1535686 allocated for an EMF tab. array
*** Blank common cells from 1535687 to 1536748 allocated for an EMF tab. array
***** EMF tabulations in blank common start at location 1473543 *****
***** end at location 1536748 (I*4 addresses) *****
*** Blank common cells from 1536749 to 1537394 allocated for an EMF tab. array
*** Blank common cells from 1537395 to 1538040 allocated for an EMF tab. array
*** Blank common cells from 1538041 to 1538686 allocated for an EMF tab. array
*** Blank common cells from 1538687 to 1539332 allocated for an EMF tab. array
*** Blank common cells from 1539333 to 1539978 allocated for an EMF tab. array
*** Blank common cells from 1539979 to 1540624 allocated for an EMF tab. array
1 Quantities/Biasing associated with each media:
AIR
Rho = 1.225000E-03 g/cm**3 Rlc= 29890.6 cm
Ae = 1.51100 MeV Ue = 11086.5 MeV
Ap = 0.333333 MeV Up = 11086.0 MeV
dE/dx fluctuations activated for this medium, level 1
below the threshold for explicit secondary electron production
(up to 2I discrete levels, up to 2 K-edges)
CARBON
Rho = 2.00000 g/cm**3 Rlc= 21.3485 cm
Ae = 1.51100 MeV Ue = 11086.5 MeV
Ap = 0.333333 MeV Up = 11086.0 MeV
dE/dx fluctuations activated for this medium, level 1
below the threshold for explicit secondary electron production
(up to 2I discrete levels, up to 2 K-edges)
GLASS
Rho = 1.17500 g/cm**3 Rlc= 21.4218 cm
Ae = 1.51100 MeV Ue = 11086.5 MeV
Ap = 0.333333 MeV Up = 11086.0 MeV
dE/dx fluctuations activated for this medium, level 1
below the threshold for explicit secondary electron production
(up to 2I discrete levels, up to 2 K-edges)
1 Correspondence of regions and EMF-FLUKA material numbers and names:
Region EMF FLUKA
1 0 VACUUM 1 BLCKHOLE
Ecut = 0.0000E+00 MeV, Pcut = 0.0000E+00 MeV, BIAS = F, Ray. = F, S(q,Z) = T, Pz(q,Z) = F
2 1 AIR 36 AIR
Ecut = 1.5110E+00 MeV, Pcut = 3.3333E-01 MeV, BIAS = F, Ray. = F, S(q,Z) = T, Pz(q,Z) = F
3 2 CARBON 6 CARBON
Ecut = 1.5110E+00 MeV, Pcut = 3.3333E-01 MeV, BIAS = F, Ray. = F, S(q,Z) = T, Pz(q,Z) = F
4 3 GLASS 37 GLASS
Ecut = 1.5110E+00 MeV, Pcut = 3.3333E-01 MeV, BIAS = F, Ray. = F, S(q,Z) = T, Pz(q,Z) = F
Starting location in blank common of binning data: 1540627
Last location in blank common of binning data: 1540631
1
=== Output before the actual run - Particle properties: ===
=== Transportable Fluka particles: ===
Particle Number Mass Mean Life Charge Baryon Discard Decay PDG id
(GeV/c**2) (s) number Flag(=1) Flag
4-HELIUM -6 3.7273803 1.000E+18 2 4 0 1 9999
3-HELIUM -5 2.8083922 1.000E+18 2 3 0 1 9999
TRITON -4 2.8089218 1.000E+18 1 3 0 1 9999
DEUTERON -3 1.8756134 1.000E+18 1 2 0 1 9999
HEAVYION -2 5.6015195 1.000E+18 3 6 0 1 9999
OPTIPHOT -1 0.0000000 1.000E+18 0 0 0 1 9999
RAY 0 0.0000000 0.00 0 0 0 1 9999
PROTON 1 0.9382723 1.000E+18 1 1 0 1 2212
APROTON 2 0.9382723 1.000E+18 -1 -1 0 1 -2212
ELECTRON 3 0.0005110 1.000E+18 -1 0 0 1 11
POSITRON 4 0.0005110 1.000E+18 1 0 0 1 -11
NEUTRIE 5 0.0000000 1.000E+18 0 0 1 1 12
ANEUTRIE 6 0.0000000 1.000E+18 0 0 1 1 -12
PHOTON 7 0.0000000 1.000E+18 0 0 0 1 22
NEUTRON 8 0.9395656 889. 0 1 0 1 2112
ANEUTRON 9 0.9395656 889. 0 -1 0 1 -2112
MUON+ 10 0.1056584 2.197E-06 1 0 0 1 -13
MUON- 11 0.1056584 2.197E-06 -1 0 0 1 13
KAONLONG 12 0.4976710 5.170E-08 0 0 0 1 130
PION+ 13 0.1395699 2.603E-08 1 0 0 1 211
PION- 14 0.1395699 2.603E-08 -1 0 0 1 -211
KAON+ 15 0.4936460 1.237E-08 1 0 0 1 321
KAON- 16 0.4936460 1.237E-08 -1 0 0 1 -321
LAMBDA 17 1.1156300 2.632E-10 0 1 0 1 3122
ALAMBDA 18 1.1156300 2.632E-10 0 -1 0 1 -3122
KAONSHRT 19 0.4976710 8.922E-11 0 0 0 1 310
SIGMA- 20 1.1974300 1.479E-10 -1 1 0 1 3112
SIGMA+ 21 1.1893700 7.990E-11 1 1 0 1 3222
SIGMAZER 22 1.1925500 7.400E-20 0 1 0 1 3212
PIZERO 23 0.1349764 8.400E-17 0 0 0 1 111
KAONZERO 24 0.4976710 0.00 0 0 0 1 311
AKAONZER 25 0.4976710 0.00 0 0 0 1 -311
RESERVED 26 0.1349764 8.400E-17 0 0 0 1 0
NEUTRIM 27 0.0000000 1.000E+18 0 0 1 1 14
ANEUTRIM 28 0.0000000 1.000E+18 0 0 1 1 -14
RESERVED 29 0.5474500 0.00 0 0 0 1 0
RESERVED 30 0.0000000 0.00 0 0 0 1 0
ASIGMA- 31 1.1893700 7.990E-11 -1 -1 0 1 -3222
ASIGMAZE 32 1.1925500 7.400E-20 0 -1 0 1 -3212
ASIGMA+ 33 1.1974300 1.479E-10 1 -1 0 1 -3112
XSIZERO 34 1.3149000 2.900E-10 0 1 0 1 3322
AXSIZERO 35 1.3149000 2.900E-10 0 -1 0 1 -3322
XSI- 36 1.3213200 1.639E-10 -1 1 0 1 3312
AXSI+ 37 1.3213200 1.639E-10 1 -1 0 1 -3312
OMEGA- 38 1.6724300 8.220E-11 -1 1 0 1 3334
AOMEGA+ 39 1.6724300 8.220E-11 1 -1 0 1 -3334
RESERVED 40 0.9395656 889. 0 1 0 1 2112
TAU+ 41 1.7770000 2.910E-13 1 0 0 1 -15
TAU- 42 1.7770000 2.910E-13 -1 0 0 1 15
NEUTRIT 43 0.0000000 1.000E+18 0 0 1 1 16
ANEUTRIT 44 0.0000000 1.000E+18 0 0 1 1 -16
D+ 45 1.8693000 1.057E-12 1 0 0 1 411
D- 46 1.8693000 1.057E-12 -1 0 0 1 -411
D0 47 1.8646000 4.150E-13 0 0 0 1 421
D0BAR 48 1.8646000 4.150E-13 0 0 0 1 -421
DS+ 49 1.9683000 4.900E-13 1 0 0 1 431
DS- 50 1.9683000 4.900E-13 -1 0 0 1 -431
LAMBDAC+ 51 2.2849000 2.060E-13 1 1 0 1 4122
XSIC+ 52 2.4663000 3.300E-13 1 1 0 1 4232
XSIC0 53 2.4718000 9.800E-14 0 1 0 1 4132
XSIPC+ 54 2.5741000 5.000E-20 1 1 0 1 4322
XSIPC0 55 2.5788000 5.000E-20 0 1 0 1 4312
OMEGAC0 56 2.7300000 6.400E-14 0 1 0 1 4332
ALAMBDC- 57 2.2849000 2.060E-13 -1 -1 0 1 -4122
AXSIC- 58 2.4663000 3.300E-13 -1 -1 0 1 -4232
AXSIC0 59 2.4718000 9.800E-14 0 -1 0 1 -4132
AXSIPC- 60 2.5741000 5.000E-20 -1 -1 0 1 -4322
AXSIPC0 61 2.5788000 5.000E-20 0 -1 0 1 -4312
AOMEGAC0 62 2.7300000 6.400E-14 0 -1 0 1 -4332
RESERVED 63 0.0000000 0.00 0 0 0 1 9999
RESERVED 64 0.0000000 0.00 0 0 0 1 9999
1
=== Generalised particles (201-233) (for scoring): ===
Generalised particle Number
ALL-PART 201
ALL-CHAR 202
ALL-NEUT 203
ALL-NEGA 204
ALL-POSI 205
NUCLEONS 206
NUC&PI+- 207
ENERGY 208
PIONS+- 209
BEAMPART 210
EM-ENRGY 211
MUONS 212
E+&E- 213
AP&AN 214
KAONS 215
STRANGE 216
KAONS+- 217
HAD-CHAR 218
FISSIONS 219
HE-FISS 220
LE-FISS 221
NEU-BALA 222
HAD-NEUT 223
KAONS0 224
C-MESONS 225
C-(A)BAR 226
CHARMED 227
BLANK 228
UNB-ENER 229
UNB-EMEN 230
X-MOMENT 231
Y-MOMENT 232
Z-MOMENT 233
ACTIVITY 234
ACTOMASS 235
1
=== Output before the actual run - Beam properties ===
Fluka incident beam properties:
Beam particle: HEAVYION Id: -2 (Fluka) 9999 (PDG) Charge: 3 Baryon n.: 6
Mass: 5.602 (GeV/c^2) Mean life: 1.0000E+18 (s) Weight: 1.000
Average beam momentum : 13.508903 (GeV/c)
Average beam kinetic energy: 9.022689 (GeV)
Momentum deviation at FWHM (gaussian): 0.0000000 (GeV/c)
Beam hit position : 0.00000000 0.00000000 -100.000000 cm
Beam direction cosines: 0.00000000 0.00000000 1.00000000
Beam spot FWHM X-width (Rectangular ): 10.000 cm
Beam spot FWHM Y-width (Rectangular ): 10.000 cm
Beam FWHM angular divergence (Rectangular ): 0.0000 (mrad)
The nominal beam position belongs to region: 2(Reg.#2 ), lattice cell: 0( )
Heavy ion beam requested, properties:
Mass number of the projectile: 6
Charge number of the projectile: 3
Isomer number of the projectile: 0
Laboratory momentum of the projectile : 2.246 (GeV/c/amu)
Laboratory kin. energy of the projectile: 1.500 (GeV/amu)
Laboratory tot. energy of the projectile: 2.431 (GeV/amu)
=== User written source: ===
The use of a user written source routine has been requested. The beam features listed above could be overridden.
=== Particle transport thresholds:
Global cut-off kinetic energy for particle transport: 1.000E-02 GeV
The cut-off kinetic energy is superseded by individual particle thresholds if set
Cut-off kinetic energy for 4-HELIUM transport: 1.000E-02 GeV
Cut-off kinetic energy for 3-HELIUM transport: 1.000E-02 GeV
Cut-off kinetic energy for TRITON transport: 1.000E-02 GeV
Cut-off kinetic energy for DEUTERON transport: 1.000E-02 GeV
Cut-off kinetic energy for PROTON transport: 1.000E-02 GeV
Cut-off kinetic energy for APROTON transport: 1.000E-02 GeV
Cut-off kinetic energy for ELECTRON transport defined in the Emfcut card
Cut-off kinetic energy for POSITRON transport defined in the Emfcut card
Cut-off kinetic energy for NEUTRIE transport: 0.000E+00 GeV
Cut-off kinetic energy for ANEUTRIE transport: 0.000E+00 GeV
Cut-off kinetic energy for PHOTON transport defined in the Emfcut card
Cut-off kinetic energy for NEUTRON transport: 1.960E-02 GeV
Cut-off kinetic energy for ANEUTRON transport: 1.000E-05 GeV
Cut-off kinetic energy for MUON+ transport: 1.000E-02 GeV
Cut-off kinetic energy for MUON- transport: 1.000E-02 GeV
Cut-off kinetic energy for KAONLONG transport: 1.000E-02 GeV
Cut-off kinetic energy for PION+ transport: 1.000E-02 GeV
Cut-off kinetic energy for PION- transport: 1.000E-02 GeV
Cut-off kinetic energy for KAON+ transport: 1.000E-02 GeV
Cut-off kinetic energy for KAON- transport: 1.000E-02 GeV
Cut-off kinetic energy for LAMBDA transport: 1.000E-02 GeV
Cut-off kinetic energy for ALAMBDA transport: 1.000E-02 GeV
Cut-off kinetic energy for KAONSHRT transport: 1.000E-02 GeV
Cut-off kinetic energy for SIGMA- transport: 1.000E-02 GeV
Cut-off kinetic energy for SIGMA+ transport: 1.000E-02 GeV
Cut-off kinetic energy for SIGMAZER transport: 1.000E-02 GeV
Cut-off kinetic energy for PIZERO transport: 1.000E-02 GeV
Cut-off kinetic energy for KAONZERO transport: 1.000E-02 GeV
Cut-off kinetic energy for AKAONZER transport: 1.000E-02 GeV
Cut-off kinetic energy for RESERVED transport: 1.000E-02 GeV
Cut-off kinetic energy for NEUTRIM transport: 0.000E+00 GeV
Cut-off kinetic energy for ANEUTRIM transport: 0.000E+00 GeV
Cut-off kinetic energy for RESERVED transport: 1.000E-02 GeV
Cut-off kinetic energy for RESERVED transport: 1.000E-02 GeV
Cut-off kinetic energy for ASIGMA- transport: 1.000E-02 GeV
Cut-off kinetic energy for ASIGMAZE transport: 1.000E-02 GeV
Cut-off kinetic energy for ASIGMA+ transport: 1.000E-02 GeV
Cut-off kinetic energy for XSIZERO transport: 1.000E-02 GeV
Cut-off kinetic energy for AXSIZERO transport: 1.000E-02 GeV
Cut-off kinetic energy for XSI- transport: 1.000E-02 GeV
Cut-off kinetic energy for AXSI+ transport: 1.000E-02 GeV
Cut-off kinetic energy for OMEGA- transport: 1.000E-02 GeV
Cut-off kinetic energy for AOMEGA+ transport: 1.000E-02 GeV
Cut-off kinetic energy for RESERVED transport: 1.000E-02 GeV
Cut-off kinetic energy for TAU+ transport: 1.000E-02 GeV
Cut-off kinetic energy for TAU- transport: 1.000E-02 GeV
Cut-off kinetic energy for NEUTRIT transport: 0.000E+00 GeV
Cut-off kinetic energy for ANEUTRIT transport: 0.000E+00 GeV
Cut-off kinetic energy for D+ transport: 1.000E-02 GeV
Cut-off kinetic energy for D- transport: 1.000E-02 GeV
Cut-off kinetic energy for D0 transport: 1.000E-02 GeV
Cut-off kinetic energy for D0BAR transport: 1.000E-02 GeV
Cut-off kinetic energy for DS+ transport: 1.000E-02 GeV
Cut-off kinetic energy for DS- transport: 1.000E-02 GeV
Cut-off kinetic energy for LAMBDAC+ transport: 1.000E-02 GeV
Cut-off kinetic energy for XSIC+ transport: 1.000E-02 GeV
Cut-off kinetic energy for XSIC0 transport: 1.000E-02 GeV
Cut-off kinetic energy for XSIPC+ transport: 1.000E-02 GeV
Cut-off kinetic energy for XSIPC0 transport: 1.000E-02 GeV
Cut-off kinetic energy for OMEGAC0 transport: 1.000E-02 GeV
Cut-off kinetic energy for ALAMBDC- transport: 1.000E-02 GeV
Cut-off kinetic energy for AXSIC- transport: 1.000E-02 GeV
Cut-off kinetic energy for AXSIC0 transport: 1.000E-02 GeV
Cut-off kinetic energy for AXSIPC- transport: 1.000E-02 GeV
Cut-off kinetic energy for AXSIPC0 transport: 1.000E-02 GeV
Cut-off kinetic energy for AOMEGAC0 transport: 1.000E-02 GeV
Cut-off kinetic energy for RESERVED transport: 1.000E-02 GeV
Cut-off kinetic energy for RESERVED transport: 1.000E-02 GeV
"Heavy" residual transport activated
=== Termination conditions: ===
Minimum cpu-time reserved for output: 80.00 sec
Maximum number of beam particles to be followed: 1
Maximum number of stars to be generated: infinite
=== Multiple Coulomb scattering: ===
Moliere Coulomb scattering for primaries: T
Moliere Coulomb scattering for secondaries: T
Hadrons/muons:
Flag for MCS check with boundary normals: F
Flag for Coulomb single scattering(s) at boundaries: F
(# of Coulomb single scattering(s) at boundaries: 1)
Flag for single scatterings below min. (Moliere) energy: F
=== Electromagnetic Showers: ===
EM showers are treated by the EMF (A.Fasso`,A.Ferrari,P.R.Sala) code
Electrons/positrons:
Flag for MCS check with boundary normals: F
Flag for Coulomb single scattering(s) at boundaries: F
(# of Coulomb single scattering(s) at boundaries: 1)
Flag for single scatterings below min. (Moliere) energy: F
1
Region Particle importances RR factor Cut off N.A. abs. N.A. abs. Bias. Dow. Bias. Dow.
number Fluka part. EM part. Low en. n. group group factor group factor
1 1.0000 1.0000 1.0000 73 72 9.500E-01 0 1.500E+00
2 1.0000 1.0000 1.0000 73 72 9.500E-01 0 1.500E+00
3 1.0000 1.0000 1.0000 73 72 9.500E-01 0 1.500E+00
4 1.0000 1.0000 1.0000 73 72 9.500E-01 0 1.500E+00
******* "usrbin" option:
Cartesian binning n. 1 " DI40 " , generalized particle n. 208
X coordinate: from -3.0000E+01 to 3.0000E+01 cm, 1 bins ( 6.0000E+01 cm wide)
Y coordinate: from -3.0000E+01 to 3.0000E+01 cm, 1 bins ( 6.0000E+01 cm wide)
Z coordinate: from 0.0000E+00 to 1.0000E+00 cm, 1 bins ( 1.0000E+00 cm wide)
data will be printed on unit 42 (unformatted if < 0)
accurate deposition along the tracks requested
normalized (per unit volume) data will be printed at the end of the run
******* "USRBDX" option:
No user bdrx defined
******* "USRTRACK" option:
No user track-length estimator defined
******* "USRCOLL" option:
No user collision density estimator defined
******* "Usryield" option:
No user yield estimator defined
******* "RESNUCLEI" option:
No residual nuclei estimator defined
******* "Detect" option:
No detector defined
=== Material compositions: ===
Material Atomic Atomic Density Inelastic Elastic Radiation Inelastic
Number&Name Number Weight Scattering Scattering Length Scattering
Length for Length for Length for
HEAVYION at HEAVYION at neutrons at
Beam energy Beam energy Threshold
Momentum
g/cm**3 cm cm cm cm
1 BLCKHOLE 0.000 0.000 0.000 0.1000E+31 0.1000E+31 0.1000E+31 0.1000E+31
2 VACUUM 0.000 0.000 0.000 0.1000E+31 0.1000E+31 0.1000E+31 0.1000E+31
3 HYDROGEN 1.000 1.008 0.8900E-04 0.1232E+06 0.1124E+35 0.7083E+06 0.6109E+09
4 HELIUM 2.000 4.003 0.1660E-03 0.1095E+06 0.6024E+34 0.5682E+06 0.6024E+34
5 BERYLLIU 4.000 9.012 1.848 12.96 0.5411E+30 35.28 17.02
6 CARBON 6.000 12.01 2.000 14.69 0.5000E+30 21.35 18.09
7 NITROGEN 7.000 14.01 0.1170E-02 0.2710E+05 0.8547E+33 0.3247E+05 0.3283E+05
8 OXYGEN 8.000 16.00 0.1330E-02 0.2517E+05 0.7519E+33 0.2574E+05 0.2979E+05
9 MAGNESIU 12.00 24.30 1.740 22.60 0.5747E+30 14.39 26.35
10 ALUMINUM 13.00 26.98 2.699 15.47 0.3705E+30 8.896 17.08
11 IRON 26.00 55.84 7.874 7.817 0.1270E+30 1.757 9.077
12 COPPER 29.00 63.55 8.960 7.343 0.1116E+30 1.436 8.381
13 SILVER 47.00 107.9 10.50 8.253 0.9524E+29 0.8543 10.58
14 SILICON 14.00 28.09 2.329 18.43 0.4294E+30 9.370 20.21
15 GOLD 79.00 197.0 19.32 6.090 0.5176E+29 0.3344 7.014
16 MERCURY 80.00 200.6 13.55 8.764 0.7382E+29 0.4752 9.847
17 LEAD 82.00 207.2 11.35 10.63 0.8811E+29 0.5612 12.05
18 TANTALUM 73.00 180.9 16.65 6.766 0.6005E+29 0.4094 7.382
19 SODIUM 11.00 22.99 0.9710 39.36 0.1030E+31 28.56 47.49
20 ARGON 18.00 39.95 0.1660E-02 0.3081E+05 0.6024E+33 0.1178E+05 0.3734E+05
21 CALCIUM 20.00 40.08 1.550 33.56 0.6452E+30 10.42 36.21
22 TIN 50.00 118.7 7.310 12.42 0.1368E+30 1.206 14.35
23 TUNGSTEN 74.00 183.8 19.30 5.885 0.5181E+29 0.3504 6.169
24 TITANIUM 22.00 47.87 4.540 12.46 0.2203E+30 3.560 15.05
25 NICKEL 28.00 58.69 8.902 7.126 0.1123E+30 1.424 8.396
26 SILICON 14.00 28.09 2.330 18.42 0.4292E+30 9.366 20.20
27 PHOSPHO 15.00 30.97 2.200 20.48 0.4545E+30 9.639 24.24
28 SULFUR 16.00 32.07 2.070 22.28 0.4831E+30 9.419 27.23
29 TITANIUM 22.00 47.87 4.500 12.57 0.2222E+30 3.592 15.18
30 CHROMIUM 24.00 52.00 7.190 8.229 0.1391E+30 2.079 9.677
31 NICKEL 28.00 58.69 7.190 8.823 0.1391E+30 1.763 10.40
32 ARGON 18.00 39.95 1.400 36.54 0.7143E+30 13.96 44.28
33 BORON 5.000 10.81 2.460 11.09 0.4065E+30 21.42 14.50
34 AUSTEN 25.70 55.12 7.620 8.028 0.1312E+30 1.825 9.355
Material Number Atom content Partial Densities
CARBON 6 0.67991E-06 0.11289E-05
SILICON 26 0.10216E-01 0.39665E-01
PHOSPHO 27 0.45960E-04 0.19680E-03
SULFUR 28 0.30629E-03 0.13578E-02
TITANIUM 29 0.40921E-02 0.27079E-01
CHROMIUM 30 0.18358 1.3196
IRON 11 0.69459 5.3625
NICKEL 31 0.10717 0.86959
35 WATER 3.333 6.005 1.000 27.22 0.1000E+31 36.08 44.61
Material Number Atom content Partial Densities
HYDROGEN 3 0.66667 0.11189
OXYGEN 8 0.33333 0.88811
36 AIR 7.262 14.55 0.1225E-02 0.2633E+05 0.8163E+33 0.2989E+05 0.3174E+05
Material Number Atom content Partial Densities
NITROGEN 7 0.78475 0.92558E-03
OXYGEN 8 0.21057 0.28370E-03
ARGON 32 0.46731E-02 0.15720E-04
37 GLASS 10.28 20.64 1.175 32.51 0.8511E+30 21.42 37.36
Material Number Atom content Partial Densities
SILICON 14 0.21856 0.34941
OXYGEN 8 0.62426 0.56853
CALCIUM 21 0.89199E-01 0.20349
MAGNESIU 9 0.15286E-01 0.21148E-01
BORON 33 0.52690E-01 0.32427E-01
=== Regions: materials and fields ===
Region N. and Name Material N. and Name Magn./El. Field (on/off)
Minimum and Maximum step size (cm)
1 Reg.#1 1 BLCKHOLE OFF 0.00000E+00 9.99852E+04
2 Reg.#2 36 AIR OFF 0.00000E+00 9.99852E+04
3 Reg.#3 6 CARBON OFF 0.00000E+00 9.99852E+04
4 Reg.#4 37 GLASS OFF 0.00000E+00 9.99852E+04
=== End of the output associated with the input ===
Total time used for initialization: 5.22 s
1NUMBER OF BEAM NUMBER OF BEAM APPROXIMATE NUMBER AVERAGE TIME USED TIME LEFT (RESERVED NUMBER OF STARS
PARTICLES HANDLED PARTICLES LEFT OF BEAM PARTICLES BY A BEAM PARTICLE 80.0 SECONDS CREATED
THAT CAN STILL BE FOR PRINTOUT)
HANDLED
NEXT SEEDS: 0 0 0 0 0 0 181CD 3039 0 0
USER-ROUTINE SOURCE last modified Markus Kirk 11.06.2007
IJBEAM = -2 AND LRDBEA = F
read data from file SDSOU from directory...
0. 8.28639E-06
0.0151515 1.05402E-05
0.030303 1.3358E-05
0.0454545 1.6867E-05
0.0606061 2.12196E-05
0.0757576 2.65976E-05
0.0909091 3.32164E-05
0.106061 4.13302E-05
0.121212 5.12374E-05
0.136364 6.32865E-05
0.151515 7.78825E-05
0.166667 9.54935E-05
0.181818 0.000116657
0.19697 0.000141989
0.212121 0.000172189
0.227273 0.000208045
0.242424 0.000250447
0.257576 0.000300385
0.272727 0.000358959
0.287879 0.000427383
0.30303 0.000506984
0.318182 0.000599206
0.333333 0.000705607
0.348485 0.000827855
0.363636 0.000967722
0.378788 0.00112707
0.393939 0.00130785
0.409091 0.00151206
0.424242 0.00174174
0.439394 0.00199896
0.454545 0.00228575
0.469697 0.0026041
0.484848 0.00295592
0.5 0.00334297
0.515152 0.00376683
0.530303 0.00422888
0.545455 0.00473019
0.560606 0.00527154
0.575758 0.0058533
0.590909 0.00647543
0.606061 0.00713743
0.621212 0.00783826
0.636364 0.00857634
0.651515 0.00934952
0.666667 0.010155
0.681818 0.0109895
0.69697 0.011849
0.712121 0.0127288
0.727273 0.0136238
0.742424 0.0145283
0.757576 0.015436
0.772727 0.0163403
0.787879 0.0172342
0.80303 0.0181103
0.818182 0.0189612
0.833333 0.0197793
0.848485 0.0205571
0.863636 0.021287
0.878788 0.0219621
0.893939 0.0225756
0.909091 0.023121
0.924242 0.0235929
0.939394 0.0239861
0.954545 0.0242964
0.969697 0.0245206
0.984848 0.024656
1. 0.0247014
1.01515 0.024656
1.0303 0.0245206
1.04545 0.0242964
1.06061 0.0239861
1.07576 0.0235929
1.09091 0.023121
1.10606 0.0225756
1.12121 0.0219621
1.13636 0.021287
1.15152 0.0205571
1.16667 0.0197793
1.18182 0.0189612
1.19697 0.0181103
1.21212 0.0172342
1.22727 0.0163403
1.24242 0.015436
1.25758 0.0145283
1.27273 0.0136238
1.28788 0.0127288
1.30303 0.011849
1.31818 0.0109895
1.33333 0.010155
1.34848 0.00934952
1.36364 0.00857634
1.37879 0.00783826
1.39394 0.00713743
1.40909 0.00647544
1.42424 0.0058533
1.43939 0.00527154
1.45455 0.00473019
1.4697 0.00422888
1.48485 0.00376683
1.5 0.00334297
...done.
MAXEPB = 0.0247014
WHASOU(2) = 1.
ENERG = 0.943765222
WEIGHT = 0.999999983
WTFLK (NPFLKA) = 1.
WEIPRI = 1.
IONID = IJHION = 300630
This archive was generated by hypermail 2.1.6 : Thu Jun 14 2007 - 10:19:46 CEST