diff --git a/script/SINSB01_phase_scan.py b/script/SINSB01_phase_scan.py deleted file mode 100755 index f5c3dd1..0000000 --- a/script/SINSB01_phase_scan.py +++ /dev/null @@ -1 +0,0 @@ -lscan(SINSB01_phase, BC1_energy, 0.0, 360.0, 10.0, 2.0) \ No newline at end of file diff --git a/script/test_didier.py b/script/test_didier.py deleted file mode 100755 index 4394e03..0000000 --- a/script/test_didier.py +++ /dev/null @@ -1,38 +0,0 @@ -################################################################################################### -# Multiple Gaussians peak search with mathutils.py -################################################################################################### - - -from mathutils import estimate_peak_indexes, fit_gaussians, create_fit_point_list - -start = 0 -end = 50 -step_size = 0.2 - -result= lscan(ao1,ai1,start,end,[step_size,]) - -readable = result.getReadable(0) -positions = result.getPositions(0) - -threshold = (min(readable) + max(readable))/2 -min_peak_distance = 5.0 - -peaks = estimate_peak_indexes(readable, positions, threshold, min_peak_distance) -print "Peak indexes: " + str(peaks) -print "Peak x: " + str(map(lambda x:positions[x], peaks)) -print "Peak y: " + str(map(lambda x:readable[x], peaks)) - - - -gaussians = fit_gaussians(readable, positions, peaks) - - -plots = plot([readable],["sin"],[positions], title="Data" ) -for i in range(len(peaks)): - peak = peaks[i] - (norm, mean, sigma) = gaussians[i] - if abs(mean - positions[peak]) < min_peak_distance: - print "Peak -> " + str(mean) - plots[0].addMarker(mean, None, "N="+str(round(norm,2)), None) - else: - print "Invalid gaussian fit: " + str(mean) diff --git a/script/test_scan_with_fit.py b/script/test_scan_with_fit.py deleted file mode 100755 index c11684d..0000000 --- a/script/test_scan_with_fit.py +++ /dev/null @@ -1,7 +0,0 @@ - - - -res = lscan(SINSB01_phase, BC1_energy, -90, 90, 21, latency = 1.0, relative = True) -y = res.getReadable(0) -x = res.getPositions(0) -fit(y, x) \ No newline at end of file