Isolating the chiral magnetic effect from backgrounds by pair invariant mass

TL;DR

This paper proposes a novel method to isolate the chiral magnetic effect in heavy ion collisions by analyzing the pair invariant mass dependence of charge correlations, aiming to reduce background interference.

Contribution

It introduces a differential measurement approach based on pair invariant mass to distinguish the CME signal from background effects in experimental data.

Findings

Model studies show the method's feasibility.

High invariant mass regions reduce background influence.

The approach enhances CME detection sensitivity.

Abstract

Topological gluon configurations in quantum chromodynamics induce quark chirality imbalance in local domains, which can result in the chiral magnetic effect (CME)--an electric charge separation along a strong magnetic field. Experimental searches for the CME in relativistic heavy ion collisions via the charge-dependent azimuthal correlator ($\Delta\gamma$) suffer from large backgrounds arising from particle correlations (e.g. due to resonance decays) coupled with the elliptic anisotropy. We propose differential measurements of the $\Delta\gamma$ as a function of the pair invariant mass ($m_{\rm inv}$), by restricting to high $m_{\rm inv}$ thus relatively background free, and by studying the $m_{\rm inv}$ dependence to separate the possible CME signal from backgrounds. We demonstrate by model studies the feasibility and effectiveness of such measurements for the CME search.

Figures

Equations

$$\gamma \equiv ⟨\cos (\alpha +\beta -2\psi )⟩,$$

$$\frac{d^{2}N}{d{p}_{T}d\eta d\varphi }=\frac{d^{2}N}{2\pi d{p}_{T}d\eta }\left(1+2v_2\cos 2\varphi +2a_1\sin \varphi \right),$$

$$\Delta \gamma =\frac{N_{\mathrm{SS}}({\gamma }_{\mathrm{CME}}+{\gamma }_{m.c.})+N_{\rho }{\gamma }_{\rho }}{N_{\mathrm{SS}}+N_{\rho }}-(-{\gamma }_{\mathrm{CME}}+{\gamma }_{m.c.})=r({\gamma }_{\rho }-{\gamma }_{m.c.})+(1-r/2)\Delta {\gamma }_{\mathrm{CME}}.$$

$$\Delta \gamma (m_{\mathrm{inv}})=r(m_{\mathrm{inv}})R(m_{\mathrm{inv}})+\Delta {\gamma }_{\mathrm{CME}}(m_{\mathrm{inv}}).$$

$$R(m_{\mathrm{inv}})\equiv ⟨f(m_{\mathrm{inv}})v_2(m_{\mathrm{inv}})⟩-{\gamma }_{m.c.}$$

$$q_2=\frac{1}{\sqrt{M}}\sum _{j=1}^N{e}^{in{\varphi }_j},$$

Full text

Chiral magnetic effect, invariant mass, QGP, QCD

I Introduction

Vacuum fluctuations in quantum chromodynamics (QCD) can result in metastable local domains of gluon fields with non-vanishing topological charges Lee:1974ma ; Kharzeev:1998kz ; Kharzeev:1999cz . Interactions with those gluon fields can, under the approximate chiral symmetry restoration, change the overall chirality of quarks in those domains Kharzeev:2007tn ; Kharzeev:2007jp . The chirality imbalance yields an electric charge separation under a strong magnetic field, a phenomenon called the chiral magnetic effect (CME) Fukushima:2008xe ; Muller:2010jd ; Liu:2011ys . Strong magnetic fields are generated at early times by the spectator protons in relativistic heavy ion collisions, raising the possibility to detect the CME in those collisions Voloshin:2004vk ; Kharzeev:2004ey . An observation of the CME would confirm a fundamental property of QCD and is therefore of great importance Kharzeev:2015znc . CME-like phenomena are not specific only to QCD and may have been observed in condense matter physics Li:2014bha .

A commonly used variable to measure the CME-induced charge separation in heavy ion collisions is the three-point correlator Voloshin:2004vk ,

[TABLE]

