! plot_twist_lrs_limit scale_xlo = 200 scale_xhi = 1200 scale_ylo = -0.1 scale_yhi = 0.1 scale scale_xlo scale_xhi 5 scale_ylo scale_yhi 4 x_mass = [scale_xlo;scale_xhi] y_zeta = [scale_ylo;scale_yhi] goto input_decay_parameter_results set linthk -6 ! use 68.27%CL (one sigma) !cl = 0.6827 ! 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 ! ! 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 ! lrs model, limit on zeta from TWIST rho; rho>=0.75, so use one-sided limit llim_rho_twist_macd = 0.75-sig_clg1*drho_twist_macd ! note this doesn't use rho_twist zeta_rho_twist_macd_1[1] = -sqrt((0.75 - llim_rho_twist_macd)*2./3.) zeta_rho_twist_macd_1[2] = zeta_rho_twist_macd_1[1] zeta_rho_twist_macd_2[1] = -zeta_rho_twist_macd_1[1] zeta_rho_twist_macd_2[2] = zeta_rho_twist_macd_2[1] ! 2008 est, for final results rho_twist_08fp = 0.7500 drho_twist_08fp = 1.e-4*sqrt(1.3**2 + 2.4**2 + 1.1**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 ! lrs model, limit on zeta from final TWIST rho llim_rho_twist_08fp = 0.75-sig_clg1*drho_twist_08fp zeta_rho_twist_1_08fp[1] = -sqrt((0.75 - llim_rho_twist_08fp)*2./3.) zeta_rho_twist_1_08fp[2] = zeta_rho_twist_1_08fp[1] zeta_rho_twist_2_08fp[1] = -zeta_rho_twist_1_08fp[1] zeta_rho_twist_2_08fp[2] = zeta_rho_twist_2_08fp[1] ! 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 ! Jamieson et al., 2006 pmuxi_twist_bj = 1.0003 dpmuxi_twist_bj = 1.e-4*sqrt(6**2 + 38**2) llim_pmuxi_twist_bj = pmuxi_twist_bj - sig_clg*dpmuxi_twist_bj ulim_pmuxi_twist_bj = pmuxi_twist_bj + sig_clg*dpmuxi_twist_bj ! estimated final result, 2008 pmuxi_twist_08fp = 1.000 dpmuxi_twist_08fp = 1.e-4*sqrt(3.**2 + 5.9**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 ! pmu*xi*delta/rho ! put in numbers from MacDonald explicitly llim_pmuxi_twist_js = 0.99524 ulim_pmuxi_twist_js = 1.00091 input_decay_parameter_results: @decay_parameter_results @pmuxi_lim llim_pmuxi_bel x_pmuxi_bel = massvar_nonmanif y_pmuxi_bel = zetavar_nonmanif x_pmuxi_bel_m = massvar_manif y_pmuxi_bel_m = zetavar_manif @pmuxi_lim llim_pmuxi_twist_jami x_pmuxi_twist_jami = massvar_nonmanif y_pmuxi_twist_jami = zetavar_nonmanif x_pmuxi_twist_jami_m = massvar_manif y_pmuxi_twist_jami_m = zetavar_manif @pmuxi_lim llim_pmuxi_twist_macd x_pmuxi_twist_macd = massvar_nonmanif y_pmuxi_twist_macd = zetavar_nonmanif x_pmuxi_twist_macd_m = massvar_manif y_pmuxi_twist_macd_m = zetavar_manif @pmuxi_lim llim_pmuxi_twist_08fp x_pmuxi_twist_08fp = massvar_nonmanif y_pmuxi_twist_08fp = zetavar_nonmanif x_pmuxi_twist_08fp_m = massvar_manif y_pmuxi_twist_08fp_m = zetavar_manif ! set up plot shape, size, etc. enable border orientation portrait window 21 0 20 100 100 !window 21 0 0 100 100 window 21 set xtica 90 set ytica -90 if (~exist(`letter_size')) then letter_size = 3 set %txthit letter_size set %xlabsz letter_size set %ylabsz letter_size set %xnumsz letter_size set %ynumsz letter_size set %xlaxis 25 set %ylaxis 25 pl_xlab = `(g<_>L<^>/g<_>R<^>)m<_>2<^> (GeV/c<^>2<_>)' pl_ylab = `(g<_>R<^>/g<_>L<^>)' pl_titl = `General LRS model' ! pl_titl = ` ' set histyp 0 label\xaxis pl_xlab label\yaxis pl_ylab clear legend frame 30 77 85 87 legend frame off !!! plot_beltrami: legend_text = - `Beltrami et al. P<_><^> ' set lintyp 4 color white legend on gra\noaxes legend_text x_pmuxi_bel y_pmuxi_bel legend off !!! plot_twist: ! plot TWIST current plot_jami: legend frame 30 21 95 37 legend_text = - `TWIST P<_><^> (Jamieson) ' set lintyp 3 color blue legend on gra\noaxes legend_text x_pmuxi_twist_jami y_pmuxi_twist_jami legend off !!! plot_twist_macd: legend_text = - `TWIST , (MacDonald) + Jodidio' set lintyp 5 color blue legend on gra\noaxes legend_text x_pmuxi_twist_macd y_pmuxi_twist_macd legend off !!! plot_08fp: ! plot TWIST expected final values as of 2008 legend_text = - `TWIST , P<_><^> estimated final ' set lintyp 1 color red legend on gra\noaxes legend_text x_pmuxi_twist_08fp y_pmuxi_twist_08fp legend off window\noconfirm 99 70 64 95 67 erasewindow window 21 color white set %xloc 90 set %yloc 57 set linthk -10 set cursor -11 text\noconfirm `allowed' graph\axesonly x_mass y_zeta