Radiative Bhabha scattering tests by Jean-Jacques Blaising 1) -Wh2.8.2, beams = e1, E1, phs_q_scale = 1e-4 GeV cuts = all E > 10.0 GeV and Theta > 10 degree and Theta < 170 degree ["A"] | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | 16 999986 1.4111753E+05 1.63E+02 0.12 1.16 1.94 -Wh2.8.2-extended, beams = e1, E1, phs_q_scale = 1e-4 GeV cuts = all E > 10.0 GeV and Theta > 10 degree and Theta < 170 degree ["A"] | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | 16 999986 1.4302628E+05 1.65E+02 0.12 1.15 2.13 =>Without beam spectrum no difference between Wh2.8.2 and Wh2.8.2-extended 2) -Wh2.8.2, beams = e1, E1 => circe2, phs_q_scale = 1e-4 GeV cuts = all E > 10.0 GeV and Theta > 10 degree and Theta < 170 degree ["A"] | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | 16 999986 6.0234723E+05 1.65E+03 0.27 2.74 0.29 Wh2.8.2-extended, beams = e1, E1 => circe2, phs_q_scale = 1e-4 GeV cuts = all E > 10.0 GeV and Theta > 10 degree and Theta < 170 degree ["A"] | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | 16 999990 5.9987900E+05 1.28E+03 0.21 2.13* 0.47 =>With beam spectrum no difference between Wh2.8.2 and Wh2.8.2-extended Wh2.8.2-quadruple is stll running but it gives thae same results as Wh2.8.2-extended and Wh2.8.2 3)The cross sections without and with beam spectrum are very different; with beam spectrum the cross section is 4.2 larger, does this make sense ? With Wh2.8.2, beams = e1, E1, at 300 GeV the cross section is ~ 200000 fb With Wh2.8.2, beams = e1, E1, at 190 GeV the cross section is ~ 380000 fb convoluting such cross sections with the beam spectrum is unlikely to lead to 600000 fb What can explain such high cross with the beam spectrum. ? More details below, the calculation does not show convergence problems. 1) Wh2.8.2, beams = e1, E1, phs_q_scale = 1e-4 GeV cuts = all E > 10.0 GeV and Theta > 10 degree and Theta < 170 degree ["A"] | Integrator: 400000 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| | VAMP: parameter mismatch, discarding grid file 'decay_proc.m1.vg' 1 399984 1.5018765E+05 2.14E+04 14.25 90.12* 0.03 2 399970 1.4130322E+05 2.02E+03 1.43 9.06* 0.17 3 399962 1.4161518E+05 5.80E+02 0.41 2.59* 0.68 4 399954 1.4138456E+05 3.48E+02 0.25 1.56* 1.06 5 399942 1.4133164E+05 2.58E+02 0.18 1.15* 2.95 6 399928 1.4059233E+05 2.41E+02 0.17 1.08* 3.94 7 399916 1.4084697E+05 2.49E+02 0.18 1.12 2.21 8 399904 1.4132581E+05 2.60E+02 0.18 1.17 2.32 9 399890 1.4123835E+05 2.50E+02 0.18 1.12* 3.05 10 399878 1.4139472E+05 2.81E+02 0.20 1.26 2.37 11 399868 1.4148723E+05 2.74E+02 0.19 1.23* 2.60 12 399852 1.4123048E+05 2.71E+02 0.19 1.21* 2.89 13 399840 1.4117716E+05 2.59E+02 0.18 1.16* 3.40 14 399830 1.4138563E+05 2.36E+02 0.17 1.06* 4.64 15 399822 1.4153760E+05 2.60E+02 0.18 1.16 2.40 |-----------------------------------------------------------------------------| 15 5998540 1.4123235E+05 7.50E+01 0.05 1.30 2.40 0.97 15 |-----------------------------------------------------------------------------| 16 999986 1.4111753E+05 1.63E+02 0.12 1.16 1.94 |=============================================================================| | Time estimate for generating 10000 events: 0d:00h:00m:41s 2)Wh2.8.2, beams = e1, E1 => circe2, phs_q_scale = 1e-4 GeV cuts = all E > 10.0 GeV and Theta > 10 degree and Theta < 170 degree ["A"] | Integrator: 400000 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 399984 3.8191926E+05 8.40E+04 21.99 139.04* 0.01 2 399972 5.8180318E+05 2.76E+04 4.75 30.03* 0.03 3 399964 6.1001290E+05 5.93E+03 0.97 6.15* 0.13 4 399950 6.0643196E+05 2.86E+03 0.47 2.98* 0.35 5 399942 6.0171595E+05 1.96E+03 0.33 2.06* 0.84 6 399934 6.0428321E+05 1.76E+03 0.29 1.84* 0.83 7 399920 6.0199620E+05 1.76E+03 0.29 1.85 0.89 8 399904 6.0242771E+05 1.87E+03 0.31 1.97 0.83 9 399892 6.0026475E+05 1.77E+03 0.29 1.86* 0.84 10 399882 6.0351273E+05 1.78E+03 0.30 1.87 0.85 11 399870 6.0429589E+05 2.03E+03 0.34 2.12 0.62 12 399860 6.0191612E+05 1.94E+03 0.32 2.04* 0.90 13 399850 5.9966567E+05 1.85E+03 0.31 1.95* 0.83 14 399836 6.1205491E+05 8.75E+03 1.43 9.04 0.11 15 399826 6.0337752E+05 2.25E+03 0.37 2.36* 0.76 15 399826 6.0337752E+05 2.25E+03 0.37 2.36* 0.76 |-----------------------------------------------------------------------------| 15 5998586 6.0256054E+05 5.78E+02 0.10 2.35 0.76 1.34 15 |-----------------------------------------------------------------------------| 16 999986 6.0234723E+05 1.65E+03 0.27 2.74 0.29 |=============================================================================| | Time estimate for generating 10000 events: 0d:00h:03m:25s 3)Wh2.8.2-extended, beams = e1, E1, phs_q_scale = 1e-4 GeV cuts = all E > 10.0 GeV and Theta > 10 degree and Theta < 170 degree ["A"] | Integrator: 400000 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| | VAMP: parameter mismatch, discarding grid file 'decay_proc.m1.vg' 1 399984 1.5519016E+05 1.36E+04 8.76 55.40* 0.04 2 399970 1.4044145E+05 1.24E+03 0.88 5.59* 0.28 3 399956 1.4261852E+05 6.87E+02 0.48 3.04* 0.40 4 399942 1.4417457E+05 5.87E+02 0.41 2.57* 0.48 5 399928 1.4317257E+05 3.42E+02 0.24 1.51* 0.98 6 399916 1.4300038E+05 2.55E+02 0.18 1.13* 2.81 7 399904 1.4315115E+05 2.52E+02 0.18 1.11* 2.26 8 399894 1.4280278E+05 2.40E+02 0.17 1.06* 4.03 9 399884 1.4302509E+05 2.37E+02 0.17 1.05* 3.73 10 399874 1.4318175E+05 2.30E+02 0.16 1.02* 4.22 11 399868 1.4322374E+05 2.91E+02 0.20 1.29 2.14 12 399858 1.4306530E+05 2.86E+02 0.20 1.27* 2.72 13 399848 1.4320031E+05 2.61E+02 0.18 1.15* 2.63 14 399834 1.4270654E+05 2.43E+02 0.17 1.08* 3.56 15 399824 1.4306788E+05 2.55E+02 0.18 1.13 2.24 |-----------------------------------------------------------------------------| 15 5998484 1.4304563E+05 7.67E+01 0.05 1.31 2.24 0.98 15 |-----------------------------------------------------------------------------| 16 999986 1.4302628E+05 1.65E+02 0.12 1.15 2.13 |=============================================================================| | Time estimate for generating 10000 events: 0d:00h:01m:34s 4)Wh2.8.2-extended, beams = e1, E1 => circe2, phs_q_scale = 1e-4 GeV cuts = all E > 10.0 GeV and Theta > 10 degree and Theta < 170 degree ["A"] | Integrator: Using VAMP channel equivalences | Integrator: 400000 initial calls, 20 bins, stratified = T | Integrator: VAMP |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| | VAMP: parameter mismatch, discarding grid file 'decay_proc.m1.vg' 1 399984 5.7732822E+05 1.26E+05 21.86 138.26* 0.02 2 399972 5.6013173E+05 2.30E+04 4.10 25.93* 0.03 3 399962 5.9923881E+05 4.94E+03 0.82 5.21* 0.23 4 399952 5.9870507E+05 2.20E+03 0.37 2.32* 0.73 5 399942 6.0016911E+05 1.92E+03 0.32 2.02* 0.63 6 399926 5.9887922E+05 2.07E+03 0.34 2.18 0.46 7 399912 6.0002644E+05 1.80E+03 0.30 1.90* 0.90 8 399900 5.9685904E+05 1.80E+03 0.30 1.91 0.73 9 399886 5.9621116E+05 1.64E+03 0.27 1.74* 1.34 10 399874 5.9724524E+05 1.62E+03 0.27 1.72* 1.16 11 399858 5.9607975E+05 1.79E+03 0.30 1.90 0.60 12 399840 5.9651138E+05 1.81E+03 0.30 1.92 0.80 13 399826 5.9643940E+05 1.67E+03 0.28 1.77* 1.17 14 399814 5.9334115E+05 1.54E+03 0.26 1.64* 1.29 15 399802 5.9963237E+05 2.10E+03 0.35 2.22 0.56 |-----------------------------------------------------------------------------| 15 5998450 5.9720839E+05 5.17E+02 0.09 2.12 0.56 1.25 15 |-----------------------------------------------------------------------------| 16 999990 5.9987900E+05 1.28E+03 0.21 2.13* 0.47 |=============================================================================| | Time estimate for generating 10000 events: 0d:00h:05m:53s