where $\alpha$ and $\beta$ are the azimuthal angles of two particles and $\psi$ is that of the reaction plane (span by the beam and impact parameter directions of the colliding nuclei). Charge separation along the magnetic field ($\vec{B}$), which is perpendicular to $\psi$ on average, would yield different values of $\gamma$ for particle pairs of same-sign (SS) and opposite-sign (OS) charges: $\gamma_{\rm SS}=-1,\gamma_{\rm OS}=+1$. However, there exist background correlations unrelated to the CME Voloshin:2004vk ; Wang:2009kd ; Bzdak:2009fc ; Liao:2010nv ; Bzdak:2010fd ; Schlichting:2010qia ; Pratt:2010zn ; Petersen:2010di ; Toneev:2012zx . For example, transverse momentum conservation induces correlations among particles enhancing back-to-back pairs Bzdak:2009fc ; Liao:2010nv ; Bzdak:2010fd ; Schlichting:2010qia ; Pratt:2010zn . This background is independent of particle charges, affecting SS and OS pairs equally and cancels in the difference, $\Delta\gamma\equiv\gamma_{\rm OS}-\gamma_{\rm SS}$. Recent experimental searches have thus focused on the $\Delta\gamma$ observable Kharzeev:2015znc ; Zhao:2018ixy ; Zhao:2018skm ; the CME would yield $\Delta\gamma>0$. There are, however, also mundane physics that differ between SS and OS pairs. One such physics is resonance/cluster decays Voloshin:2004vk ; Wang:2009kd ; Bzdak:2009fc ; Liao:2010nv ; Bzdak:2010fd ; Schlichting:2010qia ; Pratt:2010zn , more significantly affecting OS pairs than SS pairs. Backgrounds arise from the coupling of elliptical anisotropy ($v_{2}$, a common phenomenon in heavy ion collisions Heinz:2013th ) of resonances/clusters and the angular correlations between their decay daughters (nonflow) Voloshin:2004vk ; Wang:2009kd ; Bzdak:2009fc ; Schlichting:2010qia . Take $\rho\rightarrow\pi^{+}\pi^{-}$ as an example. The background is $(\Delta\gamma)_{\rho}=r_{\rho}\gamma_{\rho}$, where $r_{\rho}=N_{\rho}/(N_{\pi^{+}}N_{\pi^{-}})$ is the relative abundance of $\rho$-decay pairs over all OS pairs, and $\gamma_{\rho}\equiv\langle f_{\rho}v_{2,\rho}\rangle=\langle\cos(\alpha+\beta-2\phi_{\rho})\cos 2(\phi_{\rho}-\psi)\rangle$ quantifies the $\rho$ decay angular correlations coupled with its $v_{2}$ Voloshin:2004vk ; Wang:2016iov ; Zhao:2018ixy ; Zhao:2018skm .

Experimentally, significant positive $\Delta\gamma$ values have been observed at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC) Abelev:2009ac ; Abelev:2009ad ; Abelev:2012pa ; Adamczyk:2013hsi ; Adamczyk:2014mzf . The relative background and CME contributions are under extensive debate QM17 . The recent observations of comparable $\Delta\gamma$ in small system collisions Khachatryan:2016got ; Sirunyan:2017quh ; Zhao:2017ckp ; ZhaoJieQM17 ; Zhao:2017wck ; Zhao:2018pnk , where any CME signals would average to zero Khachatryan:2016got ; Belmont:2016oqp , challenge the CME interpretation of the measured $\Delta\gamma$ in heavy ion collisions. The major difficulty in distinguishing CME from backgrounds with the $\Delta\gamma$ observable is their similar behaviors with respect to the event multiplicity ZhaoJieQM17 ; Zhao:2017wck ; Zhao:2018pnk . This is because the magnetic field strength, to which the CME is sensitive, has a similar dependence as the $v_{2}$ (backgrounds) on the event multiplicity Kharzeev:2004ey ; Bzdak:2011yy ; Deng:2012pc ; Bloczynski:2012en ; Chatterjee:2014sea . There have been various proposals and attempts to reduce or eliminate the backgrounds Ajitanand:2010rc ; Bzdak:2011np ; Adamczyk:2013kcb ; Wang:2016iov ; Wen:2016zic ; Acharya:2017fau ; Sirunyan:2017quh . The central idea is to “hold” the magnetic field fixed (in a narrow centrality) and vary the event-by-event $v_{2}$ from statistical and dynamical fluctuations Voloshin:2010ut ; Adamczyk:2013kcb ; Chatterjee:2014sea . The first attempt was carried out by STAR Adamczyk:2013kcb where a charge asymmetry observable was analyzed as a function of the observed event-by-event $v_{2}$. A linear dependence was observed, expected from background, and the intercept was extracted representing a background-suppressed signal. ALICE Acharya:2017fau divided their data in each collision centrality according to $v_{2}$ in one phase space, and found the $\Delta\gamma$ to be approximately proportional to the $v_{2}$ in the phase space of the $\Delta\gamma$ measurement, consistent with background contributions. However, as recently pointed out by two of us Wang:2016iov , those methods suppressing background may not completely eliminate it. Another way to help search for the CME is to compare isobaric collisions Voloshin:2010ut , where the magnetic fields differ and the backgrounds are expected to be the same Deng:2016knn . However, these simple expectations may not be correct because of the non-identical isobaric nuclear structures Xu:2017zcn .

