# A Method for the Direct Absolute Measurement ofJ/ψ Decay with ψ(3686) Data Set

Authors: Fang Yan; Bo Zheng
Journal: Advances in High Energy Physics
 (2019)

---

## Abstract

To take the full advantage of theψ(3686) data set collected at e+e- collider at τ-charm energy region, a tag method is developed to directly measure the J/ψ meson decay branching fractions absolutely. The J/ψ meson decay can be measured with the J/ψ sample tagged by the two soft charged pions from the decay ψ(3686)→J/ψπ+π-. This method is illustrated by comparing the input and output branching fractions of J/ψ→γη with 106 million inclusive ψ(3686) Monte Carlo samples. The consistent result confirms the validity of the tag method.


---

## Body

## 1. Introduction

The firstly discovered charmoniumJ/ψ [1, 2], JPC=1--, is the lowest one among those which can be produced directly in e+e- annihilation. Many experiments [3–14] have been performed to study its production and decay properties. However, the summation of the measured J/ψ decay branching fraction is not more than 60% [15] without considering the interference between resonance and continuum amplitudes until now, which hampers the understanding of its properties. The precisely measured branching fractions of J/ψ decay not only provide information to understand the J/ψ properties, but also can test OZI (Okubo-Zweig-Iizuka) rule, flavor SU(3) symmetry, and perturbative QCD [16].TheJ/ψ production and decay property are studied by using the data collected from e+e- collisions at the J/ψ resonance traditionally. At this energy, the events consist mostly of e+e-→J/ψ, e+e-→l+l- (l represents e,μ), with small amounts of three-flavor continuum and other processes such as γV, where V represents vector meson. Due to some unavoidable influences, such as interference effect and undistinguishable backgrounds, some J/ψ decay channels can not be studied by using this kind of data sample. An energy scan experiment performed around the J/ψ resonance can solve the problem induced by interference effect, while it is not effective for the final states with undistinguishable backgrounds from other processes. Usually, a large data set will be taken at ψ(3686) resonance for τ-charm factory. Taking the Beijing Spectrometer III (BESIII) [17] at the Beijing Electron Positron Collider II (BEPCII) [17] as an example, a goal of 3.2 billion ψ(3686) events is set to be taken before its closure. Considering the large branching fraction of ψ(3686)→J/ψπ+π-, (34.67±0.30)% [15], this sample can naturally be used to study the J/ψ decay.In this paper, we propose a method to construct aJ/ψ sample by tagging the J/ψ meson with the two soft charged pions from ψ(3686)→J/ψπ+π-(called tag method for convenient) in the ψ(3686) data sample. The feasibility of the method is then verified by examining the input and output branching fractions of J/ψ→γη in the inclusive ψ(3686) Monte Carlo (MC) sample.


## 2. BESIII Detector and Simulation

There are four subdetectors at BESIII detector, which has been described elsewhere [17]. From the inner to the outside is Main Drift Chamber (MDC), Time of Flight (TOF), Electromagnetic Calorimeter (EMC), and Muon Counter (MUC). The information from four subdetectors is used to identify and select candidate particles. The Superconducting Magnet, between EMC and MUC, provides 1 T magnetic field.The work is performed in the framework of the BESIII Offline Software System (BOSS) [18], the GAUDI [19] based, which contains five subprojects such as framework, simulation, calibration, reconstruction, and analysis tools. Monte Carlo (MC) simulations are used to optimize the event selection and background estimation. The simulation software, the GEANT4-based, includes the geometric and material description of the BESIII detector, the detector response, running conditions, and performance. The production of ψ(3686) is simulated by the KKMC [20, 21] generator, while its decay is generated by EVTGEN [22, 23] for known decay channels with branching fractions being set to the PDG [15] values, and by LUNDCHARM [24, 25] for the remaining unknown decay. In this work, 106 million inclusive ψ(3686) events generated by data production group of BESIII are used. In addition, 100 000 exclusive ψ(3686)→J/ψπ+π- with J/ψ→γη and η→γγ events are generated with JPIPI generator [22, 23] for ψ(3686)→J/ψπ+π- and phase space generator for J/ψ→γη and η→γγ.


## 3. Introduction of the Method

