OKG4Test --compute --save --dbgseqhis 8ccccccd
8ccccd 0.439 4386 [6 ] TO BT BT BT BT SA
8ccccccd 0.424 4240 [8 ] TO BT BT BT BT BT BT SA <<< CG4 recording has extra "BT BT" cf Opticks
2016-10-01 18:25:39.730 INFO [1273721] [CRecorder::Dump@670] CSteppingAction::UserSteppingAction DONE record_id 6893
2016-10-01 18:25:39.730 INFO [1273721] [CRecorder::Dump@673] seqhis 8ccccccd TORCH BOUNDARY_TRANSMIT BOUNDARY_TRANSMIT BOUNDARY_TRANSMIT BOUNDARY_TRANSMIT BOUNDARY_TRANSMIT BOUNDARY_TRANSMIT SURFACE_ABSORB . . . . . . . .
2016-10-01 18:25:39.731 INFO [1273721] [CRecorder::Dump@677] seqmat 44343231 GdDopedLS Acrylic LiquidScintillator Acrylic MineralOil Acrylic MineralOil MineralOil - - - - - - - -
TODO: get seqmat to include the Stainless Steel on which SURFACE_ABSORB happens
Investigated issue in tlaser and seqvol : volume sequence indexing. Next: rather different AB (BULK_ABSORB).
delta:opticksnpy blyth$ tokg4.py
/Users/blyth/opticks/ana/tokg4.py
writing opticks environment to /tmp/blyth/opticks/opticks_env.bash
[2016-10-02 18:35:49,636] p48179 {/Users/blyth/opticks/ana/tokg4.py:25} INFO - tag 1 src torch det dayabay c2max 2.0
[2016-10-02 18:35:50,057] p48179 {/Users/blyth/opticks/ana/tokg4.py:39} INFO - a : dayabay/torch/ 1 : 20161002-1835 /tmp/blyth/opticks/evt/dayabay/torch/1/fdom.npy
[2016-10-02 18:35:50,057] p48179 {/Users/blyth/opticks/ana/tokg4.py:40} INFO - b : dayabay/torch/ -1 : 20161002-1835 /tmp/blyth/opticks/evt/dayabay/torch/-1/fdom.npy
A Evt( 1,"torch","dayabay","dayabay/torch/ 1 : ", seqs="[]") 20161002-1835 /tmp/blyth/opticks/evt/dayabay/torch/1
B Evt( -1,"torch","dayabay","dayabay/torch/ -1 : ", seqs="[]") 20161002-1835 /tmp/blyth/opticks/evt/dayabay/torch/-1
A:seqhis_ana 1:dayabay
8ccccd 0.434 43405 [6 ] TO BT BT BT BT SA
ccaccccccd 0.090 9009 [10] TO BT BT BT BT BT BT SR BT BT
4cccccd 0.061 6104 [7 ] TO BT BT BT BT BT AB
4d 0.061 6051 [2 ] TO AB
cccbcccccd 0.038 3822 [10] TO BT BT BT BT BT BR BT BT BT
8cccccd 0.030 2978 [7 ] TO BT BT BT BT BT SA
8cbcccccd 0.025 2511 [9 ] TO BT BT BT BT BT BR BT SA
8cccc6d 0.022 2165 [7 ] TO SC BT BT BT BT SA
8ccccccd 0.016 1565 [8 ] TO BT BT BT BT BT BT SA
4ccd 0.013 1347 [4 ] TO BT BT AB
cccc9ccccd 0.013 1326 [10] TO BT BT BT BT DR BT BT BT BT
8cccc5d 0.011 1126 [7 ] TO RE BT BT BT BT SA
cccacccccd 0.011 1103 [10] TO BT BT BT BT BT SR BT BT BT
accccccccd 0.009 861 [10] TO BT BT BT BT BT BT BT BT SR
45d 0.008 833 [3 ] TO RE AB
4ccccd 0.007 701 [6 ] TO BT BT BT BT AB
7cccccd 0.006 636 [7 ] TO BT BT BT BT BT SD
8cccc55d 0.006 609 [8 ] TO RE RE BT BT BT BT SA
4cccccccd 0.005 542 [9 ] TO BT BT BT BT BT BT BT AB
cccccccccd 0.004 450 [10] TO BT BT BT BT BT BT BT BT BT
100000 1.00
B:seqhis_ana -1:dayabay
8ccccd 0.461 46136 [6 ] TO BT BT BT BT SA
ccaccccccd 0.081 8088 [10] TO BT BT BT BT BT BT SR BT BT
4d 0.077 7687 [2 ] TO AB
c0cac0cccd 0.067 6679 [10] TO BT BT BT ?0? BT SR BT ?0? BT
cccbcccccd 0.046 4560 [10] TO BT BT BT BT BT BR BT BT BT
8cccccd 0.045 4459 [7 ] TO BT BT BT BT BT SA
8cbcccccd 0.023 2320 [9 ] TO BT BT BT BT BT BR BT SA
4ccd 0.019 1939 [4 ] TO BT BT AB
8cccc6d 0.018 1810 [7 ] TO SC BT BT BT BT SA
cccc9ccccd 0.014 1363 [10] TO BT BT BT BT DR BT BT BT BT
cac00cc0cd 0.013 1284 [10] TO BT ?0? BT BT ?0? ?0? BT SR BT
8ccccccd 0.011 1129 [8 ] TO BT BT BT BT BT BT SA
7cccccd 0.011 1110 [7 ] TO BT BT BT BT BT SD
abaccccccd 0.008 847 [10] TO BT BT BT BT BT BT SR BR SR
4ccccd 0.007 731 [6 ] TO BT BT BT BT AB
4cccccd 0.006 553 [7 ] TO BT BT BT BT BT AB
ccc0b0cccd 0.005 473 [10] TO BT BT BT ?0? BR ?0? BT BT BT
cabac0cccd 0.004 388 [10] TO BT BT BT ?0? BT SR BR SR BT
8c0cccd 0.003 308 [7 ] TO BT BT BT ?0? BT SA
cacccccc6d 0.003 300 [10] TO SC BT BT BT BT BT BT SR BT
100000 1.00
Associated to Acrylic to Acrylic steps between OAV and reflector
Problem with volume sequencing is the large number of volumes and expensive storage of int32*10 sequence but actually the number relevant to critical optical path is not so big, so judicious favoriting of 15 volumes 0x0->0xE specific to targetted AD and overflow 0xF for all others
OKG4Test --compute --save
OKG4Test --compute --save --steppingdbg ## very verbose
OKG4Test --compute --save --dbgseqhis c0cac0cccd ## just dump particular history photons, with some zero flags
lldb OKG4Test -- --compute --save
(lldb) b "OpRayleigh::GetMeanFreePath(G4Track const&, double, G4ForceCondition*)"
tokg4.py
ipython -i $(which tokg4.py)
OKTest --load # works
OKTest --load --vizg4 # was failing wrt photon buffer
OKG4Test --load --vizg4 # succeeds to load g4evt, fix index loading with indexPresentationPrep, but suspect using Opticks index with G4 evt
This was using old propagator upload, fixed by move to viz->uploadEvent following OKG4Test:
2016-09-30 20:44:59.861 INFO [1132529] [OpticksViz::uploadEvent@281] OpticksViz::uploadEvent (0) DONE
2016-09-30 20:44:59.861 INFO [1132529] [OEvent::createBuffers@62] OEvent::createBuffers genstep 1,6,4 nopstep NULL photon 100000,4,4 record 100000,10,2,4 phosel 100000,1,4 recsel 100000,10,1,4 sequence 100000,1,2 seed 0,1,1 hit 897,4,4
2016-09-30 20:44:59.861 FATAL [1132529] [OContext::createBuffer@423] OContext::createBuffer CANNOT createBufferFromGLBO as not uploaded name photon buffer_id -1
Assertion failed: (buffer_id > -1), function createBuffer, file /Users/blyth/opticks/optixrap/OContext.cc, line 427.
Abort trap: 6
2016-09-30 20:47:25.236 INFO [1133593] [OpticksEvent::importParameters@1593] OpticksEvent::importParameters mode_ COMPUTE_MODE --> COMPUTE_MODE
2016-09-30 20:47:25.238 WARN [1133593] [*Index::load@370] Index::load FAILED to load index idpath /tmp/blyth/opticks/evt/dayabay/torch/-1 itemtype Boundary_Index Source path /tmp/blyth/opticks/evt/dayabay/torch/-1/Boundary_IndexSource.json Local path /tmp/blyth/opticks/evt/dayabay/torch/-1/Boundary_IndexLocal.json
This is torch running with a point source and big bugs, position and polarization were wrong
A:seqhis_ana 1:dayabay
8ccccd 0.434 43405 [6 ] TO BT BT BT BT SA
ccaccccccd 0.090 9009 [10] TO BT BT BT BT BT BT SR BT BT
4cccccd 0.061 6104 [7 ] TO BT BT BT BT BT AB
4d 0.061 6051 [2 ] TO AB
cccbcccccd 0.038 3822 [10] TO BT BT BT BT BT BR BT BT BT
8cccccd 0.030 2978 [7 ] TO BT BT BT BT BT SA
8cbcccccd 0.025 2511 [9 ] TO BT BT BT BT BT BR BT SA
8cccc6d 0.022 2165 [7 ] TO SC BT BT BT BT SA
8ccccccd 0.016 1565 [8 ] TO BT BT BT BT BT BT SA
4ccd 0.013 1347 [4 ] TO BT BT AB
cccc9ccccd 0.013 1326 [10] TO BT BT BT BT DR BT BT BT BT
8cccc5d 0.011 1126 [7 ] TO RE BT BT BT BT SA
cccacccccd 0.011 1103 [10] TO BT BT BT BT BT SR BT BT BT
accccccccd 0.009 861 [10] TO BT BT BT BT BT BT BT BT SR
45d 0.008 833 [3 ] TO RE AB
4ccccd 0.007 701 [6 ] TO BT BT BT BT AB
7cccccd 0.006 636 [7 ] TO BT BT BT BT BT SD
8cccc55d 0.006 609 [8 ] TO RE RE BT BT BT BT SA
4cccccccd 0.005 542 [9 ] TO BT BT BT BT BT BT BT AB
cccccccccd 0.004 450 [10] TO BT BT BT BT BT BT BT BT BT
100000 1.00
B:seqhis_ana -1:dayabay
666666666d 0.141 14077 [10] TO SC SC SC SC SC SC SC SC SC
4d 0.084 8431 [2 ] TO AB
46d 0.072 7202 [3 ] TO SC AB
8d 0.070 6966 [2 ] TO SA
866d 0.065 6510 [4 ] TO SC SC SA
8666d 0.063 6344 [5 ] TO SC SC SC SA
86d 0.061 6118 [3 ] TO SC SA
466d 0.059 5900 [4 ] TO SC SC AB
86666d 0.056 5624 [6 ] TO SC SC SC SC SA
866666d 0.051 5056 [7 ] TO SC SC SC SC SC SA
4666d 0.049 4887 [5 ] TO SC SC SC AB
8666666d 0.043 4257 [8 ] TO SC SC SC SC SC SC SA
46666d 0.038 3785 [6 ] TO SC SC SC SC AB
86666666d 0.034 3375 [9 ] TO SC SC SC SC SC SC SC SA
466666d 0.030 2979 [7 ] TO SC SC SC SC SC AB
866666666d 0.028 2802 [10] TO SC SC SC SC SC SC SC SC SA
4666666d 0.024 2412 [8 ] TO SC SC SC SC SC SC AB
46666666d 0.018 1847 [9 ] TO SC SC SC SC SC SC SC AB
466666666d 0.014 1381 [10] TO SC SC SC SC SC SC SC SC AB
666666cc6d 0.000 4 [10] TO SC BT BT SC SC SC SC SC SC
Now at least in same ballpark, some zero flags to identify:
A:seqhis_ana 1:dayabay
8ccccd 0.434 43405 [6 ] TO BT BT BT BT SA
ccaccccccd 0.090 9009 [10] TO BT BT BT BT BT BT SR BT BT
4cccccd 0.061 6104 [7 ] TO BT BT BT BT BT AB
4d 0.061 6051 [2 ] TO AB
cccbcccccd 0.038 3822 [10] TO BT BT BT BT BT BR BT BT BT
8cccccd 0.030 2978 [7 ] TO BT BT BT BT BT SA
8cbcccccd 0.025 2511 [9 ] TO BT BT BT BT BT BR BT SA
8cccc6d 0.022 2165 [7 ] TO SC BT BT BT BT SA
8ccccccd 0.016 1565 [8 ] TO BT BT BT BT BT BT SA
4ccd 0.013 1347 [4 ] TO BT BT AB
cccc9ccccd 0.013 1326 [10] TO BT BT BT BT DR BT BT BT BT
8cccc5d 0.011 1126 [7 ] TO RE BT BT BT BT SA
cccacccccd 0.011 1103 [10] TO BT BT BT BT BT SR BT BT BT
accccccccd 0.009 861 [10] TO BT BT BT BT BT BT BT BT SR
45d 0.008 833 [3 ] TO RE AB
4ccccd 0.007 701 [6 ] TO BT BT BT BT AB
7cccccd 0.006 636 [7 ] TO BT BT BT BT BT SD
8cccc55d 0.006 609 [8 ] TO RE RE BT BT BT BT SA
4cccccccd 0.005 542 [9 ] TO BT BT BT BT BT BT BT AB
cccccccccd 0.004 450 [10] TO BT BT BT BT BT BT BT BT BT
100000 1.00
B:seqhis_ana -1:dayabay
8ccccccd 0.423 42297 [8 ] TO BT BT BT BT BT BT SA
cc0ccccccd 0.115 11515 [10] TO BT BT BT BT BT BT ?0? BT BT
4d 0.076 7584 [2 ] TO AB
c0c0c0cccd 0.067 6652 [10] TO BT BT BT ?0? BT ?0? BT ?0? BT
cccbcccccd 0.062 6219 [10] TO BT BT BT BT BT BR BT BT BT
cccccccccd 0.031 3065 [10] TO BT BT BT BT BT BT BT BT BT
ccbccccccd 0.023 2309 [10] TO BT BT BT BT BT BT BR BT BT
4ccd 0.019 1902 [4 ] TO BT BT AB
8cccccc6d 0.016 1552 [9 ] TO SC BT BT BT BT BT BT SA
c0c00cc0cd 0.013 1325 [10] TO BT ?0? BT BT ?0? ?0? BT ?0? BT
8ccccd 0.012 1167 [6 ] TO BT BT BT BT SA
b0ccccccd 0.008 815 [9 ] TO BT BT BT BT BT BT ?0? BR
8cbcccccd 0.007 681 [9 ] TO BT BT BT BT BT BR BT SA
4ccccd 0.007 674 [6 ] TO BT BT BT BT AB
4cccccd 0.006 570 [7 ] TO BT BT BT BT BT AB
ccc0b0cccd 0.005 481 [10] TO BT BT BT ?0? BR ?0? BT BT BT
4ccccccd 0.005 480 [8 ] TO BT BT BT BT BT BT AB
c0cccccc6d 0.004 440 [10] TO SC BT BT BT BT BT BT ?0? BT
c0b0c0cccd 0.004 385 [10] TO BT BT BT ?0? BT ?0? BR ?0? BT
cbcccccccd 0.004 384 [10] TO BT BT BT BT BT BT BT BR BT
After identify SR SURFACE_SREFLECT with SpikeReflection eliminate some zero flags, suspect remainder due to SameMaterial steps:
B:seqhis_ana -1:dayabay
8ccccccd 0.420 419905 [8 ] TO BT BT BT BT BT BT SA
ccaccccccd 0.081 81049 [10] TO BT BT BT BT BT BT SR BT BT
4d 0.078 77610 [2 ] TO AB
c0cac0cccd 0.066 66482 [10] TO BT BT BT ?0? BT SR BT ?0? BT
cccbcccccd 0.063 63079 [10] TO BT BT BT BT BT BR BT BT BT
cc9ccccccd 0.034 33940 [10] TO BT BT BT BT BT BT DR BT BT
cccccccccd 0.031 30534 [10] TO BT BT BT BT BT BT BT BT BT
ccbccccccd 0.023 23444 [10] TO BT BT BT BT BT BT BR BT BT
4ccd 0.019 19127 [4 ] TO BT BT AB
8cccccc6d 0.015 15140 [9 ] TO SC BT BT BT BT BT BT SA
cac00cc0cd 0.013 12771 [10] TO BT ?0? BT BT ?0? ?0? BT SR BT
8ccccd 0.012 12083 [6 ] TO BT BT BT BT SA
abaccccccd 0.008 8032 [10] TO BT BT BT BT BT BT SR BR SR
8cbcccccd 0.008 7512 [9 ] TO BT BT BT BT BT BR BT SA
4ccccd 0.007 7050 [6 ] TO BT BT BT BT AB
4cccccd 0.006 5645 [7 ] TO BT BT BT BT BT AB
ccc0b0cccd 0.005 4873 [10] TO BT BT BT ?0? BR ?0? BT BT BT
4ccccccd 0.005 4847 [8 ] TO BT BT BT BT BT BT AB
cbcccccccd 0.004 3735 [10] TO BT BT BT BT BT BT BT BR BT
cabac0cccd 0.004 3725 [10] TO BT BT BT ?0? BT SR BR SR BT
1000000 1.00
simon:~ blyth$ find /usr/local/opticks -name abbrev.json
/usr/local/opticks/opticksdata/export/DayaBay/GMaterialLib/abbrev.json
/usr/local/opticks/opticksdata/resource/GFlags/abbrev.json
simon:~ blyth$ cat /usr/local/opticks/opticksdata/resource/GFlags/abbrev.json
{
"CERENKOV":"CK",
"SCINTILLATION":"SI",
"TORCH":"TO",
"MISS":"MI",
"BULK_ABSORB":"AB",
"BULK_REEMIT":"RE",
"BULK_SCATTER":"SC",
"SURFACE_DETECT":"SD",
"SURFACE_ABSORB":"SA",
"SURFACE_DREFLECT":"DR",
"SURFACE_SREFLECT":"SR",
"BOUNDARY_REFLECT":"BR",
"BOUNDARY_TRANSMIT":"BT",
"NAN_ABORT":"NA"
}
OKG4Test --compute --save --dbgseqhis c0cac0cccd ## just dump particular history photons, with some zero flags
----CSteppingAction----
2016-10-01 17:53:24.033 INFO [1267418] [CRecorder::Dump@670] CSteppingAction::UserSteppingAction DONE record_id 68
2016-10-01 17:53:24.033 INFO [1267418] [CRecorder::Dump@673] seqhis c0cac0cccd TORCH BOUNDARY_TRANSMIT BOUNDARY_TRANSMIT BOUNDARY_TRANSMIT . BOUNDARY_TRANSMIT SURFACE_SREFLECT BOUNDARY_TRANSMIT . BOUNDARY_TRANSMIT . . . . . .
2016-10-01 17:53:24.033 INFO [1267418] [CRecorder::Dump@677] seqmat 233ff33231 GdDopedLS Acrylic LiquidScintillator Acrylic Acrylic Air Air Acrylic Acrylic LiquidScintillator - - - - - -
2016-10-01 17:53:24.033 INFO [1267418] [Rec::Dump@226] CSteppingAction::UserSteppingAction (Rec)DONE nstates 11
[ 0/ 11]
stepStatus [ Undefined/ GeomBoundary]
flag [ TORCH/ BOUNDARY_TRANSMIT]
bs pri/cur [ Und/ FrT]
material [ GdDopedLS/ Acrylic]
(rec state ;opticalphoton stepNum 11(tk ;opticalphoton tid 69 pid 0 nm 430 mm ori[ -1.81e+04 -8e+05-6.60e+03] pos[ -7791.08e+03-2.04e+03] )
pre d/Geometry/AD/lvIAV#pvGDS rials/GdDopedLS noProc Undefined pos[ 0 0 0] dir[ -0.238 0.33 -0.913] pol[ 0.971 0.081 -0.224] ns 0.100 nm 430.000
post d/Geometry/AD/lvLSO#pvIAV terials/Acrylic Transportation GeomBoundary pos[ -530 734-2.03e+03] dir[ -0.237 0.329 -0.914] pol[ 0.971 0.0807 -0.223] ns 11.512 nm 430.000
)
[ 1/ 11]
stepStatus [ GeomBoundary/ GeomBoundary]
flag [ BOUNDARY_TRANSMIT/ BOUNDARY_TRANSMIT]
bs pri/cur [ FrT/ FrT]
material [ Acrylic/ LiquidScintillator]
(rec state ;opticalphoton stepNum 11(tk ;opticalphoton tid 69 pid 0 nm 430 mm ori[ -1.81e+04 -8e+05-6.60e+03] pos[ -7791.08e+03-2.04e+03] )
pre d/Geometry/AD/lvLSO#pvIAV terials/Acrylic Transportation GeomBoundary pos[ -530 734-2.03e+03] dir[ -0.237 0.329 -0.914] pol[ 0.971 0.0807 -0.223] ns 11.512 nm 430.000
post d/Geometry/AD/lvOAV#pvLSO uidScintillator Transportation GeomBoundary pos[ -534 739-2.04e+03] dir[ -0.238 0.33 -0.913] pol[ 0.971 0.081 -0.224] ns 11.597 nm 430.000
)
[ 2/ 11]
stepStatus [ GeomBoundary/ GeomBoundary]
flag [ BOUNDARY_TRANSMIT/ BOUNDARY_TRANSMIT]
bs pri/cur [ FrT/ FrT]
material [ LiquidScintillator/ Acrylic]
(rec state ;opticalphoton stepNum 11(tk ;opticalphoton tid 69 pid 0 nm 430 mm ori[ -1.81e+04 -8e+05-6.60e+03] pos[ -7791.08e+03-2.04e+03] )
pre d/Geometry/AD/lvOAV#pvLSO uidScintillator Transportation GeomBoundary pos[ -534 739-2.04e+03] dir[ -0.238 0.33 -0.913] pol[ 0.971 0.081 -0.224] ns 11.597 nm 430.000
post d/Geometry/AD/lvOIL#pvOAV terials/Acrylic Transportation GeomBoundary pos[ -649 899-2.49e+03] dir[ -0.237 0.329 -0.914] pol[ 0.971 0.0807 -0.223] ns 14.082 nm 430.000
)
[ 3/ 11]
stepStatus [ GeomBoundary/ GeomBoundary]
flag [ BOUNDARY_TRANSMIT/ .]
bs pri/cur [ FrT/ SAM]
material [ Acrylic/ Acrylic]
(rec state ;opticalphoton stepNum 11(tk ;opticalphoton tid 69 pid 0 nm 430 mm ori[ -1.81e+04 -8e+05-6.60e+03] pos[ -7791.08e+03-2.04e+03] )
pre d/Geometry/AD/lvOIL#pvOAV terials/Acrylic Transportation GeomBoundary pos[ -649 899-2.49e+03] dir[ -0.237 0.329 -0.914] pol[ 0.971 0.0807 -0.223] ns 14.082 nm 430.000
post y/AD/lvOIL#pvBotReflector terials/Acrylic Transportation GeomBoundary pos[ -654 905-2.50e+03] dir[ -0.237 0.329 -0.914] pol[ 0.971 0.0807 -0.223] ns 14.184 nm 430.000
)
[ 4/ 11]
stepStatus [ GeomBoundary/ GeomBoundary]
flag [ ./ BOUNDARY_TRANSMIT]
bs pri/cur [ SAM/ FrT]
material [ Acrylic/ Air]
(rec state ;opticalphoton stepNum 11(tk ;opticalphoton tid 69 pid 0 nm 430 mm ori[ -1.81e+04 -8e+05-6.60e+03] pos[ -7791.08e+03-2.04e+03] )
pre y/AD/lvOIL#pvBotReflector terials/Acrylic Transportation GeomBoundary pos[ -654 905-2.50e+03] dir[ -0.237 0.329 -0.914] pol[ 0.971 0.0807 -0.223] ns 14.184 nm 430.000
post vBotReflector#pvBotRefGap d/Materials/Air Transportation GeomBoundary pos[ -656 909-2.51e+03] dir[ -0.357 0.494 -0.793] pol[ 0.932 0.123 -0.342] ns 14.241 nm 430.000
)
///
/// zero flags confirmed to be associated with same material (Acrylic->Acrylic) steps,
/// but the Acrylic are from different volumes OAV -> botReflector
///
[ 5/ 11]
stepStatus [ GeomBoundary/ GeomBoundary]
flag [ BOUNDARY_TRANSMIT/ SURFACE_SREFLECT]
bs pri/cur [ FrT/ SpR]
material [ Air/ ESR]
(rec state ;opticalphoton stepNum 11(tk ;opticalphoton tid 69 pid 0 nm 430 mm ori[ -1.81e+04 -8e+05-6.60e+03] pos[ -7791.08e+03-2.04e+03] )
pre vBotReflector#pvBotRefGap d/Materials/Air Transportation GeomBoundary pos[ -656 909-2.51e+03] dir[ -0.357 0.494 -0.793] pol[ 0.932 0.123 -0.342] ns 14.241 nm 430.000
post ails/lvBotRefGap#pvBotESR d/Materials/ESR Transportation GeomBoundary pos[ -656 909-2.51e+03] dir[ -0.357 0.494 0.793] pol[ -0.932 -0.123 -0.342] ns 14.241 nm 430.000
)
[ 6/ 11]
stepStatus [ GeomBoundary/ GeomBoundary]
flag [ SURFACE_SREFLECT/ NAN_ABORT]
bs pri/cur [ SpR/ STS]
material [ ESR/ Air]
(rec state ;opticalphoton stepNum 11(tk ;opticalphoton tid 69 pid 0 nm 430 mm ori[ -1.81e+04 -8e+05-6.60e+03] pos[ -7791.08e+03-2.04e+03] )
pre ails/lvBotRefGap#pvBotESR d/Materials/ESR Transportation GeomBoundary pos[ -656 909-2.51e+03] dir[ -0.357 0.494 0.793] pol[ -0.932 -0.123 -0.342] ns 14.241 nm 430.000
post vBotReflector#pvBotRefGap d/Materials/Air Transportation GeomBoundary pos[ -656 909-2.51e+03] dir[ -0.357 0.494 0.793] pol[ -0.932 -0.123 -0.342] ns 14.241 nm 430.000
)
[ 7/ 11]
stepStatus [ GeomBoundary/ GeomBoundary]
flag [ SURFACE_SREFLECT/ BOUNDARY_TRANSMIT]
bs pri/cur [ STS/ FrT]
material [ Air/ Acrylic]
(rec state ;opticalphoton stepNum 11(tk ;opticalphoton tid 69 pid 0 nm 430 mm ori[ -1.81e+04 -8e+05-6.60e+03] pos[ -7791.08e+03-2.04e+03] )
pre vBotReflector#pvBotRefGap d/Materials/Air Transportation GeomBoundary pos[ -656 909-2.51e+03] dir[ -0.357 0.494 0.793] pol[ -0.932 -0.123 -0.342] ns 14.241 nm 430.000
post y/AD/lvOIL#pvBotReflector terials/Acrylic Transportation GeomBoundary pos[ -656 909-2.51e+03] dir[ -0.237 0.329 0.914] pol[ -0.358 -0.904 0.232] ns 14.241 nm 430.000
)
[ 8/ 11]
stepStatus [ GeomBoundary/ GeomBoundary]
flag [ BOUNDARY_TRANSMIT/ .]
bs pri/cur [ FrT/ SAM]
material [ Acrylic/ Acrylic]
(rec state ;opticalphoton stepNum 11(tk ;opticalphoton tid 69 pid 0 nm 430 mm ori[ -1.81e+04 -8e+05-6.60e+03] pos[ -7791.08e+03-2.04e+03] )
pre y/AD/lvOIL#pvBotReflector terials/Acrylic Transportation GeomBoundary pos[ -656 909-2.51e+03] dir[ -0.237 0.329 0.914] pol[ -0.358 -0.904 0.232] ns 14.241 nm 430.000
post d/Geometry/AD/lvOIL#pvOAV terials/Acrylic Transportation GeomBoundary pos[ -659 913-2.50e+03] dir[ -0.237 0.329 0.914] pol[ -0.358 -0.904 0.232] ns 14.297 nm 430.000
)
[ 9/ 11]
stepStatus [ GeomBoundary/ GeomBoundary]
flag [ ./ BOUNDARY_TRANSMIT]
bs pri/cur [ SAM/ FrT]
material [ Acrylic/ LiquidScintillator]
(rec state ;opticalphoton stepNum 11(tk ;opticalphoton tid 69 pid 0 nm 430 mm ori[ -1.81e+04 -8e+05-6.60e+03] pos[ -7791.08e+03-2.04e+03] )
pre d/Geometry/AD/lvOIL#pvOAV terials/Acrylic Transportation GeomBoundary pos[ -659 913-2.50e+03] dir[ -0.237 0.329 0.914] pol[ -0.358 -0.904 0.232] ns 14.297 nm 430.000
post d/Geometry/AD/lvOAV#pvLSO uidScintillator Transportation GeomBoundary pos[ -664 919-2.49e+03] dir[ -0.238 0.33 0.913] pol[ -0.359 -0.904 0.233] ns 14.400 nm 430.000
)
[ 10/ 11]
stepStatus [ GeomBoundary/ GeomBoundary]
flag [ BOUNDARY_TRANSMIT/ BOUNDARY_TRANSMIT]
bs pri/cur [ FrT/ FrT]
material [ LiquidScintillator/ Acrylic]
(rec state ;opticalphoton stepNum 11(tk ;opticalphoton tid 69 pid 0 nm 430 mm ori[ -1.81e+04 -8e+05-6.60e+03] pos[ -7791.08e+03-2.04e+03] )
pre d/Geometry/AD/lvOAV#pvLSO uidScintillator Transportation GeomBoundary pos[ -664 919-2.49e+03] dir[ -0.238 0.33 0.913] pol[ -0.359 -0.904 0.233] ns 14.400 nm 430.000
post d/Geometry/AD/lvLSO#pvIAV terials/Acrylic Transportation GeomBoundary pos[ -7791.08e+03-2.04e+03] dir[ -0.237 0.329 0.914] pol[ -0.358 -0.904 0.232] ns 16.884 nm 430.000
)
(rec)FlagSequence TORCH BOUNDARY_TRANSMIT BOUNDARY_TRANSMIT BOUNDARY_TRANSMIT . BOUNDARY_TRANSMIT SURFACE_SREFLECT BOUNDARY_TRANSMIT . BOUNDARY_TRANSMIT . . . . . .
(rec)MaterialSequence GdDopedLS Acrylic LiquidScintillator Acrylic Acrylic Air Air Acrylic Acrylic LiquidScintillator - - - - - -
2016-10-01 17:53:24.035 INFO [1267418] [Rec::sequence@200] Rec::sequence nstep 11
Photons should all be starting from same place:
ipython -i $(which tokg4.py)
In [2]: a.rpost_(0)
Out[2]:
A()sliced
A([[ -18079.4443, -799699.4149, -6604.9499, 0.0977],
[ -18079.4443, -799699.4149, -6604.9499, 0.0977],
[ -18079.4443, -799699.4149, -6604.9499, 0.0977],
...,
[ -18079.4443, -799699.4149, -6604.9499, 0.0977],
[ -18079.4443, -799699.4149, -6604.9499, 0.0977],
[ -18079.4443, -799699.4149, -6604.9499, 0.0977]])
In [2]: a.gs
Out[2]:
A(torch,1,dayabay)-
A([[[ 0. , 0. , 0. , 0. ],
[ -18079.4531, -799699.4375, -6605. , 0.1 ],
[ 0. , 0. , 1. , 1. ],
[ 0. , 0. , 0. , 430. ],
[ 0. , 1. , 0. , 1. ],
[ 0. , 0. , 0. , 0. ]]], dtype=float32)
## huh: B photons not starting from where genstep points
## OR CTorchSource::configure _t 0.1 _radius 0 _pos -18079.4531,-799699.4375,-6605.0000
In [3]: b.rpost_(0)
Out[3]:
A()sliced
A([[ -24230.8603, -809820.8603, -0.0785, 0.0977],
[ -24230.8603, -809820.8603, -0.0785, 0.0977],
[ -24230.8603, -809820.8603, -0.0785, 0.0977],
...,
[ -24230.8603, -809820.8603, -0.0785, 0.0977],
[ -24230.8603, -809820.8603, -0.0785, 0.0977],
[ -24230.8603, -809820.8603, -0.0785, 0.0977]])
In [3]: b.gs
Out[3]:
A(torch,-1,dayabay)-
A([[[ 0. , 0. , 0. , 0. ],
[ -18079.4531, -799699.4375, -6605. , 0.1 ],
[ 0. , 0. , 1. , 1. ],
[ 0. , 0. , 0. , 430. ],
[ 0. , 1. , 0. , 1. ],
[ 0. , 0. , 0. , 0. ]]], dtype=float32)
## primaries collected from CTorchSource all at (0,0,0,0,0.1)
## where is frame setup for the default torch source done ?
In [4]: pr = np.load("/tmp/blyth/opticks/cg4/primary.npy")
In [5]: pr
Out[5]:
array([[[ 0. , 0. , 0. , 0.1],
[ 0. , 0. , 0. , 0. ],
[ 0. , 0. , 0. , 0. ],
[ 0. , 0. , 0. , 0. ]],
The gensteps have position that looks to be frame targetted:
21 const char* TorchStepNPY::DEFAULT_CONFIG =
22 "type=sphere_"
23 "frame=3153_"
24 "source=0,0,0_"
25 "target=0,0,1_"
26 "photons=100000_"
27 "material=GdDopedLS_"
28 "wavelength=430_"
29 "weight=1.0_"
30 "time=0.1_"
31 "zenithazimuth=0,1,0,1_"
32 "radius=0_" ;
33
34 // Aug 2016: change default torch wavelength from 380nm to 430nm
35 //
36 //
37 // NB time 0.f causes 1st step record rendering to be omitted, as zero is special
38 // NB the material string needs to be externally translated into a material line
Huh CTorchSource operating direct from TorchStepNPY, not the targetted NPY that it creates:
35
36 CTorchSource::CTorchSource(TorchStepNPY* torch, unsigned int verbosity)
37 :
38 CSource(verbosity),
39 m_torch(torch),
With bouncemax zero propagation is immediately terminated in both Opticks and G4, so can see initial photon position from photon buffer without the compression/decompression complications of the record buffer:
OKG4Test --save --compute --bouncemax 0
In [5]: a.ox[:,0] ## Opticks as expected
Out[5]:
A()sliced
A([[ -18079.453, -799699.438, -6605. , 0.1 ],
[ -18079.453, -799699.438, -6605. , 0.1 ],
[ -18079.453, -799699.438, -6605. , 0.1 ],
...,
[ -18079.453, -799699.438, -6605. , 0.1 ],
[ -18079.453, -799699.438, -6605. , 0.1 ],
[ -18079.453, -799699.438, -6605. , 0.1 ]], dtype=float32)
In [6]: b.ox[:,0] ## G4: real crazy position and time
Out[6]:
A()sliced
A([[ 0. , 0. , -2400000. , 8005.638],
[ 0. , 0. , -816713.875, 2724.364],
[ 0. , 0. , -1618713.875, 5399.548],
...,
[ 0. , 0. , -2062325.125, 6879.276],
[ 0. , 0. , -2400000. , 8005.638],
[ 0. , 0. , -1681468.25 , 5608.874]], dtype=float32)
After handling sphere positioning, gets a bit better:
In [2]: a.ox[:,0]
Out[2]:
A()sliced
A([[ -18079.4531, -799699.4375, -6605. , 0.1 ],
[ -18079.4531, -799699.4375, -6605. , 0.1 ],
[ -18079.4531, -799699.4375, -6605. , 0.1 ],
...,
[ -18079.4531, -799699.4375, -6605. , 0.1 ],
[ -18079.4531, -799699.4375, -6605. , 0.1 ],
[ -18079.4531, -799699.4375, -6605. , 0.1 ]], dtype=float32)
In [3]: b.ox[:,0]
Out[3]:
A()sliced
A([[ -18079.4531, -799699.4375, -8635. , 10.5231],
[ -18079.4531, -799699.4375, -6798.9727, 1.096 ],
[ -18079.4531, -799699.4375, -8635. , 10.5231],
...,
[ -18079.4531, -799699.4375, -8635. , 10.5231],
[ -18079.4531, -799699.4375, -8635. , 10.5231],
[ -18079.4531, -799699.4375, -8635. , 10.5231]], dtype=float32)
Direction should be random not all in -z dir:
In [8]: pr = np.load("cg4/primary.npy")
In [9]: pr
Out[9]:
array([[[ -18079.4531, -799699.4375, -6605. , 0.1 ],
[ 0. , 0. , -1. , 1. ],
[ 1. , 0. , 0. , 430. ],
[ 0. , 0. , 0. , 0. ]],
After setting iso get:
In [10]: pr = np.load("cg4/primary.npy")
In [11]: pr
Out[11]:
array([[[ -18079.4531, -799699.4375, -6605. , 0.1 ],
[ 0.6034, -0.673 , -0.4279, 1. ],
[ 0.7975, 0.5092, 0.3237, 430. ],
[ 0. , 0. , 0. , 0. ]],
[[ -18079.4531, -799699.4375, -6605. , 0.1 ],
[ 0.084 , -0.4561, 0.886 , 1. ],
[ 0.9965, 0.0384, -0.0747, 430. ],
[ 0. , 0. , 0. , 0. ]],
[[ -18079.4531, -799699.4375, -6605. , 0.1 ],
[ -0.2997, -0.7136, -0.6332, 1. ],
[ 0.954 , -0.2242, -0.1989, 430. ],
[ 0. , 0. , 0. , 0. ]],
Initially matseq stuck on G4:
A:seqmat_ana 1:dayabay
443231 0.441 44062 [6 ] Gd Ac LS Ac MO MO
33ff343231 0.090 9021 [10] Gd Ac LS Ac MO Ac Ai Ai Ac Ac
11 0.061 6051 [2 ] Gd Gd
aa33231 0.049 4859 [7 ] Gd Ac LS Ac Ac ES ES
3343343231 0.037 3688 [10] Gd Ac LS Ac MO Ac Ac MO Ac Ac
4432311 0.034 3351 [7 ] Gd Gd Ac LS Ac MO MO
dd43231 0.032 3164 [7 ] Gd Ac LS Ac MO Vm Vm
443343231 0.021 2116 [9 ] Gd Ac LS Ac MO Ac Ac MO MO
3323443231 0.015 1465 [10] Gd Ac LS Ac MO MO Ac LS Ac Ac
2231 0.013 1319 [4 ] Gd Ac LS LS
aa34231 0.011 1104 [7 ] Gd Ac LS MO Ac ES ES
111 0.011 1067 [3 ] Gd Gd Gd
4443231 0.009 940 [7 ] Gd Ac LS Ac MO MO MO
ff33424321 0.008 832 [10] Gd LS Ac MO LS MO Ac Ac Ai Ai
44323111 0.008 780 [8 ] Gd Gd Gd Ac LS Ac MO MO
dde43231 0.007 693 [8 ] Gd Ac LS Ac MO Py Vm Vm
334ff33231 0.007 656 [10] Gd Ac LS Ac Ac Ai Ai MO Ac Ac
4432231 0.006 559 [7 ] Gd Ac LS LS Ac MO MO
44ee43231 0.005 465 [9 ] Gd Ac LS Ac MO Py Py MO MO
44343231 0.004 423 [8 ] Gd Ac LS Ac MO Ac MO MO
100000 1.00
B:seqmat_ana -1:dayabay
11111111 0.429 42900 [8 ] Gd Gd Gd Gd Gd Gd Gd Gd
1111111111 0.399 39924 [10] Gd Gd Gd Gd Gd Gd Gd Gd Gd Gd
11 0.076 7584 [2 ] Gd Gd
111111111 0.037 3712 [9 ] Gd Gd Gd Gd Gd Gd Gd Gd Gd
1111 0.020 1952 [4 ] Gd Gd Gd Gd
111111 0.019 1900 [6 ] Gd Gd Gd Gd Gd Gd
1111111 0.012 1208 [7 ] Gd Gd Gd Gd Gd Gd Gd
111 0.005 455 [3 ] Gd Gd Gd
11111 0.004 365 [5 ] Gd Gd Gd Gd Gd
100000 1.00
A:seqhis_ana 1:dayabay
8ccccd 0.439 4386 [6 ] TO BT BT BT BT SA
ccaccccccd 0.089 893 [10] TO BT BT BT BT BT BT SR BT BT
4cccccd 0.060 603 [7 ] TO BT BT BT BT BT AB
4d 0.060 600 [2 ] TO AB
cccbcccccd 0.039 389 [10] TO BT BT BT BT BT BR BT BT BT
8cccccd 0.029 286 [7 ] TO BT BT BT BT BT SA
8cbcccccd 0.028 280 [9 ] TO BT BT BT BT BT BR BT SA
8cccc6d 0.021 206 [7 ] TO SC BT BT BT BT SA
8ccccccd 0.016 160 [8 ] TO BT BT BT BT BT BT SA
4ccd 0.014 137 [4 ] TO BT BT AB
cccc9ccccd 0.014 136 [10] TO BT BT BT BT DR BT BT BT BT
8cccc5d 0.012 118 [7 ] TO RE BT BT BT BT SA
cccacccccd 0.011 108 [10] TO BT BT BT BT BT SR BT BT BT
accccccccd 0.008 82 [10] TO BT BT BT BT BT BT BT BT SR
45d 0.008 76 [3 ] TO RE AB
7cccccd 0.006 61 [7 ] TO BT BT BT BT BT SD
8cccc55d 0.005 54 [8 ] TO RE RE BT BT BT BT SA
4ccccd 0.005 53 [6 ] TO BT BT BT BT AB
4cccccccd 0.005 49 [9 ] TO BT BT BT BT BT BT BT AB
cccccccccd 0.004 45 [10] TO BT BT BT BT BT BT BT BT BT
10000 1.00
B:seqhis_ana -1:dayabay
8ccccccd 0.424 4240 [8 ] TO BT BT BT BT BT BT SA <<< CG4 recording has more "BT BT" than Opticks does ???
ccaccccccd 0.081 813 [10] TO BT BT BT BT BT BT SR BT BT
4d 0.079 787 [2 ] TO AB
cccbcccccd 0.062 619 [10] TO BT BT BT BT BT BR BT BT BT
c0cac0cccd 0.062 616 [10] TO BT BT BT ?0? BT SR BT ?0? BT
cc9ccccccd 0.032 320 [10] TO BT BT BT BT BT BT DR BT BT
cccccccccd 0.030 296 [10] TO BT BT BT BT BT BT BT BT BT
ccbccccccd 0.023 227 [10] TO BT BT BT BT BT BT BR BT BT
4ccd 0.019 188 [4 ] TO BT BT AB
8cccccc6d 0.015 146 [9 ] TO SC BT BT BT BT BT BT SA
8ccccd 0.013 131 [6 ] TO BT BT BT BT SA
cac00cc0cd 0.013 131 [10] TO BT ?0? BT BT ?0? ?0? BT SR BT
abaccccccd 0.008 79 [10] TO BT BT BT BT BT BT SR BR SR
4ccccd 0.007 69 [6 ] TO BT BT BT BT AB
8cbcccccd 0.007 68 [9 ] TO BT BT BT BT BT BR BT SA
4ccccccd 0.006 60 [8 ] TO BT BT BT BT BT BT AB
ccc0b0cccd 0.006 59 [10] TO BT BT BT ?0? BR ?0? BT BT BT
4cccccd 0.005 46 [7 ] TO BT BT BT BT BT AB
bc9ccccccd 0.004 42 [10] TO BT BT BT BT BT BT DR BT BR
cabac0cccd 0.004 40 [10] TO BT BT BT ?0? BT SR BR SR BT
10000 1.00
A:seqmat_ana 1:dayabay
443231 0.444 4437 [6 ] Gd Ac LS Ac MO MO
33ff343231 0.089 894 [10] Gd Ac LS Ac MO Ac Ai Ai Ac Ac
11 0.060 600 [2 ] Gd Gd
aa33231 0.048 478 [7 ] Gd Ac LS Ac Ac ES ES
3343343231 0.037 373 [10] Gd Ac LS Ac MO Ac Ac MO Ac Ac
4432311 0.033 328 [7 ] Gd Gd Ac LS Ac MO MO
dd43231 0.029 295 [7 ] Gd Ac LS Ac MO Vm Vm
443343231 0.023 232 [9 ] Gd Ac LS Ac MO Ac Ac MO MO
3323443231 0.015 146 [10] Gd Ac LS Ac MO MO Ac LS Ac Ac
2231 0.013 135 [4 ] Gd Ac LS LS
aa34231 0.011 113 [7 ] Gd Ac LS MO Ac ES ES
4443231 0.011 105 [7 ] Gd Ac LS Ac MO MO MO
111 0.010 97 [3 ] Gd Gd Gd
ff33424321 0.008 81 [10] Gd LS Ac MO LS MO Ac Ac Ai Ai
334ff33231 0.008 78 [10] Gd Ac LS Ac Ac Ai Ai MO Ac Ac
44323111 0.007 70 [8 ] Gd Gd Gd Ac LS Ac MO MO
dde43231 0.007 67 [8 ] Gd Ac LS Ac MO Py Vm Vm
4432231 0.006 62 [7 ] Gd Ac LS LS Ac MO MO
44ee43231 0.005 54 [9 ] Gd Ac LS Ac MO Py Py MO MO
ff343231 0.005 47 [8 ] Gd Ac LS Ac MO Ac Ai Ai
10000 1.00
B:seqmat_ana -1:dayabay
44343231 0.425 4249 [8 ] Gd Ac LS Ac MO Ac MO MO
43ff343231 0.081 814 [10] Gd Ac LS Ac MO Ac Ai Ai Ac MO
11 0.079 787 [2 ] Gd Gd
233ff33231 0.061 607 [10] Gd Ac LS Ac Ac Ai Ai Ac Ac LS
2343343231 0.039 390 [10] Gd Ac LS Ac MO Ac Ac MO Ac LS
4344343231 0.030 298 [10] Gd Ac LS Ac MO Ac MO MO Ac MO
4343343231 0.022 216 [10] Gd Ac LS Ac MO Ac Ac MO Ac MO
ebd5e43231 0.021 215 [10] Gd Ac LS Ac MO Py Bk Vm OV Py
443231 0.020 200 [6 ] Gd Ac LS Ac MO MO
2231 0.019 188 [4 ] Gd Ac LS LS
443432311 0.015 146 [9 ] Gd Gd Ac LS Ac MO Ac MO MO
3ff3332331 0.013 131 [10] Gd Ac Ac LS Ac Ac Ac Ai Ai Ac
4e55e43231 0.013 130 [10] Gd Ac LS Ac MO Py Bk Bk Py MO
4ed5e43231 0.009 95 [10] Gd Ac LS Ac MO Py Bk Vm Py MO
ffff343231 0.008 81 [10] Gd Ac LS Ac MO Ac Ai Ai Ai Ai
443343231 0.007 70 [9 ] Gd Ac LS Ac MO Ac Ac MO MO
444343231 0.006 62 [9 ] Gd Ac LS Ac MO Ac MO MO MO
e5d5e43231 0.005 46 [10] Gd Ac LS Ac MO Py Bk Vm Bk Py
4323333231 0.004 42 [10] Gd Ac LS Ac Ac Ac Ac LS Ac MO
3ffff33231 0.004 40 [10] Gd Ac LS Ac Ac Ai Ai Ai Ai Ac
10000 1.00