Update NeutronTrack0.py

HumedadConNeutrones/codigos python/NeutronTrack0.py

@@ -2,27 +2,35 @@
 # fuente hasta que es abosorbido. 
 import matplotlib.pyplot as plt
 
+#####################################
+#Leemos los datos para extraer los puntos de interaccion y  previo a la interacion con el suelo y capa de deteccion.
 dlayer= pd.read_table('/work/detectorLayer300000.dat')
 dftrack = dlayer[dlayer['Neutron_Number'] == 0]
+##########################################
 
-#Detector Layer
+#Dimensiones de la capa de decteción.
 x = [0,200]
 y = [2,2]
 y2 = [2.5,2.5]
 
-# Ground
+# Suelo. (Interfaz suelo-aire)
 y1 = [0,0]
 
-#"Tracking del primer neutron"
+#"Tracking del primer neutron" Los puntos extraidos leyendo el dataframe 'dftrack'
 x1 = [178.324,186.75,157.713,30.1466,55.3114,61.3679,127.643,190.402,197.086,196.943,187.272,187.315]
 z1 = [32.9121,2.5,53.5416,2.5,-0.0295019,2.0,24.2086,2.5,-0.0517591,2.0,-0.022933,2.0]
 
-plt.plot(x1,z1,'black',label='Neutron Journey')
-
-plt.scatter(x1,z1,color='red',s=16)
+#graficamos el escenario (Capa de decteccion y suelo)
 plt.plot(x,y,'g',label='Detector layer')
-plt.plot(x,y1,'brown', label='Ground')
 plt.plot(x,y2,'g')
+plt.plot(x,y1,'brown', label='Ground')
+##########################################
+#graficamos las posiciones del neutron
+plt.plot(x1,z1,'black',label='Neutron Journey')
+plt.scatter(x1,z1,color='red',s=16)
+##########################################
+#######################################3##
+#
 plt.xlim(0,200)
 plt.ylim(-2,55)
 plt.xlabel('x_[m]', size=16)
@@ -32,9 +40,6 @@ plt.tick_params(labelsize=16)
 plt.text(169,35,'91Mev',fontsize=17)
 plt.grid()
 plt.legend( fontsize='x-large')
-data1.plot(y='hola', x='x_at_interface_[m]', kind="scatter", s = 3, fontsize=15, alpha=0.3)
-plt.xlim(-300,300)
-plt.ylim(-1,0)
 plt.show()
 
 plt.scatter(0,0, c='red')