Usually, theψ(3686)→J/ψπ+π- is used to study J/ψ decay experimentally by calculating the branching fraction of J/ψ→f (f denotes the studied final states) with formula(1)BJ/ψ→f=Nobsϵ×Nψ3686tot×Bψ3686→J/ψπ+π-,where B, Nobs, ϵ, and Nψ(3686)tot represent the branching fraction, the number of observed events, the detection efficiency for the whole process (ψ(3686)→J/ψπ+π-) with J/ψ→f), and the total number of ψ(3686) events. This is an indirect measurement, which replies on the input B(ψ(3686)→J/ψπ+π-). While an update measurement of B(ψ(3686)→J/ψπ+π-) is made, the B(J/ψ→f) should be updated accordingly.To solve this problem, a tag method can be employed to determine theB(J/ψ→f), which is a direct measurement method. With the dominant ψ(3686) decay channel, ψ(3686)→J/ψπ+π-, the J/ψ meson can be tagged with the two soft opposite charged pions. If the two pions are reconstructed correctly, there must be a J/ψ meson in the event, then the J/ψ decay can be studied. With this method, the branching fraction can be calculated directly by(2)BJ/ψ→f=Nobsϵ×NJ/ψtag,where NJ/ψtag is the number of tagged J/ψ mesons in ψ(3686) sample and ϵ is the detection efficiency of J/ψ→f. In the following, we take the determination of B(J/ψ→γη) as an example to illustrate and validate this method.


## 4. General Track Selection

To be accepted as a good photon candidate, a neutral electromagnetic shower in the EMC must satisfy fiducial and shower-quality requirement. The good photon candidate showers reconstructed from the barrel EMC (cos⁡θ<0.80) must have a minimum energy of 25 MeV, while those in the end caps (0.86<cos⁡θ<0.92) must have at least 50 MeV, where the energies deposited in nearby TOF counters are included. Showers in the region between the barrel and the end caps are poorly measured and excluded. To eliminate showers from charged particle, a photon candidate must be separated by at least 10∘ from any charged track. The time of EMC cluster (TEMC) requirements is used to suppress electronic noise, which is 0≤TEMC≤14 (in unit of 50 ns).Charged tracks in BESIII detector are reconstructed from MDC hits. Each charged track is required to satisfyVz<10 cm, Rxy<1 cm and cos⁡θ<0.93, where Vz and Rxy are the closest approach to the beam axis in z direction and x-y plane and θ is the polar angle.


## 5. Reconstruction of Tag Side

The two charged tracks in the tag side are required to satisfyp→<0.45 GeV/c and cos⁡θ+-<0.95, where p→ is the momentum of the candidate track and θ+- is the angle between positive and negative charged tracks. The candidate tracks are assumed to be pions. The recoil mass of π+π-, Mπ+π-rec, is calculated by(3)Mπ+π-rec=ptot-pπ+-pπ-2,where ptot, pπ+, and pπ- are the four momenta of e+e-, π+, and π-. To reject the obvious background events, we only keep the events satisfying 3.03<Mπ+π-rec<3.17 GeV/c2. The Mπ+π-rec spectra of candidate events are shown in Figure 1, which shows a clear J/ψ peak over a smooth background.Figure 1
TheMπ+π-rec spectra of candidate events from inclusive ψ(3686) MC sample (the dots with error bars). The solid line shows the best fit to the spectra and the dashed line shows the background. The J/ψ signal region is shown between the two arrows.We fit the distribution ofMπ+π-rec with the signal shape obtained from exclusive decay of ψ(3686)→J/ψπ+π-,J/ψ→μ+μ-, which is shown in Figure 2, and a polynomial background. In Figure 2, the dots with error bar show the distribution of Mπ+π-rec from pure ψ(3686)→J/ψπ+π-,J/ψ→μ+μ- MC events, and the solid line shows the extracted shape. Fitting the Mπ+π-rec spectra of candidate events from 106 million inclusive ψ(3686) MC samples with maximum likelihood method, the number of tagged J/ψ mesons can be obtained. The fitting results are shown in Figure 1, where the dots with error bar represent the numbers of events, the solid curve shows the total fitting result, and the dashed line shows the background. The fitting results given (2.2515±0.0005)×107 tagged J/ψ mesons. Considering the input B(ψ(3686)→J/ψπ+π-)=32.6%, the tag efficiency is 65.2%.Figure 2
TheMπ+π-rec spectra of the signal events from pure ψ(3686)→J/ψπ+π-,J/ψ→μ+μ- sample (the dots with error bars) and the signal shape (the solid blue line).To study theJ/ψ decay in the recoil side, the Mπ+π-rec is required to be located in J/ψ signal region, which is defined as from 3.082 to 3.112 GeV/c2 according to the resolution of Mπ+π-rec from the fitting. Integrating the signal distribution in the signal region, we obtain (2.1120±0.0005)×107J/ψ mesons.


## 6. Analysis of theJ/ψ Decay

