R_tunnel_init=28.0 # kOhm
#R_tunnel_init=1.0 # kOhm
TEC_init=5e-3 # K
#R,C,T,vol
#bounds = [(1.0,200.0),(1.0,200.0),(5.0,100.0),(0.2,200.0)]
#bounds = [(1.0,100.0),(100.0,920.0),(5.0,250.0)]
bounds = [(1.0,100.0),(230.0,240.0),(5.0,250.0)]
#bounds = None
filenames = ["data/test_meas_data.txt"]
fitters=[]
for i,filename in enumerate(filenames):
fitter=CBT_fitter(filename=filename,T_init=T_init,island_size_init=island_size_init,
R_tunnel_init=R_tunnel_init,TEC_init=TEC_init,
bounds=bounds,const_P=500e-18,
excitation=excitation,parallel_arrays=parallel_arrays,
junctions_in_series = N_junctions)
fitter.find_offset(5)
fitter.print_offset_curve()
fitter.fit_classic_curve()
fitter.plot_result_initial()
fitter.fit_full_curve()
fitter.print_elapset_time()
fitter.plot_nonlin_results()
fitters.append(fitter)
#fitter.plot_all_results()
cc = CBT_plot_data_pyx(fitters)
cc.plot_data(filename="res_8p5_50_100")
island_size_init=island_size_init,
R_tunnel_init=R_tunnel_init,
TEC_init=TEC_init,
bounds=bounds,
const_P=500e-18,
excitation=excitation,
parallel_arrays=parallel_arrays,
junctions_in_series=N_junctions)
fitter.find_offset(
5) # finds offset with 5 points from minimum each direction
#fitter.print_offset_curve()
fitter.fit_classic_curve(
) # fit analytic curve to get starting point for nonlin fit
#fitter.plot_result_initial()
fitter.fit_full_curve() # nonlin fitting
fitter.print_elapset_time()
#fitter.plot_nonlin_results()
fitters.append(fitter) # save results
#fitter.plot_all_results()
# fitting all together
multi_fitter = CBT_multi_fitter(
fitters,
result_text_filename=results_txt_file,
interpolation_table_file=interpolation_filename)
# plot results
cc = CBT_plot_data_pyx(fitters)
cc.plot_multi_data(filename=result_pdf_1)
cc.plot_data(filename=result_pdf_2)
print "it took", time.time() - start, "seconds."
interpolation_filename="41A4_interpolation_table.txt"
# separate fitting of each curve
fitters=[] # list of classes CBT_fitter that also store results
for i,filename in enumerate(filenames):
fitter=CBT_fitter(filename=filename,T_init=T_init,island_size_init=island_size_init,
R_tunnel_init=R_tunnel_init,TEC_init=TEC_init,
bounds=bounds,const_P=500e-18,
excitation=excitation,parallel_arrays=parallel_arrays,
junctions_in_series = N_junctions)
fitter.find_offset(5) # finds offset with 5 points from minimum each direction
#fitter.print_offset_curve()
fitter.fit_classic_curve() # fit analytic curve to get starting point for nonlin fit
#fitter.plot_result_initial()
fitter.fit_full_curve() # nonlin fitting
fitter.print_elapset_time()
#fitter.plot_nonlin_results()
fitters.append(fitter) # save results
#fitter.plot_all_results()
# fitting all together
multi_fitter = CBT_multi_fitter(fitters,
result_text_filename=results_txt_file,
interpolation_table_file=interpolation_filename)
# plot results
cc = CBT_plot_data_pyx(fitters)
cc.plot_multi_data(filename=result_pdf_1)
cc.plot_data(filename=result_pdf_2)
print "it took", time.time() - start, "seconds."
