mcdoc - default values for input parameters
Thomas C Hansen
hansen at ill.fr
Mon Mar 13 15:14:40 CET 2000
Dear Kristian,
I see your approach, looks good! By the way, I wrote a perl script (attached), which
handles huge command lines, as I suggested a couple of days ago. To use it, you put
all my files in a directory and type, e.g.,
ln -s mcsim.pl mcsim; chmod a+x mcsim
mcsim d20 par=ge113.par Sample_file=0
This will call 'mcrun d20', take default parameters from 'd20.par', overwrite them
with parameters eventually found in 'ge113.par', and overwrite in the end the
parameter 'Sample_file'. If any needed parameter as found in 'd20.instr' is missing,
the script will ask for it and append a corresponding entry in 'd20.par'. All these
parameters are not simply given in a command line, as this won't work for command
lines above a certain limited length, but only the values are written in the right
order in a temporary file 'd20.stdin'. Then a shell script containing 'mcrun d20 <
d20.stdin > d20.stdout' is executed. In order to suppress all messages asking for
all these parameters, the output is deviated to 'd20.stdout' and the script is
printing in the end only lines not containing these questions to the screen.
Unfortunately I haven't found yet a possibility to filter that standard output in
real time - that would be preferable. With this perl script also gscan will work,
using e.g. following command line:
gscan 10 1e5 mcsim d20.scan instr=d20 par=ge113.par Source_N_E=1,10
That's just a suggestion, it works fine for me, but maybe some future testing and
standardization has to be done to make it a secure tool for everybody handling long
lists of parameters. It will be easy to adapt the script to default values as you
suggested, just by an additional string analysis of the 'DEFINE INSTRUMENT'
following lines. In that case a default parameter file ('d20.par' in the example)
may become redundant.
Best Regards,
Thomas
Kristian Nielsen wrote:
> > Date: Mon, 13 Mar 2000 11:23:43 +0100
> > From: Thomas C Hansen <hansen at ill.fr>
>
> > Dear Kristian and Emmanuel,
> >
> > The tables in the html files coming out from mcdoc contain a column for default
> > values of input parameters. I haven't found any component alrady using this
> > possibility, so what's the foreseen syntax to fill in that column as well?
>
> For example,
>
> DEFINE COMPONENT GUIDE
> SETTING PARAMETERS(R0 = 1, ...)
>
> This would make the R0 parameter optional, with a default value of
> 1. This is already implemented in my current development version. I
> guess a similar syntax will be implemented for instrument input
> parameters.
>
> - Kristian.
--
Dr. Thomas C HANSEN, Instrument Scientist D20, Diffraction Group
Institut Max von Laue-Paul Langevin, BP 156, 38042 Grenoble Cédex 9, France
ILL4-1st floor room 100, Phone +33(0)476207044, Fax +33(0)476207648
mailto:hansen at ill.fr, http://www.ill.fr/YellowBook/D20
Private: 13, allée des étangs, 38240 Meylan, France, +33(0)476180553
-------------- next part --------------
A non-text attachment was scrubbed...
Name: d20.instr
Type: image/x-xbitmap
Size: 40898 bytes
Desc: not available
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20000313/fea05225/attachment.xbm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: mcsim.pl
Type: application/x-perl
Size: 3366 bytes
Desc: Document
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20000313/fea05225/attachment.pl>
-------------- next part --------------
Source_alpha=.1 Source_beta=.1
Source_N_E=71 Source_N_xpos=1 Source_N_xdiv=1
Source_Lmin=2.3 Source_Lmax=2.5
Filter_d=0.05
Mono_position=17.2
Mono_omega=0
Mono_foccus=0.3
Mono_mosh=10 Mono_mosv=10
Mono_a=5.65 Mono_Vc=180.362
Mono_h0=3 Mono_k0=3 Mono_l0=5
Mono_h=1 Mono_k=1 Mono_l=3
Mono_F2=536.64 Mono_B=1
Mono_xmax=0.008 Mono_ymax=0.019 Mono_dist=0.020 Mono_zmax=0.07
Mono_blades=15
Mono_harmonic=3
Mono_LIMIT=10001
Takeoff=-90
-------------- next part --------------
A non-text attachment was scrubbed...
Name: d20.par
Type: image/x-xbitmap
Size: 861 bytes
Desc: not available
URL: <http://mailman2.mcstas.org/pipermail/mcstas-users/attachments/20000313/fea05225/attachment-0001.xbm>
-------------- next part --------------
#0.700000 0.003500 0.050000 0.000100 # sample: packing radius height wall
#1.472205 0.153000 0.100000 0.053205 3.000000 # PSD: radius height pitch gap pressure
160.176 1.368 0.032 # Si powder: Vc sig_a sig_i
24
2.004 8 550.881 0.979320 1 1 1 1
3.272 12 1101.763 0.945799 1 2 2 0
4.008 8 0.000 0.919810 1 2 2 2
3.837 24 550.881 0.926239 1 3 1 1
5.043 24 550.880 0.876036 1 3 3 1
6.012 8 550.881 0.828554 1 3 3 3
4.628 6 1101.763 0.894535 1 4 0 0
5.668 24 1101.763 0.846050 1 4 2 2
6.545 12 1101.763 0.800193 1 4 4 0
6.942 24 0.000 0.778205 1 4 4 2
8.016 8 1101.763 0.715801 1 4 4 4
6.012 24 550.881 0.828554 1 5 1 1
6.845 48 550.881 0.783645 1 5 3 1
7.587 24 550.881 0.741170 1 5 3 3
8.262 24 550.882 0.700998 1 5 5 1
8.887 24 550.882 0.663003 1 5 5 3
7.317 24 1101.763 0.756822 1 6 2 0
7.674 24 0.000 0.736026 1 6 2 2
8.658 48 1101.763 0.677004 1 6 4 2
9.540 24 0.000 0.622715 1 6 4 4
8.262 24 550.882 0.700998 1 7 1 1
8.887 48 550.881 0.663003 1 7 3 1
9.470 24 550.883 0.627067 1 7 3 3
9.256 6 1101.763 0.640309 1 8 0 0
1889.494 133.139 19.677 YIG
210 YIG
0.719 12 0.000 0.998322
1.016 6 0.000 0.996636
1.245 24 203.580 0.984168
1.438 12 3016.019 0.986776
1.607 24 0.000 0.992151
1.761 8 0.000 0.989722
1.902 48 1506.029 0.978739
2.033 6 98281.828 0.986199
2.156 12 0.000 0.986071
2.156 24 0.000 0.981520
2.273 24 54691.691 0.986475
2.384 24 36367.449 0.977947
2.490 24 18260.684 0.979746
2.592 48 13085.104 0.977365
2.592 24 0.000 0.980939
2.784 48 9884.264 0.981401
2.875 12 11381.769 0.984795
2.963 24 0.000 0.972488
2.963 24 0.000 0.974189
3.049 24 0.000 0.970052
3.049 6 0.000 0.980046
3.133 48 4738.101 0.980699
3.133 24 26526.645 0.972831
3.214 24 17280.061 0.952795
3.294 48 1072.821 0.963691
3.371 24 0.000 0.964791
3.447 48 5640.986 0.945413
3.521 8 108324.977 0.971225
3.594 48 710.790 0.956927
3.594 12 0.000 0.972702
3.594 24 0.000 0.964389
3.665 24 190097.016 0.960594
3.735 24 8854.975 0.967531
3.735 24 2605.230 0.928074
3.735 48 136.468 0.844741
3.803 48 127188.094 0.952464
3.871 24 0.000 0.959452
4.002 48 48.465 1.173685
4.002 48 14212.914 0.934375
4.066 6 244666.688 0.951433
4.129 24 0.000 0.941082
4.129 48 7215.575 0.943539
4.129 24 0.000 0.937901
4.191 24 0.000 0.946819
4.191 24 17019.494 0.941878
4.252 48 8083.322 0.921666
4.313 12 734.671 1.049026
4.313 24 10628.211 0.910545
4.372 48 7646.340 0.936915
4.372 24 0.000 0.915093
4.372 48 776.047 0.936915
4.431 24 0.000 0.930654
4.489 48 8357.111 0.913423
4.546 24 99628.523 0.932912
4.602 24 0.000 0.925655
4.602 24 0.000 0.927404
4.658 48 41391.684 0.946833
4.713 24 6454.844 0.952557
4.713 48 20689.980 0.912980
4.713 48 1554.205 0.876190
4.768 24 82241.906 0.926566
4.822 48 21384.988 0.923807
4.822 48 15318.223 0.923807
4.822 24 0.000 0.925684
4.928 48 504.394 0.824814
4.928 48 19.874 1.482843
4.980 24 16359.466 0.889442
5.031 12 0.000 0.922682
5.031 48 39592.711 0.917322
5.031 48 25325.520 0.917322
5.082 24 129.762 0.915708
5.082 6 0.000 0.924715
5.133 24 1579.174 0.975029
5.133 24 20095.832 0.897369
5.183 48 673.349 1.077036
5.183 24 43752.160 0.932161
5.233 48 10300.698 0.910883
5.233 24 0.000 0.893915
5.282 8 0.000 0.903652
5.282 24 0.000 0.903433
5.330 48 11530.723 0.931612
5.330 48 131.671 0.698822
5.330 48 1564.328 0.844278
5.426 24 0.000 0.903811
5.426 24 0.000 0.911782
5.426 48 5560.977 0.904489
5.474 48 43032.797 0.926892
5.474 24 121046.172 0.917166
5.521 48 869.710 0.995956
5.521 24 1071.323 0.986376
5.567 48 100115.828 0.901113
5.614 48 1334.206 0.898139
5.614 48 1100.267 0.898139
5.614 24 0.000 0.894209
5.705 48 44.953 0.516079
5.705 48 11047.518 0.922668
5.705 48 4519.932 0.938451
5.750 12 231623.359 0.905655
5.795 24 0.000 0.898561
5.795 24 0.000 0.895735
5.839 24 0.000 0.902606
5.839 24 0.000 0.884242
5.883 24 24914.727 0.869469
5.883 48 7679.861 0.854214
5.883 48 125.314 1.182243
5.883 48 1178.122 0.802006
5.927 24 7067.095 0.950199
5.927 24 11585.631 0.936204
5.970 48 462.440 0.885574
5.970 48 5778.562 0.885574
6.014 48 22267.492 0.884016
6.056 48 486.486 1.034189
6.099 24 72231.336 0.882867
6.099 6 19129.297 0.884716
6.141 48 72.866 0.879357
6.141 48 27934.947 0.879357
6.141 48 5137.507 0.879357
6.141 24 0.000 0.883748
6.141 24 0.000 0.887918
6.183 24 66256.578 0.901881
6.225 24 1145.936 0.779438
6.225 48 2368.012 0.946966
6.225 48 2.107 0.093289
6.266 48 84037.055 0.876638
6.266 24 134249.734 0.876239
6.307 48 1321.134 0.873184
6.307 48 23028.475 0.873184
6.388 48 378.576 0.698084
6.388 48 1202.565 0.771261
6.429 24 4163.282 0.777967
6.469 24 0.000 0.847652
6.469 12 0.000 0.913683
6.469 48 1046.849 0.867055
6.469 24 0.000 0.732765
6.509 24 0.000 0.861072
6.509 24 10020.346 0.865529
6.509 48 100.722 0.865529
6.548 48 5719.824 0.913443
6.548 24 10481.625 0.900393
6.548 48 34038.199 0.884106
6.587 48 2212.532 0.734415
6.627 48 8340.317 0.860968
6.627 24 0.000 0.864604
6.627 48 1144.634 0.860968
6.627 24 0.000 0.866319
6.665 24 0.000 0.861820
6.704 48 4163.081 0.799453
6.704 48 2879.308 0.931003
6.704 48 18606.156 0.830018
6.743 24 31417.455 0.935357
6.781 24 0.000 0.851566
6.781 48 19656.873 0.854924
6.781 24 0.000 0.899097
6.819 48 154874.438 0.872072
6.819 24 45146.090 0.888126
6.856 48 1.584 6.976236
6.856 48 51.244 1.500796
6.856 48 1127.712 0.974933
6.894 48 50511.836 0.853395
6.931 48 34190.152 0.848923
6.931 48 2164.648 0.848923
6.931 48 524.398 0.848923
7.006 48 10537.984 0.806079
7.042 8 279514.156 0.860435
7.079 48 4450.304 0.842964
7.079 48 91.726 0.842964
7.079 24 0.000 0.833867
7.079 48 236.346 0.842964
7.079 48 13131.040 0.842964
7.079 24 0.000 0.818209
7.115 48 2832.459 0.841480
7.115 6 0.000 0.853769
7.152 24 55446.820 0.821886
7.152 24 10563.778 0.798793
7.152 48 29967.029 0.815409
7.152 24 871.975 0.991503
7.188 48 2430.303 0.697682
7.188 12 61323.203 0.809452
7.188 24 15699.059 0.781564
7.223 24 0.000 0.866574
7.223 48 26204.113 0.837046
7.259 24 0.000 0.837881
7.259 24 0.000 0.829693
7.295 48 11099.855 0.792550
7.295 48 1751.051 0.731513
7.295 48 1620.526 0.947878
7.295 48 5809.077 0.776945
7.295 48 3604.480 0.909587
7.330 24 73445.188 0.835414
7.365 48 15921.470 0.831170
7.365 48 5873.266 0.831170
7.400 48 76166.891 0.860376
7.400 24 54682.270 0.865962
7.435 48 19258.697 0.795594
7.435 24 2714.380 0.742726
7.470 24 165646.812 0.828709
7.470 24 74290.812 0.829656
7.470 48 151288.562 0.828798
7.504 48 10344.729 0.825335
7.504 48 5411.900 0.825336
7.504 24 0.000 0.822316
7.573 48 1254.346 0.694604
7.573 48 148.693 0.481002
7.573 48 30448.402 0.795615
7.607 48 1966.262 1.010599
7.640 24 0.000 0.827541
7.640 48 6363.631 0.819542
7.640 24 0.000 0.837935
7.674 24 0.000 0.824846
7.674 24 0.000 0.818697
160.176 1.368 0.032 Si
24
2.004 8 550.837 0.940795
3.272 12 1101.675 0.849807
4.008 8 0.000 0.783395
3.837 24 550.837 0.799495
5.043 24 550.836 0.679416
6.012 8 550.837 0.577373
4.628 6 1101.675 0.722172
5.668 24 1101.675 0.613707
6.545 12 1101.675 0.521533
6.942 24 0.000 0.480775
8.016 8 1101.675 0.376637
6.012 24 550.837 0.577373
6.845 48 550.837 0.490656
7.587 24 550.837 0.416963
8.262 24 550.838 0.354338
8.887 24 550.838 0.301119
7.317 24 1101.675 0.443202
7.674 24 0.000 0.408566
8.658 48 1101.675 0.320069
9.540 24 0.000 0.250740
8.262 24 550.837 0.354338
8.887 48 550.836 0.301119
9.470 24 550.839 0.255893
9.256 6 1101.675 0.271997
24.265 5.120 0.800 Fe
10 Fe
3.100 12 357.210 0.940960 1
4.384 6 357.210 0.885405 1
5.369 24 357.210 0.833130 1
6.200 12 357.210 0.783942 1
6.932 24 357.210 0.737658 1
7.593 8 357.210 0.694106 1
8.202 48 357.210 0.653126 1
8.768 6 357.210 0.614565 1
9.300 24 357.210 0.578281 1
9.300 12 357.210 0.578281 1
-------------- next part --------------
# 0.600000 0.003500 0.050000 0.000100 # sample: packing radius height wall
# 1.472205 0.153000 0.100000 0.053205 3.000000 # PSD: radius height pitch gap pressure
1889.494 133.139 19.677 # YIG powder: Vc sig_a sig_i
400 YIG
0.719 12 0.000 0.998322 1 1 0
1.016 6 0.000 0.996636 2 0 0
1.245 24 203.580 0.984168 2 1 1
1.438 12 3016.019 0.986776 2 2 0
1.607 24 0.000 0.992151 3 1 0
1.761 8 0.000 0.989722 2 2 2
1.902 48 1506.029 0.978739 3 2 1
2.033 6 98281.828 0.986199 4 0 0
2.156 12 0.000 0.986071 3 3 0
2.156 24 0.000 0.981520 4 1 1
2.273 24 54691.691 0.986475 4 2 0
2.384 24 36367.449 0.977947 3 3 2
2.490 24 18260.684 0.979746 4 2 2
2.592 48 13085.104 0.977365 4 3 1
2.592 24 0.000 0.980939 5 1 0
2.784 48 9884.264 0.981401 5 2 1
2.875 12 11381.769 0.984795 4 4 0
2.963 24 0.000 0.972488 4 3 3
2.963 24 0.000 0.974189 5 3 0
3.049 24 0.000 0.970052 4 4 2
3.049 6 0.000 0.980046 6 0 0
3.133 48 4738.101 0.980699 5 3 2
3.133 24 26526.645 0.972831 6 1 1
3.214 24 17280.061 0.952795 6 2 0
3.294 48 1072.821 0.963691 5 4 1
3.371 24 0.000 0.964791 6 2 2
3.447 48 5640.986 0.945413 6 3 1
3.521 8 108324.977 0.971225 4 4 4
3.594 48 710.790 0.956927 5 4 3
3.594 12 0.000 0.972702 5 5 0
3.594 24 0.000 0.964389 7 1 0
3.665 24 190097.016 0.960594 6 4 0
3.735 24 8854.975 0.967531 5 5 2
3.735 24 2605.230 0.928074 6 3 3
3.735 48 136.468 0.844741 7 2 1
3.803 48 127188.094 0.952464 6 4 2
3.871 24 0.000 0.959452 7 3 0
4.002 48 48.465 1.173685 6 5 1
4.002 48 14212.914 0.934375 7 3 2
4.066 6 244666.688 0.951433 8 0 0
4.129 24 0.000 0.941082 5 5 4
4.129 48 7215.575 0.943539 7 4 1
4.129 24 0.000 0.937901 8 1 1
4.191 24 0.000 0.946819 6 4 4
4.191 24 17019.494 0.941878 8 2 0
4.252 48 8083.322 0.921666 6 5 3
4.313 12 734.671 1.049026 6 6 0
4.313 24 10628.211 0.910545 8 2 2
4.372 48 7646.340 0.936915 7 4 3
4.372 24 0.000 0.915093 7 5 0
4.372 48 776.047 0.936915 8 3 1
4.431 24 0.000 0.930654 6 6 2
4.489 48 8357.111 0.913423 7 5 2
4.546 24 99628.523 0.932912 8 4 0
4.602 24 0.000 0.925655 8 3 3
4.602 24 0.000 0.927404 9 1 0
4.658 48 41391.684 0.946833 8 4 2
4.713 24 6454.844 0.952557 6 5 5
4.713 48 20689.980 0.912980 7 6 1
4.713 48 1554.205 0.876190 9 2 1
4.768 24 82241.906 0.926566 6 6 4
4.822 48 21384.988 0.923807 7 5 4
4.822 48 15318.223 0.923807 8 5 1
4.822 24 0.000 0.925684 9 3 0
4.928 48 504.394 0.824814 7 6 3
4.928 48 19.874 1.482843 9 3 2
4.980 24 16359.466 0.889442 8 4 4
5.031 12 0.000 0.922682 7 7 0
5.031 48 39592.711 0.917322 8 5 3
5.031 48 25325.520 0.917322 9 4 1
5.082 24 129.762 0.915708 8 6 0
5.082 6 0.000 0.924715 10 0 0
5.133 24 1579.174 0.975029 7 7 2
5.133 24 20095.832 0.897369 10 1 1
5.183 48 673.349 1.077036 8 6 2
5.183 24 43752.160 0.932161 10 2 0
5.233 48 10300.698 0.910883 9 4 3
5.233 24 0.000 0.893915 9 5 0
5.282 8 0.000 0.903652 6 6 6
5.282 24 0.000 0.903433 10 2 2
5.330 48 11530.723 0.931612 7 6 5
5.330 48 131.671 0.698822 9 5 2
5.330 48 1564.328 0.844278 10 3 1
5.426 24 0.000 0.903811 7 7 4
5.426 24 0.000 0.911782 8 5 5
5.426 48 5560.977 0.904489 8 7 1
5.474 48 43032.797 0.926892 8 6 4
5.474 24 121046.172 0.917166 10 4 0
5.521 48 869.710 0.995956 9 6 1
5.521 24 1071.323 0.986376 10 3 3
5.567 48 100115.828 0.901113 10 4 2
5.614 48 1334.206 0.898139 8 7 3
5.614 48 1100.267 0.898139 9 5 4
5.614 24 0.000 0.894209 11 1 0
5.705 48 44.953 0.516079 9 6 3
5.705 48 11047.518 0.922668 10 5 1
5.705 48 4519.932 0.938451 11 2 1
5.750 12 231623.359 0.905655 8 8 0
5.795 24 0.000 0.898561 9 7 0
5.795 24 0.000 0.895735 11 3 0
5.839 24 0.000 0.902606 8 8 2
5.839 24 0.000 0.884242 10 4 4
5.883 24 24914.727 0.869469 7 7 6
5.883 48 7679.861 0.854214 9 7 2
5.883 48 125.314 1.182243 10 5 3
5.883 48 1178.122 0.802006 11 3 2
5.927 24 7067.095 0.950199 8 6 6
5.927 24 11585.631 0.936204 10 6 0
5.970 48 462.440 0.885574 8 7 5
5.970 48 5778.562 0.885574 11 4 1
6.014 48 22267.492 0.884016 10 6 2
6.056 48 486.486 1.034189 9 6 5
6.099 24 72231.336 0.882867 8 8 4
6.099 6 19129.297 0.884716 12 0 0
6.141 48 72.866 0.879357 9 7 4
6.141 48 27934.947 0.879357 9 8 1
6.141 48 5137.507 0.879357 11 4 3
6.141 24 0.000 0.883748 11 5 0
6.141 24 0.000 0.887918 12 1 1
6.183 24 66256.578 0.901881 12 2 0
6.225 24 1145.936 0.779438 10 5 5
6.225 48 2368.012 0.946966 10 7 1
6.225 48 2.107 0.093289 11 5 2
6.266 48 84037.055 0.876638 10 6 4
6.266 24 134249.734 0.876239 12 2 2
6.307 48 1321.134 0.873184 9 8 3
6.307 48 23028.475 0.873184 12 3 1
6.388 48 378.576 0.698084 10 7 3
6.388 48 1202.565 0.771261 11 6 1
6.429 24 4163.282 0.777967 12 4 0
6.469 24 0.000 0.847652 8 7 7
6.469 12 0.000 0.913683 9 9 0
6.469 48 1046.849 0.867055 11 5 4
6.469 24 0.000 0.732765 12 3 3
6.509 24 0.000 0.861072 8 8 6
6.509 24 10020.346 0.865529 10 8 0
6.509 48 100.722 0.865529 12 4 2
6.548 48 5719.824 0.913443 9 7 6
6.548 24 10481.625 0.900393 9 9 2
6.548 48 34038.199 0.884106 11 6 3
6.587 48 2212.532 0.734415 10 8 2
6.627 48 8340.317 0.860968 9 8 5
6.627 24 0.000 0.864604 11 7 0
6.627 48 1144.634 0.860968 12 5 1
6.627 24 0.000 0.866319 13 1 0
6.665 24 0.000 0.861820 10 6 6
6.704 48 4163.081 0.799453 10 7 5
6.704 48 2879.308 0.931003 11 7 2
6.704 48 18606.156 0.830018 13 2 1
6.743 24 31417.455 0.935357 12 4 4
6.781 24 0.000 0.851566 9 9 4
6.781 48 19656.873 0.854924 12 5 3
6.781 24 0.000 0.899097 13 3 0
6.819 48 154874.438 0.872072 10 8 4
6.819 24 45146.090 0.888126 12 6 0
6.856 48 1.584 6.976236 10 9 1
6.856 48 51.244 1.500796 11 6 5
6.856 48 1127.712 0.974933 13 3 2
6.894 48 50511.836 0.853395 12 6 2
6.931 48 34190.152 0.848923 11 7 4
6.931 48 2164.648 0.848923 11 8 1
6.931 48 524.398 0.848923 13 4 1
7.006 48 10537.984 0.806079 10 9 3
7.042 8 279514.156 0.860435 8 8 8
7.079 48 4450.304 0.842964 9 8 7
7.079 48 91.726 0.842964 11 8 3
7.079 24 0.000 0.833867 12 5 5
7.079 48 236.346 0.842964 12 7 1
7.079 48 13131.040 0.842964 13 4 3
7.079 24 0.000 0.818209 13 5 0
7.115 48 2832.459 0.841480 12 6 4
7.115 6 0.000 0.853769 14 0 0
7.152 24 55446.820 0.821886 9 9 6
7.152 24 10563.778 0.798793 10 7 7
7.152 48 29967.029 0.815409 13 5 2
7.152 24 871.975 0.991503 14 1 1
7.188 48 2430.303 0.697682 10 8 6
7.188 12 61323.203 0.809452 10 10 0
7.188 24 15699.059 0.781564 14 2 0
7.223 24 0.000 0.866574 11 9 0
7.223 48 26204.113 0.837046 12 7 3
7.259 24 0.000 0.837881 10 10 2
7.259 24 0.000 0.829693 14 2 2
7.295 48 11099.855 0.792550 10 9 5
7.295 48 1751.051 0.731513 11 7 6
7.295 48 1620.526 0.947878 11 9 2
7.295 48 5809.077 0.776945 13 6 1
7.295 48 3604.480 0.909587 14 3 1
7.330 24 73445.188 0.835414 12 8 0
7.365 48 15921.470 0.831170 11 8 5
7.365 48 5873.266 0.831170 13 5 4
7.400 48 76166.891 0.860376 12 8 2
7.400 24 54682.270 0.865962 14 4 0
7.435 48 19258.697 0.795594 13 6 3
7.435 24 2714.380 0.742726 14 3 3
7.470 24 165646.812 0.828709 10 10 4
7.470 24 74290.812 0.829656 12 6 6
7.470 48 151288.562 0.828798 14 4 2
7.504 48 10344.729 0.825335 11 9 4
7.504 48 5411.900 0.825336 12 7 5
7.504 24 0.000 0.822316 13 7 0
7.573 48 1254.346 0.694604 11 10 1
7.573 48 148.693 0.481002 13 7 2
7.573 48 30448.402 0.795615 14 5 1
7.607 48 1966.262 1.010599 12 8 4
7.640 24 0.000 0.827541 9 9 8
7.640 48 6363.631 0.819542 12 9 1
7.640 24 0.000 0.837935 15 1 0
7.674 24 0.000 0.824846 10 8 8
7.674 24 0.000 0.818697 14 4 4
7.708 48 3849.914 0.896828 10 9 7
7.708 48 1094.647 0.676433 11 10 3
7.708 48 79.345 0.362349 13 6 5
7.708 48 79.048 1.454426 14 5 3
7.708 48 5521.893 0.883341 15 2 1
7.741 24 74.444 0.141779 14 6 0
7.775 48 15044.647 0.813789 11 8 7
7.775 48 11309.317 0.813789 12 9 3
7.775 48 14739.329 0.813789 13 7 4
7.775 48 3247.167 0.813789 13 8 1
7.775 24 0.000 0.828548 15 3 0
7.808 24 0.000 0.844919 10 10 6
7.808 48 11997.312 0.812357 14 6 2
7.841 48 19856.273 0.775865 11 9 6
7.841 48 13406.674 0.854643 15 3 2
7.906 12 0.000 0.816166 11 11 0
7.906 24 0.000 0.830149 12 7 7
7.906 48 7061.392 0.808076 13 8 3
7.906 48 2060.587 0.808076 15 4 1
7.939 48 93269.773 0.837780 12 8 6
7.939 24 47890.504 0.850255 12 10 0
7.971 48 7796.717 0.748176 11 10 5
7.971 24 7086.558 0.745438 11 11 2
7.971 24 2328.366 0.702413 14 5 5
7.971 48 19.962 0.066294 14 7 1
8.004 48 112340.125 0.806478 12 10 2
8.004 48 53296.637 0.807681 14 6 4
8.036 48 14131.926 0.802403 12 9 5
8.036 24 0.000 0.829866 13 9 0
8.036 48 956.852 0.802403 15 4 3
8.036 24 0.000 0.785195 15 5 0
8.100 48 125.046 1.339667 13 7 6
8.100 48 2835.213 0.703768 13 9 2
8.100 48 10.602 0.000098 14 7 3
8.100 48 1178.623 0.653589 15 5 2
8.132 6 151833.891 0.825886 16 0 0
8.163 24 0.000 0.795598 11 11 4
8.163 48 1550.868 0.796771 13 8 5
8.163 24 0.000 0.872871 16 1 1
8.195 48 242.107 0.795369 12 10 4
8.195 24 33236.344 0.795369 14 8 0
8.195 24 34660.660 0.795369 16 2 0
8.227 24 12606.899 0.842757 10 9 9
8.227 48 2381.923 0.687250 15 6 1
8.258 24 864.789 0.513735 10 10 8
8.258 48 6570.184 0.908619 14 8 2
8.258 24 6778.777 0.906760 16 2 2
8.289 48 966.176 0.791177 11 9 8
8.289 48 43581.391 0.791178 12 11 1
8.289 48 7196.764 0.791178 13 9 4
8.289 48 10670.142 0.791177 15 5 4
8.289 48 3204.440 0.791177 16 3 1
8.320 24 0.000 0.808097 14 6 6
8.351 48 200.804 1.226601 11 10 7
8.351 48 2200.392 0.910606 13 10 1
8.351 48 1091.037 0.629967 14 7 5
8.351 48 3256.225 0.888213 15 6 3
8.382 24 93462.555 0.790193 12 8 8
8.382 24 75682.148 0.790548 16 4 0
8.413 48 11445.481 0.785624 12 9 7
8.413 48 142.793 0.785624 12 11 3
8.413 24 0.000 0.807413 15 7 0
8.413 24 0.000 0.793343 16 3 3
8.443 48 66963.406 0.824903 14 8 4
8.443 48 152509.391 0.811070 16 4 2
8.474 24 991.613 0.613598 11 11 6
8.474 48 1375.104 0.637806 13 10 3
8.474 48 3873.620 0.694648 14 9 1
8.474 48 6604.764 0.853982 15 7 2
8.504 48 105229.992 0.784568 12 10 6
8.535 48 20650.426 0.780109 13 8 7
8.535 48 11068.548 0.780109 16 5 1
8.595 48 1777.985 0.920500 13 9 6
8.595 48 9026.879 0.839398 14 9 3
8.595 48 1867.975 0.649371 15 6 5
8.625 12 440.601 0.382222 12 12 0
8.625 24 4123.180 0.934647 16 4 4
8.655 48 2586.402 0.774632 12 11 5
8.655 24 0.000 0.773600 13 11 0
8.655 48 3.233 0.774635 15 7 4
8.655 48 53071.383 0.774632 15 8 1
8.655 48 3789.377 0.774632 16 5 3
8.655 24 0.000 0.795470 17 1 0
8.685 24 0.000 0.771895 12 12 2
8.685 24 3643.341 0.773269 16 6 0
8.714 48 941.143 0.591423 13 10 5
8.714 48 3285.156 0.878372 13 11 2
8.714 24 39280.148 0.801991 14 7 7
8.714 48 8748.503 0.710066 17 2 1
8.744 48 1214.614 0.512271 14 8 6
8.744 24 2282.333 0.576950 14 10 0
8.744 48 57.059 2.642583 16 6 2
8.773 48 308.901 0.769195 15 8 3
8.773 24 0.000 0.819555 17 3 0
8.803 8 0.000 0.755858 10 10 10
8.803 48 454.897 0.767840 14 10 2
8.832 48 12771.266 0.820756 11 10 9
8.832 48 6171.375 0.845140 14 9 5
8.832 48 590.266 1.034670 17 3 2
8.861 24 104101.289 0.832626 12 12 4
8.890 24 0.000 0.770600 11 11 8
8.890 24 0.000 0.675269 12 9 9
8.890 48 9594.028 0.763795 13 11 4
8.890 24 0.000 0.801700 15 9 0
8.890 24 0.000 0.763028 16 5 5
8.890 48 11653.093 0.763795 16 7 1
8.890 48 960.200 0.763795 17 4 1
8.919 48 41982.527 0.818643 12 10 8
8.919 48 84537.617 0.801842 16 6 4
8.948 48 1440.681 0.931723 15 7 6
8.948 48 6684.949 0.687764 15 9 2
8.977 48 31867.904 0.765981 14 10 4
9.006 48 3182.443 0.758433 12 11 7
9.006 48 58.134 0.758433 13 9 8
9.006 48 11570.396 0.758433 13 12 1
9.006 48 5064.502 0.758433 15 8 5
9.006 48 6890.246 0.758433 16 7 3
9.006 48 12437.841 0.758433 17 4 3
9.006 24 0.000 0.781995 17 5 0
9.063 48 2402.020 0.633288 13 10 7
9.063 48 935.103 0.974633 14 11 1
9.063 48 2552.169 0.636789 17 5 2
9.092 24 489417.750 0.772806 16 8 0
9.120 48 16024.504 0.753109 13 12 3
9.120 48 2811.926 0.753109 15 9 4
9.148 24 0.000 0.772726 12 12 6
9.148 24 0.000 0.725124 14 8 8
9.148 48 1665.530 0.751784 16 8 2
9.148 6 0.000 0.769562 18 0 0
9.176 48 13771.264 0.697036 13 11 6
9.176 48 29490.500 0.713739 14 9 7
9.176 48 15184.293 0.699537 14 11 3
9.176 48 1268.854 0.940422 15 10 1
9.176 48 229.789 0.388343 17 6 1
9.176 24 5101.780 0.663727 18 1 1
9.205 24 26858.607 0.685284 16 6 6
9.205 24 1676.267 1.033228 18 2 0
9.233 48 4989.955 0.747822 16 7 5
9.233 48 25568.219 0.747822 17 5 4
9.260 48 15964.095 0.746507 14 10 6
9.260 24 0.000 0.738002 18 2 2
9.288 48 4842.825 0.841539 15 10 3
9.288 48 59.952 1.821303 17 6 3
9.288 48 842.021 0.986068 18 3 1
9.316 48 73451.273 0.748262 16 8 4
9.344 48 1779.935 0.742573 13 12 5
9.344 12 0.000 0.750385 13 13 0
9.344 48 7446.031 0.742573 15 8 7
9.344 48 64.178 0.742571 16 9 1
9.344 24 0.000 0.721414 17 7 0
9.371 24 116010.547 0.777449 14 12 0
9.371 24 84.167 2.662763 18 4 0
9.399 24 16458.018 0.792424 11 11 10
9.399 24 25.086 2.650020 13 13 2
9.399 48 4530.413 0.841504 14 11 5
9.399 48 16069.531 0.793065 15 9 6
9.399 48 967.204 0.542545 17 7 2
9.399 24 13988.126 0.796950 18 3 3
9.426 24 37421.613 0.744935 12 10 10
9.426 48 39291.531 0.744784 14 12 2
9.426 48 52253.723 0.743969 18 4 2
9.454 48 8956.533 0.737360 12 11 9
9.454 24 0.000 0.769230 15 11 0
9.454 48 6119.978 0.737360 16 9 3
9.508 48 367.693 0.425812 13 10 9
9.508 48 205.257 1.279353 15 10 5
9.508 48 18085.365 0.685595 15 11 2
9.508 48 47413.023 0.704197 17 6 5
9.508 48 11488.113 0.673343 18 5 1
9.535 24 2087.741 0.507875 12 12 8
9.562 48 8.731 0.732185 13 11 8
9.562 24 0.000 0.730241 13 13 4
9.562 24 0.000 0.779491 16 7 7
9.562 48 7435.995 0.732184 17 7 4
9.562 48 4377.982 0.732184 17 8 1
9.589 48 21.935 0.730896 14 12 4
9.589 48 592.034 0.730895 16 8 6
9.589 24 12893.613 0.730895 16 10 0
9.589 24 0.000 0.733011 18 4 4
9.616 24 394.922 0.426086 14 9 9
9.616 48 1802.463 0.577434 18 5 3
9.643 48 986.512 1.137751 14 10 8
9.643 48 4248.782 0.563507 16 10 2
9.643 24 10720.993 0.622106 18 6 0
9.670 48 30420.504 0.727044 13 12 7
9.670 48 6393.077 0.727044 15 11 4
9.670 48 3875.861 0.727044 16 9 5
9.670 48 9838.925 0.727044 17 8 3
9.670 24 0.000 0.748096 19 1 0
9.697 48 8127.695 0.725765 18 6 2
9.723 48 770.210 0.495789 14 11 7
9.723 48 2651.060 0.595845 14 13 1
9.723 48 16128.882 0.670819 19 2 1
9.776 48 3330.684 0.721940 15 9 8
9.776 48 723.887 0.721940 15 12 1
9.776 24 0.000 0.759183 17 9 0
9.776 24 0.000 0.781093 19 3 0
9.802 48 6867.757 0.886021 16 10 4
9.829 24 35440.613 0.757362 13 13 6
9.829 48 10398.315 0.790245 14 13 3
9.829 48 6696.791 0.635780 15 10 7
9.829 48 9738.723 0.792662 17 7 6
9.829 48 285.456 1.197523 17 9 2
9.829 24 14787.541 0.662560 18 5 5
9.829 48 926.501 0.969829 18 7 1
9.829 48 3458.568 0.844357 19 3 2
9.855 48 101912.961 0.722259 14 12 6
9.855 48 70842.172 0.723083 18 6 4
9.881 48 4636.409 0.716873 15 12 3
9.881 48 871.953 0.716872 16 11 1
9.881 48 2569.428 0.716873 17 8 5
9.881 48 64.828 0.716872 19 4 1
9.933 48 2296.289 0.871588 15 11 6
9.933 48 3882.000 0.833895 18 7 3
9.959 24 429715.625 0.735588 16 8 8
9.985 24 0.000 0.728021 12 11 11
9.985 48 3829.123 0.711840 16 9 7
9.985 48 2300.313 0.711840 16 11 3
9.985 48 2157.346 0.711840 17 9 4
9.985 48 5.185 0.711839 19 4 3
9.985 24 0.000 0.731896 19 5 0
10.011 24 0.000 0.729847 12 12 10
10.011 24 8407.514 0.710588 18 8 0
10.037 48 4403.742 0.604092 13 11 10
10.037 48 0.061 333.109955 14 13 5
10.037 48 5228.630 0.612431 17 10 1
10.037 48 10274.168 0.639479 19 5 2
10.063 12 11096.196 0.827805 14 14 0
10.063 48 491.430 0.266939 16 10 6
10.063 48 354.569 0.210751 18 8 2
10.088 48 23376.588 0.706843 13 12 9
10.088 48 33480.711 0.706843 15 12 5
10.088 24 0.000 0.738399 15 13 0
10.114 24 0.000 0.727591 14 10 10
10.114 24 0.000 0.748768 14 14 2
10.114 24 0.000 0.757605 18 6 6
10.139 48 7466.300 0.621679 14 11 9
10.139 48 990.249 0.957920 15 13 2
10.139 48 1190.770 0.934025 17 10 3
10.139 48 9.645 0.329466 18 7 5
10.139 48 8873.878 0.784734 19 6 1
10.165 24 135888.438 0.706895 16 12 0
10.165 6 200452.516 0.706227 20 0 0
10.190 24 0.000 0.702737 13 13 8
10.190 48 12294.665 0.701881 16 11 5
10.190 48 5908.372 0.701881 17 8 7
10.190 48 3440.059 0.701881 19 5 4
10.190 24 0.000 0.664889 20 1 1
10.215 48 103570.594 0.744078 14 12 8
10.215 48 78204.258 0.750741 16 12 2
10.215 48 108975.664 0.742971 18 8 4
10.215 24 200123.938 0.731759 20 2 0
10.241 48 16443.350 0.642387 15 10 9
10.241 48 821.729 0.463152 17 9 6
10.241 48 183.235 0.256509 18 9 1
10.241 48 11.070 4.744914 19 6 3
10.266 24 85986.914 0.703017 14 14 4
10.266 24 176294.906 0.701557 20 2 2
10.291 48 1421.877 0.696954 15 11 8
10.291 48 10078.014 0.696954 15 13 4
10.291 24 0.000 0.705679 17 11 0
10.291 24 0.000 0.727756 19 7 0
10.291 48 7.288 0.696954 20 3 1
10.341 48 211.520 1.312292 14 13 7
10.341 48 1385.313 0.912704 17 10 5
10.341 48 16359.064 0.636597 17 11 2
10.341 48 490.017 0.394712 18 9 3
10.341 48 6772.580 0.789503 19 7 2
10.366 48 1908.914 0.434801 16 12 4
10.366 24 23541.840 0.613565 20 4 0
10.391 48 5924.918 0.692061 15 12 7
10.391 24 0.000 0.799210 16 9 9
10.391 24 0.000 0.729334 20 3 3
10.416 48 6394.384 0.690844 16 10 8
10.416 48 3046.540 0.690844 20 4 2
10.440 48 2461.563 0.543780 15 14 1
10.440 24 135.264 1.506200 18 7 7
10.440 48 19656.266 0.636037 19 6 5
10.465 48 22722.738 0.606361 18 8 6
10.465 24 17729.400 0.595911 18 10 0
10.490 48 5276.253 0.687204 16 11 7
10.490 48 448.920 0.687203 16 13 1
10.490 48 8295.166 0.687203 17 11 4
10.490 48 139.147 0.687203 19 7 4
10.490 48 0.328 0.687215 19 8 1
10.490 48 3527.839 0.687203 20 5 1
10.514 24 0.000 0.727619 14 14 6
10.514 48 48499.012 0.685994 18 10 2
10.539 48 647.314 0.413689 15 13 6
10.539 48 26527.465 0.637960 15 14 3
10.539 48 487.215 0.378284 18 9 5
10.563 8 52768.902 0.807684 12 12 12
10.563 24 92219.422 0.776506 20 4 4
10.588 48 5757.705 0.682379 13 12 11
10.588 48 9340.360 0.682379 16 13 3
10.588 48 710.295 0.682379 17 9 8
10.588 48 10436.022 0.682379 17 12 1
10.588 48 27884.703 0.682379 19 8 3
10.588 48 5.528 0.682372 20 5 3
10.612 48 62078.328 0.740588 16 12 6
10.612 24 2835.733 0.980520 20 6 0
10.637 24 336.916 1.175382 13 13 10
10.637 24 33559.977 0.637768 14 11 11
10.637 48 930.464 0.961942 17 10 7
10.661 48 98878.328 0.683607 14 12 10
10.661 48 42282.145 0.686166 18 10 4
10.661 48 68886.531 0.684564 20 6 2
10.685 48 1498.741 0.677589 17 12 3
10.685 24 0.000 0.698189 19 9 0
10.685 24 0.000 0.703940 21 1 0
10.733 48 2054.250 0.511190 14 13 9
10.733 48 13512.246 0.608545 15 11 10
10.733 48 874.002 0.433063 15 14 5
10.733 48 8341.328 0.765061 17 11 6
10.733 48 3496.692 0.547453 18 11 1
10.733 48 5957.996 0.576137 19 7 6
10.733 48 1422.699 0.480880 19 9 2
10.733 48 29.973 2.908144 21 2 1
10.781 48 1753.692 0.672832 15 12 9
10.781 12 0.000 0.664564 15 15 0
10.781 48 1753.916 0.672832 16 13 5
10.781 48 40.492 0.672836 19 8 5
10.781 24 0.000 0.690230 20 5 5
10.781 48 3651.796 0.672833 20 7 1
10.781 24 0.000 0.707799 21 3 0
10.805 24 53989.008 0.671649 16 14 0
10.805 24 0.000 0.643831 18 8 8
10.805 48 9669.600 0.671649 20 6 4
10.829 24 13108.416 0.602064 15 15 2
10.829 48 7.145 1.037184 18 9 7
10.829 48 5454.835 0.565999 18 11 3
10.829 48 11229.632 0.748505 21 3 2
10.853 24 15315.414 0.565674 14 14 8
10.853 24 22509.879 0.761307 16 10 10
10.853 48 3599.667 0.465110 16 14 2
10.877 48 18971.936 0.668109 15 13 8
10.877 48 6143.480 0.668110 16 11 9
10.877 48 6768.690 0.668109 17 12 5
10.877 24 0.000 0.688241 17 13 0
10.877 48 21970.449 0.668109 19 9 4
10.877 48 6369.703 0.668110 20 7 3
10.877 48 1052.865 0.668110 21 4 1
10.900 48 3017.939 0.666934 18 10 6
10.924 48 19252.697 0.725662 17 13 2
10.924 48 25930.057 0.616213 19 10 1
10.948 48 84530.609 0.670057 16 12 8
10.948 24 97406.680 0.669681 20 8 0
10.971 24 0.000 0.685600 15 15 4
10.971 48 10599.709 0.663419 21 4 3
10.971 24 0.000 0.643299 21 5 0
10.995 48 100173.562 0.711064 16 14 4
10.995 24 166767.672 0.699932 18 12 0
10.995 48 159852.547 0.700750 20 8 2
11.018 48 1658.682 0.878372 15 14 7
11.018 48 17370.697 0.600177 17 10 9
11.018 48 597.529 0.367525 18 11 5
11.018 48 997.503 0.426396 19 10 3
11.018 48 3084.364 0.521335 21 5 2
11.042 48 106769.031 0.664861 18 12 2
11.042 24 137048.594 0.664281 20 6 6
11.065 48 2199.089 0.658762 16 13 7
11.065 48 22348.607 0.658762 17 11 8
11.065 48 5517.560 0.658762 17 13 4
11.065 48 2181.360 0.658762 19 8 7
11.065 48 792.531 0.658762 20 7 5
11.112 48 94.408 1.790153 19 9 6
11.112 48 8604.981 0.749077 21 6 1
11.135 48 12669.628 0.783483 20 8 4
11.158 24 0.000 0.646298 13 13 12
11.158 48 9421.717 0.654138 16 15 1
11.158 48 210.036 0.654139 17 12 7
11.158 24 0.000 0.710653 19 11 0
11.158 48 92.088 0.654137 20 9 1
11.158 48 2810.931 0.654138 21 5 4
11.181 24 0.000 0.727377 14 12 12
11.181 48 5305.505 0.652987 18 12 4
11.181 6 0.000 0.654133 22 0 0
11.204 48 15691.037 0.720583 14 13 11
11.204 24 33526.344 0.606793 15 15 6
11.204 48 3536.515 0.800657 17 14 1
11.204 24 287.097 0.250511 18 9 9
11.204 48 7806.454 0.750322 19 10 5
11.204 48 812.864 0.390627 19 11 2
11.204 48 8561.631 0.745729 21 6 3
11.204 24 20528.598 0.594559 22 1 1
11.227 48 14.069 2.231752 16 14 6
11.227 48 3712.752 0.899219 18 10 8
11.227 24 759.892 1.253699 22 2 0
11.250 48 487.995 0.649546 15 12 11
11.250 48 2527.021 0.649546 16 15 3
11.250 48 4456.323 0.649545 20 9 3
11.250 24 0.000 0.652303 21 7 0
11.273 24 0.000 0.789600 14 14 10
11.273 24 0.000 0.653409 22 2 2
11.296 48 706.681 0.367562 15 13 10
11.296 48 297.195 0.248811 17 13 6
11.296 48 7711.250 0.554293 17 14 3
11.296 48 2570.974 0.490818 18 11 7
11.296 48 8067.584 0.745120 18 13 1
11.296 48 22961.691 0.704436 21 7 2
11.296 48 87.333 0.057870 22 3 1
11.342 24 0.000 0.685507 16 11 11
11.342 48 332.649 0.644986 19 11 4
11.342 24 0.000 0.746241 20 7 7
11.364 48 159977.438 0.683996 16 12 10
11.364 48 130867.500 0.688308 20 8 6
11.364 24 72803.609 0.703794 20 10 0
11.364 24 23748.637 0.750559 22 4 0
11.387 48 1212.476 0.912182 15 14 9
11.387 48 14.392 4.882132 18 13 3
11.387 48 1769.587 0.455264 21 6 5
11.387 24 51470.301 0.681020 22 3 3
11.410 48 99131.141 0.647012 18 12 6
11.410 48 29163.473 0.651607 20 10 2
11.410 48 92906.383 0.647192 22 4 2
11.432 48 10999.578 0.640458 16 13 9
11.432 48 4368.299 0.640459 16 15 5
11.432 48 31403.549 0.640459 19 9 8
11.432 48 119.824 0.640458 19 12 1
11.432 48 2748.602 0.640458 20 9 5
11.432 48 346.130 0.640457 21 7 4
11.432 48 8339.304 0.640458 21 8 1
11.478 48 5003.653 0.521834 17 11 10
11.478 48 3454.537 0.793780 17 14 5
11.478 48 4331.639 0.776329 19 10 7
11.478 48 503.373 0.311025 22 5 1
11.500 12 218894.016 0.675633 16 16 0
11.522 24 0.000 0.639138 15 15 8
11.522 48 13001.710 0.635962 17 12 9
11.522 24 0.000 0.631716 17 15 0
11.522 48 1780.457 0.635963 19 12 3
11.522 48 6706.062 0.635963 21 8 3
11.545 48 8059.045 0.634844 16 14 8
11.545 24 0.000 0.604006 16 16 2
11.545 48 8097.510 0.634843 20 10 4
11.545 24 0.000 0.657954 22 4 4
11.567 48 3729.106 0.498882 17 15 2
11.567 48 8971.104 0.729207 18 13 5
11.567 48 7212.025 0.535136 19 11 6
11.567 48 5994.438 0.526029 22 5 3
11.590 24 1666.043 1.035919 18 14 0
11.590 24 518.950 0.157468 22 6 0
11.612 48 17038.049 0.631498 17 13 8
11.612 48 2803.506 0.631498 20 11 1
11.612 24 0.000 0.667525 21 9 0
11.634 24 0.000 0.638608 18 10 10
11.634 48 22130.412 0.630387 18 14 2
11.634 48 1981.995 0.630388 22 6 2
11.656 48 116.182 0.077477 18 11 9
11.656 48 11704.309 0.550516 21 7 6
11.656 48 4335.886 0.502654 21 9 2
11.678 24 116167.578 0.633385 16 16 4
11.678 24 107535.102 0.633591 20 8 8
11.700 48 4930.005 0.627065 16 15 7
11.700 48 4708.191 0.627065 17 15 4
11.700 48 28224.625 0.627065 19 12 5
11.700 24 0.000 0.596998 19 13 0
11.700 48 2646.947 0.627066 20 9 7
11.700 48 918.285 0.627066 20 11 3
11.700 48 9310.486 0.627065 21 8 5
11.700 24 0.000 0.611802 23 1 0
11.723 48 74428.492 0.687569 18 12 8
11.745 24 36005.445 0.578920 14 13 13
11.745 48 21234.656 0.565384 17 14 7
11.745 48 5825.659 0.746672 19 13 2
11.745 24 23817.857 0.568623 22 5 5
11.745 48 2404.847 0.819144 22 7 1
11.745 48 14864.012 0.554119 23 2 1
11.766 24 117187.000 0.629019 14 14 12
11.766 48 163276.672 0.628214 18 14 4
11.766 48 51583.629 0.631695 20 10 6
11.766 48 66068.156 0.630769 22 6 4
11.788 48 2424.192 0.622663 15 13 12
11.788 48 12175.204 0.622663 21 9 4
11.788 24 0.000 0.684423 23 3 0
11.832 48 5088.439 0.752444 15 14 11
11.832 48 461.719 0.273537 18 13 7
11.832 48 9.878 6.746250 19 10 9
11.832 48 13722.126 0.546340 21 10 1
11.832 48 9360.422 0.716544 22 7 3
11.832 48 285.144 0.203293 23 3 2
11.854 24 245.206 1.926058 16 12 12
11.854 24 5207.864 0.841583 20 12 0
11.876 48 2550.054 0.618292 16 13 11
11.876 48 1977.929 0.618292 17 16 1
11.876 48 31793.830 0.618292 19 11 8
11.876 48 2590.990 0.618292 19 13 4
11.876 48 4312.363 0.618292 20 11 5
11.876 48 68.349 0.618293 23 4 1
11.897 24 0.000 0.601473 16 16 6
11.897 48 325.180 0.617205 20 12 2
11.897 24 7184.949 0.617204 22 8 0
11.919 24 2757.569 0.800199 15 15 10
11.919 48 6327.436 0.507762 17 15 6
11.919 48 4055.922 0.482407 18 15 1
11.919 48 242.952 0.172778 21 10 3
11.941 48 7888.565 0.458459 16 14 10
11.941 48 109.173 2.873094 22 8 2
11.962 48 10852.533 0.613952 17 12 11
11.962 48 16083.902 0.613952 17 16 3
11.962 48 21445.627 0.613952 19 12 7
11.962 48 1.472 0.613933 21 8 7
11.962 48 37939.918 0.613952 23 4 3
11.962 24 0.000 0.616921 23 5 0
11.984 48 1296.799 0.612872 18 14 6
11.984 24 0.000 0.712688 22 6 6
12.006 48 82.405 2.026851 17 13 10
12.006 48 18215.365 0.546157 18 15 3
12.006 48 9411.382 0.709300 19 14 1
12.006 48 2181.693 0.822372 21 9 6
12.006 48 335.092 1.210647 22 7 5
12.006 48 8659.716 0.517815 23 5 2
12.027 48 95207.109 0.718642 20 12 4
12.048 48 62.905 0.609640 16 15 9
12.048 24 0.000 0.614455 20 9 9
12.048 24 0.000 0.621012 21 11 0
12.070 48 20326.646 0.733442 20 10 8
12.070 48 105682.531 0.661900 22 8 4
12.091 48 21634.326 0.671532 17 14 9
12.091 24 972.384 0.937657 18 11 11
12.091 48 1994.939 0.830481 19 13 6
12.091 48 0.028 1184.492188 19 14 3
12.091 48 3860.352 0.764279 21 10 5
12.091 48 3519.251 0.462570 21 11 2
12.091 48 775.554 0.981970 22 9 1
12.091 48 10175.889 0.519880 23 6 1
12.113 48 76630.094 0.613261 18 12 10
12.134 48 3683.249 0.605362 17 16 5
12.134 48 3521.664 0.605362 20 11 7
12.134 48 124.058 0.605363 20 13 1
12.134 48 3462.126 0.605362 23 5 4
12.176 48 55.890 0.000408 18 13 9
12.176 48 97.423 1.905106 18 15 5
12.176 48 14498.582 0.528827 22 9 3
12.176 48 4836.832 0.477510 23 6 3
12.198 24 146829.391 0.652999 16 16 8
12.198 6 48972.145 0.691485 24 0 0
12.219 48 2249.633 0.601113 17 15 8
12.219 12 0.000 0.648104 17 17 0
12.219 48 731.749 0.601113 20 13 3
12.219 48 2019.349 0.601113 21 11 4
12.219 24 0.000 0.620825 23 7 0
12.219 24 0.000 0.585538 24 1 1
12.240 24 2309.716 0.600055 18 16 0
12.240 48 12503.156 0.600055 20 12 6
12.240 24 392.155 0.600055 24 2 0
12.261 24 8259.682 0.705990 17 17 2
12.261 48 1170.735 0.357480 19 11 10
12.261 48 3571.711 0.453139 19 14 5
12.261 24 23870.430 0.540168 22 7 7
12.261 48 11697.327 0.515888 23 7 2
12.282 48 1961.352 0.301177 18 14 8
12.282 48 2399.438 0.325265 18 16 2
12.282 48 8974.643 0.446666 22 8 6
12.282 24 20183.016 0.494621 22 10 0
12.282 24 4232.669 0.857612 24 2 2
12.303 48 8762.402 0.596893 19 12 9
12.303 24 0.000 0.652194 19 15 0
12.303 48 16475.064 0.596893 21 9 8
12.303 48 21749.035 0.596893 21 12 1
12.303 48 0.566 0.596886 24 3 1
12.324 8 0.000 0.643371 14 14 14
12.324 48 12293.235 0.595843 22 10 2
12.345 48 147.154 1.621171 15 14 13
12.345 48 1708.446 0.389242 19 15 2
12.345 48 516.791 1.079764 21 10 7
12.345 48 11336.400 0.509842 22 9 5
12.345 48 18111.604 0.527044 23 6 5
12.366 24 58389.941 0.601851 24 4 0
12.387 24 0.000 0.671713 15 15 12
12.387 24 0.000 0.735260 16 13 13
12.387 48 12095.004 0.592703 17 16 7
12.387 24 0.000 0.624602 17 17 4
12.387 48 6109.655 0.592703 19 13 8
12.387 48 1167.612 0.592703 20 13 5
12.387 48 486.571 0.592703 21 12 3
12.387 48 968.762 0.592702 23 7 4
12.387 48 2560.149 0.592703 23 8 1
12.387 24 0.000 0.590848 24 3 3
12.408 48 33698.273 0.690808 16 14 12
12.408 48 106267.023 0.646548 18 16 4
12.408 24 27260.693 0.702372 20 14 0
12.408 48 100993.477 0.647996 24 4 2
12.428 48 8679.107 0.493428 18 15 7
12.428 48 9943.497 0.499549 21 11 6
12.449 24 29955.414 0.601035 20 10 10
12.449 48 71944.438 0.596959 20 14 2
12.449 48 57275.621 0.597853 22 10 4
12.470 48 789.999 0.588542 16 15 11
12.470 48 10761.667 0.588542 17 13 12
12.470 48 26082.959 0.588543 19 15 4
12.470 48 4142.475 0.588543 20 11 9
12.470 48 3639.900 0.588543 23 8 3
12.470 48 17676.646 0.588542 24 5 1
12.511 48 236.687 0.130661 17 14 11
12.511 48 4685.903 0.455545 19 14 7
12.511 48 15084.020 0.511461 22 11 1
12.532 48 4970.170 0.384055 20 12 8
12.532 24 2797.412 0.326392 24 4 4
12.552 48 42385.359 0.584411 21 12 5
12.552 24 0.000 0.603295 21 13 0
12.552 24 0.000 0.614612 23 9 0
12.552 48 1401.835 0.584411 24 5 3
12.573 24 0.000 0.628182 16 16 10
12.573 24 0.000 0.574667 18 12 12
12.573 48 27.689 0.583384 20 14 4
12.573 24 0.000 0.593523 22 8 8
12.573 24 974.674 0.583382 24 6 0
12.594 48 7713.880 0.696704 17 15 10
12.594 24 625.189 0.261879 17 17 6
12.594 48 342.001 0.179135 18 13 11
12.594 48 2097.192 0.810675 18 17 1
12.594 48 486.338 1.098598 21 13 2
12.594 48 19185.775 0.653674 22 9 7
12.594 48 257.410 0.137906 22 11 3
12.594 48 21796.854 0.649147 23 7 6
12.594 48 2041.253 0.814043 23 9 2
12.614 48 42259.598 0.660673 18 16 6
12.614 48 1627.606 1.038560 24 6 2
12.635 48 9.570 0.580310 19 16 1
12.635 48 4582.844 0.580308 20 13 7
12.635 48 424.410 0.580309 23 8 5
12.655 48 31903.266 0.579287 18 14 10
12.655 48 5253.120 0.579288 22 10 6
12.675 48 2611.902 0.404740 18 17 3
12.675 48 13236.974 0.497375 19 15 6
12.675 48 1066.550 0.911786 21 10 9
12.716 48 6204.871 0.576235 17 16 9
12.716 48 6512.263 0.576235 19 12 11
12.716 48 7935.641 0.576235 19 16 3
12.716 48 1384.119 0.576235 20 15 1
12.716 48 6295.197 0.576235 21 11 8
12.716 48 24718.215 0.576235 21 13 4
12.716 48 5930.177 0.576235 23 9 4
12.716 24 0.000 0.586224 24 5 5
12.716 48 14408.281 0.576235 24 7 1
12.716 24 0.000 0.599632 25 1 0
12.736 24 120591.719 0.628172 22 12 0
12.736 48 34299.668 0.676419 24 6 4
12.757 48 156.663 0.060188 18 15 9
12.757 48 3943.698 0.739427 19 13 10
12.757 48 10074.289 0.481447 22 11 5
12.757 48 1050.248 0.913322 23 10 1
12.757 48 14211.491 0.658643 25 2 1
12.777 48 118851.406 0.579183 20 14 6
12.777 48 81153.352 0.580444 22 12 2
12.797 48 1636.235 0.572190 20 15 3
12.797 48 5481.227 0.572190 21 12 7
12.797 48 18053.295 0.572190 24 7 3
12.797 24 0.000 0.589848 25 3 0
12.837 48 7405.788 0.689516 18 17 5
12.837 48 2624.350 0.395211 19 14 9
12.837 48 2138.960 0.801308 21 14 1
12.837 48 10.241 1.610380 23 10 3
12.837 48 475.319 1.107143 25 3 2
12.857 24 400540.719 0.601597 24 8 0
12.878 24 0.000 0.594411 17 17 8
12.878 48 0.125 0.568303 19 16 5
12.878 24 0.000 0.668721 20 11 11
12.878 48 9156.891 0.568173 23 8 7
12.878 48 626.749 0.568173 25 4 1
12.898 48 752.168 0.567174 18 16 8
12.898 48 5026.329 0.567173 20 12 10
12.898 48 1582.216 0.567173 22 12 4
12.898 48 10010.276 0.567173 24 8 2
12.918 24 201.271 1.472532 15 15 14
12.918 48 222.695 0.098690 21 13 6
12.918 48 6281.894 0.448796 21 14 3
12.918 24 124.054 0.027292 22 9 9
12.918 48 7534.736 0.458458 23 9 6
12.938 24 6316.977 0.786949 16 14 14
12.938 12 29334.926 0.663537 18 18 0
12.938 48 1147.344 1.143397 22 10 8
12.938 24 3723.319 0.331233 24 6 6
12.958 48 19312.793 0.564184 16 15 13
12.958 48 1304.676 0.564185 19 15 8
12.958 24 0.000 0.629484 19 17 0
12.958 48 7762.704 0.564184 20 13 9
12.958 48 0.036 0.564203 20 15 5
12.958 24 0.000 0.552595 23 11 0
12.958 48 1174.496 0.564185 24 7 5
12.958 48 1864.236 0.564185 25 4 3
12.958 24 0.000 0.610466 25 5 0
12.977 24 0.000 1.007893 18 18 2
12.997 48 1257.012 0.317586 17 14 13
12.997 48 803.021 0.965953 19 17 2
12.997 48 3526.509 0.407812 22 11 7
12.997 48 33105.316 0.617943 22 13 1
12.997 48 3190.903 0.400563 23 10 5
12.997 48 2609.250 0.385052 23 11 2
12.997 48 628.391 0.236561 25 5 2
13.017 24 57828.707 0.570060 16 16 12
13.017 24 20781.746 0.576011 20 16 0
13.017 48 89788.805 0.568307 24 8 4
13.037 48 46834.953 0.560224 17 15 12
13.037 48 4458.039 0.560224 24 9 1
13.057 48 64854.797 0.633972 20 14 8
13.057 48 1893.262 1.063184 20 16 2
13.077 24 11205.765 0.468121 18 13 13
13.077 48 2285.663 0.369342 18 17 7
13.077 48 339.441 0.148501 21 11 10
13.077 48 8367.582 0.454671 21 14 5
13.077 48 7108.424 0.446360 22 13 3
13.077 48 21281.016 0.492060 25 6 1
13.096 48 68508.367 0.565478 18 14 12
13.096 24 155737.828 0.562708 18 18 4
13.096 48 103275.125 0.563951 22 12 6
13.116 48 11.861 0.556303 17 16 11
13.116 48 9.261 0.556294 19 16 7
13.116 48 858.044 0.556291 19 17 4
13.116 48 3173.471 0.556291 21 12 9
13.116 24 0.000 0.612472 21 15 0
13.116 48 3159.007 0.556291 23 11 4
13.116 48 11758.186 0.556291 24 9 3
13.116 48 817.126 0.556292 25 5 4
13.155 48 2785.415 0.760410 18 15 11
13.155 48 12620.544 0.468744 21 15 2
13.155 48 9819.387 0.660046 25 6 3
13.175 48 16619.480 0.433205 20 16 4
13.195 48 1551.486 0.552387 19 13 12
13.195 48 2451.736 0.552386 20 15 7
13.195 48 8249.820 0.552386 21 13 8
13.195 48 4248.809 0.552386 23 9 8
13.195 48 21355.518 0.552386 23 12 1
13.195 24 0.000 0.511862 24 7 7
13.195 24 0.000 0.532794 25 7 0
13.214 48 45.962 0.551413 24 8 6
13.214 24 1375.366 0.551414 24 10 0
13.214 6 0.000 0.581809 26 0 0
13.234 24 125.795 1.815142 17 17 10
13.234 48 3624.754 0.730497 19 14 11
13.234 48 20871.211 0.622388 22 13 5
13.234 48 16.691 0.846322 23 10 7
13.234 48 1982.095 0.799991 25 7 2
13.234 24 4838.143 0.704812 26 1 1
13.253 48 6534.444 0.364097 18 16 10
13.253 24 2222.547 0.955805 22 14 0
13.253 48 24731.307 0.449439 24 10 2
13.253 24 3270.319 0.298545 26 2 0
13.273 48 7338.060 0.548509 21 15 4
13.273 48 4058.442 0.548509 23 12 3
13.273 48 14392.250 0.548509 24 9 5
13.292 24 0.000 0.517475 18 18 6
13.292 24 0.000 0.554471 22 10 10
13.292 48 5628.166 0.547544 22 14 2
13.292 24 0.000 0.523200 26 2 2
13.311 48 727.746 0.986569 19 15 10
13.311 48 5154.896 0.696772 19 17 6
13.311 48 1101.667 0.894374 19 18 1
13.311 48 81.138 2.249579 21 14 7
13.311 48 2995.469 0.747331 22 11 9
13.311 48 1254.767 0.869951 23 11 6
13.311 48 2509.770 0.767480 25 6 5
13.311 48 2730.918 0.757605 26 3 1
13.331 24 43466.543 0.728675 20 12 12
13.350 48 216.962 0.544659 20 13 11
13.350 48 1647.968 0.544658 20 17 1
13.350 48 4276.428 0.544658 25 7 4
13.350 48 15238.885 0.544658 25 8 1
13.370 48 84474.062 0.610509 20 16 6
13.370 48 74722.828 0.614865 22 12 8
13.370 48 40800.441 0.641053 24 10 4
13.370 24 179990.234 0.589052 26 4 0
13.389 48 404.751 0.154578 18 17 9
13.389 48 4652.005 0.403703 19 18 3
13.389 24 391.049 1.179888 26 3 3
13.408 48 232808.609 0.546398 20 14 10
13.408 48 99927.586 0.548833 22 14 4
13.408 48 188579.734 0.546911 26 4 2
13.427 48 641.491 0.540835 19 16 9
13.427 48 3573.906 0.540835 20 17 3
13.427 48 1412.832 0.540835 21 16 1
13.427 48 2392.580 0.540836 23 12 5
13.427 24 0.000 0.583230 23 13 0
13.427 48 13539.380 0.540835 24 11 1
13.427 48 4325.955 0.540835 25 8 3
13.466 48 24231.410 0.475178 21 15 6
13.466 48 2542.877 0.355048 22 13 7
13.466 48 9143.574 0.651304 23 13 2
13.466 48 44289.496 0.491389 26 5 1
13.485 24 386159.906 0.572714 24 8 8
13.504 24 0.000 0.537275 16 15 15
13.504 48 79.871 0.537040 20 15 9
13.504 48 17695.494 0.537038 21 12 11
13.504 48 895.248 0.537037 21 16 3
13.504 48 12284.010 0.537038 24 9 7
13.504 48 2309.103 0.537038 24 11 3
13.504 24 0.000 0.591044 25 9 0
13.523 24 0.000 0.402816 16 16 14
13.523 24 0.000 0.504490 26 4 4
13.542 48 9775.623 0.436148 17 15 14
13.542 48 7637.786 0.659352 19 18 5
13.542 48 804.636 0.959318 21 13 10
13.542 48 3708.660 0.379537 22 15 1
13.542 48 7600.334 0.423640 23 10 9
13.542 48 924.706 0.926972 25 7 6
13.542 48 1905.093 0.325385 25 9 2
13.542 48 2232.918 0.339564 26 5 3
13.561 24 21222.117 0.656080 18 18 8
13.561 48 7232.998 0.354522 24 10 6
13.561 24 9536.754 0.720598 26 6 0
13.580 48 7.376 0.533264 17 16 13
13.580 48 8833.213 0.533268 19 17 8
13.580 48 507.977 0.533268 20 17 5
13.580 48 4764.306 0.533268 23 11 8
13.580 48 16286.252 0.533268 23 13 4
13.580 48 13357.467 0.533268 25 8 5
13.599 24 0.000 0.552398 18 14 14
13.599 48 240.404 0.532327 22 14 6
13.599 48 142.940 0.532329 26 6 2
13.618 48 1512.714 0.298800 18 15 13
13.618 48 13081.259 0.444836 21 14 9
13.618 48 5860.847 0.404917 22 15 3
13.637 48 68102.938 0.539222 20 16 8
13.637 24 161144.547 0.536148 24 12 0
13.656 24 0.000 0.567319 17 17 12
13.656 12 0.000 0.641853 19 19 0
13.656 48 1567.160 0.529524 21 16 5
13.656 48 1346.462 0.529525 23 12 7
13.656 48 21105.100 0.529525 24 11 5
13.656 48 24573.625 0.529525 25 9 4
13.675 48 55987.738 0.613047 18 16 12
13.675 24 294443.688 0.564651 20 18 0
13.675 48 83364.617 0.597342 24 12 2
13.694 48 21622.609 0.459374 19 14 13
13.694 24 833.602 0.229171 19 19 2
13.694 24 41057.730 0.580182 22 11 11
13.694 48 258.246 0.079186 23 14 1
13.694 48 12107.555 0.626306 25 10 1
13.694 24 5123.758 0.392501 26 5 5
13.694 48 1369.571 0.284268 26 7 1
13.713 48 123013.125 0.533304 20 18 2
13.713 48 71978.680 0.535331 22 12 10
13.713 48 60597.105 0.536107 26 6 4
13.732 48 36370.348 0.525808 19 15 12
13.732 48 6984.553 0.525808 21 15 8
13.732 24 0.000 0.556326 21 17 0
13.732 24 0.000 0.520802 27 1 0
13.769 48 805.905 0.956158 18 17 11
13.769 48 3654.744 0.365521 19 18 7
13.769 48 11134.515 0.429707 21 17 2
13.769 48 514.928 1.084913 22 13 9
13.769 48 1388.901 0.841474 22 15 5
13.769 48 5090.758 0.680468 23 13 6
13.769 48 118.760 1.919536 23 14 3
13.769 48 2065.000 0.778836 25 10 3
13.769 48 1673.389 0.810194 26 7 3
13.769 48 1752.430 0.304277 27 2 1
13.788 48 115.481 0.140374 24 12 4
13.807 48 13001.454 0.522117 19 16 11
13.807 24 0.000 0.550671 19 19 4
13.807 24 0.000 0.644168 20 13 13
13.807 48 3757.070 0.522116 20 17 7
13.807 24 0.000 0.619965 24 9 9
13.807 48 166.534 0.522118 25 8 7
13.807 24 0.000 0.642318 27 3 0
13.825 48 432.177 0.521196 20 14 12
13.825 48 12337.985 0.521198 20 18 4
13.825 24 6312.746 0.521198 22 16 0
13.825 48 1634.139 0.521198 24 10 8
13.825 24 6345.703 0.521198 26 8 0
13.844 48 1706.373 0.805092 25 9 6
13.844 48 3221.206 0.721440 27 3 2
13.863 48 3274.244 0.263069 22 14 8
13.863 48 3664.048 0.274854 22 16 2
13.863 48 4512.130 0.805788 26 8 2
13.881 48 202.873 0.518453 20 15 11
13.881 48 53.521 0.518447 21 16 7
13.881 48 164.154 0.518448 21 17 4
13.881 48 4790.123 0.518451 24 11 7
13.881 48 143.928 0.518449 24 13 1
13.881 24 0.000 0.551518 25 11 0
13.881 48 5639.001 0.518451 27 4 1
13.900 24 0.000 0.480472 18 18 10
13.900 24 0.000 0.536103 26 6 6
13.919 48 4715.988 0.374636 19 17 10
13.919 48 1152.720 0.252972 23 11 10
13.919 48 14412.105 0.608090 23 14 5
13.919 48 8711.256 0.409941 25 10 5
13.919 48 12031.894 0.425067 25 11 2
13.919 48 997.063 0.904037 26 7 5
13.956 48 1712.422 0.514812 21 13 12
13.956 48 51.106 0.514808 23 12 9
13.956 24 0.000 0.493669 23 15 0
13.956 48 3151.326 0.514811 24 13 3
13.956 48 4966.151 0.514812 27 4 3
13.956 24 0.000 0.544797 27 5 0
13.974 48 184851.469 0.560507 20 16 10
13.974 48 70229.789 0.590530 22 16 4
13.974 48 102809.844 0.576854 24 12 6
13.974 48 77585.695 0.586684 26 8 4
13.993 24 27149.818 0.579212 19 19 6
13.993 48 2431.204 0.749978 21 14 11
13.993 48 34488.859 0.457621 22 15 7
13.993 48 364.326 0.108468 23 15 2
13.993 48 1418.150 0.270548 26 9 1
13.993 48 9578.805 0.410609 27 5 2
14.011 48 46500.246 0.523157 20 18 6
14.030 48 0.084 0.511293 20 19 1
14.030 48 3199.971 0.511198 23 13 8
14.030 48 13297.949 0.511198 25 11 4
-------------- next part --------------
# 0.600000 0.003500 0.050000 0.000100 # sample: packing radius height wall
# 1.472205 0.153000 0.100000 0.053205 3.000000 # PSD: radius height pitch gap pressure
1889.494 133.139 19.677 # YIG powder: Vc sigma_abs sigma_inc
1111 YIG
0.719 12 0.000 0.998322 1 1 0
1.016 6 0.000 0.996636 2 0 0
1.245 24 203.580 0.984168 2 1 1
1.438 12 3016.019 0.986776 2 2 0
1.607 24 0.000 0.992151 3 1 0
1.761 8 0.000 0.989722 2 2 2
1.902 48 1506.029 0.978739 3 2 1
2.033 6 98281.828 0.986199 4 0 0
2.156 12 0.000 0.986071 3 3 0
2.156 24 0.000 0.981520 4 1 1
2.273 24 54691.691 0.986475 4 2 0
2.384 24 36367.449 0.977947 3 3 2
2.490 24 18260.684 0.979746 4 2 2
2.592 48 13085.104 0.977365 4 3 1
2.592 24 0.000 0.980939 5 1 0
2.784 48 9884.264 0.981401 5 2 1
2.875 12 11381.769 0.984795 4 4 0
2.963 24 0.000 0.972488 4 3 3
2.963 24 0.000 0.974189 5 3 0
3.049 24 0.000 0.970052 4 4 2
3.049 6 0.000 0.980046 6 0 0
3.133 48 4738.101 0.980699 5 3 2
3.133 24 26526.645 0.972831 6 1 1
3.214 24 17280.061 0.952795 6 2 0
3.294 48 1072.821 0.963691 5 4 1
3.371 24 0.000 0.964791 6 2 2
3.447 48 5640.986 0.945413 6 3 1
3.521 8 108324.977 0.971225 4 4 4
3.594 48 710.790 0.956927 5 4 3
3.594 12 0.000 0.972702 5 5 0
3.594 24 0.000 0.964389 7 1 0
3.665 24 190097.016 0.960594 6 4 0
3.735 24 8854.975 0.967531 5 5 2
3.735 24 2605.230 0.928074 6 3 3
3.735 48 136.468 0.844741 7 2 1
3.803 48 127188.094 0.952464 6 4 2
3.871 24 0.000 0.959452 7 3 0
4.002 48 48.465 1.173685 6 5 1
4.002 48 14212.914 0.934375 7 3 2
4.066 6 244666.688 0.951433 8 0 0
4.129 24 0.000 0.941082 5 5 4
4.129 48 7215.575 0.943539 7 4 1
4.129 24 0.000 0.937901 8 1 1
4.191 24 0.000 0.946819 6 4 4
4.191 24 17019.494 0.941878 8 2 0
4.252 48 8083.322 0.921666 6 5 3
4.313 12 734.671 1.049026 6 6 0
4.313 24 10628.211 0.910545 8 2 2
4.372 48 7646.340 0.936915 7 4 3
4.372 24 0.000 0.915093 7 5 0
4.372 48 776.047 0.936915 8 3 1
4.431 24 0.000 0.930654 6 6 2
4.489 48 8357.111 0.913423 7 5 2
4.546 24 99628.523 0.932912 8 4 0
4.602 24 0.000 0.925655 8 3 3
4.602 24 0.000 0.927404 9 1 0
4.658 48 41391.684 0.946833 8 4 2
4.713 24 6454.844 0.952557 6 5 5
4.713 48 20689.980 0.912980 7 6 1
4.713 48 1554.205 0.876190 9 2 1
4.768 24 82241.906 0.926566 6 6 4
4.822 48 21384.988 0.923807 7 5 4
4.822 48 15318.223 0.923807 8 5 1
4.822 24 0.000 0.925684 9 3 0
4.928 48 504.394 0.824814 7 6 3
4.928 48 19.874 1.482843 9 3 2
4.980 24 16359.466 0.889442 8 4 4
5.031 12 0.000 0.922682 7 7 0
5.031 48 39592.711 0.917322 8 5 3
5.031 48 25325.520 0.917322 9 4 1
5.082 24 129.762 0.915708 8 6 0
5.082 6 0.000 0.924715 10 0 0
5.133 24 1579.174 0.975029 7 7 2
5.133 24 20095.832 0.897369 10 1 1
5.183 48 673.349 1.077036 8 6 2
5.183 24 43752.160 0.932161 10 2 0
5.233 48 10300.698 0.910883 9 4 3
5.233 24 0.000 0.893915 9 5 0
5.282 8 0.000 0.903652 6 6 6
5.282 24 0.000 0.903433 10 2 2
5.330 48 11530.723 0.931612 7 6 5
5.330 48 131.671 0.698822 9 5 2
5.330 48 1564.328 0.844278 10 3 1
5.426 24 0.000 0.903811 7 7 4
5.426 24 0.000 0.911782 8 5 5
5.426 48 5560.977 0.904489 8 7 1
5.474 48 43032.797 0.926892 8 6 4
5.474 24 121046.172 0.917166 10 4 0
5.521 48 869.710 0.995956 9 6 1
5.521 24 1071.323 0.986376 10 3 3
5.567 48 100115.828 0.901113 10 4 2
5.614 48 1334.206 0.898139 8 7 3
5.614 48 1100.267 0.898139 9 5 4
5.614 24 0.000 0.894209 11 1 0
5.705 48 44.953 0.516079 9 6 3
5.705 48 11047.518 0.922668 10 5 1
5.705 48 4519.932 0.938451 11 2 1
5.750 12 231623.359 0.905655 8 8 0
5.795 24 0.000 0.898561 9 7 0
5.795 24 0.000 0.895735 11 3 0
5.839 24 0.000 0.902606 8 8 2
5.839 24 0.000 0.884242 10 4 4
5.883 24 24914.727 0.869469 7 7 6
5.883 48 7679.861 0.854214 9 7 2
5.883 48 125.314 1.182243 10 5 3
5.883 48 1178.122 0.802006 11 3 2
5.927 24 7067.095 0.950199 8 6 6
5.927 24 11585.631 0.936204 10 6 0
5.970 48 462.440 0.885574 8 7 5
5.970 48 5778.562 0.885574 11 4 1
6.014 48 22267.492 0.884016 10 6 2
6.056 48 486.486 1.034189 9 6 5
6.099 24 72231.336 0.882867 8 8 4
6.099 6 19129.297 0.884716 12 0 0
6.141 48 72.866 0.879357 9 7 4
6.141 48 27934.947 0.879357 9 8 1
6.141 48 5137.507 0.879357 11 4 3
6.141 24 0.000 0.883748 11 5 0
6.141 24 0.000 0.887918 12 1 1
6.183 24 66256.578 0.901881 12 2 0
6.225 24 1145.936 0.779438 10 5 5
6.225 48 2368.012 0.946966 10 7 1
6.225 48 2.107 0.093289 11 5 2
6.266 48 84037.055 0.876638 10 6 4
6.266 24 134249.734 0.876239 12 2 2
6.307 48 1321.134 0.873184 9 8 3
6.307 48 23028.475 0.873184 12 3 1
6.388 48 378.576 0.698084 10 7 3
6.388 48 1202.565 0.771261 11 6 1
6.429 24 4163.282 0.777967 12 4 0
6.469 24 0.000 0.847652 8 7 7
6.469 12 0.000 0.913683 9 9 0
6.469 48 1046.849 0.867055 11 5 4
6.469 24 0.000 0.732765 12 3 3
6.509 24 0.000 0.861072 8 8 6
6.509 24 10020.346 0.865529 10 8 0
6.509 48 100.722 0.865529 12 4 2
6.548 48 5719.824 0.913443 9 7 6
6.548 24 10481.625 0.900393 9 9 2
6.548 48 34038.199 0.884106 11 6 3
6.587 48 2212.532 0.734415 10 8 2
6.627 48 8340.317 0.860968 9 8 5
6.627 24 0.000 0.864604 11 7 0
6.627 48 1144.634 0.860968 12 5 1
6.627 24 0.000 0.866319 13 1 0
6.665 24 0.000 0.861820 10 6 6
6.704 48 4163.081 0.799453 10 7 5
6.704 48 2879.308 0.931003 11 7 2
6.704 48 18606.156 0.830018 13 2 1
6.743 24 31417.455 0.935357 12 4 4
6.781 24 0.000 0.851566 9 9 4
6.781 48 19656.873 0.854924 12 5 3
6.781 24 0.000 0.899097 13 3 0
6.819 48 154874.438 0.872072 10 8 4
6.819 24 45146.090 0.888126 12 6 0
6.856 48 1.584 6.976236 10 9 1
6.856 48 51.244 1.500796 11 6 5
6.856 48 1127.712 0.974933 13 3 2
6.894 48 50511.836 0.853395 12 6 2
6.931 48 34190.152 0.848923 11 7 4
6.931 48 2164.648 0.848923 11 8 1
6.931 48 524.398 0.848923 13 4 1
7.006 48 10537.984 0.806079 10 9 3
7.042 8 279514.156 0.860435 8 8 8
7.079 48 4450.304 0.842964 9 8 7
7.079 48 91.726 0.842964 11 8 3
7.079 24 0.000 0.833867 12 5 5
7.079 48 236.346 0.842964 12 7 1
7.079 48 13131.040 0.842964 13 4 3
7.079 24 0.000 0.818209 13 5 0
7.115 48 2832.459 0.841480 12 6 4
7.115 6 0.000 0.853769 14 0 0
7.152 24 55446.820 0.821886 9 9 6
7.152 24 10563.778 0.798793 10 7 7
7.152 48 29967.029 0.815409 13 5 2
7.152 24 871.975 0.991503 14 1 1
7.188 48 2430.303 0.697682 10 8 6
7.188 12 61323.203 0.809452 10 10 0
7.188 24 15699.059 0.781564 14 2 0
7.223 24 0.000 0.866574 11 9 0
7.223 48 26204.113 0.837046 12 7 3
7.259 24 0.000 0.837881 10 10 2
7.259 24 0.000 0.829693 14 2 2
7.295 48 11099.855 0.792550 10 9 5
7.295 48 1751.051 0.731513 11 7 6
7.295 48 1620.526 0.947878 11 9 2
7.295 48 5809.077 0.776945 13 6 1
7.295 48 3604.480 0.909587 14 3 1
7.330 24 73445.188 0.835414 12 8 0
7.365 48 15921.470 0.831170 11 8 5
7.365 48 5873.266 0.831170 13 5 4
7.400 48 76166.891 0.860376 12 8 2
7.400 24 54682.270 0.865962 14 4 0
7.435 48 19258.697 0.795594 13 6 3
7.435 24 2714.380 0.742726 14 3 3
7.470 24 165646.812 0.828709 10 10 4
7.470 24 74290.812 0.829656 12 6 6
7.470 48 151288.562 0.828798 14 4 2
7.504 48 10344.729 0.825335 11 9 4
7.504 48 5411.900 0.825336 12 7 5
7.504 24 0.000 0.822316 13 7 0
7.573 48 1254.346 0.694604 11 10 1
7.573 48 148.693 0.481002 13 7 2
7.573 48 30448.402 0.795615 14 5 1
7.607 48 1966.262 1.010599 12 8 4
7.640 24 0.000 0.827541 9 9 8
7.640 48 6363.631 0.819542 12 9 1
7.640 24 0.000 0.837935 15 1 0
7.674 24 0.000 0.824846 10 8 8
7.674 24 0.000 0.818697 14 4 4
7.708 48 3849.914 0.896828 10 9 7
7.708 48 1094.647 0.676433 11 10 3
7.708 48 79.345 0.362349 13 6 5
7.708 48 79.048 1.454426 14 5 3
7.708 48 5521.893 0.883341 15 2 1
7.741 24 74.444 0.141779 14 6 0
7.775 48 15044.647 0.813789 11 8 7
7.775 48 11309.317 0.813789 12 9 3
7.775 48 14739.329 0.813789 13 7 4
7.775 48 3247.167 0.813789 13 8 1
7.775 24 0.000 0.828548 15 3 0
7.808 24 0.000 0.844919 10 10 6
7.808 48 11997.312 0.812357 14 6 2
7.841 48 19856.273 0.775865 11 9 6
7.841 48 13406.674 0.854643 15 3 2
7.906 12 0.000 0.816166 11 11 0
7.906 24 0.000 0.830149 12 7 7
7.906 48 7061.392 0.808076 13 8 3
7.906 48 2060.587 0.808076 15 4 1
7.939 48 93269.773 0.837780 12 8 6
7.939 24 47890.504 0.850255 12 10 0
7.971 48 7796.717 0.748176 11 10 5
7.971 24 7086.558 0.745438 11 11 2
7.971 24 2328.366 0.702413 14 5 5
7.971 48 19.962 0.066294 14 7 1
8.004 48 112340.125 0.806478 12 10 2
8.004 48 53296.637 0.807681 14 6 4
8.036 48 14131.926 0.802403 12 9 5
8.036 24 0.000 0.829866 13 9 0
8.036 48 956.852 0.802403 15 4 3
8.036 24 0.000 0.785195 15 5 0
8.100 48 125.046 1.339667 13 7 6
8.100 48 2835.213 0.703768 13 9 2
8.100 48 10.602 0.000098 14 7 3
8.100 48 1178.623 0.653589 15 5 2
8.132 6 151833.891 0.825886 16 0 0
8.163 24 0.000 0.795598 11 11 4
8.163 48 1550.868 0.796771 13 8 5
8.163 24 0.000 0.872871 16 1 1
8.195 48 242.107 0.795369 12 10 4
8.195 24 33236.344 0.795369 14 8 0
8.195 24 34660.660 0.795369 16 2 0
8.227 24 12606.899 0.842757 10 9 9
8.227 48 2381.923 0.687250 15 6 1
8.258 24 864.789 0.513735 10 10 8
8.258 48 6570.184 0.908619 14 8 2
8.258 24 6778.777 0.906760 16 2 2
8.289 48 966.176 0.791177 11 9 8
8.289 48 43581.391 0.791178 12 11 1
8.289 48 7196.764 0.791178 13 9 4
8.289 48 10670.142 0.791177 15 5 4
8.289 48 3204.440 0.791177 16 3 1
8.320 24 0.000 0.808097 14 6 6
8.351 48 200.804 1.226601 11 10 7
8.351 48 2200.392 0.910606 13 10 1
8.351 48 1091.037 0.629967 14 7 5
8.351 48 3256.225 0.888213 15 6 3
8.382 24 93462.555 0.790193 12 8 8
8.382 24 75682.148 0.790548 16 4 0
8.413 48 11445.481 0.785624 12 9 7
8.413 48 142.793 0.785624 12 11 3
8.413 24 0.000 0.807413 15 7 0
8.413 24 0.000 0.793343 16 3 3
8.443 48 66963.406 0.824903 14 8 4
8.443 48 152509.391 0.811070 16 4 2
8.474 24 991.613 0.613598 11 11 6
8.474 48 1375.104 0.637806 13 10 3
8.474 48 3873.620 0.694648 14 9 1
8.474 48 6604.764 0.853982 15 7 2
8.504 48 105229.992 0.784568 12 10 6
8.535 48 20650.426 0.780109 13 8 7
8.535 48 11068.548 0.780109 16 5 1
8.595 48 1777.985 0.920500 13 9 6
8.595 48 9026.879 0.839398 14 9 3
8.595 48 1867.975 0.649371 15 6 5
8.625 12 440.601 0.382222 12 12 0
8.625 24 4123.180 0.934647 16 4 4
8.655 48 2586.402 0.774632 12 11 5
8.655 24 0.000 0.773600 13 11 0
8.655 48 3.233 0.774635 15 7 4
8.655 48 53071.383 0.774632 15 8 1
8.655 48 3789.377 0.774632 16 5 3
8.655 24 0.000 0.795470 17 1 0
8.685 24 0.000 0.771895 12 12 2
8.685 24 3643.341 0.773269 16 6 0
8.714 48 941.143 0.591423 13 10 5
8.714 48 3285.156 0.878372 13 11 2
8.714 24 39280.148 0.801991 14 7 7
8.714 48 8748.503 0.710066 17 2 1
8.744 48 1214.614 0.512271 14 8 6
8.744 24 2282.333 0.576950 14 10 0
8.744 48 57.059 2.642583 16 6 2
8.773 48 308.901 0.769195 15 8 3
8.773 24 0.000 0.819555 17 3 0
8.803 8 0.000 0.755858 10 10 10
8.803 48 454.897 0.767840 14 10 2
8.832 48 12771.266 0.820756 11 10 9
8.832 48 6171.375 0.845140 14 9 5
8.832 48 590.266 1.034670 17 3 2
8.861 24 104101.289 0.832626 12 12 4
8.890 24 0.000 0.770600 11 11 8
8.890 24 0.000 0.675269 12 9 9
8.890 48 9594.028 0.763795 13 11 4
8.890 24 0.000 0.801700 15 9 0
8.890 24 0.000 0.763028 16 5 5
8.890 48 11653.093 0.763795 16 7 1
8.890 48 960.200 0.763795 17 4 1
8.919 48 41982.527 0.818643 12 10 8
8.919 48 84537.617 0.801842 16 6 4
8.948 48 1440.681 0.931723 15 7 6
8.948 48 6684.949 0.687764 15 9 2
8.977 48 31867.904 0.765981 14 10 4
9.006 48 3182.443 0.758433 12 11 7
9.006 48 58.134 0.758433 13 9 8
9.006 48 11570.396 0.758433 13 12 1
9.006 48 5064.502 0.758433 15 8 5
9.006 48 6890.246 0.758433 16 7 3
9.006 48 12437.841 0.758433 17 4 3
9.006 24 0.000 0.781995 17 5 0
9.063 48 2402.020 0.633288 13 10 7
9.063 48 935.103 0.974633 14 11 1
9.063 48 2552.169 0.636789 17 5 2
9.092 24 489417.750 0.772806 16 8 0
9.120 48 16024.504 0.753109 13 12 3
9.120 48 2811.926 0.753109 15 9 4
9.148 24 0.000 0.772726 12 12 6
9.148 24 0.000 0.725124 14 8 8
9.148 48 1665.530 0.751784 16 8 2
9.148 6 0.000 0.769562 18 0 0
9.176 48 13771.264 0.697036 13 11 6
9.176 48 29490.500 0.713739 14 9 7
9.176 48 15184.293 0.699537 14 11 3
9.176 48 1268.854 0.940422 15 10 1
9.176 48 229.789 0.388343 17 6 1
9.176 24 5101.780 0.663727 18 1 1
9.205 24 26858.607 0.685284 16 6 6
9.205 24 1676.267 1.033228 18 2 0
9.233 48 4989.955 0.747822 16 7 5
9.233 48 25568.219 0.747822 17 5 4
9.260 48 15964.095 0.746507 14 10 6
9.260 24 0.000 0.738002 18 2 2
9.288 48 4842.825 0.841539 15 10 3
9.288 48 59.952 1.821303 17 6 3
9.288 48 842.021 0.986068 18 3 1
9.316 48 73451.273 0.748262 16 8 4
9.344 48 1779.935 0.742573 13 12 5
9.344 12 0.000 0.750385 13 13 0
9.344 48 7446.031 0.742573 15 8 7
9.344 48 64.178 0.742571 16 9 1
9.344 24 0.000 0.721414 17 7 0
9.371 24 116010.547 0.777449 14 12 0
9.371 24 84.167 2.662763 18 4 0
9.399 24 16458.018 0.792424 11 11 10
9.399 24 25.086 2.650020 13 13 2
9.399 48 4530.413 0.841504 14 11 5
9.399 48 16069.531 0.793065 15 9 6
9.399 48 967.204 0.542545 17 7 2
9.399 24 13988.126 0.796950 18 3 3
9.426 24 37421.613 0.744935 12 10 10
9.426 48 39291.531 0.744784 14 12 2
9.426 48 52253.723 0.743969 18 4 2
9.454 48 8956.533 0.737360 12 11 9
9.454 24 0.000 0.769230 15 11 0
9.454 48 6119.978 0.737360 16 9 3
9.508 48 367.693 0.425812 13 10 9
9.508 48 205.257 1.279353 15 10 5
9.508 48 18085.365 0.685595 15 11 2
9.508 48 47413.023 0.704197 17 6 5
9.508 48 11488.113 0.673343 18 5 1
9.535 24 2087.741 0.507875 12 12 8
9.562 48 8.731 0.732185 13 11 8
9.562 24 0.000 0.730241 13 13 4
9.562 24 0.000 0.779491 16 7 7
9.562 48 7435.995 0.732184 17 7 4
9.562 48 4377.982 0.732184 17 8 1
9.589 48 21.935 0.730896 14 12 4
9.589 48 592.034 0.730895 16 8 6
9.589 24 12893.613 0.730895 16 10 0
9.589 24 0.000 0.733011 18 4 4
9.616 24 394.922 0.426086 14 9 9
9.616 48 1802.463 0.577434 18 5 3
9.643 48 986.512 1.137751 14 10 8
9.643 48 4248.782 0.563507 16 10 2
9.643 24 10720.993 0.622106 18 6 0
9.670 48 30420.504 0.727044 13 12 7
9.670 48 6393.077 0.727044 15 11 4
9.670 48 3875.861 0.727044 16 9 5
9.670 48 9838.925 0.727044 17 8 3
9.670 24 0.000 0.748096 19 1 0
9.697 48 8127.695 0.725765 18 6 2
9.723 48 770.210 0.495789 14 11 7
9.723 48 2651.060 0.595845 14 13 1
9.723 48 16128.882 0.670819 19 2 1
9.776 48 3330.684 0.721940 15 9 8
9.776 48 723.887 0.721940 15 12 1
9.776 24 0.000 0.759183 17 9 0
9.776 24 0.000 0.781093 19 3 0
9.802 48 6867.757 0.886021 16 10 4
9.829 24 35440.613 0.757362 13 13 6
9.829 48 10398.315 0.790245 14 13 3
9.829 48 6696.791 0.635780 15 10 7
9.829 48 9738.723 0.792662 17 7 6
9.829 48 285.456 1.197523 17 9 2
9.829 24 14787.541 0.662560 18 5 5
9.829 48 926.501 0.969829 18 7 1
9.829 48 3458.568 0.844357 19 3 2
9.855 48 101912.961 0.722259 14 12 6
9.855 48 70842.172 0.723083 18 6 4
9.881 48 4636.409 0.716873 15 12 3
9.881 48 871.953 0.716872 16 11 1
9.881 48 2569.428 0.716873 17 8 5
9.881 48 64.828 0.716872 19 4 1
9.933 48 2296.289 0.871588 15 11 6
9.933 48 3882.000 0.833895 18 7 3
9.959 24 429715.625 0.735588 16 8 8
9.985 24 0.000 0.728021 12 11 11
9.985 48 3829.123 0.711840 16 9 7
9.985 48 2300.313 0.711840 16 11 3
9.985 48 2157.346 0.711840 17 9 4
9.985 48 5.185 0.711839 19 4 3
9.985 24 0.000 0.731896 19 5 0
10.011 24 0.000 0.729847 12 12 10
10.011 24 8407.514 0.710588 18 8 0
10.037 48 4403.742 0.604092 13 11 10
10.037 48 0.061 333.109955 14 13 5
10.037 48 5228.630 0.612431 17 10 1
10.037 48 10274.168 0.639479 19 5 2
10.063 12 11096.196 0.827805 14 14 0
10.063 48 491.430 0.266939 16 10 6
10.063 48 354.569 0.210751 18 8 2
10.088 48 23376.588 0.706843 13 12 9
10.088 48 33480.711 0.706843 15 12 5
10.088 24 0.000 0.738399 15 13 0
10.114 24 0.000 0.727591 14 10 10
10.114 24 0.000 0.748768 14 14 2
10.114 24 0.000 0.757605 18 6 6
10.139 48 7466.300 0.621679 14 11 9
10.139 48 990.249 0.957920 15 13 2
10.139 48 1190.770 0.934025 17 10 3
10.139 48 9.645 0.329466 18 7 5
10.139 48 8873.878 0.784734 19 6 1
10.165 24 135888.438 0.706895 16 12 0
10.165 6 200452.516 0.706227 20 0 0
10.190 24 0.000 0.702737 13 13 8
10.190 48 12294.665 0.701881 16 11 5
10.190 48 5908.372 0.701881 17 8 7
10.190 48 3440.059 0.701881 19 5 4
10.190 24 0.000 0.664889 20 1 1
10.215 48 103570.594 0.744078 14 12 8
10.215 48 78204.258 0.750741 16 12 2
10.215 48 108975.664 0.742971 18 8 4
10.215 24 200123.938 0.731759 20 2 0
10.241 48 16443.350 0.642387 15 10 9
10.241 48 821.729 0.463152 17 9 6
10.241 48 183.235 0.256509 18 9 1
10.241 48 11.070 4.744914 19 6 3
10.266 24 85986.914 0.703017 14 14 4
10.266 24 176294.906 0.701557 20 2 2
10.291 48 1421.877 0.696954 15 11 8
10.291 48 10078.014 0.696954 15 13 4
10.291 24 0.000 0.705679 17 11 0
10.291 24 0.000 0.727756 19 7 0
10.291 48 7.288 0.696954 20 3 1
10.341 48 211.520 1.312292 14 13 7
10.341 48 1385.313 0.912704 17 10 5
10.341 48 16359.064 0.636597 17 11 2
10.341 48 490.017 0.394712 18 9 3
10.341 48 6772.580 0.789503 19 7 2
10.366 48 1908.914 0.434801 16 12 4
10.366 24 23541.840 0.613565 20 4 0
10.391 48 5924.918 0.692061 15 12 7
10.391 24 0.000 0.799210 16 9 9
10.391 24 0.000 0.729334 20 3 3
10.416 48 6394.384 0.690844 16 10 8
10.416 48 3046.540 0.690844 20 4 2
10.440 48 2461.563 0.543780 15 14 1
10.440 24 135.264 1.506200 18 7 7
10.440 48 19656.266 0.636037 19 6 5
10.465 48 22722.738 0.606361 18 8 6
10.465 24 17729.400 0.595911 18 10 0
10.490 48 5276.253 0.687204 16 11 7
10.490 48 448.920 0.687203 16 13 1
10.490 48 8295.166 0.687203 17 11 4
10.490 48 139.147 0.687203 19 7 4
10.490 48 0.328 0.687215 19 8 1
10.490 48 3527.839 0.687203 20 5 1
10.514 24 0.000 0.727619 14 14 6
10.514 48 48499.012 0.685994 18 10 2
10.539 48 647.314 0.413689 15 13 6
10.539 48 26527.465 0.637960 15 14 3
10.539 48 487.215 0.378284 18 9 5
10.563 8 52768.902 0.807684 12 12 12
10.563 24 92219.422 0.776506 20 4 4
10.588 48 5757.705 0.682379 13 12 11
10.588 48 9340.360 0.682379 16 13 3
10.588 48 710.295 0.682379 17 9 8
10.588 48 10436.022 0.682379 17 12 1
10.588 48 27884.703 0.682379 19 8 3
10.588 48 5.528 0.682372 20 5 3
10.612 48 62078.328 0.740588 16 12 6
10.612 24 2835.733 0.980520 20 6 0
10.637 24 336.916 1.175382 13 13 10
10.637 24 33559.977 0.637768 14 11 11
10.637 48 930.464 0.961942 17 10 7
10.661 48 98878.328 0.683607 14 12 10
10.661 48 42282.145 0.686166 18 10 4
10.661 48 68886.531 0.684564 20 6 2
10.685 48 1498.741 0.677589 17 12 3
10.685 24 0.000 0.698189 19 9 0
10.685 24 0.000 0.703940 21 1 0
10.733 48 2054.250 0.511190 14 13 9
10.733 48 13512.246 0.608545 15 11 10
10.733 48 874.002 0.433063 15 14 5
10.733 48 8341.328 0.765061 17 11 6
10.733 48 3496.692 0.547453 18 11 1
10.733 48 5957.996 0.576137 19 7 6
10.733 48 1422.699 0.480880 19 9 2
10.733 48 29.973 2.908144 21 2 1
10.781 48 1753.692 0.672832 15 12 9
10.781 12 0.000 0.664564 15 15 0
10.781 48 1753.916 0.672832 16 13 5
10.781 48 40.492 0.672836 19 8 5
10.781 24 0.000 0.690230 20 5 5
10.781 48 3651.796 0.672833 20 7 1
10.781 24 0.000 0.707799 21 3 0
10.805 24 53989.008 0.671649 16 14 0
10.805 24 0.000 0.643831 18 8 8
10.805 48 9669.600 0.671649 20 6 4
10.829 24 13108.416 0.602064 15 15 2
10.829 48 7.145 1.037184 18 9 7
10.829 48 5454.835 0.565999 18 11 3
10.829 48 11229.632 0.748505 21 3 2
10.853 24 15315.414 0.565674 14 14 8
10.853 24 22509.879 0.761307 16 10 10
10.853 48 3599.667 0.465110 16 14 2
10.877 48 18971.936 0.668109 15 13 8
10.877 48 6143.480 0.668110 16 11 9
10.877 48 6768.690 0.668109 17 12 5
10.877 24 0.000 0.688241 17 13 0
10.877 48 21970.449 0.668109 19 9 4
10.877 48 6369.703 0.668110 20 7 3
10.877 48 1052.865 0.668110 21 4 1
10.900 48 3017.939 0.666934 18 10 6
10.924 48 19252.697 0.725662 17 13 2
10.924 48 25930.057 0.616213 19 10 1
10.948 48 84530.609 0.670057 16 12 8
10.948 24 97406.680 0.669681 20 8 0
10.971 24 0.000 0.685600 15 15 4
10.971 48 10599.709 0.663419 21 4 3
10.971 24 0.000 0.643299 21 5 0
10.995 48 100173.562 0.711064 16 14 4
10.995 24 166767.672 0.699932 18 12 0
10.995 48 159852.547 0.700750 20 8 2
11.018 48 1658.682 0.878372 15 14 7
11.018 48 17370.697 0.600177 17 10 9
11.018 48 597.529 0.367525 18 11 5
11.018 48 997.503 0.426396 19 10 3
11.018 48 3084.364 0.521335 21 5 2
11.042 48 106769.031 0.664861 18 12 2
11.042 24 137048.594 0.664281 20 6 6
11.065 48 2199.089 0.658762 16 13 7
11.065 48 22348.607 0.658762 17 11 8
11.065 48 5517.560 0.658762 17 13 4
11.065 48 2181.360 0.658762 19 8 7
11.065 48 792.531 0.658762 20 7 5
11.112 48 94.408 1.790153 19 9 6
11.112 48 8604.981 0.749077 21 6 1
11.135 48 12669.628 0.783483 20 8 4
11.158 24 0.000 0.646298 13 13 12
11.158 48 9421.717 0.654138 16 15 1
11.158 48 210.036 0.654139 17 12 7
11.158 24 0.000 0.710653 19 11 0
11.158 48 92.088 0.654137 20 9 1
11.158 48 2810.931 0.654138 21 5 4
11.181 24 0.000 0.727377 14 12 12
11.181 48 5305.505 0.652987 18 12 4
11.181 6 0.000 0.654133 22 0 0
11.204 48 15691.037 0.720583 14 13 11
11.204 24 33526.344 0.606793 15 15 6
11.204 48 3536.515 0.800657 17 14 1
11.204 24 287.097 0.250511 18 9 9
11.204 48 7806.454 0.750322 19 10 5
11.204 48 812.864 0.390627 19 11 2
11.204 48 8561.631 0.745729 21 6 3
11.204 24 20528.598 0.594559 22 1 1
11.227 48 14.069 2.231752 16 14 6
11.227 48 3712.752 0.899219 18 10 8
11.227 24 759.892 1.253699 22 2 0
11.250 48 487.995 0.649546 15 12 11
11.250 48 2527.021 0.649546 16 15 3
11.250 48 4456.323 0.649545 20 9 3
11.250 24 0.000 0.652303 21 7 0
11.273 24 0.000 0.789600 14 14 10
11.273 24 0.000 0.653409 22 2 2
11.296 48 706.681 0.367562 15 13 10
11.296 48 297.195 0.248811 17 13 6
11.296 48 7711.250 0.554293 17 14 3
11.296 48 2570.974 0.490818 18 11 7
11.296 48 8067.584 0.745120 18 13 1
11.296 48 22961.691 0.704436 21 7 2
11.296 48 87.333 0.057870 22 3 1
11.342 24 0.000 0.685507 16 11 11
11.342 48 332.649 0.644986 19 11 4
11.342 24 0.000 0.746241 20 7 7
11.364 48 159977.438 0.683996 16 12 10
11.364 48 130867.500 0.688308 20 8 6
11.364 24 72803.609 0.703794 20 10 0
11.364 24 23748.637 0.750559 22 4 0
11.387 48 1212.476 0.912182 15 14 9
11.387 48 14.392 4.882132 18 13 3
11.387 48 1769.587 0.455264 21 6 5
11.387 24 51470.301 0.681020 22 3 3
11.410 48 99131.141 0.647012 18 12 6
11.410 48 29163.473 0.651607 20 10 2
11.410 48 92906.383 0.647192 22 4 2
11.432 48 10999.578 0.640458 16 13 9
11.432 48 4368.299 0.640459 16 15 5
11.432 48 31403.549 0.640459 19 9 8
11.432 48 119.824 0.640458 19 12 1
11.432 48 2748.602 0.640458 20 9 5
11.432 48 346.130 0.640457 21 7 4
11.432 48 8339.304 0.640458 21 8 1
11.478 48 5003.653 0.521834 17 11 10
11.478 48 3454.537 0.793780 17 14 5
11.478 48 4331.639 0.776329 19 10 7
11.478 48 503.373 0.311025 22 5 1
11.500 12 218894.016 0.675633 16 16 0
11.522 24 0.000 0.639138 15 15 8
11.522 48 13001.710 0.635962 17 12 9
11.522 24 0.000 0.631716 17 15 0
11.522 48 1780.457 0.635963 19 12 3
11.522 48 6706.062 0.635963 21 8 3
11.545 48 8059.045 0.634844 16 14 8
11.545 24 0.000 0.604006 16 16 2
11.545 48 8097.510 0.634843 20 10 4
11.545 24 0.000 0.657954 22 4 4
11.567 48 3729.106 0.498882 17 15 2
11.567 48 8971.104 0.729207 18 13 5
11.567 48 7212.025 0.535136 19 11 6
11.567 48 5994.438 0.526029 22 5 3
11.590 24 1666.043 1.035919 18 14 0
11.590 24 518.950 0.157468 22 6 0
11.612 48 17038.049 0.631498 17 13 8
11.612 48 2803.506 0.631498 20 11 1
11.612 24 0.000 0.667525 21 9 0
11.634 24 0.000 0.638608 18 10 10
11.634 48 22130.412 0.630387 18 14 2
11.634 48 1981.995 0.630388 22 6 2
11.656 48 116.182 0.077477 18 11 9
11.656 48 11704.309 0.550516 21 7 6
11.656 48 4335.886 0.502654 21 9 2
11.678 24 116167.578 0.633385 16 16 4
11.678 24 107535.102 0.633591 20 8 8
11.700 48 4930.005 0.627065 16 15 7
11.700 48 4708.191 0.627065 17 15 4
11.700 48 28224.625 0.627065 19 12 5
11.700 24 0.000 0.596998 19 13 0
11.700 48 2646.947 0.627066 20 9 7
11.700 48 918.285 0.627066 20 11 3
11.700 48 9310.486 0.627065 21 8 5
11.700 24 0.000 0.611802 23 1 0
11.723 48 74428.492 0.687569 18 12 8
11.745 24 36005.445 0.578920 14 13 13
11.745 48 21234.656 0.565384 17 14 7
11.745 48 5825.659 0.746672 19 13 2
11.745 24 23817.857 0.568623 22 5 5
11.745 48 2404.847 0.819144 22 7 1
11.745 48 14864.012 0.554119 23 2 1
11.766 24 117187.000 0.629019 14 14 12
11.766 48 163276.672 0.628214 18 14 4
11.766 48 51583.629 0.631695 20 10 6
11.766 48 66068.156 0.630769 22 6 4
11.788 48 2424.192 0.622663 15 13 12
11.788 48 12175.204 0.622663 21 9 4
11.788 24 0.000 0.684423 23 3 0
11.832 48 5088.439 0.752444 15 14 11
11.832 48 461.719 0.273537 18 13 7
11.832 48 9.878 6.746250 19 10 9
11.832 48 13722.126 0.546340 21 10 1
11.832 48 9360.422 0.716544 22 7 3
11.832 48 285.144 0.203293 23 3 2
11.854 24 245.206 1.926058 16 12 12
11.854 24 5207.864 0.841583 20 12 0
11.876 48 2550.054 0.618292 16 13 11
11.876 48 1977.929 0.618292 17 16 1
11.876 48 31793.830 0.618292 19 11 8
11.876 48 2590.990 0.618292 19 13 4
11.876 48 4312.363 0.618292 20 11 5
11.876 48 68.349 0.618293 23 4 1
11.897 24 0.000 0.601473 16 16 6
11.897 48 325.180 0.617205 20 12 2
11.897 24 7184.949 0.617204 22 8 0
11.919 24 2757.569 0.800199 15 15 10
11.919 48 6327.436 0.507762 17 15 6
11.919 48 4055.922 0.482407 18 15 1
11.919 48 242.952 0.172778 21 10 3
11.941 48 7888.565 0.458459 16 14 10
11.941 48 109.173 2.873094 22 8 2
11.962 48 10852.533 0.613952 17 12 11
11.962 48 16083.902 0.613952 17 16 3
11.962 48 21445.627 0.613952 19 12 7
11.962 48 1.472 0.613933 21 8 7
11.962 48 37939.918 0.613952 23 4 3
11.962 24 0.000 0.616921 23 5 0
11.984 48 1296.799 0.612872 18 14 6
11.984 24 0.000 0.712688 22 6 6
12.006 48 82.405 2.026851 17 13 10
12.006 48 18215.365 0.546157 18 15 3
12.006 48 9411.382 0.709300 19 14 1
12.006 48 2181.693 0.822372 21 9 6
12.006 48 335.092 1.210647 22 7 5
12.006 48 8659.716 0.517815 23 5 2
12.027 48 95207.109 0.718642 20 12 4
12.048 48 62.905 0.609640 16 15 9
12.048 24 0.000 0.614455 20 9 9
12.048 24 0.000 0.621012 21 11 0
12.070 48 20326.646 0.733442 20 10 8
12.070 48 105682.531 0.661900 22 8 4
12.091 48 21634.326 0.671532 17 14 9
12.091 24 972.384 0.937657 18 11 11
12.091 48 1994.939 0.830481 19 13 6
12.091 48 0.028 1184.492188 19 14 3
12.091 48 3860.352 0.764279 21 10 5
12.091 48 3519.251 0.462570 21 11 2
12.091 48 775.554 0.981970 22 9 1
12.091 48 10175.889 0.519880 23 6 1
12.113 48 76630.094 0.613261 18 12 10
12.134 48 3683.249 0.605362 17 16 5
12.134 48 3521.664 0.605362 20 11 7
12.134 48 124.058 0.605363 20 13 1
12.134 48 3462.126 0.605362 23 5 4
12.176 48 55.890 0.000408 18 13 9
12.176 48 97.423 1.905106 18 15 5
12.176 48 14498.582 0.528827 22 9 3
12.176 48 4836.832 0.477510 23 6 3
12.198 24 146829.391 0.652999 16 16 8
12.198 6 48972.145 0.691485 24 0 0
12.219 48 2249.633 0.601113 17 15 8
12.219 12 0.000 0.648104 17 17 0
12.219 48 731.749 0.601113 20 13 3
12.219 48 2019.349 0.601113 21 11 4
12.219 24 0.000 0.620825 23 7 0
12.219 24 0.000 0.585538 24 1 1
12.240 24 2309.716 0.600055 18 16 0
12.240 48 12503.156 0.600055 20 12 6
12.240 24 392.155 0.600055 24 2 0
12.261 24 8259.682 0.705990 17 17 2
12.261 48 1170.735 0.357480 19 11 10
12.261 48 3571.711 0.453139 19 14 5
12.261 24 23870.430 0.540168 22 7 7
12.261 48 11697.327 0.515888 23 7 2
12.282 48 1961.352 0.301177 18 14 8
12.282 48 2399.438 0.325265 18 16 2
12.282 48 8974.643 0.446666 22 8 6
12.282 24 20183.016 0.494621 22 10 0
12.282 24 4232.669 0.857612 24 2 2
12.303 48 8762.402 0.596893 19 12 9
12.303 24 0.000 0.652194 19 15 0
12.303 48 16475.064 0.596893 21 9 8
12.303 48 21749.035 0.596893 21 12 1
12.303 48 0.566 0.596886 24 3 1
12.324 8 0.000 0.643371 14 14 14
12.324 48 12293.235 0.595843 22 10 2
12.345 48 147.154 1.621171 15 14 13
12.345 48 1708.446 0.389242 19 15 2
12.345 48 516.791 1.079764 21 10 7
12.345 48 11336.400 0.509842 22 9 5
12.345 48 18111.604 0.527044 23 6 5
12.366 24 58389.941 0.601851 24 4 0
12.387 24 0.000 0.671713 15 15 12
12.387 24 0.000 0.735260 16 13 13
12.387 48 12095.004 0.592703 17 16 7
12.387 24 0.000 0.624602 17 17 4
12.387 48 6109.655 0.592703 19 13 8
12.387 48 1167.612 0.592703 20 13 5
12.387 48 486.571 0.592703 21 12 3
12.387 48 968.762 0.592702 23 7 4
12.387 48 2560.149 0.592703 23 8 1
12.387 24 0.000 0.590848 24 3 3
12.408 48 33698.273 0.690808 16 14 12
12.408 48 106267.023 0.646548 18 16 4
12.408 24 27260.693 0.702372 20 14 0
12.408 48 100993.477 0.647996 24 4 2
12.428 48 8679.107 0.493428 18 15 7
12.428 48 9943.497 0.499549 21 11 6
12.449 24 29955.414 0.601035 20 10 10
12.449 48 71944.438 0.596959 20 14 2
12.449 48 57275.621 0.597853 22 10 4
12.470 48 789.999 0.588542 16 15 11
12.470 48 10761.667 0.588542 17 13 12
12.470 48 26082.959 0.588543 19 15 4
12.470 48 4142.475 0.588543 20 11 9
12.470 48 3639.900 0.588543 23 8 3
12.470 48 17676.646 0.588542 24 5 1
12.511 48 236.687 0.130661 17 14 11
12.511 48 4685.903 0.455545 19 14 7
12.511 48 15084.020 0.511461 22 11 1
12.532 48 4970.170 0.384055 20 12 8
12.532 24 2797.412 0.326392 24 4 4
12.552 48 42385.359 0.584411 21 12 5
12.552 24 0.000 0.603295 21 13 0
12.552 24 0.000 0.614612 23 9 0
12.552 48 1401.835 0.584411 24 5 3
12.573 24 0.000 0.628182 16 16 10
12.573 24 0.000 0.574667 18 12 12
12.573 48 27.689 0.583384 20 14 4
12.573 24 0.000 0.593523 22 8 8
12.573 24 974.674 0.583382 24 6 0
12.594 48 7713.880 0.696704 17 15 10
12.594 24 625.189 0.261879 17 17 6
12.594 48 342.001 0.179135 18 13 11
12.594 48 2097.192 0.810675 18 17 1
12.594 48 486.338 1.098598 21 13 2
12.594 48 19185.775 0.653674 22 9 7
12.594 48 257.410 0.137906 22 11 3
12.594 48 21796.854 0.649147 23 7 6
12.594 48 2041.253 0.814043 23 9 2
12.614 48 42259.598 0.660673 18 16 6
12.614 48 1627.606 1.038560 24 6 2
12.635 48 9.570 0.580310 19 16 1
12.635 48 4582.844 0.580308 20 13 7
12.635 48 424.410 0.580309 23 8 5
12.655 48 31903.266 0.579287 18 14 10
12.655 48 5253.120 0.579288 22 10 6
12.675 48 2611.902 0.404740 18 17 3
12.675 48 13236.974 0.497375 19 15 6
12.675 48 1066.550 0.911786 21 10 9
12.716 48 6204.871 0.576235 17 16 9
12.716 48 6512.263 0.576235 19 12 11
12.716 48 7935.641 0.576235 19 16 3
12.716 48 1384.119 0.576235 20 15 1
12.716 48 6295.197 0.576235 21 11 8
12.716 48 24718.215 0.576235 21 13 4
12.716 48 5930.177 0.576235 23 9 4
12.716 24 0.000 0.586224 24 5 5
12.716 48 14408.281 0.576235 24 7 1
12.716 24 0.000 0.599632 25 1 0
12.736 24 120591.719 0.628172 22 12 0
12.736 48 34299.668 0.676419 24 6 4
12.757 48 156.663 0.060188 18 15 9
12.757 48 3943.698 0.739427 19 13 10
12.757 48 10074.289 0.481447 22 11 5
12.757 48 1050.248 0.913322 23 10 1
12.757 48 14211.491 0.658643 25 2 1
12.777 48 118851.406 0.579183 20 14 6
12.777 48 81153.352 0.580444 22 12 2
12.797 48 1636.235 0.572190 20 15 3
12.797 48 5481.227 0.572190 21 12 7
12.797 48 18053.295 0.572190 24 7 3
12.797 24 0.000 0.589848 25 3 0
12.837 48 7405.788 0.689516 18 17 5
12.837 48 2624.350 0.395211 19 14 9
12.837 48 2138.960 0.801308 21 14 1
12.837 48 10.241 1.610380 23 10 3
12.837 48 475.319 1.107143 25 3 2
12.857 24 400540.719 0.601597 24 8 0
12.878 24 0.000 0.594411 17 17 8
12.878 48 0.125 0.568303 19 16 5
12.878 24 0.000 0.668721 20 11 11
12.878 48 9156.891 0.568173 23 8 7
12.878 48 626.749 0.568173 25 4 1
12.898 48 752.168 0.567174 18 16 8
12.898 48 5026.329 0.567173 20 12 10
12.898 48 1582.216 0.567173 22 12 4
12.898 48 10010.276 0.567173 24 8 2
12.918 24 201.271 1.472532 15 15 14
12.918 48 222.695 0.098690 21 13 6
12.918 48 6281.894 0.448796 21 14 3
12.918 24 124.054 0.027292 22 9 9
12.918 48 7534.736 0.458458 23 9 6
12.938 24 6316.977 0.786949 16 14 14
12.938 12 29334.926 0.663537 18 18 0
12.938 48 1147.344 1.143397 22 10 8
12.938 24 3723.319 0.331233 24 6 6
12.958 48 19312.793 0.564184 16 15 13
12.958 48 1304.676 0.564185 19 15 8
12.958 24 0.000 0.629484 19 17 0
12.958 48 7762.704 0.564184 20 13 9
12.958 48 0.036 0.564203 20 15 5
12.958 24 0.000 0.552595 23 11 0
12.958 48 1174.496 0.564185 24 7 5
12.958 48 1864.236 0.564185 25 4 3
12.958 24 0.000 0.610466 25 5 0
12.977 24 0.000 1.007893 18 18 2
12.997 48 1257.012 0.317586 17 14 13
12.997 48 803.021 0.965953 19 17 2
12.997 48 3526.509 0.407812 22 11 7
12.997 48 33105.316 0.617943 22 13 1
12.997 48 3190.903 0.400563 23 10 5
12.997 48 2609.250 0.385052 23 11 2
12.997 48 628.391 0.236561 25 5 2
13.017 24 57828.707 0.570060 16 16 12
13.017 24 20781.746 0.576011 20 16 0
13.017 48 89788.805 0.568307 24 8 4
13.037 48 46834.953 0.560224 17 15 12
13.037 48 4458.039 0.560224 24 9 1
13.057 48 64854.797 0.633972 20 14 8
13.057 48 1893.262 1.063184 20 16 2
13.077 24 11205.765 0.468121 18 13 13
13.077 48 2285.663 0.369342 18 17 7
13.077 48 339.441 0.148501 21 11 10
13.077 48 8367.582 0.454671 21 14 5
13.077 48 7108.424 0.446360 22 13 3
13.077 48 21281.016 0.492060 25 6 1
13.096 48 68508.367 0.565478 18 14 12
13.096 24 155737.828 0.562708 18 18 4
13.096 48 103275.125 0.563951 22 12 6
13.116 48 11.861 0.556303 17 16 11
13.116 48 9.261 0.556294 19 16 7
13.116 48 858.044 0.556291 19 17 4
13.116 48 3173.471 0.556291 21 12 9
13.116 24 0.000 0.612472 21 15 0
13.116 48 3159.007 0.556291 23 11 4
13.116 48 11758.186 0.556291 24 9 3
13.116 48 817.126 0.556292 25 5 4
13.155 48 2785.415 0.760410 18 15 11
13.155 48 12620.544 0.468744 21 15 2
13.155 48 9819.387 0.660046 25 6 3
13.175 48 16619.480 0.433205 20 16 4
13.195 48 1551.486 0.552387 19 13 12
13.195 48 2451.736 0.552386 20 15 7
13.195 48 8249.820 0.552386 21 13 8
13.195 48 4248.809 0.552386 23 9 8
13.195 48 21355.518 0.552386 23 12 1
13.195 24 0.000 0.511862 24 7 7
13.195 24 0.000 0.532794 25 7 0
13.214 48 45.962 0.551413 24 8 6
13.214 24 1375.366 0.551414 24 10 0
13.214 6 0.000 0.581809 26 0 0
13.234 24 125.795 1.815142 17 17 10
13.234 48 3624.754 0.730497 19 14 11
13.234 48 20871.211 0.622388 22 13 5
13.234 48 16.691 0.846322 23 10 7
13.234 48 1982.095 0.799991 25 7 2
13.234 24 4838.143 0.704812 26 1 1
13.253 48 6534.444 0.364097 18 16 10
13.253 24 2222.547 0.955805 22 14 0
13.253 48 24731.307 0.449439 24 10 2
13.253 24 3270.319 0.298545 26 2 0
13.273 48 7338.060 0.548509 21 15 4
13.273 48 4058.442 0.548509 23 12 3
13.273 48 14392.250 0.548509 24 9 5
13.292 24 0.000 0.517475 18 18 6
13.292 24 0.000 0.554471 22 10 10
13.292 48 5628.166 0.547544 22 14 2
13.292 24 0.000 0.523200 26 2 2
13.311 48 727.746 0.986569 19 15 10
13.311 48 5154.896 0.696772 19 17 6
13.311 48 1101.667 0.894374 19 18 1
13.311 48 81.138 2.249579 21 14 7
13.311 48 2995.469 0.747331 22 11 9
13.311 48 1254.767 0.869951 23 11 6
13.311 48 2509.770 0.767480 25 6 5
13.311 48 2730.918 0.757605 26 3 1
13.331 24 43466.543 0.728675 20 12 12
13.350 48 216.962 0.544659 20 13 11
13.350 48 1647.968 0.544658 20 17 1
13.350 48 4276.428 0.544658 25 7 4
13.350 48 15238.885 0.544658 25 8 1
13.370 48 84474.062 0.610509 20 16 6
13.370 48 74722.828 0.614865 22 12 8
13.370 48 40800.441 0.641053 24 10 4
13.370 24 179990.234 0.589052 26 4 0
13.389 48 404.751 0.154578 18 17 9
13.389 48 4652.005 0.403703 19 18 3
13.389 24 391.049 1.179888 26 3 3
13.408 48 232808.609 0.546398 20 14 10
13.408 48 99927.586 0.548833 22 14 4
13.408 48 188579.734 0.546911 26 4 2
13.427 48 641.491 0.540835 19 16 9
13.427 48 3573.906 0.540835 20 17 3
13.427 48 1412.832 0.540835 21 16 1
13.427 48 2392.580 0.540836 23 12 5
13.427 24 0.000 0.583230 23 13 0
13.427 48 13539.380 0.540835 24 11 1
13.427 48 4325.955 0.540835 25 8 3
13.466 48 24231.410 0.475178 21 15 6
13.466 48 2542.877 0.355048 22 13 7
13.466 48 9143.574 0.651304 23 13 2
13.466 48 44289.496 0.491389 26 5 1
13.485 24 386159.906 0.572714 24 8 8
13.504 24 0.000 0.537275 16 15 15
13.504 48 79.871 0.537040 20 15 9
13.504 48 17695.494 0.537038 21 12 11
13.504 48 895.248 0.537037 21 16 3
13.504 48 12284.010 0.537038 24 9 7
13.504 48 2309.103 0.537038 24 11 3
13.504 24 0.000 0.591044 25 9 0
13.523 24 0.000 0.402816 16 16 14
13.523 24 0.000 0.504490 26 4 4
13.542 48 9775.623 0.436148 17 15 14
13.542 48 7637.786 0.659352 19 18 5
13.542 48 804.636 0.959318 21 13 10
13.542 48 3708.660 0.379537 22 15 1
13.542 48 7600.334 0.423640 23 10 9
13.542 48 924.706 0.926972 25 7 6
13.542 48 1905.093 0.325385 25 9 2
13.542 48 2232.918 0.339564 26 5 3
13.561 24 21222.117 0.656080 18 18 8
13.561 48 7232.998 0.354522 24 10 6
13.561 24 9536.754 0.720598 26 6 0
13.580 48 7.376 0.533264 17 16 13
13.580 48 8833.213 0.533268 19 17 8
13.580 48 507.977 0.533268 20 17 5
13.580 48 4764.306 0.533268 23 11 8
13.580 48 16286.252 0.533268 23 13 4
13.580 48 13357.467 0.533268 25 8 5
13.599 24 0.000 0.552398 18 14 14
13.599 48 240.404 0.532327 22 14 6
13.599 48 142.940 0.532329 26 6 2
13.618 48 1512.714 0.298800 18 15 13
13.618 48 13081.259 0.444836 21 14 9
13.618 48 5860.847 0.404917 22 15 3
13.637 48 68102.938 0.539222 20 16 8
13.637 24 161144.547 0.536148 24 12 0
13.656 24 0.000 0.567319 17 17 12
13.656 12 0.000 0.641853 19 19 0
13.656 48 1567.160 0.529524 21 16 5
13.656 48 1346.462 0.529525 23 12 7
13.656 48 21105.100 0.529525 24 11 5
13.656 48 24573.625 0.529525 25 9 4
13.675 48 55987.738 0.613047 18 16 12
13.675 24 294443.688 0.564651 20 18 0
13.675 48 83364.617 0.597342 24 12 2
13.694 48 21622.609 0.459374 19 14 13
13.694 24 833.602 0.229171 19 19 2
13.694 24 41057.730 0.580182 22 11 11
13.694 48 258.246 0.079186 23 14 1
13.694 48 12107.555 0.626306 25 10 1
13.694 24 5123.758 0.392501 26 5 5
13.694 48 1369.571 0.284268 26 7 1
13.713 48 123013.125 0.533304 20 18 2
13.713 48 71978.680 0.535331 22 12 10
13.713 48 60597.105 0.536107 26 6 4
13.732 48 36370.348 0.525808 19 15 12
13.732 48 6984.553 0.525808 21 15 8
13.732 24 0.000 0.556326 21 17 0
13.732 24 0.000 0.520802 27 1 0
13.769 48 805.905 0.956158 18 17 11
13.769 48 3654.744 0.365521 19 18 7
13.769 48 11134.515 0.429707 21 17 2
13.769 48 514.928 1.084913 22 13 9
13.769 48 1388.901 0.841474 22 15 5
13.769 48 5090.758 0.680468 23 13 6
13.769 48 118.760 1.919536 23 14 3
13.769 48 2065.000 0.778836 25 10 3
13.769 48 1673.389 0.810194 26 7 3
13.769 48 1752.430 0.304277 27 2 1
13.788 48 115.481 0.140374 24 12 4
13.807 48 13001.454 0.522117 19 16 11
13.807 24 0.000 0.550671 19 19 4
13.807 24 0.000 0.644168 20 13 13
13.807 48 3757.070 0.522116 20 17 7
13.807 24 0.000 0.619965 24 9 9
13.807 48 166.534 0.522118 25 8 7
13.807 24 0.000 0.642318 27 3 0
13.825 48 432.177 0.521196 20 14 12
13.825 48 12337.985 0.521198 20 18 4
13.825 24 6312.746 0.521198 22 16 0
13.825 48 1634.139 0.521198 24 10 8
13.825 24 6345.703 0.521198 26 8 0
13.844 48 1706.373 0.805092 25 9 6
13.844 48 3221.206 0.721440 27 3 2
13.863 48 3274.244 0.263069 22 14 8
13.863 48 3664.048 0.274854 22 16 2
13.863 48 4512.130 0.805788 26 8 2
13.881 48 202.873 0.518453 20 15 11
13.881 48 53.521 0.518447 21 16 7
13.881 48 164.154 0.518448 21 17 4
13.881 48 4790.123 0.518451 24 11 7
13.881 48 143.928 0.518449 24 13 1
13.881 24 0.000 0.551518 25 11 0
13.881 48 5639.001 0.518451 27 4 1
13.900 24 0.000 0.480472 18 18 10
13.900 24 0.000 0.536103 26 6 6
13.919 48 4715.988 0.374636 19 17 10
13.919 48 1152.720 0.252972 23 11 10
13.919 48 14412.105 0.608090 23 14 5
13.919 48 8711.256 0.409941 25 10 5
13.919 48 12031.894 0.425067 25 11 2
13.919 48 997.063 0.904037 26 7 5
13.956 48 1712.422 0.514812 21 13 12
13.956 48 51.106 0.514808 23 12 9
13.956 24 0.000 0.493669 23 15 0
13.956 48 3151.326 0.514811 24 13 3
13.956 48 4966.151 0.514812 27 4 3
13.956 24 0.000 0.544797 27 5 0
13.974 48 184851.469 0.560507 20 16 10
13.974 48 70229.789 0.590530 22 16 4
13.974 48 102809.844 0.576854 24 12 6
13.974 48 77585.695 0.586684 26 8 4
13.993 24 27149.818 0.579212 19 19 6
13.993 48 2431.204 0.749978 21 14 11
13.993 48 34488.859 0.457621 22 15 7
13.993 48 364.326 0.108468 23 15 2
13.993 48 1418.150 0.270548 26 9 1
13.993 48 9578.805 0.410609 27 5 2
14.011 48 46500.246 0.523157 20 18 6
14.030 48 0.084 0.511293 20 19 1
14.030 48 3199.971 0.511198 23 13 8
14.030 48 13297.949 0.511198 25 11 4
-------------- next part --------------
/*******************************************************************************
*
* McStas, version 1.2 released February 2000
* Maintained by Kristian Nielsen and Kim Lefmann,
* Risoe National Laboratory, Roskilde, Denmark
*
* %IDENTIFICATION
*
* Author: <a href="mailto:hansen at ill.fr">Thomas C Hansen</a>
* Date: 08 March 2000
* Version: $Revision: 1.1 $
* Origin: <a href="http://www.ill.fr">ILL</a> (Dif/<a href="http://www.ill.fr/YellowBook/D20">D20</a>)
*
* General powder sample in incoherent scattering cylindrical can
*
* %DESCRIPTION
*
* This is a general powder sample in incoherent scattering cylindrical vanadium can.
* It creates elastic coherent scattering from a table of structure factors and incoherent
* background scattering from the sample itself and its container. To be efficient, the
* component is focussing on a given target, normally corresponding to a PSD or a multi-
* detector bank. Absorption is considered, as well attenuation of the beam, so the
* probability of scattering is higher closer to the incoming beam. Multiple scattering
* (and so secondary extinction) is not considered yet, nor diffuse scattering, or elastic
* coherent scattering from a sample can, nor any inelastic scattering or elastic scattering
* from amorphous materials. It is planned to implement multiple scattering and elastic
* coherent scattering from the sample can in the near future. Transmitted neutrons are
* normally not created, but only scattered neutrons leaving towards the detector target.
*
* %PARAMETERS
*
* INPUT PARAMETERS
*
* radius: (m) Radius of sample in (x,z) plane (0.005)
* h: (m) Height of sample y direction (0.05)
* pack: (1) Packing factor (1)
* Vc: (AA**3) Volume of unit cell
* sigma_a: (fm**2) Absorption cross section per unit cell at 2200 m/s
* j: (*1) Array of multiplicities
* q: (*1/AA) Array of wavevectors
* F2: (*fm**2) Array of structure factors
* DW: (*1/AA**2) Array of Debye-Waller factors
* write: (1) Output flag for writing in file sample.dat (0)
* transmission: (1) Ratio of transmitted neutron histories (0)
* nbInt: (1) Number of Bragg relections (size of the arrays)
* sigma_i: (fm**2) Incoherent scattering cross section of the sample per unit cell
* d_V: (m) Thickness of the Vanadium sample can (0.0001)
* broadening: (1) FWHM of the lattice constant variations for strain broadening (0)
* ttmin: (deg) Minimum 2theta angle of target <a href="../monitors/PSD_curved.html">PSD_curved.comp</a> (0)
* ttmax: (deg) Maximum 2theta angle of target <a href="../monitors/PSD_curved.html">PSD_curved.comp</a> (156.36)
* PSD_r: (m) Radius of curved linar PSD target <a href="../monitors/PSD_curved.html">PSD_curved.comp</a> (1.471)
* PSD_h: (m) Height of curved linear PSD target <a href="../monitors/PSD_curved.html">PSD_curved.comp</a> (0.15)
* sign: (1) Chirality of 1st diffractometer axis=sign[takeoff of <a href="../d20adapt.instr">d20adapt.instr]</a> (-1)
*
* OUTPUT PARAMETERS
*
* my_s_v2: (m/s**2) Attenuation factor due to elastic coherent scattering, multiplied by neutron velocity**2
* my_a_v: (1/s) Attenuation factor due to absorption, multiplied by neutron velocity
* q_v: (m/s) Corresponding velocity of wavevector Q
* my_i: (1/m) Attenuation factor due to elastic incoherent scattering
*
* %LINKS
* <a href="../d20adapt.instr">Source code of d20adapt.instr</a>, where this component is used
* %LINKS
* A possible target detector: <a href="../monitors/PSD_curved.html">PSD_curved.comp</a>
* %LINKS
* <A HREF="http://neutron.risoe.dk/mcstas/components/Powder1/">The original component Powder1</A>, with explanation of focusing.
* %LINKS
* <A HREF="http://neutron.risoe.dk/mcstas/components/tests/powder/">Test results from Powder1</A> (not up-to-date).
*
* %END
*
*******************************************************************************/
DEFINE COMPONENT Powder0
DEFINITION PARAMETERS (radius, h, pack, Vc, sigma_a, j, q, F2, DW,write,transmission,nbInt,sigma_i,d_V,broadening,ttmin,ttmax)
SETTING PARAMETERS (PSD_r, PSD_h, sign)
OUTPUT PARAMETERS (my_s_v2, my_a_v, q_v,my_i)
STATE PARAMETERS (x,y,z,vx,vy,vz,t,s1,s2,p)
DECLARE
%{
#if __dest_os == __mac_os
/* void Event_loop(); */
#endif
double my_i;
unsigned long end_counter=0;
double my_s_v2,my_s_v2_total, my_a_v, q_v;
double d_phi0,d_Scherrer,total_F2;
int total,sign2,i;
FILE *outfile;
double V_pack=1.0; /* Vanadium sample can */
double V_sigma_a = 5.08; /* Absorption cross section per atom (barns) */
double V_sigma_i = 4.935 ; /* Incoherent scattering cross section per atom (barns) */
double V_rho ; /* Density of atoms (AA-3) */
double V_my_s ;
double V_my_a_v,s ;
%}
INITIALIZE
%{
V_rho =(2.0 * V_pack/(3.024*3.024*3.024)); /* Density of atoms (AA-3) */
V_my_s =(V_rho * 100.0 * V_sigma_i);
V_my_a_v =(V_rho * 100.0 * V_sigma_a * 2200.0);
my_i= sigma_i * pack/Vc * 100.0 ;
outfile=fopen("sample.dat","w");
my_a_v = sigma_a * pack/Vc * 2200.0; /* Is not yet divided by v (which changes) */
my_s_v2_total=0.0;
for (i=1,my_s_v2=0.0, total_F2 = 0.0;(i<=nbInt);i++)
{
if ((i<2)||(i>(nbInt-1))) printf("i=%4d Q=%9.2f j=%4d F2=%14.2f DW=%9.2f %9.5f\n",i,q[i],j[i],F2[i],DW[i],my_a_v);
}
d_phi0 = RAD2DEG*atan(0.5*(h+ PSD_h)/PSD_r);
total = 0;
printf("Vertical divergence (sample-PSD): %lfdeg, ",d_phi0);
printf("my(incoh.=)%lg/m, sigma(incoh.)=%lgbarn\n",my_i,sigma_i);
%}
TRACE
%{
double t20, t21,t2i0,t2i1,t0, t1,ti0,ti1, v, l_full,l_powder,l_can, l, l_1,l_2,l_p,l_c, dt, d_phi, theta, my_s, my,muR,muV;
double aim_x, aim_y, aim_z, axis_x, axis_y, axis_z,my_total_powder,my_total_can;
double tmp_vx, tmp_vy, tmp_vz, vout_x, vout_y, vout_z, Choice,total_F2_1,X;
int i,coherent,powder;
if (vz<0) ABSORB;
my_s_v2_total=0.0;
v = sqrt(vx*vx + vy*vy + vz*vz); /* Very first, we should provide the neutron velocity ... */
for (i=1,my_s_v2=0.0, total_F2 = 0.0;i<=nbInt;i++) if (q[i]*K2V<(2.0*v))
{
my_s_v2_total += PI*PI*PI*pack*j[i]*F2[i]*DW[i]/(Vc*Vc*V2K*V2K*q[i]); /* all 0 < 2theta < 180 deg */
if ((fabs(q[i]*K2V) < (sin(fabs(ttmax)/2*PI/180)*2.0*v)) && (fabs(q[i]*K2V) > (sin(fabs(ttmin)/2*PI/180)*2.0*v)))
total_F2 += j[i]*F2[i]*DW[i]/q[i]; /* only ttmin < 2theta < ttmax */
}
total_F2_1=total_F2;
total++;
if (cylinder_intersect(&t0, &t1, x, y, z, vx, vy, vz, radius+d_V, h))
{
if (t0 < 0) ABSORB; /* Neutron enters at t=t0. */
l_full = v * (t1 - t0); /* Length of full path through sample AND can - if there won't be any scattering */
if (rand01() < (double)transmission) /* Transmission (IF any ...) */
{
my=my_i+my_a_v/v+my_s_v2_total/v/v; /* not exactly correct yet: considering can as sample ... */
dt = l_full/v + t0;
PROP_DT(dt);
p*=exp(-my*l_full)/(double)transmission;
}
else
{
/* Diffraction (elastic coherent AND incoherent scattering) */
/* it is sufficient to look here for choosing kind of scattering ... */
/* First, we must choose between scattering from sample or from container */
/* Therefore, we need (mu_i+mu_c)*R and mu_i(V)*d(V) and a random choice as done for coherent/incoherent */
/* To do it properly, we need the path of an unscattered neutron in Vanadium and in sample */
/* To do it even better, we may even consider 'multiple scattering', but that's valid for the sample itself beforehand ... */
if (cylinder_intersect(&ti0, &ti1, x, y, z, vx, vy, vz, radius, h))
{
if (ti0 < 0) ABSORB;
if (ti0 < t0) printf("\n%lg / %lg - %lg / %lg\n",t0,ti0,ti1,t1);
l_powder = v * (ti1 - ti0); /* Length of full path through sample AND can */
}
else
{
ti0=(t0+t1)/2.0;
ti1=ti0;
l_powder=0.0;
}
if (ti0 < t0) printf("\n---> %12.6lg / %12.6lg - %12.6lg - %12.6lg / %12.6lg ---> %8.6lg %8.6lg %8.6lg - %8.6lg\n",t0,ti0,dt,ti1,t1,x,y,z,l_powder);
l_can=l_full-l_powder;
muR=(my_s_v2_total/v/v+my_i)*l_powder;
muV=(V_my_s)*l_can;
Choice = fabs((muR+muV)*randpm1());
if (muR>=Choice) /****************** Scattering in Powder Sample ******************/
{
powder=1;
p*= (muR+muV)/muR ;
/*** Now, we must choose between coherent or incoherent scattering ***/
Choice = fabs((my_s_v2_total/v/v+my_i)*randpm1());
if ((my_s_v2_total/(v*v))>=Choice) /****************** Coherent Scattering in Powder Sample ******************/
{ coherent=1;
if (ti0 < t0) printf("\n>>>> %lg / %lg - %lg / %lg\n",t0,ti0,ti1,t1);
/**** In case of coherent scattering, choose a Bragg reflection ****/
total_F2=total_F2_1;
Choice = fabs(total_F2 * randpm1());
for (i=1,total_F2 = 0.0; i<=nbInt; i++)
{
if ((fabs(q[i]*K2V) < (sin(fabs(ttmax)/2*PI/180)*2.0*v)) && (fabs(q[i]*K2V) > (sin(fabs(ttmin)/2*PI/180)*2.0*v)))
{
total_F2 += j[i]*F2[i]*DW[i]/q[i];
if (total_F2 >= Choice)
{
break;
}
}
}
/*printf("i=%3d, F2=%10.1lf, Q=%4.1lf, p=%10.2lg, ",i,F2[i],q[i],p);*/
if (i==0) i=1;
if (i>nbInt) i=nbInt;
p*=total_F2_1/(j[i]*F2[i]*DW[i]/q[i])*(my_s_v2_total/(v*v)+my_i)/(my_s_v2_total/(v*v));
/*printf("p=%10.2lg, ",p);*/
q_v = q[i]*K2V;
/* NEW: quick and dirty strain peak broadening (TH 22/10/99, mod 25/10/99) */
do
{
s = rand01();
s = 2*s - 1;
} while(s == 0);
X =sqrt(1.0/fabs(s)-1.0)/3.1415*s/fabs(s); /* Lorentzian distributuion */
q_v=1.0/(1.0/q_v*(1.0+broadening*X));
if ((2.0*v)<q_v)
{
ABSORB;
}
my_s_v2= PI*PI*PI*pack*j[i]*F2[i]*DW[i]/(Vc*Vc*V2K*V2K*q[i]);
my=my_s_v2/v/v+my_a_v/v; /* my_s_v2 only known for coherent case */
}
else /****************** Incoherent Scattering in Powder Sample ******************/
{ coherent=0;
my=my_i+my_a_v/v;
if (ti0 < t0) printf("\n===> %lg / %lg - %lg / %lg\n",t0,ti0,ti1,t1);
/* in case of incoherent scattering choose random a scattering angle */
/* if ever possible choose it from the target solid angle ... but this will be difficult */
theta = (ttmin/2*+rand01()*(ttmax-ttmin)/2)*DEG2RAD; /* Bragg scattering law */
q_v = sin(theta)*2.0*v; /* Bragg scattering law */
/*printf("Incoherent in powder p=%10.2lg, ",p);*/
p*=(my_s_v2_total/(v*v)+my_i)/(my_i)*sin(2.0*theta);
/*printf("p=%10.2lg, ",p);*/
p*=PI/2*(-cos(fabs(ttmax/2)/180*PI)+cos(fabs(ttmin/2)/180*PI));
if (q[0]<0)
{
theta = -theta;
q_v = -q_v;
}
}
dt = 1.0/my/v * -log(rand01()*(1-exp(-my*l_powder))+exp(-my*l_powder)) + ti0;
/* dt is NOT a difference, but the absolute time of scattering */
PROP_DT(dt); /* Point of scattering */
SCATTER;
l = v*(dt -t0 ); /* Penetration in sample AND can */
l_c = v*(ti0-t0 ); /* Penetration in can */
l_p = v*(dt -ti0); /* Penetration in sample */
}
else /****************** Incoherent Scattering in Vanadium Sample Can ******************/
{
powder=0;
p*= (muR+muV)/muV ;
coherent=0; /* assume no Bragg scattering by Vanadium can - may be changed later! */
my=my_i+my_a_v/v;
/* in case of incoherent scattering choose random a scattering angle */
theta = (ttmin/2+rand01()*(ttmax-ttmin)/2)*DEG2RAD; /* Bragg scattering law */
q_v = sin(theta)*2.0*v; /* Bragg scattering law */
/* printf("Incoherent in sample can: 2theta=%lf p=%10.2lg, ",2*theta/DEG2RAD,p); */
p*=sin(2.0*theta);
/* printf("p=%10.2lg, ",p); */
p*=PI/2*(-cos(fabs(ttmax/2)/180*PI)+cos(fabs(ttmin/2)/180*PI));
/* printf("p=%10.2lg, ",p); */
/* p is correct only if there is no elastic scattering */
if (q[0]<0)
{
theta = -theta;
q_v = -q_v;
}
/* simple choice, but modify p depending on scattering in entering or outgoing wall */
dt = rand01()*((t1-t0)-(ti1-ti0))+t0;
if (dt>ti0) dt+=(ti1-ti0);
PROP_DT(dt); /* Point of scattering */
SCATTER;
l = v*(dt -t0 ); /* Penetration in sample AND can */
if (dt>ti1)
{
l_c = v*(dt - t0 - (ti1-ti0) ); /* Penetration in can */
l_p = v*(ti1-ti0);
}
else
{
l_c = v*(ti0-t0 ); /* Penetration in can */
l_p = v*(dt -ti0); /* Penetration in sample */
}
}
}
if (write>=1) fprintf(outfile,"%lf %lf %lf\n ", x,y,z);
theta = asin(q_v/(2.0*v)); /* Bragg scattering law */
d_Scherrer=2.0*PSD_r*sin(2.0*theta);
if (d_Scherrer<PSD_h) d_phi0 = 180.0; /* in this case take full cone! */
else d_phi0 = 2.0*RAD2DEG*asin(PSD_h/d_Scherrer); /* otherwise take only less than half of the cone */
d_phi = d_phi0*DEG2RAD/2.0*randpm1(); /* for very small scattering angles d_phi0 is variable ! */
p *= d_phi0/360.0;
if (coherent == 0) i=0; /* not important, but must be an existing index ... */
if (q[i]<0) sign2=sign; /* we only consider an one-sided detector (in the higher resolution part) */
else sign2=-sign;
aim_x = (double)sign2*PSD_r*sin(2*theta)-x; /* Vector pointing at target (anal./det.) */
aim_y = -y ;
/* aim_z = fabs(PSD_r*cos(2*theta)-z); <- may be obsolete now */
aim_z = PSD_r*cos(2*theta)-z; /* the new way ... */
vec_prod(axis_x, axis_y, axis_z, vx, vy, vz, aim_x, aim_y, aim_z);
rotate(tmp_vx, tmp_vy, tmp_vz, vx, vy, vz, 2*theta, axis_x, axis_y, axis_z);
rotate(vout_x, vout_y, vout_z, tmp_vx, tmp_vy, tmp_vz, d_phi, vx, vy, vz);
vx = vout_x;
vy = vout_y;
vz = vout_z;
if(!cylinder_intersect(&t2i0, &t2i1, x, y, z,vout_x, vout_y, vout_z, radius, h)) /* a modifier ... */
{
/* Do not exit here, the neutron may have hit only the wall of the sample can (TH 09/10/99) */
l_1=0.0; /* Only the sample can is hit by the neutron, there is not trace through the powder sample AFTER scattering */
l_2=0.0;
t2i0=0.0;
t2i1=0.0;
}
else
{
l_1 = v*t2i1; /* l_1 is the length of trace in the powder sample itself AFTER scattering */
if (t2i0>0)
{
l_2 = v*t2i0; /* l_3 is the length or trace in the entering wall of the sample can, after scattering */
l_1-= v*t2i0; /* l_1 should not contain the part of the trace in the sample can if scattering happens in the entering wall */
}
else l_2=0.0; /* l_3 should be zero if the powder itself scatters, or the exiting wall of the sample can */
}
if(!cylinder_intersect(&t20, &t21, x, y, z,vout_x, vout_y, vout_z, radius+d_V, h)) /* a modifier ... */
{
if (ti0 < t0) printf("\n---> %12.6lg / %12.6lg - %12.6lg - %12.6lg / %12.6lg ---> %8.6lg %8.6lg %8.6lg - %8.6lg\n",t0,ti0,dt,ti1,t1,x,y,z,l_powder);
printf("\nCoherent: %d, powder: %d\n",coherent,powder);
printf("%lg %lg %lg (%lg %lg) %lg %lg\n",t0,ti0,dt,ti1,t1,t21,t2i1);
printf("%6.1lf deg scattering angle (%lg, %d)\n", asin(q_v/(2.0*v))*2.0*RAD2DEG,q[i]*K2V,i);
printf("FATAL ERROR: Did not hit outer cylinder from inside.\n");
exit(1); /* This can really not happen ... at least, theoretically ... (TH 10/10/99) */
}
l_2 += v*(t21-t2i1);
if (powder==0)
{
my_s = V_my_s;
l_full=l_can;
}
else
{
l_full=l_powder;
if (coherent==1) my_s = my_s_v2/(v*v); /* my_s_v2 only known for coherent case */
else my_s=my_i; /* of course ... even if it does not sound logic at the first glimpse */
}
/* l_full is now the full path without scattering through the finally scattering material, powder or can (TH 10/10/99) */
my_total_powder=my_i+my_a_v/v+my_s_v2_total/(v*v);
my_total_can =V_my_s+V_my_a_v/v;
p *= fabs(l_full*my_s*exp(-my_total_powder*(l_p+l_1)-my_total_can*(l_c+l_2)))/(1.-(double)transmission);
/*printf("%10.2lg, %10.2lg, p=%10.2lg\n",my_total_powder,my_s_v2_total,p);*/
end_counter++;
/* if (p>1) exit(1); causes trouble with new source and weight>1*/
/*
if ((unsigned long)fmod(end_counter,50) == 0)
{
if (powder==0) printf("-");
if (coherent==1) printf("*");
else if (powder==1) printf("+");
fflush(stdout);
#if __dest_os == __mac_os
Event_loop();
#endif
}
*/
}
else ABSORB; /* ... if the neutron does not hit the sample ... */
%}
FINALLY
%{
fclose(outfile);
/*
printf("\n");
for (i=1;i<=nbInt;i++)
{
printf("i=%4d N=%8d <P>=%18.15f \n",i,Nsum[i],Psum[i]/Nsum[i]);
}
*/
%}
MCDISPLAY
%{
magnify("xyz");
circle("xz", 0, h/2.0, 0, radius);
circle("xz", 0, -h/2.0, 0, radius);
circle("xz", 0, h/2.0, 0, radius+d_V);
circle("xz", 0, -h/2.0, 0, radius+d_V);
line(-radius, -h/2.0, 0, -radius, +h/2.0, 0);
line(+radius, -h/2.0, 0, +radius, +h/2.0, 0);
line(0, -h/2.0, -radius, 0, +h/2.0, -radius);
line(0, -h/2.0, +radius, 0, +h/2.0, +radius);
line(-radius-d_V, -h/2.0, 0, -radius-d_V, +h/2.0, 0);
line(+radius+d_V, -h/2.0, 0, +radius+d_V, +h/2.0, 0);
line(0, -h/2.0, -radius-d_V, 0, +h/2.0, -radius-d_V);
line(0, -h/2.0, +radius+d_V, 0, +h/2.0, +radius+d_V);
%}
END
-------------- next part --------------
/*******************************************************************************
*
* McStas, version 1.2 released February 2000
* Maintained by Kristian Nielsen and Kim Lefmann,
* Risoe National Laboratory, Roskilde, Denmark
*
* %IDENTIFICATION
*
* Written by: <a href="mailto:hansen at ill.fr">Thomas C Hansen</a>
* Date: 07 March 2000
* Version: $Revision: 1.0 $
* Origin: <a href="http://www.ill.fr">ILL</a> (Dif/<a href="http://www.ill.fr/YellowBook/D20">D20</a>)
*
* Pyrolytic graphite filter
*
* %DESCRIPTION
*
* This pyrolytic graphite filter, defined by two identical rectangular opening apertures,
* is based upon empiric data and determines by
* interpolation the transmission of an incomming neutron, depending on its
* wavelength. This type of filter is e.g. used to supress higher harmonics, so
* the 1.2 AA contribution to the 2.4 AA obtained by a highly orientated
* pyrolytic graphite (HOPG) monochromator at a takeoff angle of 42 deg., as
* it is the case for the powder diffractometer D20 at the ILL Grenoble
*
* %PARAMETERS
*
* INPUT PARAMETERS:
*
* xmin: (m) Lower x bound (-0.16)
* xmax: (m) Upper x bound ( 0.16)
* ymin: (m) Lower y bound (-0.16)
* ymax: (m) Upper y bound ( 0.16)
* len: (m) Thickness of graphite plate ( 0.05)
*
* %LINKS
* <a href="../d20adapt.instr">Source code of d20adapt.instr</a>, where this component is used
*
* %END
*
*******************************************************************************/
DEFINE COMPONENT Filter_Graphite
DEFINITION PARAMETERS (xmin, xmax, ymin, ymax, len)
SETTING PARAMETERS ()
OUTPUT PARAMETERS ()
STATE PARAMETERS (x,y,z,vx,vy,vz,t,s1,s2,p)
DECLARE
%{
%}
INITIALIZE
%{
%}
TRACE
%{
double L,L1,L0,T1,T0, Filt_T;
double dt;
PROP_Z0;
L = (2*PI/V2K)/sqrt(vx*vx + vy*vy + vz*vz);
if (x<xmin || x>xmax || y<ymin || y>ymax) ABSORB;
dt = len/vz;
PROP_DT(dt);
if (x<xmin || x>xmax || y<ymin || y>ymax) ABSORB;
if (L>2.60) {L1= L ;L0=2.60;T1=0.06;T0=0.06;};
if (L<2.60) {L1=2.60;L0=2.40;T1=0.06;T0=0.05;};
if (L<2.40) {L1=2.40;L0=2.30;T1=0.05;T0=0.09;};
if (L<2.30) {L1=2.30;L0=2.20;T1=0.09;T0=0.22;};
if (L<2.20) {L1=2.20;L0=2.00;T1=0.22;T0=0.34;};
if (L<2.00) {L1=2.00;L0=1.90;T1=0.34;T0=0.61;};
if (L<1.90) {L1=1.90;L0=1.80;T1=0.61;T0=0.23;};
if (L<1.80) {L1=1.80;L0=1.60;T1=0.23;T0=0.25;};
if (L<1.60) {L1=1.60;L0=1.40;T1=0.25;T0=0.27;};
if (L<1.40) {L1=1.40;L0=1.30;T1=0.27;T0=0.53;};
if (L<1.30) {L1=1.30;L0=1.20;T1=0.53;T0=0.98;};
if (L<1.20) {L1=1.20;L0=1.10;T1=0.98;T0=0.89;};
if (L<1.10) {L1=1.10;L0=1.00;T1=0.89;T0=0.52;};
if (L<1.00) {L1=1.00;L0=0.87;T1=0.52;T0=0.47;};
if (L<0.87) {L1=0.87;L0= L ;T1=0.47;T0=0.47;};
Filt_T=(T0+(L-L0)/(L1-L0)*(T1-T0))*100.0;
Filt_T = exp(-Filt_T*len);
p*=Filt_T;
%}
FINALLY
%{
%}
MCDISPLAY
%{
magnify("xy");
multiline(5, (double)xmin, (double)ymin, 0.0,
(double)xmax, (double)ymin, 0.0,
(double)xmax, (double)ymax, 0.0,
(double)xmin, (double)ymax, 0.0,
(double)xmin, (double)ymin, 0.0);
multiline(5, (double)xmin, (double)ymin, (double)len,
(double)xmax, (double)ymin, (double)len,
(double)xmax, (double)ymax, (double)len,
(double)xmin, (double)ymax, (double)len,
(double)xmin, (double)ymin, (double)len);
line(xmin, ymin, 0.0, xmin, ymin, len);
line(xmax, ymin, 0.0, xmax, ymin, len);
line(xmin, ymax, 0.0, xmin, ymax, len);
line(xmax, ymax, 0.0, xmax, ymax, len);
%}
END
-------------- next part --------------
/*******************************************************************************
*
* McStas, version 1.2 released February 2000
* Maintained by Kristian Nielsen and Kim Lefmann,
* Risoe National Laboratory, Roskilde, Denmark
*
* %IDENTIFICATION
*
* Written by: <a href="mailto:hansen at ill.fr">Thomas C Hansen</a>
* Date: 08 March 2000
* Version: $Revision: 2.0 $
* Origin: <a href="http://www.ill.fr">ILL</a> (Dif/<a href="http://www.ill.fr/YellowBook/D20">D20</a>)
*
* Non-flat monochromator crystal with multiple scattering
*
* %DESCRIPTION
*
* Non-flat monochromator which uses a small-mosity approximation as well as
* the approximation vy^2 << vz^2 + vx^2. The crystal may be in transmission
* geometry as well as in reflection. Multiple scattering is possible. The
* individual reflectivity for each neutron is calculated from its velocity
* vector, the structure factor of the monochromator crystal and its mosaic.
* For an unrotated monochromator component, the crystal plane lies in the y-z
* plane (ie. parallel to the beam).
*
* %PARAMETERS
*
* INPUT PARAMETERS:
*
* xmin: (m) Lower x-bound of crystal
* xmax: (m) Upper x-bound of crystal
* zmin: (m) Lower z-bound of crystal
* zmax: (m) Upper z-bound of crystal
* ymin: (m) Lower y-bound of crystal
* ymax: (m) Upper y-bound of crystal
* mosh: (min) Horizontal mosaic (FWHM)
* mosv: (min) Vertical mosaic (FWHM)
* dist: (m) Distance from target (sample), not yet used (17.2)
* xw: (m) Width of target (sample), not yet used (0.01)
* write: (1) Flag (1/0) to write debugging information in a file
* F2: (fm**2) Structure factor of monochromator material
* Vc: (AA**3) Unit cell volume of monochromator material
* omega: (deg) Rocking angle difference from Bragg reflection position (0)
* Q: (1/AA) Wavevector of scattering
* harmonic: (1) Order of first higher, contamining harmonic (2 for Cu or HOPG, 3 for Ge)
* B: (1) Debye-Waller factor of monochromator material
* LIMIT: (1) Upper limit of multiple scattering events to avoid infinite loops (1001)
* blade: (1) Number of blade in a focussing monochromator of several blades (0)
* present: (1) Flag (1/0) of presence of the component in focussing monochromators with variable number of elements (1)
*
* OUTPUT PARAMETERS:
*
* %LINKS
* <a href="../d20adapt.instr">Source code of d20adapt.instr</a>, where this component is used
*
* %END
*
*******************************************************************************/
DEFINE COMPONENT Monochromator0
DEFINITION PARAMETERS (xmin, xmax, zmin, zmax, ymin, ymax, mosh, mosv, dist, xw, write, F2, Vc, omega, Q, harmonic, B, LIMIT,blade, present)
SETTING PARAMETERS ()
OUTPUT PARAMETERS ()
STATE PARAMETERS (x,y,z,vx,vy,vz,t,s1,s2,p)
DECLARE
%{
#define DIV_CUTOFF 2 /* ~ 10^-5 cutoff. */
double d0,my,tmp7,tmp8;
int z_edge,x_edge;
int z_edge0,x_edge0;
int counter;
%}
INITIALIZE
%{
/*
double my,r0,l,Mono_d,v;
*/
tmp7=0;
tmp8=0;
z_edge=0;
x_edge=0;
counter=0;
fprintf(monofile,"p0,r0, my0, l0, phi0,p1,r1,p2,vx2,vy2\n");
/*
v = sqrt(vx*vx+vy*vy+vz*vz);
Mono_d = 2*PI/Q;
my = PI*PI*PI*2.0*F2*exp(-Mono_B/4.0/Mono_d/Mono_d)/(Vc*Vc*V2K*V2K*Q)/v/v*(360*360*60*60/mosh/mosv/PI);
l=xmax-xmin;
if (l > (zmax-zmin)) l=zmax-zmin;
r0=1.0-exp(-my*l);
printf ("Monochromator: my=%lf/cm => R=%lf percent for l=%lfmm and v=%lfm/s\n",my/100,r0*100,l*1000,v);
*/
%}
TRACE
%{
double dphi,tmp1,tmp2,tmp3,tmp4,tmp5,tmp6,vratio,phi,theta0,theta,v,cs,sn;
double old_x = x, old_y = y, old_z = z, old_t = t, old_vx=vx,old_vz=vz, tmp_vx=vx,tmp_vz=vz;
double dt, dtx,dty,dtz,l, r0,q,l_total,tmp_x,tmp_y,tmp_z;
double new_x,new_y,new_z,new_t;
int j;
old_x=x; old_y=y;old_z=z;old_t=t;
if (present)
{
if (vx == 0.0) dtx=0.; else if ((xmin-x)/vx > (xmax-x)/vx) dtx = (xmax-x)/vx; else dtx = (xmin-x)/vx;
if (vy == 0.0) dty=0.; else if ((ymin-y)/vy > (ymax-y)/vy) dty = (ymax-y)/vy; else dty = (ymin-y)/vy;
if (vz == 0.0) dtz=0.; else if ((zmin-z)/vz > (zmax-z)/vz) dtz = (zmax-z)/vz; else dtz = (zmin-z)/vz;
if (dtx > dtz) dt=dtx; else dt=dtz;
if (dty > dt ) dt=dty;
new_x = x+vx*dt; new_y = y+vy*dt; new_z = z+vz*dt; new_t = t+dt;
if (new_x>=xmin-0.000001 && new_x<=xmax+0.000001 && new_z>=zmin-0.000001 && new_z<=zmax+0.000001 && new_y>=ymin-0.000001 && new_y<=ymax+0.000001)
{
x=new_x;y=new_y;z=new_z;t=new_t;
counter++;
if ((int)write == 1)
{
if (fabs(x-xmin)<.000001)
x_edge++;
if (fabs(z-zmin)<.000001)
z_edge++;
}
if (vx == 0.0) dtx=0; else if ((xmin-x)/vx < (xmax-x)/vx) dtx = (xmax-x)/vx; else dtx = (xmin-x)/vx;
if (vy == 0.0) dty=0; else if ((ymin-y)/vy < (ymax-y)/vy) dty = (ymax-y)/vy; else dty = (ymin-y)/vy;
if (vz == 0.0) dtz=0; else if ((zmin-z)/vz < (zmax-z)/vz) dtz = (zmax-z)/vz; else dtz = (zmin-z)/vz;
if (dtx < dtz) dt=dtx; else dt=dtz;
if (dty < dt ) dt=dty;
if (dt>0)
{
v = sqrt(vx*vx+vy*vy+vz*vz);
/* First: scattering in plane */
theta0 = atan2(vx,vz)-omega*PI/180; /* neutron angle to slab */
/* if(theta0 < 0) theta = -theta; - useless statement ... but, careful, why has it been written once? ... */
theta = asin(Q2V*Q/(2.0*v)); /* Bragg's law */
if (counter <0) printf("Theta0=%lf, Theta=%lf (%lf), Omega=%lf, ",theta0*180/PI,theta*180/PI,Q2V*Q/(2.0*v), omega);
tmp3 = (theta-theta0)/(MIN2RAD*mosh);
theta = (Q2V*Q*(double)harmonic/(2.0*v)); /* Bragg's law */
if (theta>1) theta=asin(1);
else theta = asin(theta); /* Bragg's law */
if (counter <0) printf("for %d.harmonic: Theta=%lf (%lf)",harmonic,theta*180/PI,Q2V*Q*(double)harmonic/(2.0*v));
if (fabs(tmp3)>fabs((theta-theta0)/(MIN2RAD*mosh))) q=Q*(double)harmonic;
else q=Q;
theta = asin(Q2V*q/(2.0*v)); /* Bragg's law */
tmp3 = (theta-theta0)/(MIN2RAD*mosh);
Mono_d = 2*PI/q;
if (counter <0) if (q!=Q) printf(" *");
if (counter <0) printf("\n");
l = v * dt;
l_total=l;
tmp4=exp(-tmp3*tmp3*4*log(2)); /*** the relative orientation goes into the calculation of my !!!(?) ******/
phi = atan2(vy,vz); /* ... but not the out-of plane angle! (? - because this will change with multiple scattering) */
my = tmp4*PI*PI*PI*2.0*F2*exp(-B/4.0/Mono_d/Mono_d)/(Vc*Vc*V2K*V2K*q)/v/v*(360*360*60*60/mosh/mosv/PI);
/*
my = tmp4*PI*PI*PI*2.0*F2*exp(-Mono_B/4.0/Mono_d/Mono_d)/(Vc*Vc*V2K*V2K*q)/v/v*(360*60/mosh/PI);
*/
/* Reflectivity r0 */
r0=1.0-exp(-my*l/**exp(-phi*phi*4*log(2)/(MIN2RAD*mosv)/(MIN2RAD*mosv))*/);
if ((counter <0)||(r0<0)) printf("F2=%lf, Vc=%lf, mosh=%lf, mosv=%lf\n",F2,Vc,mosh,mosv);
if ((counter <0)||(r0<0)) printf("v=%lfm/s, Q=%lf/AA\nMu=%lf/mm, path=%10.2lgmm, ",v,q,my/1000,l*1000);
if ((counter <0)||(r0<0)) printf("MuR=%lf, Reflectivity=%10.2lg, Transmission=%lg, time= %lf s\n",my*l,100*r0,100*(1-r0),dt);
d0 =1-r0;
d0+=rand01()*r0;
if ((counter <0)||(r0<0)) printf("Transmission at point of scattering: %lg, ",d0*100);
d0 =-log(d0);
d0/=my;
dt=d0/v;
if (counter <0) printf("\nEnter at x=%lgmm, z=%lgmm",x*1000,z*1000);
tmp_x=x; tmp_y=y; tmp_z=z;
x += vx*dt; y += vy*dt; z += vz*dt; t += dt;
l_total = sqrt((x-tmp_x)*(x-tmp_x)+(y-tmp_y)*(y-tmp_y)+(z-tmp_z)*(z-tmp_z));
fflush(monofile);
if (counter <0) printf("\nScattering at %lgmm, time=%lgs, x=%lgmm, z=%lgmm\n",d0*1000,dt,x*1000,z*1000);
if ((int)write == 1) if (counter < 1e4) fprintf(monofile,"%lf\n", p);
/* First: scattering in plane */
if(fabs(tmp3) > DIV_CUTOFF)
{
x = old_x; y = old_y; z = old_z; t = old_t;
}
else
{
SCATTER;
if (counter < 1e4) fprintf(monofile,"%12lg %12lg %12lg %12lg %12lg ",p,r0, my, l, phi);
if (r0<0) printf("p_in=%10.2lg, ",p);
incident_p+=p;
p *= r0;
first_p+=p;
if (counter < 1e4) fprintf(monofile,"%12lg ",p);
/* p *= tmp4; we consider this already in the calculation of my, don't we? */
tmp1 = 2*theta;
tmp6=1;
cs = cos(tmp1);
sn = sin(tmp1);
tmp2 = cs*vx - sn*vz;
if (counter <0) printf("vx=%lgm/s, vz=%lgm/s =>",vx,vz);
old_vz=vz; /***** NEW! *****/
old_vx=vx; /***** NEW! *****/
vy = vy;
vz = cs*vz + sn*vx;
vx = tmp2;
if (counter <0) printf("vx=%lgm/s, vz=%lgmm\n",vx,vz);
/* Second: scattering out of plane. **********************************/
/* Approximation is that Debye-Scherrer cone is a plane **************/
phi = atan2(vy,vz); /* out-of plane angle */
dphi = (MIN2RAD*mosv)/(2*sqrt(2*log(2)))*randnorm(); /* MC choice: **/
/* Vertical angle of the crystallite */
vy = vz*tan(phi+2*dphi*sin(theta));
vratio = v/sqrt(vx*vx+vy*vy+vz*vz);
vz = vz*vratio;
vy = vy*vratio; /* Renormalize v */
vx = vx*vratio;
/***************** NEW: MULTIPLE SCATTERING ************************************************************************/
j=1;
if (vx == 0.0) dtx=0; else if ((xmin-x)/vx < (xmax-x)/vx) dtx = (xmax-x)/vx; else dtx = (xmin-x)/vx;
if (vy == 0.0) dty=0; else if ((ymin-y)/vy < (ymax-y)/vy) dty = (ymax-y)/vy; else dty = (ymin-y)/vy;
if (vz == 0.0) dtz=0; else if ((zmin-z)/vz < (zmax-z)/vz) dtz = (zmax-z)/vz; else dtz = (zmin-z)/vz;
if (dtx < dtz) dt=dtx; else dt=dtz;
if (dty < dt ) dt=dty;
l = v * dt;
tmp5=rand01();
r0=1.0-exp(-my*l/**exp(-phi*phi*4*log(2)/(MIN2RAD*mosv)/(MIN2RAD*mosv))*/);
while (((tmp5<r0)||((j-2*(j/2))==0))&&(j<LIMIT))
{
j+=1;
d0 =1-r0;
if ((j-2*(j/2))==0) /* if neutron is about to leave in wrong direction, scatter it in any case and adjust probability ... */
{
p*=r0;
d0+=tmp5*r0;
}
else /* ... otherwise make a MC choice, also if limit has been passed */
{
d0=tmp5;
}
d0 =-log(d0);
d0/=my;
if (d0 > l) break;
dt=d0/v;
if (counter <0) printf("Scattering at %lgmm, time=%lgs",d0*1000,dt);
tmp_x=x; tmp_y=y; tmp_z=z;
x += vx*dt; y += vy*dt; z += vz*dt; t += dt;
SCATTER;
l_total+= sqrt((x-tmp_x)*(x-tmp_x)+(y-tmp_y)*(y-tmp_y)+(z-tmp_z)*(z-tmp_z));
if (counter <0) printf(", x=%lgmm, z=%lgmm\n",x*1000,z*1000);
tmp_vx=vx;
tmp_vz=vz;
vx=old_vx;
vz=old_vz;
old_vx=tmp_vx;
old_vz=tmp_vz;
phi = atan2(vy,vz); /* out-of plane angle */
dphi = (MIN2RAD*mosv)/(2*sqrt(2*log(2)))*randnorm(); /* MC choice: **/
/* Vertical angle of the crystallite */
vy = vz*tan(phi+2*dphi*sin(theta));
vratio = v/sqrt(vx*vx+vy*vy+vz*vz);
vz = vz*vratio;
vy = vy*vratio; /* Renormalize v */
vx = vx*vratio;
if (counter <0) printf("vx=%lgm/s, vz=%lgmm\n",vx,vz);
if (vx == 0.0) dtx=0; else if ((xmin-x)/vx < (xmax-x)/vx) dtx = (xmax-x)/vx; else dtx = (xmin-x)/vx;
if (vy == 0.0) dty=0; else if ((ymin-y)/vy < (ymax-y)/vy) dty = (ymax-y)/vy; else dty = (ymin-y)/vy;
if (vz == 0.0) dtz=0; else if ((zmin-z)/vz < (zmax-z)/vz) dtz = (zmax-z)/vz; else dtz = (zmin-z)/vz;
if (dtx < dtz) dt=dtx; else dt=dtz;
if (dty < dt ) dt=dty;
l = v * dt;
r0=1.0-exp(-my*l/**exp(-phi*phi*4*log(2)/(MIN2RAD*mosv)/(MIN2RAD*mosv))*/);
tmp5=rand01();
}
if (((vx*theta)>0)||(vx>0))
{
if (counter < 1e4) fprintf(monofile,"%12lg %12lg %12lg %12lg\n",r0,p,vx,vz);
leaving_p+=p;
ABSORB;
}
/****************** END OF MULTIPLE SCATTERING ***********************/
if (p<0) printf("r0=%10.2lg, tmp3=%10.2lg, p_out=%10.2lg\n",r0,tmp3,p);
leaving_p+=p;
if (counter < 1e4) fprintf(monofile,"%12lg %12lg %12lg %12lg\n",r0,p,vx,vz);
fflush(monofile);
}
}
else
{
/* x = old_x; y = old_y; z = old_z; t = old_t; */
}
}
else
{
/* x = old_x; y = old_y; z = old_z; t = old_t; */
}
if (counter<0) counter=0;
} /* if present */
%}
FINALLY
%{
if ((int)write == 1) printf("x edge : %d neutrons, z edge : %d neutrons\n",x_edge ,z_edge );
%}
MCDISPLAY
%{
magnify("xyz");
multiline(16,(double)xmin*present, (double)ymin*present, (double)zmin*present,
(double)xmin*present, (double)ymax*present, (double)zmin*present,
(double)xmin*present, (double)ymax*present, (double)zmax*present,
(double)xmin*present, (double)ymin*present, (double)zmax*present,
(double)xmin*present, (double)ymin*present, (double)zmin*present,
(double)xmax*present, (double)ymin*present, (double)zmin*present,
(double)xmax*present, (double)ymax*present, (double)zmin*present,
(double)xmin*present, (double)ymax*present, (double)zmin*present,
(double)xmax*present, (double)ymax*present, (double)zmin*present,
(double)xmax*present, (double)ymax*present, (double)zmax*present,
(double)xmin*present, (double)ymax*present, (double)zmax*present,
(double)xmax*present, (double)ymax*present, (double)zmax*present,
(double)xmax*present, (double)ymin*present, (double)zmax*present,
(double)xmin*present, (double)ymin*present, (double)zmax*present,
(double)xmax*present, (double)ymin*present, (double)zmax*present,
(double)xmax*present, (double)ymin*present, (double)zmin*present);
%}
END
-------------- next part --------------
/***********************************************************************
*
* McStas, version 1.2 released February 2000
* Maintained by Kristian Nielsen and Kim Lefmann,
* Risoe National Laboratory, Roskilde, Denmark
*
* %IDENTIFICATION
*
* Author: <a href="mailto:hansen at ill.fr">Thomas C Hansen</a>
* Date: 07 March 2000
* Version: $Revision: 1.0 $
* Origin: <a href="http://www.ill.fr">ILL</a> (Dif/<a href="http://www.ill.fr/YellowBook/D20">D20</a>)
*
* Curved linear 1D MSGC PSD
*
* %DESCRIPTION
*
* A curved linear 1D PSD monitor using a cylindrical projection. This detector type is commonly
* used in constant wavelength neutron powder diffraction, based either on multiwire- (MWGC) or
* microstrip- (MSGC) gaschamber technology. This implementation is mainly pushed by the MSGC
* realisation (D20) of such a PSD, as done on D20 at ILL. So some kind of 'polygonality' effect
* resulting from the polygonal arrangement of MSGC plates will be considered in the near future
* by an approximation of the electron avalanche trace in the electrical field.
*
* %PARAMETERS
*
* INPUT PARAMETERS:
*
* radius: (m) Radius of detector at center of MSGC plate (=radius in <a href="PSD_entry.html">PSD_entry</a> + gap) (1.471)
* height: (m) Height of detector (=height in <a href="PSD_entry.html">PSD_entry</a>) (0.15)
* nd: (1) Number of cells (=nd in <a href="PSD_entry.html">PSD_entry</a>) (1536)
* pitch: (deg) angular pitch (=pitch in <a href="PSD_entry.html">PSD_entry</a>) (0.1)
* gap: (m) gaschamber detection gap (0.053)
* filename: (string) Name of file in which to store the detector image (NULL)
* sign: (1) Chirality of 1st diffractometer axis=sign[takeoff of <a href="../d20adapt.instr">d20adapt.instr]</a> (-1)
* tt0: (deg) Angular position of 1st PSD cell (=tt0 in <a href="PSD_entry.html">PSD_entry</a>) (0)
* pdet: (bar) Pressure of 3He detection gas (1.2)
* pT_trace: (m) Average effective length of proton/tritium trace (0.005)
* period: (1) Number of detection cells per MSGC plate (=period in <a href="PSD_entry.html">PSD_entry</a>) (32)
* cellwidth: (m) Width of one detection cell on a MSGC plate (0.002568)
*
* OUTPUT PARAMETERS:
*
* PSD_N: (*1) Array of neutron counts
* PSD_p: (*1.0) Array of neutron weight counts
* PSD_p2: (*1.0) Array of second moments
* alpha: (rad) Angle covered by one MSGC plate
* R1: (m) PSD radius at MSGC edge
*
* %LINKS
* <a href="../d20adapt.instr">Source code of d20adapt.instr</a>, where this component is used
* %LINKS
* <a href="PSD_entry.html">Source code of PSD_entry.comp</a>, a corresponding entry window
*
* %END
*
***********************************************************************/
DEFINE COMPONENT PSD_curved
DEFINITION PARAMETERS ( radius, height, nd, pitch, gap, filename, sign, tt0, pdet, pT_trace, period, cellwidth)
SETTING PARAMETERS ()
OUTPUT PARAMETERS (PSD_N, PSD_p, PSD_p2, alpha, R1)
STATE PARAMETERS (x,y,z,vx,vy,vz,t,s1,s2,p)
DECLARE
%{
int PSD_N[3601];
double PSD_p[3601];
double PSD_p2[3601];
double alpha, R1;
%}
INITIALIZE
%{
int i;
alpha=period*pitch/360.0*PI;
R1=radius/cos(alpha);
for (i=0; i<nd; i++)
{
PSD_N [i] = 0;
PSD_p [i] = 0;
PSD_p2[i] = 0;
}
%}
TRACE
%{
double l,phi,v, l_full, t0, dt, t1, t2,t3, twotheta,eff,p_detection;
double psi,theta, ax,ay,az, tx,ty,tz, ox,oy,oz, bx,by,bz;
double x0,z0,R2,x1,z1,twotheta1,twotheta0,xp0,zp0,xp1,zp1,xx,zz,tmp;
int i0,i1,i,plate0,plate1;
if(cylinder_intersect(&t0, &t1, x, y, z, vx, vy, vz, radius-gap, height) /* && t1 > 0 */)
{
/*if(t0 < 0) t0 = t1;*/
if(t0 < t1) t0 = t1;
/*PROP_DT(t0);*/
x+=vx*t0;
y+=vy*t0;
z+=vz*t0;
t+=t0;
SCATTER; /* incomming neutron */
/* if ((double)sign*twotheta/PI*180. < tt0-pitch) ABSORB; */
/* if ((double)sign*twotheta/PI*180. > tt0+(nd+1)*pitch) ABSORB; */
/* Polygonality effect in MSGC PSD - the detection gap is not really constant ... */
/* Where would the n hit, if the outer radius is where the plates are closest to the entry */
cylinder_intersect(&t2, &t3, x, y, z, vx, vy, vz, radius, height);
if(t2 < 0) t2 = t3;
x0=x+t3*vx;
z0=z+t3*vz;
twotheta0 = -atan2(x0,z0);
i0 = floor(0.5+((double)sign*twotheta0/PI*180.0-tt0)/pitch);
plate0=i0/period;
/*
printf("\nplate: %d, cell %d, twotheta %lf\n",plate0,i0,twotheta0/PI*180.0);
*/
/* Where would the n hit, if the outer radius is where two plates join */
cylinder_intersect(&t2, &t3, x, y, z, vx, vy, vz, R1, height);
if(t2 < 0) t2 = t3;
x1=x+t3*vx;
z1=z+t3*vz;
twotheta1 = -atan2(x1,z1);
i1 = floor(0.5+((double)sign*twotheta1/PI*180.0-tt0)/pitch);
plate1=i1/period;
/*
printf("plate: %d, cell %d, twotheta %lf\n",plate1,i1,twotheta1/PI*180.0);
*/
/* We try now to get the radius R2 where the n hits the plate */
twotheta=(tt0-pitch/2.0+(plate0+0.5)*period*pitch)*PI/180.0;
/*
printf("twotheta of plate %d: %lf\n",plate0,twotheta*180.0/PI);
*/
zp0=cos(twotheta-alpha)*R1;
xp0=sin(twotheta-alpha)*R1;
zp1=cos(twotheta+alpha)*R1;
xp1=sin(twotheta+alpha)*R1;
xx = (z0-(z1 -z0 )/(x1 -x0 )*x0) - (zp0-(zp1-zp0)/(xp1-xp0)*xp0);
xx/= (zp1-zp0)/(xp1-xp0) - (z1 -z0 )/(x1 -x0 );
zz = (z0-(z1 -z0 )/(x1 -x0 )*x0) + (z1 -z0 )/(x1 -x0 ) * xx ;
/*
printf("mt=%lf;bt=%lf \n",(z1 -z0 )/(x1 -x0 ),(z0 -(z1 -z0 )/(x1 -x0 )*x0 ));
printf("mp=%lf;bp=%lf \n",(zp1-zp0)/(xp1-xp0),(zp0-(zp1-zp0)/(xp1-xp0)*xp0));
printf("xx[0]=%lf;zz[0]=%lf\n",xx,zz);
*/
R2=sqrt(xx*xx+zz*zz);
/*
printf("R0=%lf;R1=%lf;R2=%lf\n",radius,R1,R2);
printf("neutron hits plate %d (radius %lf < %lf < %lf) %lf %lf\n",plate0,radius,R2,R1,sqrt(x0*x0+z0*z0),sqrt(x1*x1+z1*z1));
printf("twotheta = %lf to %lf\n",twotheta0/PI*180.0,twotheta1/PI*180.0);
printf("plate %d (%lf,%lf) to (%lf,%lf)\n",plate0,xp0,zp0,xp1,zp1);
printf("n hits between (%lf,%lf) and (%lf,%lf)\n",x0,z0,x1,z1);
printf("%lf of plate \n",b);
*/
if (plate1 != plate0)
{
twotheta=(tt0-pitch/2.0+(plate1+1)*period*pitch)*PI/180.0;
zp0=cos(twotheta-alpha)*R1;
xp0=sin(twotheta-alpha)*R1;
zp1=cos(twotheta+alpha)*R1;
xp1=sin(twotheta+alpha)*R1;
xx = (z0-(z1 -z0 )/(x1 -x0 )*x0) - (zp0-(zp1-zp0)/(xp1-xp0)*xp0);
xx/= (zp1-zp0)/(xp1-xp0) - (z1 -z0 )/(x1 -x0 );
zz = (z0-(z1 -z0 )/(x1 -x0 )*x0) + (z1 -z0 )/(x1 -x0 ) * xx ;
tmp=sqrt(xx*xx+zz*zz);
if (tmp<R2)
{
R2=tmp;
/*
printf("neutron hits plate %d (radius %lf < %lf < %lf)\n",plate1,radius,R2,R1);
*/
}
/*
else printf("neutron hits plate %d (radius %lf < %lf < %lf)\n",plate0,radius,R2,R1);
*/
}
if (R2 > R1)
{
/*
printf("neutron hits plate %d or %d (radius %lf < %lf < %lf)\n",plate0,plate1,radius,R2,R1);
printf("##############################################################\n");
*/
R2=R1;
}
/* Now we got the virtual outer radius R2 which is in between 'radius' and 'R1' */
cylinder_intersect(&t2, &t3, x, y, z, vx, vy, vz, R2, height);
if(t2 < 0) t2 = t3;
/* The following thing is not working as exspected to trace the trajectory of a neutron through the detection gap
PROP_DT(t2);
SCATTER;
PROP_DT(-t2);
SCATTER;
*/
x+=vx*t2;
y+=vy*t2;
z+=vz*t2;
SCATTER;
x-=vx*t2;
y-=vy*t2;
z-=vz*t2;
SCATTER;
v = sqrt(vx*vx + vy*vy + vz*vz);
l_full = v * (t3); /* Length of full path through sample */
/*
printf("Trajectory in gap %lfmm (%lf microsecs), ",l_full*1000,t3*1000000);
*/
eff=1.0-exp(-0.07417*pdet*l_full*100.0*VL/v);
dt = (t3) * -log(rand01()*(eff)+1.0-eff);/* + t0*/
/*
printf("efficiency %lf percent, capture after %lfmm\n",eff*100,dt*v*1000);
*/
PROP_DT(dt);
SCATTER; /* Point of capture */
l = v*dt; /* Penetration in detector */
p_detection=p*eff;
p*=1-eff; /* a not detected (=captured) neutron may continue the flightpath!) */
/* Now, 2Theta for the place, where the neutron should have been detected ... */
twotheta = -atan2(x,z);
/* target vector b from neutron to PSD (in plane: y=0) */
bx= sin(twotheta)*radius-x;
by= 0.0;
bz= cos(twotheta)*radius-z;
/* axis a perpendicular to neutron vector v */
vec_prod(ax,ay,az, vx,vy,vz, bx,by,bz);
/* MC choice of angle component psi in a circle */
psi =2*PI*rand01();
/* rotation in circle of neutron vector v around axis a (in plane) -> temporary trace vector t */
rotate (tx,ty,tz, vx,vy,vz, psi, ax,ay,az);
/* MC choice of second angle component in a sphere */
theta=acos(randpm1());
/* rotation in sphere of temporary trace vector t (out of plane) -> output pT trace vector o */
rotate (ox,oy,oz, tx,ty,tz, theta, vx,vy,vz);
/* vx=ox;vy=oy;vz=oz; - NOT the neutron changes it's direction but the pT trace o */
NORM(ox,oy,oz);
/* neutron is virtually detected at about a quarter of the total p/T trace length from capture */
ox*=pT_trace/4.0;
oy*=pT_trace/4.0;
oz*=pT_trace/4.0;
x+=ox;y+=oy;z+=oz;
/* SHOW the part of pT trace leading to the point of virtual neutron capture detcted by the cas chamber PSD */
SCATTER;
/* NOW 2Theta and phi is calculated for the gravicenter of electron avalanche */
twotheta = -atan2(x,z);
phi = asin(y/radius);
/* i = floor(0.5+(nx-1)*( (double)sign*twotheta/(2*PI) - tt0/360.)); */
i = floor(0.5+((double)sign*twotheta/PI*180.0-tt0)/pitch);
if ((i<nd && i>=0))
{
PSD_N [i] ++;
PSD_p [i] += p_detection;
PSD_p2[i] += p_detection*p_detection;
}
else
{
ABSORB;
}
}
else
{
ABSORB;
}
/* Yes, indeed, ABSORB IS a good idea for non-detected neutrons, as they do not interest anymore, and to use Check_adapt afterwards! */
%}
FINALLY
%{
int i;
double total=0.0,histories=0.0;
char string[40];
for (i=0;i<nd;i++)
{
total+=PSD_p[i];
histories+=PSD_N[i];
}
sprintf(string,"Curved lin. PSD %9.4lgn/s (%8.3lg events)",total,histories);
DETECTOR_OUT_1D(
string,
"2theta [deg]",
"Intensity",
"twotheta", tt0, tt0+nd*pitch, nd,
&PSD_N[0],&PSD_p[0],&PSD_p2[0],
filename);
%}
MCDISPLAY
%{
double R2,x0,x1,z1,z0,twotheta,beta;
int plate,cell;
magnify("xyz");
R1=radius/cos(period*pitch/360.0*PI);
for (plate=0;plate<ceil((double)nd/(double)period);plate++)
{
twotheta=(tt0-pitch/2.0+(plate+0.5)*period*pitch)*PI/180.0;
beta=atan(period/2.0*cellwidth/radius);
R2=radius/cos(beta);
z0=cos(twotheta-beta)*R2;
x0=sin(twotheta-beta)*R2;
z1=cos(twotheta+beta)*R2;
x1=sin(twotheta+beta)*R2;
line(x0,-height/2.0, z0,x1,-height/2.0, z1);
line(x0,+height/2.0, z0,x1,+height/2.0, z1);
line(x0,-height/2.0, z0,x0,+height/2.0, z0);
line(x1,-height/2.0, z1,x1,+height/2.0, z1);
for (cell=0;cell<period;cell++)
{
circle("xz",x0+((double)cell+0.5)*(x1-x0)/(double)period,0,z0+((double)cell+0.5)*(z1-z0)/(double)period,cellwidth/2.0);
}
}
R2=radius-gap;
twotheta=(tt0-pitch/2.0)*PI/180.0;
z0=cos(twotheta)*R2;
x0=sin(twotheta)*R2;
line(x0, -height/2.0,z0,x0, height/2.0,z0);
for (cell=1;cell<=ceil((double)nd/(double)period)*period;cell++)
{
twotheta=(tt0-pitch/2.0+pitch*cell)*PI/180.0;
z1=cos(twotheta)*R2;
x1=sin(twotheta)*R2;
line(x0, -height/2.0,z0,x1, -height/2.0,z1);
line(x0, height/2.0,z0,x1, height/2.0,z1);
/* printf("%4d %10.5lf %10.5lf %10.5lf %10.5lf\n",cell ,x0, z0,x1,z1);*/
x0=x1; z0=z1;
}
line(x0, -height/2.0,z0,x0, height/2.0,z0);
%}
END
-------------- next part --------------
/***********************************************************************
*
* McStas, version 1.2 released February 2000
* Maintained by Kristian Nielsen and Kim Lefmann,
* Risoe National Laboratory, Roskilde, Denmark
*
* %IDENTIFICATION
*
* Author: <a href="mailto:hansen at ill.fr">Thomas C Hansen</a>
* Date: 07 March 2000
* Version: $Revision: 1.1 $
* Origin: <a href="http://www.ill.fr">ILL</a> (Dif/<a href="http://www.ill.fr/YellowBook/D20">D20</a>)
*
* Entry window of a curved linear gaschamber PSD
*
* %DESCRIPTION
*
* Aluminium entry window for a curved linear position sensitive detector (PSD),
* as used in constant wavelength neutron powder diffraction (ILL Grenoble: D20).
* Normally, this component is always used in front of a <a href="PSD_curved.html">PSD_curved</a> component.
*
*
* %PARAMETERS
*
* INPUT PARAMETERS:
*
* radius: (m) Outer radius of detector entry window (=radius-gap in <a href="PSD_curved.html">PSD_curved</a>) (1.418)
* height: (m) Height of detector (=height in <a href="PSD_curved.html">PSD_curved</a>) (0.15)
* nd: (1) Number of cells (=nd in <a href="PSD_curved.html">PSD_curved</a>) (1536)
* pitch: (deg) angular pitch (=pitch in <a href="PSD_curved.html">PSD_curved</a>) (0.1)
* gap: (m) Thickness of entry window (0.007)
* tt0: (deg) detection angle of first cell (=tt0 in <a href="PSD_curved.html">PSD_curved</a>) (0)
* pack: (1) Packing density of window material
* Vc: (AA**3) Unit cell volume of window material
* sigma_a: (barn) Absorption cross section
* q: (1/AA) Wavevector
* j: (1) Multiplicity
* F2: (fm**2) Structure factor
* DW: (1/AA**2) Debye-Waller factor
* period: (1) Number of cells on one MSGC plate (=period in <a href="PSD_curved.html">PSD_curved</a>) (32)
* sign: (1) Chirality of 1st diffractometer axis=sign[takeoff of <a href="../d20adapt.instr">d20adapt.instr]</a> (-1)
*
* OUTPUT PARAMETERS:
*
* my_s_v2: (m/s**2) Attenuation factor due to elastic coherent scattering, multiplied by neutron velocity**2
* my_a_v: (1/s) Attenuation factor due to absorption, multiplied by neutron velocity
* q_v: (m/s) Corresponding velocity of wavevector Q
*
* %LINKS
* <a href="../d20adapt.instr">Source code of d20adapt.instr</a>, where this component is used
* %LINKS
* The corresponding detector <a href="PSD_curved.html">PSD_curved.comp</<a>
*
* %END
*
***********************************************************************/
DEFINE COMPONENT PSD_entry
DEFINITION PARAMETERS(radius,height,nd,pitch,gap,tt0,pack,Vc,sigma_a,q,j,F2,DW,period,sign)
SETTING PARAMETERS ()
OUTPUT PARAMETERS (my_s_v2, my_a_v, q_v)
STATE PARAMETERS (x,y,z,vx,vy,vz,t,s1,s2,p)
DECLARE
%{
double my_s_v2, my_a_v, q_v;
%}
INITIALIZE
%{
my_a_v = sigma_a/Vc*2200; /* Is not yet divided by v */
my_s_v2 = PI*PI*PI*pack*j*F2*DW/(Vc*Vc*V2K*V2K*q);
q_v = q*K2V;
%}
TRACE
%{
double t0, t1, t2, t3, twotheta, d_phi, phi,eff,p_transmission;
double my, l_full, dt, aim_x, aim_y, aim_z, axis_x, axis_y, axis_z;
double l, tmp_vx, tmp_vy, tmp_vz, vout_x, vout_y, vout_z,v;
int i;
if(cylinder_intersect(&t0, &t1, x, y, z, vx, vy, vz, radius-gap, height) && t1 > 0)
{
if(t0 < 0) t0 = t1;
PROP_DT(t0);
/*SCATTER;*/
twotheta = -atan2(x,z);
/*printf("%lf\n",twotheta);*/
if (((double)sign*twotheta>=tt0-pitch/2.0) && ((double)sign*twotheta<=tt0+(double)nd*pitch+pitch/2.0))
{
/*SCATTER;*/
cylinder_intersect(&t2, &t3, x, y, z, vx, vy, vz, radius, height);
if(t2 < 0) t2 = t3;
v = sqrt(vx*vx + vy*vy + vz*vz);
l_full = v * (t2); /* Length of full path through PSD entry window of Al */
my=my_s_v2/v/v+my_a_v/v;
p_transmission=exp(-my*l_full);
if (rand01() > p_transmission)
{
eff=1.0-exp(-my_s_v2/v/v*l_full);
dt =(t2)*-log(rand01()*(eff)+1.0-eff);
PROP_DT(dt); /* Point of diffraction */
l = v*dt; /* Penetration in entry window */
p*=(eff)/(1-p_transmission); /* a not detected (=captured) neutron may continue the flightpath!) */
d_phi = 360.0*DEG2RAD/2.0*randpm1();
twotheta = 2.0*asin(q_v/(2.0*v)); /* Bragg scattering law */
aim_x = sin(twotheta)-x; /* Vector pointing at target (anal./det.) */
aim_y = -y ;/* + PSD_h/2*randpm1() ; doppelt wg. d_phi?*/
aim_z = cos(twotheta)-z;
vec_prod(axis_x, axis_y, axis_z, vx, vy, vz, aim_x, aim_y, aim_z);
rotate(tmp_vx, tmp_vy, tmp_vz, vx, vy, vz, twotheta, axis_x, axis_y, axis_z);
rotate(vout_x, vout_y, vout_z, tmp_vx, tmp_vy, tmp_vz, d_phi, vx, vy, vz);
vx = vout_x;
vy = vout_y;
vz = vout_z;
}
else
{
PROP_DT(t2);
/*SCATTER;*/
}
}
else ABSORB;
}
else ABSORB;
/* Yes, ABSORB is a good idea for non-detectedable neutrons, as they do not interest anymore, and to use Check_adapt afterwards! */
%}
FINALLY
%{
%}
MCDISPLAY
%{
double R2, twotheta,x0,z0,x1,z1;
int cell;
magnify("xyz");
R2=radius-gap;
twotheta=(tt0-pitch/2.0)*PI/180.0;
z0=cos(twotheta)*R2;
x0=sin(twotheta)*R2;
line(x0, -height/2.0,z0,x0, height/2.0,z0);
for (cell=1;cell<=ceil((double)nd/(double)period)*period;cell++)
{
twotheta=(tt0-pitch/2.0+pitch*cell)*PI/180.0;
z1=cos(twotheta)*R2;
x1=sin(twotheta)*R2;
line(x0, -height/2.0,z0,x1, -height/2.0,z1);
line(x0, height/2.0,z0,x1, height/2.0,z1);
x0=x1; z0=z1;
}
line(x0, -height/2.0,z0,x0, height/2.0,z0);
/* Circles are too roughly drawn for curved detector presentations */
/*
circle("xz", 0, height/2.0, 0, radius);
circle("xz", 0, -height/2.0, 0, radius);
line(-radius, -height/2.0, 0, -radius, +height/2.0, 0);
line(+radius, -height/2.0, 0, +radius, +height/2.0, 0);
line(0, -height/2.0, -radius, 0, +height/2.0, -radius);
line(0, -height/2.0, +radius, 0, +height/2.0, +radius);
circle("xz", 0, height/2.0, 0, (radius-gap));
circle("xz", 0, -height/2.0, 0, (radius-gap));
line(-(radius-gap), -height/2.0, 0, -(radius-gap), +height/2.0, 0);
line(+(radius-gap), -height/2.0, 0, +(radius-gap), +height/2.0, 0);
line(0, -height/2.0, -(radius-gap), 0, +height/2.0, -(radius-gap));
line(0, -height/2.0, +(radius-gap), 0, +height/2.0, +(radius-gap));
*/
%}
END
-------------- next part --------------
/***********************************************************************
*
* McStas, version 1.2 released February 2000
* Maintained by Kristian Nielsen and Kim Lefmann,
* Risoe National Laboratory, Roskilde, Denmark
*
* %IDENTIFICATION
*
* Author: <a href="mailto:hansen at ill.fr">Thomas C Hansen</a>
* Date: 07 March 2000
* Version: $Revision: 1.0 $
* Origin: <a href="http://www.ill.fr">ILL</a> (Dif/<a href="http://www.ill.fr/YellowBook/D20">D20</a>)
*
* Conversion from flat energy distribution to Maxwellian neutron spectra
*
* %DESCRIPTION
*
* The routine modifies/converts the incident flat (in energy) neutron spectra,
* as it comes out from a source like the standard adaptive source
* <a href="http://neutron.risoe.dk/mcstas/mcdoc/components/sources/Source_adapt.html">Source_adapt</a>,
* to obtain a Maxwellian neutron distribution, as it is approximatively the case for thermal neutrons
* in a reactor's moderator (e.g. the HFR of the ILL in Grenoble). So, this
* component is only a virtual one, not corresponding to a real element in a
* neutron beam.
*
* %PARAMETERS
*
* INPUT PARAMETERS:
*
* T: (K) The temperature of the source, 293 K for thermal neutrons (beam tube H11 at ILL)
* E0: (meV) Mean energy of neutrons (as for <a href="http://neutron.risoe.dk/mcstas/mcdoc/components/sources/Source_adapt.html">Source_adapt</a>)
* dE: (meV) Energy spread as for <a href="http://neutron.risoe.dk/mcstas/mcdoc/components/sources/Source_adapt.html">Source_adapt</a> (energy range is from E0-dE to E0+dE)
*
* OUTPUT PARAMETERS:
*
* Lmin: (AA) Minimum wavelength in energy range
* Lmax: (AA) Maxium wavelength in energy range
* l0: (AA) Mean wavelength of neutron spectra, about 1.8 AA for thermal neutrons
* L2P: (AA**4) Normalisation constant for Maxwell distribution
*
* %LINKS
* <a href="../d20adapt.instr">Source code of d20adapt.instr</a>, where this component is used
* %LINKS
* A source producing a flat in energy distribution: <a href="http://neutron.risoe.dk/mcstas/mcdoc/components/sources/Source_adapt.html">Source_adapt.comp</a>
*
* %END
*
***********************************************************************/
DEFINE COMPONENT Convert_FlatE_2_Maxwell
DEFINITION PARAMETERS (T,E0,dE)
SETTING PARAMETERS ()
OUTPUT PARAMETERS (Lmin,Lmax,l0,L2P)
STATE PARAMETERS (x,y,z,vx,vy,vz,t,s1,s2,p)
DECLARE
%{
double l0, L2P, Lmin,Lmax;
%}
INITIALIZE
%{
double k;
k = 1.38066e-23;
l0=1.0e10*sqrt(HBAR*HBAR*4.0*PI*PI/2.0/MNEUTRON/k/T);
L2P=2*l0*l0*l0*l0;
Lmin = sqrt(81.81/(E0+dE)); /* AAngstroem */
Lmax = sqrt(81.81/(E0-dE));
%}
TRACE
%{
double l,E, Maxwell,dE_dl;
l= (2*PI/V2K)/sqrt(vx*vx + vy*vy + vz*vz);
E=(VS2E*4.0*(PI*PI)/(V2K*V2K))/l/l;
dE_dl=E/l*2.0;
p*=2.0*dE/(Lmax - Lmin)/dE_dl;
Maxwell=L2P/l/l/l/l/l*exp(-l0*l0/l/l);
p*= Maxwell;
%}
FINALLY
%{
%}
MCDISPLAY
%{
magnify("xy");
circle("xy",0,0,0,0.1);
%}
END
-------------- next part --------------
/*******************************************************************************
*
* McStas, version 1.1, released
* Maintained by Kristian Nielsen and Kim Lefmann,
* Risoe National Laboratory, Roskilde, Denmark
*
* Component: Hdiv_monitor
*
* Written by: KL, Nov. 11, 1998
* Written by: Thomas Hhansen, December, 1999
*
* A divergence sensitive monitor. The counts are distributed in
* n pixels.
*
* INPUT PARAMETERS:
*
* xmin: Lower x bound of detector opening (m)
* xmax: Upper x bound of detector opening (m)
* ymin: Lower y bound of detector opening (m)
* ymax: Upper y bound of detector opening (m)
* nh: Number of pixel rows (1)
* h_maxdiv Maximal vertical divergence detected (degrees)
* filename: Name of file in which to store the detector image (text)
*
* OUTPUT PARAMETERS:
*
* Div_N: Array of neutron counts
* Div_p: Array of neutron weight counts
* Div_p2: Array of second moments
*
*******************************************************************************/
DEFINE COMPONENT Hdiv_monitor
DEFINITION PARAMETERS (xmin, xmax, ymin, ymax,
nh, h_maxdiv, filename)
SETTING PARAMETERS ()
OUTPUT PARAMETERS (Div_N, Div_p, Div_p2)
STATE PARAMETERS (x,y,z,vx,vy,vz,t,s1,s2,p)
DECLARE
%{
int Div_N[nh];
double Div_p[nh];
double Div_p2[nh];
%}
INITIALIZE
%{
int i;
for (i=0; i<nh; i++)
{
/* printf("HDiv_monitor: %d\n",i); */
Div_N[i] = 0;
Div_p[i] = 0;
Div_p2[i] = 0;
}
/* printf("%d %d %d\n",i,nh,h_maxdiv); */
%}
TRACE
%{
int i;
double h_div;
PROP_Z0;
if (x>xmin && x<xmax && y>ymin && y<ymax)
{
h_div = RAD2DEG*atan2(vx,vz);
if (h_div < (double)h_maxdiv && h_div > -(double)h_maxdiv)
{
i = floor((h_div + (double)h_maxdiv)*nh/(2.0*(double)h_maxdiv));
Div_N[i]++;
Div_p[i] += p;
Div_p2[i] += p*p;
}
}
%}
FINALLY
%{
DETECTOR_OUT_1D(
"horizontal divergence monitor",
"horizontal divergence [deg]",
"Intensity",
"divergence", -h_maxdiv, h_maxdiv, nh,
&Div_N[0],&Div_p[0],&Div_p2[0],
filename);
%}
MCDISPLAY
%{
magnify("xy");
multiline(5, (double)xmin, (double)ymin, 0.0,
(double)xmax, (double)ymin, 0.0,
(double)xmax, (double)ymax, 0.0,
(double)xmin, (double)ymax, 0.0,
(double)xmin, (double)ymin, 0.0);
%}
END
-------------- next part --------------
/*******************************************************************************
*
* McStas, version 1.0, released October 26, 1998
* Maintained by Kristian Nielsen and Kim Lefmann,
* Risoe National Laboratory, Roskilde, Denmark
*
* Component: PSD_bidim
*
* Written by: KL, Feb 3, 1998 (PSD_monitor.comp)
* Rewritten by: Thomas Hansen, ILL Grenoble, France, 1999
* Output modified
*
* An (n times m) pixel PSD monitor. This component may also be used as a beam
* detector.
*
* INPUT PARAMETERS:
*
* xmin: Lower x bound of detector opening (m)
* xmax: Upper x bound of detector opening (m)
* ymin: Lower y bound of detector opening (m)
* ymax: Upper y bound of detector opening (m)
* nx: Number of pixel columns (1)
* ny: Number of pixel rows (1)
* filename: Name of file in which to store the detector image (text)
*
* OUTPUT PARAMETERS:
*
* PSD_N: Array of neutron counts
* PSD_p: Array of neutron weight counts
* PSD_p2: Array of second moments
*
*******************************************************************************/
DEFINE COMPONENT PSD_bidim
DEFINITION PARAMETERS (xmin, xmax, ymin, ymax, nx, ny, filename)
SETTING PARAMETERS ()
OUTPUT PARAMETERS (PSD_N, PSD_p, PSD_p2)
STATE PARAMETERS (x,y,z,vx,vy,vz,t,s1,s2,p)
DECLARE
%{
int PSD_N[nx][ny];
double PSD_p[nx][ny];
double PSD_p2[nx][ny];
%}
INITIALIZE
%{
int i,j;
for (i=0; i<nx; i++)
for (j=0; j<ny; j++)
{
PSD_N[i][j] = 0;
PSD_p[i][j] = 0;
PSD_p2[i][j] = 0;
}
%}
TRACE
%{
int i,j;
PROP_Z0;
if (x>xmin && x<xmax && y>ymin && y<ymax)
{
i = floor((x - xmin)*nx/(xmax - xmin));
j = floor((y - ymin)*ny/(ymax - ymin));
PSD_N[i][j]++;
PSD_p[i][j] += p;
PSD_p2[i][j] += p*p;
}
%}
FINALLY
%{
int i,j,n=0;
double flux=0.0,xx,yy,A;
char string[40];
for (i=0; i<nx; i++)
{
xx=xmin+(xmax-xmin)*(double)i/(double)nx;
for (j=0; j<ny; j++)
{
yy=ymin+(ymax-ymin)*(double)j/(double)ny;
if ((yy*yy+xx*xx)<=(0.01*0.01/PI))
{ n++;
flux +=PSD_p[i][j];
}
}
}
A=(double)n*(xmax-xmin)/(double)nx*(ymax-ymin)/(double)ny*100.0*100.0; /* cm2 */
flux/=A;
/* flux/=(double)created_neutrons; */
sprintf(string,"BiDimPSD %9.4lgn/s/cm2",flux);
DETECTOR_OUT_2D(
string,
"X position [cm]",
"Y position [cm]",
xmin*100.0, xmax*100.0, ymin*100.0, ymax*100.0,
nx, ny,
&PSD_N[0][0],&PSD_p[0][0],&PSD_p2[0][0],
filename);
%}
MCDISPLAY
%{
magnify("xy");
multiline(5, (double)xmin, (double)ymin, 0.0,
(double)xmax, (double)ymin, 0.0,
(double)xmax, (double)ymax, 0.0,
(double)xmin, (double)ymax, 0.0,
(double)xmin, (double)ymin, 0.0);
%}
END
More information about the mcstas-users
mailing list