A new method to search for the CME, as we demonstrate in this article, is to eliminate the resonance background contributions using particle pair invariant mass ($m_{\rm inv}$) by (i) applying a lower cut on the $m_{\rm inv}$, (ii) fitting the low $m_{\rm inv}$ region by a two-component (TC) model, and (iii) fitting the low $m_{\rm inv}$ region by a two-component model aided by event-shape engineering (TC+ESE). We illustrate our method using the AMPT (A Multi-Phase Transport) model Zhang:1999bd and a toy Monte Carlo (MC) simulation. In both, the resonance masses are sampled from Breit-Wigner distributions Agashe:2014kda ; Lin:2014tya . We use pions within pseudorapidity $|\eta|<1$ and $0.2<p_{T}<2$ GeV/$c$, except otherwise specified.

It is worthwhile to note here, since the CME is generally a low $p_{T}$ phenomenon, that large $m_{\rm inv}$ region may not be the best place to search for it. If no CME is found at large $m_{\rm inv}$, it does not necessarily mean that CME does not exist at low $m_{\rm inv}$. However, if a finite signal is observed at large $m_{\rm inv}$, where background sources diminish, it could be a strong evidence for the possible existence of the CME that would call for further investigations. It is also worthwhile to note that our low $m_{\rm inv}$ fitting procedure is not without assumptions. Although the background shape may be assessed by the ESE method, there is still the question of the CME shape as function of $m_{\rm inv}$. In principle, if the CME and the background have the same $m_{\rm inv}$ shape, then our invariant mass method would not be able to distinguish the two contributions. Theoretical inputs are much needed in this respect. Nevertheless, our method provides a potentially powerful way to extract the CME with reasonable assumptions of the $m_{\rm inv}$ dependence of the CME, which can be gauged to some extent by the fitting quality.

The rest of the paper is organized as follows. Section II presents the AMPT and toy-model studies by applying a lower $m_{\rm inv}$ cut to extract the $\Delta\gamma$ in the high $m_{\rm inv}$ region. Section III describes a two-component model fit to the $\Delta\gamma(m_{\rm inv})$. Section IV extends the fitting method further by using ESE to determine the background $\Delta\gamma(m_{\rm inv})$ shape to help extract the CME signal. Section V puts our methods into experimental context and summarizes the paper.

II High-$m_{\rm inv}$ region

II.1 A transport model study with null CME

AMPT is a parton transport model Zhang:1999bd . It consists of a fluctuating initial condition, parton elastic scatterings, quark coalescence for hadronization, and hadronic interactions. The initial condition is taken from HIJING Gyulassy:1994ew . The string melting version Lin:2001zk is used in this study. Two-body elastic parton scatterings are treated with Zhang’s Parton Cascade Zhang:1997ej , where the parton scattering cross section is set to 3 mb. After partons stop interacting, a simple quark coalescence model is applied to describe the hadronization process that converts partons into hadrons Lin:2004en . Subsequent interactions of these formed hadrons are modeled by a hadron cascade Lin:2004en . However, it is known that this version of the hadron cascade does not conserve charge, which is critical to the charge correlation study here. The hadronic scatterings, while responsible for the majority of the $v_{2}$ mass splitting, are unimportant for the main development of $v_{2}$ Li:2016flp ; Li:2016ubw , and thus may not be critical for the CME backgrounds. We thus turn off hadronic cascade in AMPT for our study here, as was done in Ref. Ma:2011uma . AMPT has been quite successful in describing variety of heavy ion data Lin:2014tya . It reproduces approximately the measured particle yields and distributions, and therefore should approximately describe those of resonances as well, which is relevant to the CME background study here. We simulate Au+Au collisions at $\sqrt{s_{{}_{\rm NN}}}=200$ GeV of various impact parameter ($b$) ranges. For simplicity we use the known reaction plane in our analysis, fixed at $\psi=0$.

