! compilation of results for muon decay parameters

! use 90%CL
cl = 0.9

! sigma corresponding to this cl
! usual gaussian, 1.64@90%
sig_clg = sqrt(2)*aerf(cl)
! or one-sided limit, 1.28@90%
sig_clg1 = sqrt(2)*aerf(1 - 2*(1 - cl))

!!! =================== rho

! Derenzo PR 181 (1969) 1954 0.7518 pm 0.0026
rho_der = 0.7518
drho_der = 1.e-4*26.
llim_rho_der = rho_der-sig_clg*drho_der
ulim_rho_der = rho_der+sig_clg*drho_der

! macdonald et al. (2008) 0.75014 pm 0.00017 pm 0.00044 pm 0.00011
rho_twist_macd = 0.75014
drho_twist_macd = 1.e-4*sqrt(1.7**2 + 4.6**2 + 1.1**2)
llim_rho_twist_macd = rho_twist_macd-sig_clg*drho_twist_macd
ulim_rho_twist_macd = rho_twist_macd+sig_clg*drho_twist_macd

! TWIST estimated final result, 2008  SM pm 0.00013 pm 0.00024 (pm 0.00011 for eta?)
rho_twist_08fp = 0.7500
drho_twist_08fp = 1.e-4*sqrt(1.3**2 + 2.4**2)
llim_rho_twist_08fp = rho_twist_08fp-sig_clg*drho_twist_08fp
ulim_rho_twist_08fp = rho_twist_08fp+sig_clg*drho_twist_08fp

!!!=========================  delta

! Balke et al. 1988 0.7486 pm 0.0026 pm 0.0028
delta_bal = 0.7486
ddelta_bal = 1.e-3*sqrt(2.6**2 + 2.8**2)
llim_delta_bal = delta_bal - sig_clg*ddelta_bal
ulim_delta_bal = delta_bal + sig_clg*ddelta_bal

! Gaponenko et al. 2005  0.74964 pm 0.00066 pm 0.00112
delta_twist_gap = 0.74964
ddelta_twist_gap = 1.e-4*sqrt(6.6**2 + 11.2**2)
llim_delta_twist_gap = delta_twist_gap - sig_clg*ddelta_twist_gap
ulim_delta_twist_gap = delta_twist_gap + sig_clg*ddelta_twist_gap

! MacDonald et al., (2008) 0.75067 pm 0.00030 pm 0.00067
delta_twist_macd = 0.75067
ddelta_twist_macd = 1.e-4*sqrt(3.0**2 + 6.7**2)
llim_delta_twist_macd = delta_twist_macd - sig_clg*ddelta_twist_macd
ulim_delta_twist_macd = delta_twist_macd + sig_clg*ddelta_twist_macd

! TWIST estimated final result, 2008  SM pm 0.00023 pm 0.00022
delta_twist_08fp = 0.7500
ddelta_twist_08fp = 1.e-4*sqrt(2.3**2 + 3.2**2)
llim_delta_twist_08fp = delta_twist_08fp - sig_clg*ddelta_twist_08fp
ulim_delta_twist_08fp = delta_twist_08fp + sig_clg*ddelta_twist_08fp

!!!========================== pmu*xi

! Beltrami et al., PLB 194, 326 (1987)
pmuxi_bel = 1.0027
dpmuxi_bel = 1.e-4*sqrt(79**2 + 30**2)
llim_pmuxi_bel = pmuxi_bel - sig_clg*dpmuxi_bel
ulim_pmuxi_bel = pmuxi_bel + sig_clg*dpmuxi_bel

! Gaponenko et al.2005 (indirect, 90%cl)
llim_pmuxi_twist_gap = 0.9960 
ulim_pmuxi_twist_gap = 1.0040

! Jamieson et al., 2006
pmuxi_twist_jami = 1.0003
dpmuxi_twist_jami = 1.e-4*sqrt(6**2 + 38**2)
llim_pmuxi_twist_jami = pmuxi_twist_jami - sig_clg*dpmuxi_twist_jami
ulim_pmuxi_twist_jami = pmuxi_twist_jami + sig_clg*dpmuxi_twist_jami

! MacDonald et al., 2008 (indirect, 90%cl) 0.99524-1.00091
llim_pmuxi_twist_macd = 0.99524 
ulim_pmuxi_twist_macd = 1.00091

! TWIST estimated final result, 2008 SM pm 0.00030 pm 0.00060
pmuxi_twist_08fp = 1.000
dpmuxi_twist_08fp = 1.e-4*sqrt(3.0**2 + 6.0**2)
llim_pmuxi_twist_08fp = pmuxi_twist_08fp - sig_clg*dpmuxi_twist_08fp
ulim_pmuxi_twist_08fp = pmuxi_twist_08fp + sig_clg*dpmuxi_twist_08fp

!!! global analysis RH

! MacDonald et al., 2008
!QRR       = 0.000505 +- 0.000341
!QLR       = 0.000859 +- 0.000405
qrr_macd = 0.000505
dqrr_macd = 0.000341
qlr_macd = 0.000859
dqlr_macd = 0.000405
qrr_lim_macd = qrr_macd + sig_clg1*dqrr_macd
qlr_lim_macd = qlr_macd + sig_clg1*dqlr_macd
qmur_lim_macd = qrr_lim_macd + qlr_lim_macd

!qrr_lim_macd = 0.00096
!qlr_lim_macd = 0.00138
!qmur_lim_macd = 0.0024

! estimated final, 2008
qrr_08fp = 0.000327
dqrr_08fp = 0.000223
qlr_08fp = 0.000508
dqlr_08fp = 0.000259
qrr_lim_08fp = qrr_08fp + sig_clg1*dqrr_08fp
qlr_lim_08fp = qlr_08fp + sig_clg1*dqlr_08fp
qmur_lim_08fp = qrr_lim_08fp + qlr_lim_08fp