diff --git a/e2-estia.i b/e2-estia.i index 9e3fce4..dcb51f1 100644 --- a/e2-estia.i +++ b/e2-estia.i @@ -17,6 +17,7 @@ c tallies read file=geometry/estia_inbunker.i echo read file=geometry/estia_guides.i echo read file=geometry/estia_cave.i echo +read file=geometry/estia_universes.i echo c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++ c ------------------------------------------------------- @@ -31,10 +32,10 @@ c 700 s 117.0 -155.00 13.5 3.0 703 1 s 1082 -25.0 0 5.0 704 1 s 1082 -0.25 0 5.0 705 1 s 1165 0.7 0 5.0 -706 1 s 1328 5 0 12.0 +706 1 s 1320 5 0 12.0 707 1 s 2190 3.5 0 12.0 -708 s 730 -1380 0 30.0 -709 s 1020 -1170 0 30.0 +708 s 746 -1370 0 30.0 +709 s 1030 -1155 0 30.0 710 s 1328 -1756 -13 12.0 c Estia geometry read file=geometry/estia_surfaces.i echo @@ -95,6 +96,10 @@ c ------------------- TALLY CARDS --------------------------- c c read file=tallies/bunker_mesh.i echo c read file=tallies/bunker_mesh_p.i echo +c read file=tallies/bunker_mesh_np.i echo +c read file=tallies/full_mesh.i echo +c read file=tallies/full_mesh_p.i echo +c read file=tallies/full_mesh_np.i echo c read file=tallies/inbunker_f4.i echo c read file=tallies/chopper_pit_doses.i echo c read file=tallies/bunker_wall_f4.i echo @@ -122,7 +127,7 @@ ext:n 0 0 0.5V2 0 0.5V2 0 8r 0.5V2 7r 0 0.5V2 0 0.5V2 8r 0 0.5V2 0 0.5V2 0 0 0.5V2 0 0.5V2 0 0 0.5V2 0 0.5V2 0 0.5V2 0 0.5V2 0 22r 0.5V2 10r 0 4r 0.5V2 3r 0 0.5V2 13r 0 48r -0.5V1 2r 0 4r -0.1V1 9r 0 0.1V3 3r 0 11r 0.1V3 2r - 0 0.1V3 0 0.1V3 0 0.1V3 5r 0 30r + 0 0.1V3 0 0 0.1V3 0 0 0.1V3 5r 0 42r c exponential transform reference vectors: origin, selene 1, inst. shutter vect v1 0 0 0 v2 1000 -1394 -6.5 v3 1282 -1793 -13 c @@ -137,9 +142,9 @@ mode n p n h / d t s a c volume calculated for tally cells (spheres) c and set for chopper disk/spindle and stages x1/x2/z1/z2/rot vol 1 507r 10j 1 107r 3848.5 10104 1 13r 479.5 479.5 913 913 1734 - 1 18r 3.34e6 1 2r 3.5e6 1 3.5e6 1 6r 2.05e5 1 12r $ wall segments -imp:h,/,d,t,s,a 0 0 1 691r -imp:n,p 0 0 1 691r + 1 18r 3.34e6 1 2r 3.5e6 1 3.5e6 1 6r 2.05e5 1 26r $ wall segments +imp:h,/,d,t,s,a 0 0 1 705r +imp:n,p 0 0 1 705r cut:h,/,d,t,s,a 1e+08 1e-3 c cut:n j 0.99E-6 cut:n j 1e-11 @@ -158,6 +163,7 @@ c * McStas source terms for gamma doeses c read file=sources/mcstas_source_vs.i echo c read file=sources/mcstas_source_sample.i echo c read file=sources/mcstas_gammas.i echo +c read file=sources/mcstas_gammas_2.i echo prdmp 5e7 -60 1 2 5e6 print 30 32 $ 120 128 150 160 161 162 c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++ diff --git a/geometry/ess_cells.i b/geometry/ess_cells.i index 9cfb2ee..f72817d 100644 --- a/geometry/ess_cells.i +++ b/geometry/ess_cells.i @@ -5,15 +5,10 @@ c Main ESS plant geometry without the instrument 159 ) -1 ( -154 : 161 : 166 ) ( -62 : 64 : 63 ) ( -3 : -103 : -70 : -64 : 105 : 104 : 72 ) ( -121 : -65 : -63 : 122 : 101 : 124 ) - #((68 -72 8209 -70 -104 -2) - :(72 -8101 8201 8202 8203 8209 -2) - :(8709 8710 8701 8703 8706 8707 8708 8709 8710 (8704 : 8705)) - :(8709 8710 8702 8703 8711 8712) - :(8709 8809 -8311 -8712 -8701 -8313) - :(8101 -8709 8301 8302 8303 8209) 3) 600 #707 - #(72 -8709 -104 -8201 8203 8204) - #(((112 -8202):(8202 -8302 8101)) 72 -8709 8203 8204) #(8709 8720 112 -104 8703 8701) + #(-104 112 72 8710 8209 8203) + #( 104 8708 8707 8710 8703 8701) + #(68 -72 8209 -70 -104 -2) 3 118 0.122325 24 4 -19 4 0 159 -137 -21 153 -154 5 118 0.122325 37 -5 -19 @@ -496,7 +491,7 @@ c Main ESS plant geometry without the instrument -112 257 200 -3.9 ( -3 : -69 : -2 : -70 : -64 : 71 : 68 ) 3 -102 2 103 66 -67 64 -258 49 0.074981 64 -67 66 -104 101 -72 3 +258 49 0.074981 -104 101 -67 ((64 66 -72 3):(72 -610 59)) 259 49 0.074981 ( -3 : -69 : -2 : -70 : -64 : 71 : 72 ) 3 -72 103 -104 -66 70 64 260 200 -3.9 -67 66 -107 -2 -72 68 3 diff --git a/geometry/ess_surfaces.i b/geometry/ess_surfaces.i index a9574c4..1b87f39 100644 --- a/geometry/ess_surfaces.i +++ b/geometry/ess_surfaces.i @@ -100,10 +100,13 @@ c 70 pz -334.4 98 pz 341.389474 99 pz 351.126316 100 pz 360.863158 -101 p -0.906307787 -0.422618262 0.0 -5.52451529e-14 +c Building Wall E01 +c 101 p -0.906307787 -0.422618262 0.0 -5.52451529e-14 +101 p -0.866025 -0.5 0.0 80.345 102 cz 2450 103 p 0.0 -1 0.0 -80 -104 p -0.906307787 -0.422618262 0.0 80 +c 104 p -0.906307787 -0.422618262 0.0 80 +104 p -0.866025 -0.5 0.0 180.345 105 pz 375.6 106 cz 2770 107 p -0.125718017 -0.992066016 0.0 0.0 @@ -379,3 +382,6 @@ c 463 pz 13.5 572 p 0.809016994 0.587785252 0.0 17.5262 575 pz 7.4 600 s 1603.5 2265.0 13.0 3.0 +c additional surfaces from CATIA +c D01/D02 boundary +610 px 1090 diff --git a/geometry/estia_cave.i b/geometry/estia_cave.i index dafdb56..d32f953 100644 --- a/geometry/estia_cave.i +++ b/geometry/estia_cave.i @@ -1,7 +1,7 @@ c Experimental cave c First largest area: up-stream/down-stream/roof/floor/walls c then cut corner and add extended roof -8701 49 0.074981 ((8709 8710 8701 8703 8706 8707 8708 8709 8710 (8704 : 8705)) +8701 4 -1.2E-3 ((8709 8710 8701 8703 8706 8707 8708 8709 8710 (8704 : 8705)) :(8709 8710 8702 8703 8711 8712)) #((-8809 -8810 -8801 -8803 -8806 -8807 -8808 -8809 -8810 (-8804 : -8805)) @@ -13,10 +13,15 @@ c then cut corner and add extended roof (-8804 : -8805)) :(-8809 -8810 -8802 -8803 -8811 -8812) :(-8311 -8312 8709 8809 -8313 -8803)) 3 + #(8621 8622 8623 8624 8626) #6111 #8501 #8511 #9001 c sample position 9001 19 -8.65 -9001 +c 9001 4 -1.2E-3 -9001 c area around cave -9900 0 8709 8710 -104 (-8708:-8707) 8703 8701 -9901 0 8709 8710 112 ((-8704 -8705):-8706) 8703 8701 -9902 0 8720 -8710 112 -104 8703 8701 +9900 4 -1.2E-3 8709 8710 -104 (-8708:-8707) 8703 8701 +9901 4 -1.2E-3 8709 8710 112 ((-8704 -8705):-8706) 8703 8701 +9902 4 -1.2E-3 8720 -8710 112 -104 8703 8701 +c visualize Selene gamma sources (comment for real runs) +c 50001 0 -50001 +c 50002 0 -50002 diff --git a/geometry/estia_guides.i b/geometry/estia_guides.i index 91e3b5a..f8317b5 100644 --- a/geometry/estia_guides.i +++ b/geometry/estia_guides.i @@ -13,19 +13,23 @@ c -----Wall Cut-Out--- 6003 4 -1.2E-3 (6103 -72 -6601 -6602 -6603 -6604 -6605) 3 #(7103 8405 -8921) #(7102 -7103 7205 7206 7207 7208) - #(-8405 8401 8402 8403 8404) + #(8421 8422 8423 8424 8425 8426) +c #(-8405 8401 8402 8403 8404) c Vacuum before and inside Selene Guide 1 -6101 0 ((7103 8405 -8922 3) - :(-8405 8406 -8411 8412 8413 8402 3)) - #706 #8311 #8411 #8412 #8413 -6102 5 0.0582256 (7103 8405 8922 -8921 3) -c Vaccum after Selene 1, until and inside Selen 2 -6111 0 ((-8406 -8100 -8932 3) +6101 0 (7103 8435 -8922 3) #706 +6102 5 0.0582256 (7103 -8425 8922 -8921 3) +6103 0 (-8431 -8432 -8433 -8434 -8435 -8436 (-8441:-8442) (-8443:-8444)) + #8301 #8311 #8411 #8412 #8413 +c Vacuum after Selene 1, until Selen 2 +6111 0 ((8436 -8100 -8932 3) :(8100 -8101 -8501 -8502 -8503 -8504 3) - :(8101 -8121 -8932 3) - :(8121 -8122 8601 -8611 -8612 -8613 3)) - #707 #8511 -6112 5 0.0582256 (-8406 -8121 8932 -8931 3 #(8100 -8101)) + :(8101 8635 -8932 3)) + #707 +6112 5 0.0582256 (-8426 -8625 8932 -8931 3 #(8100 -8101)) +c Vacuum inside Selene 2 +6113 0 ((-8631 -8632 -8633 -8634 -8635 -8636 + #(8121 -8122 8601 8602 8603 8604)) + :(8121 -8122 8601 -8611 -8612 -8613 3)) #8511 c -----Bunker Wall Collimator (BWC)----- 7001 0 (7100 -7103 -7301 -7302 -7303 -7304 3) #(7101 -7401 7402) @@ -48,38 +52,35 @@ c ---- Outside of Bunker c -- Shielding Walls -- c concrete floor 8001 49 0.074981 ((68 -72 8209 -70 -104 -2) - :(72 -8101 8209 -59 8201 8202) - :(8101 -8709 8209 -59 8301 8302)) 3 + :(72 8710 8707 112 8209 (-8304:(-8305 -610 -101)))) + 3 -104 c ceiling outside -8002 49 0.074981 ( 72 -8101 8207 8203 8201 8202) +8002 4 -1.2E-3 ( 72 -8101 8207 8203 8201 8202) c left outer wall (E01) -8003 4 -1.2E-3 ( 72 -8100 8201 -8205 59 -8207 3) +8003 4 -1.2E-3 ( 72 -8100 8201 -8205 59 -8207 3) #(-8305 -610) c Selene 1/2 separator 8004 200 -3.9 (8100 -8101 8201 8202 59 -8207 3) #(8100 -8101 -8221 -8222 -8223 -8224) c Optical cave -8005 200 -3.9 (8101 -8709 59 8301 8302 8303 3) +8005 4 -1.2E-3 (8101 -8709 59 8201 8302 8203 3) #(8101 -8709 -8311 -8312 -8313 -8314) c air 8008 4 -1.2E-3 (72 -8100 -8211 -8212 -8213 -8214 3) - #(-8406 -8931) - #(-8405 8406 8401 8402 8403 8404) -c #8370 #8371 #8372 + #(8421 8422 8423 8424 8425 8426) + #(-8426 -8931) #(-8305 -610) 8009 4 -1.2E-3 (8101 -8709 -8311 -8312 -8313 -8314 3) #(8121 -8122 8601 8602 8603 8604) #(-8121 -8931) + #(8621 8622 8623 8624 8625 8626) #707 #710 -c B4C concrete inner layers -8011 4 -1.2E-3 (8211 8215 -8213 -8214 72 -8100 3) -8012 16 -2.4 (8205 8202 8213 -8218 72 -8100 3) -8013 4 -1.2E-3 ((-8216 8111 -8100):(72 -8111 18202)) - 8212 -8213 -8214 3 -8014 16 -2.4 (8202 -8206 8217 -8213 -8214 8111 -8100 3) -c Steel middle layers -8020 71 0.0843223 (-18202 8202 -8213 -8214 72 -8111 3) -8021 71 0.0843223 (-8215 8205 -8213 -8214 72 -8100 3) -8022 71 0.0843223 (8205 8202 8218 -8207 72 -8100 3) -8023 71 0.0843223 (8216 8202 (-8217:8206) -8213 -8214 8111 -8100 3) +c Selene 1 shielding +8011 4 -1.2E-3 (8211 8205 -8213 -8214 72 -8100 3) #(-8305 -610) +8012 16 -2.400 (8205 8202 8213 -8218 72 -8100 3) +8013 4 -1.2E-3 72 8202 8212 8217 -8213 -8214 3 #(-8305 -610) +c +8020 71 0.0843223 (-8217 -8100 8215 8202 -8213 -8214 3) +8021 71 0.0843223 (-8217 -8100 -8215 8212 -8213 -8214 3) +8022 71 0.0843223 (8205 8202 8218 -8207 72 -8100 3) c c Selene 1 carrier 8301 18 -2.4 -8405 8406 ((8401 8402 8403 8404) @@ -87,6 +88,14 @@ c Selene 1 carrier c Selene 1 glass (change material!!) 8311 12 -2.51 -8405 8406 -8903 -8904 -8518 #(-8901 -8902) #(-8515 -8517) 3 +c Selene 1 Vacuum tank (stainless steel) +8321 3 0.0876276 ( 8421 8422 8423 8424 8425 8426) + #(-8431 -8432 -8433 -8434 -8435 -8436) + #(8425 8435 -8922) + #(8426 8436 -8932) +c Lead shielding layer at top-right corner of tank +c 13 -11.34 +8422 0 -8435 -8436 -8432 -8433 ((8441 8442):(8443 8444)) c Reference sensor material for Selene 1 dose simulations c 8370 71 0.0843223 -8370 c 8371 71 0.0843223 -8371 @@ -110,9 +119,19 @@ c Selene 2 carrier c Selene 2 glass (change material!!) 8511 12 -2.51 8121 -8122 -8913 -8914 8915 #(-8911 -8912) #(-8515 -8517) 3 +c Selene 2 Vacuum tank (stainless steel) +8521 3 0.0876276 ( 8621 8622 8623 8624 8625 8626) + #(-8631 -8632 -8633 -8634 -8635 -8636) + #(8625 8635 -8932) +c #(8626 8636 -8932) c Air outside the guide shielding -8601 0 72 -8709 -104 -8201 8203 8204 +8601 4 -1.2E-3 ((72 -8709 -104 -8201 8203 8204) + :(-8701 8203 8710 8709 -104 8302)) + #(101 -610 -67) #(-8305 -610) + #(8302 8702 (8712:(8809 -8311 -8313))) #708 -8602 0 ((112 -8202) - :(8202 -8302 8101)) 72 -8709 8203 8204 +8602 4 -1.2E-3 ((112 -8202 -8709) + :(8202 -8302 8101 -8709) + :(-8701 8203 8710 8709 -8302 112)) + 72 8203 8204 #(-8305 -610) #709 diff --git a/geometry/estia_surfaces.i b/geometry/estia_surfaces.i index eb6803c..bdb70c0 100644 --- a/geometry/estia_surfaces.i +++ b/geometry/estia_surfaces.i @@ -471,7 +471,7 @@ c Sides 4 c Sides 5 6601 p -0.809017 -0.587785 0.000000 85.500 6602 p 0.809017 0.587785 0.000000 41.800 -6603 pz 76.000 +6603 pz 85.0 6604 p 0.000000 0.000000 -1.000000 81.000 6605 p -0.987688 0.156434 0.000000 -895.551 c Cu-collimator @@ -508,8 +508,8 @@ c Wall radii 8100 cz 2112.01 8101 cz 2157.17 c 8102 cz 3023.41 -8103 cz 1355.64 -8104 cz 2075.58 +c 8103 cz 1355.64 +c 8104 cz 2075.58 8111 cz 1655.0 8112 cz 1665.0 8113 cz 1710.0 @@ -521,9 +521,9 @@ c 8102 cz 3023.41 c Shielding wall geometry c first section outside 8201 p 0.813101 0.582123 0.000000 -148.964 -8202 p -0.776046 -0.630676 0.000000 1.917 -18202 p -0.776046 -0.630676 0.000000 16.917 -8203 p 0.000000 0.000000 -1.000000 -166.400 +8202 p -0.771625 -0.636078 0.000000 -21.127 +c 18202 p -0.776046 -0.630676 0.000000 16.917 +8203 p 0.000000 0.000000 -1.000000 -666.400 8204 pz -200.000 c steel/concrete interface 8205 p 0.813101 0.582123 0.000000 -98.964 @@ -533,14 +533,14 @@ c 8208 p -0.776046 -0.630676 0.000000 31.917 8209 pz -300.000 c inside 8211 p -0.813101 -0.582123 0.000000 68.794 -8212 p 0.813101 0.582123 0.000000 71.037 +8212 p 0.790690 0.612217 0.000000 0.000 8213 pz 85.900 8214 p 0.000000 0.000000 -1.000000 200.000 c B4C concrete/steel interface -8215 p 0.813101 0.582123 0.000000 -88.964 -8216 p 0.813101 0.582123 0.000000 91.037 -8217 p 0.813101 0.582123 0.000000 101.037 -8218 pz 105.900 +8215 p 0.739506 0.673150 0.000000 -102.437 +8216 pz 50.000 +8217 p -0.612217 0.790690 0.000000 -1958.974 +8218 pz 105.900 c 8219 p -0.776046 -0.630676 0.000000 111.917 c wall breakthrough 8221 p -0.813101 -0.582123 0.000000 32.866 @@ -551,8 +551,10 @@ c Optical Cave Part c Outside 8301 p 0.813101 0.582123 0.000000 -147.866 8302 p -0.813101 -0.582123 0.000000 -102.134 -8303 p 0.000000 0.000000 -1.000000 -132.500 +8303 p 0.000000 0.000000 -1.000000 -666 +c -132.500 8304 pz -200.000 +8305 pz -150.000 c Inside 8311 p -0.813101 -0.582123 0.000000 117.866 8312 p 0.813101 0.582123 0.000000 72.134 @@ -578,6 +580,26 @@ c Inside 8411 p -0.813101 -0.582123 0.000000 14.866 8412 p -0.007112 0.009934 -0.999925 -38.883 8413 p 0.007112 -0.009934 0.999925 -11.117 +c Selene 1 Vacuum +c Outside +8421 p 0.813101 0.582123 0.000000 -66.866 +8422 p -0.813101 -0.582123 0.000000 -35.134 +8423 p -0.007112 0.009934 -0.999925 -96.883 +8424 p 0.007112 -0.009934 0.999925 -52.617 +8425 p 0.582080 -0.813040 -0.012217 1345.888 +8426 p -0.582080 0.813040 0.012217 -2097.888 +c Inside +8431 p -0.813101 -0.582123 0.000000 64.166 +8432 p 0.813101 0.582123 0.000000 31.934 +8433 p 0.007112 -0.009934 0.999925 94.183 +8434 p -0.007112 0.009934 -0.999925 50.917 +8435 p -0.582080 0.813040 0.012217 -1347.388 +8436 p 0.582080 -0.813040 -0.012217 2096.388 +c Anti-gamma lead shield concept +8441 p 0.813101 0.582123 0.000000 28.934 +8442 p 0.007112 -0.009934 0.999925 44.183 +8443 p 0.007112 -0.009934 0.999925 91.183 +8444 p 0.813101 0.582123 0.000000 -9.934 c Beam through wall c Inside 8501 p -0.792291 -0.610143 -0.000426 79.405 @@ -597,19 +619,39 @@ c 8516 p -0.007112 0.009934 -0.999925 -5.400 8519 p 0.813101 0.582123 0.000000 -9.866 c Selene 2 Carrier c Outside -8601 p 0.813101 0.582123 0.000000 -1.866 -8602 p -0.813101 -0.582123 0.000000 -58.134 -8603 p -0.007112 0.009934 -0.999925 -68.883 -8604 p 0.007112 -0.009934 0.999925 -41.117 +8601 p 0.813101 0.582123 0.000000 -2.7658 +8602 p -0.813101 -0.582123 0.000000 -42.9342 +8603 p -0.007112 0.009934 -0.999925 -58.883 +8604 p 0.007112 -0.009934 0.999925 -31.117 c Inside -8611 p 0.813101 0.582123 0.000000 19.634 -8612 p 0.007112 -0.009934 0.999925 31.383 -8613 p -0.007112 0.009934 -0.999925 1.617 -c Experimental Cave +8611 p 0.813101 0.582123 0.000000 19.1342 +8612 p 0.007112 -0.009934 0.999925 38.8829 +8613 p -0.007112 0.009934 -0.999925 11.1171 +c Guide field plate +8616 p -0.813101 -0.582123 0.000000 -0.1342 +8617 p 0.813101 0.582123 0.000000 0.8658 +8618 p -0.007112 0.009934 -0.999925 -30.5829 +8619 p 0.007112 -0.009934 0.999925 -2.8171 +c Selene 2 Vacuum c Outside +8621 p 0.813101 0.582123 0.000000 -30.866 +8622 p -0.813101 -0.582123 0.000000 -71.134 +8623 p -0.007112 0.009934 -0.999925 -96.883 +8624 p 0.007112 -0.009934 0.999925 -52.617 +8625 p 0.582080 -0.813040 -0.012217 2545.388 +8626 p -0.582080 0.813040 0.012217 -3297.388 +c Inside +8631 p -0.813101 -0.582123 0.000000 28.166 +8632 p 0.813101 0.582123 0.000000 67.934 +8633 p 0.007112 -0.009934 0.999925 94.183 +8634 p -0.007112 0.009934 -0.999925 50.117 +8635 p -0.582080 0.813040 0.012217 -2546.888 +8636 p 0.582080 -0.813040 -0.012217 3295.888 +c Experimental Cave +c Outside 8701 p 0.000000 0.000000 -1.000000 -87.940 8702 p 0.000000 0.000000 -1.000000 -136.940 -8703 pz -300.000 +8703 pz -200.000 8704 p -0.813101 -0.582123 0.000000 -142.134 8705 p 0.582123 -0.813101 0.000000 3364.795 8706 p -0.736097 -0.676876 0.000000 196.513 @@ -620,7 +662,7 @@ c Outside 8711 p -0.813101 -0.582123 0.000000 -102.134 8712 p 0.813101 0.582123 0.000000 -77.866 8720 p -0.582123 0.813101 0.000000 -4464.795 -c Inside +c Inside 8801 pz 57.941 8802 pz 108.941 8803 p 0.000000 0.000000 -1.000000 166.960 @@ -652,4 +694,23 @@ c Vacuum vessels before S1/S1-S2 8932 3 c/x 0 -2 15.0 c Sphere at sample position to simulate gamma generation 9001 15 s 0 0 0 12.0 +c UNIVERSES +c Selene 1 +18301 3 px -295 +18302 3 px -247 +18311 3 c/y -291.25 1.75 1.8 +18312 3 c/y -291.25 12.35 1.8 +18313 3 c/y -291.25 18.75 1.8 +18314 3 c/y -291.25 -1.75 1.8 +18315 3 c/y -291.25 -12.35 1.8 +18316 3 c/y -291.25 -18.75 1.8 +18321 3 c/y -251.75 22.25 1.8 +18322 3 c/y -251.75 12.35 1.8 +18323 3 c/y -251.75 18.75 1.8 +18324 3 c/y -251.75 -22.25 1.8 +18325 3 c/y -251.75 -12.35 1.8 +18326 3 c/y -251.75 -18.75 1.8 +c visualize Selene gamma sources +50001 1 rpp 1339.25 2060.75 2 8 2 8 +50002 1 rpp 2539.25 3260.75 -8 -2 -8 -2 c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++ diff --git a/geometry/estia_universes.i b/geometry/estia_universes.i new file mode 100644 index 0000000..0b06a7f --- /dev/null +++ b/geometry/estia_universes.i @@ -0,0 +1,17 @@ +c UNIVERSES for Selene geometry lattices +c Selene 1 holes for adjusters +18300 0 18301 -18302 U=1 fill=2 lat=1 +18301 18 -2.4 18311 18312 18313 18314 18315 18316 + 18321 18322 18323 18324 18325 18326 U=2 +18311 0 -18311:-18314 U=2 +18312 0 -18312 U=2 +18313 0 -18313 U=2 +c 18314 0 -18314 U=2 +18315 0 -18315 U=2 +18316 0 -18316 U=2 +18321 0 -18321:-18323 U=2 +18322 0 -18322 U=2 +c 18323 0 -18323 U=2 +18324 0 -18324:-18326 U=2 +18325 0 -18325 U=2 +c 18326 0 -18326 U=2 diff --git a/plot_tally.py b/plot_tally.py index b695f50..63e0303 100644 --- a/plot_tally.py +++ b/plot_tally.py @@ -13,7 +13,7 @@ from pylab import * from matplotlib.colors import LogNorm COLORMAP='tab20c' -SCALING=1.56e16/8.7374e15 +SCALING=1.0 def read_tally(fname, use_tally=None): if fname.endswith('.gz'): @@ -68,9 +68,16 @@ def plot_xy(outfile, res, zidx=5, Imin=1e-6, Imax=1e4, cticks=None): #fig.savefig(outfile) #pc=p.get_facecolors() #pc[:,3]=where(I[:,:,zidx].T.flatten()>0, 0.9-0.25*dI[:,:,zidx].T.flatten(), 0.) - ax.add_patch(Circle((0,0), radius=5.50, color='black', fill=False, lw=2)) - ax.add_patch(Circle((0,0), radius=11.50, color='black', fill=False, lw=2)) - ax.add_patch(Circle((0,0), radius=15.00, color='black', fill=False, lw=2)) + ax.add_patch(Circle((-0.151,0.021), radius=5.50, color='black', fill=False, lw=2)) + ax.add_patch(Circle((-0.151,0.021), radius=11.50, color='black', fill=False, lw=2)) + ax.add_patch(Circle((-0.151,0.021), radius=15.00, color='black', fill=False, lw=2)) + # shielding walls outside bunker + ax.add_line(Line2D([21.12-0.151, 21.12-0.151], [-1.5,1.5], color='black', lw=2)) + ax.add_line(Line2D([14.845,19.426], [-0.541,-0.713], color='black', lw=2)) + ax.add_line(Line2D([19.426,19.405], [-0.713,-1.248], color='black', lw=2)) + ax.add_line(Line2D([19.405,21.205], [-1.248,-1.460], color='black', lw=2)) + ax.add_line(Line2D([14.723, 22.740], [2.232, 3.084], color='black', lw=2)) + ax.add_line(Line2D([14.801, 20.894], [-1.210, -1.724], color='black', lw=2)) if outfile is not None: fig.savefig(outfile, transparent=False, dpi=300) @@ -89,9 +96,12 @@ def plot_xz(outfile, res, yidx=5, Imin=1e-6, Imax=1e4, cticks=None): #fig.savefig(outfile) #pc=p.get_facecolors() #pc[:,3]=where(I[:,yidx].T.flatten()>0, 0.9-0.25*dI[:,yidx].T.flatten(), 0.) - ax.plot([sqrt(5.5**2-ypos**2), sqrt(5.5**2-ypos**2)], [z.min()/100., z.max()/100.], color='black', lw=2) - ax.plot([sqrt(11.5**2-ypos**2), sqrt(11.5**2-ypos**2)], [z.min()/100., z.max()/100.], color='black', lw=2) - ax.plot([sqrt(15.0**2-ypos**2), sqrt(15.0**2-ypos**2)], [z.min()/100., z.max()/100.], color='black', lw=2) + ax.plot([sqrt(5.5**2-ypos**2)-0.151, sqrt(5.5**2-ypos**2)-0.151], + [z.min()/100., z.max()/100.], color='black', lw=2) + ax.plot([sqrt(11.5**2-ypos**2)-0.151, sqrt(11.5**2-ypos**2)-0.151], + [z.min()/100., z.max()/100.], color='black', lw=2) + ax.plot([sqrt(15.0**2-ypos**2)-0.151, sqrt(15.0**2-ypos**2)-0.151], + [z.min()/100., z.max()/100.], color='black', lw=2) if outfile is not None: fig.savefig(outfile, transparent=False, dpi=300) @@ -112,26 +122,67 @@ def plot_yz(outfile, res, xidx=5, Imin=1e-6, Imax=1e4, cticks=None): if outfile is not None: fig.savefig(outfile, transparent=False, dpi=300) +def get_name(fname): + # extract tag-name and version info from filename + stripped=os.path.basename(fname)[9:].split('.')[0] + tag, version=stripped.split('-', 1) + return tag, version + + if __name__=='__main__': import sys, os - res=read_tally(sys.argv[1], use_tally=1) + if '-o' in sys.argv: + idx=sys.argv.index('-o') + sys.argv.pop(idx) + prefix=sys.argv.pop(idx) + else: + prefix='estia_dosemap' + print sys.argv[1] + if ':' in sys.argv[1]: + fname, tally=sys.argv[1].split(':',2) + tally=int(tally) + else: + fname=sys.argv[1] + tally=None + res=read_tally(fname, use_tally=tally) x,y,z,I,dI=res I*=SCALING - fpre=os.path.join('.', 'images', os.path.basename(sys.argv[1])[:-4]) + fpath=os.path.join('.', 'images', get_name(fname)[0]+'-'+get_name(fname)[1]) + if tally is not None: + fpath+='[%s]'%tally + if len(sys.argv)>2: + for fi in sys.argv[2:]: + print fi + if ':' in fi: + fi, tally=fi.split(':',2) + tally=int(tally) + else: + tally=None + resi=read_tally(fi, use_tally=tally) + xi,yi,zi,Ii,dIi=resi + Ii*=SCALING + I+=Ii + dI=sqrt(dI**2+dIi**2) + res=(x,y,z,I,dI) + fpath+='+'+get_name(fi)[1] + if tally is not None: + fpath+='[%s]'%tally + os.mkdir(fpath) + fpre=os.path.join(fpath, prefix) xi, yi, zi=res[3].shape for zidx in range(zi): print fpre+'_xy_%02i.png'%zidx figure(figsize=(12,8)) plot_xy(fpre+'_xy_%02i.png'%zidx, res, zidx=zidx, - Imin=1e-8, Imax=1e2, cticks=[1e-8, 1e-6, 1e-4, 1e-2, 1, 100]) + Imin=1e-7, Imax=1e-2, cticks=[1e-7, 1e-6, 1e-5, 1e-4, 1e-3, 1e-2]) for yidx in range(yi): print fpre+'_xz_%02i.png'%yidx figure(figsize=(12,8)) plot_xz(fpre+'_xz_%02i.png'%yidx, res, yidx=yidx, - Imin=1e-8, Imax=1e2, cticks=[1e-8, 1e-6, 1e-4, 1e-2, 1, 100]) + Imin=1e-7, Imax=1e-2, cticks=[1e-7, 1e-6, 1e-5, 1e-4, 1e-3, 1e-2]) for xidx in range(xi): print fpre+'_yz_%02i.png'%xidx figure(figsize=(8,10)) plot_yz(fpre+'_yz_%02i.png'%xidx, res, xidx=xidx, - Imin=1e-8, Imax=1e2, cticks=[1e-8, 1e-6, 1e-4, 1e-2, 1, 100]) + Imin=1e-7, Imax=1e-2, cticks=[1e-7, 1e-6, 1e-5, 1e-4, 1e-3, 1e-2]) diff --git a/sources/e2_2m.i b/sources/e2_2m.i index fc22312..d84ab58 100644 --- a/sources/e2_2m.i +++ b/sources/e2_2m.i @@ -1,10 +1,10 @@ c source definition based on tally_181.txt c but adapted to model geometry and actual size of Estia beam extraction c make sure to supply the correct WWINP file -WWP:n j j 5 j -1 j 1.0 +WWP:n j j 5 j -1 j 1.79 c the source definition with neutrons/s sdef dir=d1 erg=fdir d2 x=199.98 y=d23 z=d24 vec=1 0 0 tr=101 - par=n ara=165.0 wgt=3.6041775E13 + par=n ara=165.0 wgt=6.436556E13 si1 -1 0 0.5 0.866 0.9659 0.9848 0.99619 0.9975640 0.998629 0.999390 0.999847 1 sp1 0 5.01272E+12 1.04029E+12 2.01432E+12 1.39792E+12 diff --git a/sources/mcstas_gammas.i b/sources/mcstas_gammas.i index 802d2d1..a56b283 100644 --- a/sources/mcstas_gammas.i +++ b/sources/mcstas_gammas.i @@ -1,145 +1,144 @@ -c source definition based on tally_181.txt -c but adapted to model geometry and actual size of Estia beam extraction +c source definition based on McStas c make sure to supply the correct WWINP file WWP:n j j 5 j -1 j 1.0 -c the source definition with neutrons/s +c the source definition with gammas/s sdef erg=fx d1 x=d31 y=d32 z=d33 tr=1 par=p wgt=9.59833e+09 ds1 q 1339.25 100 - 1387.35 101 1435.45 102 1483.55 103 1531.65 104 1579.75 105 - 1627.85 106 1675.95 107 1724.05 108 1772.15 109 1820.25 110 + 1387.35 101 1435.45 102 1483.55 103 1531.65 104 1579.75 105 + 1627.85 106 1675.95 107 1724.05 108 1772.15 109 1820.25 110 1868.35 111 1916.45 112 1964.55 113 2012.65 114 2060.75 115 c Energy distribution for all positions along the guide si100 0 1 sp100 0 1e-20 -si101 1e-12 0.15 0.2 0.3 0.4 - 0.5 0.6 0.8 1.0 1.5 - 2.0 3.0 4.0 5.0 6.0 +si101 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 sp101 0.00000E+00 1.56557E+05 2.76591E+05 8.82267E+05 3.79148E+06 7.22100E+09 7.15782E+04 2.36360E+05 1.64987E+06 1.20477E+07 5.83592E+06 5.13734E+06 1.05131E+07 9.13771E+06 2.23406E+06 5.55688E+06 3.73177E+06 9.84371E+06 3.67345E+04 5.62062E+04 C ------------------------------------------------------------------------ -si102 1e-12 0.15 0.2 0.3 0.4 - 0.5 0.6 0.8 1.0 1.5 - 2.0 3.0 4.0 5.0 6.0 +si102 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 sp102 0.00000E+00 1.75066E+05 3.14271E+05 1.00187E+06 4.19715E+06 4.29576E+08 6.49313E+04 8.35318E+04 1.59407E+06 1.13391E+07 3.24825E+06 2.40414E+06 3.23346E+06 2.35587E+06 1.81032E+06 4.22018E+06 2.81568E+06 1.08595E+07 3.50283E+04 1.16665E+04 C ------------------------------------------------------------------------ -si103 1e-12 0.15 0.2 0.3 0.4 - 0.5 0.6 0.8 1.0 1.5 - 2.0 3.0 4.0 5.0 6.0 +si103 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 sp103 0.00000E+00 1.14282E+05 2.09752E+05 6.68141E+05 2.71818E+06 2.88507E+08 4.31570E+04 5.55428E+04 1.05765E+06 7.30837E+06 2.10495E+06 1.58164E+06 2.11537E+06 1.54343E+06 1.20473E+06 2.79310E+06 1.87709E+06 7.24741E+06 2.27791E+04 7.70741E+03 C ------------------------------------------------------------------------ -si104 1e-12 0.15 0.2 0.3 0.4 - 0.5 0.6 0.8 1.0 1.5 - 2.0 3.0 4.0 5.0 6.0 +si104 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 sp104 0.00000E+00 1.19942E+05 2.24594E+05 7.14916E+05 2.83168E+06 2.69546E+08 4.59551E+04 5.82997E+04 1.12511E+06 7.56863E+06 2.17669E+06 1.65382E+06 2.17880E+06 1.58300E+06 1.28282E+06 2.95750E+06 2.00035E+06 7.75812E+06 2.37907E+04 7.90258E+03 C ------------------------------------------------------------------------ -si105 1e-12 0.15 0.2 0.3 0.4 - 0.5 0.6 0.8 1.0 1.5 - 2.0 3.0 4.0 5.0 6.0 +si105 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 sp105 0.00000E+00 1.02808E+05 1.94757E+05 6.19690E+05 2.41655E+06 2.29099E+08 3.97544E+04 5.03268E+04 9.72514E+05 6.44006E+06 1.85566E+06 1.42096E+06 1.86323E+06 1.35351E+06 1.11024E+06 2.55211E+06 1.73259E+06 6.72701E+06 2.03449E+04 6.77894E+03 C ------------------------------------------------------------------------ -si106 1e-12 0.15 0.2 0.3 0.4 - 0.5 0.6 0.8 1.0 1.5 - 2.0 3.0 4.0 5.0 6.0 +si106 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 sp106 0.00000E+00 8.86042E+04 1.69642E+05 5.39580E+05 2.07421E+06 1.99597E+08 3.45601E+04 4.37451E+04 8.44765E+05 5.51349E+06 1.59281E+06 1.22897E+06 1.60647E+06 1.16764E+06 9.65712E+05 2.21413E+06 1.50825E+06 5.85933E+06 1.74992E+04 5.87135E+03 C ------------------------------------------------------------------------ -si107 1e-12 0.15 0.2 0.3 0.4 - 0.5 0.6 0.8 1.0 1.5 - 2.0 3.0 4.0 5.0 6.0 +si107 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 sp107 0.00000E+00 7.87877E+04 1.51840E+05 4.82849E+05 1.83971E+06 1.73871E+08 3.08863E+04 3.89901E+04 7.54666E+05 4.88100E+06 1.41072E+06 1.09310E+06 1.42351E+06 1.03396E+06 8.63179E+05 1.97565E+06 1.34859E+06 5.24414E+06 1.55375E+04 5.20473E+03 C ------------------------------------------------------------------------ -si108 1e-12 0.15 0.2 0.3 0.4 - 0.5 0.6 0.8 1.0 1.5 - 2.0 3.0 4.0 5.0 6.0 +si108 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 sp108 0.00000E+00 5.01840E+04 9.68964E+04 3.08110E+05 1.17096E+06 1.12987E+08 1.97075E+04 2.49206E+04 4.81428E+05 3.10577E+06 8.98765E+05 6.97579E+05 9.09029E+05 6.60776E+05 5.50891E+05 1.26041E+06 8.60883E+05 3.34661E+06 9.89468E+03 3.33164E+03 C ------------------------------------------------------------------------ -si109 1e-12 0.15 0.2 0.3 0.4 - 0.5 0.6 0.8 1.0 1.5 - 2.0 3.0 4.0 5.0 6.0 +si109 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 sp109 0.00000E+00 4.11446E+04 7.95509E+04 2.52944E+05 9.59525E+05 9.23981E+07 1.61749E+04 2.04454E+04 3.95097E+05 2.54402E+06 7.36328E+05 5.72028E+05 7.44921E+05 5.41445E+05 4.52162E+05 1.03416E+06 7.06655E+05 2.74751E+06 8.11001E+03 2.73082E+03 C ------------------------------------------------------------------------ -si110 1e-12 0.15 0.2 0.3 0.4 - 0.5 0.6 0.8 1.0 1.5 - 2.0 3.0 4.0 5.0 6.0 +si110 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 sp110 0.00000E+00 2.51204E+04 4.88245E+04 1.55217E+05 5.84610E+05 5.28251E+07 9.90982E+03 1.24431E+04 2.42028E+05 1.54694E+06 4.46975E+05 3.48156E+05 4.50560E+05 3.26731E+05 2.76963E+05 6.32447E+05 4.32872E+05 1.68610E+06 4.94348E+03 1.64531E+03 C ------------------------------------------------------------------------ -si111 1e-12 0.15 0.2 0.3 0.4 - 0.5 0.6 0.8 1.0 1.5 - 2.0 3.0 4.0 5.0 6.0 +si111 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 sp111 0.00000E+00 2.36608E+04 4.56443E+04 1.45144E+05 5.52264E+05 4.88854E+07 9.27597E+03 1.16381E+04 2.26683E+05 1.46399E+06 4.22049E+05 3.26897E+05 4.23743E+05 3.07060E+05 2.59131E+05 5.92778E+05 4.04758E+05 1.57629E+06 4.66257E+03 1.54109E+03 C ------------------------------------------------------------------------ -si112 1e-12 0.15 0.2 0.3 0.4 - 0.5 0.6 0.8 1.0 1.5 - 2.0 3.0 4.0 5.0 6.0 +si112 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 sp112 0.00000E+00 2.34986E+04 4.48276E+04 1.42601E+05 5.50860E+05 4.89730E+07 9.13069E+03 1.14783E+04 2.23305E+05 1.46457E+06 4.21412E+05 3.23890E+05 4.21794E+05 3.05694E+05 2.54959E+05 5.84877E+05 3.97945E+05 1.54818E+06 4.64119E+03 1.52926E+03 C ------------------------------------------------------------------------ -si113 1e-12 0.15 0.2 0.3 0.4 - 0.5 0.6 0.8 1.0 1.5 - 2.0 3.0 4.0 5.0 6.0 +si113 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 sp113 0.00000E+00 1.87764E+04 3.52661E+04 1.12245E+05 4.42775E+05 3.94292E+07 7.20582E+03 9.07992E+03 1.76424E+05 1.18187E+06 3.39132E+05 2.57888E+05 3.37896E+05 2.44903E+05 2.01079E+05 4.63087E+05 3.13513E+05 1.21806E+06 3.72001E+03 1.21944E+03 C ------------------------------------------------------------------------ -si114 1e-12 0.15 0.2 0.3 0.4 - 0.5 0.6 0.8 1.0 1.5 - 2.0 3.0 4.0 5.0 6.0 +si114 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 sp114 0.00000E+00 1.82897E+04 3.37346E+04 1.07439E+05 4.34219E+05 4.04983E+07 6.92221E+03 8.78390E+03 1.69672E+05 1.16466E+06 3.33775E+05 2.50952E+05 3.32061E+05 2.41077E+05 1.93078E+05 4.46799E+05 3.00726E+05 1.16534E+06 3.63777E+03 1.19682E+03 C ------------------------------------------------------------------------ -si115 1e-12 0.15 0.2 0.3 0.4 - 0.5 0.6 0.8 1.0 1.5 - 2.0 3.0 4.0 5.0 6.0 +si115 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 sp115 0.00000E+00 1.08301E+04 1.94946E+04 6.21412E+04 2.59394E+05 2.48374E+07 4.02189E+03 5.13569E+03 9.87454E+04 6.99903E+05 @@ -147,17 +146,17 @@ sp115 0.00000E+00 1.08301E+04 1.94946E+04 6.21412E+04 2.59394E+05 2.61025E+05 1.74300E+05 6.73544E+05 2.16453E+03 7.10515E+02 C ------------------------------------------------------------------------ c x-distribution -si31 1339.25 - 1387.35 1435.45 1483.55 1531.65 1579.75 - 1627.85 1675.95 1724.05 1772.15 1820.25 +si31 1339.25 + 1387.35 1435.45 1483.55 1531.65 1579.75 + 1627.85 1675.95 1724.05 1772.15 1820.25 1868.35 1916.45 1964.55 2012.65 2060.75 -sp31 0 - 7.29220E+09 4.79340E+08 3.21180E+08 3.03858E+08 2.58578E+08 - 2.25073E+08 1.96544E+08 1.27443E+08 1.04253E+08 6.00576E+07 - 5.56826E+07 5.57082E+07 4.47933E+07 4.57106E+07 2.79109E+07 +sp31 0 + 7.29220E+09 4.79340E+08 3.21180E+08 3.03858E+08 2.58578E+08 + 2.25073E+08 1.96544E+08 1.27443E+08 1.04253E+08 6.00576E+07 + 5.56826E+07 5.57082E+07 4.47933E+07 4.57106E+07 2.79109E+07 c y-distribution -si32 -6 -1.0 +si32 2 8 sp32 0 1 c z-distribution -si33 -6 6 +si33 2 8 sp33 0 1 diff --git a/sources/mcstas_gammas_2.i b/sources/mcstas_gammas_2.i new file mode 100644 index 0000000..1dad2c3 --- /dev/null +++ b/sources/mcstas_gammas_2.i @@ -0,0 +1,162 @@ +c source definition based on McStas +c make sure to supply the correct WWINP file +WWP:n j j 5 j -1 j 1.0 +c the source definition with gammas/s +sdef erg=fx d1 x=d31 y=d32 z=d33 tr=1 par=p wgt=1.04957e+09 +ds1 q 2539.25 100 + 2587.35 101 2635.45 102 2683.55 103 2731.65 104 2779.75 105 + 2827.85 106 2875.95 107 2924.05 108 2972.15 109 3020.25 110 + 3068.35 111 3116.45 112 3164.55 113 3212.65 114 3260.75 115 +c Energy distribution for all positions along the guide +si100 0 1 +sp100 0 1e-20 +si101 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp101 0.00000E+00 6.74946E+03 1.26349E+04 4.02190E+04 1.59434E+05 + 4.37259E+07 2.64491E+03 3.95862E+03 6.43047E+04 4.33561E+05 + 1.34669E+05 1.05571E+05 1.54428E+05 1.18439E+05 7.48297E+04 + 1.73991E+05 1.17757E+05 4.37637E+05 1.36106E+03 6.36018E+02 +C ------------------------------------------------------------------------ +si102 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp102 0.00000E+00 1.17151E+04 2.23209E+04 7.10081E+04 2.74744E+05 + 1.91981E+07 4.53666E+03 5.59356E+03 1.11043E+05 7.29343E+05 + 2.07983E+05 1.59135E+05 2.04496E+05 1.47055E+05 1.26489E+05 + 2.89983E+05 1.97214E+05 7.70669E+05 2.31033E+03 7.27223E+02 +C ------------------------------------------------------------------------ +si103 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp103 0.00000E+00 1.20232E+04 2.32095E+04 7.38021E+04 2.80546E+05 + 2.04130E+07 4.70687E+03 5.81143E+03 1.15089E+05 7.42422E+05 + 2.12512E+05 1.64201E+05 2.10375E+05 1.51470E+05 1.31339E+05 + 3.00162E+05 2.04993E+05 8.01335E+05 2.36550E+03 7.53558E+02 +C ------------------------------------------------------------------------ +si104 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp104 0.00000E+00 1.75885E+04 3.43282E+04 1.09117E+05 4.08635E+05 + 3.05406E+07 6.94852E+03 8.58466E+03 1.69754E+05 1.07843E+06 + 3.09623E+05 2.41213E+05 3.08149E+05 2.22054E+05 1.94011E+05 + 4.42225E+05 3.03071E+05 1.18520E+06 3.45334E+03 1.10997E+03 +C ------------------------------------------------------------------------ +si105 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp105 0.00000E+00 1.85638E+04 3.67221E+04 1.16674E+05 4.28978E+05 + 3.21759E+07 7.41406E+03 9.14695E+03 1.80954E+05 1.12799E+06 + 3.24748E+05 2.55490E+05 3.24699E+05 2.34011E+05 2.07133E+05 + 4.70575E+05 3.23871E+05 1.26777E+06 3.63481E+03 1.17512E+03 +C ------------------------------------------------------------------------ +si106 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp106 0.00000E+00 2.21857E+04 4.43306E+04 1.40802E+05 5.10575E+05 + 3.77016E+07 8.93175E+03 1.09930E+04 2.17852E+05 1.33861E+06 + 3.85960E+05 3.05846E+05 3.86785E+05 2.78627E+05 2.49620E+05 + 5.65666E+05 3.90547E+05 1.53035E+06 4.33437E+03 1.40296E+03 +C ------------------------------------------------------------------------ +si107 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp107 0.00000E+00 3.02306E+04 6.08550E+04 1.93239E+05 6.93589E+05 + 4.94114E+07 1.22401E+04 1.50130E+04 2.98419E+05 1.81412E+06 + 5.23227E+05 4.16740E+05 5.24421E+05 3.77374E+05 3.42124E+05 + 7.73787E+05 5.35477E+05 2.10065E+06 5.89540E+03 1.90221E+03 +C ------------------------------------------------------------------------ +si108 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp108 0.00000E+00 2.94446E+04 5.92364E+04 1.88103E+05 6.75727E+05 + 4.79016E+07 1.19154E+04 1.46108E+04 2.90519E+05 1.76765E+06 + 5.09678E+05 4.05737E+05 5.10570E+05 3.67353E+05 3.33032E+05 + 7.53323E+05 5.21216E+05 2.04477E+06 5.74269E+03 1.85095E+03 +C ------------------------------------------------------------------------ +si109 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp109 0.00000E+00 4.35594E+04 8.73170E+04 2.77304E+05 1.00114E+06 + 7.16969E+07 1.75773E+04 2.15781E+04 4.28663E+05 2.62181E+06 + 7.55655E+05 6.00004E+05 7.56553E+05 5.44496E+05 4.91231E+05 + 1.11219E+06 7.68647E+05 3.01416E+06 8.50262E+03 2.74125E+03 +C ------------------------------------------------------------------------ +si110 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp110 0.00000E+00 5.08122E+04 1.01261E+05 3.21649E+05 1.17065E+06 + 8.44071E+07 2.04082E+04 2.50830E+04 4.97895E+05 3.07096E+06 + 8.84254E+05 6.99122E+05 8.83955E+05 6.36310E+05 5.70220E+05 + 1.29292E+06 8.91913E+05 3.49561E+06 9.93105E+03 3.19801E+03 +C ------------------------------------------------------------------------ +si111 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp111 0.00000E+00 5.54824E+04 1.09734E+05 3.48652E+05 1.28218E+06 + 9.19204E+07 2.21468E+04 2.72337E+04 5.40609E+05 3.37053E+06 + 9.68856E+05 7.61698E+05 9.65731E+05 6.95046E+05 6.18581E+05 + 1.40515E+06 9.67039E+05 3.78823E+06 1.08606E+04 3.48353E+03 +C ------------------------------------------------------------------------ +si112 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp112 0.00000E+00 5.97040E+04 1.16814E+05 3.71278E+05 1.38574E+06 + 9.88383E+07 2.36237E+04 2.90770E+04 5.77107E+05 3.65346E+06 + 1.04777E+06 8.17221E+05 1.04035E+06 7.48620E+05 6.59504E+05 + 1.50209E+06 1.03023E+06 4.03280E+06 1.17128E+04 3.73774E+03 +C ------------------------------------------------------------------------ +si113 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp113 0.00000E+00 4.85803E+04 9.36078E+04 2.97674E+05 1.13437E+06 + 7.99568E+07 1.89850E+04 2.33908E+04 4.64308E+05 3.00271E+06 + 8.58300E+05 6.62033E+05 8.47386E+05 6.09533E+05 5.29619E+05 + 1.21081E+06 8.26424E+05 3.23184E+06 9.55949E+03 3.02649E+03 +C ------------------------------------------------------------------------ +si114 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp114 0.00000E+00 5.88344E+04 1.11186E+05 3.53808E+05 1.38411E+06 + 9.77202E+07 2.26369E+04 2.79595E+04 5.54391E+05 3.68221E+06 + 1.04882E+06 7.98041E+05 1.02934E+06 7.40419E+05 6.30958E+05 + 1.44953E+06 9.83220E+05 3.83908E+06 1.16223E+04 3.65352E+03 +C ------------------------------------------------------------------------ +si115 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp115 0.00000E+00 5.61276E+04 1.03236E+05 3.28823E+05 1.33382E+06 + 9.44561E+07 2.11336E+04 2.61958E+04 5.18590E+05 3.57223E+06 + 1.01277E+06 7.56511E+05 9.86081E+05 7.09333E+05 5.88350E+05 + 1.36095E+06 9.15062E+05 3.56505E+06 1.11463E+04 3.47063E+03 +C ------------------------------------------------------------------------ +c x-distribution +si31 2539.25 + 2587.35 2635.45 2683.55 2731.65 2779.75 + 2827.85 2875.95 2924.05 2972.15 3020.25 + 3068.35 3116.45 3164.55 3212.65 3260.75 +sp31 0 + 4.57688E+07 2.25345E+07 2.38501E+07 3.55841E+07 3.75155E+07 + 4.40951E+07 5.81307E+07 5.63921E+07 8.42501E+07 9.90333E+07 + 1.07862E+08 1.15949E+08 9.38290E+07 1.14450E+08 1.10325E+08 +c y-distribution +si32 -8 -2 +sp32 0 1 +c z-distribution +si33 -8 -2 +sp33 0 1 diff --git a/sources/proton.i b/sources/proton.i index 1137f5f..d267d14 100644 --- a/sources/proton.i +++ b/sources/proton.i @@ -5,7 +5,7 @@ WWP:n,p j j 10 j -1 j 8.7374e15 $ WW was generated for source scaling 1 c energy splitting for higher neutron energies and roulette below 1MeV esplt:n 0.05 1e-6 0.2 1e-3 1 10 3 100 10 3000 c the source definition with protons/s -sdef dir=1 erg=2000 par=9 vec=0.0 1 0.0 x=d1 y=-30 z=d2 wgt=8.7374e15 +sdef dir=1 erg=2000 par=9 vec=0.0 1 0.0 x=d1 y=-30 z=d2 wgt=1.560377E16 si1 A -8 -7.5 -7 -6.5 -6 -5.5 -5 -4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 sp1 D 0 0.12109 0.23438 0.33984 0.4375 0.52734 0.60938 0.68359 0.75 diff --git a/tallies/bunker_mesh.i b/tallies/bunker_mesh.i index 72e54df..58159d5 100644 --- a/tallies/bunker_mesh.i +++ b/tallies/bunker_mesh.i @@ -2,8 +2,8 @@ c large mesh tally covering the beamline from monolith to S1/2 separator tmesh rmesh1:n dose 1 trans 1 cora1 500 10i 1000 20i 1500 16i 2300 - corb1 -305 5i -55 4i 45 5i 295 - corc1 -200 1i -100 5i 100 2i 300 + corb1 -305 5i -55 4i 45 9i 495 + corc1 -200 1i -100 5i 100 5i 600 c Official ESS dose conversion factors from n/cm²/s to Sv/h c Is described in ESS-0019931 mshmf1 0.000000000 1.00000e-08 diff --git a/tallies/bunker_mesh_np.i b/tallies/bunker_mesh_np.i new file mode 100644 index 0000000..e3db20e --- /dev/null +++ b/tallies/bunker_mesh_np.i @@ -0,0 +1,80 @@ +c large mesh tally covering the beamline from monolith to S1/2 separator +tmesh + rmesh1:n dose 1 trans 1 + cora1 500 10i 1000 20i 1500 16i 2300 + corb1 -305 5i -55 4i 45 9i 495 + corc1 -200 1i -100 5i 100 5i 600 +c Official ESS dose conversion factors from n/cm²/s to Sv/h +c Is described in ESS-0019931 + mshmf1 0.000000000 1.00000e-08 + 1.00000e-09 1.11000e-08 1.00000e-08 1.28000e-08 + 2.50000e-08 1.44000e-08 1.00000e-07 1.87000e-08 + 2.00000e-07 2.11000e-08 5.00000e-07 2.37000e-08 + 1.00000e-06 2.53000e-08 2.00000e-06 2.66000e-08 + 5.00000e-06 2.78000e-08 1.00000e-05 2.82000e-08 + 2.00000e-05 2.82000e-08 5.00000e-05 2.82000e-08 + 1.00000e-04 2.80000e-08 2.00000e-04 2.78000e-08 + 5.00000e-04 2.71000e-08 1.00000e-03 2.71000e-08 + 2.00000e-03 2.74000e-08 5.00000e-03 2.87000e-08 + 1.00000e-02 3.28000e-08 2.00000e-02 4.39000e-08 + 3.00000e-02 5.65000e-08 5.00000e-02 8.28000e-08 + 7.00000e-02 1.10000e-07 1.00000e-01 1.51000e-07 + 1.50000e-01 2.18000e-07 2.00000e-01 2.84000e-07 + 3.00000e-01 4.10000e-07 5.00000e-01 6.37000e-07 + 7.00000e-01 8.35000e-07 9.00000e-01 1.00000e-06 + 1.00000e+00 1.08000e-06 1.20000e+00 1.19000e-06 + 1.50000e+00 1.31000e-06 2.00000e+00 1.47000e-06 + 3.00000e+00 1.65000e-06 4.00000e+00 1.74000e-06 + 5.00000e+00 1.78000e-06 6.00000e+00 1.79000e-06 + 7.00000e+00 1.80000e-06 8.00000e+00 1.80000e-06 + 9.00000e+00 1.80000e-06 1.00000e+01 1.80000e-06 + 1.20000e+01 1.80000e-06 1.40000e+01 1.78000e-06 + 1.50000e+01 1.77000e-06 1.60000e+01 1.76000e-06 + 1.80000e+01 1.74000e-06 2.00000e+01 1.72000e-06 + 2.10000e+01 1.71000e-06 3.00000e+01 1.63000e-06 + 5.00000e+01 1.56000e-06 7.50000e+01 1.58000e-06 + 1.00000e+02 1.60000e-06 1.30000e+02 1.61000e-06 + 1.50000e+02 1.61000e-06 1.80000e+02 1.61000e-06 + 2.00000e+02 1.61000e-06 3.00000e+02 1.67000e-06 + 4.00000e+02 1.79000e-06 5.00000e+02 1.92000e-06 + 6.00000e+02 2.05000e-06 7.00000e+02 2.16000e-06 + 8.00000e+02 2.24000e-06 9.00000e+02 2.30000e-06 + 1.00000e+03 2.35000e-06 2.00000e+03 2.76000e-06 + 5.00000e+03 3.64000e-06 1.00000e+04 4.75000e-06 +c + rmesh11:p dose 2 trans 1 + cora11 500 10i 1000 20i 1500 16i 2300 + corb11 -305 5i -55 4i 45 9i 495 + corc11 -200 1i -100 5i 100 5i 600 +c Official ESS dose conversion factors from gamma/cm²/s to Sv/h +c Is described in ESS-0019931 + mshmf2 0.000000 1.00000e-04 + 1.00E-02 2.47E-10 1.50E-02 5.62E-10 + 2.00E-02 8.1E-10 3.00E-02 1.13E-09 + 4.00E-02 1.26E-09 5.00E-02 1.33E-09 + 6.00E-02 1.4E-09 7.00E-02 1.49E-09 + 8.00E-02 1.6E-09 1.00E-01 1.87E-09 + 1.50E-01 2.69E-09 2.00E-01 3.6E-09 + 3.00E-01 5.44E-09 4.00E-01 7.2E-09 + 5.00E-01 8.89E-09 5.11E-01 9.07E-09 + 6.00E-01 1.05E-08 6.62E-01 1.14E-08 + 8.00E-01 1.34E-08 1.00E+00 1.62E-08 + 1.12E+00 1.76E-08 1.33E+00 2.01E-08 + 1.50E+00 2.2E-08 2.00E+00 2.69E-08 + 3.00E+00 3.51E-08 4.00E+00 4.21E-08 + 5.00E+00 4.82E-08 6.00E+00 5.4E-08 + 6.13E+00 5.47E-08 8.00E+00 6.7E-08 + 1.00E+01 7.92E-08 1.50E+01 1.09E-07 + 2.00E+01 1.38E-07 3.00E+01 1.85E-07 + 4.00E+01 2.23E-07 5.00E+01 2.6E-07 + 6.00E+01 2.95E-07 8.00E+01 3.52E-07 + 1.00E+02 3.96E-07 1.50E+02 4.68E-07 + 2.00E+02 5.15E-07 3.00E+02 5.8E-07 + 4.00E+02 6.19E-07 5.00E+02 6.48E-07 + 6.00E+02 6.7E-07 8.00E+02 7.02E-07 + 1.00E+03 7.24E-07 1.50E+03 7.63E-07 + 2.00E+03 7.92E-07 3.00E+03 8.35E-07 + 4.00E+03 8.75E-07 5.00E+03 9.04E-07 + 6.00E+03 9.29E-07 8.00E+03 9.65E-07 + 1.00E+04 9.94E-07 +endmd diff --git a/tallies/bunker_mesh_p.i b/tallies/bunker_mesh_p.i index 3b8510a..517f570 100644 --- a/tallies/bunker_mesh_p.i +++ b/tallies/bunker_mesh_p.i @@ -2,8 +2,8 @@ c large mesh tally covering the beamline from monolith to S1/2 separator tmesh rmesh1:p dose 1 trans 1 cora1 500 10i 1000 20i 1500 16i 2300 - corb1 -305 5i -55 4i 45 5i 295 - corc1 -200 1i -100 5i 100 2i 300 + corb1 -305 5i -55 4i 45 9i 495 + corc1 -200 1i -100 5i 100 5i 600 c Official ESS dose conversion factors from gamma/cm²/s to Sv/h c Is described in ESS-0019931 mshmf1 0.000000 1.00000e-04 diff --git a/tallies/full_mesh.i b/tallies/full_mesh.i new file mode 100644 index 0000000..2ef99fd --- /dev/null +++ b/tallies/full_mesh.i @@ -0,0 +1,44 @@ +c large and detailed mesh tally covering the beamline from monolith to full cave +tmesh + rmesh1:n dose 1 trans 1 + cora1 500 10i 1000 20i 1500 60i 4500 + corb1 -705 13i -55 20i 55 13i 705 + corc1 -200 2i -50 30i 100 2i 300 +c Official ESS dose conversion factors from n/cm²/s to Sv/h +c Is described in ESS-0019931 + mshmf1 0.000000000 1.00000e-08 + 1.00000e-09 1.11000e-08 1.00000e-08 1.28000e-08 + 2.50000e-08 1.44000e-08 1.00000e-07 1.87000e-08 + 2.00000e-07 2.11000e-08 5.00000e-07 2.37000e-08 + 1.00000e-06 2.53000e-08 2.00000e-06 2.66000e-08 + 5.00000e-06 2.78000e-08 1.00000e-05 2.82000e-08 + 2.00000e-05 2.82000e-08 5.00000e-05 2.82000e-08 + 1.00000e-04 2.80000e-08 2.00000e-04 2.78000e-08 + 5.00000e-04 2.71000e-08 1.00000e-03 2.71000e-08 + 2.00000e-03 2.74000e-08 5.00000e-03 2.87000e-08 + 1.00000e-02 3.28000e-08 2.00000e-02 4.39000e-08 + 3.00000e-02 5.65000e-08 5.00000e-02 8.28000e-08 + 7.00000e-02 1.10000e-07 1.00000e-01 1.51000e-07 + 1.50000e-01 2.18000e-07 2.00000e-01 2.84000e-07 + 3.00000e-01 4.10000e-07 5.00000e-01 6.37000e-07 + 7.00000e-01 8.35000e-07 9.00000e-01 1.00000e-06 + 1.00000e+00 1.08000e-06 1.20000e+00 1.19000e-06 + 1.50000e+00 1.31000e-06 2.00000e+00 1.47000e-06 + 3.00000e+00 1.65000e-06 4.00000e+00 1.74000e-06 + 5.00000e+00 1.78000e-06 6.00000e+00 1.79000e-06 + 7.00000e+00 1.80000e-06 8.00000e+00 1.80000e-06 + 9.00000e+00 1.80000e-06 1.00000e+01 1.80000e-06 + 1.20000e+01 1.80000e-06 1.40000e+01 1.78000e-06 + 1.50000e+01 1.77000e-06 1.60000e+01 1.76000e-06 + 1.80000e+01 1.74000e-06 2.00000e+01 1.72000e-06 + 2.10000e+01 1.71000e-06 3.00000e+01 1.63000e-06 + 5.00000e+01 1.56000e-06 7.50000e+01 1.58000e-06 + 1.00000e+02 1.60000e-06 1.30000e+02 1.61000e-06 + 1.50000e+02 1.61000e-06 1.80000e+02 1.61000e-06 + 2.00000e+02 1.61000e-06 3.00000e+02 1.67000e-06 + 4.00000e+02 1.79000e-06 5.00000e+02 1.92000e-06 + 6.00000e+02 2.05000e-06 7.00000e+02 2.16000e-06 + 8.00000e+02 2.24000e-06 9.00000e+02 2.30000e-06 + 1.00000e+03 2.35000e-06 2.00000e+03 2.76000e-06 + 5.00000e+03 3.64000e-06 1.00000e+04 4.75000e-06 +endmd diff --git a/tallies/full_mesh_np.i b/tallies/full_mesh_np.i new file mode 100644 index 0000000..e09639f --- /dev/null +++ b/tallies/full_mesh_np.i @@ -0,0 +1,80 @@ +c large and detailed mesh tally covering the beamline from monolith to full cave +tmesh + rmesh1:n dose 1 trans 1 + cora1 500 10i 1000 20i 1500 60i 4500 + corb1 -705 13i -55 20i 55 13i 705 + corc1 -200 2i -50 30i 100 2i 300 +c Official ESS dose conversion factors from n/cm²/s to Sv/h +c Is described in ESS-0019931 + mshmf1 0.000000000 1.00000e-08 + 1.00000e-09 1.11000e-08 1.00000e-08 1.28000e-08 + 2.50000e-08 1.44000e-08 1.00000e-07 1.87000e-08 + 2.00000e-07 2.11000e-08 5.00000e-07 2.37000e-08 + 1.00000e-06 2.53000e-08 2.00000e-06 2.66000e-08 + 5.00000e-06 2.78000e-08 1.00000e-05 2.82000e-08 + 2.00000e-05 2.82000e-08 5.00000e-05 2.82000e-08 + 1.00000e-04 2.80000e-08 2.00000e-04 2.78000e-08 + 5.00000e-04 2.71000e-08 1.00000e-03 2.71000e-08 + 2.00000e-03 2.74000e-08 5.00000e-03 2.87000e-08 + 1.00000e-02 3.28000e-08 2.00000e-02 4.39000e-08 + 3.00000e-02 5.65000e-08 5.00000e-02 8.28000e-08 + 7.00000e-02 1.10000e-07 1.00000e-01 1.51000e-07 + 1.50000e-01 2.18000e-07 2.00000e-01 2.84000e-07 + 3.00000e-01 4.10000e-07 5.00000e-01 6.37000e-07 + 7.00000e-01 8.35000e-07 9.00000e-01 1.00000e-06 + 1.00000e+00 1.08000e-06 1.20000e+00 1.19000e-06 + 1.50000e+00 1.31000e-06 2.00000e+00 1.47000e-06 + 3.00000e+00 1.65000e-06 4.00000e+00 1.74000e-06 + 5.00000e+00 1.78000e-06 6.00000e+00 1.79000e-06 + 7.00000e+00 1.80000e-06 8.00000e+00 1.80000e-06 + 9.00000e+00 1.80000e-06 1.00000e+01 1.80000e-06 + 1.20000e+01 1.80000e-06 1.40000e+01 1.78000e-06 + 1.50000e+01 1.77000e-06 1.60000e+01 1.76000e-06 + 1.80000e+01 1.74000e-06 2.00000e+01 1.72000e-06 + 2.10000e+01 1.71000e-06 3.00000e+01 1.63000e-06 + 5.00000e+01 1.56000e-06 7.50000e+01 1.58000e-06 + 1.00000e+02 1.60000e-06 1.30000e+02 1.61000e-06 + 1.50000e+02 1.61000e-06 1.80000e+02 1.61000e-06 + 2.00000e+02 1.61000e-06 3.00000e+02 1.67000e-06 + 4.00000e+02 1.79000e-06 5.00000e+02 1.92000e-06 + 6.00000e+02 2.05000e-06 7.00000e+02 2.16000e-06 + 8.00000e+02 2.24000e-06 9.00000e+02 2.30000e-06 + 1.00000e+03 2.35000e-06 2.00000e+03 2.76000e-06 + 5.00000e+03 3.64000e-06 1.00000e+04 4.75000e-06 +c + rmesh11:p dose 2 trans 1 + cora11 500 10i 1000 20i 1500 60i 4500 + corb11 -705 13i -55 20i 55 13i 705 + corc11 -200 2i -50 30i 100 2i 300 +c Official ESS dose conversion factors from gamma/cm²/s to Sv/h +c Is described in ESS-0019931 + mshmf2 0.000000 1.00000e-04 + 1.00E-02 2.47E-10 1.50E-02 5.62E-10 + 2.00E-02 8.1E-10 3.00E-02 1.13E-09 + 4.00E-02 1.26E-09 5.00E-02 1.33E-09 + 6.00E-02 1.4E-09 7.00E-02 1.49E-09 + 8.00E-02 1.6E-09 1.00E-01 1.87E-09 + 1.50E-01 2.69E-09 2.00E-01 3.6E-09 + 3.00E-01 5.44E-09 4.00E-01 7.2E-09 + 5.00E-01 8.89E-09 5.11E-01 9.07E-09 + 6.00E-01 1.05E-08 6.62E-01 1.14E-08 + 8.00E-01 1.34E-08 1.00E+00 1.62E-08 + 1.12E+00 1.76E-08 1.33E+00 2.01E-08 + 1.50E+00 2.2E-08 2.00E+00 2.69E-08 + 3.00E+00 3.51E-08 4.00E+00 4.21E-08 + 5.00E+00 4.82E-08 6.00E+00 5.4E-08 + 6.13E+00 5.47E-08 8.00E+00 6.7E-08 + 1.00E+01 7.92E-08 1.50E+01 1.09E-07 + 2.00E+01 1.38E-07 3.00E+01 1.85E-07 + 4.00E+01 2.23E-07 5.00E+01 2.6E-07 + 6.00E+01 2.95E-07 8.00E+01 3.52E-07 + 1.00E+02 3.96E-07 1.50E+02 4.68E-07 + 2.00E+02 5.15E-07 3.00E+02 5.8E-07 + 4.00E+02 6.19E-07 5.00E+02 6.48E-07 + 6.00E+02 6.7E-07 8.00E+02 7.02E-07 + 1.00E+03 7.24E-07 1.50E+03 7.63E-07 + 2.00E+03 7.92E-07 3.00E+03 8.35E-07 + 4.00E+03 8.75E-07 5.00E+03 9.04E-07 + 6.00E+03 9.29E-07 8.00E+03 9.65E-07 + 1.00E+04 9.94E-07 +endmd diff --git a/tallies/full_mesh_p.i b/tallies/full_mesh_p.i new file mode 100644 index 0000000..2ceb94f --- /dev/null +++ b/tallies/full_mesh_p.i @@ -0,0 +1,38 @@ +c large and detailed mesh tally covering the beamline from monolith to full cave +tmesh + rmesh1:p dose 1 trans 1 + cora1 500 10i 1000 20i 1500 60i 4500 + corb1 -705 13i -55 20i 55 13i 705 + corc1 -200 2i -50 30i 100 2i 300 +c Official ESS dose conversion factors from gamma/cm²/s to Sv/h +c Is described in ESS-0019931 + mshmf1 0.000000 1.00000e-04 + 1.00E-02 2.47E-10 1.50E-02 5.62E-10 + 2.00E-02 8.1E-10 3.00E-02 1.13E-09 + 4.00E-02 1.26E-09 5.00E-02 1.33E-09 + 6.00E-02 1.4E-09 7.00E-02 1.49E-09 + 8.00E-02 1.6E-09 1.00E-01 1.87E-09 + 1.50E-01 2.69E-09 2.00E-01 3.6E-09 + 3.00E-01 5.44E-09 4.00E-01 7.2E-09 + 5.00E-01 8.89E-09 5.11E-01 9.07E-09 + 6.00E-01 1.05E-08 6.62E-01 1.14E-08 + 8.00E-01 1.34E-08 1.00E+00 1.62E-08 + 1.12E+00 1.76E-08 1.33E+00 2.01E-08 + 1.50E+00 2.2E-08 2.00E+00 2.69E-08 + 3.00E+00 3.51E-08 4.00E+00 4.21E-08 + 5.00E+00 4.82E-08 6.00E+00 5.4E-08 + 6.13E+00 5.47E-08 8.00E+00 6.7E-08 + 1.00E+01 7.92E-08 1.50E+01 1.09E-07 + 2.00E+01 1.38E-07 3.00E+01 1.85E-07 + 4.00E+01 2.23E-07 5.00E+01 2.6E-07 + 6.00E+01 2.95E-07 8.00E+01 3.52E-07 + 1.00E+02 3.96E-07 1.50E+02 4.68E-07 + 2.00E+02 5.15E-07 3.00E+02 5.8E-07 + 4.00E+02 6.19E-07 5.00E+02 6.48E-07 + 6.00E+02 6.7E-07 8.00E+02 7.02E-07 + 1.00E+03 7.24E-07 1.50E+03 7.63E-07 + 2.00E+03 7.92E-07 3.00E+03 8.35E-07 + 4.00E+03 8.75E-07 5.00E+03 9.04E-07 + 6.00E+03 9.29E-07 8.00E+03 9.65E-07 + 1.00E+04 9.94E-07 +endmd