Links

Content Skeleton

This Page

Nopstep failing to appear : viz debuggΒΆ

As usual when adding a vizualization : nothing appears and no error messages.

ggv-;ggv-g4gun
      # geant4 particle gun simulation within default DYB geometry, loaded from GDML

ggv-;ggv-g4gun --load
      # visualize the geant4 propagation, with GGeoView

Gun particles start from mid GdLS (volume 3153):

0.543      0.840      0.000 -18079.453
-0.840      0.543      0.000 -799699.438
0.000      0.000      1.000  -6605.000
0.000      0.000      0.000      1.000

At save they look to be in the correct region:

[2016-May-25 16:41:06.388981]:info: NumpyEvt::save (nopstep) (144,4,4)

(  0)  -18079.453  -799699.438   -6605.000       0.000
(  0)       0.000       0.000       1.000       1.000
(  0)       0.543      -0.840       0.000   10000.000
(  0)       0.000       0.000       0.000       0.000
(  1)  -18079.453  -799699.438   -6601.928       0.025
(  1)      -0.007       0.073       0.997       1.000
(  1)       0.543      -0.840       0.000    7611.374

With a 10 MeV muon the range of the nopstep is very large due to some neutrinos that exit the world.:

Rdr::upload mvn name rpos type  numbytes 9216 stride 64 offset 0 count 144 extent 2400000.000000
          m_low vec3  -735007.625 -799902.000 -2400000.000
         m_high vec3  560972.188   548652.250  2400000.000
   m_dimensions vec3  1295979.750 1348554.250  4800000.000
       m_center vec3  -87017.719 -125624.875      0.000
m_model_to_world mat4
2400000.000      0.000      0.000 -87017.719
     0.000 2400000.000      0.000 -125624.875
     0.000      0.000 2400000.000      0.000
     0.000      0.000      0.000      1.000

In [19]: (a[:,0,0].min(),a[:,0,1].min(),a[:,0,2].min(),a[:,0,3].min())
Out[19]: (-735007.62, -799902.0, -2400000.0, 0.0)

In [20]: (a[:,0,0].max(),a[:,0,1].max(),a[:,0,2].max(),a[:,0,3].max())
Out[20]: (560972.19, 548652.25, 2400000.0, 12151.727)

Switch to 10 MeV e+ (later 100 MeV), much more compact:

Rdr::upload mvn name rpos type  numbytes 13568 stride 64 offset 0 count 212 extent 440.272217
          m_low vec3  -18345.670 -799975.375  -7250.606
         m_high vec3  -18067.615 -799503.625  -6370.062
   m_dimensions vec3     278.055    471.750    880.544
       m_center vec3  -18206.643 -799739.500  -6810.334

Adding a simple “dbg” point shader succeeds to find the gross problems.