We choose the channelJ/ψ→γη to study the validation of this method. The η meson is reconstructed with η→γγ and therefore there are three photons in the final states, one energetic radiative photon (denote by γ1) and other two relative soft photons (denote by γ2 and γ3). Exactly three photons are required in each candidate event. For the energetic radiative photon, the deposited energy in EMC is required to be greater than 1.0 GeV in J/ψ rest frame. For the other two photons, the invariant mass of them (Mγ2γ3) is required to satisfy 0.45<Mγ2γ3<0.65 GeV/c2. To exclude the background from ψ(3686)→γη′,η′→ηπ+π-, the invariant mass of γ2γ3 and two soft pions (Mγ2γ3π+π-) is required to be greater than 1.0 GeV/c2. A four-momentum constraint kinematic fit is performed to the candidate tracks and the χ2 of the kinematic fit is required to be less than 40. The survived events are treated as J/ψ→γη candidates. The Mγ2γ3 spectra are examined to determine the number of signal events.There are backgrounds from other decay channels. The MC truth information is used to study the background. Table1 shows the decay chain of background channels for J/ψ→γη,η→γγ, and the number of background events from each channel is listed in the last column. The Mγ2γ3 spectra of background events are shown in Figure 3, where the signal region is between two arrows. The figure shows that there is no peak background contribution to J/ψ→γη; the γ2γ3 spectra from background events can be described with smooth function.Table 1
The background ofJ/ψ→γη,η→γγ. N denotes the number of events passed from each decay chain.
decay chain N ψ ( 3686 ) → J / ψ π + π - , J / ψ → γ η ′ , η ′ → γ γ 51 ψ(3686)→J/ψπ+π-,J/ψ→γπ0 20 ψ(3686)→J/ψπ+π-,J/ψ→γf4,f4→π0π0 7 ψ(3686)→γFSRe+e-γFSR 1 ψ(3686)→J/ψπ+π-,J/ψ→γηπ0,η→γγ 1 ψ(3686)→J/ψπ+π-,J/ψ→γf2,f2→π0π0 1 ψ(3686)→γη′,η′→π+π-η,η→γγ 1 ψ(3686)→π+π-π0η,η→γγ 1 ψ(3686)→π0h1,h1→ρ+π-,ρ+→π+π0 1 ψ(3686)→ρ+a0-,ρ+→π+π0,a0-→π-η,η→γγ 1Figure 3
TheMγ2γ3 spectra of background events. The signal region of Mγ2γ3 is shown between the two arrows.


## 7. Results

TheMγ2γ3 distribution from signal MC sample is shown in Figure 4 in dots with error bar, and the solid curve in the figure shows the probability density function extracted from the shape. The Mγ2γ3 spectra for the candidate events from inclusive ψ(3686) sample are shown in Figure 5. We fit the spectra with the shape of Mγ2γ3 from signal MC sample to describe the signal and 1st order Chebychev polynomial to describe the background. The fitting results are shown in Figure 5, where the dots with error bars represent the number of events, the solid curve shows the total fit results, and the dashed line shows the background. The detection efficiency ϵ can be determined with signal MC sample, which is listed in the Table 2. Inserting the numbers into (2), considering the input B(η→γγ)=0.3925, we can calculate B(J/ψ→γη),(4)BJ/ψ→γη=3399±600.4149×2.1120×107×0.3925=9.9±0.2×10-4which is consistent with the input branching fraction 9.8×10-4.Table 2
The input and output results.
input output Ntag  ( 2.1120 ± 0.0005 ) × 1 0 7 Nobs  3399 ± 60 ϵ  ( 41.49 ± 0.39 ) % B(J/ψ→γη) 9.8 × 1 0 - 4 ( 9.9 ± 0.2 ) × 1 0 - 4Figure 4
TheMγ2γ3 spectra of signal events from J/ψ→γη MC sample (dots with error bar) and the shape used to fit the Mγ2γ3 (solid blue line).Figure 5
TheMγ2γ3 spectra of candidate events from 106 million inclusive ψ(3686) MC samples (dots with error bar). The solid line shows the best fitting to the spectra and the dashed line shows the background.When considering the systematic uncertainty in the measurement ofJ/ψ→γη, three sources, which are from the total number of ψ(3686), the tracking efficiency of two soft pions, and the B(ψ(3686)→J/ψπ+π-), should be considered when the traditional method is used, but will not present in the tag method. The numbers for theses sources are 0.7% [26], 2.0% [27] in the BESIII experiment, and 0.9% [15], respectively. Therefore the total systematic uncertainty of the branching fraction measurement of this channel by using the tag method will be less than that of the traditional one.


## 8. Conclusion

In conclusion, we have developed a method to study theJ/ψ decay with the ψ(3686) data set. The two soft opposite charged pions from ψ(3686)→J/ψπ+π- are used to tag the J/ψ meson. With the tagged J/ψ mesons from 106 million inclusive ψ(3686) MC samples, the output branching fraction of J/ψ→γη is in good agreement with the input one, which gives a solid validation of this method. By employing this method, the J/ψ decay branching fractions can be precisely measured absolutely with the large data set accumulated with BESIII and other detectors which will be run in charm energy region in the future [28] for the channels affected by interference effect or with undistinguishable background.



---

*Source: 2567070-2019-01-17.xml*