diff --git a/e2-estia.i b/e2-estia.i index b71ebe6..70b43bc 100644 --- a/e2-estia.i +++ b/e2-estia.i @@ -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 ++++++++++++++++++++++++++++ diff --git a/mcstas_reference/VS_L_total.dat b/mcstas_reference/VS_L_total.dat new file mode 100644 index 0000000..aace003 --- /dev/null +++ b/mcstas_reference/VS_L_total.dat @@ -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 diff --git a/mcstas_reference/convert_source.py b/mcstas_reference/convert_source.py new file mode 100644 index 0000000..40a4b3a --- /dev/null +++ b/mcstas_reference/convert_source.py @@ -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 diff --git a/mcstas_source_vs.i b/mcstas_source_vs.i new file mode 100644 index 0000000..d36fc3b --- /dev/null +++ b/mcstas_source_vs.i @@ -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