From 4cefa58ab3ddd851e0325ec1770dce6f0102cd83 Mon Sep 17 00:00:00 2001 From: Artur Glavic Date: Fri, 31 May 2019 20:20:38 +0200 Subject: [PATCH] Generate weight window mesh for NBPI window --- e2-estia.i | 21 +++++++++++---------- sources/proton.i | 7 ++++--- 2 files changed, 15 insertions(+), 13 deletions(-) diff --git a/e2-estia.i b/e2-estia.i index 85c6900..c46da38 100644 --- a/e2-estia.i +++ b/e2-estia.i @@ -93,16 +93,17 @@ c ------------------------------------------------------- c --------------- WEIGHT CARDS -------------------------- c ------------------------------------------------------- c weight window generator, can be run with previous weight window files -c WWGE:n 1e-3 1e-1 1 10 100 3e3 +WWGE:n 1e-3 1e-1 1 10 100 3e3 c mesh based weight windows, can be reused after change of cell structure -c WWG 55 0 0 +WWG 514 0 0 c MESH GEOM=CYL ORIGIN=0 0 -500 REF=133.7 -176.7 13.5 AXS=0 0 1 -c IMESH=200 550 1150 1500 2112.01 2157.17 4600 -c IINTS= 2 5 1 10 1 1 4 -c JMESH=100 300 460 500 540 700 900 1200 -c JINTS= 2 4 4 8 8 4 4 6 -c KMESH=298 303 309 315 360 -c KINTS= 14 1 20 1 2 +MESH GEOM=CYL ORIGIN=0 0 -500 REF=0 0 13.4 AXS=0 0 1 + IMESH=200 550 1150 1500 2112.01 2157.17 4600 + IINTS= 2 5 1 10 1 1 4 + JMESH=100 300 460 500 540 700 900 1200 + JINTS= 2 4 4 8 8 4 4 6 + KMESH=298 303 309 315 360 + KINTS= 14 1 20 1 2 c generate cell based weight windows, needs to be updated on geometry change c WWG 55 1042 0 c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++ @@ -160,7 +161,7 @@ c ------------------------------------------------------- c --------------- PHYSICS CARDS -------------------------- c ------------------------------------------------------- c dbcn 375642321 j j j j j j j j j j j 15291711 -nps 1e5 +nps 1e7 c histp mode n p h / d $ deuteron @@ -204,7 +205,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 5e8 -60 1 2 5e7 +prdmp 1e6 -60 1 2 1e5 print 30 32 40 $ 120 128 150 160 161 162 c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++ diff --git a/sources/proton.i b/sources/proton.i index 650f280..219e758 100644 --- a/sources/proton.i +++ b/sources/proton.i @@ -1,11 +1,12 @@ c Original definition of source from proton beam -cut:n j 0.99e-06 +c cut:n j 0.99e-06 +cut:n j 0 cut:p 1e+08 0.5 c weight windows from external file generated by WWG c make sure to supply the correct WWINP file c 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 -c esplt:n 0.05 1e-6 0.2 1e-3 1 10 3 100 10 3000 +c energy splitting for higher neutron energies and roulette below 1keV +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=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