question about USRYIELD

From: Fan, Ph.D., Jiajin (Jiajin.Fan@fccc.edu)
Date: Mon Nov 27 2006 - 04:25:27 CET

  • Next message: Vadim Talanov: "Messages in the ERR file"

    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
    *******************************************
     







  • Next message: Vadim Talanov: "Messages in the ERR file"

    This archive was generated by hypermail 2.1.6 : Mon Nov 27 2006 - 11:00:42 CET