diff --git a/HumedadConNeutrones/codigos python/NeutronTrack0.py b/HumedadConNeutrones/codigos python/NeutronTrack0.py index a8826b9f31cb35b0207cb6f29f26de308728f12b..f3ed2130ccae7f86d0133120eec3138e6f43891c 100644 --- a/HumedadConNeutrones/codigos python/NeutronTrack0.py +++ b/HumedadConNeutrones/codigos python/NeutronTrack0.py @@ -1,34 +1,42 @@ -#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)