Visualisation of CP -violation effects in decay-time-dependent analyses of multibody B -meson decays

Decay-time-dependent CP -violation effects in transitions of neutral B mesons to CP -eigenstates can be visualised by oscillations in the asymmetry, as a function of decay time, between decay yields from mesons tagged as initially having B or B flavour. Such images, for example for B 0 → J/ψK 0S decays where the magnitude of the oscillation is proportional to sin(2 β ) with β being an angle of the Cabibbo– Kobayashi–Maskawa Unitarity Triangle, provide a straightforward illustration of the underlying physics. Until now there has been no comparable method to provide visualisation for the case of decays to multibody final states that are not CP - eigenstates, where interference between CP -even and -odd amplitudes provides additional physics sensitivity. A method is proposed to weight the data so that the terms of interest can be projected out and used to obtain asymmetries that visualise the relevant effects. Application of the weighting to B 0 s decays, where effects due to non-zero width difference are not negligible, provides a novel method to observe CP violation in interference between mixing and decay without tagging the production flavour.


Introduction
The topic of CP violation is the subject of numerous investigations in contemporary particle physics.This includes many studies of decay-time-dependent CP violation in transitions of neutral B mesons to both CP -eigenstates and non-CP -eigenstates.A salient example is that of B 0 → J/ψK 0 S decays, where the asymmetry in rates for neutral B mesons, tagged as either B 0 or B 0 at time t = 0 and decaying at time t, can be written as where ∆m d is the difference in masses between the heavier and lighter mass eigenstates of the B 0 -B 0 system.In the limit that the B 0 → J/ψK 0 S transition is dominated by the so-called tree amplitude, involving the Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing matrix [1,2] elements V cb V * cs , it can be shown that S = sin(2β) and C = 0 to an excellent approximation [3,4], where β is an angle of the CKM Unitarity Triangle.Plots of this asymmetry in bins of decay-time were crucial for the BaBar and Belle experiments to demonstrate the existence of CP violation in the B-meson sector, with the final results from those experiments shown in Fig. 1 [5,6].
There are also numerous examples of B meson decays to multibody final states that are not CP -eigenstates, which are of interest to probe CP violation effects.In such cases, interference between CP -even and -odd amplitudes provides additional physics sensitivity.For example, the decay-time-dependent analysis of B 0 → J/ψK * (892), with K * (892) → K 0 S π 0 , 1 decays can be used to determine cos(2β) as well as sin(2β), hence resolving a trigonometric ambiguity in the possible solutions for 2β.Such analyses have been performed by BaBar and Belle [7,8], although the results have large uncertainties.Decay-time-dependent amplitude analyses with sensitivity to both cos (2β) and sin (2β) have also been carried out with the modes B 0 → D ( * ) h 0 with h 0 = π 0 , η or ω, D * → Dπ 0 and D → K 0 S π + π − [9], B 0 → J/ψπ + π − [10], B 0 → K + K − K 0 S [11,12], and B 0 → K 0 S π + π − [13,14], although in some of these cases the interpretation of the results in terms of CKM angles involves significant theoretical uncertainties.In these publications, however, any visualisation of the data in terms of decay-time-dependent asymmetries is done by either integrating over the full phase-space of the decay or by selecting regions dominated by particular resonances. 2 This tends to dilute, and can even completely remove, the dependence of the asymmetry on the physics parameters of interest.There are also additional methods, yet to be implemented in analysis of experimental data, that 1 The notation K * (892) is used to denote an admixture of the K * (892) 0 and K * (892) 0 states, which is obtained since the final state, K 0 S π 0 , does not determine the strangeness.Similarly, the notation D ( * ) will be used to denote a neutral D ( * ) meson that can be any admixture of D ( * )0 and D ( * )0 states. 2 An exception to this appears in Ref. [7], where Fig. 12 shows the asymmetry of data weighted by the moment of C, where C is a function of angular variables and parameters in the decay B 0 → J/ψK * (892) with K * (892) → K 0 S π 0 .This weighting allows the dependence of the data on cos(2β) to be visualised, and can be considered a specific implementation of the more general weighting approach proposed in this paper.S and similar decays from the (left) BaBar [5] and (right) Belle [6] experiments.In both cases the data are separated into (top for BaBar; left for Belle) CP -odd final states such as J/ψK 0 S and (bottom for BaBar; right for Belle) CP -even final states such as J/ψK 0 L .Both experiments exploit B production via the Υ(4S) → BB process, and hence use the time variable ∆t, which is the difference in decay times of the signal and tagging B meson in the BB pair.The data are separated into cases where the tagging B meson is identified as having B 0 (blue for BaBar; red for Belle) or B 0 (red for BaBar; blue for Belle) flavour; due to entanglement of the pair the signal B meson is known to have the opposite flavour at the time of the tagging B meson decay.The BaBar data include background contributions, indicated by hatched regions, while the Belle data are background-subtracted.In both cases the asymmetry is diluted by imperfect flavour-tagging capability (the dilution is smaller for Belle since the figure includes only data with good quality tagging information).
provide sensitivity to both cos(2β) and sin(2β), such as one using B 0 → D CP π + π − decays, where D CP indicates that the neutral D meson is reconstructed in a CP eigenstate [15,16].There are further multibody B 0 meson decay modes that provide sensitivity to other CP -violation parameters, including B 0 → π + π − π 0 [17][18][19][20].Moreover, there are many B 0 s decays of interest for similar reasons, including , ϕϕ [26] and K * (892) 0 K * (892) 0 [27]. 3As such, it is of interest to develop methods that enable the dependence of the data on the underlying physical parameters to be clearly visualised, with minimal dilution of the sensitivity.
In the remainder of this paper, a method to achieve this goal is described.In Sec. 2 the underlying methodology is set out, and an illustration of the method is given in Sec. 3. The B 0 → D CP π + π − process is used as an example, but the method is general to any self-conjugate final state that is a mixture of CP -even and CP -odd contributions.This can be any multibody final state that contains a particle-antiparticle pair (e.g.π + π − , K + K − , D 0 D 0 and so on), since final states in which all particles are spin-0 CP eigenstates are themselves CP eigenstates [28].Approaches to handle experimental effects are discussed in Sec. 4. Uses of the method in the case of non-zero width difference, as is the case for B 0 s -meson decays, are described in Sec. 5.The impact of CP violation in decay is considered in Sec.6, and a summary concludes the paper in Sec. 7.

