diff --git a/HumedadConNeutrones/codigos python/NeutronTrack0.py b/HumedadConNeutrones/codigos python/NeutronTrack0.py
index 3359655f0f9c15cedccc11611926b822f5ec5e0d..a8826b9f31cb35b0207cb6f29f26de308728f12b 100644
--- a/HumedadConNeutrones/codigos python/NeutronTrack0.py	
+++ b/HumedadConNeutrones/codigos python/NeutronTrack0.py	
@@ -1,38 +1,38 @@
-# -*- coding: utf-8 -*-
-"""
-Created on Sun Dec 11 20:12:15 2022
-
-@author: Usuario
-"""
-
+#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
 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%
+
+#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)
+
+#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)
+plt.grid()
+plt.show()
+###########
+#Graficamos vista  transversal del suelo
+# Es decir el plano x-z con y=0
+#########
+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.scatter(0,0, c='red')
+plt.xlabel('x(m)', size=18)
+plt.ylabel('z(m)', size=18)
+plt.grid()
+plt.show()
 
-# Track del neutron 0
-
-#Detector Layer
-x = [0,200]
-y = [2,2]
-y2 = [2.5,2.5]
-
-# 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.legend( fontsize='x-large')
\ No newline at end of file