From: Lutz Moritz (moritz@triumf.ca)
Date: Thu Oct 07 2004 - 18:41:04 CEST
Now that the summer holidays are well out of the way, I wonder if
someone could help me with the question I asked last April.
<http://www.fluka.org/web_archive/earchive/new-fluka-discuss/0154.html>
I have a set up a very simple problem: a cylindrical (almost) target
made of UC2 bombarded by 500 MeV protons (input /file isac_uc2.inp/
attached). We are trying to estimate the target radioactivity inventory
and also to optimize target and beam parameters for this and more
complicated target configurations. I use RESNUC to obtain the residual
nuclei distributions and these give reasonable looking results. BUT, in
trying to score the density of fissions in the target (total fission,
low-energy fission and high-energy fission, particle types 219, 221 and
220) using USRBIN, I get only /zeros/, despite the fact that in the
normal output listing of the statistics of the run there are totals for
both high and low energy fission:
7.1909E-02 (100.%) HIGH ENERGY FISSIONS GENERATED 'PER BEAM PARTICLE
OUT OF WHICH
6.8742E-02 (95.6%) PROTON
0.0000E+00 ( 0.0%) APROTON
0.0000E+00 ( 0.0%) ELECTRON
0.0000E+00 ( 0.0%) POSITRON
0.0000E+00 ( 0.0%) NEUTRIE
0.0000E+00 ( 0.0%) ANEUTRIE
0.0000E+00 ( 0.0%) PHOTON
2.8202E-03 ( 3.9%) NEUTRON
0.0000E+00 ( 0.0%) ANEUTRON
2.2848E-01 (100.%) SECONDARIES CREATED PER BEAM PARTICLE IN LOW
ENERGY NEUTRONS INTERACTIONS OUT OF WHICH
1.3210E-01 (57.8%) NEUTRON
1.5782E-02 ( 6.9%) FISSION NEUTRON
0.0000E+00 ( 0.0%) PROTON
8.0592E-02 (35.3%) PHOTON
Why can't I score the fission density with USRBIN?
Thanks,
Lutz
-- Lutz Moritz Manager, EH&S TRIUMF 4004 Wesbrook Mall Vancouver, B.C. V6T 2A3 Canada ------------------------------------------------------------------------ Tel: 604 222 7372 E-mail: lutz.moritz@triumf.ca <mailto:lutz.moritz@triumf.ca> Website: http://www.triumf.ca/EHS
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
* This is the setup for an ISAC UC2 target
* Lutz Moritz, June 2004
TITLE
ISAC UC2 direct target
* *
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
******************* GEOMETRY DEFINITION *********************************
GEOBEGIN COMBINAT
0 0 ISAC UC2 Target
*
* X-axis up, Y-axis right, Z-axis into the page.
* Origin is at the centre of the target
*
* Body Definitions
*
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
*
* External void
*
RPP 1 -999999.9 999999.9 -999999.9 999999.9 -999999.9 999999.9
*
* Target cylinder
*
ZCC 2 0.0 0.0 0.9
*
* Surrounding vacuum
*
RCC 3 0.0 0.0 -5.0 0.0 0.0 23.0
15.0
*
* Scoring planes along target
*
XYP 4 0.0
XYP 5 4.0
*
* Plane for cut-out for ion streaming
*
YZP 6 0.5
*
END
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
*
* Region Definitions
*
**REG.....OR....1OR....2OR....3OR....4OR....5OR....6OR....7OR....8OR....9
*
* 1) Black hole around geometry
*
BH1 +1 -3
*
* 2) Target
*
TA2 +2 +6 +5 -4
*
* 3) Surrounding vacuum
*
VC3 +3 +4
VC4 +3 -2 +5 -4
VC5 +3 -5
VC6 +2 -6 +5 -4
*
END
**REG.....OR....1OR....2OR....3OR....4OR....5OR....6OR....7OR....8OR....9
*
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
*
GEOEND
*
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
*GEOEND 400.0 400.0 1000.0 -400.0 -400.0 -500.0DEBUG
*GEOEND 80.0 80.0 200.0 &
*STOP
*
******************* MATERIAL DEFINITIONS ********************************
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
*
MATERIAL 6.0 12.011 2.0 3.0 0.0 0.0CARBON
MATERIAL 92.0 238.000 19.07 4.0 0.0 0.0URANIUM
*
* Note: density of UC2 assumed to be 45% of actual because target is made of
* stack of thin disks
*
MATERIAL 0.0 0.0 5.076 5.0 0.0 0.0UC2
COMPOUND 1.0 4.0 2.0 3.0 0.0 0.0UC2
*
LOW-MAT 3.0 6.0 -2.0 293.0
LOW-MAT 4.0 92.0 238.0 293.0
*
************* Material to region correspondence *************************
*
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
ASSIGNMA 1.0 1.0 BLACKHOLE
ASSIGNMA 5.0 2.0 UC2
ASSIGNMA 2.0 3.0 6.0 VACUUM
**
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
******************* BEAM SPECIFICATION **********************************
*
BEAM -0.500 0.0 0.0 -0.25 -0.25 0.0PROTON
*
BEAMPOS 0.0 0.0 -0.1 0.0 0.0 0.0
*
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
******************* EVENT GENERATOR, PARTICLES AND FIELDS ***************
*
EVENTYPE 0.0 0.0 1.0 0.0 0.0 0.0EVAP
PHYSICS 2.0 0.0 0.0 0.0 0.0 0.0EVAPORAT
*
* Activate to discard neutrinos and antineutrinos,
DISCARD 5.0 6.0
*
* Activate to discard pi-zeros, electrons and photons
DISCARD 23.0 3.0 4.0 7.0
*
* Activate to discard muons
DISCARD 10.0 11.0
*
* Set threshold to 19.6 MeV for neutrons
*
PART-THR -0.0196 8.0
*
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
******************* LOW ENERGY NEUTRONS *********************************
*
* Switch on low-energy neutron transport
*
LOW-NEUT 72.0 22.0 0.0196 1.0
*
*
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
******************* OUTPUT AND SCORING **********************************
*
EMF EMF-OFF
*
RESNUCLE 3.0 72.0 0.0 0.0 2.0 1.0UC2
*
SCORE 208.0 201.0
*
* Score beam particles
*
USRBIN 11.0 210.0 73.0 5.0 0.0 5.0BEAM
USRBIN 0.0 0.0 -1.0 20.0 1.0 24.0&
*
* Score fission density
*
USRBIN 11.0 219.0 74.0 5.0 0.0 5.0ALLFISS
USRBIN 0.0 0.0 -1.0 20.0 1.0 24.0&
*
* Score high-energy fission density
*
USRBIN 11.0 220.0 75.0 5.0 0.0 5.0HIFISS
USRBIN 0.0 0.0 -1.0 20.0 1.0 24.0&
*
* Score low-energy fission density
*
USRBIN 11.0 221.0 76.0 5.0 0.0 5.0LOFISS
USRBIN 0.0 0.0 -1.0 20.0 1.0 24.0&
*
*
*************** START SIMULATION ************************************
* Read seeds from unit 1
RANDOMIZE 1.0
* Start the calculation
* /No. protons/No. stars/Print time/Inc.beam//default interactive tlim./
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
*STOP
START 1.0E+07
STOP
This archive was generated by hypermail 2.1.6 : Thu Oct 07 2004 - 23:02:10 CEST