Shoot horizontal laser in X direction (vertical geometry too involved):
19 tlaser--(){
21 local msg="=== $FUNCNAME :"
23 local det=$(tlaser-det)
24 local tag=$(tlaser-tag)
25
26 local torch_config=(
27 type=point
28 frame=3153
29 source=0,0,0
30 target=1,0,0
31 photons=10000
32 material=GdDopedLS
33 wavelength=430
34 weight=1.0
35 time=0.1
36 zenithazimuth=0,1,0,1
37 radius=0
38 )
40 op.sh \
41 $* \
42 --animtimemax 15 \
43 --timemax 15 \
44 --eye 0,1,0 \
45 --torch --torchconfig "$(join _ ${torch_config[@]})" \
46 --torchdbg \
47 --save --tag $tag --cat $det
51 }
tlaser- ; tlaser-- --okg4 --compute
Proceed from this in SURFACE_ABSORB
In general the progressive mask totals show good step-by-step agreement, discrepancies coming in only at last step (AB or SA).
[2016-10-28 11:16:32,771] p43831 {/Users/blyth/opticks/ana/tlaser.py:48} INFO - a : laser/torch/ 1 : 20161028-1116 /tmp/blyth/opticks/evt/laser/torch/1/fdom.npy
[2016-10-28 11:16:32,772] p43831 {/Users/blyth/opticks/ana/tlaser.py:49} INFO - b : laser/torch/ -1 : 20161028-1116 /tmp/blyth/opticks/evt/laser/torch/-1/fdom.npy
seqhis_ana 1:laser -1:laser c2 ab ba
8ccccd 813163 813761 0.22 1.00 +- 0.00 1.00 +- 0.00 [6 ] TO BT BT BT BT SA
4d 45622 45617 0.00 1.00 +- 0.00 1.00 +- 0.00 [2 ] TO AB
cccc9ccccd 27443 27012 3.41 1.02 +- 0.01 0.98 +- 0.01 [10] TO BT BT BT BT DR BT BT BT BT
8cccc6d 15516 18592 277.41 0.83 +- 0.01 1.20 +- 0.01 [7 ] TO SC BT BT BT BT SA ## ~20% final SA
4ccd 10975 11210 2.49 0.98 +- 0.01 1.02 +- 0.01 [4 ] TO BT BT AB
4ccccd 9002 8820 1.86 1.02 +- 0.01 0.98 +- 0.01 [6 ] TO BT BT BT BT AB
8cccc5d 8433 8284 1.33 1.02 +- 0.01 0.98 +- 0.01 [7 ] TO RE BT BT BT BT SA
8cc6ccd 3370 3943 44.90 0.85 +- 0.01 1.17 +- 0.02 [7 ] TO BT BT SC BT BT SA ## ~20% final SA
cacccccc6d 3345 2435 143.27 1.37 +- 0.02 0.73 +- 0.01 [10] TO SC BT BT BT BT BT BT SR BT ## trunc
cccccc6ccd 2930 2396 53.54 1.22 +- 0.02 0.82 +- 0.02 [10] TO BT BT SC BT BT BT BT BT BT ## trunc
86ccccd 2554 2707 4.45 0.94 +- 0.02 1.06 +- 0.02 [7 ] TO BT BT BT BT SC SA ## ~20% final SA
45d 2436 2490 0.59 0.98 +- 0.02 1.02 +- 0.02 [3 ] TO RE AB
4ccccc6d 2431 78 2206.70 31.17 +- 0.63 0.03 +- 0.00 [8 ] TO SC BT BT BT BT BT AB ## drastic AB discrep
tlaser-v shows to be associated with specific geometry in viscinity of bottom reflector
tlaser-vg4 cannot show the 78 as does not make it into the top chart
8cccc55d 2180 2119 0.87 1.03 +- 0.02 0.97 +- 0.02 [8 ] TO RE RE BT BT BT BT SA
89ccccd 2011 2152 4.78 0.93 +- 0.02 1.07 +- 0.02 [7 ] TO BT BT BT BT DR SA ## final SA
cccc6ccccd 2068 1750 26.49 1.18 +- 0.03 0.85 +- 0.02 [10] TO BT BT BT BT SC BT BT BT BT ## trunc
4cccd 2065 1990 1.39 1.04 +- 0.02 0.96 +- 0.02 [5 ] TO BT BT BT AB
8ccccc6d 991 1985 332.00 0.50 +- 0.02 2.00 +- 0.04 [8 ] TO SC BT BT BT BT BT SA ## final SA (OK is half of G4)
8cc5ccd 1898 1964 1.13 0.97 +- 0.02 1.03 +- 0.02 [7 ] TO BT BT RE BT BT SA
ccbccccc6d 1621 1309 33.22 1.24 +- 0.03 0.81 +- 0.02 [10] TO SC BT BT BT BT BT BR BT BT ## trunc
1000000 1000000 37.28
Progressive mask development of the 20% discrepant 8cccc6d shows problem to be
all in final SURFACE_ABSORB SA step, with G4 absorbing 20% more than OK.
Note that top line SA is in agreement, but 2nd step SC means are going in a
random direction, indicating an issue with the "average" absorbing surface
that is not present with the direct surface pointed at by the laser.
6d 36156 35863 1.19 1.01 +- 0.01 0.99 +- 0.01 [2 ] TO SC
c6d 32422 32101 1.60 1.01 +- 0.01 0.99 +- 0.01 [3 ] TO SC BT
cc6d 32333 32014 1.58 1.01 +- 0.01 0.99 +- 0.01 [4 ] TO SC BT BT
ccc6d 31049 30857 0.60 1.01 +- 0.01 0.99 +- 0.01 [5 ] TO SC BT BT BT
cccc6d 30884 30721 0.43 1.01 +- 0.01 0.99 +- 0.01 [6 ] TO SC BT BT BT BT
8cccc6d 15516 18592 277.41 0.83 +- 0.01 1.20 +- 0.01 [7 ] TO SC BT BT BT BT SA
Same again issue with final SA.
cd 892640 893243 0.20 1.00 +- 0.00 1.00 +- 0.00 [2 ] TO BT
ccd 891267 891910 0.23 1.00 +- 0.00 1.00 +- 0.00 [3 ] TO BT BT
6ccd 9025 9035 0.01 1.00 +- 0.01 1.00 +- 0.01 [4 ] TO BT BT SC
c6ccd 8675 8640 0.07 1.00 +- 0.01 1.00 +- 0.01 [5 ] TO BT BT SC BT
cc6ccd 8446 8392 0.17 1.01 +- 0.01 0.99 +- 0.01 [6 ] TO BT BT SC BT BT
8cc6ccd 3370 3943 44.90 0.85 +- 0.01 1.17 +- 0.02 [7 ] TO BT BT SC BT BT SA
## before fix
seqhis_ana_4 1:laser -1:laser c2 ab ba
cccd 86472 86511 0.01 1.00 +- 0.00 1.00 +- 0.00 [4 ] TO BT BT BT
4d 4577 4623 0.23 0.99 +- 0.01 1.01 +- 0.01 [2 ] TO AB
cc6d 3246 3187 0.54 1.02 +- 0.02 0.98 +- 0.02 [4 ] TO SC BT BT
cc5d 1589 1467 4.87 1.08 +- 0.03 0.92 +- 0.02 [4 ] TO RE BT BT
4ccd 1070 1101 0.44 0.97 +- 0.03 1.03 +- 0.03 [4 ] TO BT BT AB
6ccd 890 899 0.05 0.99 +- 0.03 1.01 +- 0.03 [4 ] TO BT BT SC
5ccd 606 616 0.08 0.98 +- 0.04 1.02 +- 0.04 [4 ] TO BT BT RE
c55d 416 418 0.00 1.00 +- 0.05 1.00 +- 0.05 [4 ] TO RE RE BT
45d 240 239 0.00 1.00 +- 0.06 1.00 +- 0.06 [3 ] TO RE AB
46d 195 160 3.45 1.22 +- 0.09 0.82 +- 0.06 [3 ] TO SC AB
555d 175 193 0.88 0.91 +- 0.07 1.10 +- 0.08 [4 ] TO RE RE RE
4cd 126 101 2.75 1.25 +- 0.11 0.80 +- 0.08 [3 ] TO BT AB
c66d 124 107 1.25 1.16 +- 0.10 0.86 +- 0.08 [4 ] TO SC SC BT
c56d 57 74 2.21 0.77 +- 0.10 1.30 +- 0.15 [4 ] TO SC RE BT
455d 57 60 0.08 0.95 +- 0.13 1.05 +- 0.14 [4 ] TO RE RE AB
c65d 55 53 0.04 1.04 +- 0.14 0.96 +- 0.13 [4 ] TO RE SC BT
c6d ## 0 40 40.00 0.00 +- 0.00 0.00 +- 0.00 [3 ] TO SC BT <<< fixed by SAM to BT
c5d ## 0 37 37.00 0.00 +- 0.00 0.00 +- 0.00 [3 ] TO RE BT <<< fixed by SAM to BT
556d 24 21 0.20 1.14 +- 0.23 0.88 +- 0.19 [4 ] TO SC RE RE
655d 21 15 1.00 1.40 +- 0.31 0.71 +- 0.18 [4 ] TO RE RE SC
## after fix, cleans up the nasty zeros
seqhis_ana_4 1:laser -1:laser c2 ab ba
cccd 86472 86512 0.01 1.00 +- 0.00 1.00 +- 0.00 [4 ] TO BT BT BT
4d 4577 4623 0.23 0.99 +- 0.01 1.01 +- 0.01 [2 ] TO AB
cc6d 3246 3228 0.05 1.01 +- 0.02 0.99 +- 0.02 [4 ] TO SC BT BT
cc5d 1589 1504 2.34 1.06 +- 0.03 0.95 +- 0.02 [4 ] TO RE BT BT
4ccd 1070 1101 0.44 0.97 +- 0.03 1.03 +- 0.03 [4 ] TO BT BT AB
6ccd 890 899 0.05 0.99 +- 0.03 1.01 +- 0.03 [4 ] TO BT BT SC
5ccd 606 616 0.08 0.98 +- 0.04 1.02 +- 0.04 [4 ] TO BT BT RE
c55d 416 418 0.00 1.00 +- 0.05 1.00 +- 0.05 [4 ] TO RE RE BT
45d 240 239 0.00 1.00 +- 0.06 1.00 +- 0.06 [3 ] TO RE AB
46d 195 160 3.45 1.22 +- 0.09 0.82 +- 0.06 [3 ] TO SC AB
555d 175 193 0.88 0.91 +- 0.07 1.10 +- 0.08 [4 ] TO RE RE RE
4cd 126 101 2.75 1.25 +- 0.11 0.80 +- 0.08 [3 ] TO BT AB
c66d 124 107 1.25 1.16 +- 0.10 0.86 +- 0.08 [4 ] TO SC SC BT
c56d 57 74 2.21 0.77 +- 0.10 1.30 +- 0.15 [4 ] TO SC RE BT
455d 57 60 0.08 0.95 +- 0.13 1.05 +- 0.14 [4 ] TO RE RE AB
c65d 55 53 0.04 1.04 +- 0.14 0.96 +- 0.13 [4 ] TO RE SC BT
556d 24 21 0.20 1.14 +- 0.23 0.88 +- 0.19 [4 ] TO SC RE RE
655d 21 15 1.00 1.40 +- 0.31 0.71 +- 0.18 [4 ] TO RE RE SC
c6cd 16 14 0.00 1.14 +- 0.29 0.88 +- 0.23 [4 ] TO BT SC BT
466d 1 10 0.00 0.10 +- 0.10 10.00 +- 3.16 [4 ] TO SC SC AB
100000 100000 0.84
seqhis_ana 1:laser -1:laser c2 ab ba
8ccccd 81381 81398 0.00 1.00 +- 0.00 1.00 +- 0.00 [6 ] TO BT BT BT BT SA
4d 4577 4623 0.23 0.99 +- 0.01 1.01 +- 0.01 [2 ] TO AB
cccc9ccccd 2627 2685 0.63 0.98 +- 0.02 1.02 +- 0.02 [10] TO BT BT BT BT DR BT BT BT BT
8cccc6d 1570 1845 22.14 0.85 +- 0.02 1.18 +- 0.03 [7 ] TO SC BT BT BT BT SA
4ccd 1070 1101 0.44 0.97 +- 0.03 1.03 +- 0.03 [4 ] TO BT BT AB
4ccccd 859 846 0.10 1.02 +- 0.03 0.98 +- 0.03 [6 ] TO BT BT BT BT AB
8cccc5d 815 808 0.03 1.01 +- 0.04 0.99 +- 0.03 [7 ] TO RE BT BT BT BT SA
8cc6ccd 356 400 2.56 0.89 +- 0.05 1.12 +- 0.06 [7 ] TO BT BT SC BT BT SA
cacccccc6d 312 236 10.54 1.32 +- 0.07 0.76 +- 0.05 [10] TO SC BT BT BT BT BT BT SR BT
cccccc6ccd 291 226 8.17 1.29 +- 0.08 0.78 +- 0.05 [10] TO BT BT SC BT BT BT BT BT BT
86ccccd 268 267 0.00 1.00 +- 0.06 1.00 +- 0.06 [7 ] TO BT BT BT BT SC SA
4ccccc6d 257 4 245.25 64.25 +- 4.01 0.02 +- 0.01 [8 ] TO SC BT BT BT BT BT AB
45d 240 239 0.00 1.00 +- 0.06 1.00 +- 0.06 [3 ] TO RE AB
8cccc55d 222 224 0.01 0.99 +- 0.07 1.01 +- 0.07 [8 ] TO RE RE BT BT BT BT SA
4cccd 220 198 1.16 1.11 +- 0.07 0.90 +- 0.06 [5 ] TO BT BT BT AB
cccc6ccccd 216 174 4.52 1.24 +- 0.08 0.81 +- 0.06 [10] TO BT BT BT BT SC BT BT BT BT
8ccccc6d 100 211 39.62 0.47 +- 0.05 2.11 +- 0.15 [8 ] TO SC BT BT BT BT BT SA
89ccccd 203 201 0.01 1.01 +- 0.07 0.99 +- 0.07 [7 ] TO BT BT BT BT DR SA
8cc5ccd 187 200 0.44 0.94 +- 0.07 1.07 +- 0.08 [7 ] TO BT BT RE BT BT SA
46d 195 160 3.45 1.22 +- 0.09 0.82 +- 0.06 [3 ] TO SC AB
100000 100000 13.72
seqhis_ana_1 1:laser -1:laser c2 ab ba
d 100000 100000 0.00 1.00 +- 0.00 1.00 +- 0.00 [1 ] TO
100000 100000 0.00
seqhis_ana_2 1:laser -1:laser c2 ab ba
cd 89182 89249 0.03 1.00 +- 0.00 1.00 +- 0.00 [2 ] TO BT
4d 4577 4623 0.23 0.99 +- 0.01 1.01 +- 0.01 [2 ] TO AB
6d 3674 3629 0.28 1.01 +- 0.02 0.99 +- 0.02 [2 ] TO SC
5d 2566 2498 0.91 1.03 +- 0.02 0.97 +- 0.02 [2 ] TO RE
bd 1 1 0.00 1.00 +- 1.00 1.00 +- 1.00 [2 ] TO BR
100000 100000 0.36
seqhis_ana_3 1:laser -1:laser c2 ab ba
ccd 89038 89128 0.05 1.00 +- 0.00 1.00 +- 0.00 [3 ] TO BT BT
4d 4577 4623 0.23 0.99 +- 0.01 1.01 +- 0.01 [2 ] TO AB
c6d 3253 3237 0.04 1.00 +- 0.02 1.00 +- 0.02 [3 ] TO SC BT
c5d 1593 1508 2.33 1.06 +- 0.03 0.95 +- 0.02 [3 ] TO RE BT
55d 669 686 0.21 0.98 +- 0.04 1.03 +- 0.04 [3 ] TO RE RE
45d 240 239 0.00 1.00 +- 0.06 1.00 +- 0.06 [3 ] TO RE AB
46d 195 160 3.45 1.22 +- 0.09 0.82 +- 0.06 [3 ] TO SC AB
66d 134 122 0.56 1.10 +- 0.09 0.91 +- 0.08 [3 ] TO SC SC
4cd 126 101 2.75 1.25 +- 0.11 0.80 +- 0.08 [3 ] TO BT AB
56d 92 108 1.28 0.85 +- 0.09 1.17 +- 0.11 [3 ] TO SC RE
65d 63 65 0.03 0.97 +- 0.12 1.03 +- 0.13 [3 ] TO RE SC
6cd 18 18 0.00 1.00 +- 0.24 1.00 +- 0.24 [3 ] TO BT SC
6d 0 1 0.00 0.00 +- 0.00 0.00 +- 0.00 [2 ] TO SC
b6d 0 1 0.00 0.00 +- 0.00 0.00 +- 0.00 [3 ] TO SC BR
cbd 1 1 0.00 1.00 +- 1.00 1.00 +- 1.00 [3 ] TO BR BT
cd 0 1 0.00 0.00 +- 0.00 0.00 +- 0.00 [2 ] TO BT
b5d 1 0 0.00 0.00 +- 0.00 0.00 +- 0.00 [3 ] TO RE BR
bcd 0 1 0.00 0.00 +- 0.00 0.00 +- 0.00 [3 ] TO BT BR
100000 100000 0.91
tlaser-tfx()
{
tlaser-t --fxabconfig 10000 --fxab --fxscconfig 10000 --fxsc --fxreconfig 0.5 --fxre $*
}
seqhis_ana__ 1:laser -1:laser c2 ab ba
8ccccd 61164 60977 0.29 1.00 +- 0.00 1.00 +- 0.00 [6 ] TO BT BT BT BT SA
4d 6614 6749 1.36 0.98 +- 0.01 1.02 +- 0.01 [2 ] TO AB
8cccc6d 4284 5178 84.47 0.83 +- 0.01 1.21 +- 0.02 [7 ] TO SC BT BT BT BT SA
8cccc5d 2430 2404 0.14 1.01 +- 0.02 0.99 +- 0.02 [7 ] TO RE BT BT BT BT SA
4ccccd 1712 1714 0.00 1.00 +- 0.02 1.00 +- 0.02 [6 ] TO BT BT BT BT AB
cccc9ccccd 1691 1646 0.61 1.03 +- 0.02 0.97 +- 0.02 [10] TO BT BT BT BT DR BT BT BT BT
4ccd 1416 1440 0.20 0.98 +- 0.03 1.02 +- 0.03 [4 ] TO BT BT AB
8cc6ccd 964 1103 9.35 0.87 +- 0.03 1.14 +- 0.03 [7 ] TO BT BT SC BT BT SA
cacccccc6d 951 597 80.95 1.59 +- 0.05 0.63 +- 0.03 [10] TO SC BT BT BT BT BT BT SR BT <<<<
46d 893 866 0.41 1.03 +- 0.03 0.97 +- 0.03 [3 ] TO SC AB
86ccccd 717 775 2.25 0.93 +- 0.03 1.08 +- 0.04 [7 ] TO BT BT BT BT SC SA
cccccc6ccd 667 511 20.66 1.31 +- 0.05 0.77 +- 0.03 [10] TO BT BT SC BT BT BT BT BT BT <<<<
8cccc66d 575 651 4.71 0.88 +- 0.04 1.13 +- 0.04 [8 ] TO SC SC BT BT BT BT SA
cccc6ccccd 576 461 12.75 1.25 +- 0.05 0.80 +- 0.04 [10] TO BT BT BT BT SC BT BT BT BT <<<<
8cc5ccd 540 542 0.00 1.00 +- 0.04 1.00 +- 0.04 [7 ] TO BT BT RE BT BT SA
cccccccc6d 522 53 382.54 9.85 +- 0.43 0.10 +- 0.01 [10] TO SC BT BT BT BT BT BT BT BT <<<< TRUNCATION BEHAVIOUR MISMATCH ???
8ccccc6d 281 505 63.84 0.56 +- 0.03 1.80 +- 0.08 [8 ] TO SC BT BT BT BT BT SA
45d 455 412 2.13 1.10 +- 0.05 0.91 +- 0.04 [3 ] TO RE AB
ccbccccc6d 429 349 8.23 1.23 +- 0.06 0.81 +- 0.04 [10] TO SC BT BT BT BT BT BR BT BT
cacccccc5d 393 347 2.86 1.13 +- 0.06 0.88 +- 0.05 [10] TO RE BT BT BT BT BT BT SR BT
100000 100000 20.83
seqhis_ana_1 1:laser -1:laser c2 ab ba
d 100000 100000 0.00 1.00 +- 0.00 1.00 +- 0.00 [1 ] TO
100000 100000 0.00
seqhis_ana_2 1:laser -1:laser c2 ab ba
cd 73222 73372 0.15 1.00 +- 0.00 1.00 +- 0.00 [2 ] TO BT
6d 13499 13327 1.10 1.01 +- 0.01 0.99 +- 0.01 [2 ] TO SC
4d 6614 6749 1.36 0.98 +- 0.01 1.02 +- 0.01 [2 ] TO AB
5d 6664 6552 0.95 1.02 +- 0.01 0.98 +- 0.01 [2 ] TO RE
bd 1 0 0.00 0.00 +- 0.00 0.00 +- 0.00 [2 ] TO BR
100000 100000 0.89
seqhis_ana_3 1:laser -1:laser c2 ab ba
ccd 73075 73211 0.13 1.00 +- 0.00 1.00 +- 0.00 [3 ] TO BT BT
c6d 9922 9829 0.44 1.01 +- 0.01 0.99 +- 0.01 [3 ] TO SC BT
4d 6614 6749 1.36 0.98 +- 0.01 1.02 +- 0.01 [2 ] TO AB
c5d 4943 4893 0.25 1.01 +- 0.01 0.99 +- 0.01 [3 ] TO RE BT
66d 1784 1744 0.45 1.02 +- 0.02 0.98 +- 0.02 [3 ] TO SC SC
56d 897 885 0.08 1.01 +- 0.03 0.99 +- 0.03 [3 ] TO SC RE
46d 893 866 0.41 1.03 +- 0.03 0.97 +- 0.03 [3 ] TO SC AB
65d 830 843 0.10 0.98 +- 0.03 1.02 +- 0.03 [3 ] TO RE SC
45d 455 412 2.13 1.10 +- 0.05 0.91 +- 0.04 [3 ] TO RE AB
55d 436 404 1.22 1.08 +- 0.05 0.93 +- 0.05 [3 ] TO RE RE
4cd 81 86 0.15 0.94 +- 0.10 1.06 +- 0.11 [3 ] TO BT AB
6cd 66 70 0.12 0.94 +- 0.12 1.06 +- 0.13 [3 ] TO BT SC
b6d 3 3 0.00 1.00 +- 0.58 1.00 +- 0.58 [3 ] TO SC BR
cd 0 3 0.00 0.00 +- 0.00 0.00 +- 0.00 [2 ] TO BT
bcd 0 2 0.00 0.00 +- 0.00 0.00 +- 0.00 [3 ] TO BT BR
cbd 1 0 0.00 0.00 +- 0.00 0.00 +- 0.00 [3 ] TO BR BT
100000 100000 0.57
seqhis_ana_4 1:laser -1:laser c2 ab ba
cccd 67354 67347 0.00 1.00 +- 0.00 1.00 +- 0.00 [4 ] TO BT BT BT
cc6d 9893 9645 3.15 1.03 +- 0.01 0.97 +- 0.01 [4 ] TO SC BT BT
4d 6614 6749 1.36 0.98 +- 0.01 1.02 +- 0.01 [2 ] TO AB
cc5d 4930 4747 3.46 1.04 +- 0.01 0.96 +- 0.01 [4 ] TO RE BT BT
6ccd 2877 2966 1.36 0.97 +- 0.02 1.03 +- 0.02 [4 ] TO BT BT SC
5ccd 1428 1456 0.27 0.98 +- 0.03 1.02 +- 0.03 [4 ] TO BT BT RE
4ccd 1416 1440 0.20 0.98 +- 0.03 1.02 +- 0.03 [4 ] TO BT BT AB
c66d 1349 1346 0.00 1.00 +- 0.03 1.00 +- 0.03 [4 ] TO SC SC BT
46d 893 866 0.41 1.03 +- 0.03 0.97 +- 0.03 [3 ] TO SC AB
c56d 692 666 0.50 1.04 +- 0.04 0.96 +- 0.04 [4 ] TO SC RE BT
c65d 635 629 0.03 1.01 +- 0.04 0.99 +- 0.04 [4 ] TO RE SC BT
45d 455 412 2.13 1.10 +- 0.05 0.91 +- 0.04 [3 ] TO RE AB
c55d 344 297 3.45 1.16 +- 0.06 0.86 +- 0.05 [4 ] TO RE RE BT
666d 204 210 0.09 0.97 +- 0.07 1.03 +- 0.07 [4 ] TO SC SC SC
c6d 0 138 138.00 0.00 +- 0.00 0.00 +- 0.00 [3 ] TO SC BT ## whats this, different from above ???
c5d 0 130 130.00 0.00 +- 0.00 0.00 +- 0.00 [3 ] TO RE BT ## again diff to above ??? maybe bad abbr zeros effect
^^^^^^^^^^^^^^^^ maybe bad abbr zero : missing flags effect ?? ^^^^^^^^^^^^^^^^^^^^^^^^
566d 116 88 3.84 1.32 +- 0.12 0.76 +- 0.08 [4 ] TO SC SC RE
656d 116 104 0.65 1.12 +- 0.10 0.90 +- 0.09 [4 ] TO SC RE SC
466d 115 99 1.20 1.16 +- 0.11 0.86 +- 0.09 [4 ] TO SC SC AB
665d 90 114 2.82 0.79 +- 0.08 1.27 +- 0.12 [4 ] TO RE SC SC
100000 100000 10.22
seqhis_ana_5 1:laser -1:laser c2 ab ba
ccccd 67087 67078 0.00 1.00 +- 0.00 1.00 +- 0.00 [5 ] TO BT BT BT BT
ccc6d 8882 8656 2.91 1.03 +- 0.01 0.97 +- 0.01 [5 ] TO SC BT BT BT
4d 6614 6749 1.36 0.98 +- 0.01 1.02 +- 0.01 [2 ] TO AB
ccc5d 4419 4256 3.06 1.04 +- 0.02 0.96 +- 0.01 [5 ] TO RE BT BT BT
c6ccd 2617 2614 0.00 1.00 +- 0.02 1.00 +- 0.02 [5 ] TO BT BT SC BT
4ccd 1416 1440 0.20 0.98 +- 0.03 1.02 +- 0.03 [4 ] TO BT BT AB
cc66d 1344 1313 0.36 1.02 +- 0.03 0.98 +- 0.03 [5 ] TO SC SC BT BT
c5ccd 1244 1288 0.76 0.97 +- 0.03 1.04 +- 0.03 [5 ] TO BT BT RE BT
46d 893 866 0.41 1.03 +- 0.03 0.97 +- 0.03 [3 ] TO SC AB
cc56d 689 651 1.08 1.06 +- 0.04 0.94 +- 0.04 [5 ] TO SC RE BT BT
cc65d 634 609 0.50 1.04 +- 0.04 0.96 +- 0.04 [5 ] TO RE SC BT BT
6cc6d 509 486 0.53 1.05 +- 0.05 0.95 +- 0.04 [5 ] TO SC BT BT SC
45d 455 412 2.13 1.10 +- 0.05 0.91 +- 0.04 [3 ] TO RE AB
cc55d 342 289 4.45 1.18 +- 0.06 0.85 +- 0.05 [5 ] TO RE RE BT BT
6cc5d 269 268 0.00 1.00 +- 0.06 1.00 +- 0.06 [5 ] TO RE BT BT SC
4cc6d 259 249 0.20 1.04 +- 0.06 0.96 +- 0.06 [5 ] TO SC BT BT AB
5cc6d 225 254 1.76 0.89 +- 0.06 1.13 +- 0.07 [5 ] TO SC BT BT RE
66ccd 131 187 9.86 0.70 +- 0.06 1.43 +- 0.10 [5 ] TO BT BT SC SC
c666d 161 160 0.00 1.01 +- 0.08 0.99 +- 0.08 [5 ] TO SC SC SC BT
4cccd 142 118 2.22 1.20 +- 0.10 0.83 +- 0.08 [5 ] TO BT BT BT AB
100000 100000 7.81
Possible causes of discrep
seqhis_ana 1:laser -1:laser c2 ab ba
8ccccd 76521 81427 152.38 0.94 +- 0.00 1.06 +- 0.00 [6 ] TO BT BT BT BT SA
4d 5573 4758 64.29 1.17 +- 0.02 0.85 +- 0.01 [2 ] TO AB
cccc9ccccd 2428 2700 14.43 0.90 +- 0.02 1.11 +- 0.02 [10] TO BT BT BT BT DR BT BT BT BT
8cccc6d 1590 1863 21.58 0.85 +- 0.02 1.17 +- 0.03 [7 ] TO SC BT BT BT BT SA
4ccd 1194 1133 1.60 1.05 +- 0.03 0.95 +- 0.03 [4 ] TO BT BT AB
8cccc5d 1074 750 57.55 1.43 +- 0.04 0.70 +- 0.03 [7 ] TO RE BT BT BT BT SA
4ccccd 822 828 0.02 0.99 +- 0.03 1.01 +- 0.04 [6 ] TO BT BT BT BT AB
45d 754 216 298.40 3.49 +- 0.13 0.29 +- 0.02 [3 ] TO RE AB
8cccc55d 561 211 158.68 2.66 +- 0.11 0.38 +- 0.03 [8 ] TO RE RE BT BT BT BT SA
8cc6ccd 366 382 0.34 0.96 +- 0.05 1.04 +- 0.05 [7 ] TO BT BT SC BT BT SA
455d 345 47 226.54 7.34 +- 0.40 0.14 +- 0.02 [4 ] TO RE RE AB
cacccccc6d 325 228 17.01 1.43 +- 0.08 0.70 +- 0.05 [10] TO SC BT BT BT BT BT BT SR BT
86ccccd 291 268 0.95 1.09 +- 0.06 0.92 +- 0.06 [7 ] TO BT BT BT BT SC SA
cccccc6ccd 291 239 5.10 1.22 +- 0.07 0.82 +- 0.05 [10] TO BT BT SC BT BT BT BT BT BT
4ccccc6d 263 5 248.37 52.60 +- 3.24 0.02 +- 0.01 [8 ] TO SC BT BT BT BT BT AB
46d 244 165 15.26 1.48 +- 0.09 0.68 +- 0.05 [3 ] TO SC AB
8cccc555d 243 56 116.95 4.34 +- 0.28 0.23 +- 0.03 [9 ] TO RE RE RE BT BT BT BT SA
8cc5ccd 236 191 4.74 1.24 +- 0.08 0.81 +- 0.06 [7 ] TO BT BT RE BT BT SA
cccc6ccccd 227 164 10.15 1.38 +- 0.09 0.72 +- 0.06 [10] TO BT BT BT BT SC BT BT BT BT
8ccccc6d 116 223 33.77 0.52 +- 0.05 1.92 +- 0.13 [8 ] TO SC BT BT BT BT BT SA
100000 100000 32.58
seqhis_ana_1 1:laser -1:laser c2 ab ba
d 100000 100000 0.00 1.00 +- 0.00 1.00 +- 0.00 [1 ] TO
100000 100000 0.00
seqhis_ana_2 1:laser -1:laser c2 ab ba
cd 84925 89281 108.92 0.95 +- 0.00 1.05 +- 0.00 [2 ] TO BT
4d 5573 4758 64.29 1.17 +- 0.02 0.85 +- 0.01 [2 ] TO AB
5d 5457 2348 1238.42 2.32 +- 0.03 0.43 +- 0.01 [2 ] TO RE
6d 4044 3612 24.38 1.12 +- 0.02 0.89 +- 0.01 [2 ] TO SC
bd 1 1 0.00 1.00 +- 1.00 1.00 +- 1.00 [2 ] TO BR
100000 100000 359.00
seqhis_ana_3 1:laser -1:laser c2 ab ba
ccd 84790 89153 109.44 0.95 +- 0.00 1.05 +- 0.00 [3 ] TO BT BT
4d 5573 4758 64.29 1.17 +- 0.02 0.85 +- 0.01 [2 ] TO AB
c6d 3406 3217 5.39 1.06 +- 0.02 0.94 +- 0.02 [3 ] TO SC BT
55d 2595 631 1195.69 4.11 +- 0.08 0.24 +- 0.01 [3 ] TO RE RE
c5d 2034 1428 106.08 1.42 +- 0.03 0.70 +- 0.02 [3 ] TO RE BT
45d 754 216 298.40 3.49 +- 0.13 0.29 +- 0.02 [3 ] TO RE AB
46d 244 165 15.26 1.48 +- 0.09 0.68 +- 0.05 [3 ] TO SC AB
56d 230 100 51.21 2.30 +- 0.15 0.43 +- 0.04 [3 ] TO SC RE
66d 164 128 4.44 1.28 +- 0.10 0.78 +- 0.07 [3 ] TO SC SC
4cd 116 100 1.19 1.16 +- 0.11 0.86 +- 0.09 [3 ] TO BT AB
65d 74 73 0.01 1.01 +- 0.12 0.99 +- 0.12 [3 ] TO RE SC
6cd 19 26 1.09 0.73 +- 0.17 1.37 +- 0.27 [3 ] TO BT SC
bcd 0 1 0.00 0.00 +- 0.00 0.00 +- 0.00 [3 ] TO BT BR
b6d 0 1 0.00 0.00 +- 0.00 0.00 +- 0.00 [3 ] TO SC BR
cbd 1 1 0.00 1.00 +- 1.00 1.00 +- 1.00 [3 ] TO BR BT
cd 0 1 0.00 0.00 +- 0.00 0.00 +- 0.00 [2 ] TO BT
6d 0 1 0.00 0.00 +- 0.00 0.00 +- 0.00 [2 ] TO SC
100000 100000 154.37
seqhis_ana 1:laser -1:laser c2 ab ba
8ccccd 76521 81336 146.87 0.94 1.06 [6 ] TO BT BT BT BT SA
4d 5573 4699 74.36 1.19 0.84 [2 ] TO AB
cccc9ccccd 2428 2661 10.67 0.91 1.10 [10] TO BT BT BT BT DR BT BT BT BT
8cccc6d 1980 1899 1.69 1.04 0.96 [7 ] TO SC BT BT BT BT SA
4ccd 1194 1161 0.46 1.03 0.97 [4 ] TO BT BT AB
8cccc5d 1074 753 56.40 1.43 0.70 [7 ] TO RE BT BT BT BT SA
4ccccd 822 858 0.77 0.96 1.04 [6 ] TO BT BT BT BT AB
45d 754 211 305.54 3.57 0.28 [3 ] TO RE AB
8cccc55d 561 230 138.51 2.44 0.41 [8 ] TO RE RE BT BT BT BT SA
8cc6ccd 413 403 0.12 1.02 0.98 [7 ] TO BT BT SC BT BT SA
455d 345 67 187.58 5.15 0.19 [4 ] TO RE RE AB
86ccccd 299 263 2.31 1.14 0.88 [7 ] TO BT BT BT BT SC SA
cccccc6ccd 262 198 8.90 1.32 0.76 [10] TO BT BT SC BT BT BT BT BT BT
8cccc555d 243 66 101.39 3.68 0.27 [9 ] TO RE RE RE BT BT BT BT SA
8cc5ccd 236 190 4.97 1.24 0.81 [7 ] TO BT BT RE BT BT SA
cccc6ccccd 229 164 10.75 1.40 0.72 [10] TO BT BT BT BT SC BT BT BT BT
89ccccd 191 218 1.78 0.88 1.14 [7 ] TO BT BT BT BT DR SA
46d 217 141 16.13 1.54 0.65 [3 ] TO SC AB
4cccd 209 207 0.01 1.01 0.99 [5 ] TO BT BT BT AB
cacccccc6d 205 208 0.02 0.99 1.01 [10] TO SC BT BT BT BT BT BT SR BT
100000 100000 29.77
seqhis_ana_1 1:laser -1:laser c2 ab ba
d 100000 100000 0.00 1.00 1.00 [1 ] TO
100000 100000 0.00
seqhis_ana_2 1:laser -1:laser c2 ab ba
cd 84925 89211 105.49 0.95 1.05 [2 ] TO BT <<< G4 5% more get to boundary without AB RE or SC happening
4d 5573 4699 74.36 1.19 0.84 [2 ] TO AB <<< Opticks 20% more AB
5d 5457 2411 1179.22 2.26 0.44 [2 ] TO RE <<< Opticks 2.2x RE
6d 4044 3678 17.35 1.10 0.91 [2 ] TO SC <<< Opticks 10% more SC
bd 1 1 0.00 1.00 1.00 [2 ] TO BR
100000 100000 344.11
Given tpmt excellent agreement (PMTInBox of mineral oil) suspect issue with scintillator
try to confirm by tpmt with scintillator...
seqhis_ana 1:laser -1:laser c2 ab ba
8ccccd 76521 81336 146.87 0.94 +- 0.00 1.06 +- 0.00 [6 ] TO BT BT BT BT SA
4d 5573 4699 74.36 1.19 +- 0.02 0.84 +- 0.01 [2 ] TO AB
cccc9ccccd 2428 2661 10.67 0.91 +- 0.02 1.10 +- 0.02 [10] TO BT BT BT BT DR BT BT BT BT
8cccc6d 1980 1899 1.69 1.04 +- 0.02 0.96 +- 0.02 [7 ] TO SC BT BT BT BT SA
4ccd 1194 1161 0.46 1.03 +- 0.03 0.97 +- 0.03 [4 ] TO BT BT AB
8cccc5d 1074 753 56.40 1.43 +- 0.04 0.70 +- 0.03 [7 ] TO RE BT BT BT BT SA
4ccccd 822 858 0.77 0.96 +- 0.03 1.04 +- 0.04 [6 ] TO BT BT BT BT AB
45d 754 211 305.54 3.57 +- 0.13 0.28 +- 0.02 [3 ] TO RE AB
8cccc55d 561 230 138.51 2.44 +- 0.10 0.41 +- 0.03 [8 ] TO RE RE BT BT BT BT SA
8cc6ccd 413 403 0.12 1.02 +- 0.05 0.98 +- 0.05 [7 ] TO BT BT SC BT BT SA
455d 345 67 187.58 5.15 +- 0.28 0.19 +- 0.02 [4 ] TO RE RE AB
86ccccd 299 263 2.31 1.14 +- 0.07 0.88 +- 0.05 [7 ] TO BT BT BT BT SC SA
cccccc6ccd 262 198 8.90 1.32 +- 0.08 0.76 +- 0.05 [10] TO BT BT SC BT BT BT BT BT BT
8cccc555d 243 66 101.39 3.68 +- 0.24 0.27 +- 0.03 [9 ] TO RE RE RE BT BT BT BT SA
8cc5ccd 236 190 4.97 1.24 +- 0.08 0.81 +- 0.06 [7 ] TO BT BT RE BT BT SA
cccc6ccccd 229 164 10.75 1.40 +- 0.09 0.72 +- 0.06 [10] TO BT BT BT BT SC BT BT BT BT
89ccccd 191 218 1.78 0.88 +- 0.06 1.14 +- 0.08 [7 ] TO BT BT BT BT DR SA
46d 217 141 16.13 1.54 +- 0.10 0.65 +- 0.05 [3 ] TO SC AB
4cccd 209 207 0.01 1.01 +- 0.07 0.99 +- 0.07 [5 ] TO BT BT BT AB
cacccccc6d 205 208 0.02 0.99 +- 0.07 1.01 +- 0.07 [10] TO SC BT BT BT BT BT BT SR BT
100000 100000 29.77
seqhis_ana_1 1:laser -1:laser c2 ab ba
d 100000 100000 0.00 1.00 +- 0.00 1.00 +- 0.00 [1 ] TO
100000 100000 0.00
seqhis_ana_2 1:laser -1:laser c2 ab ba
cd 84925 89211 105.49 0.95 +- 0.00 1.05 +- 0.00 [2 ] TO BT
4d 5573 4699 74.36 1.19 +- 0.02 0.84 +- 0.01 [2 ] TO AB
5d 5457 2411 1179.22 2.26 +- 0.03 0.44 +- 0.01 [2 ] TO RE
6d 4044 3678 17.35 1.10 +- 0.02 0.91 +- 0.01 [2 ] TO SC
bd 1 1 0.00 1.00 +- 1.00 1.00 +- 1.00 [2 ] TO BR
100000 100000 344.11
seqhis_ana_3 1:laser -1:laser c2 ab ba
ccd 84790 89093 106.48 0.95 +- 0.00 1.05 +- 0.00 [3 ] TO BT BT
4d 5573 4699 74.36 1.19 +- 0.02 0.84 +- 0.01 [2 ] TO AB
c6d 3440 3320 2.13 1.04 +- 0.02 0.97 +- 0.02 [3 ] TO SC BT
55d 2595 704 1083.93 3.69 +- 0.07 0.27 +- 0.01 [3 ] TO RE RE
c5d 2034 1431 104.94 1.42 +- 0.03 0.70 +- 0.02 [3 ] TO RE BT
45d 754 211 305.54 3.57 +- 0.13 0.28 +- 0.02 [3 ] TO RE AB
46d 217 141 16.13 1.54 +- 0.10 0.65 +- 0.05 [3 ] TO SC AB
56d 211 93 45.80 2.27 +- 0.16 0.44 +- 0.05 [3 ] TO SC RE
66d 176 123 9.39 1.43 +- 0.11 0.70 +- 0.06 [3 ] TO SC SC
4cd 116 89 3.56 1.30 +- 0.12 0.77 +- 0.08 [3 ] TO BT AB
65d 74 63 0.88 1.17 +- 0.14 0.85 +- 0.11 [3 ] TO RE SC
6cd 19 28 1.72 0.68 +- 0.16 1.47 +- 0.28 [3 ] TO BT SC
b5d 0 2 0.00 0.00 +- 0.00 0.00 +- 0.00 [3 ] TO RE BR
bcd 0 1 0.00 0.00 +- 0.00 0.00 +- 0.00 [3 ] TO BT BR
b6d 0 1 0.00 0.00 +- 0.00 0.00 +- 0.00 [3 ] TO SC BR
cbd 1 1 0.00 1.00 +- 1.00 1.00 +- 1.00 [3 ] TO BR BT
100000 100000 146.24
seqhis_ana_4 1:laser -1:laser c2 ab ba
cccd 81407 86458 151.98 0.94 +- 0.00 1.06 +- 0.00 [4 ] TO BT BT BT
4d 5573 4699 74.36 1.19 +- 0.02 0.84 +- 0.01 [2 ] TO AB
cc6d 3433 3254 4.79 1.06 +- 0.02 0.95 +- 0.02 [4 ] TO SC BT BT
cc5d 2028 1393 117.87 1.46 +- 0.03 0.69 +- 0.02 [4 ] TO RE BT BT
555d 1241 185 782.00 6.71 +- 0.19 0.15 +- 0.01 [4 ] TO RE RE RE
5ccd 1239 590 230.29 2.10 +- 0.06 0.48 +- 0.02 [4 ] TO BT BT RE
4ccd 1194 1161 0.46 1.03 +- 0.03 0.97 +- 0.03 [4 ] TO BT BT AB
c55d 966 434 202.16 2.23 +- 0.07 0.45 +- 0.02 [4 ] TO RE RE BT
6ccd 950 882 2.52 1.08 +- 0.03 0.93 +- 0.03 [4 ] TO BT BT SC
45d 754 211 305.54 3.57 +- 0.13 0.28 +- 0.02 [3 ] TO RE AB
455d 345 67 187.58 5.15 +- 0.28 0.19 +- 0.02 [4 ] TO RE RE AB
46d 217 141 16.13 1.54 +- 0.10 0.65 +- 0.05 [3 ] TO SC AB
c66d 153 108 7.76 1.42 +- 0.11 0.71 +- 0.07 [4 ] TO SC SC BT
4cd 116 89 3.56 1.30 +- 0.12 0.77 +- 0.08 [3 ] TO BT AB
556d 112 16 72.00 7.00 +- 0.66 0.14 +- 0.04 [4 ] TO SC RE RE
c56d 71 66 0.18 1.08 +- 0.13 0.93 +- 0.11 [4 ] TO SC RE BT
c65d 59 51 0.58 1.16 +- 0.15 0.86 +- 0.12 [4 ] TO RE SC BT
c6d 0 59 59.00 0.00 +- 0.00 0.00 +- 0.00 [3 ] TO SC BT
655d 43 18 10.25 2.39 +- 0.36 0.42 +- 0.10 [4 ] TO RE RE SC
c5d 0 36 36.00 0.00 +- 0.00 0.00 +- 0.00 [3 ] TO RE BT
100000 100000 107.88
tlaser-;tlaser-t
tlaser.py
seqhis_ana 1:laser -1:laser c2
8ccccd 763501 813497 1585.04 [6 ] TO BT BT BT BT SA
cccc9ccccd 25263 26200 17.06 [10] TO BT BT BT BT DR BT BT BT BT
In [2]: 25263./(763501.+25263.)
Out[2]: 0.03202859156858072
In [3]: 26200./(813497.+26200.)
Out[3]: 0.0312017311006232
In [1]: 813497./763501. TODO: include the ratio in the output (expected reflectivity is ballpark 4%)
Out[1]: 1.0654825599442568
4d 55825 47634 648.49 [2 ] TO AB
8cccc6d 19707 18533 36.04 [7 ] TO SC BT BT BT BT SA
4ccd 12576 11563 42.51 [4 ] TO BT BT AB
8cccc5d 11183 7742 625.65 [7 ] TO RE BT BT BT BT SA
4ccccd 8554 8756 2.36 [6 ] TO BT BT BT BT AB
45d 7531 2208 2909.37 [3 ] TO RE AB
8cccc55d 5362 2116 1409.00 [8 ] TO RE RE BT BT BT BT SA
8cc6ccd 4109 4155 0.26 [7 ] TO BT BT SC BT BT SA
455d 3588 621 2091.49 [4 ] TO RE RE AB
86ccccd 2836 2743 1.55 [7 ] TO BT BT BT BT SC SA
cccccc6ccd 2674 1919 124.11 [10] TO BT BT SC BT BT BT BT BT BT
8cccc555d 2524 610 1168.92 [9 ] TO RE RE RE BT BT BT BT SA
8cc5ccd 2359 1866 57.53 [7 ] TO BT BT RE BT BT SA
89ccccd 1880 2221 28.35 [7 ] TO BT BT BT BT DR SA
cacccccc6d 2210 2127 1.59 [10] TO SC BT BT BT BT BT BT SR BT
46d 2118 1569 81.75 [3 ] TO SC AB
cccc6ccccd 2060 1752 24.89 [10] TO BT BT BT BT SC BT BT BT BT
4cccd 1940 1981 0.43 [5 ] TO BT BT BT AB
1000000 1000000 106.82
Dump top line, RSOilSurface as dielectric_metal when its MO/Ac ?:
tlaser-;tlaser-t --dbgseqhis 8ccccd
----CRecorder::compare---- record_id 5 --dindex
2016-10-25 20:11:36.056 INFO [3525262] [CRecorder::Dump@847] CRecorder::compare (rdr-dump)DONE record_id 5
2016-10-25 20:11:36.056 INFO [3525262] [CRecorder::Dump@850] seqhis 8ccccd TORCH BOUNDARY_TRANSMIT BOUNDARY_TRANSMIT BOUNDARY_TRANSMIT BOUNDARY_TRANSMIT SURFACE_ABSORB . . . . . . . . . .
2016-10-25 20:11:36.056 INFO [3525262] [CRecorder::Dump@854] seqmat 343231 GdDopedLS Acrylic LiquidScintillator Acrylic MineralOil Acrylic - - - - - - - - - -
2016-10-25 20:11:36.056 INFO [3525262] [CRec::dump@40] crec record_id 5 nstp 5 Ori[ -18079.453-799699.438-6605.000]
0[ 0](Stp ;opticalphoton stepNum 1513010768(tk ;opticalphoton tid 6 pid 0 nm 430 mm ori[ -18079.453-799699.438-6605.000] pos[ 1255.240-1878.345 0.000] )
pre d/Geometry/AD/lvIAV#pvGDS rials/GdDopedLS noProc Undefined pos[ 0.000 0.000 0.000] dir[ 0.556 -0.831 0.000] pol[ -1.000 0.023 0.000] ns 0.100 nm 430.000
post d/Geometry/AD/lvLSO#pvIAV terials/Acrylic Transportation GeomBoundary pos[ 861.221 -1288.733 0.000] dir[ 0.556 -0.831 0.000] pol[ -1.000 0.023 0.000] ns 8.059 nm 430.000
)
1[ 1](Stp ;opticalphoton stepNum 1513010768(tk ;opticalphoton tid 6 pid 0 nm 430 mm ori[ -18079.453-799699.438-6605.000] pos[ 1255.240-1878.345 0.000] )
pre d/Geometry/AD/lvLSO#pvIAV terials/Acrylic Transportation GeomBoundary pos[ 861.221 -1288.733 0.000] dir[ 0.556 -0.831 0.000] pol[ -1.000 0.023 0.000] ns 8.059 nm 430.000
post d/Geometry/AD/lvOAV#pvLSO uidScintillator Transportation GeomBoundary pos[ 866.777 -1297.048 0.000] dir[ 0.556 -0.831 0.000] pol[ -1.000 0.023 0.000] ns 8.110 nm 430.000
)
2[ 2](Stp ;opticalphoton stepNum 1513010768(tk ;opticalphoton tid 6 pid 0 nm 430 mm ori[ -18079.453-799699.438-6605.000] pos[ 1255.240-1878.345 0.000] )
pre d/Geometry/AD/lvOAV#pvLSO uidScintillator Transportation GeomBoundary pos[ 866.777 -1297.048 0.000] dir[ 0.556 -0.831 0.000] pol[ -1.000 0.023 0.000] ns 8.110 nm 430.000
post d/Geometry/AD/lvOIL#pvOAV terials/Acrylic Transportation GeomBoundary pos[ 1101.250 -1647.913 0.000] dir[ 0.556 -0.831 0.000] pol[ -1.000 0.023 0.000] ns 10.301 nm 430.000
)
3[ 3](Stp ;opticalphoton stepNum 1513010768(tk ;opticalphoton tid 6 pid 0 nm 430 mm ori[ -18079.453-799699.438-6605.000] pos[ 1255.240-1878.345 0.000] )
pre d/Geometry/AD/lvOIL#pvOAV terials/Acrylic Transportation GeomBoundary pos[ 1101.250 -1647.913 0.000] dir[ 0.556 -0.831 0.000] pol[ -1.000 0.023 0.000] ns 10.301 nm 430.000
post d/Geometry/AD/lvSST#pvOIL ials/MineralOil Transportation GeomBoundary pos[ 1111.251 -1662.879 0.000] dir[ 0.556 -0.831 0.000] pol[ -1.000 0.023 0.000] ns 10.393 nm 430.000
)
4[ 4](Stp ;opticalphoton stepNum 1513010768(tk ;opticalphoton tid 6 pid 0 nm 430 mm ori[ -18079.453-799699.438-6605.000] pos[ 1255.240-1878.345 0.000] )
pre d/Geometry/AD/lvSST#pvOIL ials/MineralOil Transportation GeomBoundary pos[ 1111.251 -1662.879 0.000] dir[ 0.556 -0.831 0.000] pol[ -1.000 0.023 0.000] ns 10.393 nm 430.000
post D/lvOIL#pvRadialShield:20 terials/Acrylic Transportation GeomBoundary pos[ 1255.240 -1878.345 0.000] dir[ 0.556 -0.831 0.000] pol[ -1.000 0.023 0.000] ns 11.738 nm 430.000
)
2016-10-25 20:11:36.057 INFO [3525262] [*DsG4OpBoundaryProcess::PostStepDoIt@442] OpticalSurface name RSOilSurface thePhotonMomentum (eV) 2.88335 theReflectivity 0.0409174 theEfficiency 0. dielectric_metal ground - m1 /dd/Materials/MineralOil m2 /dd/Materials/Acrylic
2016-10-25 20:11:36.057 INFO [3525262] [*DsG4OpBoundaryProcess::PostStepDoIt@442] OpticalSurface name RSOilSurface thePhotonMomentum (eV) 2.88335 theReflectivity 0.0409174 theEfficiency 0. dielectric_metal ground - m1 /dd/Materials/MineralOil m2 /dd/Materials/Acrylic
2016-10-25 20:11:31.336 INFO [3525262] [GSurLib::dump@196] GGeo::loadFromCache GSurLib::dump
0 S( 0 NearPoolCoverSurface) nlv 1 npvp 1 [ obnd 3:Air/NearPoolCoverSurface//PPE]
1 B( 1 NearDeadLinerSurface) nlv 1 npvp 1 [ obnd 13:DeadWater/NearDeadLinerSurface//Tyvek]
2 B( 2 NearOWSLinerSurface) nlv 1 npvp 1 [ ibnd 14:Tyvek//NearOWSLinerSurface/OwsWater]
3 B( 3 NearIWSCurtainSurface) nlv 1 npvp 1 [ ibnd 16:Tyvek//NearIWSCurtainSurface/IwsWater]
4 B( 4 SSTWaterSurfaceNear1) nlv 1 npvp 1 [ obnd 18:IwsWater/SSTWaterSurfaceNear1//StainlessSteel]
5 B( 5 SSTOilSurface) nlv 1 npvp 2 [ ibnd 19:StainlessSteel//SSTOilSurface/MineralOil]
6 S( 6 lvPmtHemiCathodeSensorSurface) nlv 1 npvp 672 [ obnd 29:Vacuum/lvPmtHemiCathodeSensorSurface//Bialkali]
7 S( 7 lvHeadonPmtCathodeSensorSurface) nlv 1 npvp 12 [ obnd 34:Vacuum/lvHeadonPmtCathodeSensorSurface//Bialkali]
8 S( 8 RSOilSurface) nlv 1 npvp 64 [ obnd 37:MineralOil/RSOilSurface//Acrylic] <-- FLIPPED ???
9 B( 9 ESRAirSurfaceTop) nlv 1 npvp 2 [ obnd 39:Air/ESRAirSurfaceTop//ESR]
10 B( 10 ESRAirSurfaceBot) nlv 1 npvp 2 [ obnd 40:Air/ESRAirSurfaceBot//ESR]
11 S( 11 AdCableTraySurface) nlv 1 npvp 2 [ obnd 76:IwsWater/AdCableTraySurface//UnstStainlessSteel]
12 B( 12 SSTWaterSurfaceNear2) nlv 1 npvp 1 [ obnd 80:IwsWater/SSTWaterSurfaceNear2//StainlessSteel]
op --surf 8 ## type 0, is dielectric_metal ... TODO: trace this
2016-10-25 20:29:13.727 INFO [3530462] [GSurfaceLib::dump@717] ( 8, 0, 3,100) GPropertyMap<T>:: 8 surface s: GOpticalSurface type 0 model 1 finish 3 value 1 RSOilSurface k:detect absorb reflect_specular reflect_diffuse extra_x extra_y extra_z extra_w RSOilSurface
domain detect absorb reflect_specular reflect_diffuse extra_x
60 0 0.827 0 0.173 -1
80 0 0.827015 0 0.172985 -1
100 0 0.85649 0 0.14351 -1
120 0 0.885965 0 0.114035 -1
140 0 0.897743 0 0.102257 -1
160 0 0.909501 0 0.0904994 -1
180 0 0.921258 0 0.0787423 -1
200 0 0.933007 0 0.0669933 -1
220 0 0.938282 0 0.0617179 -1
240 0 0.943557 0 0.0564426 -1
260 0 0.947648 0 0.0523518 -1
280 0 0.95055 0 0.0494499 -1
300 0 0.953451 0 0.0465491 -1
320 0 0.954789 0 0.0452105 -1
340 0 0.956128 0 0.043872 -1
Other than perfect additions all surfaces are type=dielectric_metal with finish ground (other than ESRAir.. which is polished)
Looks to be a surface type bug.
Hmm the perfect surfaces listed as finish: polishedfrontpainted
61 enum G4OpticalSurfaceFinish
62 {
63 polished, // smooth perfectly polished surface
64 polishedfrontpainted, // smooth top-layer (front) paint
65 polishedbackpainted, // same is 'polished' but with a back-paint
66
67 ground, // rough surface
68 groundfrontpainted, // rough top-layer (front) paint
69 groundbackpainted, // same as 'ground' but with a back-paint
65 enum G4SurfaceType
66 {
67 dielectric_metal, // dielectric-metal interface
68 dielectric_dielectric, // dielectric-dielectric interface
69 dielectric_LUT, // dielectric-Look-Up-Table interface
70 dielectric_dichroic, // dichroic filter interface
71 firsov, // for Firsov Process
72 x_ray // for x-ray mirror process
73 };
74
75 /////////////////////
76 // Class Definition
77 /////////////////////
78
79 class G4SurfaceProperty
80 {
op --surf
2016-10-25 20:54:23.188 INFO [3537695] [GSurfaceLib::Summary@137] GSurfaceLib::dump NumSurfaces 48 NumFloat4 2
2016-10-25 20:54:23.189 INFO [3537695] [GSurfaceLib::dump@651] (index,type,finish,value)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] NearPoolCoverSurface ( 0, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] NearDeadLinerSurface ( 1, 0, 3, 20)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] NearOWSLinerSurface ( 2, 0, 3, 20)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] NearIWSCurtainSurface ( 3, 0, 3, 20)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] SSTWaterSurfaceNear1 ( 4, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] SSTOilSurface ( 5, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] lvPmtHemiCathodeSensorSurface ( 6, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] lvHeadonPmtCathodeSensorSurface ( 7, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] RSOilSurface ( 8, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] ESRAirSurfaceTop ( 9, 0, 0, 0)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] ESRAirSurfaceBot ( 10, 0, 0, 0)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] AdCableTraySurface ( 11, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] SSTWaterSurfaceNear2 ( 12, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] PmtMtTopRingSurface ( 13, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] PmtMtBaseRingSurface ( 14, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] PmtMtRib1Surface ( 15, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] PmtMtRib2Surface ( 16, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] PmtMtRib3Surface ( 17, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] LegInIWSTubSurface ( 18, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] TablePanelSurface ( 19, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] SupportRib1Surface ( 20, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] SupportRib5Surface ( 21, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] SlopeRib1Surface ( 22, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] SlopeRib5Surface ( 23, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] ADVertiCableTraySurface ( 24, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] ShortParCableTraySurface ( 25, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] NearInnInPiperSurface ( 26, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] NearInnOutPiperSurface ( 27, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] LegInOWSTubSurface ( 28, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] UnistrutRib6Surface ( 29, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] UnistrutRib7Surface ( 30, 0, 3,100)
2016-10-25 20:54:23.189 WARN [3537695] [GSurfaceLib::dump@658] UnistrutRib3Surface ( 31, 0, 3,100)
2016-10-25 20:54:23.190 WARN [3537695] [GSurfaceLib::dump@658] UnistrutRib5Surface ( 32, 0, 3,100)
2016-10-25 20:54:23.190 WARN [3537695] [GSurfaceLib::dump@658] UnistrutRib4Surface ( 33, 0, 3,100)
2016-10-25 20:54:23.190 WARN [3537695] [GSurfaceLib::dump@658] UnistrutRib1Surface ( 34, 0, 3,100)
2016-10-25 20:54:23.190 WARN [3537695] [GSurfaceLib::dump@658] UnistrutRib2Surface ( 35, 0, 3,100)
2016-10-25 20:54:23.190 WARN [3537695] [GSurfaceLib::dump@658] UnistrutRib8Surface ( 36, 0, 3,100)
2016-10-25 20:54:23.190 WARN [3537695] [GSurfaceLib::dump@658] UnistrutRib9Surface ( 37, 0, 3,100)
2016-10-25 20:54:23.190 WARN [3537695] [GSurfaceLib::dump@658] TopShortCableTraySurface ( 38, 0, 3,100)
2016-10-25 20:54:23.190 WARN [3537695] [GSurfaceLib::dump@658] TopCornerCableTraySurface ( 39, 0, 3,100)
2016-10-25 20:54:23.190 WARN [3537695] [GSurfaceLib::dump@658] VertiCableTraySurface ( 40, 0, 3,100)
2016-10-25 20:54:23.190 WARN [3537695] [GSurfaceLib::dump@658] NearOutInPiperSurface ( 41, 0, 3,100)
2016-10-25 20:54:23.190 WARN [3537695] [GSurfaceLib::dump@658] NearOutOutPiperSurface ( 42, 0, 3,100)
2016-10-25 20:54:23.190 WARN [3537695] [GSurfaceLib::dump@658] LegInDeadTubSurface ( 43, 0, 3,100)
2016-10-25 20:54:23.190 WARN [3537695] [GSurfaceLib::dump@658] perfectDetectSurface ( 44, 1, 1,100)
2016-10-25 20:54:23.190 WARN [3537695] [GSurfaceLib::dump@658] perfectAbsorbSurface ( 45, 1, 1,100)
2016-10-25 20:54:23.190 WARN [3537695] [GSurfaceLib::dump@658] perfectSpecularSurface ( 46, 1, 1,100)
2016-10-25 20:54:23.190 WARN [3537695] [GSurfaceLib::dump@658] perfectDiffuseSurface ( 47, 1, 1,100)
248 void G4DAEWriteStructure::
249 OpticalSurfaceWrite(xercesc::DOMElement* targetElement,
250 const G4OpticalSurface* const surf)
251 {
252 xercesc::DOMElement* optElement = NewElement("opticalsurface");
253 G4OpticalSurfaceModel smodel = surf->GetModel();
254 G4double sval = (smodel==glisur) ? surf->GetPolish() : surf->GetSigmaAlpha();
255
256 optElement->setAttributeNode(NewNCNameAttribute("name", surf->GetName()));
257 optElement->setAttributeNode(NewAttribute("model", smodel));
258 optElement->setAttributeNode(NewAttribute("finish", surf->GetFinish()));
259 optElement->setAttributeNode(NewAttribute("type", surf->GetType()));
260 optElement->setAttributeNode(NewAttribute("value", sval));
261
262 G4MaterialPropertiesTable* ptable = surf->GetMaterialPropertiesTable();
263 PropertyWrite( optElement, ptable );
264
265 targetElement->appendChild(optElement);
266 }
delta:ana blyth$ tlaser.py ## apply seqhis selection to pick the most common seqs for A and B
- A:seqhis_ana noname
- 8ccccd 1.000 7673 [6 ] TO BT BT BT BT SA
- 7673 1.00
- B:seqhis_ana noname
- 8c0cc0cd 1.000 7030 [8 ] TO BT ?0? BT BT ?0? BT SA
- 7030 1.00
Huh looks like laser going in different directions:
In [6]: a.rpost_(slice(0,6)) ## heading in some combination of X and Y direction
Out[6]:
A()sliced
A([[[ -18079.4443, -799699.4149, -6604.9499, 0.0998],
[ -17219.8321, -800985.8917, -6604.9499, 7.8266],
[ -17214.1845, -800994.1278, -6604.9499, 7.8765],
[ -16980.2796, -801344.2792, -6604.9499, 9.98 ],
[ -16970.161 , -801359.3395, -6604.9499, 10.0702],
[ -16826.3825, -801575.3603, -6604.9499, 11.3474]],
In [13]: b.rpost_(slice(0,6)) ## huh heading in -Z direction
Out[13]:
A()sliced
A([[[ -18079.4443, -799699.4149, -6604.9499, 0.0998],
[ -18079.4443, -799699.4149, -8635.0278, 10.5229],
[ -18079.4443, -799699.4149, -8650.0881, 10.6008],
[ -18079.4443, -799699.4149, -8850.1073, 11.639 ],
[ -18079.4443, -799699.4149, -8895.0528, 11.8702],
[ -18079.4443, -799699.4149, -9092.013 , 12.8928]],
OKTest --load --vizg4 --cat laser
OKG4Test --load --vizg4 --cat laser
Gensteps are same by construction, suspect CTorchSource not reading it:
In [3]: a.gs
Out[3]:
A(torch,1,laser)-
A([[[ 0. , 0. , 0. , 0. ],
[ -18079.4531, -799699.4375, -6605. , 0.1 ],
[ 0.5556, -0.8314, 0. , 1. ],
[ 0. , 0. , 0. , 430. ],
[ 0. , 1. , 0. , 1. ],
[ 0. , 0. , 0. , 0. ]]], dtype=float32)
In [4]: b.gs
Out[4]:
A(torch,-1,laser)-
A([[[ 0. , 0. , 0. , 0. ],
[ -18079.4531, -799699.4375, -6605. , 0.1 ],
[ 0.5556, -0.8314, 0. , 1. ],
[ 0. , 0. , 0. , 430. ],
[ 0. , 1. , 0. , 1. ],
[ 0. , 0. , 0. , 0. ]]], dtype=float32)
Restricting to top seq:
A:seqhis_ana noname
8ccccd 1.000 7673 [6 ] TO BT BT BT BT SA
7673 1.00
B:seqhis_ana noname
8ccccccd 1.000 7500 [8 ] TO BT BT BT BT BT BT SA
7500 1.00
tlaser- ; tlaser-- --okg4 --compute --dbgseqhis 8ccccccd
In [8]: a.rpost_(slice(0,9))[0]
Out[8]:
A()sliced
A([[ -18079.4443, -799699.4149, -6604.9499, 0.0998],
[ -17219.8321, -800985.8917, -6604.9499, 7.8266],
[ -17214.1845, -800994.1278, -6604.9499, 7.8765],
[ -16980.2796, -801344.2792, -6604.9499, 9.98 ],
[ -16970.161 , -801359.3395, -6604.9499, 10.0702],
[ -16826.3825, -801575.3603, -6604.9499, 11.3474],
[ -16520. , -802110. , -7125. , 0. ], << decompression dummies
[ -16520. , -802110. , -7125. , 0. ],
[ -16520. , -802110. , -7125. , 0. ]])
In [14]: a.ox[:,0] # final position photons, no compression
Out[14]:
A()sliced
A([[ -16826.3945, -801575.375 , -6605. , 11.3472],
[ -16826.3945, -801575.375 , -6605. , 11.3472],
[ -16826.3945, -801575.375 , -6605. , 11.3472],
...,
[ -16826.3945, -801575.375 , -6605. , 11.3472],
[ -16826.3945, -801575.375 , -6605. , 11.3472],
[ -16826.3945, -801575.375 , -6605. , 11.3472]], dtype=float32)
In [9]: b.rpost_(slice(0,9))[0]
Out[9]:
A()sliced
A([[ -18079.4443, -799699.4149, -6604.9499, 0.0998],
[ -17218.1849, -800988.2449, -6604.9499, 8.0587],
[ -17212.7726, -800996.481 , -6604.9499, 8.1104],
[ -16978.1618, -801347.3383, -6604.9499, 10.2771],
[ -16968.2785, -801362.3986, -6604.9499, 10.3705],
[ -16824.2646, -801577.7134, -6604.9499, 11.6829],
[ -16822.6174, -801580.3019, -6604.9499, 11.6985],
[ -16696.9582, -801768.0847, -6604.9499, 12.842 ],
[ -16520. , -802110. , -7125. , 0. ]])
In [15]: b.ox[:,0]
Out[15]:
A()sliced
A([[ -16697.0586, -801768.0625, -6605. , 12.842 ],
[ -16697.0586, -801768.0625, -6605. , 12.842 ],
[ -16697.0586, -801768.0625, -6605. , 12.842 ],
...,
[ -16697.0586, -801768.0625, -6605. , 12.842 ],
[ -16697.0586, -801768.0625, -6605. , 12.842 ],
[ -16697.0586, -801768.0625, -6605. , 12.842 ]], dtype=float32)
In [17]: a.ox[:7500,0] - b.ox[:,0]
Out[17]:
A()sliced
A([[-129.3359, 192.6875, 0. , -1.4948],
[-129.3359, 192.6875, 0. , -1.4948],
[-129.3359, 192.6875, 0. , -1.4948],
...,
[-129.3359, 192.6875, 0. , -1.4948],
[-129.3359, 192.6875, 0. , -1.4948],
[-129.3359, 192.6875, 0. , -1.4948]], dtype=float32)
In [1]: a.rpost_(slice(0,9))[0]
Out[1]:
A()sliced
A([[ -18079.4443, -799699.4149, -6604.9499, 0.0998],
[ -17219.8321, -800985.8917, -6604.9499, 7.8266],
[ -17214.1845, -800994.1278, -6604.9499, 7.8765],
[ -16980.2796, -801344.2792, -6604.9499, 9.98 ],
[ -16970.161 , -801359.3395, -6604.9499, 10.0702],
[ -16826.3825, -801575.3603, -6604.9499, 11.3474],
[ -16520. , -802110. , -7125. , 0. ],
[ -16520. , -802110. , -7125. , 0. ],
[ -16520. , -802110. , -7125. , 0. ]])
In [2]: b.rpost_(slice(0,9))[0]
Out[2]:
A()sliced
A([[ -18079.4443, -799699.4149, -6604.9499, 0.0998],
[ -17218.1849, -800988.2449, -6604.9499, 8.0587],
[ -17212.7726, -800996.481 , -6604.9499, 8.1104],
[ -16978.1618, -801347.3383, -6604.9499, 10.2771],
[ -16968.2785, -801362.3986, -6604.9499, 10.3705],
[ -16824.2646, -801577.7134, -6604.9499, 11.6829],
[ -16520. , -802110. , -7125. , 0. ],
[ -16520. , -802110. , -7125. , 0. ],
[ -16520. , -802110. , -7125. , 0. ]])
In [4]: a.ox[:,0]
Out[4]:
A()sliced
A([[ -16826.3945, -801575.375 , -6605. , 11.3472],
[ -16826.3945, -801575.375 , -6605. , 11.3472],
[ -16826.3945, -801575.375 , -6605. , 11.3472],
...,
[ -16826.3945, -801575.375 , -6605. , 11.3472],
[ -16826.3945, -801575.375 , -6605. , 11.3472],
[ -16826.3945, -801575.375 , -6605. , 11.3472]], dtype=float32)
In [5]: b.ox[:,0]
Out[5]:
A()sliced
A([[ -16824.2129, -801577.8125, -6605. , 11.6829],
[ -16824.2129, -801577.8125, -6605. , 11.6829],
[ -16824.2129, -801577.8125, -6605. , 11.6829],
...,
[ -16824.2129, -801577.8125, -6605. , 11.6829],
[ -16824.2129, -801577.8125, -6605. , 11.6829],
[ -16824.2129, -801577.8125, -6605. , 11.6829]], dtype=float32)
In [8]: a.ox[:,0] - b.ox[:763501,0] ## few mm presumably tesselation effect
Out[8]:
A()sliced
A([[-2.1816, 2.4375, 0. , -0.3357],
[-2.1816, 2.4375, 0. , -0.3357],
[-2.1816, 2.4375, 0. , -0.3357],
...,
[-2.1816, 2.4375, 0. , -0.3357],
[-2.1816, 2.4375, 0. , -0.3357],
[-2.1816, 2.4375, 0. , -0.3357]], dtype=float32)
Time shift is smaller than I recall the groupvel issue being:
In [30]: 0.33/11.
Out[30]: 0.030
134 #define FLAGS(p, s, prd) \
135 { \
136 p.flags.i.x = prd.boundary ; \
137 p.flags.u.y = s.identity.w ; \
138 p.flags.u.w |= s.flag ; \
139 } \
( 37) om: MineralOil os: RSOilSurface is: im: Acrylic
(signed boundaries are 1-based, as 0 means miss : so subtract 1 for the 0-based op --bnd)
GSurLib::pushBorderSurfaces does not list it, so it should be isur/osur duped in order to be relevant in both directions ???
WHAT IS THE CG4 8? just the slot
HUH : -ve boundary corresponds to inward going photons ???
In [21]: a.ox[:,3].view(np.int32)
Out[21]:
A()sliced
A([[ -38, 0, 67305984, 6272],
[ -38, 0, 67305984, 6272],
[ -38, 0, 67305984, 6272],
...,
[ -38, 0, 67305984, 6272],
[ -38, 0, 67305984, 6272],
[ -38, 0, 67305984, 6272]], dtype=int32)
In [22]: b.ox[:,3].view(np.int32)
Out[22]:
A()sliced
A([[ 8, 0, 67305984, 6272],
[ 8, 0, 67305984, 6272],
[ 8, 0, 67305984, 6272],
...,
[ 8, 0, 67305984, 6272],
[ 8, 0, 67305984, 6272],
[ 8, 0, 67305984, 6272]], dtype=int32)
586 void CRecorder::RecordPhoton(const G4Step* step)
587 {
588 // gets called at last step (eg absorption) or when truncated
...
609 target->setUInt(target_record_id, 3, 0, 0, m_slot );
610 target->setUInt(target_record_id, 3, 0, 1, 0u );
611 target->setUInt(target_record_id, 3, 0, 2, m_c4.u );
612 target->setUInt(target_record_id, 3, 0, 3, m_mskhis );
613
z is c4:
309 // initial quadrant
310 uifchar4 c4 ;
311 c4.uchar_.x =
312 ( p.position.x > 0.f ? QX : 0u )
313 |
314 ( p.position.y > 0.f ? QY : 0u )
315 |
316 ( p.position.z > 0.f ? QZ : 0u )
317 ;
318
319 c4.uchar_.y = 2u ; // 3-bytes up for grabs
320 c4.uchar_.z = 3u ;
321 c4.uchar_.w = 4u ;
322
323 p.flags.f.z = c4.f ;
In [28]: a.c4
Out[28]:
rec.array([(0, 2, 3, 4), (0, 2, 3, 4), (0, 2, 3, 4), ..., (0, 2, 3, 4), (0, 2, 3, 4), (0, 2, 3, 4)],
dtype=[('x', 'u1'), ('y', 'u1'), ('z', 'u1'), ('w', 'u1')])
In [29]: b.c4
Out[29]:
rec.array([(0, 2, 3, 4), (0, 2, 3, 4), (0, 2, 3, 4), ..., (0, 2, 3, 4), (0, 2, 3, 4), (0, 2, 3, 4)],
dtype=[('x', 'u1'), ('y', 'u1'), ('z', 'u1'), ('w', 'u1')])
2016-10-02 16:51:37.006 INFO [1411044] [CBorderSurfaceTable::init@21] CBorderSurfaceTable::init nsurf 11
0 NearDeadLinerSurface pv1 /dd/Geometry/Sites/lvNearHallBot#pvNearPoolDead #0 pv2 /dd/Geometry/Pool/lvNearPoolDead#pvNearPoolLiner #0
1 NearOWSLinerSurface pv1 /dd/Geometry/Pool/lvNearPoolLiner#pvNearPoolOWS #0 pv2 /dd/Geometry/Pool/lvNearPoolDead#pvNearPoolLiner #0
2 NearIWSCurtainSurface pv1 /dd/Geometry/Pool/lvNearPoolCurtain#pvNearPoolIWS #0 pv2 /dd/Geometry/Pool/lvNearPoolOWS#pvNearPoolCurtain #0
3 SSTWaterSurfaceNear1 pv1 /dd/Geometry/Pool/lvNearPoolIWS#pvNearADE1 #0 pv2 /dd/Geometry/AD/lvADE#pvSST #0
4 SSTOilSurface pv1 /dd/Geometry/AD/lvSST#pvOIL #0 pv2 /dd/Geometry/AD/lvADE#pvSST #0
5 SSTOilSurface pv1 /dd/Geometry/AD/lvSST#pvOIL #0 pv2 /dd/Geometry/AD/lvADE#pvSST #0
6 ESRAirSurfaceTop pv1 /dd/Geometry/AdDetails/lvTopReflector#pvTopRefGap #0 pv2 /dd/Geometry/AdDetails/lvTopRefGap#pvTopESR #0
7 ESRAirSurfaceTop pv1 /dd/Geometry/AdDetails/lvTopReflector#pvTopRefGap #0 pv2 /dd/Geometry/AdDetails/lvTopRefGap#pvTopESR #0
8 ESRAirSurfaceBot pv1 /dd/Geometry/AdDetails/lvBotReflector#pvBotRefGap #0 pv2 /dd/Geometry/AdDetails/lvBotRefGap#pvBotESR #0
9 ESRAirSurfaceBot pv1 /dd/Geometry/AdDetails/lvBotReflector#pvBotRefGap #0 pv2 /dd/Geometry/AdDetails/lvBotRefGap#pvBotESR #0
10 SSTWaterSurfaceNear2 pv1 /dd/Geometry/Pool/lvNearPoolIWS#pvNearADE2 #0 pv2 /dd/Geometry/AD/lvADE#pvSST #0
2016-10-02 16:51:37.006 INFO [1411044] [CBorderSurfaceTable::dump@47] CGeometryTest CBorderSurfaceTable
2016-10-02 16:51:37.006 INFO [1411044] [CSkinSurfaceTable::init@22] CSkinSurfaceTable::init nsurf 36
0 NearPoolCoverSurface lv NULL
1 lvPmtHemiCathodeSensorSurface lv NULL
2 lvHeadonPmtCathodeSensorSurface lv NULL
3 RSOilSurface lv NULL
4 AdCableTraySurface lv NULL
5 PmtMtTopRingSurface lv NULL
6 PmtMtBaseRingSurface lv NULL
7 PmtMtRib1Surface lv NULL
8 PmtMtRib2Surface lv NULL
9 PmtMtRib3Surface lv NULL
10 LegInIWSTubSurface lv NULL
11 TablePanelSurface lv NULL
12 SupportRib1Surface lv NULL
13 SupportRib5Surface lv NULL
14 SlopeRib1Surface lv NULL
15 SlopeRib5Surface lv NULL
16 ADVertiCableTraySurface lv NULL
17 ShortParCableTraySurface lv NULL
18 NearInnInPiperSurface lv NULL
19 NearInnOutPiperSurface lv NULL
20 LegInOWSTubSurface lv NULL
21 UnistrutRib6Surface lv NULL
22 UnistrutRib7Surface lv NULL
23 UnistrutRib3Surface lv NULL
24 UnistrutRib5Surface lv NULL
25 UnistrutRib4Surface lv NULL
26 UnistrutRib1Surface lv NULL
27 UnistrutRib2Surface lv NULL
28 UnistrutRib8Surface lv NULL
29 UnistrutRib9Surface lv NULL
30 TopShortCableTraySurface lv NULL
31 TopCornerCableTraySurface lv NULL
32 VertiCableTraySurface lv NULL
33 NearOutInPiperSurface lv NULL
34 NearOutOutPiperSurface lv NULL
35 LegInDeadTubSurface lv NULL
Steps looking rather similar now, next issue more BULK_ABSORB AB in CG4 than OK.
A:seqhis_ana 1:laser
8ccccd 0.764 763501 [6 ] TO BT BT BT BT SA
4d 0.056 55825 [2 ] TO AB
cccc9ccccd 0.025 25263 [10] TO BT BT BT BT DR BT BT BT BT
8cccc6d 0.020 19707 [7 ] TO SC BT BT BT BT SA
4ccd 0.013 12576 [4 ] TO BT BT AB
8cccc5d 0.011 11183 [7 ] TO RE BT BT BT BT SA
4ccccd 0.009 8554 [6 ] TO BT BT BT BT AB
45d 0.008 7531 [3 ] TO RE AB
8cccc55d 0.005 5362 [8 ] TO RE RE BT BT BT BT SA
8cc6ccd 0.004 4109 [7 ] TO BT BT SC BT BT SA
455d 0.004 3588 [4 ] TO RE RE AB
86ccccd 0.003 2836 [7 ] TO BT BT BT BT SC SA
cccccc6ccd 0.003 2674 [10] TO BT BT SC BT BT BT BT BT BT
8cccc555d 0.003 2524 [9 ] TO RE RE RE BT BT BT BT SA
8cc5ccd 0.002 2359 [7 ] TO BT BT RE BT BT SA
cacccccc6d 0.002 2210 [10] TO SC BT BT BT BT BT BT SR BT
46d 0.002 2118 [3 ] TO SC AB
cccc6ccccd 0.002 2060 [10] TO BT BT BT BT SC BT BT BT BT
4cccd 0.002 1940 [5 ] TO BT BT BT AB
89ccccd 0.002 1880 [7 ] TO BT BT BT BT DR SA
1000000 1.00
B:seqhis_ana -1:laser
8ccccd 0.813 813472 [6 ] TO BT BT BT BT SA
4d 0.072 71523 [2 ] TO AB
cccc9ccccd 0.027 27170 [10] TO BT BT BT BT DR BT BT BT BT
4ccd 0.017 17386 [4 ] TO BT BT AB
8cccc6d 0.015 15107 [7 ] TO SC BT BT BT BT SA
4ccccd 0.009 8842 [6 ] TO BT BT BT BT AB
cacccccc6d 0.004 3577 [10] TO SC BT BT BT BT BT BT SR BT
8cc6ccd 0.003 3466 [7 ] TO BT BT SC BT BT SA
46d 0.003 2515 [3 ] TO SC AB
86ccccd 0.002 2476 [7 ] TO BT BT BT BT SC SA
cac0ccc6d 0.002 2356 [9 ] TO SC BT BT BT ?0? BT SR BT
cccccc6ccd 0.002 2157 [10] TO BT BT SC BT BT BT BT BT BT
89ccccd 0.002 2127 [7 ] TO BT BT BT BT DR SA
4cccd 0.002 1977 [5 ] TO BT BT BT AB
cccc6ccccd 0.002 1949 [10] TO BT BT BT BT SC BT BT BT BT
8ccccc6d 0.002 1515 [8 ] TO SC BT BT BT BT BT SA
ccbccccc6d 0.001 1429 [10] TO SC BT BT BT BT BT BR BT BT
4cc9ccccd 0.001 1215 [9 ] TO BT BT BT BT DR BT BT AB
4cd 0.001 1077 [3 ] TO BT AB
4cc6d 0.001 802 [5 ] TO SC BT BT AB
1000000 1.00
A:seqhis_ana 1:laser
8ccccd 0.767 7673 [6 ] TO BT BT BT BT SA
4d 0.055 553 [2 ] TO AB
cccc9ccccd 0.024 242 [10] TO BT BT BT BT DR BT BT BT BT
8cccc6d 0.019 188 [7 ] TO SC BT BT BT BT SA
4ccd 0.012 122 [4 ] TO BT BT AB
8cccc5d 0.012 121 [7 ] TO RE BT BT BT BT SA
45d 0.006 65 [3 ] TO RE AB
4ccccd 0.006 63 [6 ] TO BT BT BT BT AB
8cccc55d 0.005 52 [8 ] TO RE RE BT BT BT BT SA
8cc6ccd 0.004 39 [7 ] TO BT BT SC BT BT SA
455d 0.003 34 [4 ] TO RE RE AB
cccccc6ccd 0.003 34 [10] TO BT BT SC BT BT BT BT BT BT
8cc5ccd 0.003 27 [7 ] TO BT BT RE BT BT SA
86ccccd 0.003 27 [7 ] TO BT BT BT BT SC SA
8cccc555d 0.003 26 [9 ] TO RE RE RE BT BT BT BT SA
4cccd 0.003 25 [5 ] TO BT BT BT AB
cacccccc5d 0.002 22 [10] TO RE BT BT BT BT BT BT SR BT
46d 0.002 21 [3 ] TO SC AB
cccc6ccccd 0.002 20 [10] TO BT BT BT BT SC BT BT BT BT
4ccccc5d 0.002 19 [8 ] TO RE BT BT BT BT BT AB
10000 1.00
B:seqhis_ana -1:laser
8ccccccd 0.750 7500 [8 ] TO BT BT BT BT BT BT SA
4d 0.074 741 [2 ] TO AB
cc9ccccccd 0.043 433 [10] TO BT BT BT BT BT BT DR BT BT
cb9ccccccd 0.027 271 [10] TO BT BT BT BT BT BT DR BR BT
4ccd 0.018 175 [4 ] TO BT BT AB
8cccccc6d 0.014 138 [9 ] TO SC BT BT BT BT BT BT SA
4ccccd 0.009 88 [6 ] TO BT BT BT BT AB
4c9ccccccd 0.008 78 [10] TO BT BT BT BT BT BT DR BT AB
4ccccccd 0.007 70 [8 ] TO BT BT BT BT BT BT AB
cacccccc6d 0.004 35 [10] TO SC BT BT BT BT BT BT SR BT
8cc6ccccd 0.003 25 [9 ] TO BT BT BT BT SC BT BT SA
cccc6ccccd 0.002 22 [10] TO BT BT BT BT SC BT BT BT BT
89ccccccd 0.002 22 [9 ] TO BT BT BT BT BT BT DR SA
ccbccccc6d 0.002 22 [10] TO SC BT BT BT BT BT BR BT BT
4cccd 0.002 21 [5 ] TO BT BT BT AB
8cccc6ccd 0.002 21 [9 ] TO BT BT SC BT BT BT BT SA
cac0ccc6d 0.002 21 [9 ] TO SC BT BT BT ?0? BT SR BT
46d 0.002 18 [3 ] TO SC AB
cccccc6ccd 0.002 17 [10] TO BT BT SC BT BT BT BT BT BT
bc9ccccccd 0.002 16 [10] TO BT BT BT BT BT BT DR BT BR
10000 1.00
ipython -i $(which tokg4.py) -- --det laser
/Users/blyth/opticks/ana/tokg4.py --det laser
writing opticks environment to /tmp/blyth/opticks/opticks_env.bash
[2016-10-02 11:10:22,331] p22488 {/Users/blyth/opticks/ana/tokg4.py:25} INFO - tag 1 src torch det laser c2max 2.0
[2016-10-02 11:10:22,397] p22488 {/Users/blyth/opticks/ana/tokg4.py:36} INFO - a : laser/torch/ 1 : 20161002-1106 /tmp/blyth/opticks/evt/laser/torch/1/fdom.npy
[2016-10-02 11:10:22,397] p22488 {/Users/blyth/opticks/ana/tokg4.py:37} INFO - b : laser/torch/ -1 : 20161002-1106 /tmp/blyth/opticks/evt/laser/torch/-1/fdom.npy
A Evt( 1,"torch","laser","laser/torch/ 1 : ", seqs="[]") 20161002-1106 /tmp/blyth/opticks/evt/laser/torch/1
B Evt( -1,"torch","laser","laser/torch/ -1 : ", seqs="[]") 20161002-1106 /tmp/blyth/opticks/evt/laser/torch/-1
A:seqhis_ana 1:laser
8ccccd 0.767 7673 [6 ] TO BT BT BT BT SA
4d 0.055 553 [2 ] TO AB
cccc9ccccd 0.024 242 [10] TO BT BT BT BT DR BT BT BT BT
8cccc6d 0.019 188 [7 ] TO SC BT BT BT BT SA
4ccd 0.012 122 [4 ] TO BT BT AB
8cccc5d 0.012 121 [7 ] TO RE BT BT BT BT SA
45d 0.006 65 [3 ] TO RE AB
4ccccd 0.006 63 [6 ] TO BT BT BT BT AB
8cccc55d 0.005 52 [8 ] TO RE RE BT BT BT BT SA
8cc6ccd 0.004 39 [7 ] TO BT BT SC BT BT SA
455d 0.003 34 [4 ] TO RE RE AB
cccccc6ccd 0.003 34 [10] TO BT BT SC BT BT BT BT BT BT
8cc5ccd 0.003 27 [7 ] TO BT BT RE BT BT SA
86ccccd 0.003 27 [7 ] TO BT BT BT BT SC SA
8cccc555d 0.003 26 [9 ] TO RE RE RE BT BT BT BT SA
4cccd 0.003 25 [5 ] TO BT BT BT AB
cacccccc5d 0.002 22 [10] TO RE BT BT BT BT BT BT SR BT
46d 0.002 21 [3 ] TO SC AB
cccc6ccccd 0.002 20 [10] TO BT BT BT BT SC BT BT BT BT
4ccccc5d 0.002 19 [8 ] TO RE BT BT BT BT BT AB
10000 1.00
B:seqhis_ana -1:laser
8c0cc0cd 0.703 7030 [8 ] TO BT ?0? BT BT ?0? BT SA
4d 0.090 899 [2 ] TO AB
4c9c0cc0cd 0.030 301 [10] TO BT ?0? BT BT ?0? BT DR BT AB
cb9c0cc0cd 0.029 285 [10] TO BT ?0? BT BT ?0? BT DR BR BT
4cc0cd 0.022 217 [6 ] TO BT ?0? BT BT AB
40cd 0.020 201 [4 ] TO BT ?0? AB
8cccccc6d 0.015 152 [9 ] TO SC BT BT BT BT BT BT SA
4c0cc0cd 0.015 145 [8 ] TO BT ?0? BT BT ?0? BT AB
bb9c0cc0cd 0.011 105 [10] TO BT ?0? BT BT ?0? BT DR BR BR
cac0ccc6d 0.005 52 [9 ] TO SC BT BT BT ?0? BT SR BT
46d 0.005 49 [3 ] TO SC AB
cc0b0ccc6d 0.004 44 [10] TO SC BT BT BT ?0? BR ?0? BT BT
cc9c0cc0cd 0.004 43 [10] TO BT ?0? BT BT ?0? BT DR BT BT
cacccccc6d 0.004 40 [10] TO SC BT BT BT BT BT BT SR BT
4c6c0cc0cd 0.004 39 [10] TO BT ?0? BT BT ?0? BT SC BT AB
cccc6cc0cd 0.002 21 [10] TO BT ?0? BT BT SC BT BT BT BT
4cd 0.002 20 [3 ] TO BT AB
bc6c0cc0cd 0.002 17 [10] TO BT ?0? BT BT ?0? BT SC BT BR
c9cccccc6d 0.002 17 [10] TO SC BT BT BT BT BT BT DR BT
cccccccc6d 0.002 17 [10] TO SC BT BT BT BT BT BT BT BT
10000 1.00
A:seqmat_ana 1:laser
443231 0.774 7736 [6 ] Gd Ac LS Ac MO MO
11 0.055 553 [2 ] Gd Gd
4432311 0.031 314 [7 ] Gd Gd Ac LS Ac MO MO
3323443231 0.026 265 [10] Gd Ac LS Ac MO MO Ac LS Ac Ac
2231 0.012 122 [4 ] Gd Ac LS LS
111 0.009 86 [3 ] Gd Gd Gd
44323111 0.007 72 [8 ] Gd Gd Gd Ac LS Ac MO MO
4432231 0.007 71 [7 ] Gd Ac LS LS Ac MO MO
4443231 0.005 46 [7 ] Gd Ac LS Ac MO MO MO
fff3432311 0.004 39 [10] Gd Gd Ac LS Ac MO Ac Ai Ai Ai
3323132231 0.004 39 [10] Gd Ac LS LS Ac Gd Ac LS Ac Ac
1111 0.004 35 [4 ] Gd Gd Gd Gd
4433432311 0.003 33 [10] Gd Gd Ac LS Ac MO Ac Ac MO MO
443231111 0.003 31 [9 ] Gd Gd Gd Gd Ac LS Ac MO MO
aa332311 0.003 26 [8 ] Gd Gd Ac LS Ac Ac ES ES
33231 0.003 25 [5 ] Gd Ac LS Ac Ac
11111 0.002 20 [5 ] Gd Gd Gd Gd Gd
dd432311 0.002 20 [8 ] Gd Gd Ac LS Ac MO Vm Vm
44322231 0.002 17 [8 ] Gd Ac LS LS LS Ac MO MO
331 0.001 14 [3 ] Gd Ac Ac
10000 1.00
B:seqmat_ana -1:laser
44332331 0.718 7175 [8 ] Gd Ac Ac LS Ac Ac MO MO
11 0.090 899 [2 ] Gd Gd
ff44332331 0.034 340 [10] Gd Ac Ac LS Ac Ac MO MO Ai Ai
3444332331 0.026 264 [10] Gd Ac Ac LS Ac Ac MO MO MO Ac
332331 0.022 217 [6 ] Gd Ac Ac LS Ac Ac
3331 0.020 201 [4 ] Gd Ac Ac Ac
443432311 0.015 154 [9 ] Gd Gd Ac LS Ac MO Ac MO MO
4444332331 0.013 134 [10] Gd Ac Ac LS Ac Ac MO MO MO MO
33ff332311 0.005 52 [10] Gd Gd Ac LS Ac Ac Ai Ai Ac Ac
f344332331 0.005 51 [10] Gd Ac Ac LS Ac Ac MO MO Ac Ai
111 0.005 49 [3 ] Gd Gd Gd
3233332311 0.004 43 [10] Gd Gd Ac LS Ac Ac Ac Ac LS Ac
3ff3432311 0.004 40 [10] Gd Gd Ac LS Ac MO Ac Ai Ai Ac
3344332331 0.003 29 [10] Gd Ac Ac LS Ac Ac MO MO Ac Ac
f444332331 0.003 29 [10] Gd Ac Ac LS Ac Ac MO MO MO Ai
331 0.002 20 [3 ] Gd Ac Ac
444332331 0.002 19 [9 ] Gd Ac Ac LS Ac Ac MO MO MO
3443432311 0.002 17 [10] Gd Gd Ac LS Ac MO Ac MO MO Ac
3232332331 0.002 16 [10] Gd Ac Ac LS Ac Ac LS Ac LS Ac
3433432311 0.002 15 [10] Gd Gd Ac LS Ac MO Ac Ac MO Ac
10000 1.00