Re: Neutron problem

From: Sebastien WURTH (wurth@ipno.in2p3.fr)
Date: Wed Sep 20 2006 - 10:29:16 CEST

  • Next message: Stefan Roesler: "Re: Neutron problem"

    Hello Nisy,

    If you are using a routine which converts fluence to ambiant dose
    equivalent by means of USERWEIG card (maybe some deq99c if I figure well
    your USERTRACK SDUMs) you will indeeed get results in pSv/primary,
    pSv/neutrons here.
    I cant' see anything wrong in your work, maybe the FLUKA experts will.
    But in terms of radiation safety, a difference of a factor 6 doesn't
    seem to be too important. Of course, it depends on what accuracy you are
    looking at...
    It depends also on what you are calculating in your "simple calculation
    of dose transmitted", if you do not use the same conversion factors,
    what method you use, is it an ambiant dose equivalent calculation, etc,
    you can rapidly become a 6 factor difference in your two results without
    anything being false or wrong.

    Hope it helps a little.
    Regards.
    Sebastien.

    Nisy Ipe a écrit :

    >Hello all,
    >
    >I have a simple neutron problem - a 252Cf source inside a tunnel with
    >concrete walls but no roof. The californium source has a dose rate
    >profile and so a source routine was written (attached). Ambient dose
    >equivalent for all particles, neutrons and photons were scored at
    >various locations using USRTRACK and compiling with ustsuw.f. Input and
    >results are attached. The results (calif. dose) read p/cmq/pr but I
    >think in my case this would be pSv/n? Multiplying the mass by 2.4 E12
    >n/g of 252Cf and the results by the number of neutrons, I obtained the
    >ambient dose equivalents at various locations.
    >
    >If I do a simple calculation of dose transmitted through the concrete
    >wall I get results which are 6 x higher tthan the FLUKA results. So I
    >think I am doing something wrong but can't figure it out. Could the
    >experts please help?
    >
    >Thanks and Regards,
    >
    >Nisy
    >
    >
    >------------------------------------------------------------------------
    >
    >*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
    >TITLE
    >Dose rates from 252Cf, h=4.35 m
    >*9/17/06
    >*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
    >*
    >BEAM -2.E-03 NEUTRON
    >SOURCE
    >*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
    >*TNA at 1.85 m (4.35 m high) above center,11 m in y 1.65 m from tunnel center (z)
    >BEAMPOS 185.0 1100.0 -165.0 0. 0.
    >GEOBEGIN COMBINAT
    > 0 -10 Beam strikes object
    >*Black Box
    >*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
    > RPP 1 -5000000.0 +5000000.0 -5000000.0
    > +5000000.0 -5000000.0 +5000000.0
    >* VAC Box
    > RPP 2 -50000.0 +50000.0 -50000.0
    > +50000.0 -50000.0 +50000.0
    >* AIR Box
    > RPP 3 -310.0 +30250.0 -3250.0
    > +2750.0 -1250.0 +2750.0
    >* ConcreteO
    > RPP 4 -290.0 +250.0 -1700.0
    > +1700.0 -650.0 +550.0
    >* AIR Inside
    > RPP 5 -250.0 +250.0 -1675.0
    > +1675.0 -625.0 +500.0
    >* Office Concrete
    > RPP 6 -290.0 +250.0 -1900.0
    > +1900.0 -175.0 +175.0
    >* Office Air
    > RPP 7 -250.0 +250.0 -1900.0
    > +1900.0 -150.0 +150.0
    >* RPP for object which is 30 cm x 30 cm x 30 cm
    >*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
    > RPP 8 +135.0 +165.0 +1085.0
    > +1115.0 -15.0 +15.0
    >*Plane for soil
    > YZP 9 -270.0
    >* Detectors in Check out 1 foot from wall at 0.5,1 and 4.35 m heights,30cmSphere
    >*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
    >* Check Out centered at 1730, -205, h=0.5, 1.0 and 4.35 m
    > SPH 10 -200.0 1730.0 -205.0
    > 15.0
    > SPH 11 -150.0 1730.0 -205.0
    > 15.0
    > SPH 12 185.0 1730.0 -205.0
    > 15.0
    >*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
    >*Beyond Check Out Wall centered at 1930, -205, h=0.5, 1.,4.35 m
    > SPH 13 -200.0 1930.0 -205.0
    > 15.0
    > SPH 14 -150.0 1930.0 -205.0
    > 15.0
    > SPH 15 185.0 1930.0 -205.0
    > 15.0
    >*Outside Rolling Door centered at 1930, 0, h=0.5, 1.,4.35 m
    > SPH 16 -200.0 1930.0 0.0
    > 15.0
    > SPH 17 -150.0 1930.0 0.0
    > 15.0
    > SPH 18 185.0 1930.0 0.0
    > 15.0
    >*Beyond Plant wall centered at 1930, 205, h=0.5, 1.,4.35 m
    > SPH 19 100.0 1930.0 205.0
    > 15.0
    > SPH 20 150.0 1930.0 205.0
    > 15.0
    > SPH 21 200.0 1930.0 205.0
    > 15.0
    >*Inside Plant Room centered at 100, 1730, 205, h=1, 1.5, 2 m
    > SPH 22 100.0 1730.0 205.0
    > 15.0
    > SPH 23 150.0 1730.0 205.0
    > 15.0
    > SPH 24 200.0 1730.0 205.0
    > 15.0
    >*Inside Plant Room centered at 100, 1730, 305, h=1, 1.5, 2 m
    > SPH 25 100.0 1730.0 305.0
    > 15.0
    > SPH 26 150.0 1730.0 305.0
    > 15.0
    > SPH 27 200.0 1730.0 305.0
    > 15.0
    >*Inside Plant Room centered at 100, 1730, 405, h=1, 1.5, 2 m
    > SPH 28 100.0 1730.0 405.0
    > 15.0
    > SPH 29 150.0 1730.0 405.0
    > 15.0
    > SPH 30 200.0 1730.0 405.0
    > 15.0
    >*Outside 25 cm wall
    > SPH 31 100.0 1100.0 -680.0
    > 15.0
    > SPH 32 150.0 1100.0 -680.0
    > 15.0
    > SPH 33 200.0 1100.0 -680.0
    > 15.0
    >*Exclusion Zone beyond 25 cm wall
    > SPH 34 100.0 1100.0 -1220.0
    > 15.0
    > SPH 35 150.0 1100.0 -1220.0
    > 15.0
    > SPH 36 200.0 1100.0 -1220.0
    > 15.0
    >*Outside 50 cm wall
    > SPH 37 100.0 1100.0 580.0
    > 15.0
    > SPH 38 150.0 1100.0 580.0
    > 15.0
    > SPH 39 200.0 1100.0 580.0
    > 15.0
    >*Near Exclusion Zone Beyond 50 cm wall
    > SPH 40 100.0 1100.0 2720.0
    > 15.0
    > SPH 41 150.0 1100.0 2720.0
    > 15.0
    > SPH 42 200.0 1100.0 2720.0
    > 15.0
    >*Near Exclusion Zone Beyond Small Tunnel
    > SPH 43 100.0 2720.0 0.0
    > 15.0
    > SPH 44 150.0 2720.0 0.0
    > 15.0
    > SPH 45 200.0 2720.0 0.0
    > 15.0
    > END
    >* black hole
    >*.........or.....or.....or.....or.....or.....or.....or.....or.....or.....
    > BH1 5 +1 -2
    >* Vacuum around
    >*.........or.....or.....or.....or.....or.....or.....or.....or.....or.....
    > VA2 5 +2 -3
    >* Exclusion zone
    >*.........or.....or.....or.....or.....or.....or.....or.....or.....or.....
    > AI3 5 +3 -4 -6 -9 -10 -11 -12 -13
    > -14 -15 -16 -17 -18 -19 -20 -21
    > -22 -23 -24 -25 -26 -27 -28 -29
    > -30 -31 -32 -33 -34 -35 -36 -37
    > -38 -39 -40 -41 -42 -43 -44 -45
    >*.........or.....or.....or.....or.....or.....or.....or.....or.....or.....
    >* Air Inside small tunnel
    > AI4 5 +7 -5
    >* Air Inside Large tunnel
    > AI5 5 +5 -8
    >* Concrete for small tunnel
    > CO6 5 +6 -7 -4
    >* Concrete for large tunnel
    > CO7 5 +4 -5 -7
    >*.........or.....or.....or.....or.....or.....or.....or.....or.....or.....
    >*Soil below
    > SO8 5 +3 -4 -6 +9
    >*.........or.....or.....or.....or.....or.....or.....or.....or.....or.....
    >*OBJECT
    > OB9 5 +8
    >*Check Out
    > D10 5 +10
    > D11 5 +11
    > D12 5 +12
    >*Beyond Check Out Wall
    > D13 5 +13
    > D14 5 +14
    > D15 5 +15
    >*Outside Rolling Door
    > D16 5 +16
    > D17 5 +17
    > D18 5 +18
    >*Beyond Plant Wall
    > D19 5 +19
    > D20 5 +20
    > D21 5 +21
    >*Inside Plant Room at z = 205
    > D22 5 +22
    > D23 5 +23
    > D24 5 +24
    >*Inside Plant Wall at z =305
    > D25 5 +25
    > D26 5 +26
    > D27 5 +27
    >*Inside Plant Wall at z =405
    > D28 5 +28
    > D29 5 +29
    > D30 5 +30
    >*Outside 25 cm Wall
    > D31 5 +31
    > D32 5 +32
    > D33 5 +33
    >*Near Exclusion Zone Outside 25 cm Wall
    > D34 5 +34
    > D35 5 +35
    > D36 5 +36
    >*Outside 50 cm Wall
    > D37 5 +37
    > D38 5 +38
    > D39 5 +39
    >*Near Exclusion Zone Outside 50 cm Wall
    > D40 5 +40
    > D41 5 +41
    > D42 5 +42
    >*Near Exclusion Zone Outside Small Tunnel
    > D43 5 +43
    > D44 5 +44
    > D45 5 +45
    > END
    >*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
    >*GEOEND 1000. 0. 1000. -1000. 0. -1000. DEBUG
    >*GEOEND 200. 1. 200. &
    >*STOP
    >GEOEND
    >MATERIAL 1. 1.0079 0.000089 3. 0.0 0. HYDROGEN
    >MATERIAL 6. 12.011 2.265 6. 0.0 0. CARBON
    >MATERIAL 8. 15.9994 0.001429 8. 0.0 0. OXYGEN
    >MATERIAL 12. 24.305 1.738 9. 0.0 0. MAGNESIU
    >MATERIAL 13. 26.982 2.6989 10. 0.0 0. ALUMINUM
    >MATERIAL 26. 55.847 7.87 11. 0.0 0. IRON
    >MATERIAL 29. 63.546 8.96 12. 0.0 0. COPPER
    >MATERIAL 22. 47.867 6.8281 13. 0.0 0. TITANIUM
    >MATERIAL 18.0 39.948 1.662E-3 23. 0.0 0. ARGON
    >MATERIAL 0.0012048 24. 0.0 0. AIR
    >COMPOUND -1.24E-4 6.0 -0.755267 7.0 -0.231781 8. AIR
    >COMPOUND -0.012827 23.0 0.0 0.0 0.0 0. AIR
    >MATERIAL 14. 28.0855 2.33 26. 0.0 0. SILICON
    >MATERIAL 20. 40.078 1.55 27. 0.0 0. CALCIUM
    >MATERIAL 82. 207.2 11.35 28. 0.0 0. LEAD
    >MATERIAL 11. 22.989768 0.971 29. 0.0 0. SODIUM
    >MATERIAL 19. 39.0983 0.862 30. 0.0 0. POTASSIU
    >* H C O Na Mg Al Si K Ca Fe
    >MATERIAL 0.0 0.0 2.35 31. 0.0 0. CONCRETE
    >COMPOUND -0.01 3. -0.001 6.0 -0.529107 8. CONCRETE
    >COMPOUND -0.016 29. -0.002 9.0 -0.033872 10. CONCRETE
    >COMPOUND -0.337021 26. -0.013 30.0 -0.044 27. CONCRETE
    >COMPOUND -0.014 11. CONCRETE
    >*
    >MATERIAL 0.0 0.0 1.0 32.0 TISSUE
    >COMPOUND -76.2 8.0 -11.1 6.0 -10.1 3. TISSUE
    >COMPOUND -2.6 7.0 TISSUE
    >* H,O, Al, Si
    >MATERIAL 0.0 0.0 1.7 33. 0.0 0. SOIL
    >COMPOUND +9.77E+24 3.0 3.48E+22 8.0 4.88E+22 10. SOIL
    >COMPOUND +1.16E+22 26.0 SOIL
    >*
    >MATERIAL 0.0 0.0 1.0 34. 0.0 0. WATER
    >COMPOUND 2.0 3.0 1.0 8.0 WATER
    >*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
    >* External Black Hole
    >ASSIGNMAT 1.0 1.0
    >ASSIGNMAT 2.0 2.0 ! vacuum
    >ASSIGNMAT 24.0 3.0 5.0 1.0 ! air
    >* Concrete
    >ASSIGNMAT 31.0 6.0 7.0 1.0! concrete
    >*Soil
    >ASSIGNMAT 33.0 8.0 ! soil
    >*Object changed to air
    >ASSIGNMAT 24.0 9.0 ! air
    >*Water Object
    >*ASSIGNMAT 34.0 9.0 ! Water
    >*Steel Object
    >*ASSIGNMAT 11.0 9.0 ! Iron
    >*Detector (air)
    >ASSIGNMAT 24.0 10.0 45.0 1.0 ! air
    >*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
    >*
    >*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
    >SCORE 208.0 211.0 201. 206.
    >*All particles
    >*Check Out
    >USRTRACK 1.0 201.0 -21. 10.14137.1669 1. AMB74-A10
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 201.0 -21. 11.14137.1669 1. AMB74-A11
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 201.0 -21. 12.14137.1669 1. AMB74-A12
    >USRTRACK 0.0 0.0 &
    >*Beyond Check Out Wall
    >USRTRACK 1.0 201.0 -21. 13.14137.1669 1. AMB74-A13
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 201.0 -21. 14.14137.1669 1. AMB74-A14
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 201.0 -21. 15.14137.1669 1. AMB74-A15
    >USRTRACK 0.0 0.0 &
    >*Outside Rolling Door
    >USRTRACK 1.0 201.0 -21. 16.14137.1669 1. AMB74-A16
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 201.0 -21. 17.14137.1669 1. AMB74-A17
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 201.0 -21. 18.14137.1669 1. AMB74-A18
    >USRTRACK 0.0 0.0 &
    >*Beyond Plant Wall
    >USRTRACK 1.0 201.0 -21. 19.14137.1669 1. AMB74-A19
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 201.0 -21. 20.14137.1669 1. AMB74-A20
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 201.0 -21. 21.14137.1669 1. AMB74-A21
    >USRTRACK 0.0 0.0 &
    >* Inside Plant Wall at 205
    >USRTRACK 1.0 201.0 -21. 22.14137.1669 1. AMB74-A22
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 201.0 -21. 23.14137.1669 1. AMB74-A23
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 201.0 -21. 24.14137.1669 1. AMB74-A24
    >USRTRACK 0.0 0.0
    >* Inside Plant Wall at 305
    >USRTRACK 1.0 201.0 -21. 25.14137.1669 1. AMB74-A25
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 201.0 -21. 26.14137.1669 1. AMB74-A26
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 201.0 -21. 27.14137.1669 1. AMB74-A27
    >USRTRACK 0.0 0.0 &
    >* Inside Plant Wall at 405
    >USRTRACK 1.0 201.0 -21. 28.14137.1669 1. AMB74-A28
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 201.0 -21. 29.14137.1669 1. AMB74-A29
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 201.0 -21. 30.14137.1669 1. AMB74-A30
    >USRTRACK 0.0 0.0 &
    >* Outside 25 cm
    >USRTRACK 1.0 201.0 -21. 31.14137.1669 1. AMB74-A31
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 201.0 -21. 32.14137.1669 1. AMB74-A32
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 201.0 -21. 33.14137.1669 1. AMB74-A33
    >USRTRACK 0.0 0.0 &
    >*Near Exclusion Zone Outside 25 cm
    >USRTRACK 1.0 201.0 -21. 34.14137.1669 1. AMB74-A34
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 201.0 -21. 35.14137.1669 1. AMB74-A35
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 201.0 -21. 36.14137.1669 1. AMB74-A36
    >USRTRACK 0.0 0.0 &
    >*
    >*Outside 50 cm
    >USRTRACK 1.0 201.0 -21. 37.14137.1669 1. AMB74-A37
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 201.0 -21. 38.14137.1669 1. AMB74-A38
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 201.0 -21. 39.14137.1669 1. AMB74-A39
    >USRTRACK 0.0 0.0 &
    >*Exclusion Zone Outside 50 cm
    >USRTRACK 1.0 201.0 -21. 40.14137.1669 1. AMB74-A40
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 201.0 -21. 41.14137.1669 1. AMB74-A41
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 201.0 -21. 42.14137.1669 1. AMB74-A42
    >USRTRACK 0.0 0.0 &
    >*Exclusion Zone Outside Small Tunnel
    >USRTRACK 1.0 201.0 -21. 43.14137.1669 1. AMB74-A43
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 201.0 -21. 44.14137.1669 1. AMB74-A44
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 201.0 -21. 45.14137.1669 1. AMB74-A45
    >USRTRACK 0.0 0.0 &
    >* Neutrons
    >*Check Out
    >USRTRACK 1.0 8.0 -22. 10.14137.1669 1. AMB74-N10
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -22. 11.14137.1669 1. AMB74-N11
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -22. 12.14137.1669 1. AMB74-N12
    >USRTRACK 0.0 0.0 &
    >*Beyond Check Out Wall
    >USRTRACK 1.0 8.0 -22. 13.14137.1669 1. AMB74-N13
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -22. 14.14137.1669 1. AMB74-N14
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -22. 15.14137.1669 1. AMB74-N15
    >USRTRACK 0.0 0.0 &
    >*Outside Rolling Door
    >USRTRACK 1.0 8.0 -22. 16.14137.1669 1. AMB74-N16
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -22. 17.14137.1669 1. AMB74-N17
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -22. 18.14137.1669 1. AMB74-N18
    >USRTRACK 0.0 0.0 &
    >*Beyond Plant Wall
    >USRTRACK 1.0 8.0 -22. 19.14137.1669 1. AMB74-N19
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -22. 20.14137.1669 1. AMB74-N20
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -22. 21.14137.1669 1. AMB74-N21
    >USRTRACK 0.0 0.0 &
    >*Inside Plant Room at z=205
    >USRTRACK 1.0 8.0 -22. 22.14137.1669 1. AMB74-N22
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -22. 23.14137.1669 1. AMB74-N23
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -22. 24.14137.1669 1. AMB74-N24
    >USRTRACK 0.0 0.0 &
    >*Inside Plant Room at z=305
    >USRTRACK 1.0 8.0 -22. 25.14137.1669 1. AMB74-N25
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -22. 26.14137.1669 1. AMB74-N26
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -22. 27.14137.1669 1. AMB74-N27
    >USRTRACK 0.0 0.0 &
    >*Inside Plant Room at z=405
    >USRTRACK 1.0 8.0 -22. 28.14137.1669 1. AMB74-N28
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -22. 29.14137.1669 1. AMB74-N29
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -22. 30.14137.1669 1. AMB74-N30
    >USRTRACK 0.0 0.0 &
    >*Outside 25 cm Wall
    >USRTRACK 1.0 8.0 -22. 31.14137.1669 1. AMB74-N31
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -22. 32.14137.1669 1. AMB74-N32
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -22. 33.14137.1669 1. AMB74-N33
    >USRTRACK 0.0 0.0 &
    >*Near Exclusion Zone Outside 25 cm Wall
    >USRTRACK 1.0 8.0 -22. 34.14137.1669 1. AMB74-N34
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -22. 35.14137.1669 1. AMB74-N35
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -22. 36.14137.1669 1. AMB74-N36
    >USRTRACK 0.0 0.0 &
    >*Outside 50 cm Wall
    >USRTRACK 1.0 8.0 -22. 37.14137.1669 1. AMB74-N37
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -22. 38.14137.1669 1. AMB74-N38
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -22. 39.14137.1669 1. AMB74-N39
    >USRTRACK 0.0 0.0 &
    >*Near Exclusion Zone Outside 50 cm Wall
    >USRTRACK 1.0 8.0 -22. 40.14137.1669 1. AMB74-N40
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -22. 41.14137.1669 1. AMB74-N41
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -22. 42.14137.1669 1. AMB74-N42
    >USRTRACK 0.0 0.0 &
    >* Near Exclusion Zone Beyond Small Tunnel
    >USRTRACK 1.0 8.0 -22. 43.14137.1669 1. AMB74-N43
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -22. 44.14137.1669 1. AMB74-N44
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -22. 45.14137.1669 1. AMB74-N45
    >USRTRACK 0.0 0.0 &
    >*Photons
    >*Check Out
    >USRTRACK 1.0 7.0 -23. 10.14137.1669 1. AMB74-P10
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 7.0 -23. 11.14137.1669 1. AMB74-P11
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -23. 12.14137.1669 1. AMB74-P12
    >USRTRACK 0.0 0.0 &
    >*Beyond Check Out Wall
    >USRTRACK 1.0 8.0 -23. 13.14137.1669 1. AMB74-P13
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -23. 14.14137.1669 1. AMB74-P14
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -23. 15.14137.1669 1. AMB74-P15
    >USRTRACK 0.0 0.0 &
    >*Outside Rolling Door
    >USRTRACK 1.0 8.0 -23. 16.14137.1669 1. AMB74-P16
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -23. 17.14137.1669 1. AMB74-P17
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -23. 18.14137.1669 1. AMB74-P18
    >USRTRACK 0.0 0.0 &
    >*Beyond Plant Wall
    >USRTRACK 1.0 8.0 -23. 19.14137.1669 1. AMB74-P19
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -23. 20.14137.1669 1. AMB74-P20
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -23. 21.14137.1669 1. AMB74-P21
    >USRTRACK 0.0 0.0 &
    >*Inside Plant Room at z = 205 cm
    >USRTRACK 1.0 8.0 -23. 22.14137.1669 1. AMB74-P22
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -23. 23.14137.1669 1. AMB74-P23
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -23. 24.14137.1669 1. AMB74-P24
    >USRTRACK 0.0 0.0 &
    >*Inside Plant Room at z =305 cm
    >USRTRACK 1.0 8.0 -23. 25.14137.1669 1. AMB74-P25
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -23. 26.14137.1669 1. AMB74-P26
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -23. 27.14137.1669 1. AMB74-P27
    >USRTRACK 0.0 0.0 &
    >*Inside Plant Room at z =405 cm
    >USRTRACK 1.0 8.0 -23. 28.14137.1669 1. AMB74-P28
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -23. 29.14137.1669 1. AMB74-P29
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -23. 30.14137.1669 1. AMB74-P30
    >USRTRACK 0.0 0.0 &
    >*Outside 25 cm Wall
    >USRTRACK 1.0 8.0 -23. 31.14137.1669 1. AMB74-P31
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -23. 32.14137.1669 1. AMB74-P32
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -23. 33.14137.1669 1. AMB74-P33
    >USRTRACK 0.0 0.0 &
    >*Near Exclusion Zone Outside 25 cm Wall
    >USRTRACK 1.0 8.0 -23. 34.14137.1669 1. AMB74-P34
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -23. 35.14137.1669 1. AMB74-P35
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -23. 36.14137.1669 1. AMB74-P36
    >USRTRACK 0.0 0.0 &
    >*Outside 50 cm Wall
    >USRTRACK 1.0 8.0 -23. 37.14137.1669 1. AMB74-P37
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -23. 38.14137.1669 1. AMB74-P38
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -23. 39.14137.1669 1. AMB74-P39
    >USRTRACK 0.0 0.0 &
    >*Near Exclusion Zone Outside 50 cm Wall
    >USRTRACK 1.0 8.0 -23. 40.14137.1669 1. AMB74-P40
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -23. 41.14137.1669 1. AMB74-P41
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -23. 42.14137.1669 1. AMB74-P42
    >USRTRACK 0.0 0.0 &
    >*Near Exclusion Zone Outside Small Tunnel
    >USRTRACK 1.0 8.0 -23. 43.14137.1669 1. AMB74-P43
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -23. 44.14137.1669 1. AMB74-P44
    >USRTRACK 0.0 0.0 &
    >USRTRACK 1.0 8.0 -23. 45.14137.1669 1. AMB74-P45
    >USRTRACK 0.0 0.0 &
    >*
    >USERWEIG 0.0 0.0 1.
    >RANDOMIZE 1.0
    >*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
    >START 500000.0
    >STOP
    >
    >
    >
    >------------------------------------------------------------------------
    >
    >*$ CREATE SOURCE.FOR
    >*COPY SOURCE
    >*
    >*=== source ===========================================================*
    >*
    > SUBROUTINE SOURCE ( NOMORE )
    >
    > INCLUDE '(DBLPRC)'
    > INCLUDE '(DIMPAR)'
    > INCLUDE '(IOUNIT)'
    >*----------------------------------------------------------------------*
    >* Version for isotropic Cf source.
    >*
    > DIMENSION CFCUM(0:72), ENEDGE(73), ANG(0:40), CUMANG(0:40)
    >*----------------------------------------------------------------------*
    >*
    > INCLUDE '(BEAMCM)'
    > INCLUDE '(CASLIM)'
    > INCLUDE '(FHEAVY)'
    > INCLUDE '(FLKSTK)'
    > INCLUDE '(IOIOCM)'
    > INCLUDE '(LTCLCM)'
    > INCLUDE '(PAPROP)'
    > INCLUDE '(SOURCM)'
    > INCLUDE '(SUMCOU)'
    >*
    > LOGICAL LFIRST
    > SAVE LFIRST
    > DATA LFIRST / .TRUE. /
    >*----------------------------------------------------------------------*
    >* Neutron energy group boundaries
    > DATA ENEDGE /
    > & 1.9600E-02, 1.7500E-02, 1.4918E-02, 1.3499E-02, 1.2214E-02,
    > & 1.1052E-02, 1.0000E-02, 9.0484E-03, 8.1873E-03, 7.4082E-03,
    > & 6.7032E-03, 6.0653E-03, 5.4881E-03, 4.9659E-03, 4.4933E-03,
    > & 4.0657E-03, 3.6788E-03, 3.3287E-03, 3.0119E-03, 2.7253E-03,
    > & 2.4660E-03, 2.2313E-03, 2.0190E-03, 1.8268E-03, 1.6530E-03,
    > & 1.4957E-03, 1.3534E-03, 1.2246E-03, 1.1080E-03, 1.0026E-03,
    > & 9.0718E-04, 8.2085E-04, 7.4274E-04, 6.0810E-04, 4.9787E-04,
    > & 4.0762E-04, 3.3373E-04, 2.7324E-04, 2.2371E-04, 1.8316E-04,
    > & 1.4996E-04, 1.2277E-04, 8.6517E-05, 5.2475E-05, 3.1828E-05,
    > & 2.1852E-05, 1.5034E-05, 1.0332E-05, 7.1018E-06, 4.8809E-06,
    > & 3.3546E-06, 2.3054E-06, 1.5846E-06, 1.0446E-06, 6.8871E-07,
    > & 4.5400E-07, 2.7537E-07, 1.6702E-07, 1.0130E-07, 6.1442E-08,
    > & 3.7267E-08, 2.2603E-08, 1.5535E-08, 1.0677E-08, 7.3375E-09,
    > & 5.0435E-09, 3.4662E-09, 2.3824E-09, 1.6374E-09, 1.1254E-09,
    > & 6.8257E-10, 4.1400E-10, 1.0000E-14/
    >* Cumulative Cf spectrum
    > DATA CFCUM / 0.D0,
    >*...+....1....+....2....+....3....+....4....+....5....+....6....+....7..
    > & 1.919922147E-05,1.076421528E-04,2.985396423E-04,6.750280793E-04,
    > & 1.460097913E-03,2.919855883E-03,5.440437967E-03,9.517761940E-03,
    > & 1.572873215E-02,2.550738486E-02,3.904978251E-02,5.635569494E-02,
    > & 7.365912720E-02,9.966894989E-02,1.279101323E-01,1.650685484E-01,
    > & 2.022289923E-01,2.393916485E-01,2.836610364E-01,3.279253392E-01,
    > & 3.733559139E-01,4.193675039E-01,4.654041752E-01,5.064859139E-01,
    > & 5.475773136E-01,5.886592151E-01,6.297539805E-01,6.708701087E-01,
    > & 6.974131472E-01,7.239697428E-01,7.505235450E-01,7.770757115E-01,
    > & 8.301848328E-01,8.661373176E-01,8.849342051E-01,9.037310414E-01,
    > & 9.225258728E-01,9.413224153E-01,9.601180671E-01,9.756795138E-01,
    > & 9.815263394E-01,9.879075224E-01,9.948263962E-01,9.971427050E-01,
    > & 9.985910051E-01,9.990648305E-01,9.995399469E-01,9.996648243E-01,
    > & 9.997897420E-01,9.999146525E-01,9.999395164E-01,9.999643698E-01,
    > & 9.999838018E-01,9.999909490E-01,9.999946297E-01,9.999962702E-01,
    > & 9.999979108E-01,9.999995515E-01,9.999997034E-01,9.999998553E-01,
    > & 9.999999446E-01,9.999999646E-01,9.999999846E-01,9.999999885E-01,
    > & 9.999999924E-01,9.999999963E-01,9.999999981E-01,9.999999988E-01,
    > & 9.999999995E-01,9.999999997E-01,9.999999999E-01,1.000000000E+00/
    >* Angle cosines
    > DATA ANG/ -1.0D0, -0.95D0, -0.9D0, -0.85D0, -0.8D0, -0.75D0,
    > & -0.7D0, -0.65D0, -0.6D0, -0.55D0, -0.5D0, -0.45D0,
    > & -0.4D0, -0.35D0, -0.3D0, -0.25D0, -0.2D0, -0.15D0,
    > & -0.1D0, -5.D-02, 0.0D0, 5.D-02, 0.1D0, 0.15D0,
    > & 0.2D0, 0.25D0, 0.3D0, 0.35D0, 0.4D0, 0.45D0,
    > & 0.5D0, 0.55D0, 0.6D0, 0.65D0, 0.7D0, 0.75D0,
    > & 0.8D0, 0.85D0, 0.9D0, 0.95D0, 1.0D0 /
    >* Cumulative angular distribution
    > DATA CUMANG / 0.0D0, 4.3787904871D-4, 9.309108266D-4,
    > & 1.5080318037D-3, 2.1835838416D-3, 2.9743548826D-3,
    > & 3.8909576476D-3, 4.8923191479D-3, 5.9862770681D-3,
    > & 7.1813938468D-3, 8.4870236975D-3, 9.9133858209D-3,
    > & 1.1471644393D-2, 1.3097707545D-2, 1.4776914040D-2,
    > & 1.6511000723D-2, 1.8301761200D-2, 2.0151047694D-2,
    > & 2.2060772967D-2, 2.4032912292D-2, 2.6069505497D-2,
    > & 2.9084203728D-2, 3.3546756876D-2, 4.0152519402D-2,
    > & 4.9930800378D-2, 6.4405250301D-2, 8.5831276094D-2,
    > & 0.1175474784D0, 0.1569148680D0, 0.1981130258D0,
    > & 0.2412270912D0, 0.2863461629D0, 0.3335634830D0,
    > & 0.3829766299D0, 0.4346877199D0, 0.4919529577D0,
    > & 0.5559660349D0, 0.6275220820D0, 0.7075099237D0,
    > & 0.8189694928D0, 1.0D0 /
    >*----------------------------------------------------------------------*
    >*...+....1....+....2....+....3....+....4....+....5....+....6....+....7..
    >*----------------------------------------------------------------------*
    > NOMORE = 0
    >* +-------------------------------------------------------------------*
    >* | First call initializations:
    > IF ( LFIRST ) THEN
    >* | *** The following 3 cards are mandatory ***
    > TKESUM = ZERZER
    > LFIRST = .FALSE.
    > LUSSRC = .TRUE.
    >* | *** User initialization ***
    > write(lunout,*)
    > write(lunout,'(a,132a)') ("*",i=1,132)
    > write(lunout,*)
    > write(lunout,*)
    > write(lunout,'(a)') "Californium source. Version corrected"
    > write(lunout,*)
    > write(lunout,*)
    > write(lunout,'(a,132a)') ("*",i=1,132)
    > write(lunout,*)
    > END IF
    >* |
    >* +-------------------------------------------------------------------*
    >* Sample the energy group
    > XI = FLRNDM(DUMMY)
    > DO 500 K = 1, 72
    > IF(XI .LE. CFCUM(K)) THEN
    > NEUGRP = K
    > ENERGY = ENEDGE(K) -
    > & (XI-CFCUM(K-1))*(ENEDGE(K)-ENEDGE(K+1))/(CFCUM(K)-CFCUM(K-1))
    > GO TO 600
    > END IF
    > 500 CONTINUE
    > STOP ' Failed to sample the energy group'
    > 600 CONTINUE
    > kount=kount+1
    >*
    >* Sample the cosine of the polar angle
    > XI = FLRNDM(XI)
    > DO 300 K = 1, 40
    > IF(XI .LE. CUMANG(K)) THEN
    > COSTHE = ANG(K-1) +
    > & (XI - CUMANG(K-1))*(ANG(K)-ANG(K-1))/(CUMANG(K)-CUMANG(K-1))
    > GO TO 400
    > END IF
    > 300 CONTINUE
    > STOP ' Failed to sample the polar angle cosine'
    > 400 CONTINUE
    > SINTHE = SQRT(ONEONE - COSTHE**2)
    > IF(FLRNDM(COSTHE) .LE. HLFHLF) SINTHE = -SINTHE
    >* Samples the azimuthal angle
    > PHI = TWOPIP * FLRNDM(COSTHE)
    > COSPHI = COS(PHI)
    >* Npflka is the stack counter: of course any time source is called it
    >* must be =0
    > NPFLKA = NPFLKA + 1
    >* Wt is the weight of the particle
    > WTFLK (NPFLKA) = ONEONE
    > WEIPRI = WEIPRI + WTFLK (NPFLKA)
    >* Particle type (1=proton.....). Ijbeam is the type set by the BEAM
    >* card
    >* +-------------------------------------------------------------------*
    >* | Heavy ion:
    > IF ( IJBEAM .EQ. -2 ) THEN
    > IJHION = IPROZ * 1000 + IPROA
    > IJHION = IJHION * 100 + KXHEAV
    > IONID = IJHION
    > CALL DCDION ( IONID )
    > CALL SETION ( IONID )
    > ILOFLK (NPFLKA) = IJHION
    >* |
    >* +-------------------------------------------------------------------*
    >* | Normal hadron:
    > ELSE
    > IONID = IJBEAM
    > ILOFLK (NPFLKA) = IJBEAM
    > END IF
    >* |
    >* +-------------------------------------------------------------------*
    >* From this point .....
    >* Particle generation (1 for primaries)
    > LOFLK (NPFLKA) = 1
    >* User dependent flag:
    > LOUSE (NPFLKA) = 0
    >* User dependent spare variables:
    > DO 100 ISPR = 1, MKBMX1
    > SPAREK (ISPR,NPFLKA) = ZERZER
    > 100 CONTINUE
    >* User dependent spare flags:
    > DO 200 ISPR = 1, MKBMX2
    > ISPARK (ISPR,NPFLKA) = 0
    > 200 CONTINUE
    >* Save the track number of the stack particle:
    > ISPARK (MKBMX2,NPFLKA) = NPFLKA
    > NPARMA = NPARMA + 1
    > NUMPAR (NPFLKA) = NPARMA
    > NEVENT (NPFLKA) = 0
    > DFNEAR (NPFLKA) = +ZERZER
    >* ... to this point: don't change anything
    >* Particle age (s)
    > AGESTK (NPFLKA) = +ZERZER
    > AKNSHR (NPFLKA) = -TWOTWO
    >* Group number for "low" energy neutrons
    > IGROUP (NPFLKA) = NEUGRP
    >* Kinetic energy of the particle (GeV)
    > TKEFLK (NPFLKA) = ENERGY
    >* Particle momentum
    > PMOFLK (NPFLKA) = SQRT ( TKEFLK (NPFLKA) * ( TKEFLK (NPFLKA)
    > * + TWOTWO * AM (ILOFLK(NPFLKA)) ) )
    >* Cosines (tx,ty,tz) (make sure they are properly normalized)
    > TZFLK (NPFLKA) = COSTHE
    > TYFLK (NPFLKA) = SINTHE * COSPHI
    > TXFLK (NPFLKA) = SQRT(ONEONE-TZFLK(NPFLKA)**2-TYFLK(NPFLKA)**2)
    > IF(FLRNDM(COSPHI) .LE. HLFHLF) TXFLK (NPFLKA) = -TXFLK (NPFLKA)
    >* Polarization cosines:
    > TXPOL (NPFLKA) = -TWOTWO
    > TYPOL (NPFLKA) = +ZERZER
    > TZPOL (NPFLKA) = +ZERZER
    >* Particle coordinates
    > XFLK (NPFLKA) = XBEAM
    > YFLK (NPFLKA) = YBEAM
    > ZFLK (NPFLKA) = ZBEAM
    >* Calculate the total kinetic energy of the primaries: don't change
    > IF ( ILOFLK (NPFLKA) .EQ. -2 .OR. ILOFLK (NPFLKA) .GT. 100000 )
    > & THEN
    > TKESUM = TKESUM + TKEFLK (NPFLKA) * WTFLK (NPFLKA)
    > ELSE IF ( ILOFLK (NPFLKA) .NE. 0 ) THEN
    > TKESUM = TKESUM + ( TKEFLK (NPFLKA) + AMDISC (ILOFLK(NPFLKA)) )
    > & * WTFLK (NPFLKA)
    > ELSE
    > TKESUM = TKESUM + TKEFLK (NPFLKA) * WTFLK (NPFLKA)
    > END IF
    >* Flag this is prompt radiation
    > LRADDC (NPFLKA) = .FALSE.
    > RADDLY (NPFLKA) = ZERZER
    >* Here we ask for the region number of the hitting point.
    >* NREG (NPFLKA) = ...
    >* The following line makes the starting region search much more
    >* robust if particles are starting very close to a boundary:
    > CALL GEOCRS ( TXFLK (NPFLKA), TYFLK (NPFLKA), TZFLK (NPFLKA) )
    > CALL GEOREG ( XFLK (NPFLKA), YFLK (NPFLKA), ZFLK (NPFLKA),
    > & NRGFLK(NPFLKA), IDISC )
    >* Do not change these cards:
    > CALL GEOHSM ( NHSPNT (NPFLKA), 1, -11, MLATTC )
    > NLATTC (NPFLKA) = MLATTC
    > CMPATH (NPFLKA) = ZERZER
    > CALL SOEVSV
    > RETURN
    >*=== End of subroutine Source =========================================*
    > END
    >
    >
    >------------------------------------------------------------------------
    >
    > Tot. response (p/cmq/pr) 7.1284106E-08 +/- 40.45634 %
    > Tot. response (p/cmq/pr) 4.7086409E-08 +/- 15.20756 %
    > Tot. response (p/cmq/pr) 7.0881796E-08 +/- 14.72943 %
    > Tot. response (p/cmq/pr) 3.0243843E-08 +/- 20.67986 %
    > Tot. response (p/cmq/pr) 7.8054548E-08 +/- 17.66477 %
    > Tot. response (p/cmq/pr) 5.3089430E-08 +/- 13.35963 %
    > Tot. response (p/cmq/pr) 1.5073347E-07 +/- 3.240386 %
    > Tot. response (p/cmq/pr) 2.0100804E-07 +/- 5.698418 %
    > Tot. response (p/cmq/pr) 1.7647442E-07 +/- 6.313897 %
    > Tot. response (p/cmq/pr) 7.5248231E-08 +/- 19.95282 %
    > Tot. response (p/cmq/pr) 6.5417119E-08 +/- 15.56184 %
    > Tot. response (p/cmq/pr) 7.6103909E-08 +/- 18.11728 %
    > Tot. response (p/cmq/pr) 2.1436688E-07 +/- 7.714584 %
    > Tot. response (p/cmq/pr) 1.6674799E-07 +/- 9.129667 %
    > Tot. response (p/cmq/pr) 1.5370897E-07 +/- 10.97864 %
    > Tot. response (p/cmq/pr) 1.4343240E-07 +/- 10.75971 %
    > Tot. response (p/cmq/pr) 1.3844644E-07 +/- 8.064213 %
    > Tot. response (p/cmq/pr) 1.3238373E-07 +/- 7.113124 %
    > Tot. response (p/cmq/pr) 1.1347294E-07 +/- 4.055789 %
    > Tot. response (p/cmq/pr) 1.3479203E-07 +/- 7.975585 %
    > Tot. response (p/cmq/pr) 9.9289736E-08 +/- 11.40696 %
    > Tot. response (p/cmq/pr) 2.7271266E-08 +/- 15.22840 %
    > Tot. response (p/cmq/pr) 5.1523276E-08 +/- 24.18129 %
    > Tot. response (p/cmq/pr) 2.7501205E-08 +/- 19.45839 %
    > Tot. response (p/cmq/pr) 8.6526759E-09 +/- 58.75982 %
    > Tot. response (p/cmq/pr) 8.0905922E-09 +/- 22.45524 %
    > Tot. response (p/cmq/pr) 3.1863443E-08 +/- 51.06798 %
    > Tot. response (p/cmq/pr) 1.3314099E-07 +/- 15.35631 %
    > Tot. response (p/cmq/pr) 1.0214023E-07 +/- 18.10196 %
    > Tot. response (p/cmq/pr) 1.0933596E-07 +/- 9.578140 %
    > Tot. response (p/cmq/pr) 7.5091338E-09 +/- 16.40095 %
    > Tot. response (p/cmq/pr) 1.1314577E-08 +/- 23.84999 %
    > Tot. response (p/cmq/pr) 2.0866164E-09 +/- 33.64861 %
    > Tot. response (p/cmq/pr) 3.0393661E-08 +/- 12.93056 %
    > Tot. response (p/cmq/pr) 2.0285617E-08 +/- 19.32285 %
    > Tot. response (p/cmq/pr) 3.0324259E-08 +/- 12.96468 %
    >
    >


  • Next message: Stefan Roesler: "Re: Neutron problem"

    This archive was generated by hypermail 2.1.6 : Wed Sep 20 2006 - 13:29:38 CEST