Pion Form Factor Measurement at BESIII Martin Ripka on behalf of the BESIII colaboration EINN - November, 2015 1 / 17
Motivation: Why Form Factor Measurements at BESIII? Muon anomalous magnetic moment a µ = (g µ 2)/2 Experimental measurement at BNL: aµ exp = 11659208.9(5.4)(3.3) 10 10 [PRD 73 072003 (2006)] Theoretical calculation: aµ theo = 11659182.8(4.9) 10 10 [J. Phys. G 38, 085003 (2011)] Theory and experiment not in agreement: aµ exp aµ theo (30 ± 8) 10 10 3 4 σ deviation M. Davier, A. Hoecker, B. Malaescu and Z. Zhang, Eur. Phys. J. C 71 1515 (2011) 2 / 17
Theoretical calculation of a µ aµ theo = aµ QED + aµ weak + aµ QCD aµ QED = (11658471.8 ± 0.0) 10 10 [PRL 109, 111808 (2012)] aµ weak = (15.36 ± 0.1) 10 10 [PRD 88, 053005 (2013)] aµ QCD = aµ LbL + aµ VP,LO + aµ VP,HO aµ VP,LO = (694.9 ± 4.2) 10 10 [J. Phys. G 38, 085003 (2011)] aµ VP,HO = ( 98.4 ± 0.7) 10 10 [J. Phys. G 38, 085003 (2011)] aµ LbL = (11.6 ± 3.9) 10 10 [Phys. Rept. 477, 1 (2009)] B-field B-field B-field had had 3 / 17
The Vacuum Polarisation Contribution to a QCD µ Loop can not be calculated for low momentum hadrons Optical theorem connects VP amplitude with hadronic cross sections: σ(s) e + e hadrons = 4πα s Im Π γ (s) 2 had had photon self-energy function aµ VP,LO = 1 4π 3 0 ds K(s)σ e + e hadrons(s) Hadronic contributions to a QCD µ : hadronic cross-section M. Davier, A. Hoecker, B. Malaescu and Z. Zhang, Eur. Phys. J. C 71 1515 (2011) 4 / 17
Hadronic Final States contributing to a VP,LO µ D. Bernard [BaBar Collaboration], PoS Hadron 2013, 126 (2013) [arxiv:1402.0618 [hep-ex]]. Most important channels: π + π, KK, π + π π 0, π + π 2π 0 Largest contribution to uncertainty: π + π, π + π 2π 0, KK 5 / 17
π + π at Babar and Kloe M. Davier, A. Hoecker, B. Malaescu and Z. Zhang, Eur. Phys. J. C 71 1515 (2011) Babar and Kloe each claim sub-percent precision Measurements do not agree with each other Another high precision measurement needed BESIII 6 / 17
BEPCII τ-charm factory Energy range: 2-4.6 GeV Design luminosity: 10 33 cm 2 s 1 (at 3.77 GeV) Linac + double storage ring 7 / 17
BESIII Detector Multilayer Drift Chamber (MDC) Time of Flight system (ToF) Electromagnetic Calorimeter (EMC) Super Conducting magnet 1 Tesla (SC) Resistive Plate Chamber (RPC) for muon detection 8 / 17
Initial State Radiation Technique I Need σ had (s) in the entire energy range where pqcd fails had Initial State Radiation (ISR) reduces the effective CMS-energy of the collision: m 2 had = E 2 CMS 2E CMSE ISR Non radiative cross-section can be obtained by dσ (had+γ) dm had = 2m had W (s, E ISR, θ ISR )σ had s Radiator-function W (s, E ISR, θ ISR ) gives the amplitude to emit an ISR photon 9 / 17
Initial State Radiation Technique II Emission of ISR photons is suppressed by α/π High integrated luminosity needed for precision measurements Untagged analysis possible above 1 GeV Tagged analysis Untagged analysis 10 / 17
Event Selection and Particle Identification (PID) TMVA overtraining check for classifier: CFMlpANN dx (1/N) dn/ Signal (test sample) Signal (training sample) 22 Background (test sample) Background (training sample) 20 Kolmogorov-Smirnov test: signal (background) probability = 0.735 (0.455) 18 16 14 12 10 8 6 4 2 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 CFMlpANN response U/O-flow (S,B): (0.0, 0.0)% / (0.0, 0.0)% Kinematic Fit for π + π γ ISR final state Standard BESIII PID system for electron rejection Artificial Neuronal Network for muon-pion separation Before ANN After ANN events / 20 MeV 25000 20000 - π + π γ MC - µ + µ γ MC data events / 20 MeV 10000 8000 - π + π γ MC - µ + µ γ MC data 15000 6000 10000 4000 5000 2000 0 0 0.5 1 1.5 2 2.5 3 3.5 4 m 2π [GeV] 0 0 0.5 1 1.5 2 2.5 3 3.5 4 m 2π [GeV] 11 / 17
Systematic Uncertainties Source Uncertainty (%) Photon efficiency correction 0.2 Pion tracking efficiency correction 0.3 Pion ANN efficiency correction 0.2 Pion e-pid efficiency correction 0.2 ANN negl. Angular acceptance 0.1 Background subtraction 0.1 Unfolding 0.2 FSR correction δ FSR 0.2 Vacuum polarisation correction δ vac 0.2 Radiator function 0.5 Luminosity L 0.5 Sum 0.9 12 / 17
π + π Cross Section - 1400 )) [nb] FSR 1200 1000 π + π - (γ e + σ bare (e 800 600 400 200 0.6 0.65 0.7 0.75 0.8 0.85 0.9 s' [GeV] σ bare ( s ) = 1 2 s W (s,x)ɛ( dn s s )Lδ vacδ FSR d s ρ-ω interference clearly visible 13 / 17
π + π Form Factor (Gounaris-Sakurai Parametrisation) 45 BESIII fit 40 BESIII 35 2 π F 30 25 20 15 10 0.6 0.65 0.7 0.75 0.8 0.85 0.9 s' [GeV] parameter BESIII value PDG 2014 m ρ [MeV/c 2 ] 776.0 ± 0.4 775.26 ± 0.25 Γ ρ [MeV] 151.7 ± 0.7 147.8 ± 0.9 m ω [MeV/c 2 ] 782.2 ± 0.6 782.65 ± 0.12 Γ ω [MeV] fixed to PDG 8.49 ± 0.08 c ρ [10 3 ] 1.7 ± 0.2 - φ ω [rad] 0.04 ± 0.13-14 / 17
Comparison to Other π + π Measurements 2 2 / BESIII fit - 1 0.15 0.1 0.05 0 BESIII fit KLOE 08 KLOE 10 KLOE 12 / BESIII fit - 1 0.15 0.1 0.05 0 BESIII fit BaBar BESIII π F -0.05 π F -0.05-0.1-0.1-0.15 0.6 0.65 0.7 0.75 0.8 0.85 0.9 s' [GeV] -0.15 0.6 0.65 0.7 0.75 0.8 0.85 0.9 s' [GeV] 0.15 BESIIIMfit M/MBESIIIMfitM-M1 2 π F 0.1 0.05 0-0.05-0.1 SND CMD2-2006 CMD2-2004 New BESIII measurement agrees with KLOE and BaBar Small shift wrt. BaBar above ρ-ω interference -0.15 0.6 0.65 0.7 0.75 0.8 0.85 0.9 s [GeV] 15 / 17
Final Result: Contribution to a VP,LO µ Precision competitive with previous measurements BESIII measurement between BaBar and KLOE a ππ,lo µ (600 900 MeV) = (370.0 ± 2.5 stat ± 3.3 sys ) 10 10 Confirms deviation of 3.4σ between experiment and theory arxiv:1507.08188 and submitted to PLB 16 / 17
Outlook Extend tagged π + π ISR study to threshold region Use Untagged ISR technique for π + π cross section at higher energies Analyse π + π form factor from R-scan data (130 points, L 1.3fb 1 ) Ongoing investigation of π + π π 0 and π + π π 0 π 0 Thank You 17 / 17