#m=manager.ProductManager.from_config('/home/zahra/PIPELINE/ex_3by3_kl/prod_params_custom.yaml')
m = manager.ProductManager.from_config(
'/home/zahra/PIPELINE/example_7by7_lmax_mmax_250/prod_params_custom.yaml')
t = m.telescope
x = t.feedpositions[:,
0] #these are x and y positions not x and y polarizations
y = t.feedpositions[:, 1]
#klobj=m.kltransforms['kl']
time_channel = 400
Ndish = 49
correlation_arr, sum_counts, corr_counts = Bls_counts(m)
arr_pt01 = Visibilities_grid(m, ts_2, rg_2, ag_2, gn_2_50_1, time_channel, 50,
1)
#arr_pt1=Visibilities_grid(m,ts_2,rg_2,ag_2,gn_2_500_1,time_channel,Tsys=500,ndays=1)
#arr_pt5=Visibilities_grid(m,ts_2,rg_2,ag_2,gn_2_750_1,time_channel,Tsys=750,ndays=1)
#arr_10=Visibilities_grid(m,ts_2,rg_2,ag_2,gn_2_5e4_1,time_channel,Tsys=5e4,ndays=1)
arr = arr_pt01
#vis=arr[2][:,time_channel] #2 is gnoise, 0 is manual
sigma = arr[-1] #sigma_gn
sigma = 1
vis = arr[1][:, time_channel]
print(arr[1][:, time_channel].min(), arr[1][:, time_channel].max(),
'meas vis no noise')
print(vis.min(), vis.max(), 'meas vis noise')
'/home/zahra/PIPELINE/example_7by7_lmax_mmax_250/draco_synthesis_pt1_gainfluc/rand_gains/sig_pt12.h5',
'r')
gn_2 = h5py.File(
'/home/zahra/PIPELINE/example_7by7_lmax_mmax_250/draco_synthesis_pt1_gainfluc/rand_gains/ts_final/app_gain_noise_2.h5',
'r')
m = manager.ProductManager.from_config(
'/home/zahra/PIPELINE/example_7by7_lmax_mmax_250/prod_params_custom.yaml')
#klobj=m.kltransforms['kl']
time_channel = 400
Ndish = 49
correlation_arr, sum_counts, corr_counts = Bls_counts(m)
arr = Visibilities_grid(m, ts_2, rg_2, ag_2, gn_2, time_channel)
vis = arr[0][:, time_channel] #2 is gnoise, 0 is manual
print(arr[1][:, time_channel].min(), arr[1][:, time_channel].max(),
'meas vis no noise')
print(vis.min(), vis.max(), 'meas vis noise')
sigma = arr[-2] #-1 is sigma_gn, -2 is sigma for manual
#sigma=1.e-6
data = np.append(vis.real, vis.imag)
#data=np.ones(72)
data = np.zeros(2 * vis.size)
data[0::2] = vis.real
data[1::2] = vis.imag
'''
sum_full_real=np.zeros((102,102))
for i in range(201):