Figure 1(a) shows the $m_{\rm inv}$ distribution of the excess OS over SS pion pairs ($N\equiv N_{\rm OS}-N_{\rm SS}$), with $b=6.8$-8.2 fm (corresponding to the 20-30% centrality of Au+Au collisions Abelev:2008ab , and average pion multiplicities $N_{\pi^{+}}\approx N_{\pi^{-}}\approx 210$ within $|\eta|<1$). The $\rho$ peak is evident; the lower mass peaks are from Dalitz decays of $\eta$ and $\omega$ mesons (the $K_{S}$ is kept stable in AMPT). Figure 1(b) shows the $\Delta\gamma$ as a function of $m_{\rm inv}$. The $\rho$ contribution is clearly seen in the $\rho$ mass region. Since no CME is present in AMPT, the finite $\Delta\gamma$ at $m_{\rm inv}\lesssim 2$ GeV/$c^{2}$ must be due to correlations from resonance decays, or generally, correlated pion pairs. This has been observed before Ma:2011uma ; ZhaoJieQM17 . For m_{\rm inv}$$>2 GeV/$c^{2}$ where resonance contribution to the OS over SS excess is small, the $\Delta\gamma$ value is essentially zero, as expected.

Figure 2 shows the $\Delta\gamma$ in AMPT from all pairs and $\Delta\gamma$(m_{\rm inv}$$>2 GeV/$c^{2}$) from pairs with m_{\rm inv}$$>2 GeV/$c^{2}$. The positive $\Delta\gamma$ is due to backgrounds; in $\Delta\gamma$(m_{\rm inv}$$>2 GeV/$c^{2}$) this background is essentially eliminated, and as expected the result is consistent with zero. With the $11\times 10^{6}$ AMPT events simulated for 200 GeV Au+Au collisions with $b=6.6$-8.2 fm, the inclusive $\Delta\gamma$ value is $(8.1\pm 0.1)\times 10^{-5}$, and $\Delta\gamma$(m_{\rm inv}$$>2 GeV/$c^{2}$)=$(-0.6\pm 0.8)\times 10^{-5}$. This represents a null signal with an upper limit of 20% of the inclusive $\Delta\gamma$ with 98% confidence level (CL).

II.2 A toy model study with finite CME

In light of the AMPT results, we propose to apply a lower $m_{\rm inv}$ cut in real data analysis to search for the CME. We illustrate this point further by using a toy MC with input CME signal. Our toy model generates primordial $\pi^{\pm}$, $K_{S}$, and resonances ($\rho,\eta,\omega$), and decays the $K_{S}$ and resonances (via both two- and three-body decays Agashe:2014kda ). Particle kinematics are sampled according to

[TABLE]

where $a_{1}$ is the CME signal parameter Voloshin:2004vk . The particle $dN/dy$, $p_{T}$ spectra, and $v_{2}(p_{T})$ correspond to the 40-50% centrality of Au+Au collisions; they are as same as those used in Ref. Adamczyk:2015lme except that the primordial pion $p_{T}$ spectra are parameterized here with a better agreement with data at high $p_{T}$, and we have added $K_{S}$. The pion multiplicities within $|\eta|<1$ are $N_{\pi^{+}}\approx N_{\pi^{-}}\approx 100$ and those of primordial pions are $N^{\rm prim}_{\pi^{+}}\approx N^{\rm prim}_{\pi^{-}}\approx 60$. The $K_{S}$ multiplicity is taken to be 1/5 of the measured $\rho$’s Adams:2003cc , because some $K_{S}$’s would have both their decay pions reconstructed as primary particles in experiments (such as STAR). We generate $200\times 10^{6}$ events with an input CME signal of overall strength $a_{1}=\pm 0.008$ for primordial $\pi^{\pm}$; for $K_{S}$ and resonances $a_{1}=0$. Our input CME is independent of the particle $p_{T}$. This is supported by a recent theoretical study Shi:2017cpu , where the CME is insensitive to $p_{T}$ once $p_{T}$ is above 0.2 GeV/$c$.

