def get_SCAT(self):
if (
len(set(self.y_Arai_segment)) == 1
): # prevents divide by zero, i.e. if all y values in segment are the same [1,1,1]
self.pars['SCAT'] = 0 #float('nan')
self.pars['fail_arai_beta_box_scatter'] = 0 # float('nan')
self.pars["fail_ptrm_beta_box_scatter"] = 0 #float('nan')
self.pars["fail_tail_beta_box_scatter"] = 0 #float('nan')
self.pars['scat_bounding_line_high'] = 0 # float('nan')
self.pars['scat_bounding_line_low'] = 0 #float('nan')
return 0
slope = self.pars['specimen_b']
x_mean, y_mean = self.x_Arai_mean, self.y_Arai_mean
x_Arai_segment, y_Arai_segment = self.x_Arai_segment, self.y_Arai_segment
if 'specimen_b_beta' in self.acceptance_criteria:
box = lib_arai.get_SCAT_box(
slope,
x_mean,
y_mean,
beta_threshold=self.acceptance_criteria['specimen_b_beta']
['value'])
else:
box = lib_arai.get_SCAT_box(slope, x_mean, y_mean)
low_bound, high_bound, x_max, y_max, low_line, high_line = box[0], box[
1], box[2], box[3], box[4], box[5]
# getting SCAT points
x_Arai_segment, y_Arai_segment = self.x_Arai_segment, self.y_Arai_segment
tmin, tmax = self.tmin, self.tmax
ptrm_checks_temps, ptrm_checks_starting_temps = self.ptrm_checks_temperatures, self.ptrm_checks_starting_temperatures
x_ptrm_check, y_ptrm_check = self.x_ptrm_check, self.y_ptrm_check
tail_checks_temps, tail_checks_starting_temps = self.tail_checks_temperatures, self.tail_checks_starting_temperatures
x_tail_check, y_tail_check = self.x_tail_check, self.y_tail_check
points, fancy_points = lib_arai.get_SCAT_points(
x_Arai_segment, y_Arai_segment, tmin, tmax, ptrm_checks_temps,
ptrm_checks_starting_temps, x_ptrm_check, y_ptrm_check,
tail_checks_temps, tail_checks_starting_temps, x_tail_check,
y_tail_check)
# checking each point
SCAT = lib_arai.get_SCAT(points, low_bound, high_bound, x_max, y_max)
fancy_SCAT, SCATs = lib_arai.fancy_SCAT(fancy_points, low_bound,
high_bound, x_max, y_max)
#'SCAT_arai': False, 'SCAT_tail': True, 'SCAT_ptrm': True})
if fancy_SCAT.upper() == 'PASS':
self.pars['SCAT'] = 't'
else:
self.pars['SCAT'] = 'f'
#self.pars['SCAT'] = fancy_SCAT
# set self.pars['SCAT'] to 'g' or 'b' here?
self.pars['fail_arai_beta_box_scatter'] = SCATs['SCAT_arai']
self.pars["fail_ptrm_beta_box_scatter"] = SCATs['SCAT_ptrm']
self.pars["fail_tail_beta_box_scatter"] = SCATs['SCAT_tail']
self.pars['scat_bounding_line_high'] = high_line # [y_int, slope]
self.pars['scat_bounding_line_low'] = low_line # [y_int, slope]
return fancy_SCAT
def get_SCAT(self):
slope = self.pars['specimen_b']
x_mean, y_mean = self.x_Arai_mean, self.y_Arai_mean
x_Arai_segment, y_Arai_segment = self.x_Arai_segment, self.y_Arai_segment
box = lib_arai.get_SCAT_box(slope, x_mean, y_mean)
low_bound, high_bound, x_max, y_max, low_line, high_line = box[0], box[1], box[2], box[3], box[4], box[5]
# getting SCAT points
x_Arai_segment, y_Arai_segment = self.x_Arai_segment, self.y_Arai_segment
tmin, tmax = self.tmin, self.tmax
ptrm_checks_temps, ptrm_checks_starting_temps = self.ptrm_checks_temperatures, self.ptrm_checks_starting_temperatures
x_ptrm_check, y_ptrm_check = self.x_ptrm_check, self.y_ptrm_check
tail_checks_temps, tail_checks_starting_temps = self.tail_checks_temperatures, self.tail_checks_starting_temperatures
x_tail_check, y_tail_check = self.x_tail_check, self.y_tail_check
points, fancy_points = lib_arai.get_SCAT_points(x_Arai_segment, y_Arai_segment, tmin, tmax,
ptrm_checks_temps, ptrm_checks_starting_temps,
x_ptrm_check, y_ptrm_check, tail_checks_temps,
tail_checks_starting_temps, x_tail_check, y_tail_check)
# checking each point
SCAT = lib_arai.get_SCAT(points, low_bound, high_bound, x_max, y_max)
fancy_SCAT, SCATs = lib_arai.fancy_SCAT(fancy_points, low_bound, high_bound, x_max, y_max)
#'SCAT_arai': False, 'SCAT_tail': True, 'SCAT_ptrm': True})
self.pars['SCAT'] = fancy_SCAT
self.pars['fail_arai_beta_box_scatter'] = SCATs['SCAT_arai']
self.pars["fail_ptrm_beta_box_scatter"] = SCATs['SCAT_ptrm']
self.pars["fail_tail_beta_box_scatter"] = SCATs['SCAT_tail']
self.pars['scat_bounding_line_high'] = high_line # [y_int, slope]
self.pars['scat_bounding_line_low'] = low_line # [y_int, slope]
return fancy_SCAT
def get_SCAT(self):
if (len(set(self.y_Arai_segment)) == 1): # prevents divide by zero, i.e. if all y values in segment are the same [1,1,1]
self.pars['SCAT'] = 0 #float('nan')
self.pars['fail_arai_beta_box_scatter'] = 0# float('nan')
self.pars["fail_ptrm_beta_box_scatter"] = 0#float('nan')
self.pars["fail_tail_beta_box_scatter"] = 0#float('nan')
self.pars['scat_bounding_line_high'] = 0# float('nan')
self.pars['scat_bounding_line_low'] = 0#float('nan')
return 0
slope = self.pars['specimen_b']
x_mean, y_mean = self.x_Arai_mean, self.y_Arai_mean
x_Arai_segment, y_Arai_segment = self.x_Arai_segment, self.y_Arai_segment
if 'specimen_b_beta' in self.acceptance_criteria:
box = lib_arai.get_SCAT_box(slope, x_mean, y_mean, beta_threshold=self.acceptance_criteria['specimen_b_beta']['value'])
else:
box = lib_arai.get_SCAT_box(slope, x_mean, y_mean)
low_bound, high_bound, x_max, y_max, low_line, high_line = box[0], box[1], box[2], box[3], box[4], box[5]
# getting SCAT points
x_Arai_segment, y_Arai_segment = self.x_Arai_segment, self.y_Arai_segment
tmin, tmax = self.tmin, self.tmax
ptrm_checks_temps, ptrm_checks_starting_temps = self.ptrm_checks_temperatures, self.ptrm_checks_starting_temperatures
x_ptrm_check, y_ptrm_check = self.x_ptrm_check, self.y_ptrm_check
tail_checks_temps, tail_checks_starting_temps = self.tail_checks_temperatures, self.tail_checks_starting_temperatures
x_tail_check, y_tail_check = self.x_tail_check, self.y_tail_check
points, fancy_points = lib_arai.get_SCAT_points(x_Arai_segment, y_Arai_segment, tmin, tmax,
ptrm_checks_temps, ptrm_checks_starting_temps,
x_ptrm_check, y_ptrm_check, tail_checks_temps,
tail_checks_starting_temps, x_tail_check, y_tail_check)
# checking each point
SCAT = lib_arai.get_SCAT(points, low_bound, high_bound, x_max, y_max)
fancy_SCAT, SCATs = lib_arai.fancy_SCAT(fancy_points, low_bound, high_bound, x_max, y_max)
#'SCAT_arai': False, 'SCAT_tail': True, 'SCAT_ptrm': True})
if fancy_SCAT.upper() == 'PASS':
self.pars['SCAT'] = 't'
else:
self.pars['SCAT'] = 'f'
#self.pars['SCAT'] = fancy_SCAT
# set self.pars['SCAT'] to 'g' or 'b' here?
self.pars['fail_arai_beta_box_scatter'] = SCATs['SCAT_arai']
self.pars["fail_ptrm_beta_box_scatter"] = SCATs['SCAT_ptrm']
self.pars["fail_tail_beta_box_scatter"] = SCATs['SCAT_tail']
self.pars['scat_bounding_line_high'] = high_line # [y_int, slope]
self.pars['scat_bounding_line_low'] = low_line # [y_int, slope]
return fancy_SCAT
def test_SCAT_box(self):
# def get_SCAT_box(slope, slope_err, x_mean, y_mean, beta_threshold = .1):
ref_y_max = 6.5
ref_x_max = 7.
ref_low_bound_result = 1.8846153846153846
ref_high_bound_result = 4.642857142857142
result = lib_arai.get_SCAT_box(self.b, self.x_mean, self.y_mean, self.beta_threshold)
self.assertAlmostEqual(ref_low_bound_result, result[0](2.))
self.assertAlmostEqual(ref_high_bound_result, result[1](2.))
self.assertAlmostEqual(ref_x_max, result[2])
self.assertAlmostEqual(ref_y_max, result[3])
def test_in_SCAT_box(self):
low_bound, high_bound, x_max, y_max, low_line, high_line = lib_arai.get_SCAT_box(self.b,
self.x_mean, self.y_mean, self.beta_threshold)
good = [(2., 1.9), (1., 5.), (5., 1.8)]
bad = [(1., 2.), (3., 4), (7.1, .2)]
# def in_SCAT_box(x, y, low_bound, high_bound, x_max, y_max):
for xy in good:
result = lib_arai.in_SCAT_box(xy[0], xy[1], low_bound, high_bound, x_max, y_max)
self.assertTrue(result)
for xy in bad:
result = lib_arai.in_SCAT_box(xy[0], xy[1], low_bound, high_bound, x_max, y_max)
self.assertFalse(result)
