diff --git a/e2-estia.i b/e2-estia.i index a0a2ea5..9e2b30e 100644 --- a/e2-estia.i +++ b/e2-estia.i @@ -149,7 +149,8 @@ c read file=tallies/bunker_mesh_np.i echo c read file=tallies/bunker_mesh_act.i echo c read file=tallies/full_mesh.i echo c read file=tallies/full_mesh_p.i echo -read file=tallies/full_mesh_np.i echo +c read file=tallies/full_mesh_np.i echo +read file=tallies/vs_fine_mesh.i echo c read file=tallies/boundary_meshes.i echo c read file=tallies/boundary_wall_meshes.i echo c read file=tallies/selene_source_mesh_p.i echo @@ -204,8 +205,8 @@ c --------------- PHYSICS CARDS -------------------------- c ------------------------------------------------------- c dbcn 375642321 j j j j j j j j j j j 15291711 c stop calculation either at XeY histories, 6.5 days or tally 15 error <0.001 -c stop nps 392426460 CTME 4e5 $f31 0.01 -stop nps 1e8 CTME 4e5 $f31 0.01 +stop nps 39242646 CTME 4e5 $f31 0.01 +c stop nps 1e8 CTME 4e5 $f31 0.01 rand gen=2 seed=29976355590901 $ default 19073486328125 c histp mode n p h / z $ neutron, photon, proton, pi+, pi0 @@ -237,18 +238,18 @@ c * Source for simulation from proton beam (needs different WW and ESplit) 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 -read file=sources/e2_2m_mod.i echo +c read file=sources/e2_2m_mod.i echo c * Estimated source within bunker wall to optimize collimator efficiency c read file=sources/bunker_wall.i echo $ approximate, not for official results c * McStas source terms for gamma doeses -c read file=sources/mcstas_rssa.i echo +read file=sources/mcstas_rssa.i echo c read file=sources/mcstas_source_vs.i echo c read file=sources/mcstas_source_sample.i echo c read file=sources/mcstas_source_mf.i echo c read file=sources/mcstas_gammas.i echo c read file=sources/mcstas_gammas_2.i echo c read file=sources/mcstas_source_shutter.i echo -prdmp 5e6 -30 1 2 5e5 +prdmp 5e5 -30 1 2 1e5 print 30 32 40 c 120 $O importance function quality c 128 $O universe map diff --git a/geometry/estia_surfaces.i b/geometry/estia_surfaces.i index 33f1c86..65166d9 100644 --- a/geometry/estia_surfaces.i +++ b/geometry/estia_surfaces.i @@ -771,8 +771,8 @@ c VS Blades (open) 224502 p -0.2588190294771241 0.9659264084517616 0.0 -290.42 224503 p -0.2588190294771241 0.9659264084517616 0.0 -290.52 224504 p 0.2588190294771241 -0.9659264084517616 0.0 292.92 -221501 pz 5.061 -221502 pz 1.061 +221501 pz 5.011 +221502 pz 1.011 221503 pz -4.739 221504 pz -0.739 224505 p 0.9659264084517616 0.2588190294771241 0.0 1073.54 @@ -782,7 +782,7 @@ c VS Blades (open) 224513 p -0.2588190294771241 0.9659264084517616 0.0 -292.856 224514 p -0.2588190294771241 0.9659264084517616 0.0 -290.906 221511 pz 5.261 -221512 pz 0.539 +221512 pz -0.489 221513 pz 1.761 221514 pz -4.539 224515 p -0.9659264084517616 -0.2588190294771241 0.0 -1078.921 diff --git a/tallies/vs_fine_mesh.i b/tallies/vs_fine_mesh.i new file mode 100644 index 0000000..1f1bf5f --- /dev/null +++ b/tallies/vs_fine_mesh.i @@ -0,0 +1,53 @@ +c large and detailed mesh tally covering the beamline from monolith to full cave +c plus one very detailed tally to get the source details around VS +tmesh + rmesh1:n dose 1 trans 101 + cora1 500 9i 1000 22i 1460 + corb1 -700 11i -100 9i 100 3i 300 + corc1 -200 1i -100 9i 100 5i 400 + rmesh11:n dose 1 trans 101 + cora11 1460 151i 4500 + corb11 -400 59i 800 + corc11 -200 39i 600 + rmesh21:n dose 1 trans 101 + cora21 1000 399i 1400 + corb21 -50 99i 50 + corc21 -50 99i 50 +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