Figure 3(a) shows the relative OS pair excess, $r(m_{\rm inv})\equiv(N_{\rm OS}-N_{\rm SS})/N_{\rm OS}$ as a function of $m_{\rm inv}$ from the toy MC. The $K_{S}$ and $\rho$ peaks are evident. Figure 3(b) shows the $\Delta\gamma(m_{\rm inv})$; the $K_{S}$ and $\rho$ contributions are clear. The inclusive $\Delta\gamma$ from Fig. 3(b) is $(24.5\pm 0.1)\times 10^{-5}$; our input CME signal of $2a_{1}^{2}$, diluted by $(N^{\rm prim}_{\pi}/N_{\pi})^{2}$, is $4.6\times 10^{-5}$, about 20% of the inclusive $\Delta\gamma$ value. The $\Delta\gamma(m_{\rm inv})$ distribution in Fig. 3(b) has a pedestal corresponding to the input CME signal. The pedestal extends to high $m_{\rm inv}$ (not shown) where resonance backgrounds vanish. A lower $m_{\rm inv}$ cut removes backgrounds to $r(m_{\rm inv})$ but not the CME signal. The value $\Delta\gamma$(m_{\rm inv}$$>2 GeV/$c^{2}$) $=(4.5\pm 0.8)\times 10^{-5}$ is consistent with the input CME signal, and it would present a $5\sigma$ measurement.

The CME is generally believed to be a low-$p_{T}$ phenomenon Kharzeev:2007jp , and would thus be more prominent in the low $m_{\rm inv}$ region. With a $m_{\rm inv}>2$ GeV/$c^{2}$ cut we used here, the particle average $p_{T}$ is typically 1.2 GeV/$c$ WangQM18 . This is not very high and the CME may still be present above such a mass cut. Moreover, a recent study Shi:2017cpu indicates that the CME signal is rather independent of $p_{T}$ at $p_{T}>0.2$ GeV/$c$, suggesting that the signal can persist to high $m_{\rm inv}$. Nevertheless, our proposal to apply a lower $m_{\rm inv}$ cut will eliminate resonance contributions to $\Delta\gamma$; any measured remaining positive $\Delta\gamma$ would point to the interesting possibility of the existence of the CME. A null measurement at high $m_{\rm inv}$, on the other hand, does not necessarily mean null CME also at low $m_{\rm inv}$.

III Two-component model fit in low-$m_{\rm inv}$ region

In what follows, we illustrate a fit method to potentially identify the possible CME at low $m_{\rm inv}$. Still use $\rho\to\pi^{+}\pi^{-}$ as an example, and consider the event to be composed of primordial pions containing CME signals ($\gamma_{\rm CME}$) and common (charge-independent) backgrounds, such as momentum conservation ($\gamma_{\rm m.c.}$) Bzdak:2010fd ; Pratt:2010zn , and decay pions containing correlations from the decay Voloshin:2004vk ; Schlichting:2010qia ; Wang:2016iov . We have

[TABLE]

(If one normalized $\gamma_{\rm OS}$ by $N_{\rm SS}$ instead, then Eq. (3) would become simpler, $\Delta\gamma^{\prime}=r^{\prime}\gamma_{\rho}+\Delta\gamma_{\rm CME}$ with $r^{\prime}=(N_{\rm OS}-N_{\rm SS})/N_{\rm SS}$.) Considering Eq. (3), the $\Delta\gamma$ can be expressed by two terms:

[TABLE]

The first term is resonance contributions, where the response function

[TABLE]

