diff --git a/HumedadConNeutrones/codigos python/DetectorHits.py b/HumedadConNeutrones/codigos python/DetectorHits.py index f5ef192420ab79efc5ae0a26df382b1e560fda9b..6fc99f35b7c8ec49f588b9e7373fae3e9aea83ca 100644 --- a/HumedadConNeutrones/codigos python/DetectorHits.py +++ b/HumedadConNeutrones/codigos python/DetectorHits.py @@ -1,9 +1,5 @@ -# -*- coding: utf-8 -*- -""" -Spyder Editor - -This is a temporary script file. -""" +#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 @@ -11,12 +7,13 @@ 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)] -data1=datareal.assign(hola=(datareal['maximum_Depth_[m]']*(-1))) +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) @@ -27,6 +24,7 @@ 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) @@ -36,7 +34,3 @@ plt.xlabel('x(m)', size=18) plt.ylabel('z(m)', size=18) plt.grid() plt.show() - - - -