diff --git a/e2-estia.i b/e2-estia.i index 51101ae..e3f8378 100644 --- a/e2-estia.i +++ b/e2-estia.i @@ -108,10 +108,10 @@ c read file=tallies/bunker_mesh_np.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/inbunker_f4.i echo -c read file=tallies/chopper_pit_doses.i echo -c read file=tallies/bunker_wall_f4.i echo -c read file=tallies/e1_e2_f4.i echo +read file=tallies/inbunker_f4.i echo +read file=tallies/chopper_pit_doses.i echo +read file=tallies/bunker_wall_f4.i echo +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 c read file=tallies/shutter_f5.i echo @@ -144,7 +144,7 @@ c ------------------------------------------------------- c --------------- PHYSICS CARDS -------------------------- c ------------------------------------------------------- c dbcn 375642321 j j j j j j j j j j j 15291711 -nps 1e6 +nps 2e8 c histp mode n p n h / d t s a c volume calculated for tally cells (spheres) @@ -154,9 +154,6 @@ vol 1 506r 11j 1 183r 3848.5 10104 1j 1 7r 479.5 479.5 913 913 1734 imp:h,/,d,t,s,a 0 0 1 796r imp:n,p 0 0 1 796r cut:h,/,d,t,s,a 1e+08 1e-3 -c cut:n j 0.99E-6 -cut:n j 1e-11 -cut:p 1e+08 0.01 phys:n 3000 0 phys:p 100 phys:/,d,t,s,a 3000 @@ -170,9 +167,9 @@ 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 c * McStas source terms for gamma doeses -c read file=sources/mcstas_source_vs.i echo +read file=sources/mcstas_source_vs.i echo c read file=sources/mcstas_source_sample.i echo -read file=sources/mcstas_source_mf.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 prdmp 1e8 -60 1 2 1e7 diff --git a/geometry/estia_surfaces.i b/geometry/estia_surfaces.i index 09c2945..23e62cb 100644 --- a/geometry/estia_surfaces.i +++ b/geometry/estia_surfaces.i @@ -705,7 +705,7 @@ c Vacuum vessels before S1/S1-S2 8922 4 c/x 2.500 0.000 15.0 8931 3 c/x 0 -2 16.0 8932 3 c/x 0 -2 15.0 -8933 3 c/x 0 -2 17.5 +8933 3 c/x 0 -2 16.5 c Sphere at sample position to simulate gamma generation 9001 15 s 0 0 0 5.0 c Sphere at MF position to simulate gamma generation from mask @@ -896,50 +896,50 @@ c VS Carry Rods 25421 c/z 669.042 -892.275 3.000 25422 c/z 630.343 -896.070 3.000 25423 c/z 652.978 -927.687 3.000 -c c VS Blades (closed) -c 24501 p -0.860742 -0.509041 0.000000 -97.298 -c 24502 p -0.860742 -0.509041 0.000000 -99.298 -c 24503 p -0.860742 -0.509041 0.000000 -99.398 -c 24504 p 0.860742 0.509041 0.000000 101.798 -c 21501 pz 4.261 -c 21502 pz 0.261 -c 21503 pz -5.539 -c 21504 pz -1.539 -c 24505 p 0.509041 -0.860742 0.000000 1110.096 -c 24506 p 0.509041 -0.860742 0.000000 1110.496 -c c -c 24511 p -0.860742 -0.509041 0.000000 -96.807 -c 24512 p -0.860742 -0.509041 0.000000 -99.257 -c 24513 p -0.860742 -0.509041 0.000000 -101.307 -c 24514 p -0.860742 -0.509041 0.000000 -99.357 -c 21511 pz 6.061 -c 21512 pz 0.261 -c 21513 pz 2.561 -c 21514 pz -3.739 -c 24515 p -0.509041 0.860742 0.000000 -1110.595 -c 24516 p -0.509041 0.860742 0.000000 -1110.895 -c VS Blades (open) -24501 p -0.936060 -0.351842 0.000000 -288.420 -24502 p -0.936060 -0.351842 0.000000 -290.420 -24503 p -0.936060 -0.351842 0.000000 -290.520 -24504 p 0.936060 0.351842 0.000000 292.920 -21501 pz 5.061 -21502 pz 1.061 -21503 pz -4.739 -21504 pz -0.739 -24505 p 0.351842 -0.936060 0.000000 1073.540 -24506 p 0.351842 -0.936060 0.000000 1073.940 +c VS Blades (closed) +24501 p -0.860742 -0.509041 0.000000 -97.298 +24502 p -0.860742 -0.509041 0.000000 -99.298 +24503 p -0.860742 -0.509041 0.000000 -99.398 +24504 p 0.860742 0.509041 0.000000 101.798 +21501 pz 4.261 +21502 pz 0.261 +21503 pz -5.539 +21504 pz -1.539 +24505 p 0.509041 -0.860742 0.000000 1110.096 +24506 p 0.509041 -0.860742 0.000000 1110.496 c -24511 p -0.936060 -0.351842 0.000000 -288.356 -24512 p -0.936060 -0.351842 0.000000 -290.806 -24513 p -0.936060 -0.351842 0.000000 -292.856 -24514 p -0.936060 -0.351842 0.000000 -290.906 -21511 pz 5.261 -21512 pz 0.539 -21513 pz 1.761 -21514 pz -4.539 -24515 p -0.351842 0.936060 0.000000 -1078.921 -24516 p -0.351842 0.936060 0.000000 -1079.221 +24511 p -0.860742 -0.509041 0.000000 -96.807 +24512 p -0.860742 -0.509041 0.000000 -99.257 +24513 p -0.860742 -0.509041 0.000000 -101.307 +24514 p -0.860742 -0.509041 0.000000 -99.357 +21511 pz 6.061 +21512 pz 0.261 +21513 pz 2.561 +21514 pz -3.739 +24515 p -0.509041 0.860742 0.000000 -1110.595 +24516 p -0.509041 0.860742 0.000000 -1110.895 +c c VS Blades (open) +c 24501 p -0.936060 -0.351842 0.000000 -288.420 +c 24502 p -0.936060 -0.351842 0.000000 -290.420 +c 24503 p -0.936060 -0.351842 0.000000 -290.520 +c 24504 p 0.936060 0.351842 0.000000 292.920 +c 21501 pz 5.061 +c 21502 pz 1.061 +c 21503 pz -4.739 +c 21504 pz -0.739 +c 24505 p 0.351842 -0.936060 0.000000 1073.540 +c 24506 p 0.351842 -0.936060 0.000000 1073.940 +c c +c 24511 p -0.936060 -0.351842 0.000000 -288.356 +c 24512 p -0.936060 -0.351842 0.000000 -290.806 +c 24513 p -0.936060 -0.351842 0.000000 -292.856 +c 24514 p -0.936060 -0.351842 0.000000 -290.906 +c 21511 pz 5.261 +c 21512 pz 0.539 +c 21513 pz 1.761 +c 21514 pz -4.539 +c 24515 p -0.351842 0.936060 0.000000 -1078.921 +c 24516 p -0.351842 0.936060 0.000000 -1079.221 c visualize Selene gamma sources 50001 1 rpp 1339.25 2060.75 2 8 2 8 50002 1 rpp 2539.25 3260.75 -8 -2 -8 -2 diff --git a/geometry/materials.i b/geometry/materials.i index b8d14b8..b77fdf1 100644 --- a/geometry/materials.i +++ b/geometry/materials.i @@ -3,7 +3,7 @@ c --------------- MATERIAL CARDS ------------------------ c ------------------------------------------------------- c c Gd3Ga5O12 -7.08 (Gadolinium Gallium Garnet - GGG) -m2 6000.60c 0.6 31000.60c 0.25 64000.35c 0.15 +m2 6000.70c 0.6 31000.60c 0.25 64000.35c 0.15 c c Material : Stainless304 rho=0.0876276 m3 6000.70c 0.00031864 14028.70c 0.00170336 15031.70c 6.95038e-05 @@ -48,7 +48,7 @@ m9 2004.70c 2.45e-05 pnlib=70u c c BorAL 1100 alloy (2.7 g/cc) m10 5010.70c -0.04 5011.70c -0.01 13027.70c -0.9288 & - 14000.70c -0.01 25055.70c -5.0E-3 26000.50c -5.0E-3 & + 14028.70c -0.01 25055.70c -5.0E-3 26000.50c -5.0E-3 & 29000.50c -1.0E-3 30000.70c -2.0E-4 c c Material : H2O rho=0.100283 @@ -57,7 +57,7 @@ m11 1001.70c 0.0668456 1002.70c 1.003e-05 8016.70c 0.0334278 mt11 lwtr.01t c c bor-float glass 13% B2O3 2.51 g/cm3 -m12 8016.70c -0.4078 11023.70c -0.096441 14000.70c -0.336563 & +m12 8016.70c -0.4078 11023.70c -0.096441 14028.70c -0.336563 & 20000.60c -0.107205 5010.70c -0.0104 5011.70c -0.0416 c lead (11.34 g/cm3) m13 82000.50c 1.0 @@ -69,7 +69,7 @@ c 6000.70c -5.23e-3 c 8016.70c -0.08511 c 12000.60c -7e-4 c 13027.70c -4.71e-3 -c 14000.70c -1.983e-2 +c 14028.70c -1.983e-2 c 16000.70c -1.368e-2 c 20000.60c -2.425e-2 c 24000.50c -0.13896 @@ -88,7 +88,7 @@ m16 1001.70c -6.09e-3 8016.70c -0.404 12000.60c -0.0322 13027.70c -0.0115 - 14000.70c -3.3906e-2 + 14028.70c -3.3906e-2 20000.60c -0.3222 26054.70c -4.22e-4 26056.70c -6.92e-3 @@ -111,7 +111,7 @@ m18 1001.70c -0.05 8016.70c -0.2325 12000.60c -0.007 13027.70c -0.023 - 14000.70c -0.38 + 14028.70c -0.38 26000.50c -0.02 28000.50c -0.004 5010.70c -0.0047 @@ -205,5 +205,6 @@ m200 01001.70c 1.3431e-01 01002.70c 1.5448e-05 08016.70c 5.2962e-01 mt200 lwtr.01t c c Ar 2.7bar, CO2 0.1bar, 3He 0.2bar -m300 2003.70c 0.02381 18000.70c 0.01190 - 6000.70c 0.32143 8000.70c 0.64286 +m300 2003.70c 0.02381 18036.70C 0.00004046 + 18038.70C 0.0000075208 18040.70C 0.0118524 + 6000.70c 0.32143 8016.70c 0.64286 diff --git a/sources/e2_2m.i b/sources/e2_2m.i index d84ab58..8339fe9 100644 --- a/sources/e2_2m.i +++ b/sources/e2_2m.i @@ -1,4 +1,6 @@ c source definition based on tally_181.txt +cut:n j 1e-11 +cut:p 1e+08 0.01 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 5 j -1 j 1.79 diff --git a/sources/mcstas_gammas.i b/sources/mcstas_gammas.i index a56b283..a48e014 100644 --- a/sources/mcstas_gammas.i +++ b/sources/mcstas_gammas.i @@ -1,4 +1,6 @@ c source definition based on McStas +cut:n j 1e-11 +cut:p 1e+08 0.01 c make sure to supply the correct WWINP file WWP:n j j 5 j -1 j 1.0 c the source definition with gammas/s diff --git a/sources/mcstas_gammas_2.i b/sources/mcstas_gammas_2.i index 1dad2c3..b4b84cb 100644 --- a/sources/mcstas_gammas_2.i +++ b/sources/mcstas_gammas_2.i @@ -1,4 +1,6 @@ c source definition based on McStas +cut:n j 1e-11 +cut:p 1e+08 0.01 c make sure to supply the correct WWINP file WWP:n j j 5 j -1 j 1.0 c the source definition with gammas/s diff --git a/sources/mcstas_source_mf.i b/sources/mcstas_source_mf.i index d41e790..37ca49b 100644 --- a/sources/mcstas_source_mf.i +++ b/sources/mcstas_source_mf.i @@ -1,4 +1,8 @@ c cold neutron source from McStas simulation for sample position +c activation (delayed particles) +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=-5.1 y=d23 z=d24 vec=1 0 0 tr=13 @@ -6,48 +10,122 @@ sdef dir=1 erg=d11 x=-5.1 y=d23 z=d24 vec=1 0 0 tr=13 c scaling wgt is to flux from McStas nps si1 -1 0.99976203 1 sp1 0 0 1 -si11 0.00000e+00 3.99924e-04 4.03684e-04 4.07480e-04 4.11347e-04 - 4.15289e-04 4.19307e-04 4.23404e-04 4.27582e-04 4.31842e-04 - 4.36189e-04 4.40624e-04 4.45150e-04 4.49770e-04 4.54487e-04 - 4.59303e-04 4.64224e-04 4.69250e-04 4.74387e-04 4.79637e-04 - 4.85005e-04 4.90494e-04 4.96110e-04 5.01855e-04 5.07734e-04 - 5.13754e-04 5.19917e-04 5.26231e-04 5.32699e-04 5.39329e-04 - 5.46125e-04 5.53095e-04 5.60246e-04 5.67583e-04 5.75116e-04 - 5.82851e-04 5.90797e-04 5.98962e-04 6.07357e-04 6.15990e-04 - 6.24872e-04 6.34014e-04 6.43428e-04 6.53125e-04 6.63119e-04 - 6.73424e-04 6.84054e-04 6.95025e-04 7.06354e-04 7.18058e-04 - 7.30157e-04 7.42670e-04 7.55620e-04 7.69030e-04 7.82924e-04 - 7.97329e-04 8.12274e-04 8.27791e-04 8.43912e-04 8.60673e-04 - 8.78114e-04 8.96276e-04 9.15205e-04 9.34952e-04 9.55569e-04 - 9.77117e-04 9.99658e-04 1.02326e-03 1.04801e-03 1.07399e-03 - 1.10128e-03 1.13000e-03 1.16026e-03 1.19219e-03 1.22592e-03 - 1.26161e-03 1.29945e-03 1.33962e-03 1.38237e-03 1.42792e-03 - 1.47659e-03 1.52869e-03 1.58460e-03 1.64476e-03 1.70967e-03 - 1.77992e-03 1.85619e-03 1.93929e-03 2.03018e-03 2.13002e-03 - 2.24019e-03 2.36239e-03 2.49870e-03 2.65173e-03 2.82474e-03 - 3.02193e-03 3.24876e-03 3.51246e-03 3.82285e-03 4.19354e-03 - 4.64408e-03 5.14032e-03 -sp11 0.00000e+00 0.00000e+00 3.85593e+04 4.20865e+04 4.58640e+04 - 5.17596e+04 5.61391e+04 6.11107e+04 6.63908e+04 7.21455e+04 - 8.09129e+04 8.94953e+04 9.82641e+04 1.09883e+05 1.15230e+05 - 1.29230e+05 1.38110e+05 1.54438e+05 1.69936e+05 1.91156e+05 - 2.06688e+05 2.27873e+05 2.47183e+05 2.76754e+05 2.99902e+05 - 3.25928e+05 3.72806e+05 4.08736e+05 4.29439e+05 4.81832e+05 - 5.27307e+05 5.61325e+05 6.25344e+05 6.99932e+05 7.59986e+05 - 8.23687e+05 9.09114e+05 9.92801e+05 1.08652e+06 1.22423e+06 - 1.31508e+06 1.42046e+06 1.55452e+06 1.73926e+06 1.95625e+06 - 2.08401e+06 2.32159e+06 2.53700e+06 2.79880e+06 3.01669e+06 - 3.34322e+06 3.66624e+06 4.05774e+06 4.40861e+06 4.82292e+06 - 5.45131e+06 6.05936e+06 6.43475e+06 7.37179e+06 7.70977e+06 - 8.68124e+06 9.22544e+06 1.06175e+07 1.13642e+07 1.25066e+07 - 1.41004e+07 1.53690e+07 1.63458e+07 1.90828e+07 2.02059e+07 - 2.28912e+07 2.63899e+07 2.77627e+07 3.17059e+07 3.37525e+07 - 3.78514e+07 4.18388e+07 4.58109e+07 5.03699e+07 5.57842e+07 - 6.23562e+07 6.85334e+07 7.62412e+07 8.32321e+07 8.98834e+07 - 1.00885e+08 1.14599e+08 1.23979e+08 1.40106e+08 1.47047e+08 - 1.73630e+08 1.90716e+08 2.03589e+08 2.22960e+08 2.43119e+08 - 2.54674e+08 2.27755e+08 1.14912e+08 3.88323e+06 6.21135e+01 - 1.34560e-07 1.83476e-18 +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 -1.1 -1.085859e+00 -1.035354e+00 -9.848485e-01 -9.343434e-01 -8.838384e-01 diff --git a/sources/mcstas_source_sample.i b/sources/mcstas_source_sample.i index 61bf4b2..b29a8f8 100644 --- a/sources/mcstas_source_sample.i +++ b/sources/mcstas_source_sample.i @@ -1,4 +1,8 @@ c cold neutron source from McStas simulation for sample position +c activation (delayed particles) +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=-5.1 y=d23 z=d24 vec=1 0 0 tr=15 @@ -6,48 +10,121 @@ sdef dir=1 erg=d11 x=-5.1 y=d23 z=d24 vec=1 0 0 tr=15 c scaling wgt is to flux from McStas nps si1 -1 0.99976203 1 sp1 0 0 1 -si11 0.00000e+00 3.99924e-04 4.03684e-04 4.07480e-04 4.11347e-04 - 4.15289e-04 4.19307e-04 4.23404e-04 4.27582e-04 4.31842e-04 - 4.36189e-04 4.40624e-04 4.45150e-04 4.49770e-04 4.54487e-04 - 4.59303e-04 4.64224e-04 4.69250e-04 4.74387e-04 4.79637e-04 - 4.85005e-04 4.90494e-04 4.96110e-04 5.01855e-04 5.07734e-04 - 5.13754e-04 5.19917e-04 5.26231e-04 5.32699e-04 5.39329e-04 - 5.46125e-04 5.53095e-04 5.60246e-04 5.67583e-04 5.75116e-04 - 5.82851e-04 5.90797e-04 5.98962e-04 6.07357e-04 6.15990e-04 - 6.24872e-04 6.34014e-04 6.43428e-04 6.53125e-04 6.63119e-04 - 6.73424e-04 6.84054e-04 6.95025e-04 7.06354e-04 7.18058e-04 - 7.30157e-04 7.42670e-04 7.55620e-04 7.69030e-04 7.82924e-04 - 7.97329e-04 8.12274e-04 8.27791e-04 8.43912e-04 8.60673e-04 - 8.78114e-04 8.96276e-04 9.15205e-04 9.34952e-04 9.55569e-04 - 9.77117e-04 9.99658e-04 1.02326e-03 1.04801e-03 1.07399e-03 - 1.10128e-03 1.13000e-03 1.16026e-03 1.19219e-03 1.22592e-03 - 1.26161e-03 1.29945e-03 1.33962e-03 1.38237e-03 1.42792e-03 - 1.47659e-03 1.52869e-03 1.58460e-03 1.64476e-03 1.70967e-03 - 1.77992e-03 1.85619e-03 1.93929e-03 2.03018e-03 2.13002e-03 - 2.24019e-03 2.36239e-03 2.49870e-03 2.65173e-03 2.82474e-03 - 3.02193e-03 3.24876e-03 3.51246e-03 3.82285e-03 4.19354e-03 - 4.64408e-03 5.14032e-03 -sp11 0.00000e+00 0.00000e+00 3.85593e+04 4.20865e+04 4.58640e+04 - 5.17596e+04 5.61391e+04 6.11107e+04 6.63908e+04 7.21455e+04 - 8.09129e+04 8.94953e+04 9.82641e+04 1.09883e+05 1.15230e+05 - 1.29230e+05 1.38110e+05 1.54438e+05 1.69936e+05 1.91156e+05 - 2.06688e+05 2.27873e+05 2.47183e+05 2.76754e+05 2.99902e+05 - 3.25928e+05 3.72806e+05 4.08736e+05 4.29439e+05 4.81832e+05 - 5.27307e+05 5.61325e+05 6.25344e+05 6.99932e+05 7.59986e+05 - 8.23687e+05 9.09114e+05 9.92801e+05 1.08652e+06 1.22423e+06 - 1.31508e+06 1.42046e+06 1.55452e+06 1.73926e+06 1.95625e+06 - 2.08401e+06 2.32159e+06 2.53700e+06 2.79880e+06 3.01669e+06 - 3.34322e+06 3.66624e+06 4.05774e+06 4.40861e+06 4.82292e+06 - 5.45131e+06 6.05936e+06 6.43475e+06 7.37179e+06 7.70977e+06 - 8.68124e+06 9.22544e+06 1.06175e+07 1.13642e+07 1.25066e+07 - 1.41004e+07 1.53690e+07 1.63458e+07 1.90828e+07 2.02059e+07 - 2.28912e+07 2.63899e+07 2.77627e+07 3.17059e+07 3.37525e+07 - 3.78514e+07 4.18388e+07 4.58109e+07 5.03699e+07 5.57842e+07 - 6.23562e+07 6.85334e+07 7.62412e+07 8.32321e+07 8.98834e+07 - 1.00885e+08 1.14599e+08 1.23979e+08 1.40106e+08 1.47047e+08 - 1.73630e+08 1.90716e+08 2.03589e+08 2.22960e+08 2.43119e+08 - 2.54674e+08 2.27755e+08 1.14912e+08 3.88323e+06 6.21135e+01 - 1.34560e-07 1.83476e-18 +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 1.008025e+04 1.044972e+04 1.080195e+04 1.099551e+04 1.142110e+04 + 1.181709e+04 1.225534e+04 1.265653e+04 1.318484e+04 1.370127e+04 + 1.398173e+04 1.444153e+04 1.484467e+04 1.555741e+04 1.604779e+04 + 1.670984e+04 1.726319e+04 1.751394e+04 1.843601e+04 1.875144e+04 + 1.950671e+04 1.983854e+04 2.102559e+04 2.130872e+04 2.256400e+04 + 2.323411e+04 2.422479e+04 2.454622e+04 2.576871e+04 2.617331e+04 + 2.714874e+04 2.822352e+04 2.861572e+04 3.009640e+04 3.096840e+04 + 3.189114e+04 3.310325e+04 3.435932e+04 3.561572e+04 3.715973e+04 + 3.800929e+04 3.925349e+04 4.016106e+04 4.122701e+04 4.330927e+04 + 4.504885e+04 4.625217e+04 4.773717e+04 4.855331e+04 5.005798e+04 + 5.340498e+04 5.366403e+04 5.580346e+04 5.827668e+04 5.980446e+04 + 6.156877e+04 6.441039e+04 6.589969e+04 6.842061e+04 7.177448e+04 + 7.391729e+04 7.531638e+04 7.798407e+04 7.931816e+04 8.240248e+04 + 8.578720e+04 8.981466e+04 9.303574e+04 9.575726e+04 9.891288e+04 + 1.020687e+05 1.050539e+05 1.107794e+05 1.124304e+05 1.173669e+05 + 1.221830e+05 1.226280e+05 1.273947e+05 1.349663e+05 1.382150e+05 + 1.403006e+05 1.466681e+05 1.508025e+05 1.548290e+05 1.643740e+05 + 1.678966e+05 1.738443e+05 1.772489e+05 1.829497e+05 1.894190e+05 + 1.975155e+05 2.020876e+05 2.114231e+05 2.185840e+05 2.248148e+05 + 2.285961e+05 2.435891e+05 2.526266e+05 2.590725e+05 2.670361e+05 + 2.739711e+05 2.823097e+05 2.920449e+05 3.042509e+05 3.131360e+05 + 3.243324e+05 3.358441e+05 3.518691e+05 3.590528e+05 3.757659e+05 + 3.840672e+05 3.964595e+05 4.109730e+05 4.229469e+05 4.303120e+05 + 4.536527e+05 4.679909e+05 4.863041e+05 4.946692e+05 5.112674e+05 + 5.337627e+05 5.489484e+05 5.651495e+05 5.796258e+05 6.096617e+05 + 6.229919e+05 6.476939e+05 6.690880e+05 7.059023e+05 7.304116e+05 + 7.422505e+05 7.690611e+05 8.046053e+05 8.331856e+05 8.555723e+05 + 8.867643e+05 9.045859e+05 9.523897e+05 9.669171e+05 1.002110e+06 + 1.061782e+06 1.060013e+06 1.105509e+06 1.165756e+06 1.199584e+06 + 1.220154e+06 1.283135e+06 1.326009e+06 1.358438e+06 1.416615e+06 + 1.452558e+06 1.498097e+06 1.565796e+06 1.607302e+06 1.653762e+06 + 1.746285e+06 1.805085e+06 1.891255e+06 1.915608e+06 1.960642e+06 + 2.055185e+06 2.096000e+06 2.207035e+06 2.274029e+06 2.296175e+06 + 2.390864e+06 2.510334e+06 2.527698e+06 2.642242e+06 2.778158e+06 + 2.870501e+06 2.989092e+06 3.037544e+06 3.206918e+06 3.252058e+06 + 3.413884e+06 3.516998e+06 3.760520e+06 3.754380e+06 3.866644e+06 + 4.163490e+06 4.217836e+06 4.427571e+06 4.646819e+06 4.748602e+06 + 4.832024e+06 4.999095e+06 5.251591e+06 5.435756e+06 5.662045e+06 + 5.797489e+06 6.035072e+06 6.314328e+06 6.399095e+06 6.620023e+06 + 7.003929e+06 7.339109e+06 7.522915e+06 7.817125e+06 8.045639e+06 + 8.372159e+06 8.604745e+06 8.862432e+06 8.919529e+06 9.447464e+06 + 9.895489e+06 1.024763e+07 1.042107e+07 1.081681e+07 1.123139e+07 + 1.148303e+07 1.208429e+07 1.223955e+07 1.297488e+07 1.337273e+07 + 1.367657e+07 1.417224e+07 1.480856e+07 1.524734e+07 1.591191e+07 + 1.634047e+07 1.695270e+07 1.762327e+07 1.816128e+07 1.892174e+07 + 1.968160e+07 2.032581e+07 2.101324e+07 2.173065e+07 2.262567e+07 + 2.353282e+07 2.458781e+07 2.508849e+07 2.624155e+07 2.765789e+07 + 2.814961e+07 2.864851e+07 3.048953e+07 3.158131e+07 3.205461e+07 + 3.385151e+07 3.460774e+07 3.673969e+07 3.788368e+07 3.829032e+07 + 3.991199e+07 4.123477e+07 4.299071e+07 4.371126e+07 4.684197e+07 + 4.798410e+07 4.913580e+07 5.042680e+07 5.164860e+07 5.405405e+07 + 5.348658e+07 5.486452e+07 5.594011e+07 5.675180e+07 5.753822e+07 + 5.614008e+07 5.634224e+07 5.035996e+07 4.056106e+07 2.424032e+07 + 8.199239e+06 1.213338e+06 5.954076e+04 8.846304e+02 3.862136e+00 + 4.347077e-03 2.869644e-06 1.210559e-09 2.320094e-13 2.105981e-17 + 4.665008e-22 6.516805e-29 6.214276e-38 0.000000e+00 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.75 -0.73232323 -0.68181818 -0.63131313 -0.58080808 -0.53030303 diff --git a/sources/mcstas_source_vs.i b/sources/mcstas_source_vs.i index d37a702..5300726 100644 --- a/sources/mcstas_source_vs.i +++ b/sources/mcstas_source_vs.i @@ -1,4 +1,8 @@ c cold neutron source from McStas simulation for VS position +c activation (delayed particles) +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=-26.25 y=d23 z=d24 vec=1 0 0 tr=4 @@ -6,48 +10,122 @@ sdef dir=1 erg=d11 x=-26.25 y=d23 z=d24 vec=1 0 0 tr=4 c scaling wgt is to flux from McStas nps si1 -1 0.99976203 1 sp1 0 0 1 -si11 0.00000e+00 3.99922e-04 4.03763e-04 4.07641e-04 4.11594e-04 - 4.15624e-04 4.19734e-04 4.23927e-04 4.28204e-04 4.32568e-04 - 4.37022e-04 4.41568e-04 4.46211e-04 4.50951e-04 4.55794e-04 - 4.60742e-04 4.65799e-04 4.70967e-04 4.76252e-04 4.81657e-04 - 4.87185e-04 4.92842e-04 4.98632e-04 5.04560e-04 5.10630e-04 - 5.16849e-04 5.23220e-04 5.29751e-04 5.36446e-04 5.43314e-04 - 5.50359e-04 5.57589e-04 5.65012e-04 5.72635e-04 5.80467e-04 - 5.88516e-04 5.96791e-04 6.05302e-04 6.14060e-04 6.23074e-04 - 6.32358e-04 6.41922e-04 6.51780e-04 6.61945e-04 6.72433e-04 - 6.83258e-04 6.94438e-04 7.05989e-04 7.17932e-04 7.30285e-04 - 7.43071e-04 7.56313e-04 7.70035e-04 7.84264e-04 7.99030e-04 - 8.14362e-04 8.30294e-04 8.46861e-04 8.64104e-04 8.82063e-04 - 9.00785e-04 9.20319e-04 9.40719e-04 9.62044e-04 9.84358e-04 - 1.00773e-03 1.03224e-03 1.05798e-03 1.08503e-03 1.11350e-03 - 1.14350e-03 1.17517e-03 1.20864e-03 1.24407e-03 1.28164e-03 - 1.32156e-03 1.36404e-03 1.40935e-03 1.45776e-03 1.50963e-03 - 1.56532e-03 1.62528e-03 1.69002e-03 1.76013e-03 1.83631e-03 - 1.91940e-03 2.01036e-03 2.11037e-03 2.22086e-03 2.34358e-03 - 2.48066e-03 2.63479e-03 2.80936e-03 3.00874e-03 3.23862e-03 - 3.50661e-03 3.82304e-03 4.20239e-03 4.66558e-03 5.24394e-03 - 5.98680e-03 6.82707e-03 -sp11 0.00000e+00 0.00000e+00 4.73317e+05 5.20458e+05 5.71528e+05 - 6.28807e+05 6.92966e+05 7.71888e+05 8.44652e+05 9.36465e+05 - 1.00925e+06 1.11561e+06 1.22453e+06 1.35942e+06 1.50732e+06 - 1.61141e+06 1.78013e+06 1.97405e+06 2.18967e+06 2.41912e+06 - 2.62402e+06 2.90044e+06 3.19088e+06 3.52382e+06 3.90203e+06 - 4.20625e+06 4.69078e+06 5.14001e+06 5.62174e+06 6.21576e+06 - 6.86171e+06 7.54765e+06 8.23914e+06 9.04558e+06 9.85563e+06 - 1.10350e+07 1.19012e+07 1.32544e+07 1.46700e+07 1.59696e+07 - 1.78449e+07 1.93294e+07 2.11647e+07 2.34912e+07 2.62256e+07 - 2.83634e+07 3.12576e+07 3.46387e+07 3.76342e+07 4.14516e+07 - 4.51023e+07 5.08468e+07 5.49930e+07 6.13677e+07 6.70328e+07 - 7.41103e+07 8.05813e+07 9.01005e+07 9.81536e+07 1.09147e+08 - 1.20167e+08 1.33930e+08 1.45450e+08 1.58317e+08 1.79033e+08 - 1.94626e+08 2.16941e+08 2.41613e+08 2.65787e+08 2.97746e+08 - 3.32680e+08 3.73719e+08 4.09773e+08 4.63621e+08 5.10589e+08 - 5.68804e+08 6.48527e+08 7.26139e+08 8.17297e+08 9.35034e+08 - 1.04309e+09 1.18888e+09 1.36129e+09 1.55131e+09 1.79901e+09 - 2.09372e+09 2.40061e+09 2.78880e+09 3.30983e+09 3.77655e+09 - 4.46366e+09 5.25730e+09 6.13953e+09 7.09792e+09 8.02417e+09 - 8.04221e+09 5.96527e+09 2.49280e+09 3.64107e+08 1.47476e+07 - 4.45917e+03 3.90976e-04 +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 + 7.314213e+04 7.534749e+04 7.770641e+04 8.071539e+04 8.293751e+04 + 8.594698e+04 8.901406e+04 9.176514e+04 9.558657e+04 9.804128e+04 + 1.019017e+05 1.049198e+05 1.082036e+05 1.119810e+05 1.161518e+05 + 1.204111e+05 1.244043e+05 1.279849e+05 1.325754e+05 1.373407e+05 + 1.417398e+05 1.447561e+05 1.520370e+05 1.556986e+05 1.621075e+05 + 1.667239e+05 1.719704e+05 1.796481e+05 1.842749e+05 1.915267e+05 + 1.959106e+05 2.038436e+05 2.091362e+05 2.173824e+05 2.250166e+05 + 2.322025e+05 2.405742e+05 2.484920e+05 2.551448e+05 2.651063e+05 + 2.738210e+05 2.836766e+05 2.930146e+05 3.016496e+05 3.134208e+05 + 3.241464e+05 3.347430e+05 3.453339e+05 3.523625e+05 3.677829e+05 + 3.832107e+05 3.940389e+05 4.068824e+05 4.229935e+05 4.368460e+05 + 4.510380e+05 4.691019e+05 4.810431e+05 4.960731e+05 5.143496e+05 + 5.290191e+05 5.481243e+05 5.673071e+05 5.877119e+05 6.059504e+05 + 6.268908e+05 6.490568e+05 6.736975e+05 6.936052e+05 7.134860e+05 + 7.438258e+05 7.645571e+05 7.920251e+05 8.170812e+05 8.452803e+05 + 8.741586e+05 8.954765e+05 9.299146e+05 9.672485e+05 9.964529e+05 + 1.031042e+06 1.066913e+06 1.091902e+06 1.127949e+06 1.177567e+06 + 1.211853e+06 1.250921e+06 1.297242e+06 1.333577e+06 1.390353e+06 + 1.425827e+06 1.476786e+06 1.544936e+06 1.583806e+06 1.622074e+06 + 1.675392e+06 1.745505e+06 1.815346e+06 1.866180e+06 1.927402e+06 + 1.987326e+06 2.061633e+06 2.121317e+06 2.193759e+06 2.262659e+06 + 2.344964e+06 2.434513e+06 2.520442e+06 2.594007e+06 2.697823e+06 + 2.759750e+06 2.855345e+06 2.968563e+06 3.052496e+06 3.157187e+06 + 3.259902e+06 3.382436e+06 3.491134e+06 3.641182e+06 3.730913e+06 + 3.850323e+06 3.980167e+06 4.119078e+06 4.256675e+06 4.408269e+06 + 4.566608e+06 4.709554e+06 4.861482e+06 5.034308e+06 5.194571e+06 + 5.381048e+06 5.568007e+06 5.730777e+06 5.921100e+06 6.144306e+06 + 6.339086e+06 6.604239e+06 6.841518e+06 7.031530e+06 7.267075e+06 + 7.535507e+06 7.766925e+06 8.005327e+06 8.275115e+06 8.603228e+06 + 8.904485e+06 9.191018e+06 9.515403e+06 9.804721e+06 1.013352e+07 + 1.051150e+07 1.090381e+07 1.128360e+07 1.168800e+07 1.196222e+07 + 1.235378e+07 1.279619e+07 1.341781e+07 1.379545e+07 1.419294e+07 + 1.482206e+07 1.516191e+07 1.582630e+07 1.631618e+07 1.680888e+07 + 1.737588e+07 1.805657e+07 1.866849e+07 1.946390e+07 2.004097e+07 + 2.074718e+07 2.146861e+07 2.227579e+07 2.306144e+07 2.391859e+07 + 2.459553e+07 2.569143e+07 2.675689e+07 2.745817e+07 2.841127e+07 + 2.960418e+07 3.042320e+07 3.174439e+07 3.306165e+07 3.419572e+07 + 3.548158e+07 3.652932e+07 3.812443e+07 3.928993e+07 4.075788e+07 + 4.244747e+07 4.399839e+07 4.544803e+07 4.759183e+07 4.911309e+07 + 5.113628e+07 5.335918e+07 5.506053e+07 5.701091e+07 5.942090e+07 + 6.202103e+07 6.403282e+07 6.708073e+07 6.934622e+07 7.219234e+07 + 7.491512e+07 7.812413e+07 8.089741e+07 8.465862e+07 8.787138e+07 + 9.175372e+07 9.626414e+07 9.888579e+07 1.033984e+08 1.083084e+08 + 1.134365e+08 1.168989e+08 1.234029e+08 1.282685e+08 1.328512e+08 + 1.389779e+08 1.463445e+08 1.522384e+08 1.605338e+08 1.675247e+08 + 1.758617e+08 1.832638e+08 1.927654e+08 2.003541e+08 2.118512e+08 + 2.234568e+08 2.342977e+08 2.453942e+08 2.572761e+08 2.729111e+08 + 2.862787e+08 3.023340e+08 3.187229e+08 3.298927e+08 3.513120e+08 + 3.685151e+08 3.904085e+08 4.115645e+08 4.337881e+08 4.588002e+08 + 4.840627e+08 5.104945e+08 5.377827e+08 5.674106e+08 6.003238e+08 + 6.381702e+08 6.705987e+08 7.056909e+08 7.455516e+08 7.926358e+08 + 8.282402e+08 8.779378e+08 9.219139e+08 9.752022e+08 1.022091e+09 + 1.058257e+09 1.113515e+09 1.117680e+09 1.125860e+09 1.072965e+09 + 9.951232e+08 8.589051e+08 7.065817e+08 5.474185e+08 3.787102e+08 + 2.327046e+08 1.220900e+08 5.516838e+07 2.114656e+07 6.556119e+06 + 1.700399e+06 2.348269e+05 1.308566e+04 3.135864e+02 1.675656e+00 + 8.780346e-04 3.404936e-06 2.293624e-09 1.910527e-12 c approximate size of beam at VS is 40mmx20mm si23 -2.2 -2.196970e+00 -2.146465e+00 -2.095960e+00 -2.045455e+00 -1.994949e+00 diff --git a/sources/proton.i b/sources/proton.i index d267d14..509caa4 100644 --- a/sources/proton.i +++ b/sources/proton.i @@ -1,4 +1,6 @@ c Original definition of source from proton beam +cut:n j 0.99e-06 +cut:p 1e+08 0.5 c weight windows from external file generated by WWG c make sure to supply the correct WWINP file WWP:n,p j j 10 j -1 j 8.7374e15 $ WW was generated for source scaling 1