Improved Event Mixing for Resonance Yield Extraction †

Resonance particles, such as the K*(892) meson, are reconstructed from the invariant mass (Minv) distribution of possible particle pairs. To extract the yield with the highest precision, the combinatorial background must be determined as precisely as possible. An event-mixed Minv distribution is often used to describe the combinatorial background. However, this distribution will not contain the mini-jet-like structures present inside an event. This analysis introduces a new re-weighing scheme, where two Minv distributions of like-sign particles in the same-event and in mixed-events are used to correct the mixed-event background estimate for the mini-jet-like structure. Using PYTHIA 8.2 generated proton-proton collisions at √ sNN = 5.02 TeV, it is shown that the new method can be used to more accurately describe the combinatorial background.


Introduction
Resonance particles, such as the K*(892) meson, are produced as intermediary states in high-energy hadronic collisions.It is not possible to directly detect resonance particles due to their short life-times.Instead, the resonance particle has to be reconstructed by measuring the invariant mass (M inv ) of its decay constituents.As it is not possible to distinguish whether a single track has decayed from a resonance, the M inv has to be calculated for every possible track pair (tracks identified as possible resonance decay daughters).The resulting M inv signal distribution will contain a peak representing the yield of the resonance particle, on top of a combinatorial background.
Event mixing is a technique commonly used to estimate the combinatorial background in the signal distribution.An M inv distribution is calculated using track pairs from different events, and is then compared to the signal distribution.The event-mixed distribution contains track pairs that are completely disjoint in time, and thus are fully uncorrelated.The event-mixed distribution is then used to describe the combinatorial background in the signal distribution.This method is often preferred over using a like-sign background estimation, as the signal and background distribution will contain track pairs that have the same acceptance bias.However, this method is not able to completely describe the combinatorial background, which has been observed in K*(892) decays from A Large Ion Collider Experiment (ALICE), shown in both Figures 1 and 2.Not only is it important to understand why the event-mixed background distribution does not fully describe the background in the same-event signal distribution, but a better description of the combinatorial background will also give a more accurate estimation of the yield for a given resonance particle.This analysis will try to explain the discrepancy between the signal and the event-mixed distribution by looking at the difference in topology between the two cases.By understanding the difference in topology, a correction can be made for the event-mixed distribution.It is not presently clear which type of correction would be best suited to describe the difference in topology.This analysis will apply a correction that accounts for the difference in angular distance between track pairs for mixed and normal events.

Materials and Methods
The results are based on p-p collisions at √ s NN = 5.02 TeV, generated using PYTHIA 8.2 [3].
The analyzed data set contains 50 M events, generated using the Monash 2013 tune.The kinematic cuts used in this analysis consist of a pseudorapidity cut at |η| < 0.8, and only tracks within the transverse momentum interval 0.5 GeV < p T < 6.0 GeV are considered.The various M inv distributions that are presented throughout the report are defined in Table 1.The difference in azimuthal angle between two reconstructed tracks is defined as ∆φ = (φ 2 − φ 1 ).The ∆φ distribution is used to illustrate the different topologies of same-event and mixed-event distributions.The M inv for two tracks (with masses m 1,2 , energies E 1,2 and momentum p 1,2 ), in terms of the angle between the two track vectors θ, is written in Equation (1): The measurement ∆R = (φ 2 − φ 1 ) 2 + (η 2 − η 1 ) 2 describes the angular distance between two tracks and is used to correct the difference between same-event and mixed-event distributions.The ∆R correction is constructed by the ratio of the ∆R distributions for LSS and LSM, given a fixed M inv and p T interval.This ratio is then applied as a weight to re-weigh the USM distribution.In practice, this means that ∆R has to be stored alongside the M inv and p T for each distribution.

Results & Discussion
The ∆φ distribution between track pairs for both same-event and mixed-event distributions are shown in Figure 3.The two ∆φ distributions are normalized by their corresponding entries.
It is shown that there is a clear difference in ∆φ between mixed-events and same-events.The structure observed in the same-event distribution is due to the di-jet topology found in many p-p collisions.The nearside peak at ∆φ ≈ 0 comes from the main-leading jet; a large concentration of track pairs are focused towards the same angular direction, which is the case for a particle jet.Likewise, the small peak at ∆φ ≈ π showcases the more sub-leading jet, which has a larger spread.This structure is clearly lost in the event-mixed distribution.There is no preferred direction for the jets in respect to φ, so the jet-like topology will vanish if two or more events are mixed with perpendicular di-jets.The resulting mixed-event distribution will consist of "events" that have a more isotropic spread across the ∆φ spectrum.Thus, mixed-event and same-event distributions have completely different topologies in respect to ∆φ.The same is also true for ∆η.In accordance to Equation ( 1), the increased ∆φ value in mixed events will lead to a larger θ, which in turn will boost the value of M inv .This effect can be seen in Figure 4.   Figures 4 and 5 highlight that, for simulated particle collisions in PYTHIA, normal event mixing is not able to accurately describe the combinatorial background seen in the M inv distribution for K*(892) resonance decay, much like in real data.Applying a simple correction in respect to ∆R, a measurement that only takes into account the ∆φ and ∆η between two particle tracks, has a very clear impact on the shape for the mixed-event background estimation.The impact of the ∆R correction varies for different M inv and p T intervals.

Conclusions
Although it is clear that there is a difference in topology between mixed-events and same-events (displayed in Figure 3), and that correcting for the difference in topology gives a more accurate estimation of the combinatorial background in the signal distribution (as observed in Figures 4 and 5), it is not yet clear exactly how this correction should be performed.While ∆R re-weighing offers a very simple solution that outperforms a non-reweighed event-mixed distribution, one can note that the reweighed estimation in Figure 5 does not fully describe the combinatorial background.Adding a more complex correction scheme (taking into account both p T and ∆R) could improve the accuracy of the event-mixed background estimation, at the cost of increased complexity.

Figure 3 .
Figure 3.The ∆φ distribution for track pairs from both mixed-events and same-events.

Figure 4 .
Figure 4. Comparison between the K*(892) resonance decay and an event-mixed background distribution.The pink distribution showcases the true signal distribution.

Figure 4
Figure 4 contains the resonance decay K * → K +− π −+ together with mixed-event distributions with and without ∆R reweighing.The background only consists of combinatorial pairs and pairs from non-K* sources.Both event-mixed backgrounds are normalized to the signal distribution at M inv = 1.1-1.15GeV/c 2 .The reduced signal-background spectra for the two different backgrounds are shown in Figure 5.

Figure 5 .
Figure 5. Signal-Background, for both reweighed and non-reweighed event-mixed distributions (the same shown in Figure 4).

Table 1 .
Table defining all the M inv distributions used in this analysis.
inv of oppositely charged track pairs from a single event.LSS Like-Sign Same-Event M inv of like-charged track pairs from a single event.USM Unlike-Sign Mixed-Event M inv of oppositely charged track pairs from two different events.LSM Like-Sign Mixed-Event M inv of like-charged track pairs from two different events.