The Dalitz plot analysis technique is used to study the resonant substructures of $B^{} \to D^{+} \pi^{} \pi^{}$ decays in a data sample corresponding to 3.0 ${\rm fb}^{1}$ of $pp$ collision data recorded by the LHCb experiment during 2011 and 2012. A modelindependent analysis of the angular moments demonstrates the presence of resonances with spins 1, 2 and 3 at high $D^{+}\pi^{}$ mass. The data are fitted with an amplitude model composed of a quasimodelindependent function to describe the $D^{+}\pi^{}$ Swave together with virtual contributions from the $D^{*}(2007)^{0}$ and $B^{*0}$ states, and components corresponding to the $D^{*}_{2}(2460)^{0}$, $D^{*}_{1}(2680)^{0}$, $D^{*}_{3}(2760)^{0}$ and $D^{*}_{2}(3000)^{0}$ resonances. The masses and widths of these resonances are determined together with the branching fractions for their production in $B^{} \to D^{+} \pi^{} \pi^{}$ decays. The $D^{+}\pi^{}$ Swave has phase motion consistent with that expected due to the presence of the $D^{*}_{0}(2400)^{0}$ state. These results constitute the first observations of the $D^{*}_{3}(2760)^{0}$ and $D^{*}_{2}(3000)^{0}$ resonances.
Results of the fit to the $B$ candidate invariant mass distribution shown with (left) linear and (right) logarithmic $y$axis scales. Contributions are as described in the legend. 
Fig1a.pdf [31 KiB] HiDef png [252 KiB] Thumbnail [84 KiB] *.C file 

Fig1b.pdf [31 KiB] HiDef png [290 KiB] Thumbnail [97 KiB] *.C file 

Distribution of $ B ^ \rightarrow D ^+ \pi ^ \pi ^ $ candidates in the signal region over (left) the DP and (right) the SDP. 
Fig2a.pdf [626 KiB] HiDef png [341 KiB] Thumbnail [48 KiB] *.C file 

Fig2b.pdf [342 KiB] HiDef png [520 KiB] Thumbnail [69 KiB] *.C file 

The first seven unnormalised angular moments for backgroundsubtracted and efficiencycorrected data (black points) as a function of $m( D ^+ \pi ^ )$ in the range $2.0$$4.0\mathrm{\,GeV} $. The blue line shows the result of the DP fit described in Sec. 7. 
Fig3a.pdf [19 KiB] HiDef png [161 KiB] Thumbnail [57 KiB] *.C file 

Fig3b.pdf [20 KiB] HiDef png [186 KiB] Thumbnail [64 KiB] *.C file 

Fig3c.pdf [19 KiB] HiDef png [163 KiB] Thumbnail [56 KiB] *.C file 

Fig3d.pdf [21 KiB] HiDef png [183 KiB] Thumbnail [65 KiB] *.C file 

Fig3e.pdf [20 KiB] HiDef png [174 KiB] Thumbnail [61 KiB] *.C file 

Fig3f.pdf [22 KiB] HiDef png [179 KiB] Thumbnail [68 KiB] *.C file 

Fig3g.pdf [23 KiB] HiDef png [175 KiB] Thumbnail [67 KiB] *.C file 

Unnormalised angular moments $\left\langle P_7 \right\rangle$ and $\left\langle P_8 \right\rangle$ for backgroundsubtracted and efficiencycorrected data (black points) as a function of $m( D ^+ \pi ^ )$ in the range $2.0$$4.0\mathrm{\,GeV} $. The blue line shows the result of the DP fit described in Sec. 7. 
Fig4a.pdf [24 KiB] HiDef png [183 KiB] Thumbnail [70 KiB] *.C file 

Fig4b.pdf [25 KiB] HiDef png [179 KiB] Thumbnail [70 KiB] *.C file 

Zoomed views of the fourth and sixth unnormalised angular moments for backgroundsubtracted and efficiencycorrected data (black points) as a function of $m( D ^+ \pi ^ )$. The blue line shows the result of the DP fit described in Sec. 7. 
Fig5a.pdf [17 KiB] HiDef png [158 KiB] Thumbnail [55 KiB] *.C file 

Fig5b.pdf [19 KiB] HiDef png [164 KiB] Thumbnail [62 KiB] *.C file 

Signal efficiency across the SDP for $ B ^ \rightarrow D ^+ \pi ^ \pi ^ $ decays. The relative uncertainty at each point is typically $5\,\%$. 
Fig6.pdf [59 KiB] HiDef png [381 KiB] Thumbnail [150 KiB] *.C file 

