# -*- coding: utf-8 -*-
"""
Created on Sun Sep 16 12:17:38 2018

@author: 
"""

import numpy as np
import matplotlib.pyplot as plt
import os


Path = "carpeta con las mediciones"     
filelist = os.listdir(Path)
#levantamos los dt y tmedio de las distintas mediciones

dt=[]
tmedio=[]
for i in filelist:
    if i.endswith(".dat"):  
        #abrimos los archivos    
        archivo=np.loadtxt(i, dtype=float)
        tmedio.append(archivo[:,0])
        dt.append(archivo[:,1])
        
#calculamos la media y el desvio por franja de todas las mediciones
tmedio_mean = np.mean(tmedio,axis=0)
tmedio_std = np.std(tmedio,axis=0)
        
dt_mean = np.mean(dt,axis=0)
dt_std = np.std(dt,axis=0)
dt_error = dt_std/(len(tmedio))**(1/2)

#graficamos el dt para cada franja con su error
plt.figure()
plt.errorbar(tmedio_mean,dt_mean, xerr=None,yerr=dt_error,fmt="o")
plt.grid('on');
plt.xlabel('tmedio');
plt.ylabel('dt [1/s]');
plt.show()

#calculamos la velocidad para cada franja con su error
d = 5
vel = d/dt_mean
dt_errel= dt_std/dt_mean
vel_error = dt_errel*vel 

#graficamos la velocidad para cada franja con su error
plt.figure()
plt.errorbar(tmedio_mean,vel, xerr=None,yerr=vel_error,fmt="o")
plt.grid('on');
plt.xlabel('tmedio x franja');
plt.ylabel('Vel media [cm/s]');
plt.show()