Test run with SSA source using mesh tally

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

rssa

@@ -0,0 +1 @@
+results/e2-estia_baseline_v09-0-gfab54af.ssa

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 

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