diff --git a/e2-estia.i b/e2-estia.i index a32290e..6abb7d6 100644 --- a/e2-estia.i +++ b/e2-estia.i @@ -139,7 +139,7 @@ c KMESH=298 303 309 315 360 c KINTS= 14 1 20 1 2 c generate cell based weight windows, needs to be updated on geometry change c WWG 34 61003 0 -TF34 2 $ set second tally cell (behind BWC) as reference +c TF34 2 $ set second tally cell (behind BWC) as reference c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++ c ----------------------------------------------------------- c ------------------- TALLY CARDS --------------------------- @@ -148,7 +148,7 @@ 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/bunker_mesh_act.i echo -c read file=tallies/full_mesh.i echo +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/vs_fine_mesh.i echo @@ -158,12 +158,12 @@ c read file=tallies/selene_source_mesh_p.i echo c read file=tallies/NBPI_f4.i echo c read file=tallies/BBG_f4.i echo c read file=tallies/BBG_f5.i echo -read file=tallies/inbunker_f4.i echo +c read file=tallies/inbunker_f4.i echo c read file=tallies/CPC1_f5.i echo c read file=tallies/CPC2_f5.i echo -read file=tallies/chopper_pit_doses.i echo +c read file=tallies/chopper_pit_doses.i echo c read file=tallies/bunker_feedthrough_f1.i echo -read file=tallies/bunker_wall_f4.i echo +c read file=tallies/bunker_wall_f4.i echo c read file=tallies/bunker_wall_f5.i echo c read file=tallies/pre_selene1_mesh.i read file=tallies/e1_e2_f4.i echo @@ -209,9 +209,10 @@ 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 39242646 CTME 4e5 $f31 0.01 -stop nps 1e9 CTME 4e5 $f31 0.01 +c stop nps 1e9 CTME 4e5 $f31 0.01 +nps 14e10 $ works with rssa file, value is scaled to source particles c write particle tracks -ssw 206505 pty=n $ record all neutron tracks crossing outward +c ssw 206505 pty=n $ record all neutron tracks crossing outward c random number generator options rand gen=2 seed=29976355590901 $ default 19073486328125 c histp @@ -251,8 +252,8 @@ 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/rssa.i echo +c read file=sources/e2_2m_mod.i echo +read file=sources/rssa.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 @@ -263,7 +264,7 @@ 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 5e7 -30 1 2 1e6 +prdmp -120 -30 1 2 -1 print 30 32 40 c 120 $O importance function quality c 128 $O universe map diff --git a/results/e2-estia_baseline_v09-0-gfab54af.out.gz b/results/e2-estia_baseline_v09-0-gfab54af.out.gz new file mode 100644 index 0000000..d4c5e3c Binary files /dev/null and b/results/e2-estia_baseline_v09-0-gfab54af.out.gz differ diff --git a/results/e2-estia_baseline_v09-0-gfab54af.ssa b/results/e2-estia_baseline_v09-0-gfab54af.ssa new file mode 100644 index 0000000..9c6bd0f Binary files /dev/null and b/results/e2-estia_baseline_v09-0-gfab54af.ssa differ diff --git a/rssa b/rssa new file mode 120000 index 0000000..835dd5c --- /dev/null +++ b/rssa @@ -0,0 +1 @@ +results/e2-estia_baseline_v09-0-gfab54af.ssa \ No newline at end of file diff --git a/sources/rssa.i b/sources/rssa.i index 7b7b6c6..43612c3 100644 --- a/sources/rssa.i +++ b/sources/rssa.i @@ -1,8 +1,7 @@ c source definition based on surface source generated from 2m source term -cut:n j 0 +cut:n j j -0.2 -0.1 $ put weight cut-off for 0-WW to 10% of source weight c cut all photons that can't produce neutrons cut:p 1e+08 1.667 -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 4 j -1 j 1.0 1j 2.0E+13 1j +WWP:n j j 4 j -1 j 1.0 1j 6.0E+14 1j ssr diff --git a/tallies/full_mesh.i b/tallies/full_mesh.i index b3de97c..eefeacc 100644 --- a/tallies/full_mesh.i +++ b/tallies/full_mesh.i @@ -1,11 +1,7 @@ c large and detailed mesh tally covering the beamline from monolith to full cave tmesh - rmesh1:n dose 1 - 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 - cora11 1460 151i 4500 + cora11 1300 159i 4500 corb11 -400 59i 800 corc11 -200 39i 600 c Official ESS dose conversion factors from n/cm²/s to Sv/h