Square Dalitz plot distributions for (left) combinatorial background and (right) $ B ^ \rightarrow D^{(*)+} K ^ \pi ^ $ decays. 
Fig7a.pdf [48 KiB] HiDef png [148 KiB] Thumbnail [54 KiB] *.C file 

Fig7b.pdf [50 KiB] HiDef png [154 KiB] Thumbnail [58 KiB] *.C file 

Real and imaginary parts of the Swave amplitude, shown in an Argand diagram. The knots are shown with statistical uncertainties only, connected by the cubic spline interpolation used in the fit. The leftmost point is that at the lowest value of $m( D ^+ \pi ^ )$, with mass increasing along the connected points. Each point, labelled 113, corresponds to the position of a knot in the spline, at values of $m( D ^+ \pi ^ ) = \{ 2.01, 2.10, 2.20, 2.30, 2.40, 2.50, 2.60, 2.70, 2.80, 2.90, 3.10, 4.10, 5.14 \} \,\mathrm{\,GeV} $. The points at $(0.5,0.0)$ and $(0.0,0.0)$ are fixed. The anticlockwise rotation of the phase at low $m( D ^+ \pi ^ )$ is as expected due to the presence of the $ D ^{*}_{0}(2400)^{0}$ resonance. 
Fig8.pdf [55 KiB] HiDef png [106 KiB] Thumbnail [27 KiB] *.C file 

Differences between the SDP distribution of the data and fit model, in terms of the normalised residual in each bin. No bin lies outside the $z$axis limits. 
Fig9.pdf [53 KiB] HiDef png [167 KiB] Thumbnail [64 KiB] *.C file 

Projections of the data and amplitude fit onto (top) $m( D ^+ \pi ^ )_{\rm min}$, (middle) $m( D ^+ \pi ^ )_{\rm max}$ and (bottom) $m(\pi ^ \pi ^ )$, with the same projections shown (right) with a logarithmic $y$axis scale. Components are described in the legend. 
Fig10a.pdf [57 KiB] HiDef png [227 KiB] Thumbnail [67 KiB] *.C file 

Fig10b.pdf [57 KiB] HiDef png [366 KiB] Thumbnail [96 KiB] *.C file 

Fig10c.pdf [59 KiB] HiDef png [308 KiB] Thumbnail [80 KiB] *.C file 

Fig10d.pdf [58 KiB] HiDef png [427 KiB] Thumbnail [106 KiB] *.C file 

Fig10e.pdf [56 KiB] HiDef png [277 KiB] Thumbnail [73 KiB] *.C file 

Fig10f.pdf [54 KiB] HiDef png [346 KiB] Thumbnail [89 KiB] *.C file 

Fig10g.pdf [12 KiB] HiDef png [79 KiB] Thumbnail [24 KiB] *.C file 

Projections of the data and amplitude fit onto (left) $m( D ^+ \pi ^ )$ and (right) the cosine of the helicity angle for the $ D ^+ \pi ^ $ system in (top to bottom) the low mass threshold region, the $ D ^{*}_{2}(2460)^{0}$ region, the $ D ^{*}_{1}(2680)^{0}$  $ D ^{*}_{3}(2760)^{0}$ region and the $ D ^{*}_{2}(3000)^{0}$ region. Components are as shown in Fig. 10. 
Fig11a.pdf [52 KiB] HiDef png [200 KiB] Thumbnail [67 KiB] *.C file 

Fig11b.pdf [66 KiB] HiDef png [210 KiB] Thumbnail [69 KiB] *.C file 

Fig11c.pdf [55 KiB] HiDef png [228 KiB] Thumbnail [72 KiB] *.C file 

Fig11d.pdf [69 KiB] HiDef png [241 KiB] Thumbnail [72 KiB] *.C file 

Fig11e.pdf [56 KiB] HiDef png [259 KiB] Thumbnail [85 KiB] *.C file 

Fig11f.pdf [71 KiB] HiDef png [325 KiB] Thumbnail [92 KiB] *.C file 

Fig11g.pdf [55 KiB] HiDef png [271 KiB] Thumbnail [79 KiB] *.C file 

Fig11h.pdf [70 KiB] HiDef png [309 KiB] Thumbnail [83 KiB] *.C file 

