diff --git a/HumedadConNeutrones/codigos python/NeutronTrack0.py b/HumedadConNeutrones/codigos python/NeutronTrack0.py index f3ed2130ccae7f86d0133120eec3138e6f43891c..ab75e60105adc009d1b5b63d2f43f9dfd678e568 100644 --- a/HumedadConNeutrones/codigos python/NeutronTrack0.py +++ b/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')