ESS-Estia/Estia-MCNP
1
0
Fork 0
Code Actions Packages Releases Activity

First version of VS source with energy deposition detectors for VS stages and chopper

Browse Source
This commit is contained in:
glavic_a
2018-10-03 21:24:40 +02:00
parent ffb3910bef
commit 59fead20ab
4 changed files with 584 additions and 79 deletions
Download Patch File Download Diff File Expand all files Collapse all files

128
e2-estia.i
View File

@@ -709,7 +709,7 @@ c -------------------------------------------------------
c --------------- WEIGHT CARDS --------------------------
c -------------------------------------------------------
c weight windows from external file generated by WWG
WWP:n,p j j 10 j -1 j 0.075 $ WW was generated for higher source scaling
c WWP:n,p j j 10 j -1 j 0.075 $ WW was generated for higher source scaling
c weight window generator, can be run with previous weight window files
c WWGE:n,p 1e-3 1e-1 1e1 3e3
c WWG:n,p 5 0 0
@@ -724,26 +724,8 @@ c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
c -----------------------------------------------------------
c ------------------- TALLY CARDS ---------------------------
c
fc14 neutron flux at BBG
f14:n 701
e14 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10
1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9
1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8
1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7
1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6
1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5
1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4
1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3
1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2
1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1
1e-0 1.58 2.51 3.98 6.31
1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1
1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2
1e+3 1.58e+3 2.51e+3
fq14 e f
fm14 8.7374e15
fc24 neutron flux before CPC1 and behind CPC1 in beam and next to beam
f24:n 702 704 703
fc24 neutron flux at VS
f24:n 704
e24 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10
1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9
1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8
@@ -760,61 +742,44 @@ e24 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10
1e+3 1.58e+3 2.51e+3
fq24 e f
fm24 8.7374e15
fc34 neutron flux behind CPC2 and BWC
f34:n 705 706
e34 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10
1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9
1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8
1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7
1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6
1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5
1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4
1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3
1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2
1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1
1e-0 1.58 2.51 3.98 6.31
1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1
1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2
1e+3 1.58e+3 2.51e+3
fq34 e f
fm34 8.7374e15
fc44 neutron flux outside bunker wall E01 and E03
f44:n 708 709
e44 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10
1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9
1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8
1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7
1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6
1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5
1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4
1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3
1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2
1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1
1e-0 1.58 2.51 3.98 6.31
1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1
1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2
1e+3 1.58e+3 2.51e+3
fq44 e f
fm44 8.7374e15
fc54 neutron flux behind instrument shutter
f54:n 707
e54 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10
1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9
1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8
1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7
1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6
1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5
1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4
1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3
1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2
1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1
1e-0 1.58 2.51 3.98 6.31
1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1
1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2
1e+3 1.58e+3 2.51e+3
fq54 e f
fm54 8.7374e15
c
fc106 deposited energy at chopper disk and spindle
+f106 4121 4131
e106 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10
1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9
1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8
1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7
1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6
1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5
1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4
1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3
1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2
1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1
1e-0 1.58 2.51 3.98 6.31
1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1
1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2
1e+3 1.58e+3 2.51e+3
fq106 e f
fm106 8.7374e15
c
fc6 deposited energy at VS stages x/z/rot
+f6 (5711 5712) (5721 5722) (5731)
e6 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10
1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9
1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8
1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7
1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6
1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5
1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4
1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3
1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2
1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1
1e-0 1.58 2.51 3.98 6.31
1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1
1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2
1e+3 1.58e+3 2.51e+3
fq6 e f
fm6 8.7374e15
c ==========================================
c
c ==========================================
@@ -844,11 +809,12 @@ c -------------------------------------------------------
c --------------- PHYSICS CARDS --------------------------
c -------------------------------------------------------
c dbcn 375642321 j j j j j j j j j j j 15291711
nps 1e9
nps 1e8
c histp
mode n p n h / d t s a
c calculate volume for tally cells (spheres)
vol 1 507r 10j 1 171r
c volume calculated for tally cells (spheres)
c and set for chopper disk/spindle and stages x1/x2/z1/z2/rot
vol 1 507r 10j 1 107r 3848.5 10104 1 13r 479.5 479.5 913 913 1734 1 42r
imp:h,/,d,t,s,a 0 0 1 687r
imp:n,p 0 0 1 687r
cut:h,/,d,t,s,a 1e+08 1e-3
@@ -862,8 +828,12 @@ phys:/,d,t,s,a 3000
phys:h 3000
lca 2 1 1 0023 1 1 0 1 1 0
lea 1 4 1 0 1 0 0 1
c * Source for simulation from proton beam (needs different WW and ESplit)
c read file=source_proton.i echo
read file=e2-source_2m.i echo
c * Source Term from Valentina for 2m adapted to size
c read file=e2-source_2m.i echo
c * McStas source terms for gamma doeses (needs no WW)
read file=mcstas_source_vs.i echo
prdmp 5e7 -60 1 2 5e6
print 30 32 $ 120 128 150 160 161 162
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++

