code: def get_SCAT_points(x_Arai_segment, y_Arai_segment, tmin, tmax, ptrm_checks_temperatures,
                      ptrm_checks_starting_temperatures, x_ptrm_check, y_ptrm_check,
                      tail_checks_temperatures, tail_checks_starting_temperatures,
                      x_tail_check, y_tail_check):
      
      points = []
      points_arai = []
      points_ptrm = []
      points_tail = []
      for i in range(len(x_Arai_segment)): # uses only the best_fit segment, so no need for further selection
          x = x_Arai_segment[i]
          y = y_Arai_segment[i]
          points.append((x, y))
          points_arai.append((x,y))
  
      for num, temp in enumerate(ptrm_checks_temperatures): #
          if temp >= tmin and temp <= tmax: # if temp is within selected range
              if (ptrm_checks_starting_temperatures[num] >= tmin and
                      ptrm_checks_starting_temperatures[num] <= tmax): # and also if it was not done after an out-of-range temp
                  x = x_ptrm_check[num]
                  y = y_ptrm_check[num]
                  points.append((x, y))
                  points_ptrm.append((x,y))
  
      for num, temp in enumerate(tail_checks_temperatures):
          if temp >= tmin and temp <= tmax:
              if (tail_checks_starting_temperatures[num] >= tmin and
                      tail_checks_starting_temperatures[num] <= tmax):
                  x = x_tail_check[num]
                  y = y_tail_check[num]
                  points.append((x, y))
                  points_tail.append((x,y))
             # print "points (tail checks added)", points
      fancy_points = {'points_arai': points_arai, 'points_ptrm': points_ptrm, 'points_tail': points_tail}
      return points, fancy_points
docs: returns relevant points for a SCAT test