A search for timeintegrated $CP$ violation in the Cabibbosuppressed decay $ D^0\rightarrow\pi^+\pi^\pi^+\pi^$ is performed using an unbinned, modelindependent technique known as the energy test. This is the first application of the energy test in fourbody decays. The search is performed for $P$even $CP$ asymmetries and, for the first time, is extended to probe the $P$odd case. Using protonproton collision data corresponding to an integrated luminosity of 3.0 fb$^{1}$ collected by the LHCb detector at centreofmass energies of $\sqrt{s}=$7 TeV and 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 $(4.6\pm0.5)\%$ in the $P$even case, and marginally consistent with a $p$value of $(0.6\pm0.2)\%$ in the $P$odd case, corresponding to a significance for $CP$ nonconservation of 2.7 standard deviations.
Distribution of $\Delta m$ with fit overlaid for the selected $ D ^{*+}$ candidates in the 2012 data. The data points and the contributions from signal, background, and their total obtained from the fit are shown. 
Fig1.pdf [26 KiB] HiDef png [268 KiB] Thumbnail [230 KiB] *.C file 

(a,b) Distribution of permutation $T$values fitted with a GEV function for the simulated sample and showing the measured $T$value as a vertical line, and (c,d,e,f) local asymmetry significances. Left column plots are for a $P$even $ C P$ violation test with a $3^\circ$ phase $ C P$ violation introduced in the $a_1(1260)^+$ resonance (see text), projected onto the (c) $m(\pi_1 \pi_2 \pi_3)$ and (e) $m(\pi_1\pi_2)$ axes. Right column plots are for a $P$odd $ C P$ violation test with $3^\circ$ phase $ C P$ violation introduced in the Pwave $\rho ^0 (770) \rho ^0 (770)$ resonance projected onto the same axes. In plots (c,d,e,f) the grey area corresponds to candidates with a contribution to the $T$value of less than one standard deviation. The pink (blue) area corresponds to candidates with a positive (negative) contribution to the $T$value. Light, medium or dark shades of pink and blue correspond to between one and two, two and three, and more than three standarddeviation contributions, respectively. 
Fig2a.pdf [28 KiB] HiDef png [163 KiB] Thumbnail [125 KiB] *.C file 

Fig2b.pdf [28 KiB] HiDef png [157 KiB] Thumbnail [126 KiB] *.C file 

Fig2c.pdf [17 KiB] HiDef png [230 KiB] Thumbnail [210 KiB] *.C file 

Fig2d.pdf [16 KiB] HiDef png [198 KiB] Thumbnail [193 KiB] *.C file 

Fig2e.pdf [17 KiB] HiDef png [233 KiB] Thumbnail [216 KiB] *.C file 

Fig2f.pdf [16 KiB] HiDef png [204 KiB] Thumbnail [201 KiB] *.C file 

(a,b) Distribution of permutation $T$values fitted with a GEV function and showing the $T$value of the data tests as a vertical line, and (c,d,e,f) local asymmetry significances. Left column plots are for the $P$even $ C P$ violation test, projected onto the (c) $m(\pi_1 \pi_2 \pi_3)$ and (e) $m(\pi_1\pi_2)$ axes. Right column plots are for the $P$odd $ C P$ violation test projected onto the same axes. In plots (c,d,e,f) the grey area correspond to candidates with a contribution to the $T$value of less than one standard deviation. In the $P$even $ C P$ violation test the positive (negative) asymmetry significance is set for the $ D ^0$ candidates having positive (negative) contribution to the measured $T$ value. In the $P$odd $ C P$ violation test the positive (negative) asymmetry significance is set for sample $I+IV$ having positive (negative) contribution to the measured $T$ value (see Sect. 5). The pink (blue) area corresponds to candidates with a positive (negative) contribution to the $T$value. Light, medium or dark shades of pink and blue correspond to between one and two, two and three, and more than three standard deviation contributions, respectively. 
Fig3a.pdf [29 KiB] HiDef png [185 KiB] Thumbnail [140 KiB] *.C file 

Fig3b.pdf [29 KiB] HiDef png [183 KiB] Thumbnail [140 KiB] *.C file 

Fig3c.pdf [16 KiB] HiDef png [176 KiB] Thumbnail [172 KiB] *.C file 

Fig3d.pdf [16 KiB] HiDef png [182 KiB] Thumbnail [177 KiB] *.C file 

Fig3e.pdf [16 KiB] HiDef png [187 KiB] Thumbnail [182 KiB] *.C file 

Fig3f.pdf [16 KiB] HiDef png [192 KiB] Thumbnail [187 KiB] *.C file 

Animated gif made out of all figures. 
PAPER2016044.gif Thumbnail 
The yields of signal events in the four samples that obtained from fits to the $\Delta m$ distribution. 
Table_1.pdf [21 KiB] HiDef png [23 KiB] Thumbnail [12 KiB] tex code 

Overview of sensitivities to various $ C P$ violation scenarios in simulation. $\Delta A$ and $\Delta \phi$ denote, respectively, the relative change in magnitude and the change in phase of the amplitude of the resonance $R$. The Pwave $\rho ^0 (770)\rho ^0 (770)$ is a $P$odd component. The phase change in this resonance is tested with the $P$odd $ C P$ violation test. Results for all other scenarios are given with the standard $P$even test. 
Table_2.pdf [54 KiB] HiDef png [134 KiB] Thumbnail [61 KiB] tex code 

Results for the $P$even and $P$odd $ C P$ violation tests for three different values of the effective radius $\delta$. The $p$values obtained with both the counting and GEV fitting methods are given (see text). The counting method is the default method. 
Table_3.pdf [38 KiB] HiDef png [56 KiB] Thumbnail [28 KiB] tex code 
Created on 17 August 2019.Citation count from INSPIRE on 20 August 2019.