363
mcstas_reference/VS_L_total.dat Normal file
View File

@@ -0,0 +1,363 @@
# Format: McCode with text headers
# URL: http://www.mccode.org
# Creator: McStas 2.4.1 - Jun. 26, 2017
# Instrument: Estia_vs.instr
# Ncount: 16666666
# Trace: no
# Gravitation: yes
# Seed: 1516871856
# Directory: .//vs_data
# Nodes: 6
# Param: omegaa=2
# Param: sample=1
# Param: sample_length=0.01
# Param: sample_height=0.015
# Param: operationmode=0
# Param: over_illumination=0.0025
# Param: theta_resolution=0.04
# Param: chopper_diameter=0.7
# Param: enable_gravity=1
# Param: enable_chopper=0
# Param: source_power=5
# Date: Thu Jan 25 09:17:20 2018 (1516868240)
# type: array_1d(325)
# Source: ESS_reflectometer_Estia (Estia_vs.instr)
# component: L_total
# position: 0.0967598 -0.134387 11.0536
# title: Wavelength monitor
# Ncount: 99999996
# filename: L_total.dat
# statistics: X0=4.88633; dX=2.39665;
# signal: Min=0; Max=1.58358e+09; Mean=1.31337e+08;
# values: 4.26845e+10 9.31529e+07 4.15394e+06
# xvar: L
# yvar: (I,I_err)
# xlabel: Wavelength [AA]
# ylabel: Intensity
# xlimits: 0.5 32
# variables: L I I_err N
0.5484615385 0 0 0
0.6453846154 0 0 0
0.7423076923 0 0 0
0.8392307692 0 0 0
0.9361538462 0 0 0
1.033076923 0 0 0
1.13 0 0 0
1.226923077 0 0 0
1.323846154 0 0 0
1.420769231 0 0 0
1.517692308 1.804512377e-05 2.017505087e-05 4
1.614615385 9.546958703e-05 8.51397802e-05 82
1.711538462 5.013169782 3.629310286 456
1.808461538 234464.456 234324.6096 1539
1.905384615 354331.3535 353390.386 3302
2.002307692 1223193.492 1110226.631 5467
2.099230769 1484227.412 302940.4747 7995
2.196153846 13171446.41 1723612.817 10562
2.293076923 46726118.4 3454748.866 12386
2.39 133957359.1 5722774.014 13168
2.486923077 315212259.4 10005157.26 13587
2.583846154 562677491.6 14275612.98 13593
2.680769231 800377598.8 17098275.54 13563
2.777692308 1035081129 19341168.31 13434
2.874615385 1259570984 21053730.4 13458
2.971538462 1409912146 22034532.9 13423
3.068461538 1495261941 22565956.83 13335
3.165384615 1573230291 22936567.89 13454
3.262307692 1583582562 22738005.34 13522
3.359230769 1509221671 21827106.04 13619
3.456153846 1475495546 21160709.81 13409
3.553076923 1382382491 20002563.19 13382
3.65 1331112277 19260473.4 13227
3.746923077 1276145569 18362368.71 13371
3.843846154 1219868995 17533764.13 13580
3.940769231 1147185641 16614885.41 13432
4.037692308 1120247827 16043864.82 13521
4.134615385 1034122132 14999725.05 13357
4.231538462 1001708082 14385972.41 13667
4.328461538 923554146.2 13394902.32 13424
4.425384615 893874143.7 12862809.87 13468
4.522307692 845849809.3 12204325.25 13439
4.619230769 804141962.1 11599866.12 13521
4.716153846 766057906.9 11012375.33 13637
4.813076923 712500227.9 10332163.85 13374
4.91 697131874.2 9993009.591 13596
5.006923077 655879369.8 9435641.498 13570
5.103846154 614556430.2 8882949.471 13447
5.200769231 584079194 8438571.767 13335
5.297692308 558368549.9 8035180.438 13367
5.394615385 532304244.4 7621451.539 13508
5.491538462 497197305.9 7167403.37 13693
5.588461538 472443879.6 6826597.21 13407
5.685384615 457712143.1 6539333.465 13602
5.782307692 425617993.6 6143139.75 13387
5.879230769 407089719.7 5869802.72 13523
5.976153846 379135775.6 5486205.957 13383
6.073076923 369476759.5 5296713.381 13449
6.17 347812875.4 4999787.809 13423
6.266923077 332011750.3 4757612.204 13488
6.363846154 316539008.4 4540009.846 13459
6.460769231 300497733.7 4315623.312 13559
6.557692308 284523951.9 4099798.897 13371
6.654615385 272899417.9 3911436.102 13490
6.751538462 261028003.2 3745580.661 13406
6.848461538 251518437.8 3583356.368 13726
6.945384615 236805438.1 3382684.865 13485
7.042307692 221159145.1 3186442.796 13412
7.139230769 214848354.8 3077509.463 13438
7.236153846 208514953.7 2958998.267 13537
7.333076923 197310772.1 2819397.109 13332
7.43 184724393.9 2653300.833 13313
7.526923077 176455888.4 2530669.755 13405
7.623846154 173926179.3 2465778.674 13715
7.720769231 163103438.7 2333325.338 13585
7.817692308 157040527.7 2234565.888 13584
7.914615385 146667489.4 2116783.877 13465
8.011538462 142259456.4 2039013.235 13421
8.108461538 137731731 1961563.54 13568
8.205384615 133504509.8 1892657.872 13609
8.302307692 127706055.4 1810973.615 13492
8.399230769 120802382.3 1729508.626 13413
8.496153846 116908561.5 1657161.346 13708
8.593076923 114338370.8 1609299.882 13596
8.69 108419265.2 1537199.827 13394
8.786923077 103211271.5 1462866.763 13509
8.883846154 101409323.7 1419545.058 13730
8.980769231 95750883.29 1353259.142 13581
9.077692308 90316761.8 1287512.019 13417
9.174615385 87246182 1239118.38 13365
9.271538462 84203001.53 1195451.766 13598
9.368461538 79251206.86 1139511.076 13434
9.465384615 78594592.56 1111453.56 13589
9.562307692 76495512.17 1075665.077 13562
9.659230769 73206986.16 1031624.696 13555
9.756153846 68731771.99 977101.1135 13459
9.853076923 66492574.31 946598.6289 13485
9.95 64154145.52 908354.8248 13544
10.04692308 63392045.89 891798.4722 13427
10.14384615 60317459.7 851494.5704 13585
10.24076923 58082928.84 820154.4409 13499
10.33769231 55413521.27 783249.9319 13564
10.43461538 54643187.74 764739.0768 13552
10.53153846 52773060.88 738571.5165 13471
10.62846154 49510113.68 701284.0428 13412
10.72538462 48486031.33 681628.6742 13556
10.82230769 45554800.63 646532.6445 13497
10.91923077 45830763.32 638044.7468 13569
11.01615385 43273784.38 608874.8894 13632
11.11307692 42048983.12 590388.6354 13487
11.21 40692300.74 569089.9319 13573
11.30692308 38655927.62 544667.9187 13518
11.40384615 37589965.88 527298.7197 13501
11.50076923 36750034.98 512095.7428 13707
11.59769231 35017707.07 490536.0385 13565
11.69461538 33511260.67 471667.6426 13350
11.79153846 32677826.49 457623.4359 13519
11.88846154 32024303.47 444981.3022 13445
11.98538462 30800464.19 429259.3647 13615
12.08230769 29514690.98 412980.511 13397
12.17923077 28591971.21 398182.1534 13471
12.27615385 27500647.62 384594.9943 13386
12.37307692 26916838.07 373857.8956 13603
12.47 25893346.53 359028.9158 13556
12.56692308 25425857.32 350603.6677 13439
12.66384615 23906836.3 332767.1321 13436
12.76076923 22926163.15 321878.708 13374
12.85769231 22060736.78 309707.4093 13421
12.95461538 21914899.08 303331.1483 13666
13.05153846 20501040.16 287051.8798 13487
13.14846154 19952767.71 278957.1389 13606
13.24538462 19744884.71 272600.9484 13530
13.34230769 19125160.51 264389.1844 13568
13.43923077 18057963.12 252369.417 13551
13.53615385 17675751.94 245141.4043 13319
13.63307692 17313987.5 238133.0457 13437
13.73 16347924.74 228088.5084 13383
13.82692308 15809819.73 220479.1942 13409
13.92384615 15468214.86 214242.3431 13449
14.02076923 15088768.34 207733.7091 13464
14.11769231 14604781.01 200940.1822 13299
14.21461538 14432425.04 196921.054 13729
14.31153846 13556495.44 187159.9898 13454
14.40846154 13328598.12 182657.2332 13486
14.50538462 12416671.88 172958.6497 13201
14.60230769 12518858.7 170232.6943 13632
14.69923077 12078971.44 164760.63 13536
14.79615385 11514590.62 157890.2981 13485
14.89307692 10962977.36 151230.8524 13466
14.99 10642630.08 147166.7606 13328
15.08692308 10341635.11 142209.3324 13476
15.18384615 10198039.5 138952.4148 13440
15.28076923 9658780.342 132682.6168 13367
15.37769231 9530428.308 130125.6605 13431
15.47461538 9207514.57 125723.4519 13429
15.57153846 8727791.436 119708.7048 13456
15.66846154 8596093.358 116913.8047 13668
15.76538462 8341021.616 113905.4179 13326
15.86230769 7981486.025 109133.5651 13399
15.95923077 7658072.028 105046.7035 13348
16.05615385 7507606.68 102424.6973 13482
16.15307692 7168083.236 98084.90439 13339
16.25 7242919.067 97024.93607 13611
16.34692308 6901561.553 93436.69443 13539
16.44384615 6569550.847 89521.97164 13368
16.54076923 6315867.76 86208.22502 13322
16.63769231 6221300.992 84145.5868 13338
16.73461538 6110233.169 82061.848 13547
16.83153846 5811312.475 78651.83064 13406
16.92846154 5660316.421 76331.34518 13504
17.02538462 5367257.755 72657.25011 13397
17.12230769 5201057.357 70773.14065 13442
17.21923077 5209116.043 69723.80539 13535
17.31615385 4969869.589 66761.78422 13559
17.41307692 4753536.465 64385.47605 13368
17.51 4618794.253 62135.77284 13516
17.60692308 4529987.481 60385.16488 13510
17.70384615 4271586.194 57692.45643 13329
17.80076923 4140433.773 55983.81375 13347
17.89769231 4129327.282 55118.80383 13339
17.99461538 3917627.186 52594.45637 13484
18.09153846 3729351.496 50376.64782 13389
18.18846154 3699360.32 49375.63681 13576
18.28538462 3591068.312 47951.18975 13406
18.38230769 3432132.574 46088.35001 13402
18.47923077 3426350.963 45270.42397 13555
18.57615385 3255159.222 43330.64086 13434
18.67307692 3115990.534 41527.65854 13481
18.77 3012187.322 40290.09584 13311
18.86692308 3054811.72 39929.55336 13640
18.96384615 2919317.099 38393.54192 13622
19.06076923 2790313.088 36828.24656 13502
19.15769231 2683526.613 35584.8971 13223
19.25461538 2625177.847 34616.33038 13460
19.35153846 2562774.756 33730.96816 13502
19.44846154 2449751.425 32206.95976 13451
19.54538462 2392355.287 31315.4148 13557
19.64230769 2277481.283 30064.82252 13447
19.73923077 2223152.592 29216.91869 13488
19.83615385 2079841.101 27723.64558 13155
19.93307692 2085460.513 27308.57127 13321
20.03 1970190.08 26048.08118 13170
20.12692308 1980347.413 25742.57316 13608
20.22384615 1882627.7 24667.69534 13354
20.32076923 1817009.141 23777.22562 13522
20.41769231 1807441.12 23360.47865 13658
20.51461538 1706702.921 22296.45736 13199
20.61153846 1649154.343 21582.46472 13215
20.70846154 1667076.695 21342.28833 13612
20.80538462 1562164.675 20297.54083 13631
20.90230769 1503963.599 19582.66911 13399
20.99923077 1448677.954 18804.22763 13474
21.09615385 1418235.915 18383.2995 13273
21.19307692 1395873.016 17985.35596 13384
21.29 1333594.425 17176.46053 13536
21.38692308 1267142.947 16451.12761 13410
21.48384615 1265824.135 16258.041 13485
21.58076923 1207974.954 15569.3802 13478
21.67769231 1164105.159 14968.10352 13470
21.77461538 1150336.803 14694.78973 13582
21.87153846 1110245.877 14179.86852 13496
21.96846154 1065093.724 13597.56349 13591
22.06538462 1033722.225 13273.49859 13323
22.16230769 972894.7138 12612.55039 13317
22.25923077 975556.4121 12402.60992 13589
22.35615385 925296.8367 11884.34829 13366
22.45307692 905485.3545 11557.42969 13384
22.55 885135.8873 11239.36717 13474
22.64692308 836832.1725 10731.89106 13310
22.74384615 820065.8859 10441.01435 13347
22.84076923 802850.5156 10175.72807 13457
22.93769231 773362.4282 9802.088391 13413
23.03461538 742132.1636 9446.173369 13418
23.13153846 718454.7572 9108.492921 13479
23.22846154 706050.235 8909.601316 13485
23.32538462 679485.8762 8561.140972 13488
23.42230769 652907.8325 8269.925772 13342
23.51923077 625421.0913 7927.097707 13284
23.61615385 616881.4791 7771.453672 13441
23.71307692 597108.4489 7508.073062 13532
23.81 580739.5517 7282.420796 13462
23.90692308 573420.4545 7125.835139 13554
24.00384615 534828.6059 6766.990189 13287
24.10076923 525901.5747 6590.446442 13308
24.19769231 519036.7816 6433.556877 13488
24.29461538 490622.926 6132.888965 13496
24.39153846 487145.0244 6052.357343 13369
24.48846154 461260.7931 5754.260081 13319
24.58538462 450918.8304 5584.427351 13534
24.68230769 430875.9691 5375.769424 13365
24.77923077 416443.2015 5184.989158 13331
24.87615385 405325.09 5013.108405 13515
24.97307692 391759.2889 4869.45276 13415
25.07 378413.9978 4676.736413 13517
25.16692308 369923.9737 4573.92567 13438
25.26384615 355845.703 4399.55204 13346
25.36076923 343742.7395 4250.889918 13376
25.45769231 339096.2494 4151.145219 13536
25.55461538 323742.7291 3991.48857 13235
25.65153846 311333.5074 3845.712867 13504
25.74846154 302109.189 3716.452675 13460
25.84538462 296900.1209 3641.456955 13577
25.94230769 287158.2403 3509.968117 13531
26.03923077 274175.5989 3364.586429 13344
26.13615385 268861.9965 3282.206403 13352
26.23307692 259135.8612 3161.671723 13387
26.33 251491.5497 3059.31828 13434
26.42692308 240858.7348 2944.204132 13249
26.52384615 231263.9592 2834.642339 13188
26.62076923 225185.0987 2746.764317 13282
26.71769231 225445.1254 2711.216558 13592
26.81461538 210910.0721 2569.097634 13363
26.91153846 203880.813 2477.175834 13483
27.00846154 199737.5831 2423.686653 13275
27.10538462 193350.6478 2343.215068 13336
27.20230769 187818.0766 2270.569118 13349
27.29923077 180296.4979 2187.36211 13282
27.39615385 178551.6445 2141.613564 13574
27.49307692 169688.0193 2058.202983 13178
27.59 163364.8238 1973.152364 13294
27.68692308 159877.9009 1919.247436 13386
27.78384615 156204.9412 1863.80961 13525
27.88076923 150040.0263 1790.270891 13501
27.97769231 143458.8879 1727.854504 13211
28.07461538 139613.3969 1671.10395 13442
28.17153846 137115.3989 1635.453271 13413
28.26846154 129672.5058 1556.369133 13103
28.36538462 122456.7296 1486.213632 13088
28.46230769 123025.6294 1468.560716 13342
28.55923077 117978.0165 1409.412827 13263
28.65615385 116472.601 1381.403111 13382
28.75307692 111112.9787 1322.125644 13314
28.85 110374.8368 1295.052148 13511
28.94692308 105545.4719 1249.271239 13262
29.04384615 102928.5056 1205.931833 13463
29.14076923 98959.6816 1166.079371 13391
29.23769231 93987.60836 1114.585717 13247
29.33461538 92653.30819 1089.110693 13411
29.43153846 90332.55801 1059.403694 13490
29.52846154 86648.38788 1013.91626 13487
29.62538462 84091.65297 990.0600405 13209
29.72230769 81958.60285 959.3004982 13315
29.81923077 78728.06541 926.1938328 13255
29.91615385 75932.37514 886.7225739 13399
30.01307692 74483.56798 868.5464541 13413
30.11 72052.104 838.0950472 13376
30.20692308 68441.16793 799.0188685 13290
30.30384615 67552.78454 783.4258878 13505
30.40076923 65049.37826 755.1546101 13325
30.49769231 63414.61361 733.9388814 13390
30.59461538 60985.0808 707.3681865 13394
30.69153846 58368.97466 678.5894295 13263
30.78846154 58132.49906 667.9489551 13438
30.88538462 54865.40704 637.3996461 13328
30.98230769 53132.93129 615.1598875 13232
31.07923077 52931.75978 606.5759982 13543
31.17615385 48957.05887 568.9599385 13130
31.27307692 47735.88313 554.3077828 13210
31.37 47385.90987 542.6063204 13336
31.46692308 45773.82776 524.9407021 13271
31.56384615 44882.07817 511.2777484 13476
31.66076923 41891.15942 484.0135041 12997
31.75769231 41681.72512 475.2583696 13420
31.85461538 39295.6672 452.1571257 13155
31.95153846 37199.78632 432.5416776 12719

