def railway_ST(C, start, end, railway_r):
delta = atan2(*(end - C)[0:2]) - atan2(*(start - C)[0:2])
dbg('railway_ST C', C, 'start', start, 'end', end, 'delta', delta)
def railway_ST(C, start, end, railway_r):
delta = atan2(*(end - C)[0:2]) - atan2(*(start - C)[0:2])
dbg('railway_ST C', C, 'start', start, 'end', end, 'delta', delta)