

In [ ]:

    
import sys
sys.path.append('../..')

import pyotc



In [ ]:

    
basename = 'B01_'
directory = '../exampleData/height_calibration_single_psds/'

hc = pyotc.HeightCalibration()
hc.load_hc_data(basename, directory=directory)



In [ ]:

    
hc.plot_pc_results()

exclude data points



In [ ]:

    
hc.exclude_heights_outside(0.8, 4, reset=True)

determine focal shift



In [ ]:

    
hc.determine_focal_shift(idx_slice=slice(16,-4))

set the focal shift



In [ ]:

    
hc.focal_shift = 0.68

fit relative drag



In [ ]:

    
hc.fit_rel_drag(method='radius', verbose=True, h0=-0.5, corr=1.2)

Plot relative drag fit result



In [ ]:

    
hc.plot_rel_drag_fit_result()

fit height data globally



In [ ]:

    
hc.fit_height_data(fit_drag_first=False, method='radius', fit_dissens_osci=True, plot_fit=True, fit_focal_shift=True)



In [ ]:

    
hc.write_results_to_file()



In [ ]: