Update gamma source for selene 1+2 and add fine mesh to test source

This commit is contained in:
2019-10-16 15:58:29 +02:00
parent d42f0d43cf
commit cedc4cee19
4 changed files with 456 additions and 62 deletions
+7 -6
View File
@@ -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
+1 -1
View File
@@ -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
+410 -55
View File
@@ -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
+38
View File
@@ -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