from plot import qtplot
import json
import math as ma
import numpy as np

path='/cloud/Public/_data/neuropercolation/4lay/cons=27-knight_steps=1000100_diares/dim=09/batch=0/allpos_ei_samples=10000/dt=1/'
eps_sp=np.linspace(0.005,0.2,40)
eps_ep=eps_sp

extremes=10000
strength=100000

savepath=path+f'extremes={extremes}_bootstrength={strength}/'

with open(path+f"phasediff_means.txt", 'r', encoding='utf-8') as f:
    means=json.load(f)
with open(path+f"phasediff_stds.txt", 'r', encoding='utf-8') as f:
    stds=json.load(f)
        

hirim=[means[i]+stds[i] for i in range(len(means))]
lorim=[means[i]-stds[i] for i in range(len(means))]

   
qtplot(f'Reduction of phase deviation to noise level',[[0]+list(eps_sp[:40])]*3,[[0]+lorim,[0]+means,[0]+hirim],['mean-std','mean','mean+std'],colors=['r','g','r'],y_tag='phase difference',y_log=False)