Methodology
The B 0 → D CP π + π − process [15,16] is used to illustrate the method.If a neutral B meson is tagged as B 0 at time t = 0, and evolves at a later time t into an admixture of B 0 and B 0 that decays to a position (m 2 + , m 2 − ) ≡ (m 2 (Dπ + ), m 2 (Dπ − )) in the D CP π + π − Dalitz plot, the corresponding amplitude is given by: Here, ∆m = m B 0 H − m B 0 L is the mass difference between the two physical eigenstates of the B 0 -B 0 system, 4 and it has been assumed that there is no CP violation in mixing so that, in terms of the usual mixing parameters p and q of the B 0 -B 0 system, q p has unit magnitude (see, for example, Refs.[29,30]).The amplitudes A(m 2 + , m 2 − ) and A(m 2 + , m 2 − ) are those for B 0 → D 0 π + π − and B 0 → D 0 π + π − decays, respectively, to the position (m 2 + , m 2 − ) in the Dalitz plot.These will be denoted subsequently by A f and A f , respectively.When A f and A f each include contributions with only one set of CKM matrix elements, then there is no

CP violation in decay,
For the B 0 → D CP π + π − case, neglecting contributions from doubly-Cabibbo-suppressed B decay amplitudes, i.e. from the decay B 0 → D 0 π + π − and the charge-conjugate process, ϕ mix+dec = −2β.In Eq. ( 2) the weak phase factors in the amplitudes are absorbed into the e −i2β term.In addition, the width difference in the B 0 -B 0 system is taken to be negligible, i.e. ∆Γ ≈ 0. A similar expression to Eq. (2), with A ↔ A and −2β ↔ +2β, holds for a neutral B meson tagged as B 0 at time t = 0.The symmetry in the Dalitz plot around the line m 2 + = m 2 − means that for decay to any CP eigenstate A(x, y) = ±A(y, x); the convention that the positive (negative) sign is used for CP -even (CP -odd) eigenstates is used in this paper.Thus, similar expressions can be used for any B 0 -meson decay to a self-conjugate final state, with the amplitudes expressed as a function of appropriately symmetric co-ordinates.
Returning to the specific B 0 → D CP π + π − case, Eq. (2) assumes that the neutral D meson is reconstructed in a CP -even eigenstate; for CP -odd D final states, the + sign between the two terms of Eq. ( 2) should be replaced by a − sign.Effects due to charm mixing and CP violation are neglected.Squaring Eq. ( 2), and the corresponding equation for a neutral B meson tagged as B 0 at time t = 0, expressions for the time-dependent decay rates are obtained: The normalisation is omitted as it is not relevant for the discussion here.It is assumed that t ranges over non-negative values, as is the case for B mesons produced in LHC collisions.For e + e − → Υ(4S) → BB production both positive and negative values of t are possible, as the tagging information comes from the other decay in the BB pair; the exponential factor in this case should instead be e −|t|/τ B 0 .The coefficient of the sinusoidal term of Eq. ( 3) is given by In case A * f A f has no imaginary component, as would be true if the entire D CP π + π − Dalitz plot were dominated by a single CP -eigenstate, there is no sensitivity to cos (2β).It is the interference in the final-state phase-space between different resonances with different CP compositions -for example, interference between S-and P-wave π + π − componentsthat generates the imaginary part of A * f A f and hence allows cos(2β) to be measured.Without loss of generality, the amplitudes A f and A f can be decomposed into CP -even (A + f ) and CP -odd (A − f ) parts, Note that here the CP eigenvalues are those of the decaying neutral B meson rather than of the final state, though in the absence of CP violation in decay, as assumed here, the two are equivalent.The impact of non-negligible CP violation in decay will be considered in Sec. 6.The previously mentioned phase convention implies that the B meson CP eigenstates are given by , which is used in Eq. ( 6) but does not affect the final results.
In the D CP π + π − final state with a CP -even D meson, since the π + π − system must also be CP -even, the overall CP eigenvalue is equal to (−1) l where l is the angular momentum between the D meson and the π + π − system.By angular momentum conservation in the B-meson decay, l is also the angular momentum in the π + π − system.Hence, in the absence of CP violation in decay, A + f (A − f ) corresponds to amplitudes where the π + π − system is in an even (odd) partial wave.This means that A + f (A − f ) is symmetric (antisymmetric) about the Dalitz-plot symmetry line m 2 + = m 2 − .Since any amplitude can be decomposed into a sum of partial waves, with the sum potentially running to infinity, this formalism includes not only amplitudes corresponding to ππ resonances, but also those corresponding to Dπ resonances.The interference term between the B 0 and B 0 decay amplitudes can therefore be written as and thus ) is symmetric (antisymmetric) about the Dalitz-plot symmetry line (i.e.under interchange of m 2 + and m 2 − ).This is illustrated in Fig. 2 using the model for B 0 → D CP π + π − decays described in the next section.Thus, from Eq. ( 5) it can be seen that a simple integration over the Dalitz plot will retain sensitivity to sin(2β), through the Re(A * f A f ) term, though diluted by cancellation of CP -even and CP -odd contributions.The sensitivity to cos(2β) will be completely cancelled out.
Weighting functions, which can be applied to the data to enable visualisation of the sensitivity to sin(2β) and cos (2β), can be proposed with the following logic.In order to visualise sin(2β) (cos(2β)), the function should be symmetric (antisymmetric) about the symmetry line, so that the relevant terms are retained when integrating over the Dalitz plot.The function should also carry the same sign as to avoid cancellation of Dalitz-plot regions where these terms have opposite signs.Finally, in the limit of a pure CP -eigenstate amplitude, the weighting function should be constant across the Dalitz plot, as in this case any variation in the weighting would dilute the asymmetry.Functions that satisfy these criteria, for the visualisation of the sensitivity to sin(2β) and cos (2β), are respectively.The normalisation is chosen such that these functions are each in the range [−1, +1], but this does not affect the outcome.The decay-time asymmetry is defined in terms of the time-dependent decay rates of Eqs. ( 3) and (4), When forming weighted decay-time asymmetry plots, each Γ term is integrated over the Dalitz plot.(It is also integrated over the range of each decay-time bin, but that is not with the Dalitz-plot model for B 0 → D CP π + π − decays from Ref. [31], illustrating the symmetry (antisymmetry . The z-axis scale is arbitrary, but common for the four plots. relevant to the discussion here.)The impact of the weights on all terms of Eqs. ( 3) and (4) must therefore be considered.It can be seen that is symmetric about the Dalitz-plot symmetry line, while is antisymmetric.Therefore, at least in the case of an antisymmetric weighting function, it is better not to weight the expressions in the denominator of the asymmetry to avoid a potential "divide by zero" problem.Weighted decay-time asymmetries are hence defined as where w Re/Im is one of the weighting functions from Eq. ( 11) and dΩ represents an element in the Dalitz-plot phase-space.For the first weighting function, the weighted asymmetry reduces to where while for the second weighting function instead where and Hence, A w−Re CP (t) and A w−Im CP (t) can be used to visualise the determination of sin(2β) and cos (2β), respectively, both appearing as sine waves.The asymmetry A w−Im CP (t) will also show a cosine oscillation, with magnitude independent of weak phase factors.The magnitudes of the oscillations are modulated by hadronic factors, but these constants are calculable assuming that the Dalitz-plot model is known.In the pure CP -eigenstate limit, the hadronic factor in Eq. ( 17) is equal to −1, while those in Eqs.(19) and (20) are equal to zero.

Illustration of the method
To illustrate the method, it is necessary to specify the amplitudes A f and A f .As previously, the B 0 → D CP π + π − decay is taken as an example, with amplitude model taken from the analysis by the LHCb collaboration [31], and implemented in the Laura ++ software package [32].Specifically, the results of the isobar model from Ref. [31] are used.The largest components in the model are Dρ 0 , D(π + π − ) S wave (including Df 0 (500)) and Df 2 (1270) together with D * 2 (2460) − π + , D * 0 (2400) − π + and a (Dπ − ) P wave π + term.The D * (2010) − π + component is vetoed with the requirement m(Dπ ± ) > 2.1 GeV/c 2 , since it is too narrow to contribute significant interference effects that would enhance the sensitivity to sin(2β) and cos (2β).
The hadronic factors in Eqs. ( 17), (19) and ( 20) are evaluated numerically with this model, and are found to be Distributions of A w−Re CP (t) and A w−Im CP (t), as defined in Eqs. ( 16) and ( 18), obtained with this model are shown in Fig. 3. Here, sin(2β) and cos(2β) are set at their world average values (corresponding to β = 22.2 • [33], so 0.700 and 0.714, respectively).These correspond to the distributions that would be obtained in the limit of infinite statistics with no experimental effects such as imperfect flavour tagging.
Samples of pseudoexperiments are generated in order to approximate what the weighted decay-time asymmetries might look like in an experiment with currently plausible precision.The B 0 → Dπ + π − yield obtained by the LHCb collaboration in Ref. [31], of ∼ 10 000, is scaled by 0.14 to account for the difference in branching fractions between the D → K + π − channel used in that analysis and the combined D → K + K − and D → π + π − decays that are expected to provide the largest D CP samples.Further scaling factors, indicated in parentheses, are applied to account for: the difference in sample size available to the LHCb collaboration in Ref. [31] to the total recorded in the LHC Run 1 and 2 operation periods (5); the facts that in Ref. [31] only a subset of possible trigger lines was used (2) and that very tight selection requirements were imposed (1.5); the typical LHCb flavour tagging power on similar decays (0.08) [34][35][36][37][38]. Rounding up optimistically, an estimate of 2000 is  [31], and with sample size corresponding to the equivalent perfectly tagged yield expected to be available in the LHCb data sample.The statistical uncertainties are evaluated with a bootstrap method [39].Results of the fits described in the text are also shown.
obtained for the equivalent yield of perfectly tagged B 0 → D CP π + π − decays available in the Run 1 and 2 LHCb data sample.
Neglecting other experimental effects, such as background, efficiency variations across the Dalitz plot or with decay time, and resolution, samples of this size are generated and then weighted with the functions of Eq. (11).Summing over the range of each decay-time bin and the whole Dalitz plot, the weighted decay-time asymmetries shown in Fig. 4 are constructed.The unweighted decay-time asymmetry is also shown for comparison.The utility of the weighting functions in visualising the physics effects of interest is evident.
It may be noted that, once these weighted decay-time asymmetries have been formed, values of sin(2β) and cos(2β) can be obtained by simple one-dimensional fits to them, assuming that the hadronic factors are fixed.Such fits would not be expected to give values as precise as can be obtained from the full unbinned decay-time-dependent amplitude analysis, but could be used to obtain a cross-check of the result.The values of sin(2β) and cos(2β) obtained by fitting the weighted asymmetries, shown in Fig. 4, are sin(2β) = 0.762 ± 0.064 and cos(2β) = 0.749 ± 0.068; for comparison the values obtained from the full decay-time-dependent amplitude fit to these samples are sin(2β) = 0.742 ± 0.045 and cos(2β) = 0.735 ± 0.045.
Alternatively, if the values of sin(2β) and cos(2β) are fixed, the weighted decay-time asymmetries can be used to check that the values of the hadronic factors in data are consistent with those obtained from the model.In either case, the method provides a useful new way to visualise the effects being measured, and to check the consistency of the data with the fit result.

Experimental effects
The decay-time distributions measured in experiment are not given only by the expressions of Eqs. ( 3) and ( 4), but are modified by contributions from background processes, asymmetries in the B 0 -B 0 production rate or detection probability, efficiency variations across the decay phase-space or with decay time, smearing of the distributions due to resolution, and imperfect flavour tagging.For most of these it is possible to apply corrections to the data, as discussed below, so that the asymmetry curve still represents that of the underlying physics.This is, however, not essential -the asymmetries measured in data without further corrections would still be useful to visualise the sensitivity to sin(2β) and cos(2β), although the impact of the experimental effects would need to be assessed on a case-by-case basis.
Background can be subtracted statistically using one of a number of methods.The sPlot method and its variants [40,41], based on the application of a signal weight for each candidate in the selected sample, provide an attractive option for many decays.Weights can also be used to correct for efficiency variation across the phase-space, as has been done in numerous Dalitz-plot analyses (see, e.g., Ref. [42]), and to correct for detection asymmetries, assuming that the size of these effects is known from appropriately calibrated simulation.Variation with decay-time of the efficiency will not affect the asymmetry as long as the bin size is chosen such that variations within each bin are small; if this is not the case, corrections can be applied.Similarly, effects due to imperfect flavour tagging can be corrected for either on a candidate-by-candidate basis or collectively, as the size of the required correction will be known from data control samples.An example of a decay-time asymmetry being visualised with candidate-by-candidate flavour-tagging corrections applied can be found in Ref. [43].A non-zero asymmetry in the production rates of B 0 and B 0 mesons will introduce a decay-time-independent offset, that can also be corrected for assuming that the size of the effect is known.
The only remaining effects to consider are those due to resolution.Non-negligible effects of resolution in the decay-time variable itself will result in dilution of the oscillations, that must be accounted for when interpreting the visualisation of the data.This is the case for the asymmetries obtained by the BaBar and Belle experiments, including those shown in Fig. 1.In the LHCb experiment the decay-time resolution is typically negligible compared to the period of B 0 -B 0 oscillations, but cannot be overlooked when considering B 0 s -B 0 s oscillations, as for example in Refs.[21,22,44].Smearing of the determination of the position in phase-space of each B decay would mean that the weight applied to each candidate will be calculated at the measured position, rather than the true position.As such, this could have non-negligible effects on the weighting procedure if the resolution is not much smaller than the scale on which the amplitudes A f and A f are varying, which for a Dalitz plot can be characterised by the width of the narrowest contributing resonance.The ability to impose kinematic constraints when determining the phase-space position will often mean this is not a concern.In cases where the effect must be accounted for, it will result in a modification of the hadronic factors of Eqs. ( 17), ( 19) and ( 20) that will be calculable for a known amplitude model and resolution function.

Non-zero width difference
To be applied to B 0 s decays, the effect of non-zero decay width difference must be accounted for.This results in modifications to Eqs. ( 3) and ( 4): in addition to replacing term is multiplied by cosh(∆Γt/2), and an additional term −2Re e i2βs A * f A f sinh(∆Γt/2) contributes to both B 0 s and B 0 s decay rates. 6Since this term appears with the same sign in both rates it does not contribute to the numerator of the asymmetry, and since the denominator is unweighted the replacements for Eqs. ( 16)-( 20) are where where and In practice, the modification of the asymmetries due to the non-zero value of ∆Γ will make little difference.The value of ∆Γ/Γ in the B 0 s -B 0 s system, 0.126 ± 0.007 [33] means that after five B 0 s lifetimes the values of (cosh(∆Γt/2), sinh(∆Γt/2)) are (1.05,0.32) -a relatively modest change from the values (1, 0) at t = 0.At very large numbers of lifetimes the extra term can significantly enhance or suppress the oscillations, depending on its sign, but there will be little data available to observe this effect.Thus it has become common in experimental analyses of B 0 s decays to show oscillations folded by the period 2π/(∆m) (see, for example, Refs.[21][22][23][24]), ignoring the modulation caused by the extra term at high decay times.This approach can also be used for the weighted asymmetries.
Nonetheless, the extra terms involving ∆Γ can have a noticeable effect in the untagged decay-time distribution, characterised by the effective lifetime [45].For decay to a CP eigenstate, in the absence of CP violation in decay, the effective lifetime is the inverse of Γ ± cos(2β s )∆Γ/2 where the sign is the same as the CP eigenvalue.It is thus interesting to ask if weighting the untagged data changes the decay-time distribution.To investigate this, it is useful to expand Assuming absence of CP violation in decay, the term in Eq. ( 29) that is symmetric (antisymmetric) across the phase-space is sensitive to cos(2β s ) (sin(2β s )).This is in contrast to Eq. ( 5), where it is the symmetric (antisymmetric) term that is sensitive to the sine (cosine) of the relevant weak phase ϕ mix+dec .The unweighted, untagged decay-time distribution is the sum of the B 0 s and B 0 s decay rates, integrated across the phase space, A Taylor expansion in t can be used to identify the coefficient of the first-order term as the inverse of the effective lifetime, where in the last line only the term that is symmetric across the phase space, and hence does not vanish on integration, is retained.Here, and in similar Taylor expansions below, terms of O(t 2 ) and higher within each set of parentheses are omitted.In the case of a pure CP eigenstate, Eq. ( 31) gives the same effective lifetime as discussed above.For a final state that is not a pure CP eigenstate, the deviation of the effective lifetime from Γ will be diluted by the hadronic factor in the last line of Eq. ( 31), which is bound in the range [−1, +1]. 7onsider now the untagged decay-time distribution obtained after weighting with the w Re function of Eq. ( 11), integrating across the phase-space, and Taylor-expanding in t, In the case of a pure CP eigenstate this again gives the same result,8 but in general a modification of the unweighted effective lifetime can be anticipated.This modification can be significant, since the w Re function can be negative as well as positive.Consequently, the coefficient of the cosh ∆Γt 2 term in the first line of Eq. ( 32) can also be either positive or negative, and can also be equal to zero in the limit that the final state is an equal mixture of CP -even and CP -odd (i.e. has net-CP of zero).In this limit, the division in the last line of Eq. ( 32) cannot be made and the distribution cannot be approximated by an exponential function.Nonetheless, the dependence of the w Re -weighted decay-time distribution of Eq. ( 32) on the weak phase factor is only through cos(2β s ), as in Eq. (31).
Turning now to the untagged decay-time distribution obtained after weighting with the w Im function of Eq. (11), integrating across the phase-space, and Taylor-expanding in t, In the last line, the hadronic factor has been written in a way that makes it independent of the normalisation of A f and A f , as done previously in Eqs.(31) and (32).Due to the cancellation of the cosh(∆Γ/2) term, the lowest-order term in the Taylor expansion in Eq. ( 33) is linear in t, rather than constant, and thus this expression cannot be approximated by a single exponential function. 9The magnitude of Γ w−Im untagged (t) is proportional to sin(2β s ), and thus corresponds to CP violation in the interference between mixing and decay in the B 0 s system.Although the w Im function can be both positive and negative, all hadronic factors in Eq. ( 33) are positive, and since ∆Γ is also known to be positive the sign of Γ w−Im untagged (t) depends only on the sign of sin(2β s ).Effectively, the weighted untagged decay-time distribution of Eq. ( 33) allows visualisation of the emergence, as a function of decay time, of an asymmetry across the Dalitz-plot.It therefore shares some similarity with methods of observing CP violation in decay through (decay-time-integrated) Dalitz-plot asymmetries [49,50], but is novel in that it allows CP violation in the interference between mixing and decay to be observed.As this method requires neither flavour tagging nor ability to resolve the fast B 0 s -B 0 s oscillations, it may enable some experiments to achieve sensitivity to sin(2β s ) that would otherwise be impossible. 10 Normalisation has been omitted in Eqs. ( 30)-( 33) for simplicity, even though the argument made in Sec. 2 that it is irrelevant due to cancellation in the asymmetries no longer holds.For a useful visualisation of data, with fit curve superimposed, the key issue is that the same normalisation should be used for both.In practice it is likely to be convenient to normalise Eq. ( 30) to unity, i.e. to treat it as a probability density function, and to use the same normalisation factor also in Eqs.(32) and (33).Since both of these weighted untagged decay-time distributions could, in principle, be equal to zero, they cannot in general be normalised individually.
To illustrate the concept of this weighted untagged decay-time distribution, the model for B 0 → D CP π + π − decays described in Sec. 3 is modified to introduce a decay-width difference.The value of the decay-width difference is taken to be that of the B 0 s -B 0 s system, specifically ∆Γ = 0.083 ps −1 [33], but all other features of the model remain the same (e.g. the mass, lifetime, mass difference and weak phase factor are those of the B 0 -B 0 system).This has the advantage that the larger magnitude of sin(2β) compared to sin(2β s ) provides a more easily visible effect in Eq. (33).The hadronic factors of Eqs. ( 31)-( 33) 9 Clearly the lowest order Taylor expansion is only a good approximation at values of t satisfying Γt ≪ 1.
It is nonetheless useful to illustrate the differences between the unweighted and weighted untagged decay-time distributions. 10The potential to determine sin(2β s ) from an untagged analysis of B 0 s → J/ψϕ decays was noted in Ref. [51], and the method was implemented by the D0 collaboration in Refs.[52,53], but the use of weighting to visualise the effect has not been described previously in the literature. [ps] [ps] t  [31], with a non-zero width difference introduced as described in the text.Curves corresponding to the expected distributions without approximations, from Eqs. (30), (32) and (33), are also shown and compared to exponential distributions with different decay widths.
are evaluated with this decay model to be where the last is the same as given in Eq. ( 22).The untagged decay-rate distributions obtained from a sample generated with this model, both without weighting and with the w Re and w Im weighting schemes are shown in Fig. 5.A large sample of 100 000 events is used to illustrate the impact of the method.
Since the net CP of the B 0 → D CP π + π − model, Eq. ( 34), is close to zero, the unweighted untagged decay-time distribution is well approximated by an exponential decay with the B 0 lifetime.The hadronic factor of Eq. ( 35) is, on the other hand, significantly different from zero, causing the w Re -weighted distribution to have a different effective lifetime (as well as being flipped negative).The modest value of the hadronic factor of Eq. ( 36) allows the term corresponding to CP -violation in the interference of mixing and decay to be seen.
There are a range of B 0 s decays where this visualisation method may be useful in future.The most extensively studied to date is B 0 s → J/ψK + K − [21-24], for which it will be important to be able to see in the data the oscillations caused by a non-zero value of sin(2β s ) as the sensitivity approaches the level at which this is likely to be observed.As noted above, the w Im -weighted untagged decay-time distribution may be particularly useful for experiments where the decay-time resolution is not sufficient to resolve the fast B 0 s -B 0 s oscillations.There could also be interesting applications of the method in the study of B 0 s → DK + K − decays, where the D meson can be reconstructed in a number of different final states, which is of interest to determine the CKM angle γ [54][55][56][57].Another potentially interesting mode is B 0 s → K 0 S π + π − , which provides an alternative way to determine γ [58,59].Consideration of the impact of CP violation in decay may be important in these cases, however.

Impact of CP violation in decay
The visualisation method discussed so far is applicable for B 0 -meson decays to multibody, self-conjugate final states where there is no CP violation in decay.This is the case, to a good approximation, for a number of decay modes of interest, including B 0 → J/ψK * (892) with K * (892) → K 0 S π 0 , B 0 → D ( * ) h 0 with D → K 0 S π + π − , B 0 → J/ψπ + π − and B 0 → D CP π + π − .In these cases it may nonetheless become important, as sample sizes increase, to be able to assess the impact of non-zero CP violation effects.Moreover, it would be of interest to apply the method also to decays where large CP violation effects are either known or expected, including B 0 → K + K − K 0 S , B 0 → K 0 S π + π − and B 0 → π + π − π 0 .The presence of CP violation in decay implies that the amplitude for decay to a final state that is, say, CP -even must include a contribution from the CP -odd B state.This breaks the interpretation of Eqs. ( 9) and (10) in terms of symmetries across the phase space -it means that A + f (A − f ) is no longer symmetric (antisymmetric) about the Dalitz-plot symmetry line m 2 + = m 2 − .Thus, application of the w Re and w Im weighting functions defined in Eq. (11) will no longer lead to the cancellation of terms that simplifies the asymmetries given in Eqs. ( 16)- (20).Additionally, the presence of CP violation in decay means that weak phase factors cannot in general be separated from hadronic factors.While in Eq. ( 2) the weak phase (−2β) could be separated from the amplitudes A f and A f , the existence of CP violation in decay means that A f and A f must be composed from two contributions with different weak phases (as well as a non-trivial strong phase difference).Since the relative magnitudes of these two contributions would be expected to vary across the phase-space, due to the presence of resonances, the weak phase difference between A f and A f will also vary across the phase-space and therefore cannot in general be factored out as a constant term.
The above observations are related to the well-known facts that for decays to CP eigenstates the presence of CP violation in decay introduces a cosine term to the decaytime-dependent asymmetry, and precludes a straightforward relation of the coefficient of the sine term to a weak phase.Nonetheless, as is the case for decays to CP eigenstates, the weighted asymmetries may provide an interesting and useful visualisation method, as it will still be possible to compare the effects observed in data to those expected in a model, including when the model is the result of a fit to data.

Summary
A new method to weight data in order to visualise decay-time-dependent CP violation effects in neutral B meson decays to multibody final states, that are not CP eigenstates, is proposed.In the common case that the weak phase difference between B 0 and B 0 decay amplitudes, also accounting for the phase from B 0 -B 0 mixing, is −2β and that effects of CP violation in decay can be neglected, the method allows separate visualisation of asymmetries proportional to sin(2β) and cos (2β).The method is based on fundamental symmetries of the decay amplitudes, and is applicable to a wide range of decay channels.Weighted untagged decay-time distributions can also be formed, and provide novel potential for observation of effects due to CP violation in the interference between mixing and decay in the B 0 s system, exploiting the non-zero decay-width difference.

Figure 1 :
Figure1: Decay-time distributions and asymmetries for B 0 → J/ψK 0 S and similar decays from the (left) BaBar[5] and (right) Belle[6] experiments.In both cases the data are separated into (top for BaBar; left for Belle) CP -odd final states such as J/ψK 0 S and (bottom for BaBar; right for Belle) CP -even final states such as J/ψK 0 L .Both experiments exploit B production via the Υ(4S) → BB process, and hence use the time variable ∆t, which is the difference in decay times of the signal and tagging B meson in the BB pair.The data are separated into cases where the tagging B meson is identified as having B 0 (blue for BaBar; red for Belle) or B 0 (red for BaBar; blue for Belle) flavour; due to entanglement of the pair the signal B meson is known to have the opposite flavour at the time of the tagging B meson decay.The BaBar data include background contributions, indicated by hatched regions, while the Belle data are background-subtracted.In both cases the asymmetry is diluted by imperfect flavour-tagging capability (the dilution is smaller for Belle since the figure includes only data with good quality tagging information).

2 mFigure 2 :
Figure 2: Dalitz-plot distributions of (top left) Re(A* f A f ) for Re(A * f A f ) > 0, (top right) −Re(A * f A f ) for Re(A * f A f ) < 0, (bottom left) Im(A * f A f ) for Im(A * f A f ) > 0,and (bottom right) −Im(A *f A f ) for Im(A * f A f ) < 0, with the Dalitz-plot model for B 0 → D CP π + π − decays from Ref.[31], illustrating the symmetry (antisymmetry) of Re(A * f A f ) (Im(A * f A f )).The z-axis scale is arbitrary, but common for the four plots.

Figure 3 :
Figure 3: Decay-time dependence of A w−Re CP (t) and A w−Im CP (t) in B 0 → D CP π + π − decays with the Dalitz-plot model from Ref. [31].

Figure 4 :
Figure 4: Distributions of time-dependent asymmetry with the time-dependent decay rates in the numerator (top) unweighted, (bottom left) weighted by w Re , and (bottom right) weighted by w Im .These are obtained from pseudoexperiments generated with the Dalitz-plot model from Ref.[31], and with sample size corresponding to the equivalent perfectly tagged yield expected to be available in the LHCb data sample.The statistical uncertainties are evaluated with a bootstrap method[39].Results of the fits described in the text are also shown.

Figure 5 :
Figure5: Distributions of untagged decay-rate distributions (top) without weighting, and with (bottom left) w Re and (bottom right) w Im weights applied.These are obtained from pseudoexperiments generated with the Dalitz-plot model from Ref.[31], with a non-zero width difference introduced as described in the text.Curves corresponding to the expected distributions without approximations, from Eqs.(30),(32) and(33), are also shown and compared to exponential distributions with different decay widths.