contains $v_{2}$ of various resonances and is likely a smooth function of $m_{\rm inv}$ while $r(m_{\rm inv})$ contains resonance spectral profile (Fig. 3(a)). Consequently, the first term is not smooth but a peaked function of $m_{\rm inv}$. The second term in Eq. (4) is the CME signal which should be a smooth function of $m_{\rm inv}$ (note we have dropped the negligible $r/2$). However, the exact functional form of $\Delta\gamma_{\rm CME}(m_{\rm inv})$ is presently unknown and needs theoretical input. It is possible that the CME may possess some resonance shapes; in the extreme case where the CME has the same $m_{\rm inv}$ shape as the background, then our method would obviously fail. Nevertheless, the $m_{\rm inv}$ dependences of the two terms are likely different, and this can be exploited to identify CME signals at low $m_{\rm inv}$. This is illustrated in Fig. 3(c) where the ratio of $\Delta\gamma/r$ is depicted. If CME signal is present, as is the case in our toy MC, $\Delta\gamma/r$ should not be smooth, but with a deviation resembling the inverse shape of $r$ in Fig. 3(a). This is clearly seen in Fig. 3(c) in the $\rho$ mass region, although not as clear in the $K_{S}$ mass region.

In Eq. (4), $\Delta\gamma(m_{\rm inv})$ and $r(m_{\rm inv})$ are measured, and $R(m_{\rm inv})$ results from known physics and can in principle be obtained from models. AMPT indicates that $R(m_{\rm inv})$ is a first-order polynomial. We can thus take a step further to fit the $\Delta\gamma(m_{\rm inv})$ in Fig. 3(b) by Eq. (4) taking $R(m_{\rm inv})$ as a first-order polynomial fit function, treating CME as a $m_{\rm inv}$-independent fit parameter (our input CME signal is $p_{T}$ independent). The fit result is superimposed as the red histogram in Fig. 3(b), and in Fig. 3(c) after divided by $r(m_{\rm inv})$ from Fig. 3(a). The straight line in blue in Fig. 3(c) is the fit result for $R(m_{\rm inv})$. The difference between the fit red histogram and the blue line is $\Delta\gamma_{\rm CME}/r(m_{\rm inv})$, which shows the inverse shape of $r(m_{\rm inv})$. It is found, with the simulated statistics, that the inverse-shape feature becomes hard to identify when the CME input signal is smaller than 10% of the inclusive $\Delta\gamma$. The fit parameters are written in Fig. 3(b). The fit parameter for CME is $\Delta\gamma_{\rm CME}=(4.2\pm 0.2)\times 10^{-5}$, not far away from the input CME signal of $4.6\times 10^{-5}$. The fit $\chi^{2}/{\rm ndf}=108/75$ is not ideal because of the approximation for the $m_{\rm inv}$-dependence of $R(m_{\rm inv})$, but it presents a potentially viable way to extract CME signals from data even at low $m_{\rm inv}$.

Theoretically, the $m_{\rm inv}$ dependence of the CME is unknown. The likely sphaleron or instanton mechanism for transitions between QCD vacuum states Kharzeev:1998kz ; Kharzeev:2007jp ; Fukushima:2008xe ; Kharzeev:2004ey , leading to the CME, might yield a broad $m_{\rm inv}$ distribution around $m_{\rm inv}\sim 1$ GeV/$c$. With such a CME signal, assuming a $m_{\rm inv}$-independent CME in our fit would probably still yield a reasonable average CME signal. To illustrate this point, we simulate sphalerons/instantons by a broad Gaussian mass distribution at 1 GeV/$c^{2}$ with width 0.5 GeV/$c^{2}$. Each sphaleron/instanton is at rest and decays into a $\pi^{+}\pi^{-}$ pair where the decay polar angle is uniform in $\sin\theta$ and the azimuthal angles of $\pi^{\pm}$ are sampled according to $dN/d\phi\propto 1\pm\sin\phi$. The number of sphalerons/instantons is Poisson and on average 0.7% (0.9% within $|\eta|<1$) of the event single charge pion multiplicity. The pion azimuthal distributions in the form of Eq. (2) are $dN_{\pi^{\pm}}/d\phi\propto 1+2v_{2}\cos 2\phi+0.009\times(1\pm\sin\phi)$. This gives an effective CME signal of $2a_{1}^{2}\approx 4.0\times 10^{-5}$. We apply our fit method, still assuming a constant CME, and obtain $\Delta\gamma_{\rm CME}=(3.7\pm 0.2)\times 10^{-5}$, consistent with the input signal, with $\chi^{2}/{\rm ndf}=187/75$.