27
mcstas_reference/convert_source.py Normal file
View File

@@ -0,0 +1,27 @@
#!/usr/bin/env python
from numpy import *
import sys
def convert(fname):
data=loadtxt(fname)[::-1].T
E=0.01239842/data[0] # convert wavelength into MeV energy
I=hstack([0., 0., data[1]])
Ebins=hstack([0, E[0]-(E[1]-E[0])/2., E[1:]-(E[1:]-E[:-1])/2., E[-1]+(E[-1]-E[-2])/2.])
out='si11 '
for i,Ei in enumerate(Ebins):
out+='%.5e'%Ei
if i%5==4:
out+='\n '
else:
out+=' '
out=out.strip()+'\nsp11 '
for i,Ii in enumerate(I):
out+='%.5e'%Ii
if i%5==4:
out+='\n '
else:
out+=' '
return out

145
mcstas_source_vs.i Normal file
View File

@@ -0,0 +1,145 @@
c cold neutron source from McStas simulation for VS position
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=5
par=n ara=165.0 wgt=4.83e-7
c scaling wgt is to proton source equivalent (McStas nps*1.13e-17)
si1 -1 0.99976203 1
sp1 0 0 1
si11 0.00000e+00 3.87448e-04 3.88629e-04 3.89813e-04 3.91004e-04
3.92203e-04 3.93409e-04 3.94623e-04 3.95844e-04 3.97073e-04
3.98309e-04 3.99553e-04 4.00805e-04 4.02065e-04 4.03333e-04
4.04609e-04 4.05892e-04 4.07184e-04 4.08485e-04 4.09793e-04
4.11110e-04 4.12436e-04 4.13770e-04 4.15113e-04 4.16464e-04
4.17824e-04 4.19193e-04 4.20572e-04 4.21959e-04 4.23355e-04
4.24761e-04 4.26176e-04 4.27601e-04 4.29035e-04 4.30479e-04
4.31932e-04 4.33396e-04 4.34869e-04 4.36353e-04 4.37846e-04
4.39350e-04 4.40864e-04 4.42389e-04 4.43924e-04 4.45470e-04
4.47027e-04 4.48594e-04 4.50173e-04 4.51763e-04 4.53364e-04
4.54976e-04 4.56600e-04 4.58236e-04 4.59884e-04 4.61543e-04
4.63214e-04 4.64898e-04 4.66593e-04 4.68301e-04 4.70022e-04
4.71756e-04 4.73502e-04 4.75261e-04 4.77033e-04 4.78819e-04
4.80618e-04 4.82431e-04 4.84257e-04 4.86097e-04 4.87951e-04
4.89820e-04 4.91702e-04 4.93600e-04 4.95512e-04 4.97439e-04
4.99381e-04 5.01338e-04 5.03310e-04 5.05298e-04 5.07302e-04
5.09322e-04 5.11358e-04 5.13411e-04 5.15479e-04 5.17565e-04
5.19668e-04 5.21787e-04 5.23925e-04 5.26079e-04 5.28252e-04
5.30442e-04 5.32651e-04 5.34878e-04 5.37124e-04 5.39389e-04
5.41673e-04 5.43976e-04 5.46299e-04 5.48643e-04 5.51006e-04
5.53389e-04 5.55794e-04 5.58219e-04 5.60666e-04 5.63134e-04
5.65624e-04 5.68136e-04 5.70671e-04 5.73228e-04 5.75808e-04
5.78412e-04 5.81039e-04 5.83691e-04 5.86366e-04 5.89066e-04
5.91791e-04 5.94542e-04 5.97318e-04 6.00120e-04 6.02949e-04
6.05805e-04 6.08687e-04 6.11597e-04 6.14536e-04 6.17502e-04
6.20497e-04 6.23522e-04 6.26576e-04 6.29660e-04 6.32775e-04
6.35921e-04 6.39098e-04 6.42307e-04 6.45548e-04 6.48823e-04
6.52130e-04 6.55472e-04 6.58848e-04 6.62259e-04 6.65705e-04
6.69188e-04 6.72707e-04 6.76263e-04 6.79857e-04 6.83490e-04
6.87162e-04 6.90873e-04 6.94624e-04 6.98417e-04 7.02251e-04
7.06128e-04 7.10047e-04 7.14010e-04 7.18018e-04 7.22071e-04
7.26170e-04 7.30316e-04 7.34510e-04 7.38752e-04 7.43043e-04
7.47384e-04 7.51776e-04 7.56221e-04 7.60718e-04 7.65269e-04
7.69874e-04 7.74536e-04 7.79254e-04 7.84030e-04 7.88865e-04
7.93761e-04 7.98717e-04 8.03735e-04 8.08817e-04 8.13964e-04
8.19176e-04 8.24456e-04 8.29804e-04 8.35222e-04 8.40711e-04
8.46273e-04 8.51909e-04 8.57621e-04 8.63410e-04 8.69277e-04
8.75225e-04 8.81254e-04 8.87367e-04 8.93566e-04 8.99852e-04
9.06227e-04 9.12693e-04 9.19251e-04 9.25905e-04 9.32656e-04
9.39506e-04 9.46457e-04 9.53512e-04 9.60673e-04 9.67942e-04
9.75323e-04 9.82816e-04 9.90426e-04 9.98154e-04 1.00600e-03
1.01398e-03 1.02208e-03 1.03031e-03 1.03868e-03 1.04718e-03
1.05582e-03 1.06461e-03 1.07355e-03 1.08263e-03 1.09187e-03
1.10127e-03 1.11084e-03 1.12057e-03 1.13047e-03 1.14055e-03
1.15081e-03 1.16126e-03 1.17190e-03 1.18273e-03 1.19377e-03
1.20502e-03 1.21648e-03 1.22816e-03 1.24006e-03 1.25220e-03
1.26458e-03 1.27721e-03 1.29009e-03 1.30323e-03 1.31665e-03
1.33034e-03 1.34432e-03 1.35860e-03 1.37318e-03 1.38808e-03
1.40331e-03 1.41888e-03 1.43479e-03 1.45107e-03 1.46772e-03
1.48475e-03 1.50219e-03 1.52004e-03 1.53832e-03 1.55705e-03
1.57623e-03 1.59590e-03 1.61606e-03 1.63674e-03 1.65795e-03
1.67972e-03 1.70208e-03 1.72503e-03 1.74861e-03 1.77285e-03
1.79776e-03 1.82339e-03 1.84976e-03 1.87690e-03 1.90485e-03
1.93365e-03 1.96332e-03 1.99393e-03 2.02550e-03 2.05809e-03
2.09175e-03 2.12652e-03 2.16248e-03 2.19966e-03 2.23815e-03
2.27801e-03 2.31932e-03 2.36215e-03 2.40659e-03 2.45274e-03
2.50070e-03 2.55056e-03 2.60246e-03 2.65651e-03 2.71285e-03
2.77164e-03 2.83303e-03 2.89720e-03 2.96435e-03 3.03468e-03
3.10843e-03 3.18586e-03 3.26724e-03 3.35289e-03 3.44316e-03
3.53842e-03 3.63910e-03 3.74568e-03 3.85869e-03 3.97874e-03
4.10649e-03 4.24273e-03 4.38832e-03 4.54426e-03 4.71169e-03
4.89194e-03 5.08653e-03 5.29726e-03 5.52620e-03 5.77584e-03
6.04912e-03 6.34955e-03 6.68141e-03 7.04990e-03 7.46144e-03
7.92406e-03 8.44791e-03 9.04600e-03 9.73537e-03 1.05387e-02
1.14868e-02 1.26227e-02 1.40088e-02 1.57380e-02 1.79567e-02
2.09084e-02 2.43033e-02
sp11 0.00000e+00 0.00000e+00 3.71998e+04 3.92957e+04 4.16817e+04
4.18912e+04 4.48821e+04 4.57738e+04 4.73859e+04 4.77359e+04
4.89571e+04 5.29318e+04 5.31329e+04 5.48654e+04 5.81325e+04
5.83690e+04 6.09851e+04 6.34146e+04 6.50494e+04 6.75528e+04
6.84412e+04 7.20521e+04 7.44836e+04 7.59324e+04 7.87281e+04
8.19586e+04 8.40917e+04 8.66484e+04 9.03326e+04 9.26533e+04
9.39876e+04 9.89597e+04 1.02929e+05 1.05545e+05 1.10375e+05
1.11113e+05 1.16473e+05 1.17978e+05 1.23026e+05 1.22457e+05
1.29673e+05 1.37115e+05 1.39613e+05 1.43459e+05 1.50040e+05
1.56205e+05 1.59878e+05 1.63365e+05 1.69688e+05 1.78552e+05
1.80296e+05 1.87818e+05 1.93351e+05 1.99738e+05 2.03881e+05
2.10910e+05 2.25445e+05 2.25185e+05 2.31264e+05 2.40859e+05
2.51492e+05 2.59136e+05 2.68862e+05 2.74176e+05 2.87158e+05
2.96900e+05 3.02109e+05 3.11334e+05 3.23743e+05 3.39096e+05
3.43743e+05 3.55846e+05 3.69924e+05 3.78414e+05 3.91759e+05
4.05325e+05 4.16443e+05 4.30876e+05 4.50919e+05 4.61261e+05
4.87145e+05 4.90623e+05 5.19037e+05 5.25902e+05 5.34829e+05
5.73420e+05 5.80740e+05 5.97108e+05 6.16881e+05 6.25421e+05
6.52908e+05 6.79486e+05 7.06050e+05 7.18455e+05 7.42132e+05
7.73362e+05 8.02851e+05 8.20066e+05 8.36832e+05 8.85136e+05
9.05485e+05 9.25297e+05 9.75556e+05 9.72895e+05 1.03372e+06
1.06509e+06 1.11025e+06 1.15034e+06 1.16411e+06 1.20797e+06
1.26582e+06 1.26714e+06 1.33359e+06 1.39587e+06 1.41824e+06
1.44868e+06 1.50396e+06 1.56216e+06 1.66708e+06 1.64915e+06
1.70670e+06 1.80744e+06 1.81701e+06 1.88263e+06 1.98035e+06
1.97019e+06 2.08546e+06 2.07984e+06 2.22315e+06 2.27748e+06
2.39236e+06 2.44975e+06 2.56277e+06 2.62518e+06 2.68353e+06
2.79031e+06 2.91932e+06 3.05481e+06 3.01219e+06 3.11599e+06
3.25516e+06 3.42635e+06 3.43213e+06 3.59107e+06 3.69936e+06
3.72935e+06 3.91763e+06 4.12933e+06 4.14043e+06 4.27159e+06
4.52999e+06 4.61879e+06 4.75354e+06 4.96987e+06 5.20912e+06
5.20106e+06 5.36726e+06 5.66032e+06 5.81131e+06 6.11023e+06
6.22130e+06 6.31587e+06 6.56955e+06 6.90156e+06 7.24292e+06
7.16808e+06 7.50761e+06 7.65807e+06 7.98149e+06 8.34102e+06
8.59609e+06 8.72779e+06 9.20751e+06 9.53043e+06 9.65878e+06
1.01980e+07 1.03416e+07 1.06426e+07 1.09630e+07 1.15146e+07
1.20790e+07 1.25189e+07 1.24167e+07 1.33286e+07 1.35565e+07
1.44324e+07 1.46048e+07 1.50888e+07 1.54682e+07 1.58098e+07
1.63479e+07 1.73140e+07 1.76758e+07 1.80580e+07 1.91252e+07
1.97449e+07 1.99528e+07 2.05010e+07 2.19149e+07 2.20607e+07
2.29262e+07 2.39068e+07 2.54259e+07 2.58933e+07 2.69168e+07
2.75006e+07 2.85920e+07 2.95147e+07 3.08005e+07 3.20243e+07
3.26778e+07 3.35113e+07 3.50177e+07 3.67500e+07 3.75900e+07
3.86559e+07 4.06923e+07 4.20490e+07 4.32738e+07 4.58308e+07
4.55548e+07 4.84860e+07 4.95101e+07 5.27731e+07 5.46432e+07
5.54135e+07 5.80829e+07 6.03175e+07 6.33920e+07 6.41541e+07
6.64926e+07 6.87318e+07 7.32070e+07 7.64955e+07 7.85946e+07
7.92512e+07 8.42030e+07 8.72462e+07 9.03168e+07 9.57509e+07
1.01409e+08 1.03211e+08 1.08419e+08 1.14338e+08 1.16909e+08
1.20802e+08 1.27706e+08 1.33505e+08 1.37732e+08 1.42259e+08
1.46667e+08 1.57041e+08 1.63103e+08 1.73926e+08 1.76456e+08
1.84724e+08 1.97311e+08 2.08515e+08 2.14848e+08 2.21159e+08
2.36805e+08 2.51518e+08 2.61028e+08 2.72899e+08 2.84524e+08
3.00498e+08 3.16539e+08 3.32012e+08 3.47813e+08 3.69477e+08
3.79136e+08 4.07090e+08 4.25618e+08 4.57712e+08 4.72444e+08
4.97197e+08 5.32304e+08 5.58369e+08 5.84079e+08 6.14556e+08
6.55879e+08 6.97132e+08 7.12500e+08 7.66058e+08 8.04142e+08
8.45850e+08 8.93874e+08 9.23554e+08 1.00171e+09 1.03412e+09
1.12025e+09 1.14719e+09 1.21987e+09 1.27615e+09 1.33111e+09
1.38238e+09 1.47550e+09 1.50922e+09 1.58358e+09 1.57323e+09
1.49526e+09 1.40991e+09 1.25957e+09 1.03508e+09 8.00378e+08
5.62677e+08 3.15212e+08 1.33957e+08 4.67261e+07 1.31714e+07
1.48423e+06 1.22319e+06 3.54331e+05 2.34464e+05 5.01317e+00
9.54696e-05 1.80451e-05 0.00000e+00 0.00000e+00 0.00000e+00
0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00 0.00000e+00
0.00000e+00 0.00000e+00
c approximate size of beam at VS is 40mmx20mm
si23 -2 2
sp23 0 1
si24 -1 1
sp24 0 1