Update NeutronTrack0.py

HumedadConNeutrones/codigos python/NeutronTrack0.py

-#Este programa grafica la ubicacion previa de los neutrones en el suelo
-#antes de ser detectados por un detector fisico simulado en URANOS
-import numpy as np
+#Este programa genera una grafica de el viaje, que hizo el neutron 0, desde que viene de la
+# fuente hasta que es abosorbido. 
 import matplotlib.pyplot as plt
-import pandas as pd
 
-#Dectector hits.
-#Los datos corresponden a 10% de humedad en el suelo con porosidad del 50%
+dlayer= pd.read_table('/work/detectorLayer300000.dat')
+dftrack = dlayer[dlayer['Neutron_Number'] == 0]
 
-#Leemos los datos y los almacenamos en un dataframe
-dlayer= pd.read_table('/work/SALIDASdetectorNeutronHitData.dat')
-datareal=dlayer[(dlayer['maximum_Depth_[m]']>0)]#Filtramos para graficar la profundidad que alcanzan las particulas en el suelo humedo
-data1=datareal.assign(hola=(datareal['maximum_Depth_[m]']*(-1)))# Es necesario invertir el signo (ya que las distancias por debajo del suelo tienen numero postivo, y es mas intuitivo trabajarlas como distancias negativas)
+#Detector Layer
+x = [0,200]
+y = [2,2]
+y2 = [2.5,2.5]
 
-#Graficamos vista aerea de el suelo
-# es decir, el plano x-y con z = 0
-datareal.plot(x='x_at_interface_[m]', y='y_at_Interface_[m]', kind="scatter", s = 4, fontsize=15, alpha=0.3)
-plt.xlabel('x(m)', size=18)
-plt.ylabel('y(m)', size=18)
-plt.scatter(0,0, c='red')
-plt.xlim(-300,300)
-plt.ylim(-300,300)
+# Ground
+y1 = [0,0]
+
+#"Tracking del primer neutron"
+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)
+plt.plot(x,y,'g',label='Detector layer')
+plt.plot(x,y1,'brown', label='Ground')
+plt.plot(x,y2,'g')
+plt.xlim(0,200)
+plt.ylim(-2,55)
+plt.xlabel('x_[m]', size=16)
+plt.ylabel('Depth_[m]', size = 16)
+plt.tick_params(labelsize=16)
+
+plt.text(169,35,'91Mev',fontsize=17)
 plt.grid()
-plt.show()
-###########
-#Graficamos vista  transversal del suelo
-# Es decir el plano x-z con y=0
-#########
+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')
 plt.xlabel('x(m)', size=18)
 plt.ylabel('z(m)', size=18)