IV Use ESE to determine background shape in two-component fit

In the TC fit in Sect. III, one needs the functional form of $R(m_{\rm inv})$ as input. The linear form used in the Fig. 3 fit was motivated by AMPT results, but is strongly model dependent. To lift this model dependency, one may resort to event-shape engineering (ESE) Schukraft:2012ah ; Acharya:2017fau ; Sirunyan:2017quh where events in each narrow centrality bin are divided into two classes according to the flow vector $q_{2}$:

[TABLE]

where $M$ is the particle multiplicity and $\phi_{j}$ is the $j$th particle azimuth. The two $q_{2}$ event classes will have different average $v_{2}$ exploiting dynamical fluctuations of $v_{2}$ Schukraft:2012ah ; Acharya:2017fau ; Sirunyan:2017quh . (In our simulation we included dynamical fluctuations of 40% of the average $v_{2}$ Wang:2016iov .) Since the magnetic fields are approximately equal for the two classes while the $v_{2}$ backgrounds differ, the difference in $\Delta\gamma$ between the two event classes is a good representation of the background shape, and can therefore be used to add to the TC fit. We refer to such a fit as TC+ESE fit. Note that in this fit model the background is not required to be strictly proportional to $v_{2}$ WangQM18 .

We calculate $q_{2}$ using particles from $-1<\eta<-0.05$ and $\Delta\gamma$ using particles from $0.05<\eta<1$, and vice versa (which we refer to as the subevent method). Figure 4 left panel shows the $m_{\rm inv}$ distributions of $\Delta\gamma_{A}-\Delta\gamma_{B}$ (scaled by a factor of 5) together with the inclusive $\Delta\gamma$. Note that, because of the different $\eta$ acceptances, the $\Delta\gamma$’s in Fig. 4 and Fig. 3 are numerically different. The inclusive $\Delta\gamma$ contains both background and CME. With the background shape given by $\Delta\gamma_{A}-\Delta\gamma_{B}$, the CME can be extracted from a fit $\Delta\gamma=k(\Delta\gamma_{A}-\Delta\gamma_{B})+\Delta\gamma_{\rm CME}$. Since the same data are used in $\Delta\gamma$ and $\Delta\gamma_{A}-\Delta\gamma_{B}$, their statistical errors are somewhat correlated. To properly handle statistical errors, one can simply fit the independent measurements of $\Delta\gamma_{A}$ versus $\Delta\gamma_{B}$, namely $\Delta\gamma_{A}=b\Delta\gamma_{B}+(1-b)\Delta\gamma_{\rm CME}$ where $b$ and $\Delta\gamma_{\rm CME}$ are the fit parameters. Figure 4 right panel shows such a fit to our toy model simulation.

Figure 5 shows the $\Delta\gamma_{\rm CME}$ parameters from the TC fit and the TC+ESE fit as functions of the input CME strength. The TC+ESE fit uncertainties on the $\Delta\gamma_{\rm CME}$ parameter are relatively large. This is because $\Delta\gamma_{A}-\Delta\gamma_{B}$ is small (see Fig. 4 left panel), corresponding to a relative difference in $v_{2}$ on the order of 10% in our toy model simulation representing mid-central Au+Au collisions at RHIC. At the LHC, this difference is larger, approximately 20% Acharya:2017fau ; Sirunyan:2017quh , because of the larger event multiplicities than at RHIC. As seen from Fig. 5, the fits appear to faithfully reproduce the input CME signal.

V Discussion and summary

