From 4486b279d6ad2fc785e796f0cabda12e17e12bff Mon Sep 17 00:00:00 2001 From: Artur Glavic Date: Sat, 4 May 2019 17:18:03 +0200 Subject: [PATCH] Run closed shutter mesh with MF derived source --- e2-estia.i | 9 ++- geometry/estia_cave.i | 1 + geometry/estia_surfaces.i | 2 + sources/mcstas_source_shutter.i | 133 ++++++++++++++++++++++++++++++++ 4 files changed, 141 insertions(+), 4 deletions(-) create mode 100644 sources/mcstas_source_shutter.i diff --git a/e2-estia.i b/e2-estia.i index 702eb28..8a1d1fe 100644 --- a/e2-estia.i +++ b/e2-estia.i @@ -107,7 +107,7 @@ 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 -c read file=tallies/full_mesh_p.i echo +read file=tallies/full_mesh_p.i echo c read file=tallies/full_mesh_np.i echo c read file=tallies/inbunker_f4.i echo c read file=tallies/chopper_pit_doses.i echo @@ -115,7 +115,7 @@ c read file=tallies/bunker_wall_f4.i echo c read file=tallies/e1_e2_f4.i echo c read file=tallies/e1_e2_f5.i echo c read file=tallies/shutter_f4.i echo -read file=tallies/shutter_f5.i echo +c read file=tallies/shutter_f5.i echo c read file=tallies/shutter_sample_f5.i echo c read file=tallies/selene1_walls_f4.i echo c read file=tallies/optical_cave_f5.i echo @@ -147,7 +147,7 @@ c ------------------------------------------------------- c --------------- PHYSICS CARDS -------------------------- c ------------------------------------------------------- c dbcn 375642321 j j j j j j j j j j j 15291711 -nps 1e10 +nps 2e8 c histp mode n p n h / d t s a c volume calculated for tally cells (spheres) @@ -168,7 +168,7 @@ c read file=sources/geometry_test.i echo 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 -read file=sources/e2_2m.i echo +c read file=sources/e2_2m.i echo c read file=sources/e2_2m_act.i echo c * McStas source terms for gamma doeses c read file=sources/mcstas_source_vs.i echo @@ -176,6 +176,7 @@ 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 +read file=sources/mcstas_source_shutter.i echo prdmp 1e8 -60 1 2 1e7 print 30 32 40 $ 120 128 150 160 161 162 c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++ diff --git a/geometry/estia_cave.i b/geometry/estia_cave.i index 2dea7a1..fe34128 100644 --- a/geometry/estia_cave.i +++ b/geometry/estia_cave.i @@ -39,3 +39,4 @@ c area around cave c visualize Selene gamma sources (comment for real runs) c 50001 0 -50001 c 50002 0 -50002 +c 50003 0 -50003 diff --git a/geometry/estia_surfaces.i b/geometry/estia_surfaces.i index 5f1f2b1..6d37846 100644 --- a/geometry/estia_surfaces.i +++ b/geometry/estia_surfaces.i @@ -957,4 +957,6 @@ c c visualize Selene gamma sources c 50001 1 rpp 1339.25 2060.75 2 8 2 8 c 50002 1 rpp 2539.25 3260.75 -8 -2 -8 -2 +c visualize Shutter Source +c 50003 3 rpp -170.5 -169.5 0 4 -6 6 c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++ diff --git a/sources/mcstas_source_shutter.i b/sources/mcstas_source_shutter.i new file mode 100644 index 0000000..153fd25 --- /dev/null +++ b/sources/mcstas_source_shutter.i @@ -0,0 +1,133 @@ +c cold neutron source from McStas simulation for middle focus position +c activation (delayed particles) +c act nonfiss=p,a,e dg=lines +cut:n j 1e-11 0 0 +cut:p j 0.01 +c as this is primary for gamma shielding purposes, +c no directional dependence has been considered +sdef dir=1 erg=d11 x=-170. y=d23 z=d24 vec=1 0 0 tr=3 + par=n ara=72.0 wgt=3.33085e+09 +c scaling wgt is to flux from McStas nps +si1 -1 0.99976203 1 +sp1 0 0 1 +si11 0.00000e+00 + 9.009970e-11 9.070338e-11 9.131316e-11 9.192910e-11 9.255130e-11 + 9.317984e-11 9.381480e-11 9.445628e-11 9.510435e-11 9.575912e-11 + 9.642068e-11 9.708911e-11 9.776452e-11 9.844700e-11 9.913665e-11 + 9.983357e-11 1.005379e-10 1.012497e-10 1.019690e-10 1.026961e-10 + 1.034309e-10 1.041737e-10 1.049245e-10 1.056835e-10 1.064507e-10 + 1.072263e-10 1.080104e-10 1.088032e-10 1.096047e-10 1.104151e-10 + 1.112345e-10 1.120630e-10 1.129009e-10 1.137482e-10 1.146051e-10 + 1.154716e-10 1.163481e-10 1.172346e-10 1.181312e-10 1.190382e-10 + 1.199556e-10 1.208837e-10 1.218226e-10 1.227725e-10 1.237336e-10 + 1.247059e-10 1.256898e-10 1.266854e-10 1.276929e-10 1.287124e-10 + 1.297441e-10 1.307884e-10 1.318453e-10 1.329150e-10 1.339978e-10 + 1.350939e-10 1.362035e-10 1.373269e-10 1.384642e-10 1.396157e-10 + 1.407816e-10 1.419621e-10 1.431576e-10 1.443683e-10 1.455943e-10 + 1.468361e-10 1.480938e-10 1.493677e-10 1.506582e-10 1.519654e-10 + 1.532897e-10 1.546315e-10 1.559909e-10 1.573683e-10 1.587640e-10 + 1.601784e-10 1.616118e-10 1.630645e-10 1.645369e-10 1.660293e-10 + 1.675421e-10 1.690757e-10 1.706304e-10 1.722067e-10 1.738049e-10 + 1.754255e-10 1.770688e-10 1.787354e-10 1.804256e-10 1.821399e-10 + 1.838787e-10 1.856425e-10 1.874319e-10 1.892472e-10 1.910890e-10 + 1.929579e-10 1.948543e-10 1.967788e-10 1.987320e-10 2.007144e-10 + 2.027266e-10 2.047692e-10 2.068428e-10 2.089481e-10 2.110857e-10 + 2.132563e-10 2.154605e-10 2.176991e-10 2.199727e-10 2.222822e-10 + 2.246282e-10 2.270115e-10 2.294330e-10 2.318934e-10 2.343937e-10 + 2.369345e-10 2.395170e-10 2.421418e-10 2.448101e-10 2.475227e-10 + 2.502806e-10 2.530849e-10 2.559366e-10 2.588368e-10 2.617865e-10 + 2.647870e-10 2.678393e-10 2.709447e-10 2.741045e-10 2.773198e-10 + 2.805921e-10 2.839226e-10 2.873127e-10 2.907640e-10 2.942778e-10 + 2.978556e-10 3.014992e-10 3.052100e-10 3.089897e-10 3.128401e-10 + 3.167629e-10 3.207600e-10 3.248331e-10 3.289844e-10 3.332158e-10 + 3.375293e-10 3.419271e-10 3.464115e-10 3.509846e-10 3.556489e-10 + 3.604068e-10 3.652608e-10 3.702136e-10 3.752677e-10 3.804261e-10 + 3.856916e-10 3.910671e-10 3.965558e-10 4.021609e-10 4.078857e-10 + 4.137336e-10 4.197082e-10 4.258131e-10 4.320522e-10 4.384294e-10 + 4.449489e-10 4.516148e-10 4.584317e-10 4.654042e-10 4.725369e-10 + 4.798348e-10 4.873031e-10 4.949472e-10 5.027725e-10 5.107849e-10 + 5.189904e-10 5.273952e-10 5.360058e-10 5.448290e-10 5.538719e-10 + 5.631419e-10 5.726465e-10 5.823937e-10 5.923920e-10 6.026500e-10 + 6.131767e-10 6.239817e-10 6.350749e-10 6.464664e-10 6.581673e-10 + 6.701887e-10 6.825426e-10 6.952412e-10 7.082975e-10 7.217250e-10 + 7.355381e-10 7.497515e-10 7.643809e-10 7.794428e-10 7.949542e-10 + 8.109333e-10 8.273991e-10 8.443715e-10 8.618716e-10 8.799214e-10 + 8.985442e-10 9.177645e-10 9.376082e-10 9.581026e-10 9.792762e-10 + 1.001160e-09 1.023785e-09 1.047186e-09 1.071398e-09 1.096460e-09 + 1.122412e-09 1.149296e-09 1.177158e-09 1.206046e-09 1.236010e-09 + 1.267104e-09 1.299387e-09 1.332920e-09 1.367768e-09 1.404000e-09 + 1.441691e-09 1.480921e-09 1.521774e-09 1.564341e-09 1.608719e-09 + 1.655013e-09 1.703335e-09 1.753804e-09 1.806549e-09 1.861711e-09 + 1.919438e-09 1.979892e-09 2.043248e-09 2.109694e-09 2.179435e-09 + 2.252693e-09 2.329707e-09 2.410739e-09 2.496073e-09 2.586020e-09 + 2.680917e-09 2.781136e-09 2.887082e-09 2.999198e-09 3.117975e-09 + 3.243949e-09 3.377716e-09 3.519930e-09 3.671319e-09 3.832688e-09 + 4.004935e-09 4.189061e-09 4.386182e-09 4.597551e-09 4.824575e-09 + 5.068840e-09 5.332136e-09 5.616492e-09 5.924217e-09 6.257943e-09 + 6.620683e-09 7.015902e-09 7.447595e-09 7.920393e-09 8.439684e-09 + 9.011772e-09 9.644063e-09 1.034531e-08 1.112592e-08 1.199833e-08 + 1.297753e-08 1.408168e-08 1.533297e-08 1.675878e-08 1.839312e-08 + 2.027874e-08 2.246991e-08 2.503641e-08 2.806911e-08 +sp11 0.00000e+00 + 1.214840e+04 1.256345e+04 1.302044e+04 1.326152e+04 1.379966e+04 + 1.421472e+04 1.475672e+04 1.528101e+04 1.583387e+04 1.642523e+04 + 1.689077e+04 1.734300e+04 1.778793e+04 1.864549e+04 1.934852e+04 + 2.014583e+04 2.083637e+04 2.110292e+04 2.219619e+04 2.249020e+04 + 2.357387e+04 2.393510e+04 2.534934e+04 2.582712e+04 2.715238e+04 + 2.786449e+04 2.887503e+04 2.949026e+04 3.083067e+04 3.151190e+04 + 3.245292e+04 3.405610e+04 3.445021e+04 3.619378e+04 3.735710e+04 + 3.828247e+04 3.996155e+04 4.109754e+04 4.278032e+04 4.438702e+04 + 4.546359e+04 4.713878e+04 4.830520e+04 4.981591e+04 5.213248e+04 + 5.434103e+04 5.570910e+04 5.739441e+04 5.832051e+04 6.047995e+04 + 6.407619e+04 6.473340e+04 6.728596e+04 6.979223e+04 7.170640e+04 + 7.417202e+04 7.763219e+04 7.943879e+04 8.224429e+04 8.616590e+04 + 8.848352e+04 9.056017e+04 9.319461e+04 9.511508e+04 9.947364e+04 + 1.036847e+05 1.084932e+05 1.116687e+05 1.146604e+05 1.193212e+05 + 1.225996e+05 1.269869e+05 1.324391e+05 1.357706e+05 1.401152e+05 + 1.464706e+05 1.468557e+05 1.539999e+05 1.613889e+05 1.656908e+05 + 1.693471e+05 1.760142e+05 1.812553e+05 1.873262e+05 1.984413e+05 + 2.015337e+05 2.106013e+05 2.119435e+05 2.196635e+05 2.285301e+05 + 2.368854e+05 2.422217e+05 2.543786e+05 2.633250e+05 2.706326e+05 + 2.750603e+05 2.910239e+05 3.045884e+05 3.130142e+05 3.188584e+05 + 3.285789e+05 3.399626e+05 3.507469e+05 3.637172e+05 3.741806e+05 + 3.904184e+05 4.008069e+05 4.254911e+05 4.293372e+05 4.523024e+05 + 4.622288e+05 4.807719e+05 4.938701e+05 5.078877e+05 5.231608e+05 + 5.412427e+05 5.589387e+05 5.817409e+05 5.971860e+05 6.122883e+05 + 6.425329e+05 6.618982e+05 6.809655e+05 7.039152e+05 7.354101e+05 + 7.549776e+05 7.766671e+05 8.026105e+05 8.429900e+05 8.740523e+05 + 8.908045e+05 9.191396e+05 9.662563e+05 1.001041e+06 1.030041e+06 + 1.060946e+06 1.089684e+06 1.150674e+06 1.164706e+06 1.200898e+06 + 1.279452e+06 1.271447e+06 1.333980e+06 1.390647e+06 1.440428e+06 + 1.471303e+06 1.539177e+06 1.586709e+06 1.634066e+06 1.695469e+06 + 1.741327e+06 1.806084e+06 1.869876e+06 1.933938e+06 1.979604e+06 + 2.093888e+06 2.166930e+06 2.265165e+06 2.294742e+06 2.354508e+06 + 2.464106e+06 2.520649e+06 2.670378e+06 2.722620e+06 2.768233e+06 + 2.867300e+06 3.003317e+06 3.028797e+06 3.179778e+06 3.309525e+06 + 3.465320e+06 3.576304e+06 3.655310e+06 3.846793e+06 3.923541e+06 + 4.092392e+06 4.248305e+06 4.518226e+06 4.487231e+06 4.671564e+06 + 4.996738e+06 5.108565e+06 5.315782e+06 5.570328e+06 5.744651e+06 + 5.872509e+06 5.987757e+06 6.341752e+06 6.470877e+06 6.772215e+06 + 6.946503e+06 7.278125e+06 7.640013e+06 7.714839e+06 8.039751e+06 + 8.494502e+06 8.869999e+06 9.033183e+06 9.488911e+06 9.766928e+06 + 1.022764e+07 1.049437e+07 1.083413e+07 1.103215e+07 1.157572e+07 + 1.212642e+07 1.268534e+07 1.294905e+07 1.330881e+07 1.407757e+07 + 1.441326e+07 1.519522e+07 1.535392e+07 1.634973e+07 1.694721e+07 + 1.747633e+07 1.814262e+07 1.910661e+07 1.957293e+07 2.046833e+07 + 2.118353e+07 2.226533e+07 2.303698e+07 2.394613e+07 2.531498e+07 + 2.605995e+07 2.715415e+07 2.846069e+07 2.953247e+07 3.073770e+07 + 3.225625e+07 3.393882e+07 3.470952e+07 3.674492e+07 3.897454e+07 + 4.005168e+07 4.103969e+07 4.364481e+07 4.538396e+07 4.664388e+07 + 4.986418e+07 5.144487e+07 5.514519e+07 5.669231e+07 5.826391e+07 + 6.109010e+07 6.432554e+07 6.770208e+07 7.007861e+07 7.545150e+07 + 7.780215e+07 8.083321e+07 8.390743e+07 8.773106e+07 9.210111e+07 + 9.338437e+07 9.794151e+07 1.017470e+08 1.055494e+08 1.079075e+08 + 1.089422e+08 1.128917e+08 1.072069e+08 9.669178e+07 7.435091e+07 + 4.410439e+07 1.706532e+07 4.223000e+06 6.719128e+05 5.827995e+04 + 2.447921e+03 4.470446e+01 6.579453e-01 1.203861e-03 1.592778e-06 + 8.984822e-10 3.613184e-14 2.100864e-19 3.859123e-29 0.000000e+00 + 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 +c approximate beam size from x and y projections +si23 0 4 +sp23 0 1 +si24 -6 6 +sp24 0 1