fixpol beam down onto a cube of scintillator within mineral oil box.
27 tboxlaser--(){
28 local cmdline=$*
29 local tag=$(tboxlaser-tag)
30 local photons=500000
31
32 #local nm=380
33 local nm=480
34
35 local m1=MineralOil
36 #local m2=Pyrex
37 local m2=GdDopedLS
38
39 ## beam in -Z direction, fixpol +Y
40
41 local torch_config=(
42 type=point
43 mode=fixpol
44 polarization=0,1,0
45 source=0,0,299
46 target=0,0,0
47 photons=$photons
48 material=$m1
49 wavelength=$nm
50 weight=1.0
51 time=0.1
52 zenithazimuth=0,1,0,1
53 radius=0
54 )
55 local test_config=(
56 mode=BoxInBox
57 analytic=1
58
59 shape=box
60 boundary=Rock//perfectAbsorbSurface/$m1
61 parameters=0,0,0,300
62
63 shape=box
64 boundary=$m1///$m2
65 parameters=0,0,0,100
66 )
67
68 op.sh \
69 --test --testconfig "$(join _ ${test_config[@]})" \
70 --torch --torchconfig "$(join _ ${torch_config[@]})" \
71 --animtimemax 10 \
72 --timemax 10 \
73 --cat boxlaser --tag $tag --save \
74 --eye 0.5,0.5,0.0 \
75 $*
76 }
77 tboxlaser-args(){ echo --tag $(tboxlaser-tag) --det boxlaser --src torch ; }
78 tboxlaser-a(){ tbox.py $(tboxlaser-args) $* ; }
79 tboxlaser-i(){ ipython -i $(which tbox.py) -- $(tboxlaser-args) $* ; }
80 tboxlaser-t()
81 {
82 tboxlaser-
83 tboxlaser-- --okg4 --compute $*
84 }
Check with constants for properties:
tboxlaser-tfx()
{
tboxlaser-t --fxabconfig 10000 --fxab --fxscconfig 10000 --fxsc --fxreconfig 0.5 --fxre $*
}
[2016-10-26 20:28:48,331] p85486 {/Users/blyth/opticks/ana/tboxlaser.py:40} INFO - a : boxlaser/torch/ 1 : 20161026-2028 /tmp/blyth/opticks/evt/boxlaser/torch/1/fdom.npy
[2016-10-26 20:28:48,331] p85486 {/Users/blyth/opticks/ana/tboxlaser.py:41} INFO - b : boxlaser/torch/ -1 : 20161026-2028 /tmp/blyth/opticks/evt/boxlaser/torch/-1/fdom.npy
seqhis_ana 1:boxlaser -1:boxlaser c2 ab ba
8ccd 443899 443465 0.21 1.00 +- 0.00 1.00 +- 0.00 [4 ] TO BT BT SA
4d 9637 9777 1.01 0.99 +- 0.01 1.01 +- 0.01 [2 ] TO AB
4ccd 8934 8976 0.10 1.00 +- 0.01 1.00 +- 0.01 [4 ] TO BT BT AB
8c6cd 8610 8410 2.35 1.02 +- 0.01 0.98 +- 0.01 [5 ] TO BT SC BT SA
86d 7090 7478 10.33 0.95 +- 0.01 1.05 +- 0.01 [3 ] TO SC SA
86ccd 6426 6918 18.14 0.93 +- 0.01 1.08 +- 0.01 [5 ] TO BT BT SC SA
4cd 4762 4767 0.00 1.00 +- 0.01 1.00 +- 0.01 [3 ] TO BT AB
8c5cd 4383 4381 0.00 1.00 +- 0.02 1.00 +- 0.02 [5 ] TO BT RE BT SA
8cc6d 2051 1728 27.61 1.19 +- 0.03 0.84 +- 0.02 [5 ] TO SC BT BT SA
8cc6ccd 1968 1603 37.31 1.23 +- 0.03 0.81 +- 0.02 [7 ] TO BT BT SC BT BT SA
46ccd 159 261 24.77 0.61 +- 0.05 1.64 +- 0.10 [5 ] TO BT BT SC AB
4c6cd 226 249 1.11 0.91 +- 0.06 1.10 +- 0.07 [5 ] TO BT SC BT AB
46d 195 222 1.75 0.88 +- 0.06 1.14 +- 0.08 [3 ] TO SC AB
86c6cd 172 211 3.97 0.82 +- 0.06 1.23 +- 0.08 [6 ] TO BT SC BT SC SA
866ccd 160 197 3.83 0.81 +- 0.06 1.23 +- 0.09 [6 ] TO BT BT SC SC SA
866d 156 196 4.55 0.80 +- 0.06 1.26 +- 0.09 [4 ] TO SC SC SA
4c5cd 99 109 0.48 0.91 +- 0.09 1.10 +- 0.11 [5 ] TO BT RE BT AB
8cb6cd 64 102 8.70 0.63 +- 0.08 1.59 +- 0.16 [6 ] TO BT SC BR BT SA
8c66cd 93 96 0.05 0.97 +- 0.10 1.03 +- 0.11 [6 ] TO BT SC SC BT SA
86c5cd 81 95 1.11 0.85 +- 0.09 1.17 +- 0.12 [6 ] TO BT RE BT SC SA
500000 500000 4.78
Pretty good, after adopt fixpol:
[2016-10-26 19:20:50,882] p77606 {/Users/blyth/opticks/ana/tboxlaser.py:40} INFO - a : boxlaser/torch/ 1 : 20161026-1920 /tmp/blyth/opticks/evt/boxlaser/torch/1/fdom.npy
[2016-10-26 19:20:50,882] p77606 {/Users/blyth/opticks/ana/tboxlaser.py:41} INFO - b : boxlaser/torch/ -1 : 20161026-1920 /tmp/blyth/opticks/evt/boxlaser/torch/-1/fdom.npy
seqhis_ana 1:boxlaser -1:boxlaser c2 ab ba
8ccd 483237 483157 0.01 1.00 +- 0.00 1.00 +- 0.00 [4 ] TO BT BT SA
4d 4070 4143 0.65 0.98 +- 0.02 1.02 +- 0.02 [2 ] TO AB
4ccd 3953 3964 0.02 1.00 +- 0.02 1.00 +- 0.02 [4 ] TO BT BT AB
4cd 3121 3080 0.27 1.01 +- 0.02 0.99 +- 0.02 [3 ] TO BT AB
8c6cd 1622 1608 0.06 1.01 +- 0.03 0.99 +- 0.02 [5 ] TO BT SC BT SA
86ccd 1264 1375 4.67 0.92 +- 0.03 1.09 +- 0.03 [5 ] TO BT BT SC SA
86d 1303 1327 0.22 0.98 +- 0.03 1.02 +- 0.03 [3 ] TO SC SA
8cc6d 392 338 3.99 1.16 +- 0.06 0.86 +- 0.05 [5 ] TO SC BT BT SA
8cc6ccd 375 313 5.59 1.20 +- 0.06 0.83 +- 0.05 [7 ] TO BT BT SC BT BT SA
8c5cd 245 369 25.04 0.66 +- 0.04 1.51 +- 0.08 [5 ] TO BT RE BT SA
8c55cd 72 82 0.65 0.88 +- 0.10 1.14 +- 0.13 [6 ] TO BT RE RE BT SA
45cd 56 20 17.05 2.80 +- 0.37 0.36 +- 0.08 [4 ] TO BT RE AB
8c555cd 52 9 30.31 5.78 +- 0.80 0.17 +- 0.06 [7 ] TO BT RE RE RE BT SA
8bd 32 40 0.89 0.80 +- 0.14 1.25 +- 0.20 [3 ] TO BR SA
8cbcd 35 25 1.67 1.40 +- 0.24 0.71 +- 0.14 [5 ] TO BT BR BT SA
46ccd 12 24 4.00 0.50 +- 0.14 2.00 +- 0.41 [5 ] TO BT BT SC AB
8c5555cd 21 1 0.00 21.00 +- 4.58 0.05 +- 0.05 [8 ] TO BT RE RE RE RE BT SA
4c6cd 19 20 0.03 0.95 +- 0.22 1.05 +- 0.24 [5 ] TO BT SC BT AB
8cb6cd 14 18 0.50 0.78 +- 0.21 1.29 +- 0.30 [6 ] TO BT SC BR BT SA
455cd 17 4 0.00 4.25 +- 1.03 0.24 +- 0.12 [5 ] TO BT RE RE AB
500000 500000 5.31
seqhis_ana_1 1:boxlaser -1:boxlaser c2 ab ba
d 500000 500000 0.00 1.00 +- 0.00 1.00 +- 0.00 [1 ] TO
500000 500000 0.00
seqhis_ana_2 1:boxlaser -1:boxlaser c2 ab ba
cd 494185 494126 0.00 1.00 +- 0.00 1.00 +- 0.00 [2 ] TO BT
4d 4070 4143 0.65 0.98 +- 0.02 1.02 +- 0.02 [2 ] TO AB
6d 1713 1691 0.14 1.01 +- 0.02 0.99 +- 0.02 [2 ] TO SC
bd 32 40 0.89 0.80 +- 0.14 1.25 +- 0.20 [2 ] TO BR
500000 500000 0.42
seqhis_ana_3 1:boxlaser -1:boxlaser c2 ab ba
ccd 488862 488848 0.00 1.00 +- 0.00 1.00 +- 0.00 [3 ] TO BT BT
4d 4070 4143 0.65 0.98 +- 0.02 1.02 +- 0.02 [2 ] TO AB
4cd 3121 3080 0.27 1.01 +- 0.02 0.99 +- 0.02 [3 ] TO BT AB
6cd 1667 1670 0.00 1.00 +- 0.02 1.00 +- 0.02 [3 ] TO BT SC
86d 1303 1327 0.22 0.98 +- 0.03 1.02 +- 0.03 [3 ] TO SC SA
5cd 499 503 0.02 0.99 +- 0.04 1.01 +- 0.04 [3 ] TO BT RE
c6d 401 345 4.20 1.16 +- 0.06 0.86 +- 0.05 [3 ] TO SC BT
8bd 32 40 0.89 0.80 +- 0.14 1.25 +- 0.20 [3 ] TO BR SA
bcd 36 25 1.98 1.44 +- 0.24 0.69 +- 0.14 [3 ] TO BT BR
46d 7 12 0.00 0.58 +- 0.22 1.71 +- 0.49 [3 ] TO SC AB
66d 2 7 0.00 0.29 +- 0.20 3.50 +- 1.32 [3 ] TO SC SC
500000 500000 0.91
seqhis_ana_4 1:boxlaser -1:boxlaser c2 ab ba
8ccd 483237 483157 0.01 1.00 +- 0.00 1.00 +- 0.00 [4 ] TO BT BT SA
4d 4070 4143 0.65 0.98 +- 0.02 1.02 +- 0.02 [2 ] TO AB
4ccd 3953 3964 0.02 1.00 +- 0.02 1.00 +- 0.02 [4 ] TO BT BT AB
4cd 3121 3080 0.27 1.01 +- 0.02 0.99 +- 0.02 [3 ] TO BT AB
6ccd 1671 1727 0.92 0.97 +- 0.02 1.03 +- 0.02 [4 ] TO BT BT SC
c6cd 1646 1638 0.02 1.00 +- 0.02 1.00 +- 0.02 [4 ] TO BT SC BT
86d 1303 1327 0.22 0.98 +- 0.03 1.02 +- 0.03 [3 ] TO SC SA
cc6d 398 342 4.24 1.16 +- 0.06 0.86 +- 0.05 [4 ] TO SC BT BT
c5cd 247 377 27.08 0.66 +- 0.04 1.53 +- 0.08 [4 ] TO BT RE BT
55cd 192 101 28.26 1.90 +- 0.14 0.53 +- 0.05 [4 ] TO BT RE RE
45cd 56 20 17.05 2.80 +- 0.37 0.36 +- 0.08 [4 ] TO BT RE AB
8bd 32 40 0.89 0.80 +- 0.14 1.25 +- 0.20 [3 ] TO BR SA
cbcd 35 25 1.67 1.40 +- 0.24 0.71 +- 0.14 [4 ] TO BT BR BT
b6cd 14 18 0.50 0.78 +- 0.21 1.29 +- 0.30 [4 ] TO BT SC BR
46d 7 12 0.00 0.58 +- 0.22 1.71 +- 0.49 [3 ] TO SC AB
46cd 4 11 0.00 0.36 +- 0.18 2.75 +- 0.83 [4 ] TO BT SC AB
b5cd 2 5 0.00 0.40 +- 0.28 2.50 +- 1.12 [4 ] TO BT RE BR
866d 2 5 0.00 0.40 +- 0.28 2.50 +- 1.12 [4 ] TO SC SC SA
66cd 2 3 0.00 0.67 +- 0.47 1.50 +- 0.87 [4 ] TO BT SC SC
4c6d 0 3 0.00 0.00 +- 0.00 0.00 +- 0.00 [4 ] TO SC BT AB
500000 500000 5.84
seqhis_ana_5 1:boxlaser -1:boxlaser c2 ab ba
8ccd 483237 483157 0.01 1.00 +- 0.00 1.00 +- 0.00 [4 ] TO BT BT SA
4d 4070 4143 0.65 0.98 +- 0.02 1.02 +- 0.02 [2 ] TO AB
4ccd 3953 3964 0.02 1.00 +- 0.02 1.00 +- 0.02 [4 ] TO BT BT AB
4cd 3121 3080 0.27 1.01 +- 0.02 0.99 +- 0.02 [3 ] TO BT AB
8c6cd 1622 1608 0.06 1.01 +- 0.03 0.99 +- 0.02 [5 ] TO BT SC BT SA
86ccd 1264 1375 4.67 0.92 +- 0.03 1.09 +- 0.03 [5 ] TO BT BT SC SA
86d 1303 1327 0.22 0.98 +- 0.03 1.02 +- 0.03 [3 ] TO SC SA
8cc6d 392 338 3.99 1.16 +- 0.06 0.86 +- 0.05 [5 ] TO SC BT BT SA
c6ccd 384 320 5.82 1.20 +- 0.06 0.83 +- 0.05 [5 ] TO BT BT SC BT
8c5cd 245 369 25.04 0.66 +- 0.04 1.51 +- 0.08 [5 ] TO BT RE BT SA
555cd 98 13 65.09 7.54 +- 0.76 0.13 +- 0.04 [5 ] TO BT RE RE RE
c55cd 76 83 0.31 0.92 +- 0.11 1.09 +- 0.12 [5 ] TO BT RE RE BT
45cd 56 20 17.05 2.80 +- 0.37 0.36 +- 0.08 [4 ] TO BT RE AB
8bd 32 40 0.89 0.80 +- 0.14 1.25 +- 0.20 [3 ] TO BR SA
8cbcd 35 25 1.67 1.40 +- 0.24 0.71 +- 0.14 [5 ] TO BT BR BT SA
46ccd 12 24 4.00 0.50 +- 0.14 2.00 +- 0.41 [5 ] TO BT BT SC AB
4c6cd 19 20 0.03 0.95 +- 0.22 1.05 +- 0.24 [5 ] TO BT SC BT AB
cb6cd 14 18 0.50 0.78 +- 0.21 1.29 +- 0.30 [5 ] TO BT SC BR BT
455cd 17 4 0.00 4.25 +- 1.03 0.24 +- 0.12 [5 ] TO BT RE RE AB
46d 7 12 0.00 0.58 +- 0.22 1.71 +- 0.49 [3 ] TO SC AB
500000 500000 7.24