def dxi(xi, n1, n2):
""" sigma(xi) = sqrt(
(dxi / dAfb * sigma(Alr)) ** 2 +
(dxi / ddPhi * sigma(dPhi)) ** 2 + ...
2*(dxi / dAfb) * (dxi / ddPhi) * cov(Afb, dPhi) + ...
)
"""
data = readFitRes('upol', 0., False)
p0, p1 = sigmaCoefs(data)
return np.sqrt((p0 + xi**2) / n1 + p1 * xi**2 / n2)
def main():
""" Unit test """
xi = 1.
print(' Abf: {:.3f}'.format(dxidAfb()))
print(' dphi: {:.3f}'.format(dxidDphi(xi)))
print('alpha: {:.3f}'.format(dxidAlpha(xi)))
print('alph1: {:.3f}'.format(dxidAlph1(xi)))
print(' dxi: {:.3f}'.format(dxi(xi, 1, 1)))
data = readFitRes('upol', 0., False)
p0, p1 = sigmaCoefs(data)
print('p0: {:.3f}, p1: {:.3f}'.format(p0, p1))
xi = [0.4, 0.6, 0.8]
asymPrecision('fb', xi)
def main():
""" Unit test """
xi = 1.
print(' Alr: {:.3f}'.format(dxidAlr()))
print(' dphi: {:.3f}'.format(dxidDphi(xi)))
print('alpha: {:.3f}'.format(dxidAlpha(xi)))
print('alph1: {:.3f}'.format(dxidAlph1(xi)))
print(' dxi: {:.3f}'.format(dxi(xi, 1, 1)))
data = readFitRes('upol', 0., False)
p0, p1 = sigmaCoefs(data)
print('p0: {:.3f}, p1: {:.3f}'.format(p0, p1))
import matplotlib.pyplot as plt
plt.rc('font', size=14)
xi = [0.4, 0.6, 0.8]
asymPrecision('lr', xi)
def loadData(xi):
""" """
return readFitRes('ss3d', xi, False)
def loadData():
""" """
return readFitRes('upol', 0., False)
def loadData(xi):
""" """
return readFitRes('full', xi, False)