Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/CaTS/scripts/hit.py

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

  1 #!/usr/bin/env python
  2 # Python program explaining  
  3 # load() function  
  4   
  5 import numpy as geek 
  6 import matplotlib.pyplot as plt 
  7   
  8 b = geek.load('hits.npy') 
  9   
 10 print("size of b is:",b.size)
 11 print("shape of b is:",b.shape)
 12 print("ndim of b is:",b.ndim)
 13 print(b)
 14 print(b[0,0,0])
 15 print(b[0,3,0])
 16 print("b is printed from hits.npy")
 17 x = b[:,0,0]
 18 y = b[:,0,1]
 19 z = b[:,0,2]
 20 lamb = b[:,2,3]
 21 print(x.shape)
 22 print(y.shape)
 23 print(lamb.shape)
 24 
 25 #axs = plt.subplots(2, 2, figsize=(5, 5))
 26 plt.subplot(131)
 27 plt.ylabel('#Photons')
 28 plt.xlabel('x-position/mm')
 29 plt.hist(x,50,(-200,200.))
 30 plt.subplot(132)
 31 plt.xlabel('x-position/mm')
 32 plt.ylabel('y-position/mm')
 33 plt.scatter(x, y)
 34 plt.subplot(133)
 35 plt.xlabel('lambda/nm')
 36 plt.ylabel('#Photons')
 37 plt.hist(lamb,50,(50,800.))
 38 #plt.hist(z)
 39 plt.show()