From cedc4cee19dde1a2d85f6bd063e3753abda7f874 Mon Sep 17 00:00:00 2001 From: Artur Glavic Date: Wed, 16 Oct 2019 15:58:29 +0200 Subject: [PATCH] Update gamma source for selene 1+2 and add fine mesh to test source --- e2-estia.i | 13 +- sources/mcstas_gammas.i | 2 +- sources/mcstas_gammas_2.i | 465 +++++++++++++++++++++++++++++---- tallies/selene_source_mesh_p.i | 38 +++ 4 files changed, 456 insertions(+), 62 deletions(-) create mode 100644 tallies/selene_source_mesh_p.i diff --git a/e2-estia.i b/e2-estia.i index 1988b28..a0afbc4 100644 --- a/e2-estia.i +++ b/e2-estia.i @@ -142,20 +142,21 @@ 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_p.i echo +c read file=tallies/full_mesh_p.i echo c read file=tallies/full_mesh_np.i echo c read file=tallies/boundary_meshes.i echo c read file=tallies/boundary_wall_meshes.i echo +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 c read file=tallies/inbunker_f4.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 +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 c read file=tallies/shutter_f5.i echo @@ -197,7 +198,7 @@ c --------------- PHYSICS CARDS -------------------------- 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 -stop nps 2e8 CTME 4e5 $f31 0.01 +stop nps 1e6 CTME 4e5 $f31 0.01 rand gen=2 seed=29976355590901 $ default 19073486328125 c histp mode n p h / z $ neutron, photon, proton, pi+, pi0 @@ -236,8 +237,8 @@ c * McStas source terms for gamma doeses c read file=sources/mcstas_source_vs.i echo c read file=sources/mcstas_source_sample.i echo c read file=sources/mcstas_source_mf.i echo -read file=sources/mcstas_gammas.i echo -c read file=sources/mcstas_gammas_2.i echo +c read file=sources/mcstas_gammas.i echo +read file=sources/mcstas_gammas_2.i echo c read file=sources/mcstas_source_shutter.i echo prdmp 5e6 -30 1 2 5e5 print 30 32 40 diff --git a/sources/mcstas_gammas.i b/sources/mcstas_gammas.i index 1e03931..bbff4ab 100644 --- a/sources/mcstas_gammas.i +++ b/sources/mcstas_gammas.i @@ -1,4 +1,4 @@ -c source definition based on McStas +c source definition for Selene 1 gammas based on McStas cut:n j 1e-11 cut:p 1e+08 0.01 c make sure to supply the correct WWINP file diff --git a/sources/mcstas_gammas_2.i b/sources/mcstas_gammas_2.i index f72d076..7fef527 100644 --- a/sources/mcstas_gammas_2.i +++ b/sources/mcstas_gammas_2.i @@ -1,164 +1,519 @@ -c source definition based on McStas +c source definition for Selene 1+2 gammas 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 -sdef erg=fx d1 x=d31 y=d32 z=d33 tr=1 par=p wgt=1.04957e+09 -ds1 q 2539.25 100 - 2587.35 101 2635.45 102 2683.55 103 2731.65 104 2779.75 105 - 2827.85 106 2875.95 107 2924.05 108 2972.15 109 3020.25 110 - 3068.35 111 3116.45 112 3164.55 113 3212.65 114 3260.75 115 -c Energy distribution for all positions along the guide -si100 0 1 -sp100 0 1e-20 -si101 1e-12 0.15 0.2 0.3 0.4 +c the source definition with total gammas/s weight of Selene 1 and 2 +sdef pos=d1 $ dummy for X-position and horizontal/vertical+guide selection + x=fpos d11 y=fpos d12 z=fpos d13 + erg=fpos d14 + tr=1 par=p wgt=1.06479e+10 +si1 l 1 0 0 2 0 0 3 0 0 4 0 0 5 0 0 + 6 0 0 7 0 0 8 0 0 9 0 0 10 0 0 + 11 0 0 12 0 0 13 0 0 14 0 0 15 0 0 + 1 1 0 2 1 0 3 1 0 4 1 0 5 1 0 + 6 1 0 7 1 0 8 1 0 9 1 0 10 1 0 + 11 1 0 12 1 0 13 1 0 14 1 0 15 1 0 + 1 0 1 2 0 1 3 0 1 4 0 1 5 0 1 + 6 0 1 7 0 1 8 0 1 9 0 1 10 0 1 + 11 0 1 12 0 1 13 0 1 14 0 1 15 0 1 + 1 1 1 2 1 1 3 1 1 4 1 1 5 1 1 + 6 1 1 7 1 1 8 1 1 9 1 1 10 1 1 + 11 1 1 12 1 1 13 1 1 14 1 1 15 1 1 +sp1 7.29220E+09 4.79340E+08 3.21180E+08 3.03858E+08 2.58578E+08 $ S1V + 2.25073E+08 1.96544E+08 1.27443E+08 1.04253E+08 6.00576E+07 + 5.56826E+07 5.57082E+07 4.47933E+07 4.57106E+07 2.79109E+07 + 7.29220E+09 4.79340E+08 3.21180E+08 3.03858E+08 2.58578E+08 $ S1H + 2.25073E+08 1.96544E+08 1.27443E+08 1.04253E+08 6.00576E+07 + 5.56826E+07 5.57082E+07 4.47933E+07 4.57106E+07 2.79109E+07 + 4.57688E+07 2.25345E+07 2.38501E+07 3.55841E+07 3.75155E+07 $ S2V + 4.40951E+07 5.81307E+07 5.63921E+07 8.42501E+07 9.90333E+07 + 1.07862E+08 1.15949E+08 9.38290E+07 1.14450E+08 1.10325E+08 + 4.57688E+07 2.25345E+07 2.38501E+07 3.55841E+07 3.75155E+07 $ S2H + 4.40951E+07 5.81307E+07 5.63921E+07 8.42501E+07 9.90333E+07 + 1.07862E+08 1.15949E+08 9.38290E+07 1.14450E+08 1.10325E+08 +c setting distributions, x and E depend on x-position +c y and z have two sets but symmetric around center +ds11 s 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 + 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 + 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 + 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 +ds12 s 201 202 203 204 205 206 207 208 207 206 205 204 203 202 201 + 211 212 213 214 215 216 217 218 217 216 215 214 213 212 211 + 221 222 223 224 225 226 227 228 227 226 225 224 223 222 221 + 231 232 233 234 235 236 237 238 237 236 235 234 233 232 231 +ds13 s 301 302 303 304 305 306 307 308 307 306 305 304 303 302 301 + 311 312 313 314 315 316 317 318 317 316 315 314 313 312 311 + 321 322 323 324 325 326 327 328 327 326 325 324 323 322 321 + 331 332 333 334 335 336 337 338 337 336 335 334 333 332 331 +ds14 s 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 + 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 + 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 + 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 +c x-ranges for each segment +c Selene 1 +si101 1340 1388 +si102 1388 1436 +si103 1436 1484 +si104 1484 1532 +si105 1532 1580 +si106 1578 1628 +si107 1628 1676 +si108 1676 1724 +si109 1724 1772 +si110 1772 1820 +si111 1820 1868 +si112 1868 1916 +si113 1916 1964 +si114 1964 2012 +si115 2012 2060 +c Selene 2 +si121 2540 2588 +si122 2588 2636 +si123 2636 2684 +si124 2684 2732 +si125 2732 2780 +si126 2780 2828 +si127 2828 2876 +si128 2876 2924 +si129 2924 2972 +si130 2972 3020 +si131 3020 3068 +si132 3068 3116 +si133 3116 3164 +si134 3164 3212 +si135 3212 3260 +sp101 0 1 +sp102 0 1 +sp103 0 1 +sp104 0 1 +sp105 0 1 +sp106 0 1 +sp107 0 1 +sp108 0 1 +sp109 0 1 +sp110 0 1 +sp111 0 1 +sp112 0 1 +sp113 0 1 +sp114 0 1 +sp115 0 1 +sp121 0 1 +sp122 0 1 +sp123 0 1 +sp124 0 1 +sp125 0 1 +sp126 0 1 +sp127 0 1 +sp128 0 1 +sp129 0 1 +sp130 0 1 +sp131 0 1 +sp132 0 1 +sp133 0 1 +sp134 0 1 +sp135 0 1 +c y-ranges, first the ellipse distance for each segment, +c then beam coverage on horizontal mirror +si201 8.400 8.969 +si202 8.969 9.429 +si203 9.429 9.796 +si204 9.796 10.080 +si205 10.080 10.288 +si206 10.288 10.424 +si207 10.424 10.492 +si208 10.492 10.493 +si211 2.304 8.684 +si212 2.723 9.199 +si213 3.142 9.612 +si214 3.561 9.938 +si215 3.979 10.184 +si216 4.398 10.356 +si217 4.817 10.458 +si218 5.236 10.492 +si221 -8.969 -8.400 +si222 -9.429 -8.969 +si223 -9.796 -9.429 +si224 -10.080 -9.796 +si225 -10.288 -10.080 +si226 -10.424 -10.288 +si227 -10.492 -10.424 +si228 -10.493 -10.492 +si231 -8.684 -2.304 +si232 -9.199 -2.723 +si233 -9.612 -3.142 +si234 -9.938 -3.561 +si235 -10.184 -3.979 +si236 -10.356 -4.398 +si237 -10.458 -4.817 +si238 -10.492 -5.236 +sp201 0 1 +sp202 0 1 +sp203 0 1 +sp204 0 1 +sp205 0 1 +sp206 0 1 +sp207 0 1 +sp208 0 1 +sp211 0 1 +sp212 0 1 +sp213 0 1 +sp214 0 1 +sp215 0 1 +sp216 0 1 +sp217 0 1 +sp218 0 1 +sp221 0 1 +sp222 0 1 +sp223 0 1 +sp224 0 1 +sp225 0 1 +sp226 0 1 +sp227 0 1 +sp228 0 1 +sp231 0 1 +sp232 0 1 +sp233 0 1 +sp234 0 1 +sp235 0 1 +sp236 0 1 +sp237 0 1 +sp238 0 1 +c z-ranges, first beam coverage on horizontal mirror, +c then the ellipse distance for each segment, +si301 -8.684 -2.304 +si302 -9.199 -2.723 +si303 -9.612 -3.142 +si304 -9.938 -3.561 +si305 -10.184 -3.979 +si306 -10.356 -4.398 +si307 -10.458 -4.817 +si308 -10.492 -5.236 +si311 -8.969 -8.400 +si312 -9.429 -8.969 +si313 -9.796 -9.429 +si314 -10.080 -9.796 +si315 -10.288 -10.080 +si316 -10.424 -10.288 +si317 -10.492 -10.424 +si318 -10.493 -10.492 +si321 2.304 8.684 +si322 2.723 9.199 +si323 3.142 9.612 +si324 3.561 9.938 +si325 3.979 10.184 +si326 4.398 10.356 +si327 4.817 10.458 +si328 5.236 10.492 +si331 8.400 8.969 +si332 8.969 9.429 +si333 9.429 9.796 +si334 9.796 10.080 +si335 10.080 10.288 +si336 10.288 10.424 +si337 10.424 10.492 +si338 10.492 10.493 +sp301 0 1 +sp302 0 1 +sp303 0 1 +sp304 0 1 +sp305 0 1 +sp306 0 1 +sp307 0 1 +sp308 0 1 +sp311 0 1 +sp312 0 1 +sp313 0 1 +sp314 0 1 +sp315 0 1 +sp316 0 1 +sp317 0 1 +sp318 0 1 +sp321 0 1 +sp322 0 1 +sp323 0 1 +sp324 0 1 +sp325 0 1 +sp326 0 1 +sp327 0 1 +sp328 0 1 +sp331 0 1 +sp332 0 1 +sp333 0 1 +sp334 0 1 +sp335 0 1 +sp336 0 1 +sp337 0 1 +sp338 0 1 +c Energy distribution for all positions along the guide for Selene 1 +si401 1e-12 0.15 0.2 0.3 0.4 0.5 0.6 0.8 1.0 1.5 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 -sp101 0.00000E+00 6.74946E+03 1.26349E+04 4.02190E+04 1.59434E+05 +sp401 0.00000E+00 1.56557E+05 2.76591E+05 8.82267E+05 3.79148E+06 + 7.22100E+09 7.15782E+04 2.36360E+05 1.64987E+06 1.20477E+07 + 5.83592E+06 5.13734E+06 1.05131E+07 9.13771E+06 2.23406E+06 + 5.55688E+06 3.73177E+06 9.84371E+06 3.67345E+04 5.62062E+04 +C ------------------------------------------------------------------------ +si402 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp402 0.00000E+00 1.75066E+05 3.14271E+05 1.00187E+06 4.19715E+06 + 4.29576E+08 6.49313E+04 8.35318E+04 1.59407E+06 1.13391E+07 + 3.24825E+06 2.40414E+06 3.23346E+06 2.35587E+06 1.81032E+06 + 4.22018E+06 2.81568E+06 1.08595E+07 3.50283E+04 1.16665E+04 +C ------------------------------------------------------------------------ +si403 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp403 0.00000E+00 1.14282E+05 2.09752E+05 6.68141E+05 2.71818E+06 + 2.88507E+08 4.31570E+04 5.55428E+04 1.05765E+06 7.30837E+06 + 2.10495E+06 1.58164E+06 2.11537E+06 1.54343E+06 1.20473E+06 + 2.79310E+06 1.87709E+06 7.24741E+06 2.27791E+04 7.70741E+03 +C ------------------------------------------------------------------------ +si404 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp404 0.00000E+00 1.19942E+05 2.24594E+05 7.14916E+05 2.83168E+06 + 2.69546E+08 4.59551E+04 5.82997E+04 1.12511E+06 7.56863E+06 + 2.17669E+06 1.65382E+06 2.17880E+06 1.58300E+06 1.28282E+06 + 2.95750E+06 2.00035E+06 7.75812E+06 2.37907E+04 7.90258E+03 +C ------------------------------------------------------------------------ +si405 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp405 0.00000E+00 1.02808E+05 1.94757E+05 6.19690E+05 2.41655E+06 + 2.29099E+08 3.97544E+04 5.03268E+04 9.72514E+05 6.44006E+06 + 1.85566E+06 1.42096E+06 1.86323E+06 1.35351E+06 1.11024E+06 + 2.55211E+06 1.73259E+06 6.72701E+06 2.03449E+04 6.77894E+03 +C ------------------------------------------------------------------------ +si406 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp406 0.00000E+00 8.86042E+04 1.69642E+05 5.39580E+05 2.07421E+06 + 1.99597E+08 3.45601E+04 4.37451E+04 8.44765E+05 5.51349E+06 + 1.59281E+06 1.22897E+06 1.60647E+06 1.16764E+06 9.65712E+05 + 2.21413E+06 1.50825E+06 5.85933E+06 1.74992E+04 5.87135E+03 +C ------------------------------------------------------------------------ +si407 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp407 0.00000E+00 7.87877E+04 1.51840E+05 4.82849E+05 1.83971E+06 + 1.73871E+08 3.08863E+04 3.89901E+04 7.54666E+05 4.88100E+06 + 1.41072E+06 1.09310E+06 1.42351E+06 1.03396E+06 8.63179E+05 + 1.97565E+06 1.34859E+06 5.24414E+06 1.55375E+04 5.20473E+03 +C ------------------------------------------------------------------------ +si408 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp408 0.00000E+00 5.01840E+04 9.68964E+04 3.08110E+05 1.17096E+06 + 1.12987E+08 1.97075E+04 2.49206E+04 4.81428E+05 3.10577E+06 + 8.98765E+05 6.97579E+05 9.09029E+05 6.60776E+05 5.50891E+05 + 1.26041E+06 8.60883E+05 3.34661E+06 9.89468E+03 3.33164E+03 +C ------------------------------------------------------------------------ +si409 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp409 0.00000E+00 4.11446E+04 7.95509E+04 2.52944E+05 9.59525E+05 + 9.23981E+07 1.61749E+04 2.04454E+04 3.95097E+05 2.54402E+06 + 7.36328E+05 5.72028E+05 7.44921E+05 5.41445E+05 4.52162E+05 + 1.03416E+06 7.06655E+05 2.74751E+06 8.11001E+03 2.73082E+03 +C ------------------------------------------------------------------------ +si410 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp410 0.00000E+00 2.51204E+04 4.88245E+04 1.55217E+05 5.84610E+05 + 5.28251E+07 9.90982E+03 1.24431E+04 2.42028E+05 1.54694E+06 + 4.46975E+05 3.48156E+05 4.50560E+05 3.26731E+05 2.76963E+05 + 6.32447E+05 4.32872E+05 1.68610E+06 4.94348E+03 1.64531E+03 +C ------------------------------------------------------------------------ +si411 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp411 0.00000E+00 2.36608E+04 4.56443E+04 1.45144E+05 5.52264E+05 + 4.88854E+07 9.27597E+03 1.16381E+04 2.26683E+05 1.46399E+06 + 4.22049E+05 3.26897E+05 4.23743E+05 3.07060E+05 2.59131E+05 + 5.92778E+05 4.04758E+05 1.57629E+06 4.66257E+03 1.54109E+03 +C ------------------------------------------------------------------------ +si412 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp412 0.00000E+00 2.34986E+04 4.48276E+04 1.42601E+05 5.50860E+05 + 4.89730E+07 9.13069E+03 1.14783E+04 2.23305E+05 1.46457E+06 + 4.21412E+05 3.23890E+05 4.21794E+05 3.05694E+05 2.54959E+05 + 5.84877E+05 3.97945E+05 1.54818E+06 4.64119E+03 1.52926E+03 +C ------------------------------------------------------------------------ +si413 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp413 0.00000E+00 1.87764E+04 3.52661E+04 1.12245E+05 4.42775E+05 + 3.94292E+07 7.20582E+03 9.07992E+03 1.76424E+05 1.18187E+06 + 3.39132E+05 2.57888E+05 3.37896E+05 2.44903E+05 2.01079E+05 + 4.63087E+05 3.13513E+05 1.21806E+06 3.72001E+03 1.21944E+03 +C ------------------------------------------------------------------------ +si414 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp414 0.00000E+00 1.82897E+04 3.37346E+04 1.07439E+05 4.34219E+05 + 4.04983E+07 6.92221E+03 8.78390E+03 1.69672E+05 1.16466E+06 + 3.33775E+05 2.50952E+05 3.32061E+05 2.41077E+05 1.93078E+05 + 4.46799E+05 3.00726E+05 1.16534E+06 3.63777E+03 1.19682E+03 +C ------------------------------------------------------------------------ +si415 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp415 0.00000E+00 1.08301E+04 1.94946E+04 6.21412E+04 2.59394E+05 + 2.48374E+07 4.02189E+03 5.13569E+03 9.87454E+04 6.99903E+05 + 1.99988E+05 1.48066E+05 1.98045E+05 1.43923E+05 1.12085E+05 + 2.61025E+05 1.74300E+05 6.73544E+05 2.16453E+03 7.10515E+02 +C ------------------------------------------------------------------------ +c Energy distribution for all positions along the guide for Selene 2 +si421 1e-12 0.15 0.2 0.3 0.4 + 0.5 0.6 0.8 1.0 1.5 + 2.0 3.0 4.0 5.0 6.0 + 7.0 8.0 9.0 10.0 11.0 +sp421 0.00000E+00 6.74946E+03 1.26349E+04 4.02190E+04 1.59434E+05 4.37259E+07 2.64491E+03 3.95862E+03 6.43047E+04 4.33561E+05 1.34669E+05 1.05571E+05 1.54428E+05 1.18439E+05 7.48297E+04 1.73991E+05 1.17757E+05 4.37637E+05 1.36106E+03 6.36018E+02 C ------------------------------------------------------------------------ -si102 1e-12 0.15 0.2 0.3 0.4 +si422 1e-12 0.15 0.2 0.3 0.4 0.5 0.6 0.8 1.0 1.5 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 -sp102 0.00000E+00 1.17151E+04 2.23209E+04 7.10081E+04 2.74744E+05 +sp422 0.00000E+00 1.17151E+04 2.23209E+04 7.10081E+04 2.74744E+05 1.91981E+07 4.53666E+03 5.59356E+03 1.11043E+05 7.29343E+05 2.07983E+05 1.59135E+05 2.04496E+05 1.47055E+05 1.26489E+05 2.89983E+05 1.97214E+05 7.70669E+05 2.31033E+03 7.27223E+02 C ------------------------------------------------------------------------ -si103 1e-12 0.15 0.2 0.3 0.4 +si423 1e-12 0.15 0.2 0.3 0.4 0.5 0.6 0.8 1.0 1.5 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 -sp103 0.00000E+00 1.20232E+04 2.32095E+04 7.38021E+04 2.80546E+05 +sp423 0.00000E+00 1.20232E+04 2.32095E+04 7.38021E+04 2.80546E+05 2.04130E+07 4.70687E+03 5.81143E+03 1.15089E+05 7.42422E+05 2.12512E+05 1.64201E+05 2.10375E+05 1.51470E+05 1.31339E+05 3.00162E+05 2.04993E+05 8.01335E+05 2.36550E+03 7.53558E+02 C ------------------------------------------------------------------------ -si104 1e-12 0.15 0.2 0.3 0.4 +si424 1e-12 0.15 0.2 0.3 0.4 0.5 0.6 0.8 1.0 1.5 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 -sp104 0.00000E+00 1.75885E+04 3.43282E+04 1.09117E+05 4.08635E+05 +sp424 0.00000E+00 1.75885E+04 3.43282E+04 1.09117E+05 4.08635E+05 3.05406E+07 6.94852E+03 8.58466E+03 1.69754E+05 1.07843E+06 3.09623E+05 2.41213E+05 3.08149E+05 2.22054E+05 1.94011E+05 4.42225E+05 3.03071E+05 1.18520E+06 3.45334E+03 1.10997E+03 C ------------------------------------------------------------------------ -si105 1e-12 0.15 0.2 0.3 0.4 +si425 1e-12 0.15 0.2 0.3 0.4 0.5 0.6 0.8 1.0 1.5 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 -sp105 0.00000E+00 1.85638E+04 3.67221E+04 1.16674E+05 4.28978E+05 +sp425 0.00000E+00 1.85638E+04 3.67221E+04 1.16674E+05 4.28978E+05 3.21759E+07 7.41406E+03 9.14695E+03 1.80954E+05 1.12799E+06 3.24748E+05 2.55490E+05 3.24699E+05 2.34011E+05 2.07133E+05 4.70575E+05 3.23871E+05 1.26777E+06 3.63481E+03 1.17512E+03 C ------------------------------------------------------------------------ -si106 1e-12 0.15 0.2 0.3 0.4 +si426 1e-12 0.15 0.2 0.3 0.4 0.5 0.6 0.8 1.0 1.5 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 -sp106 0.00000E+00 2.21857E+04 4.43306E+04 1.40802E+05 5.10575E+05 +sp426 0.00000E+00 2.21857E+04 4.43306E+04 1.40802E+05 5.10575E+05 3.77016E+07 8.93175E+03 1.09930E+04 2.17852E+05 1.33861E+06 3.85960E+05 3.05846E+05 3.86785E+05 2.78627E+05 2.49620E+05 5.65666E+05 3.90547E+05 1.53035E+06 4.33437E+03 1.40296E+03 C ------------------------------------------------------------------------ -si107 1e-12 0.15 0.2 0.3 0.4 +si427 1e-12 0.15 0.2 0.3 0.4 0.5 0.6 0.8 1.0 1.5 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 -sp107 0.00000E+00 3.02306E+04 6.08550E+04 1.93239E+05 6.93589E+05 +sp427 0.00000E+00 3.02306E+04 6.08550E+04 1.93239E+05 6.93589E+05 4.94114E+07 1.22401E+04 1.50130E+04 2.98419E+05 1.81412E+06 5.23227E+05 4.16740E+05 5.24421E+05 3.77374E+05 3.42124E+05 7.73787E+05 5.35477E+05 2.10065E+06 5.89540E+03 1.90221E+03 C ------------------------------------------------------------------------ -si108 1e-12 0.15 0.2 0.3 0.4 +si428 1e-12 0.15 0.2 0.3 0.4 0.5 0.6 0.8 1.0 1.5 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 -sp108 0.00000E+00 2.94446E+04 5.92364E+04 1.88103E+05 6.75727E+05 +sp428 0.00000E+00 2.94446E+04 5.92364E+04 1.88103E+05 6.75727E+05 4.79016E+07 1.19154E+04 1.46108E+04 2.90519E+05 1.76765E+06 5.09678E+05 4.05737E+05 5.10570E+05 3.67353E+05 3.33032E+05 7.53323E+05 5.21216E+05 2.04477E+06 5.74269E+03 1.85095E+03 C ------------------------------------------------------------------------ -si109 1e-12 0.15 0.2 0.3 0.4 +si429 1e-12 0.15 0.2 0.3 0.4 0.5 0.6 0.8 1.0 1.5 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 -sp109 0.00000E+00 4.35594E+04 8.73170E+04 2.77304E+05 1.00114E+06 +sp429 0.00000E+00 4.35594E+04 8.73170E+04 2.77304E+05 1.00114E+06 7.16969E+07 1.75773E+04 2.15781E+04 4.28663E+05 2.62181E+06 7.55655E+05 6.00004E+05 7.56553E+05 5.44496E+05 4.91231E+05 1.11219E+06 7.68647E+05 3.01416E+06 8.50262E+03 2.74125E+03 C ------------------------------------------------------------------------ -si110 1e-12 0.15 0.2 0.3 0.4 +si430 1e-12 0.15 0.2 0.3 0.4 0.5 0.6 0.8 1.0 1.5 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 -sp110 0.00000E+00 5.08122E+04 1.01261E+05 3.21649E+05 1.17065E+06 +sp430 0.00000E+00 5.08122E+04 1.01261E+05 3.21649E+05 1.17065E+06 8.44071E+07 2.04082E+04 2.50830E+04 4.97895E+05 3.07096E+06 8.84254E+05 6.99122E+05 8.83955E+05 6.36310E+05 5.70220E+05 1.29292E+06 8.91913E+05 3.49561E+06 9.93105E+03 3.19801E+03 C ------------------------------------------------------------------------ -si111 1e-12 0.15 0.2 0.3 0.4 +si431 1e-12 0.15 0.2 0.3 0.4 0.5 0.6 0.8 1.0 1.5 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 -sp111 0.00000E+00 5.54824E+04 1.09734E+05 3.48652E+05 1.28218E+06 +sp431 0.00000E+00 5.54824E+04 1.09734E+05 3.48652E+05 1.28218E+06 9.19204E+07 2.21468E+04 2.72337E+04 5.40609E+05 3.37053E+06 9.68856E+05 7.61698E+05 9.65731E+05 6.95046E+05 6.18581E+05 1.40515E+06 9.67039E+05 3.78823E+06 1.08606E+04 3.48353E+03 C ------------------------------------------------------------------------ -si112 1e-12 0.15 0.2 0.3 0.4 +si432 1e-12 0.15 0.2 0.3 0.4 0.5 0.6 0.8 1.0 1.5 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 -sp112 0.00000E+00 5.97040E+04 1.16814E+05 3.71278E+05 1.38574E+06 +sp432 0.00000E+00 5.97040E+04 1.16814E+05 3.71278E+05 1.38574E+06 9.88383E+07 2.36237E+04 2.90770E+04 5.77107E+05 3.65346E+06 1.04777E+06 8.17221E+05 1.04035E+06 7.48620E+05 6.59504E+05 1.50209E+06 1.03023E+06 4.03280E+06 1.17128E+04 3.73774E+03 C ------------------------------------------------------------------------ -si113 1e-12 0.15 0.2 0.3 0.4 +si433 1e-12 0.15 0.2 0.3 0.4 0.5 0.6 0.8 1.0 1.5 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 -sp113 0.00000E+00 4.85803E+04 9.36078E+04 2.97674E+05 1.13437E+06 +sp433 0.00000E+00 4.85803E+04 9.36078E+04 2.97674E+05 1.13437E+06 7.99568E+07 1.89850E+04 2.33908E+04 4.64308E+05 3.00271E+06 8.58300E+05 6.62033E+05 8.47386E+05 6.09533E+05 5.29619E+05 1.21081E+06 8.26424E+05 3.23184E+06 9.55949E+03 3.02649E+03 C ------------------------------------------------------------------------ -si114 1e-12 0.15 0.2 0.3 0.4 +si434 1e-12 0.15 0.2 0.3 0.4 0.5 0.6 0.8 1.0 1.5 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 -sp114 0.00000E+00 5.88344E+04 1.11186E+05 3.53808E+05 1.38411E+06 +sp434 0.00000E+00 5.88344E+04 1.11186E+05 3.53808E+05 1.38411E+06 9.77202E+07 2.26369E+04 2.79595E+04 5.54391E+05 3.68221E+06 1.04882E+06 7.98041E+05 1.02934E+06 7.40419E+05 6.30958E+05 1.44953E+06 9.83220E+05 3.83908E+06 1.16223E+04 3.65352E+03 C ------------------------------------------------------------------------ -si115 1e-12 0.15 0.2 0.3 0.4 +si435 1e-12 0.15 0.2 0.3 0.4 0.5 0.6 0.8 1.0 1.5 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 -sp115 0.00000E+00 5.61276E+04 1.03236E+05 3.28823E+05 1.33382E+06 +sp435 0.00000E+00 5.61276E+04 1.03236E+05 3.28823E+05 1.33382E+06 9.44561E+07 2.11336E+04 2.61958E+04 5.18590E+05 3.57223E+06 1.01277E+06 7.56511E+05 9.86081E+05 7.09333E+05 5.88350E+05 1.36095E+06 9.15062E+05 3.56505E+06 1.11463E+04 3.47063E+03 C ------------------------------------------------------------------------ -c x-distribution -si31 2539.25 - 2587.35 2635.45 2683.55 2731.65 2779.75 - 2827.85 2875.95 2924.05 2972.15 3020.25 - 3068.35 3116.45 3164.55 3212.65 3260.75 -sp31 0 - 4.57688E+07 2.25345E+07 2.38501E+07 3.55841E+07 3.75155E+07 - 4.40951E+07 5.81307E+07 5.63921E+07 8.42501E+07 9.90333E+07 - 1.07862E+08 1.15949E+08 9.38290E+07 1.14450E+08 1.10325E+08 -c y-distribution -si32 -8 -2 -sp32 0 1 -c z-distribution -si33 2 8 -sp33 0 1 diff --git a/tallies/selene_source_mesh_p.i b/tallies/selene_source_mesh_p.i new file mode 100644 index 0000000..34fb0ba --- /dev/null +++ b/tallies/selene_source_mesh_p.i @@ -0,0 +1,38 @@ +c high resolution mesh around neutron mirrors to check gamma source +tmesh + rmesh1:p dose 1 trans 1 + cora1 1300 31i 2100 3i 2500 31i 3300 + corb1 -30 59i 30 + corc1 -30 59i 30 +c Official ESS dose conversion factors from gamma/cm²/s to Sv/h +c Is described in ESS-0019931 + mshmf1 0.000000 1.00000e-04 + 1.00E-02 2.47E-10 1.50E-02 5.62E-10 + 2.00E-02 8.1E-10 3.00E-02 1.13E-09 + 4.00E-02 1.26E-09 5.00E-02 1.33E-09 + 6.00E-02 1.4E-09 7.00E-02 1.49E-09 + 8.00E-02 1.6E-09 1.00E-01 1.87E-09 + 1.50E-01 2.69E-09 2.00E-01 3.6E-09 + 3.00E-01 5.44E-09 4.00E-01 7.2E-09 + 5.00E-01 8.89E-09 5.11E-01 9.07E-09 + 6.00E-01 1.05E-08 6.62E-01 1.14E-08 + 8.00E-01 1.34E-08 1.00E+00 1.62E-08 + 1.12E+00 1.76E-08 1.33E+00 2.01E-08 + 1.50E+00 2.2E-08 2.00E+00 2.69E-08 + 3.00E+00 3.51E-08 4.00E+00 4.21E-08 + 5.00E+00 4.82E-08 6.00E+00 5.4E-08 + 6.13E+00 5.47E-08 8.00E+00 6.7E-08 + 1.00E+01 7.92E-08 1.50E+01 1.09E-07 + 2.00E+01 1.38E-07 3.00E+01 1.85E-07 + 4.00E+01 2.23E-07 5.00E+01 2.6E-07 + 6.00E+01 2.95E-07 8.00E+01 3.52E-07 + 1.00E+02 3.96E-07 1.50E+02 4.68E-07 + 2.00E+02 5.15E-07 3.00E+02 5.8E-07 + 4.00E+02 6.19E-07 5.00E+02 6.48E-07 + 6.00E+02 6.7E-07 8.00E+02 7.02E-07 + 1.00E+03 7.24E-07 1.50E+03 7.63E-07 + 2.00E+03 7.92E-07 3.00E+03 8.35E-07 + 4.00E+03 8.75E-07 5.00E+03 9.04E-07 + 6.00E+03 9.29E-07 8.00E+03 9.65E-07 + 1.00E+04 9.94E-07 +endmd