Animated gif made out of all figures. 
PAPER2016026.gif Thumbnail 
Measured properties of neutral excited charm states. World averages are given for the 1P resonances (top part), while all measurements are listed for the heavier states (bottom part). Where two uncertainties are given, the first is statistical and second systematic; where a third is given, it is due to model uncertainty. The uncertainties on the averages for the $ D ^{*}_{0}(2400)^{0}$ mass and the $D_1(2420)^0$ and $ D ^{*}_{2}(2460)^{0}$ masses and widths are inflated by scale factors to account for inconsistencies between measurements. The quoted $ D ^{*}_{2}(2460)^{0}$ averages do not include the recent result from Ref. [12]. 
Table_1.pdf [53 KiB] HiDef png [99 KiB] Thumbnail [16 KiB] tex code 

Yields of the various components in the fit to $ B ^ \rightarrow D ^+ \pi ^ \pi ^ $ candidate invariant mass distribution. Note that the yields in the signal region are scaled from the full mass range. 
Table_2.pdf [46 KiB] HiDef png [48 KiB] Thumbnail [8 KiB] tex code 

Signal contributions to the fit model, where parameters and uncertainties are taken from Ref. [19]. States labelled with subscript $v$ are virtual contributions. The model "MIPW" refers to the quasimodelindependent partial wave approach. 
Table_3.pdf [55 KiB] HiDef png [65 KiB] Thumbnail [9 KiB] tex code 

Masses and widths determined in the fit to data, with statistical uncertainties only. 
Table_4.pdf [44 KiB] HiDef png [80 KiB] Thumbnail [12 KiB] tex code 

Complex coefficients and fit fractions determined from the Dalitz plot fit. Uncertainties are statistical only. 
Table_5.pdf [56 KiB] HiDef png [76 KiB] Thumbnail [10 KiB] tex code 

Breakdown of experimental systematic uncertainties on the fit fractions (%) and masses and widths $(\mathrm{MeV} )$. 
Table_6.pdf [62 KiB] HiDef png [132 KiB] Thumbnail [21 KiB] tex code 

Breakdown of model uncertainties on the fit fractions (%) and masses and widths $(\mathrm{MeV} )$. 
Table_7.pdf [61 KiB] HiDef png [119 KiB] Thumbnail [17 KiB] tex code 

Results for the complex amplitudes. The three quoted errors are statistical, experimental systematic and model uncertainties. 
Table_8.pdf [55 KiB] HiDef png [168 KiB] Thumbnail [26 KiB] tex code 

Results for the $ D ^+ \pi ^ $ Swave amplitude at the spline knots. The three quoted errors are statistical, experimental systematic and model uncertainties. 
Table_9.pdf [47 KiB] HiDef png [254 KiB] Thumbnail [42 KiB] tex code 

Results for the fit fractions. The three quoted errors are statistical, experimental systematic and model uncertainties. 
Table_10.pdf [52 KiB] HiDef png [132 KiB] Thumbnail [20 KiB] tex code 

Results for the product branching fractions ${\cal B}( B ^ \rightarrow R\pi ^ ) \times {\cal B}(R \rightarrow D ^+ \pi ^ )$. The four quoted errors are statistical, experimental systematic, model and inclusive branching fraction uncertainties. 
Table_11.pdf [52 KiB] HiDef png [122 KiB] Thumbnail [19 KiB] tex code 

Interference fit fractions (%) and statistical uncertainties. The amplitudes are: ($A_0$) $ D ^{*}_{v}(2007)^{0}$ , ($A_1$) $ D ^+ \pi ^ $ Swave , ($A_2$) $ D ^{*}_{2}(2460)^{0}$ , ($A_3$) $ D ^{*}_{1}(2680)^{0}$ , ($A_4$) $B^{*0}_v$, ($A_5$) $ D ^{*}_{3}(2760)^{0}$ , ($A_6$) $ D ^{*}_{2}(3000)^{0}$ . The diagonal elements are the same as the conventional fit fractions. 
Table_12.pdf [33 KiB] HiDef png [37 KiB] Thumbnail [6 KiB] tex code 

(Top) Experimental and (bottom) model systematic uncertainties on the interference fit fractions (%). The amplitudes are: ($A_0$) $ D ^{*}_{v}(2007)^{0}$ , ($A_1$) $ D ^+ \pi ^ $ Swave , ($A_2$) $ D ^{*}_{2}(2460)^{0}$ , ($A_3$) $ D ^{*}_{1}(2680)^{0}$ , ($A_4$) $B^{*0}_{v}$, ($A_5$) $ D ^{*}_{3}(2760)^{0}$ , ($A_6$) $ D ^{*}_{2}(3000)^{0}$ . The diagonal elements are the same as the conventional fit fractions. 
Table_13.pdf [27 KiB] HiDef png [46 KiB] Thumbnail [7 KiB] tex code 
Created on 09 December 2018.Citation count from INSPIRE on 09 December 2018.