def calb(bdeg, sigb, ldeg, sigl, dkpc, sigd, Har):
b = bdeg * par.degtorad
l = ldeg * par.degtorad
zkpc = dkpc * math.sin(b)
if zkpc < 0.0:
zkpcm = -zkpc
else:
zkpcm = zkpc
adrcold = aplold(dkpc, b, l) * math.cos(b) #s^-1
errDT91 = err_DT91(bdeg, sigb, ldeg, sigl, dkpc, sigd) #s^-1
azbchfh = fhigh(zkpc) * math.sin(b) * 1.08e-19 #s^-1
azbchfl = flow(zkpc) * math.sin(b) * 1.08e-19 #s^-1
errhi = errHFhi(bdeg, sigb, dkpc, sigd) #s^-1
errlo = errHFlo(bdeg, sigb, dkpc, sigd) #s^-1
if Har == 1:
if zkpcm <= 1.5:
print("Excess_parallel_DT91, Excess_z_HF04fit = ", adrcold, ", ",
azbchfl)
else:
print("Excess_parallel_DT91, Excess_z_HF04fit = ", adrcold, ", ",
azbchfh)
else:
if zkpcm <= 1.5:
print("Excess_parallel_DT91, Excess_z_HF04fit = ", adrcold, "+/-",
errDT91, ", ", azbchfl, "+/-", errlo)
else:
print("Excess_parallel_DT91, Excess_z_HF04fit = ", adrcold, "+/-",
errDT91, ", ", azbchfh, "+/-", errhi)
return None
def caleb(bdeg, sigb, ldeg, sigl, dkpc, sigd, Har):
b = bdeg * par.degtorad
l = ldeg * par.degtorad
zkpc = dkpc * math.sin(b)
adrcold = aplold(dkpc, b, l) * math.cos(b) #s^-1
errDT91 = err_DT91(bdeg, sigb, ldeg, sigl, dkpc, sigd) #s^-1
Excz = MWBHZfo(bdeg, ldeg, dkpc) #s^-1
print("Excess_parallel_DT91, Excess_z_galpy(with BH) = ", adrcold, ", ",
Excz)
return None
def calea(bdeg, sigb, ldeg, sigl, dkpc, sigd, Har):
global excpl, exz
b = bdeg * par.degtorad
l = ldeg * par.degtorad
zkpc = dkpc * math.sin(b)
adrcold = aplold(dkpc, b, l) * math.cos(b) #s^-1
errDT91 = err_DT91(bdeg, sigb, ldeg, sigl, dkpc, sigd) #s^-1
Excz = MWZfo(bdeg, ldeg, dkpc)
print("Excess_parallel_DT91, Excess_z_galpy(without BH) = ", adrcold, ", ",
Excz)
excpl = adrcold
exz = Excz
return None
def cala(bdeg, sigb, ldeg, sigl, dkpc, sigd, Har):
b = bdeg * par.degtorad
l = ldeg * par.degtorad
zkpc = dkpc * math.sin(b)
adrcold = aplold(dkpc, b, l) * math.cos(b) #s^-1
errDT91 = err_DT91(bdeg, sigb, ldeg, sigl, dkpc, sigd) #s^-1
azbcnt = g(zkpc) * math.sin(b) #s^-1
errnt = errNT(bdeg, sigb, dkpc, sigd) #s^-1
if Har == 1:
print("Excess_parallel_DT91, Excess_z_NT95 = ", adrcold, ", ", azbcnt)
else:
print("Excess_parallel_DT91, Excess_z_NT95 = ", adrcold, "+/-",
errDT91, ", ", azbcnt, "+/-", errnt)
return None