Re: About usrtrack and usrbin

From: Alberto Fasso' (fasso@SLAC.Stanford.EDU)
Date: Tue Dec 16 2003 - 19:02:14 CET

  • Next message: Alfredo Ferrari: "FLUKA-2003.1 release"

    Dear Zhen Wu,

    your geometry has at least an error, and maybe more:
    > RPP 6 -5.0 5.0 -5.0 5.0 124.0 114.0
    Here, Zmax (114.0) is smaller than Zmin (124.)

    > RCC 4 0.0 0.0 18.85 0.0 0.0 19.1
    > 0.5
    > RCC 5 0.0 0.0 19.1 0.0 0.0 19.2
    > 0.5

    Maybe it is correct, but it seems unlikely to me. The above input defines
    two cylinders, one 0.25 cm long, and the other 19.2 cm long: is that what
    you wanted? Or rather the second cylinder was intended to be 0.1 cm long?
    In any case, it is a bad practice to define bodies "touching", i.e. with one
    face in common. Better let them overlap and "cut" them with planes.

    Alberto

    On Sat, 13 Dec 2003, [gb2312] Zhen Wu wrote:

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

    -- 
Alberto Fassò
SLAC-RP, P.O. Box 4349, MS 48, Stanford CA 94309
Phone: (1 650) 926 4762   Fax: (1 650) 926 3569
fasso@slac.stanford.edu
  • Next message: Alfredo Ferrari: "FLUKA-2003.1 release"

    This archive was generated by hypermail 2.1.6 : Wed Dec 17 2003 - 09:15:26 CET