arxiv: v1 [hep-ex] 23 Jul 2014

Transcription

1 EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH CERN-PH-EP July 1, 014 Search for exclusive photoproduction of Z ± c (3900) at COMPASS The COMPASS Collaboration arxiv: v1 [hep-ex] 3 Jul 014 Abstract A search for the exclusive production of the Z ± c (3900) hadron by virtual photons has been performed in the channel Z ± c (3900) J/ψπ ±. The data cover the range from 7 GeV to 19 GeV in the centreof-mass energy of the photon-nucleon system. The full set of the COMPASS data set collected with a muon beam between 00 and 011 has been used. An upper limit for the ratio BR(Z ± c (3900) J/ψπ ± ) σ γ N Z ± c (3900) N /σ γ N J/ψ N of has been established at the confidence level of 90%. Keywords: COMPASS, Z c (3900), photoproduction, tetraquark. (to be submitted to Phys. Lett. B)

2 The COMPASS Collaboration C. Adolph 8, R. Akhunzyanov 7, M.G. Alexeev 7, G.D. Alexeev 7, A. Amoroso 7,9, V. Andrieux, V. Anosov 7, A. Austregesilo,17, B. Badełek 31, F. Balestra 7,9, J. Barth 4, G. Baum 1, R. Beck 3, Y. Bedfer, A. Berlin, J. Bernhard 13, K. Bicker,17, E. R. Bielert, J. Bieling 4, R. Birsa 5, J. Bisplinghoff 3, M. Bodlak 19, M. Boer, P. Bordalo 1,a, F. Bradamante 4,5, C. Braun 8, A. Bressan 4,5, M. Büchele 9, E. Burtin, L. Capozza, M. Chiosso 7,9, S.U. Chung 17,b, A. Cicuttin 6,5, M.L. Crespo 6,5, Q. Curiel, S. Dalla Torre 5, S.S. Dasgupta 6, S. Dasgupta 5, O.Yu. Denisov 9, S.V. Donskov 1, N. Doshita 33, V. Duic 4, W. Dünnweber 16, M. Dziewiecki 3, A. Efremov 7, C. Elia 4,5, P.D. Eversheim 3, W. Eyrich 8, M. Faessler 16, A. Ferrero, A. Filin 1, M. Finger 19, M. Finger jr. 19, H. Fischer 9, C. Franco 1, N. du Fresne von Hohenesche 13,, J.M. Friedrich 17, V. Frolov, F. Gautheron, O.P. Gavrichtchouk 7, S. Gerassimov 15,17, R. Geyer 16, I. Gnesi 7,9, B. Gobbo 5, S. Goertz 4, M. Gorzellik 9, S. Grabmüller 17, A. Grasso 7,9, B. Grube 17, T. Grussenmeyer 9, A. Guskov 7, F. Haas 17, D. von Harrach 13, D. Hahne 4, R. Hashimoto 33, F.H. Heinsius 9, F. Herrmann 9, F. Hinterberger 3, Ch. Höppner 17, N. Horikawa 18,d, N. d Hose, S. Huber 17, S. Ishimoto 33,e, A. Ivanov 7, Yu. Ivanshin 7, T. Iwata 33, R. Jahn 3, V. Jary 0, P. Jasinski 13, P. Jörg 9, R. Joosten 3, E. Kabuß 13, B. Ketzer 17,f, G.V. Khaustov 1, Yu.A. Khokhlov 1,g, Yu. Kisselev 7, F. Klein 4, K. Klimaszewski 30, J.H. Koivuniemi, V.N. Kolosov 1, K. Kondo 33, K. Königsmann 9, I. Konorov 15,17, V.F. Konstantinov 1, A.M. Kotzinian 7,9, O. Kouznetsov 7, M. Krämer 17, Z.V. Kroumchtein 7, N. Kuchinski 7, F. Kunne, K. Kurek 30, R.P. Kurjata 3, A.A. Lednev 1, A. Lehmann 8, M. Levillain, S. Levorato 5, J. Lichtenstadt 3, A. Maggiora 9, A. Magnon, N. Makke 4,5, G.K. Mallot, C. Marchand, A. Martin 4,5, J. Marzec 3, J. Matousek 19, H. Matsuda 33, T. Matsuda 14, G. Meshcheryakov 7, W. Meyer, T. Michigami 33, Yu.V. Mikhailov 1, Y. Miyachi 33, A. Nagaytsev 7, T. Nagel 17, F. Nerling 13, S. Neubert 17, D. Neyret, V.I. Nikolaenko 1, J. Novy 0, W.-D. Nowak 9, A.S. Nunes 1, A.G. Olshevsky 7, I. Orlov 7, M. Ostrick 13, R. Panknin 4, D. Panzieri 8,9, B. Parsamyan 7,9, S. Paul 17, D.V. Peshekhonov 7, S. Platchkov, J. Pochodzalla 13, V.A. Polyakov 1, J. Pretz 4,h, M. Quaresma 1, C. Quintans 1, S. Ramos 1,a, C. Regali 9, G. Reicherz, E. Rocco, N.S. Rossiyskaya 7, D.I. Ryabchikov 1, A. Rychter 3, V.D. Samoylenko 1, A. Sandacz 30, S. Sarkar 6, I.A. Savin 7, G. Sbrizzai 4,5, P. Schiavon 4,5, C. Schill 9, T. Schlüter 16, K. Schmidt 9,c, H. Schmieden 4, K. Schönning, S. Schopferer 9, M. Schott, O.Yu. Shevchenko 7,*, L. Silva 1, L. Sinha 6, S. Sirtl 9, M. Slunecka 7, S. Sosio 7,9, F. Sozzi 5, A. Srnka 5, L. Steiger 5, M. Stolarski 1, M. Sulc 11, R. Sulej 30, H. Suzuki 33,d, A. Szabelski 30, T. Szameitat 9,c, P. Sznajder 30, S. Takekawa 7,9, J. ter Wolbeek 9,c, S. Tessaro 5, F. Tessarotto 5, F. Thibaud, S. Uhl 17, I. Uman 16, M. Virius 0, L. Wang, T. Weisrock 13, M. Wilfert 13, R. Windmolders 4, H. Wollny, K. Zaremba 3, M. Zavertyaev 15, E. Zemlyanichkina 7 and M. Ziembicki 3, A. Zink 8 1 Universität Bielefeld, Fakultät für Physik, Bielefeld, Germany i Universität Bochum, Institut für Experimentalphysik, Bochum, Germany ip 3 Universität Bonn, Helmholtz-Institut für Strahlen- und Kernphysik, Bonn, Germany i 4 Universität Bonn, Physikalisches Institut, Bonn, Germany i 5 Institute of Scientific Instruments, AS CR, 6164 Brno, Czech Republic j 6 Matrivani Institute of Experimental Research & Education, Calcutta , India k 7 Joint Institute for Nuclear Research, Dubna, Moscow region, Russia l 8 Universität Erlangen Nürnberg, Physikalisches Institut, 954 Erlangen, Germany i 9 Universität Freiburg, Physikalisches Institut, 794 Freiburg, Germany ip CERN, 111 Geneva 3, Switzerland 11 Technical University in Liberec, Liberec, Czech Republic j 1 LIP, Lisbon, Portugal m 13 Universität Mainz, Institut für Kernphysik, Mainz, Germany i 14 University of Miyazaki, Miyazaki , Japan n

