def gettraj(t, t0, u0, te, piex, piey, alpha, delta, t0par, qlat, qlong, hjds, ras, decs, diss, SPITZ):
tau = (t - t0) / te
beta = u0
qn, qe = getpar.geta(t, alpha, delta, t0par)
if SPITZ == True:
qnp, qep, qr = getpar.gets(t, hjds, ras, decs, diss, alpha, delta)
else:
qnp, qep = getpar.gett(t, qlat, qlong, alpha, delta)
qn += qnp
qe += qep
dtau = piex * qn + piey * qe
dbeta = piex * qe - piey * qn
return tau + dtau, beta + dbeta
def getqnqe(dates, alpha, delta, t0par, qlat, qlong, SPITZ):
qns, qes = [], []
qnqe = []
for t in dates:
qn, qe = getpar.geta(t, alpha, delta, t0par)
if SPITZ == True:
qnp, qep, qr = getpar.gets(t, hjds, ras, decs, diss, alpha, delta)
else:
qnp, qep = getpar.gett(t, qlat, qlong, alpha, delta)
qn += qnp
qe += qep
qnqe.append([qn, qe])
qnqe = np.array(qnqe).T
return qnqe