@@ -1050,28 +1050,15 @@ def qc3_bookend(readname, savename, time_strs, args):
10501050 save_fname = f"{ well } _{ savename } _before0.pdf"
10511051
10521052 # Plot subtraction
1053- # get_leak_corrected(first_before_current,
1054- # voltage, times,
1055- # *ramp_bounds,
1056- # save_fname=save_fname,
1057- # output_dir=output_directory)
1058-
1059- before_traces_first [well ] = get_leak_corrected (first_before_current ,
1060- voltage , times ,
1061- * ramp_bounds ,
1062- save_fname = save_fname ,
1063- output_dir = output_directory )
1064-
1065- before_traces_last [well ] = get_leak_corrected (last_before_current ,
1066- voltage , times ,
1067- * ramp_bounds )
1068-
1069- after_traces_first [well ] = get_leak_corrected (first_after_current ,
1070- voltage , times ,
1071- * ramp_bounds )
1072- after_traces_last [well ] = get_leak_corrected (last_after_current ,
1073- voltage , times ,
1074- * ramp_bounds )
1053+ before_traces_first [well ] = get_leak_corrected (
1054+ first_before_current , voltage , times , * ramp_bounds ,
1055+ save_fname = save_fname ,output_dir = output_directory )
1056+ before_traces_last [well ] = get_leak_corrected (
1057+ last_before_current , voltage , times , * ramp_bounds )
1058+ after_traces_first [well ] = get_leak_corrected (
1059+ first_after_current , voltage , times , * ramp_bounds )
1060+ after_traces_last [well ] = get_leak_corrected (
1061+ last_after_current , voltage , times , * ramp_bounds )
10751062
10761063 # Store subtracted traces
10771064 first_processed [well ] = before_traces_first [well ] - after_traces_first [well ]
@@ -1094,21 +1081,13 @@ def qc3_bookend(readname, savename, time_strs, args):
10941081 ax = fig .subplots ()
10951082 for well in args .wells :
10961083 trace1 = hergqc .filter_capacitive_spikes (
1097- first_processed [well ], times , voltage_steps
1098- ).flatten ()
1099-
1084+ first_processed [well ], times , voltage_steps ).flatten ()
11001085 trace2 = hergqc .filter_capacitive_spikes (
1101- last_processed [well ], times , voltage_steps
1102- ).flatten ()
1103-
1086+ last_processed [well ], times , voltage_steps ).flatten ()
11041087 passed = hergqc .qc3 (trace1 , trace2 )[0 ]
1105-
11061088 res_dict [well ] = passed
1107-
1108- save_fname = os .path .join (args .output_dir ,
1109- 'debug' ,
1110- f"debug_{ well } _{ savename } " ,
1111- 'qc3_bookend' )
1089+ save_fname = os .path .join (
1090+ args .output_dir , 'debug' , f"debug_{ well } _{ savename } " , 'qc3_bookend' )
11121091
11131092 ax .plot (times , trace1 )
11141093 ax .plot (times , trace2 )
0 commit comments