diff --git a/e2-estia.i b/e2-estia.i index a26b63f..4bc727e 100644 --- a/e2-estia.i +++ b/e2-estia.i @@ -94,14 +94,16 @@ c ----------------------------------------------------------- c ------------------- TALLY CARDS --------------------------- c c read file=tallies/bunker_mesh.i echo -read file=tallies/inbunker_f4.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 c read file=tallies/e1_e2_f4.i echo -read file=tallies/e1_e2_f5.i echo +c read file=tallies/e1_e2_f5.i echo c read file=tallies/shutter_f4.i echo c read file=tallies/shutter_f5.i echo -read file=tallies/optical_cave_f5.i echo +c read file=tallies/optical_cave_f5.i echo +c read file=tallies/experimental_cave_f5.i echo +read file=tallies/cave_mesh.i echo c c ========================================== c @@ -129,14 +131,14 @@ c ------------------------------------------------------- c --------------- PHYSICS CARDS -------------------------- c ------------------------------------------------------- c dbcn 375642321 j j j j j j j j j j j 15291711 -nps 1e9 +nps 3e7 c histp 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 43r -imp:h,/,d,t,s,a 0 0 1 688r -imp:n,p 0 0 1 688r +vol 1 507r 10j 1 107r 3848.5 10104 1 13r 479.5 479.5 913 913 1734 1 46r +imp:h,/,d,t,s,a 0 0 1 691r +imp:n,p 0 0 1 691r cut:h,/,d,t,s,a 1e+08 1e-3 c cut:n j 0.99E-6 cut:n j 1e-11 @@ -152,9 +154,9 @@ c read file=sources/proton.i echo c * Source Term from Valentina for 2m adapted to size c read file=sources/e2_2m.i echo c * McStas source terms for gamma doeses -read file=sources/mcstas_source_vs.i echo -c read file=sources/mcstas_source_sample.i echo -prdmp 5e7 -60 1 2 5e6 +c read file=sources/mcstas_source_vs.i echo +read file=sources/mcstas_source_sample.i echo +prdmp 5e6 -60 1 2 1e6 print 30 32 $ 120 128 150 160 161 162 c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++ diff --git a/geometry/ess_cells.i b/geometry/ess_cells.i index 2bbe7b5..fbcdf21 100644 --- a/geometry/ess_cells.i +++ b/geometry/ess_cells.i @@ -13,6 +13,7 @@ c Main ESS plant geometry without the instrument :(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) 3 118 0.122325 24 4 -19 4 0 159 -137 -21 153 -154 5 118 0.122325 37 -5 -19 diff --git a/geometry/estia_cave.i b/geometry/estia_cave.i index 010eac3..dafdb56 100644 --- a/geometry/estia_cave.i +++ b/geometry/estia_cave.i @@ -16,3 +16,7 @@ c then cut corner and add extended roof #6111 #8501 #8511 #9001 c sample position 9001 19 -8.65 -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 diff --git a/geometry/estia_guides.i b/geometry/estia_guides.i index d0b87a6..5ebf61e 100644 --- a/geometry/estia_guides.i +++ b/geometry/estia_guides.i @@ -105,8 +105,8 @@ c Selene 2 glass (change material!!) 8511 12 -2.51 8121 -8122 -8913 -8914 8915 #(-8911 -8912) #(-8515 -8517) 3 c Air outside the guide shielding -8601 4 -1.2E-3 72 -8709 -104 -8201 8203 8204 +8601 0 72 -8709 -104 -8201 8203 8204 #708 -8602 4 -1.2E-3 ((112 -8202) +8602 0 ((112 -8202) :(8202 -8302 8101)) 72 -8709 8203 8204 #709 diff --git a/geometry/estia_surfaces.i b/geometry/estia_surfaces.i index 998ef93..0b4f04f 100644 --- a/geometry/estia_surfaces.i +++ b/geometry/estia_surfaces.i @@ -614,6 +614,7 @@ c Outside 8710 p -0.582123 0.813101 0.000000 -4364.795 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 8801 pz 57.941 8802 pz 108.941 diff --git a/tallies/bunker_mesh.i b/tallies/bunker_mesh.i index 55f0707..72e54df 100644 --- a/tallies/bunker_mesh.i +++ b/tallies/bunker_mesh.i @@ -1,8 +1,44 @@ -c large mesh tally covering the beamline from monolith +c large mesh tally covering the beamline from monolith to S1/2 separator tmesh - rmesh1:n dose 31 1 2 1 trans 1 + 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 +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 -fm1 8.7374e15 diff --git a/tallies/cave_mesh.i b/tallies/cave_mesh.i new file mode 100644 index 0000000..8775ca5 --- /dev/null +++ b/tallies/cave_mesh.i @@ -0,0 +1,80 @@ +c large mesh tally covering the experimental cave area +tmesh + rmesh1:n dose 1 trans 1 + cora1 2800 170i 4500 + corb1 -370 109i 720 + corc1 -200 45i 250 +c Official ESS dose conversion factors from n/cm²/s to Sv/h +c Is described in ESS-0019931 + mshmf1 0.000000000 1.11000e-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 2800 170i 4500 + corb11 -370 109i 720 + corc11 -200 45i 250 +c Official ESS dose conversion factors from gamma/cm²/s to Sv/h +c Is described in ESS-0019931 + mshmf2 0.000000000 2.47000e-10 + 1.00000e-02 2.47000e-10 1.50000e-02 5.62000e-10 + 2.00000e-02 8.10000e-10 3.00000e-02 1.13000e-09 + 4.00000e-02 1.26000e-09 5.00000e-02 1.33000e-09 + 6.00000e-02 1.40000e-09 7.00000e-02 1.49000e-09 + 8.00000e-02 1.60000e-09 1.00000e-01 1.87000e-09 + 1.50000e-01 2.69000e-09 2.00000e-01 3.60000e-09 + 3.00000e-01 5.44000e-09 4.00000e-01 7.20000e-09 + 5.00000e-01 8.89000e-09 5.11000e-01 9.07000e-09 + 6.00000e-01 1.05000e-08 6.62000e-01 1.14000e-08 + 8.00000e-01 1.34000e-08 1.00000e+00 1.62000e-08 + 1.12000e+00 1.76000e-08 1.33000e+00 2.01000e-08 + 1.50000e+00 2.20000e-08 2.00000e+00 2.69000e-08 + 3.00000e+00 3.51000e-08 4.00000e+00 4.21000e-08 + 5.00000e+00 4.82000e-08 6.00000e+00 5.40000e-08 + 6.13000e+00 5.47000e-08 8.00000e+00 6.70000e-08 + 1.00000e+01 7.92000e-08 1.50000e+01 1.09000e-07 + 2.00000e+01 1.38000e-07 3.00000e+01 1.85000e-07 + 4.00000e+01 2.23000e-07 5.00000e+01 2.60000e-07 + 6.00000e+01 2.95000e-07 8.00000e+01 3.52000e-07 + 1.00000e+02 3.96000e-07 1.50000e+02 4.68000e-07 + 2.00000e+02 5.15000e-07 3.00000e+02 5.80000e-07 + 4.00000e+02 6.19000e-07 5.00000e+02 6.48000e-07 + 6.00000e+02 6.70000e-07 8.00000e+02 7.02000e-07 + 1.00000e+03 7.24000e-07 1.50000e+03 7.63000e-07 + 2.00000e+03 7.92000e-07 3.00000e+03 8.35000e-07 + 4.00000e+03 8.75000e-07 5.00000e+03 9.04000e-07 + 6.00000e+03 9.29000e-07 8.00000e+03 9.65000e-07 + 1.00000e+04 9.94000e-07 +endmd diff --git a/tallies/experimental_cave_f5.i b/tallies/experimental_cave_f5.i new file mode 100644 index 0000000..bc200b5 --- /dev/null +++ b/tallies/experimental_cave_f5.i @@ -0,0 +1,39 @@ +c f5 detector tallies for dose outside of experimental cave +c +fc205 neutron flux outside shielding closest to sample +f205:n 2207 -2670 -29 0.0 +e205 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10 + 1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9 + 1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8 + 1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7 + 1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6 + 1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5 + 1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4 + 1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3 + 1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2 + 1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1 + 1e-0 1.58 2.51 3.98 6.31 + 1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1 + 1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2 + 1e+3 1.58e+3 2.51e+3 +fq205 e f +fc215 photon flux outside shielding closest to sample +f215:p 2207 -2670 -29 0.0 +e215 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10 + 1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9 + 1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8 + 1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7 + 1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6 + 1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5 + 1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4 + 1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3 + 1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2 + 1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1 + 1e-0 1.58 2.51 3.98 6.31 + 1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1 + 1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2 + 1e+3 1.58e+3 2.51e+3 +fq215 e f +c 1-255 256 257 268 259-265 266-268 269-508 509-end +c pd205 0.001 254r 0.1 0.001 1 0.001 6r 1 2r 0.001 239r 0.1 110r 1 71r +c pd215 0.001 254r 0.1 0.001 1 0.001 6r 1 2r 0.001 239r 0.1 110r 1 71r