def calib(bdeg, sigb, ldeg, sigl, dkpc, sigd, Har):
global excpl, exz
b = bdeg * par.degtorad
l = ldeg * par.degtorad
zkpc = dkpc * math.sin(b)
if zkpc < 0.0:
zkpcm = -zkpc
else:
zkpcm = zkpc
azbchfh = fhigh(zkpc) * math.sin(b) * 1.08100761142e-19 #s^-1
azbchfl = flow(zkpc) * math.sin(b) * 1.08100761142e-19 #s^-1
errhi = errHFhi(bdeg, sigb, dkpc, sigd) #s^-1
errlo = errHFlo(bdeg, sigb, dkpc, sigd) #s^-1
Excplbh = MWBHpl(bdeg, ldeg, dkpc) #s^-1
if zkpcm <= 1.5:
print(
"Excess_parallel_galpy(from Vp/Vs, with BH), Excess_z_HF04fit = ",
Excplbh, ", ", azbchfl)
excpl = Excplbh
exz = azbchfl
else:
print(
"Excess_parallel_galpy(from Vp/Vs, with BH), Excess_z_HF04fit = ",
Excplbh, ", ", azbchfh)
excpl = Excplbh
exz = azbchfh
return None
def calia(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
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
Excpl = MWpl(bdeg, ldeg, dkpc) #s^-1
if zkpcm <= 1.5:
print(
"Excess_parallel_galpy(from Vp/Vs, without BH), Excess_z_HF04fit = ",
Excpl, ", ", azbchfl)
else:
print(
"Excess_parallel_galpy(from Vp/Vs, without BH), Excess_z_HF04fit = ",
Excpl, ", ", azbchfh)
return None
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 cald(bdeg, sigb, ldeg, sigl, dkpc, sigd, Har):
global excpl, exz, errpl, errz
b = bdeg * par.degtorad
l = ldeg * par.degtorad
zkpc = dkpc * math.sin(b)
if zkpc < 0.0:
zkpcm = -zkpc
else:
zkpcm = zkpc
adrc = aplmod(dkpc, b, l) * math.cos(b) #s^-1
errReid = err_Reid14(bdeg, sigb, ldeg, sigl, dkpc, sigd) #s^-1
azbchfh = fhigh(zkpc) * math.sin(b) * 1.08100761142e-19 #s^-1
azbchfl = flow(zkpc) * math.sin(b) * 1.08100761142e-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_Reid2014, Excess_z_HF04fit = ", adrc, ", ",
azbchfl)
print("Vp/Vs= ", Vprat(Rpkpcfunc(dkpc, b, l, par.Rskpc)))
excpl = adrc
exz = azbchfl
else:
print("Excess_parallel_Reid2014, Excess_z_HF04fit = ", adrc, ", ",
azbchfh)
print("Vp/Vs= ", Vprat(Rpkpcfunc(dkpc, b, l, par.Rskpc)))
excpl = adrc
exz = azbchfh
else:
if zkpcm <= 1.5:
print("Excess_parallel_Reid2014, Excess_z_HF04fit = ", adrc, "+/-",
errReid, ", ", azbchfl, "+/-", errlo)
print("Vp/Vs= ", Vprat(Rpkpcfunc(dkpc, b, l, par.Rskpc)))
excpl = adrc
exz = azbchfl
errpl = errReid
errz = errlo
else:
print("Excess_parallel_Reid2014, Excess_z_HF04fit = ", adrc, "+/-",
errReid, ", ", azbchfh, "+/-", errhi)
print("Vp/Vs= ", Vprat(Rpkpcfunc(dkpc, b, l, par.Rskpc)))
excpl = adrc
exz = azbchfh
errpl = errReid
errz = errhi
return None
def caljb(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
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
ExcRforceBH = MWBHRfo(bdeg,ldeg,dkpc) #s^-1
if zkpcm<=1.5:
print ("Excess_parallel(galpy-Rforce,with BH), Excess_z_HF04fit = ", ExcRforceBH,", ", azbchfl)
else:
print ("Excess_parallel(galpy-Rforce,with BH), Excess_z_HF04fit = ", ExcRforceBH,", ", azbchfh)
return None;