From: Fan, Ph.D., Jiajin (Jiajin.Fan@fccc.edu)
Date: Mon Nov 27 2006 - 04:25:27 CET
Hi All,
I tried to calculate the neutron production per incident proton in a
10x10x10 cm3 cubic tungsten brick with unit of neutrons/incident proton. But
the result is much smaller than I expected. I don't know that's the problem
of input file or something wrong in post processing. I use USYUSW to
analysis result. Attached is my input file and result, please help to check.
Thank you.
TITLE
Neutron generation from primary collimator
*---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+-
---8
DEFAULTS
HADROTHE
*---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+-
---8
* 0.0 x 0.0 cm2 field size, maximum energy 500 MeV, using a spectrum of
300MeV
BEAM -0.10 0.0 0.0 0.0 0.0 1.0PROTON
BEAMPOS 0.0 0.0 -49.5 0.0 0.0
POSITIVE
* SOURCE
*---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+-
---8
GEOBEGIN
COMBINAT
Disc
SPH 1 0.0 0.0 0.0 +800.0
SPH 2 0.0 0.0 0.0 +600.0
* Tungsten collimator X: -5.0~5.0; Y: -5.0~5.0; Z: -49 - -39
RPP 3 -5.0 5.0 -5.0 5.0 -49.0 -39.0
END
*---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+-
---8
*-reg-----or-----or-----or-----or-----or-----or-----or-----or-----or-----
* Black hole at the beginning
1 3 1 -2
* Vacuum
2 3 2 -3
* Tungsten collimator
3 3 3
END
GEOEND
*---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+-
---8
MATERIAL 1.00 20.0 WATER
COMPOUND 2.0 3.0 1.0 8.0 WATER
*---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+-
---8
MATERIAL 0.001295 24.0 AIR
COMPOUND -0.78084 7.0 -0.20946 8.0 AIR
*---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+-
---8
MATERIAL 0.92 26.0 POLY
COMPOUND 1.0 6.0 2.0 3.0 POLY
*---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+-
---8
MATERIAL 5.0 10.811 2.34 27.0 BORON
*---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+-
---8
MATERIAL 0.93 28.0 BPE
COMPOUND -0.612 6.0 -0.116 3.0 -0.05 27.0BPE
COMPOUND -0.222 8.0 BPE
*---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+-
---8
*Blackhole
ASSIGNMAT 1.0 1.0
* Vacuum
ASSIGNMAT 2.0 2.0
* Tungsten collimator
ASSIGNMAT 23.0 3.0
*---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+-
---8
EMF
*---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+-
---8
* to calculate dose, need to activate fluscw.f
* USERWEIG 0. 0. 4. 0. 0. 0.
*---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+-
---8
USRYIELD 1401 8.0 32.0 -1.0 -2.0
1.0Neutron
USRYIELD 0.3 0.0 1.03.14159265 0.0 0.0&
*---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+-
---8
USRYIELD 1401 8.0 -42.0 -1.0 -2.0 1.0Neu3-2
USRYIELD 0.3 0.0 100.03.14159265 0.0 0.0&
*---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+-
---8
RANDOMIZE 1.0
*---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+-
---8
START 10000.
STOP
*******************************************
Jiajin Fan, Ph.D.
Research Associate, Radiation Oncology
Fox Chase Cancer Center
432 Rhawn Street, Philadelphia, PA 19111
Tel: (215) 214-1711
Fax:(215) 728-1148
Email: jiajin.fan@fccc.edu
*******************************************
This archive was generated by hypermail 2.1.6 : Mon Nov 27 2006 - 11:00:42 CET