EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH (CERN) CERN-PH-EP-2014-251 LHCb-PAPER-2014-054 14 October 2014 Search for CP violation in D0 ! π−π+π0 decays with the energy test The LHCb collaborationy Abstract A search for time-integrated CP violation in the Cabibbo-suppressed decay D0 ! π−π+π0 is performed using for the first time an unbinned model-independent technique known as the energy test. Using proton-proton collision data, corresponding to an integrated luminosity of 2:0 fb−1 collected by the p LHCb detector at a centre-of-mass energy of s = 8 TeV, the world's best sensitivity to CP violation in this decay is obtained. The data are found to be consistent with the hypothesis of CP symmetry with a p-value of (2:6 ± 0:5)%. arXiv:1410.4170v2 [hep-ex] 11 Dec 2014 Submitted to Phys. Lett. B c CERN on behalf of the LHCb collaboration, license CC-BY-4.0. yAuthors are listed at the end of this paper. ii 1 Introduction ging of the D mesons is performed through the measurement of the soft pion (πs) charge in the The decay D0 ! π−π+π0 (charge conjugate decays ∗+ 0 + D ! D πs decay. are implied unless stated otherwise) proceeds via a The method exploited in this Letter, called the singly Cabibbo-suppressed c! dud transition with energy test [6, 7], verifies the compatibility of the a possible admixture from a penguin amplitude. observed data with CP symmetry. It is sensitive The interference of these amplitudes may give rise to local CP violation in the Dalitz plot and not to to a violation of the charge-parity symmetry (CP global asymmetries. The unbinned technique ap- violation), which may be observed as an asymmetry plied here is used for the first time. A visualisation in the total rates, or in the distribution of events method is also used that allows identification of over the Dalitz plot. Contributions from particles regions of the Dalitz plot in which CP violation that are not described in the Standard Model (SM) is observed. As this model-independent method and participate in the loops of the penguin ampli- cannot identify which amplitudes contribute to the tude can enhance the O(10−3) CP violation effects observed asymmetry, a model-dependent analysis expected within the SM [1]. Therefore CP violation would be required in the case of a signal for a non- in D0 ! π−π+π0 decays provides sensitivity to such zero CP asymmetry. non-SM physics. In addition to this direct CP violation, time- integrated CP asymmetry in D0 ! π−π+π0 decays 2 Detector and reconstruc- can also receive an indirect contribution arising tion from either the D0-D0 mixing or interference in decays following mixing. While direct CP asym- The LHCb detector [8] is a single-arm forward metry depends on the decay mode, indirect CP spectrometer covering the pseudorapidity range violation is expected to be the same for all CP 2 < η < 5, designed for the study of particles eigenstates. Recent time-dependent measurements containing b or c quarks. The detector includes 0 − + − + of D ! π π , K K decays constrain the indi- a high-precision tracking system consisting of a −3 rect CP asymmetry to the O(10 ) level [2]. silicon-strip vertex detector surrounding the pp in- 0 − + 0 The decay D ! π π π is dominated by the teraction region, a large-area silicon-strip detector ρ(770) resonances, with the ρ meson decaying into located upstream of a dipole magnet with a bend- 0 0 a pair of pions, which lead to the final states ρ π , ing power of about 4 Tm, and three stations of + − − + ρ π , and ρ π . Higher-mass ρ and f0 resonances, silicon-strip detectors and straw drift tubes placed as well as f2(1270) and σ(400) particles only con- downstream of the magnet. The combined tracking tribute with fractions at the percent level or less [3]. system provides a measurement of momentum, p, Singly Cabibbo-suppressed charm decays have with a relative uncertainty that varies from 0.4% recently received a significant attention in the liter- at low momentum to 0.6% at 100 GeV=c. The mini- ature (see e.g. Ref. [1]). A particular interest in the mum distance of a track to a primary vertex, the im- 0 − + 0 decay D ! π π π was pointed out in Ref. [4], pact parameter (IP), is measured with a resolution where the authors derived isospin relations between of (15 + 29=pT) µm, where pT is the component of p the different D ! ρπ amplitudes and discussed the transverse to the beam, in GeV=c. Charged hadrons possibility of identifying contributions from non-SM are identified using two ring-imaging Cherenkov de- physics using CP violation measurements. tectors. Photon, electron and hadron candidates Previously, the most sensitive search for CP viola- are identified by a calorimeter system consisting of tion in this decay was performed by the BaBar col- scintillating-pad and preshower detectors, an elec- laboration [5]. Their result excluded CP -violating tromagnetic calorimeter and a hadronic calorimeter. effects larger than a few percent. The results pre- Muons are identified by a system composed of al- sented here are based on a signal sample that is ternating layers of iron and multiwire proportional about eight times larger and have higher precision. chambers. This is the first CP violation analysis performed The trigger consists of a hardware stage, based on at the LHCb experiment with decays involving π0 high-pT signatures from the calorimeter and muon mesons. systems, followed by a two-level software stage, − + 0 As π π π is a self-conjugate final state and which applies partial event reconstruction. The 0 0 accessible to both D and D decays, flavour tag- software trigger at its first level requires at least 1 one good quality track associated with a particle 3 Event selection 2 having high pT and high χIP, defined as the differ- ence in χ2 of the primary pp interaction vertex (PV) The offline selection is split into a pre-selection, reconstructed with and without this particle. which follows the trigger selection, and a selection A dedicated second-level trigger partially recon- based on a boosted decision tree (BDT) [10, 11]. structs D0 ! π−π+π0 candidates coming from All the D∗+ candidates are required to pass both ∗+ 0 + levels of the software trigger. In addition, the pre- D ! D πs decays using only information from the charged particles. These requirements ensure selection requires more stringent kinematic crite- the suppression of combinatorial background with- ria than those applied in the trigger; in particular 0 out distorting the acceptance in the decay phase pT > 0:5 GeV=c is required for all the D decay space. The π−π+ pair is combined with a pion to products to reduce the combinatorial background. ∗+ 0 form a D candidate, which is accepted if it has pT For resolved π mesons, the di-photon invariant greater than 2:5 GeV=c and a difference of invariant mass has to be within 15 MeV=c2 of the known π0 − + + − + 2 mass, this corresponds to about three times the masses m(π π πs )−m(π π ) < 285 MeV=c . All the charged-particle tracks used must have good m(γγ) resolution. The invariant mass of merged quality, pT > 0:3 GeV=c and p > 3:0 GeV=c, while photons, due to its lower resolution, is required to the charged pions from the D0 decays must also be within the range of 75−195 MeV=c2. The purity 2 − + of the merged π0 sample is nevertheless significantly have high χIP. The π π combination is required to form a good quality secondary vertex significantly higher with respect to the resolved π0, as it bene- displaced from the PV, while the soft pion must fits from large transverse energy and much lower originate from the PV. combinatorial background. The inclusive D∗+ trigger was introduced at the Cross-feed from D0 ! K−π+π0 decays, with a start of 2012. The data collected by the LHCb kaon misidentified as a pion, is reduced with re- experiment in 2012 corresponds to an integrated quirements on the π± particle identification based luminosity of 2 fb−1 of pp collisions collected at a on the ring-imaging Cherenkov detectors. centre-of-mass energy of 8 TeV. The magnetic field The D0 candidates satisfying the above crite- polarity is reversed regularly during the data taking ria and having invariant mass m(π−π+π0) within with approximately half the data collected at each 40(60) MeV=c2 of the known D0 mass are accepted polarity to reduce the overall effect of any charge- in the resolved (merged) sample; this range corre- dependent detection and reconstruction efficiency. sponds to approximately four times the m(π−π+π0) In this analysis, the two π0 categories recon- resolution. The D∗+ candidates, formed with the 0 + structed in LHCb are exploited [9]. These are pions D and πs candidate, have their entire decay chain for which both final state photons are reconstructed refitted requiring that the D∗+ candidate originates separately (resolved pions), as well as pions that from the corresponding PV, and π0 and D0 can- have higher momentum (typically pT > 2 GeV=c) didates have their nominal masses. This improves and thus a smaller opening angle of the two pho- the D∗+ mass resolution and the resolution in the tons (merged pions). These π0 mesons are detected D0 ! π−π+π0 Dalitz plot, while a requirement put in the calorimeter as one merged cluster which is on the fit quality efficiently suppresses the back- further split into two subclusters based on the ex- ground.