The STAR experiment at RHIC has accumulated Au+Au minimum bias data samples of $15\times 10^{6}$ events from Run-4 Abelev:2009ac ; Abelev:2009ad , $57\times 10^{6}$ events from Run-7 Adamczyk:2014mzf , and $500\times 10^{6}$ events from Run-11 Skokov:2016yrj . If the CME signal is 1/3 of the measured inclusive $\Delta\gamma$ Abelev:2009ac ; Abelev:2009ad and independent of $m_{\rm inv}$, these data samples (with the 20-60% centrality) could yield, based on our toy MC study, a better than $5\sigma$ measurement of $\Delta\gamma$(m_{\rm inv}$$>2 GeV/$c^{2}$). If the CME signal is unobservable, then our analysis method could set, based on our AMPT result, an upper limit with 98% CL on the CME of 5% at m_{\rm inv}$$>2 GeV/$c^{2}$ relative to the measured inclusive $\Delta\gamma$. It is likely that the CME contribution decreases with $m_{\rm inv}$ and, depending on the detailed physics mechanism, may become difficult to observe at m_{\rm inv}$$>2 GeV/$c^{2}$. Our fit methods can be used to explore and extract CME signals at low $m_{\rm inv}$. The methods rely on the rather robust assumption of different $m_{\rm inv}$ dependences of peaked resonance contributions and smooth CME signal, lifting the major difficulty of similar dependences of the background and CME on experimental variables thus far. The ESE method can further help determine the background dependence on $m_{\rm inv}$, thus improving the accuracy in the extraction of the CME signal.

In summary, topological charge fluctuations, resulting in the chiral magnetic effect (CME) and charge separation in relativistic heavy ion collisions, are fundamental properties of QCD. Experimental charge separation measurements by the azimuthal correlator ($\Delta\gamma$) suffer from major backgrounds from resonance decays (generally, local charge conservation) coupled with elliptic anisotropy. In this article, we propose to measure the $\Delta\gamma$ differentially as a function of the particle pair invariant mass ($m_{\rm inv}$). By using the AMPT (A Multi-Phase Transport) model, we demonstrate that one can essentially eliminate resonance decay backgrounds to $\Delta\gamma$ by applying a lower cut on $m_{\rm inv}$. With a m_{\rm inv}$$>2 GeV/$c^{2}$ cut, an upper limit on the CME of 20% of the inclusive $\Delta\gamma$ can be achieved with $11\times 10^{6}$ AMPT events of 200 GeV Au+Au collisions with impact parameter $b=6.6$-8.2 fm. By using a toy Monte Carlo simulation with realistic resonance distributions and a $p_{T}$-independent input CME signal, we show that the resonance decay backgrounds are eliminated by the m_{\rm inv}$$>2 GeV/$c^{2}$ cut and the CME signal remains. With input CME signal of $a_{1}=0.008$ (20% of the total $\Delta\gamma$) and $200\times 10^{6}$ events corresponding to the 40-50% centrality of Au+Au collisions, a $5\sigma$ CME measurement could be achieved at m_{\rm inv}$$>2 GeV/$c^{2}$. We further show that one may be able to separate the presumably smooth CME signals from peaked resonance decay backgrounds in the low $m_{\rm inv}$ region by exploiting their different $m_{\rm inv}$ dependences. We show this by the toy MC with a $p_{T}$-independent CME signal as well as the signal from a broad sphaleron/instanton mass distribution. We demonstrate that one may use the ESE to determine the background shape in $m_{\rm inv}$, thus help the extraction of the CME signal.

We note that our lower $m_{\rm inv}$ cut may also remove the CME signal together with the background. A negative conclusion at large $m_{\rm inv}$ does therefore not necessarily mean null CME also at low $m_{\rm inv}$. A positive signal at large $m_{\rm inv}$, however, would be a good indication of the possible CME. We also note that our two-component model fit in the low $m_{\rm inv}$ region is not without assumptions. One uncertainty is the unknown $m_{\rm inv}$ dependence of the CME, where theoretical guidance is much needed. Nevertheless, our proposed invariant mass method provides a valuable tool and should help the ongoing experimental search for the CME at RHIC and the LHC.

Acknowledgments

We thank Prof. Zi-Wei Lin for useful discussions. This work was supported in part by the National Natural Science Foundation of China Grant Nos. 11647306 and 11747312 and the U.S. Department of Energy Grant No. DE-SC0012910. HL acknowledges financial support from the China Scholarship Council.