3 15 Lebedev Physical Institute, Moscow, Russia 16 Ludwig-Maximilians-Universität München, Department für Physik, Munich, Germany io 17 Technische Universität München, Physik Department, Garching, Germany io 18 Nagoya University, 464 Nagoya, Japan n 19 Charles University in Prague, Faculty of Mathematics and Physics, Prague, Czech Republic j 0 Czech Technical University in Prague, Prague, Czech Republic j 1 State Scientific Center Institute for High Energy Physics of National Research Center Kurchatov Institute, 1481 Protvino, Russia CEA IRFU/SPhN Saclay, Gif-sur-Yvette, France p 3 Tel Aviv University, School of Physics and Astronomy, Tel Aviv, Israel q 4 University of Trieste, Department of Physics, 3417 Trieste, Italy 5 Trieste Section of INFN, 3417 Trieste, Italy 6 Abdus Salam ICTP, Trieste, Italy 7 University of Turin, Department of Physics, 15 Turin, Italy 8 University of Eastern Piedmont, 150 Alessandria, Italy 9 Torino Section of INFN, 15 Turin, Italy 30 National Centre for Nuclear Research, Warsaw, Poland r 31 University of Warsaw, Faculty of Physics, Warsaw, Poland r 3 Warsaw University of Technology, Institute of Radioelectronics, Warsaw, Poland r 33 Yamagata University, Yamagata, Japan n a Also at Instituto Superior Técnico, Universidade de Lisboa, Lisbon, Portugal b Also at Department of Physics, Pusan National University, Busan , Republic of Korea and at Physics Department, Brookhaven National Laboratory, Upton, NY 11973, U.S.A. c Supported by the DFG Research Training Group Programme 1 Physics at Hadron Accelerators d Also at Chubu University, Kasugai, Aichi, Japan n e Also at KEK, 1-1 Oho, Tsukuba, Ibaraki, Japan f Present address: Universität Bonn, Helmholtz-Institut für Strahlen- und Kernphysik, Bonn, Germany g Also at Moscow Institute of Physics and Technology, Moscow Region, , Russia h present address: RWTH Aachen University, III. Physikalisches Institut, 5056 Aachen, Germany i Supported by the German Bundesministerium für Bildung und Forschung j Supported by Czech Republic MEYS Grants ME49 and LA4 k Supported by SAIL (CSR), Govt. of India l Supported by CERN-RFBR Grants and m Supported by the Portuguese FCT - Fundação para a Ciência e Tecnologia, COMPETE and QREN, Grants CERN/FP/933/009, CERN/FP/116376/0 and CERN/FP/13600/011 n Supported by the MEXT and the JSPS under the Grants No , No and No ; Daiko Foundation and Yamada Foundation o Supported by the DFG cluster of excellence Origin and Structure of the Universe ( p Supported by EU FP7 (HadronPhysics3, Grant Agreement number 8386) q Supported by the Israel Science Foundation, founded by the Israel Academy of Sciences and Humanities r Supported by the Polish NCN Grant DEC-011/01/M/ST/0350 * Deceased

4 The Z ± c (3900) state was recently discovered by the BES-III and Belle Collaborations in e + e π + π J/ψ reactions at s = 4.6 GeV [1, ] via the decay channel Z ± c (3900) J/ψπ ±. (1) It has been interpreted as a tetraquark state [3 6], although other explanations like a molecular state [7 11], a cusp effect [1] and an initial-single-pion-emission mechanism [13] were also proposed. According to the vector meson dominance (VMD) model, a photon may behave like a J/ψ so that a Z ± c (3900) can be produced by the interaction of an incoming photon with a virtual charged pion provided by the target nucleon γ N Z ± c (3900) N. () The corresponding diagram is shown in Fig. 1a. Based on the VMD model, the authors of Ref. [14] predict a sizable cross section of the reaction in Eq. () for s γn GeV. Under the assumption that the decay channel of Eq. (1) is dominant and that the total width Γtot of the Z ± c (3900) particle is 46 MeV/c, as measured by BES-III, the cross section reaches a maximum value of 50 nb to 0 nb at s γn = 7 GeV. The J/ψ production in photonnucleon interactions at COMPASS covers the range s γn from 7 GeV to 19 GeV and thus can be used to also study Z ± c (3900) production and to estimate the partial width Γ J/ψπ of the decay channel Z ± c (3900) J/ψπ ±. γ J/ψ + Z c J/ψ γ J/ψ J/ψ p π + π + n p pomeron π + n (a) (b) Fig. 1: Diagrams for (a) Z c + (3900) production via virtual π + exchange and (b) J/ψπ + production via pomeron exchange. The COMPASS experiment [15] is situated at the M beam line of the CERN Super Proton Synchrotron. The data used in the present analysis were obtained scattering positive muons of 160 GeV/c (00-0) or 00 GeV/c momentum (011) off solid 6 LiD (00-004) or NH 3 targets ( ). The longitudinally or transversely polarised targets consisted of two (00 004) or three ( ) cylindrical cells placed along the beam direction. Polarisation effects were canceled out by combining data with opposite polarisation orientations. Particle tracking and identification were performed in a two-stage spectrometer, covering a wide kinematic range. The trigger system comprises hodoscope counters and hadron calorimeters. Beam halo was rejected by veto counters upstream of the target. In the analysis presented in this Letter, the reaction µ + N µ + Z ± c (3900) N µ + J/ψπ ± N µ + µ + µ π ± N (3) was searched for. In order to select samples of exclusive µ + J/ψπ ± events, a reconstructed vertex in the target region with an incoming beam track and three outgoing muon tracks (two positive and one

5 counts/0.01 GeV/c 4 all exclusive 3 counts/1 GeV 3 µ + J/ψ µ + J/ψπ ± M µ + µ [GeV/c ] (a) E [GeV] (b) Fig. : (a) The dimuon mass distribution for all dimuons produced in muon-nucleon scattering (blue, upper curve), and for exclusively produced dimuons (yellow, lower curve). (b) Distribution for the energy balance E in the reactions Eq. (7) (yellow, upper curve) and Eq. (3) (green, lower curve). negative) is required. Tracks are attributed to muons if they cross more than 15 radiation lengths of material. Only the events with exactly three muons and one pion in the final state were selected. A pair of muons is treated as J/ψ candidate if the difference between its reconstructed mass M µ + µ (Fig. a) and the nominal J/ψ mass is less than 150 MeV/c. In case both µ + µ combinations satisfy this condition, the event is rejected. Except for the tiny recoil of the target nucleon, the sum of the scattered muon energy, E µ, and the energies of produced J/ψ and π ± mesons, E J/ψ and E π ±, should be equal to the beam energy E b for the exclusive reaction of Eq. (3). The distribution of events as a function of the energy balance E = E µ + E J/ψ + E π ± E b is presented in Fig. b. With the experimental energy resolution of about 3 GeV, the energy balance is required to be E < GeV. The distribution of the negative squared four-momentum transfer Q = (P b P µ ) is shown in Fig. 3a. Here P µ and P b are four-momenta of the scattered and incident muons, respectively. The momentum of the produced pion is required to be larger than GeV/c in order to reduce the background of exclusive events with a J/ψ and a π ± in the final state produced via pomeron exchange (Fig. 1b). The total number of selected µ + J/ψπ + and µ + J/ψπ events is 565 and 405, respectively. The distribution of the centre-of-mass energy of the photon-nucleon system s γn is shown in Fig. 3b. The mass spectrum for J/ψπ ± events is shown in Fig. 4a. It does not exhibit any statistically significant resonant structure around 3.9 GeV/c. The observed number of events N J/ψπ in the signal range 3.84 GeV/c < M J/ψπ + < 3.96 GeV/c is treated as consisting of an a priori unknown Z c ± (3900) signal N Zc and a background contribution N bkg. According to the method described in Ref. [16], the probability density function g(n Zc ) is given by g(n Zc ) = n 0 e (N Zc +Nbkg) (N Zc + N bkg ) N Zc f(n bkg )dn bkg, (4) N Zc! where n is a normalization constant and the probability density function f(n bkg ), assumed to be Gaussian, describes the background contribution in the signal interval. The mean value and the Gaussian width of f(n bkg ) are estimated by fitting a sum of two exponential functions (A e am J/ψπ + B e bm J/ψπ) to the J/ψπ ± mass spectrum in the range 3.3 GeV/c < M J/ψπ + < 6.0 GeV/c excluding the signal region. The fitted function is shown as a line in Fig. 4a. The number of expected background events in the signal region is 49.7 ± 3.4 while 51 is observed. The upper limit NZ UL c for the number of produced Z c ± (3900) events corresponding to a confidence level of CL = 90% is then determined from the expression N UL Zc 0 g(n Zc ) = 0.9 (5)

6 to be NZ UL c = 15.1 events. For the absolute normalization of the Z ± c (3900) production rate we estimated for the same data sample the number of exclusively produced J/ψ mesons from incoherent exclusive production in γ N J/ψ N, (6) the cross section of which is known for our range of s γn [17]. The same selection criteria are applied for the exclusive production of the J/ψ mesons µ + N µ + J/ψ N, (7) and Z c ± (3900) hadrons. To separate J/ψ production and non-resonant production of dimuons, the dimuon mass spectrum is fitted by a function consisting of three Gaussian (two to describe the J/ψ peak and one for the ψ(s) peak) and an exponential background under the peaks (see Fig. a). Finally events of exclusive J/ψ production remain in the sample. The distribution of the squared transverse momentum p T of the J/ψ (Fig. 4b) for the exclusive sample is fitted by a sum of two exponential functions in order to separate the contributions from exclusive coherent production on the target nuclei and exclusive production on (quasi-)free target nucleons. The contribution from coherent production is found to be 30.3% for the 6 LiD target and 38.9% for NH 3 target (36.1% averaged over the sample). The amount of nonexclusive events in the exclusive incoherent sample is estimated to be about 30 ± %. Since only the charged pion distinguishes the final state of process in Eq. () from the final state of the process in Eq. (6), the ratio R a of their acceptances is in a first approximation equal to the acceptance for this pion. Based on previous COMPASS measurements and Monte Carlo simulations this ratio is about R a = 0.5, averaged over all setup and target configurations. Thus we obtain the result BR(Z ± c (3900) J/ψπ ± ) σ γ N Z ± c (3900) N σ γ N J/ψ N sγn =13.8 GeV < 3.7 3, (8) where BR denotes the branching ratio for the Z c ± (3900) J/ψπ ± decay channel. Assuming σ γ N J/ψ N = 14.0±1.6 stat ±.5 syst nb as measured by the NA14 Collaboration for s γn = 13.7 GeV [17], the result can be presented as BR(Z c ± (3900) J/ψπ ± ) σ γ N Z ± c (3900) N < 5 pb. (9) sγn =13.8 GeV The upper limits for the ratio of the cross sections in intervals of s γn are presented in Table 1. The result shown in Eq. (9) can be converted into an upper limit for the partial width Γ J/ψπ of the decay in Eq. (1) based on the VMD model. According to Ref. [14] the cross section for the reaction in Eq. (), averaged over the measured s γn distribution for J/ψπ ± events is about Γ J/ψπ 430 pb/mev for Λ π = 0.6 GeV, a free parameter of the πnn vertex, yielding Γ J/ψπ Γ tot σ γ N Z ± c (3900) N = Γ J/ψπ 430 pb/mev < 5 pb. () Γ tot Assuming Γtot = 46 MeV/c, we obtain an upper limit Γ J/ψπ <.4 MeV/c. We estimate the systematic uncertainty of the result in Eq. (8) to be about 30%, where the main contributions come from limited knowledge of the acceptance ratio a, since the energy spectrum of the expected Z ± c (3900) events is unknown, from systematic effects in estimation of nonexclusive contamination in the reference J/ψ sample and from the background description in the signal range of the J/ψπ spectrum. The uncertainty of σ γ N J/ψ N measurement by NA14 contributes to Eq. (9), so the total systematic uncertainty of this result is about 35%.

7 counts/0.1 (GeV/c) 4 3 µ + J/ψ µ + J/ψπ ± counts/0.1 GeV µ + J/ψ µ + J/ψπ ± Q [(GeV/c) ] (a) s N [GeV] 0 γ (b) Fig. 3: Kinematic distributions for the reactions Eq. (7) (yellow, upper curves) and Eq. (3) (green, lower curves) (a) Q, (b) s γn. counts/0.0 GeV/c counts/0.0 (GeV/c) 3 NH 3 target 6 LiD target [GeV/c ] (a) M J/ψπ p [(GeV/c) ] (b) T Fig. 4: (a) Mass spectrum of the J/ψπ ± state. The fitted function is shown as a line. (b) p T exclusively produced J/ψ mesons off the 6 LiD (blue, lower) and NH 3 (red, upper) targets. distributions for Table 1: Upper limits for Z ± c (3900) production rate for intervals of s γn. Interval s γn, GeV BR(J/ψπ) σ Zc /σ J/ψ, 3 Full sγn < 1.3 GeV GeV < s γn < 14.1 GeV GeV < s γn < 15.4 GeV GeV < s γn

8 No signal of exclusive photoproduction of the Z ± c (3900) state and its decay into J/ψπ ± was found. Therefore an upper limit was determined for the product of the cross section of this process and the relative Z ± c (3900) J/ψπ ± decay probability normalized to the cross section of incoherent exclusive photoproduction of J/ψ mesons. In case Z ± c (3900) is a real hadron state, the decay channel Z ± c (3900) J/ψπ ± can not be the dominant one. This result is a significant input to clarify the nature of the Z ± c (3900) state. We gratefully acknowledge the support of the CERN management and staff as well as the skills and efforts of the technicians of the collaborating institutions. References [1] M. Ablikim, et al., Phys. Rev. Lett. 1, 5001 (013) [arxiv: [hep-ph]]. [] Z. Q. Liu, et al., Phys. Rev. Lett. 1, 500 (013). [3] L. Maiani et al., Phys. Rev. D87, 111 (013) [arxiv: [hep-ph]]. [4] J. M. Dias et al., Phys. Rev. D88, (013) [arxiv: [hep-ph]]. [5] E. Braaten, Phys. Rev. Lett. 111, (013) [arxiv: [hep-ph]]. [6] C.-F. Qiao and L. Tang, [arxiv: [hep-ph]]. [7] Q. Wang, C. Hanhart and Q. Zhao, Phys. Rev. Lett. 111, (013) [arxiv: [hepph]]. [8] C.-Y. Cui et al., [arxiv: [hep-ph]]. [9] J.-R.Zhang, Phys. Rev. D87, (013)[arXiv: [hep-ph]]. [] H.-W. Ke, Z.-T. Wei and X.-Q. Li, Eur. Phys. J. C 73, 561 (013) [arxiv: [hep-ph]]. [11] Y. Dong et al., Phys. Rev. D88, (013) [arxiv: [hep-ph]]. [1] X.-H. Liu and G. Li, Phys. Rev. D88, (013) [arxiv: [hep-ph]]. [13] D.-Y. Chen, X. Liu and T. Matsuki, Phys. Rev. D88, (013) [arxiv: [hep-ph]]. [14] Q.-Y. Lin et al., Phys. Rev. D 88, (013), [arxiv: [hep-ph]]. [15] P. Abbon et al. (COMPASS Collaboration), Nucl. Instrum. Meth. A577, 455 (007), [arxiv: [hep-ex]]. [16] O. Helene, NIM A300, 13 (1991). [17] R. Barate et al., Z.Phys. C33, 505 (1987).