About usrtrack and usrbin

From: Îäìõ (wu_zhen01@mails.tsinghua.edu.cn)
Date: Sat Dec 13 2003 - 10:24:39 CET

  • Next message: Alberto Fasso': "Re: About usrtrack and usrbin"

    Dear all,
            I have a compound target by tungsten and gold. It is hitted by a 9 MeV electron beam. I want to get the fluence in the tungsten, the gold and 1 meter away ahead of the target. I use usrtrack to get the energy spectrum of fluence, and use usrbin to get total fluence in the correspond point. In my result. the fluence got by usrtrack in the gold and 1 meter away ahead of the target is zero, but the total fluence got by usrbin is not zero. at the same time, the fluence in the tungsten got by usrtrack and usrbin are consistent. I dont know where the error is. please help me!
            Thanks a lot.

                                                                                                                                                                        Zhen Wu
                                                                                                                                                                    12/13/2003
            The following is my input file and my result:
    1. input file:
    TITLE
      Testing fluka: Beam of electron through compound target
    * by Zhen Wu 12/10/03
    *---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
    DEFAULTS NEW-DEFA
    *---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----*
    * Electron, 9 MeV Ek, circle 0.075 cm radius,
    * beam starts at (0,0,191.0) toward +Z axis
    BEAM -9.0E-3 0.0 0.0 0.075 0.0 -1.0ELECTRON
    BEAMPOS 0.0 0.0 18.5 0.0 0.0
    *---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----*
    GEOBEGIN 0.0 COMBINAT
                  outer rectangular and tungsten+gold target and shielding
    *
      RPP 1 -20000.0 20000.0 -20000.0 20000.0 -20000.0 20000.0
      RPP 2 -15000.0 15000.0 -15000.0 15000.0 -15000.0 15000.0
      RPP 3 -10000.0 10000.0 -10000.0 10000.0 -10000.0 10000.0
    * Compound Target
    * W
      RCC 4 0.0 0.0 18.85 0.0 0.0 19.1
                     0.5
    * Au
      RCC 5 0.0 0.0 19.1 0.0 0.0 19.2
                     0.5
    *detector
      RPP 6 -5.0 5.0 -5.0 5.0 124.0 114.0
      END
    *-reg-----or-----or-----or-----or-----or-----or-----or-----or-----or-----
      BL1 4 +1 -2
      VA2 4 +2 -4 -5 -6
      TU3 4 +4
      AU4 4 +5
    *-reg-----or-----or-----or-----or-----or-----or-----or-----or-----or-----
      VA5 4 +6
      END
    GEOEND
    *---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----*
    MATERIAL 79.0 196.967 19.32 15.0 0.0 0.0GOLD
    MATERIAL 82.0 207.2 11.35 17.0 0.0 0.0LEAD
    MATERIAL 74.0 184.0 19.3 23.0 0.0 0.0TUNGSTEN
    *---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----*
    ASSIGNMA 1.0 1.0
    ASSIGNMA 2.0 2.0
    ASSIGNMA 23.0 3.0
    ASSIGNMA 15.0 4.0
    ASSIGNMA 2.0 5.0
    *---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----*
    PHYSICS 2.0 EVAPORAT
    * all photonuclear phenomena are activated in all materials
    PHOTONUC 1111. 15.0
    PHOTONUC 1111. 17.0
    PHOTONUC 1111. 23.0
    * reduce inelastic interaction length of photons by 0.02 for all materials
    LAM-BIAS 0.0 0.02 0.0 7.0
    * non-analog transport of thermal neutrons in all regions
    * LOW-BIAS 73. 72. 0.85 1.0 6.0
    EMFFIX 23.0 0.05
    EMFFIX 15.0 0.05
    EMFFIX 17.0 0.05
    EMFCUT -5.0E-3 1.0E-6 0.0 1.0 8.0 1.0
    *---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----*
    OUTLEVEL 1.0 7.0
    *---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----*
    * neutron fluence in W target
    USRTRACK -1. 8. 41. 3. 1. 72.neutron
    USRTRACK 0.0196 1.E-14 &
    *---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----*
    * neutron fluence in Au target
    USRTRACK -1. 8. 41. 4. 1. 72.neutron
    USRTRACK 0.0196 1.E-14 &
    *---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----*
    * neutron total fluence in 0 degree,1 meter
    USRTRACK -1. 8. 41. 5. 1. 72.neutron
    USRTRACK 0.0196 1.E-14 &
    *---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----*
    * neutron total fluence in 0 degree,1 meter
    USRBIN 10. 8.0 61. 5.0 5.0 124.0neutron
    USRBIN -5.0 -5.0 114.0 1.0 1.0 1.0 &
    *---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----*
    * neutron total fluence in W
    USRBIN 11. 8.0 61. 0.5 0.0 19.1neutron
    USRBIN 0.0 0.0 18.85 1.0 1.0 5.0 &
    *---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----*
    * neutron total fluence in Au
    USRBIN 11. 8.0 61. 0.5 0.0 19.2neutron
    USRBIN 0.0 0.0 19.1 1.0 1.0 2.0 &
    RANDOMIZ 1.0
    START 1.0E+7
    STOP

    2. result
    usertrack result:
    1

       Track n. 1 "neutron " , generalized particle n. 8, region n. 3
          detector volume: 1.0000E+00 cm**3
          low energy neutrons scored from group 1 to group 72
          logar. energy binning from 1.9600E-02 to 1.9600E-02 GeV, 0 bins (ratio : 1.4816E+00)
          Data follow in a vector A(ie), format (1(5x,1p,10(1x,e11.4)))

          Low energy neutron data from group 1 to group 72 follow in a vector A(ig), format (1(5x,1p,10(1x,e11.4)))

           0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
           0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
           0.0000E+00 7.7668E-06 1.6520E-05 1.2455E-04 1.9992E-04 2.2520E-04 2.7827E-04 5.8891E-04 8.0825E-04 1.1837E-03
           2.0049E-03 2.2987E-03 3.3557E-03 4.5607E-03 4.8638E-03 6.7612E-03 8.1212E-03 7.7079E-03 8.1627E-03 9.6815E-03
           7.2599E-03 8.0053E-03 9.0499E-03 9.6502E-03 8.3322E-03 8.4891E-03 7.1115E-03 7.0571E-03 1.2297E-02 5.3844E-03
           3.0171E-03 7.3935E-03 3.9022E-03 1.9973E-02 2.1038E-02 5.8089E-03 2.0952E-02 0.0000E+00 0.0000E+00 0.0000E+00
           0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
           0.0000E+00 0.0000E+00
    1

       Track n. 2 "neutron " , generalized particle n. 8, region n. 4
          detector volume: 1.0000E+00 cm**3
          low energy neutrons scored from group 1 to group 72
          logar. energy binning from 1.9600E-02 to 1.9600E-02 GeV, 0 bins (ratio : 1.4816E+00)
          Data follow in a vector A(ie), format (1(5x,1p,10(1x,e11.4)))

          Low energy neutron data from group 1 to group 72 follow in a vector A(ig), format (1(5x,1p,10(1x,e11.4)))

           0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
           0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
           0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
           0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
           0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
           0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
           0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
           0.0000E+00 0.0000E+00
    1

       Track n. 3 "neutron " , generalized particle n. 8, region n. 5
          detector volume: 1.0000E+00 cm**3
          low energy neutrons scored from group 1 to group 72
          logar. energy binning from 1.9600E-02 to 1.9600E-02 GeV, 0 bins (ratio : 1.4816E+00)
          Data follow in a vector A(ie), format (1(5x,1p,10(1x,e11.4)))

          Low energy neutron data from group 1 to group 72 follow in a vector A(ig), format (1(5x,1p,10(1x,e11.4)))

           0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
           0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
           0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
           0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
           0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
           0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
           0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
           0.0000E+00 0.0000E+00
    usrbin result
    1
       Cartesian binning n. 1 "neutron " , generalized particle n. 8
          X coordinate: from -5.0000E+00 to 5.0000E+00 cm, 1 bins ( 1.0000E+01 cm wide)
          Y coordinate: from -5.0000E+00 to 5.0000E+00 cm, 1 bins ( 1.0000E+01 cm wide)
          Z coordinate: from 1.1400E+02 to 1.2400E+02 cm, 1 bins ( 1.0000E+01 cm wide)
          Data follow in a matrix A(ix,iy,iz), format (1(5x,1p,10(1x,e11.4)))

          accurate deposition along the tracks requested
          this is a track-length binning
           0.0000E+00
    1
       R - Z binning n. 2 "neutron " , generalized particle n. 8
          R coordinate: from 0.0000E+00 to 5.0000E-01 cm, 1 bins ( 5.0000E-01 cm wide)
          Z coordinate: from 1.8850E+01 to 1.9100E+01 cm, 5 bins ( 5.0000E-02 cm wide)
          axis coordinates: X = 0.0000E+00, Y = 0.0000E+00 cm
          Data follow in a matrix A(ir,iz), format (1(5x,1p,10(1x,e11.4)))

          accurate deposition along the tracks requested
          this is a track-length binning
           4.0541E-06 5.1506E-06 5.5572E-06 5.8179E-06 5.7916E-06
    1
       R - Z binning n. 3 "neutron " , generalized particle n. 8
          R coordinate: from 0.0000E+00 to 5.0000E-01 cm, 1 bins ( 5.0000E-01 cm wide)
          Z coordinate: from 1.9100E+01 to 1.9200E+01 cm, 2 bins ( 5.0000E-02 cm wide)
          axis coordinates: X = 0.0000E+00, Y = 0.0000E+00 cm
          Data follow in a matrix A(ir,iz), format (1(5x,1p,10(1x,e11.4)))

          accurate deposition along the tracks requested
          this is a track-length binning
           5.7035E-06 5.5627E-06


  • Next message: Alberto Fasso': "Re: About usrtrack and usrbin"

    This archive was generated by hypermail 2.1.6 : Sat Dec 13 2003 - 16:00:19 CET