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 ]

Diff markup

Differences between /examples/advanced/CaTS/scripts/hit.py (Version 11.3.0) and /examples/advanced/CaTS/scripts/hit.py (Version 10.7.p1)


  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()