Measurement of the fragmentation fraction ratio fs/fd

Zurich Open Repository and
Archive
University of Zurich
Main Library
Strickhofstrasse 39
CH-8057 Zurich
www.zora.uzh.ch
Year: 2013
Measurement of the fragmentation fraction ratio fs/fd and its dependence on
B meson kinematics
LHCb Collaboration; et al; Bernet, R; Müller, K; Steinkamp, O; Straumann, U; Vollhardt, A
Abstract: The relative production rate of B0s and B 0 mesons is determined with the hadronic decays
B0s→D−s�+ and B 0 → D − K +. The measurement uses data corresponding to 1.0 fb−1 of pp collisions
at a centre-of-mass energy of s√=7 TeV recorded in the forward region with the LHCb experiment. The
ratio of production rates, f s /f d , is measured to be 0.238 ± 0.004 ± 0.015 ± 0.021, where the first
uncertainty is statistical, the second systematic, and the third theoretical. This is combined with a
previous LHCb measurement to obtain f s /f d = 0.256 ± 0.020. The dependence of f s /f d on the
transverse momentum and pseudorapidity of the B meson is determined using the decays B0s→D−s�+
and B 0 → D −�+. There is evidence for a decrease with increasing transverse momentum, whereas the
ratio remains constant as a function of pseudorapidity. In addition, the ratio of branching fractions of
the decays B 0 → D − K + and B 0 → D −�+ is measured to be 0.0822 ± 0.0011 (stat) ± 0.0025 (syst).
DOI: https://doi.org/10.1007/JHEP04(2013)001
Posted at the Zurich Open Repository and Archive, University of Zurich
ZORA URL: https://doi.org/10.5167/uzh-92008
Published Version
Originally published at:
LHCb Collaboration; et al; Bernet, R; Müller, K; Steinkamp, O; Straumann, U; Vollhardt, A (2013).
Measurement of the fragmentation fraction ratio fs/fd and its dependence on B meson kinematics. Journal
of High Energy Physics, 2013:1.
DOI: https://doi.org/10.1007/JHEP04(2013)001
Published for SISSA by
Springer
Received: January 23, 2013
Accepted: March 18, 2013
Published: April 2, 2013
The LHCb collaboration
E-mail: [email protected]
Abstract: The relative production rate of Bs0 and B 0 mesons is determined with the
hadronic decays Bs0 → Ds− π + and B 0 → D− K + . The measurement uses data corre√
sponding to 1.0 fb−1 of pp collisions at a centre-of-mass energy of s = 7 TeV recorded
in the forward region with the LHCb experiment. The ratio of production rates, fs /fd ,
is measured to be 0.238 ± 0.004 ± 0.015 ± 0.021, where the first uncertainty is statistical,
the second systematic, and the third theoretical. This is combined with a previous LHCb
measurement to obtain fs /fd = 0.256 ± 0.020. The dependence of fs /fd on the transverse
momentum and pseudorapidity of the B meson is determined using the decays Bs0 → Ds− π +
and B 0 → D− π + . There is evidence for a decrease with increasing transverse momentum,
whereas the ratio remains constant as a function of pseudorapidity. In addition, the ratio
of branching fractions of the decays B 0 → D− K + and B 0 → D− π + is measured to be
0.0822 ± 0.0011 (stat) ± 0.0025 (syst).
Keywords: Hadron-Hadron Scattering, Branching fraction, B physics, Flavor physics
ArXiv ePrint: 1301.5286
Open Access, Copyright CERN,
for the benefit of the LHCb collaboration
doi:10.1007/JHEP04(2013)001
JHEP04(2013)001
Measurement of the fragmentation fraction ratio
fs/fd and its dependence on B meson kinematics
Contents
1
2 Detector and software
2
3 Event selection
3
4 Event yields
4
5 Systematic uncertainties
7
6 Results
8
7 Conclusions
10
The LHCb collaboration
13
1
Introduction
The ratio of fragmentation fractions fs /fd quantifies the relative production rate of Bs0
mesons with respect to B 0 mesons. Knowledge of this quantity is essential when determining any Bs0 branching fraction at the LHC. The measurement of the branching fraction
of the rare decay Bs0 → µ+ µ− [1] is the prime example where a precise measurement of
fs /fd is crucial for reaching the highest sensitivity in the search for physics beyond the
Standard Model. The branching fractions of a large number of B 0 and B + decays have
been measured to high precision at the B factories [2], but no Bs0 branching fraction is yet
known with sufficiently high precision to be used as a normalisation channel.
The relative production rates of b hadrons are determined by the fragmentation fractions fu , fd , fs , fc and fΛ , which describe the probability that a b quark will hadronize into
a Bq meson (where q = u, d, s, c), or a b baryon, respectively1 . The ratio of fragmentation
fractions fs /fd has been previously measured at LHCb with hadronic [3] and semileptonic
decays [4], and the resulting values were combined [4].
In this paper, the ratio of fragmentation fractions fs /fd is determined using Bs0 →
Ds− π + and B 0 → D− K + decays collected in pp collisions at a centre-of-mass energy of
√
s = 7 TeV, with data corresponding to an integrated luminosity of 1.0 fb−1 recorded
with the LHCb detector. Since the framework of factorization is well applicable to these
decays [5], their ratio of branching fractions is theoretically well understood [6] and their
1
Charge conjugation is implied throughout this paper.
–1–
JHEP04(2013)001
1 Introduction
relative decay rates can be used to determine the ratio of fragmentation fractions for Bs0
and B 0 mesons through
fs
B(B 0 → D− K + ) DK NDs π
=
fd
B(Bs0 → Ds− π + ) Ds π NDK
Vus 2 fK 2 τB 0 1 B(D− → K + π − π − ) DK NDs π
= ΦPS ,
Vud fπ
τBs0 Na NF B(Ds− → K + K − π − ) Ds π NDK
(1.1)
2
Detector and software
The LHCb detector [12] is a single-arm forward spectrometer covering the pseudorapidity
range 2 < η < 5, designed for the study of particles containing b or c quarks. The detector
includes a high precision tracking system consisting of a silicon-strip vertex detector surrounding the pp interaction region, a large-area silicon-strip detector located upstream of
a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip
detectors and straw drift tubes placed downstream. Data are taken with both magnet
polarities. The combined tracking system has momentum resolution ∆p/p that varies from
–2–
JHEP04(2013)001
where N corresponds to a signal yield, corresponds to a total efficiency, τBs0 /τB 0 =
0.984 ± 0.011 [7] corresponds to the ratio of lifetimes and B(D− → K + π − π − ) = (9.14 ±
−
0.20)% [8] and B(Ds− → K + K − π − ) = (5.50 ± 0.27)% [9] correspond to the D(s)
meson
branching fractions. The factor Na = 1.00 ± 0.02 accounts for the ratio of non-factorizable
0 → D − form factors [11], and
corrections [10], NF = 1.092 ± 0.093 for the ratio of B(s)
(s)
ΦPS = 0.971 for the difference in phase space due to the mass differences of the initial and
final state particles. The numerical values used for the CKM matrix elements are |Vus | =
0.2252, |Vud | = 0.97425, and for the decay constants are fπ = 130.41 MeV, fK = 156.1 MeV,
with negligible uncertainties, below 1% [2]. The measurement is not statistically limited by
the size of the B 0 → D− K + sample , and therefore the theoretically less clean B 0 → D− π +
decays, where exchange diagrams contribute to the total amplitude, do not contribute to
the knowledge of fs /fd .
The ratio of fragmentation fractions can depend on the centre-of-mass energy, as well
0 meson, as was studied previously at LHCb with partially
as on the kinematics of the B(s)
reconstructed B decays [4]. The dependence of the ratio of fragmentation fractions on
0 meson is determined using
the transverse momentum pT and pseudorapidity η of the B(s)
fully reconstructed B 0 → D− π + and Bs0 → Ds− π + decays. Since it is only the dependence
that is of interest here, the more abundant B 0 → D− π + decay is used rather than the
B 0 → D− K + decay. The B 0 → D− K + and B 0 → D− π + decays are also used to determine
their ratio of branching fractions, which can be used to quantify non-factorizable effects in
such heavy-to-light decays [10].
The paper is organised as follows: the detector is described in section 2, followed
by the event selection and the relative selection efficiencies in section 3. The fit to the
mass distributions and the determination of the signal yields are discussed in section 4.
The systematic uncertainties are presented in section 5, and the final results are given in
section 6.
3
Event selection
The three decay modes, B 0 → D− π + , B 0 → D− K + and Bs0 → Ds− π + , are topologically
very similar and can therefore be selected using the same event selection criteria, thus
0 candidates are reconstructed
minimizing efficiency differences between the modes. The B(s)
−
from a D(s)
candidate and an additional pion or kaon (the “bachelor” particle), with the
−
D(s) meson decaying to K + π − π − (K + K − π − ).
0 and D − masses,
After the trigger selection, a loose preselection is made using the B(s)
(s)
lifetimes and vertex qualities. A boosted decision tree (BDT) [21] is used to further separate
signal from background. The BDT is trained on half the Bs0 → Ds− π + data sample. The
0 impact parameter χ2 , the pointing angle of the
most discriminating variables are the B(s)
0
B(s) candidate to the primary vertex and the pT of the tracks. A cut value for the BDT
output variable was chosen to optimally reduce the number of combinatorial background
events, retaining approximately 84% of the signal events.
−
The D(s)
candidates are identified by requiring the invariant mass under the K + π − π −
(K + K − π − ) hypothesis to fall within the selection window 1844–1890 (1944–1990) MeV/c2 .
The relative efficiency of the selection procedure is evaluated for all decay modes using
simulated events, generated with the appropriate Dalitz plot structures [22, 23]. Since the
analysis is only sensitive to relative efficiencies, the impact of any discrepancy between
data and simulation is small.
2
Impact parameter (IP) is defined as the transverse distance of closest approach between the track and
a primary interaction.
–3–
JHEP04(2013)001
0.4% at 5 GeV/c to 0.6% at 100 GeV/c, and impact parameter2 resolution of 20 µm for tracks
with high transverse momentum. Charged hadrons are identified using two ring-imaging
Cherenkov detectors.
The trigger [13] consists of a hardware stage, based on information from the calorimeter
and muon systems, followed by a software stage which applies a full event reconstruction.
The events used in this analysis are selected at the hardware stage by requiring a cluster in
the calorimeters with transverse energy larger than 3.6 GeV. The software stage requires
a two-, three- or four-track secondary vertex with a high sum of the pT of the tracks and
a significant displacement from the primary pp interaction vertices (PVs). At least one
track should have pT greater than 1.7 GeV/c, track fit χ2 over the number of degrees of
freedom less than two, and IP χ2 with respect to the associated primary interaction greater
than sixteen. The IP χ2 is defined as the difference between the χ2 from the vertex fit of
the associated PV reconstructed with and without the considered track. A multivariate
algorithm is used for the identification of the secondary vertices consistent with the decay
of a b hadron.
In the simulation, pp collisions are generated using Pythia 6.4 [14] with a specific
LHCb configuration [15]. Decays of hadronic particles are described by EvtGen [16],
whilst final state radiation is generated using Photos [17]. The interaction of the generated
particles with the detector and its response are implemented using the Geant4 toolkit [18,
19] as described in ref. [20].
4
Event yields
The relative yields of the three decay modes are determined from unbinned extended maxi0 candidates as shown
mum likelihood fits to the mass distributions of the reconstructed B(s)
in figure 1. In order to achieve the highest sensitivity, the sample is separated according
to the two magnet polarities, allowing for possible differences in PID performance and in
running conditions. A simultaneous fit to the two magnet polarities is performed for each
decay mode, with the peak position and width of each signal shared between the two.
The signal mass shape is described by a Gaussian distribution with power-law tails on
either side to model the radiative tail and non-Gaussian detector effects. It consists of a
Crystal Ball function [24]

(m−µ)2
m−µ


for
> −α
 e− 2σ2 ,
σ
−n
fleft (m, α, n, µ, σ) = N · n n
n
m−µ
m−µ
2


· e−|α| /2 ·
−|α|−
, for
≤ −α

|α|
|α|
σ
σ
(4.1)
and a second, similar but mirrored, function to describe the right tail, resulting in the
signal mass shape f2CB (m) = fleft (m) + fright (m). The parameters of the tails are obtained
from simulated events. The mean µ and the width σ of the Gaussian distribution are equal
in both Crystal Ball functions, and are allowed to vary in the fit. The parameter N is a
normalisation factor.
Three classes of background are considered in the fit: fully reconstructed decays where
at least one track is misidentified, partially reconstructed decays with or without misidentified tracks and combinatorial background. The shapes of the invariant mass distributions
for the partially reconstructed decays are taken from large samples of simulated events. The
main sources are B 0 → D− ρ+ and B 0 → D∗− π + (K + ) for the B 0 → D− π + (K + ) sample,
and Bs0 → Ds− ρ+ and Bs0 → Ds∗− π + for the Bs0 → Ds− π + sample.
–4–
JHEP04(2013)001
The final Bs0 → Ds− π + , B 0 → D− π + and B 0 → D− K + event samples are obtained
after particle identification (PID) criteria, based on the difference in log-likelihood between
the kaon and pion hypotheses (DLL). The PID performance as a function of pT and η of
the track is estimated from data using a calibration sample of approximately 27 million
D∗− → D0 (K + π − )π − decays, which are selected using kinematic criteria only. A cut
on the bachelor particle is placed at DLL(K − π)< 0 to select the Bs0 → Ds− π + and
B 0 → D− π + event samples and at DLL(K − π)> 5 to select the B 0 → D− K + sample.
These requirements have an average efficiency of 85.5% and 73.0% respectively with a
misidentification probability of 8.81% and 2.77%. The Ds− → K + K − π − decay is further
distinguished from D− → K + π − π − decays by imposing DLL(K − π)> 5 on the kaon
candidate with the same charge as the D meson, whilst the DLL criteria for the π − and
−
K + are identical between D− and Ds− and are used to discriminate D(s)
decays from
−
background. The total (PID and invariant mass) efficiency to select the D (Ds− ) particle
is 84.6% (78.5%) with a misidentification probability of 4.57% (0.77%).
Candidates / ( 5 MeV/ c2)
12000
8000
−
B 0→ D K +
0
−
Λb → Λc π +
*
−
B 0 → D π+
−
B 0→ D ρ+
Combinatorial
6000
4000
2000
1000
5200
5400
5600
5800
m(D−π+) [MeV/ c2]
LHCb
(b)
−
B 0→ D K +
800
*−
B 0→ D K +
600
B 0→ D K *+
400
B 0→ D ρ+
−
−
−
B 0 → D π+
200
0
2000
Combinatorial
5200
5400
5600
5800
−
m(D K +) [MeV/ c2]
LHCb
(c)
−
B 0s → Ds π+
0
−
B → Ds π+
1500
−
B 0→ D π +
0
Λb→
B 0s →
B 0s →
1000
500
−
Λc π+
*−
Ds π+
−
Ds ρ+
Combinatorial
0
5200
5400
5600
5800
m(D−s π+) [MeV/ c2]
Figure 1. Invariant mass distributions of (a) B 0 → D− π + (b) B 0 → D− K + and (c) Bs0 → Ds− π +
candidates. The solid line is the result of the fit and the dotted line indicates the signal. The
stacked background shapes follow the same top-to-bottom order in the legend and the plot. The Bs0
0
and Λb backgrounds in the B 0 → D− π + mass distribution are invisibly small. The resulting signal
yields are listed in table 1. For illustration purposes the figures include events from both magnet
polarities, although they are fitted separately as described in the text.
–5–
JHEP04(2013)001
Candidates / ( 5 MeV/ c2)
−
B 0 → D π+
−
B 0s → Ds π+
10000
0
5000
Candidates / ( 5 MeV/ c2)
LHCb
(a)
Signal
B0 →
Yield
D− π+
106 197 ± 344
B 0 → D− K +
Bs0 →
Ds− π +
7 664 ± 99
17 419 ± 155
Table 1. Yields obtained from the fits to the invariant mass distributions.
–6–
JHEP04(2013)001
The invariant mass distributions of the misidentified decays are affected by the PID
criteria. The shapes are obtained from simulated events, with the appropriate mass hypothesis applied. The distribution is then reweighted in a data-driven way, according to the
particle identification cut efficiency obtained from the calibration sample, which is strongly
dependent on the momentum of the particle.
Despite the small π → K misidentification probability of 2.8%, the largest misidentified
background in the B 0 → D− K + sample originates from Cabibbo-favoured B 0 → D− π +
decays where the bachelor pion is misidentified as a kaon. The shape of this particular
misidentified decay is determined from data using a high purity sample of B 0 → D− π +
decays (see figure 1(a)), obtained by selecting events in a narrow mass window 5200–
5340 MeV/c2 . The yield of this prominent peaking background is allowed to vary in the
fit and is found to be consistent with the expected yield based on the B 0 → D− π + signal
yield and the misidentification probability. The contamination of B 0 → D− π + events
in the Bs0 → Ds− π + sample can be caused by the misidentification of either pion from
the D− decay. The misidentification probability is 2.0% (3.2%) for the higher (lower)
pT pion. After selecting the Ds− candidate within the mass window around the nominal
Ds− mass [2], the number of misidentified pions is reduced to 0.75% (0.02%). The yield
of this background is constrained in the fit, based on the B 0 → D− π + signal yield, the
misidentification probability and their associated uncertainties.
0
+
The yield of Λb → Λ−
c π decays is allowed to vary in the fit. The cross-feeds from
B 0 → D− K + and Bs0 → Ds− π + events in the B 0 → D− π + signal is small, and are constrained to their respective predicted yields. In addition, a contribution from the rare
B 0 → Ds− π + decay is expected with a yield of 3.3% compared to the Bs0 → Ds− π + signal,
and is accounted for accordingly.
The combinatorial background consists of events with random pions and kaons, forming
a fake D− or Ds− candidate, as well as real D− or Ds− mesons, that combine with a random
pion or kaon. The combinatorial background is modelled with an exponential shape.
The results of the fits are presented in figure 1, and the corresponding signal yields
are listed in table 1. The total yields of the decays B 0 → D− π + and B 0 → D− K +
are used to determine the ratio of their branching fractions, while the event yields of
the decays Bs0 → Ds− π + and B 0 → D− K + are used to measure the average ratio of
fragmentation fractions.
The dependence of the relative b-hadron production fractions as a function of the
0 meson is studied in the ranges 2.0 <
transverse momentum and pseudorapidity of the B(s)
η < 5.0 and 1.5 < pT < 40 GeV/c, using B 0 → D− π + and Bs0 → Ds− π + decays. The
Source
(%)
0.7
2.0
1.2
1.0
1.0
0.7
0.5
0.8
0.4
3.1
Bs0→Ds− π +
B 0→D− K +
0.7
2.0
1.1
1.0
1.5
1.0
0.6
1.0
—
3.4
(%)
Bs0→Ds− π +
B 0→D− π +
(%)
2.0 − 2.9
0.8
1.2
1.5
1.1
0.8
[correl.]
[correl.]
[correl.]
3.2 − 3.8
Table 2. Systematic uncertainties for the measurement of the corrected ratio of event yields used
for the measurements of fs /fd and the relative branching fraction of B 0 → D− K + . The systematic
uncertainty in pT and η bins is shown as a range in the last column, and the total systematic
uncertainty is the quadratic sum of the uncorrelated uncertainties. The systematic uncertainties on
the ratio of B 0 → D− π + and Bs0 → Ds− π + yields that are correlated among the bins do not affect
the dependence on pT or η, and are not accounted for in the total systematic uncertainty.
event sample is subdivided in 20 bins in pT and 10 bins in η, with the bin sizes chosen to
obtain approximately equal number of events per bin. The fitting model for each bin is
the same as that for the integrated samples, apart from the treatment of the exponent of
the combinatorial background distribution, which is fixed to the value obtained from the
fits to the integrated sample.
5
Systematic uncertainties
The systematic uncertainties on the measurement of the relative event yields of the B 0 →
D− π + , B 0 → D− K + and Bs0 → Ds− π + decay modes are related to trigger and offline
selection efficiency corrections, particle identification calibration and the fit model.
The response to charged pions and kaons of the hadronic calorimeter used at the
hardware trigger level has been investigated. As the hardware trigger mostly triggers on
the high-pT bachelor, a systematic uncertainty of 2% is assigned to the ratio of trigger
efficiencies for the decays B 0 → D− K + and B 0 → D− π + , estimated from dedicated studies
with D∗− → D0 (K + π − )π − decays. This uncertainty is assumed to be uncorrelated between
the individual bins in the binned analysis.
The relative selection efficiencies from simulation are studied by varying the BDT
criterion, changing the signal yields by about ±25%. The variation of the relative efficiency
is 1.0% which is assigned as systematic uncertainty.
The uncertainty on the PID efficiencies is estimated by comparing, in simulated events,
the results obtained using the D∗− calibration sample to the true simulated PID perfor-
–7–
JHEP04(2013)001
Detector acceptance
and reconstruction
Hardware trigger efficiency
Offline selection
BDT cut
PID selection
Comb. background
Signal shape (tails)
Signal shape (core)
Charmless background
Total
B 0→D− π +
B 0→D− K +
6
Results
The relative signal yields of the decays B 0 → D− π + , B 0 → D− K + and Bs0 → Ds− π +
are used to determine the branching fraction of the decay B 0 → D− K + , and the ratio of
fragmentation fractions fs /fd .
The efficiency corrected ratio of B 0 → D− K + and B 0 → D− π + signal yields results in
the ratio of branching fractions
B B 0 → D− K +
= 0.0822 ± 0.0011 (stat) ± 0.0025 (syst).
B (B 0 → D− π + )
–8–
JHEP04(2013)001
mance on the signal decays. The corresponding uncertainty ranges from 1.0% to 1.5% for
the different measurements.
The exponent of the combinatorial background distribution is allowed to vary in the fits
to the B 0 → D− π + and Bs0 → Ds− π + mass distributions. By studying D− π − and D− K −
combinations, it is suggested that the value of the exponent is smaller for the B 0 → D− K +
decays than for the B 0 → D− π + decays, and therefore in the fit to the B 0 → D− K +
candidates the exponent is fixed to half the value found in the fit to the B 0 → D− π + sample.
The uncertainty on the signal yields due to the shape of the combinatorial background
is estimated by reducing the exponent to half its value in the fits to the B 0 → D− π +
and Bs0 → Ds− π + mass distributions, and by taking a flat background for the fit to the
B 0 → D− K + mass distribution. An uncertainty of 1.0% (0.7%) is assigned to the relative
B 0 → D− K + and Bs0 → Ds− π + (B 0 → D− π + ) yields.
The tails of the signal distributions are fixed from simulation due to the presence of
large amounts of partially reconstructed decays in the lower sidebands. The uncertainty on
the signal yield is estimated by varying the parameters that describe the tails by 10%. The
uncertainty from the shape of the central peak is taken from a fit allowing for two different
widths for the Crystal Ball functions in eq. 4.1, leading to a 1.0% (0.8%) uncertainty on
the relative B 0 → D− K + and Bs0 → Ds− π + (B 0 → D− π + ) yields.
The contribution of charmless B decays without an intermediate D meson is ignored
in the fit. To evaluate the systematic uncertainty due to these decays, the B mass spectra
for candidates in the sidebands of the D mass distribution are examined. A contribution
of 0.4% relative to the signal yield is found in the B 0 → D− π + decay mode, and no
contribution is seen in the other modes. For the B 0 → D− π + decay mode no correction is
applied and the full size is taken as an uncertainty. No systematic uncertainty is assigned
for the other decay modes.
The various sources of the systematic uncertainty that contribute to the uncertainties
on the ratios of signal yields are listed in table 2. No uncertainty is associated to the
0
+
Λb → Λ−
c π background, as the yield is allowed to vary in the fit. Other cross checks, like
0
+
varying the B 0 → Ds− π + yield in the Bs0 → Ds− π + fit or including Λb → Λ−
c K in the
0
−
+
B → D K fit, show a negligible effect on the signal yields.
All systematic variations are also performed in bins, and the corresponding relative
changes in the ratio of yields have been quantified. Variations showing correlated behaviour
do not affect the slope and are therefore not considered further.
0.35
d
(a)
fs f
LHCb
/
/
fs f
d
0.35
0.3
0.3
0.25
0.25
0.2
0
10000
20000
30000
0.2
40000
p (B) [MeV/ c]
(b)
2
LHCb
3
4
5
η(B)
Figure 2. Ratio of fragmentation fractions fs /fd as functions of (a) pT and (b) η. The errors on the
data points are the statistical and uncorrelated systematic uncertainties added in quadrature. The
solid line is the result of a linear fit, and the dashed line corresponds to the fit for the no-dependence
hypothesis. The average value of pT or η is determined for each bin and used as the center of the
bin. The horizontal error bars indicate the bin size. Note that the scale is zero suppressed.
This is combined with the world average branching fraction B B 0 → D− π + = (26.8 ±
1.3) × 10−4 [2], to give
B B 0 → D− K + = (2.20 ± 0.03 ± 0.07 ± 0.11) × 10−4 ,
where the first uncertainty is statistical, the second is systematic and the last is due to the
uncertainty on the B 0 → D− π + branching fraction.
The ratio of fragmentation fractions is determined from the efficiency corrected event
yields. The ratio of efficiencies is 0.913 ± 0.027. This results in
fs
1
= (0.261 ± 0.004 ± 0.017) ×
fd
Na NF
= 0.238 ± 0.004 ± 0.015 ± 0.021 ,
where the first uncertainty is statistical, the second is systematic containing the sources
listed in table 2 as well as errors from external measurements, and the third is theoretical,
due to the knowledge of Na and NF . The last source is dominated by the uncertainty on
the form factor ratio.
This measurement supersedes and is in agreement with the previous determination with
hadronic decays [3]. It also agrees with the previous measurement based on semileptonic
decays [4]. The two independent results are combined taking into account the various
−
sources of correlated systematic uncertainties, notably the D(s)
branching fractions and
0
B(s) lifetimes, to give
fs
= 0.256 ± 0.020,
fd
which supersedes the previous measurement from LHCb.
–9–
(6.1)
JHEP04(2013)001
T
The value of fs /fd in bins of pT or η is determined using the Bs0 → Ds− π + and B 0 →
D− π + decay modes and is presented in figure 2. A linear χ2 fit gives
fs /fd (pT ) = (0.256 ± 0.020) + (−2.0 ± 0.6) × 10−3 / GeV/c × (pT − hpT i)
fs /fd (η) = (0.256 ± 0.020) + (0.005 ± 0.006) × (η − hηi),
7
Conclusions
The relative production rate of Bs0 and B 0 mesons is determined using the hadronic decays
Bs0 → Ds− π + and B 0 → D− K + resulting in fs /fd = 0.238 ± 0.004(stat) ± 0.015(syst) ±
0.021(theo). This value is consistent with a previous LHCb measurement based on semileptonic decays, with which it is averaged to obtain fs /fd = 0.256 ± 0.020. The ratio of
fragmentation fractions fs /fd is determined as a function of the transverse momentum
0 meson, and a variation consistent with a linear depenand pseudorapidity of the B(s)
0 meson is observed, with a significance
dence on the transverse momentum of the the B(s)
of three standard deviations. In addition, the ratio of branching fractions of the decays
B 0 → D− K + and B 0 → D− π + is measured to be 0.0822 ± 0.0011 (stat) ± 0.0025 (syst).
Acknowledgments
We express our gratitude to our colleagues in the CERN accelerator departments for the
excellent performance of the LHC. We thank the technical and administrative staff at the
LHCb institutes. We acknowledge support from CERN and from the national agencies:
CAPES, CNPq, FAPERJ and FINEP (Brazil); NSFC (China); CNRS/IN2P3 and Region
Auvergne (France); BMBF, DFG, HGF and MPG (Germany); SFI (Ireland); INFN (Italy);
FOM and NWO (The Netherlands); SCSR (Poland); ANCS/IFA (Romania); MinES,
Rosatom, RFBR and NRC “Kurchatov Institute” (Russia); MinECo, XuntaGal and GENCAT (Spain); SNSF and SER (Switzerland); NAS Ukraine (Ukraine); STFC (United Kingdom); NSF (USA). We also acknowledge the support received from the ERC under FP7.
The Tier1 computing centres are supported by IN2P3 (France), KIT and BMBF (Germany), INFN (Italy), NWO and SURF (The Netherlands), PIC (Spain), GridPP (United
Kingdom). We are thankful for the computing resources put at our disposal by Yandex
LLC (Russia), as well as to the communities behind the multiple open source software
packages that we depend on.
Open Access. This article is distributed under the terms of the Creative Commons
Attribution License which permits any use, distribution and reproduction in any medium,
provided the original author(s) and source are credited.
– 10 –
JHEP04(2013)001
with hpT i = 10.4 GeV/c and hηi = 3.28. The data points are normalised with a scale factor
to match the average value of 0.256. The uncertainty associated to this parameter is taken
from eq. 6.1, whilst the error from the fit is 0.003 for both pT and η.
The p-value for this linear fit is found to be 0.16 (0.87) for pT (η). The observed slope
0 meson deviates from zero with
for the dependence on the transverse momentum of the B(s)
a significance of three standard deviations. No indication of a dependence on η(B) is found.
References
[1] LHCb collaboration, First evidence for the decay Bs → µ+ µ− , Phys. Rev. Lett. 110 (2013)
021801 [arXiv:1211.2674] [INSPIRE].
[2] Particle Data Group collaboration, J. Beringer et al., Review of particle physics, Phys.
Rev. D 86 (2012) 010001 [INSPIRE].
[3] LHCb collaboration, Determination of fs /fd for 7 TeV pp collisions and a measurement of
the branching fraction of the decay Bd → D− K + ”, Phys. Rev. Lett. 107 (2011) 211801
[arXiv:1106.4435] [INSPIRE].
[5] M. Beneke, G. Buchalla, M. Neubert and C.T. Sachrajda, QCD factorization for exclusive,
nonleptonic B meson decays: General arguments and the case of heavy light final states,
Nucl. Phys. B 591 (2000) 313 [hep-ph/0006124] [INSPIRE].
[6] R. Fleischer, N. Serra and N. Tuning, A new strategy for Bs branching ratio measurements
and the search for new physics in Bs0 → µ+ µ− , Phys. Rev. D 82 (2010) 034038
[arXiv:1004.3982] [INSPIRE].
[7] Heavy Flavor Averaging Group collaboration, Y. Amhis et al., Averages of b-hadron,
c-hadron and τ -lepton properties as of early 2012, arXiv:1207.1158 [INSPIRE].
[8] CLEO collaboration, S. Dobbs et al., Measurement of absolute hadronic branching fractions
of D mesons and e+ e− → DD̄ cross-sections at the ψ(3770), Phys. Rev. D 76 (2007) 112001
[arXiv:0709.3783] [INSPIRE].
[9] CLEO collaboration, J. Alexander et al., Absolute measurement of hadronic branching
fractions of the Ds+ meson, Phys. Rev. Lett. 100 (2008) 161804 [arXiv:0801.0680]
[INSPIRE].
[10] R. Fleischer, N. Serra and N. Tuning, Tests of Factorization and SU(3) Relations in B Decays
into Heavy-Light Final States, Phys. Rev. D 83 (2011) 014017 [arXiv:1012.2784] [INSPIRE].
[11] J.A. Bailey et al., Bs → Ds /B → D semileptonic form-factor ratios and their application to
BR(Bs0 → µ+ µ− ), Phys. Rev. D 85 (2012) 114502 [Erratum ibid. D 86 (2012) 039904]
[arXiv:1202.6346] [INSPIRE].
[12] LHCb collaboration, The LHCb detector at the LHC, 2008 JINST 3 S08005 [INSPIRE].
[13] R. Aaij et al., The LHCb trigger and its performance, arXiv:1211.3055 [INSPIRE].
[14] T. Sjöstrand, S. Mrenna and P.Z. Skands, PYTHIA 6.4 physics and manual, JHEP 05
(2006) 026 [hep-ph/0603175] [INSPIRE].
[15] I. Belyaev et al., Handling of the generation of primary events in Gauss, the LHCb
simulation framework, IEEE Nucl. Sci. Symp. Conf. Rec. (2010) 1155.
[16] D. Lange, The EvtGen particle decay simulation package, Nucl. Instrum. Meth. A 462
(2001) 152 [INSPIRE].
[17] P. Golonka and Z. Was, PHOTOS Monte Carlo: a precision tool for QED corrections in Z
and W decays, Eur. Phys. J. C 45 (2006) 97 [hep-ph/0506026] [INSPIRE].
[18] GEANT4 collaboration, J. Allison et al., GEANT4 developments and applications, IEEE
Trans. Nucl. Sci. 53 (2006) 270.
– 11 –
JHEP04(2013)001
[4] LHCb collaboration, Measurement of b-hadron production fractions in 7 TeVpp collisions,
Phys. Rev. D 85 (2012) 032008 [arXiv:1111.2357] [INSPIRE].
[19] GEANT4 collaboration, S. Agostinelli et al., GEANT4: a simulation toolkit, Nucl. Instrum.
Meth. A 506 (2003) 250 [INSPIRE].
[20] M. Clemencic et al., The LHCb simulation application, Gauss: design, evolution and
experience, J. Phys. Conf. Ser. 331 (2011) 032023.
[21] L. Breiman, J.H. Friedman, R.A. Olshen and C.J. Stone, Classification and regression trees,
Wadsworth international group, Belmont U.S.A. (1984).
[22] BABAR collaboration, P. del Amo Sanchez et al., Dalitz plot analysis of Ds+ → K + K − π + ,
Phys. Rev. D 83 (2011) 052001 [arXiv:1011.4190] [INSPIRE].
[24] T. Skwarnicki, A study of the radiative cascade transitions between the Υ0 and Υ resonances,
Ph.D. thesis, Institute of Nuclear Physics, Krakow, Poland (1986) [INSPIRE].
– 12 –
JHEP04(2013)001
[23] CLEO collaboration, G. Bonvicini et al., Dalitz plot analysis of the D+ → K − π + π + decay,
Phys. Rev. D 78 (2008) 052001 [arXiv:0802.4214] [INSPIRE].
The LHCb collaboration
– 13 –
JHEP04(2013)001
R. Aaij38 , C. Abellan Beteta33,n , A. Adametz11 , B. Adeva34 , M. Adinolfi43 , C. Adrover6 ,
A. Affolder49 , Z. Ajaltouni5 , J. Albrecht35 , F. Alessio35 , M. Alexander48 , S. Ali38 , G. Alkhazov27 ,
P. Alvarez Cartelle34 , A.A. Alves Jr22,35 , S. Amato2 , Y. Amhis7 , L. Anderlini17,f , J. Anderson37 ,
R. Andreassen57 , R.B. Appleby51 , O. Aquines Gutierrez10 , F. Archilli18 , A. Artamonov 32 ,
M. Artuso53 , E. Aslanides6 , G. Auriemma22,m , S. Bachmann11 , J.J. Back45 , C. Baesso54 ,
V. Balagura28 , W. Baldini16 , R.J. Barlow51 , C. Barschel35 , S. Barsuk7 , W. Barter44 , A. Bates48 ,
Th. Bauer38 , A. Bay36 , J. Beddow48 , I. Bediaga1 , S. Belogurov28 , K. Belous32 , I. Belyaev28 ,
E. Ben-Haim8 , M. Benayoun8 , G. Bencivenni18 , S. Benson47 , J. Benton43 , A. Berezhnoy29 ,
R. Bernet37 , M.-O. Bettler44 , M. van Beuzekom38 , A. Bien11 , S. Bifani12 , T. Bird51 ,
A. Bizzeti17,h , P.M. Bjørnstad51 , T. Blake35 , F. Blanc36 , C. Blanks50 , J. Blouw11 , S. Blusk53 ,
A. Bobrov31 , V. Bocci22 , A. Bondar31 , N. Bondar27 , W. Bonivento15 , S. Borghi51 , A. Borgia53 ,
T.J.V. Bowcock49 , E. Bowen37 , C. Bozzi16 , T. Brambach9 , J. van den Brand39 , J. Bressieux36 ,
D. Brett51 , M. Britsch10 , T. Britton53 , N.H. Brook43 , H. Brown49 , A. Büchler-Germann37 ,
I. Burducea26 , A. Bursche37 , J. Buytaert35 , S. Cadeddu15 , O. Callot7 , M. Calvi20,j ,
M. Calvo Gomez33,n , A. Camboni33 , P. Campana18,35 , A. Carbone14,c , G. Carboni21,k ,
R. Cardinale19,i , A. Cardini15 , H. Carranza-Mejia47 , L. Carson50 , K. Carvalho Akiba2 , G. Casse49 ,
M. Cattaneo35 , Ch. Cauet9 , M. Charles52 , Ph. Charpentier35 , P. Chen3,36 , N. Chiapolini37 ,
M. Chrzaszcz 23 , K. Ciba35 , X. Cid Vidal34 , G. Ciezarek50 , P.E.L. Clarke47 , M. Clemencic35 ,
H.V. Cliff44 , J. Closier35 , C. Coca26 , V. Coco38 , J. Cogan6 , E. Cogneras5 , P. Collins35 ,
A. Comerma-Montells33 , A. Contu15 , A. Cook43 , M. Coombes43 , G. Corti35 , B. Couturier35 ,
G.A. Cowan36 , D. Craik45 , S. Cunliffe50 , R. Currie47 , C. D’Ambrosio35 , P. David8 ,
P.N.Y. David38 , I. De Bonis4 , K. De Bruyn38 , S. De Capua51 , M. De Cian37 , J.M. De Miranda1 ,
L. De Paula2 , W. De Silva57 , P. De Simone18 , D. Decamp4 , M. Deckenhoff9 , H. Degaudenzi36,35 ,
L. Del Buono8 , C. Deplano15 , D. Derkach14 , O. Deschamps5 , F. Dettori39 , A. Di Canto11 ,
J. Dickens44 , H. Dijkstra35 , P. Diniz Batista1 , M. Dogaru26 , F. Domingo Bonal33,n , S. Donleavy49 ,
F. Dordei11 , A. Dosil Suárez34 , D. Dossett45 , A. Dovbnya40 , F. Dupertuis36 , R. Dzhelyadin32 ,
A. Dziurda23 , A. Dzyuba27 , S. Easo46,35 , U. Egede50 , V. Egorychev28 , S. Eidelman31 ,
D. van Eijk38 , S. Eisenhardt47 , U. Eitschberger9 , R. Ekelhof9 , L. Eklund48 , I. El Rifai5 ,
Ch. Elsasser37 , D. Elsby42 , A. Falabella14,e , C. Färber11 , G. Fardell47 , C. Farinelli38 , S. Farry12 ,
V. Fave36 , D. Ferguson47 , V. Fernandez Albor34 , F. Ferreira Rodrigues1 , M. Ferro-Luzzi35 ,
S. Filippov30 , C. Fitzpatrick35 , M. Fontana10 , F. Fontanelli19,i , R. Forty35 , O. Francisco2 ,
M. Frank35 , C. Frei35 , M. Frosini17,f , S. Furcas20 , E. Furfaro21 , A. Gallas Torreira34 , D. Galli14,c ,
M. Gandelman2 , P. Gandini52 , Y. Gao3 , J. Garofoli53 , P. Garosi51 , J. Garra Tico44 , L. Garrido33 ,
C. Gaspar35 , R. Gauld52 , E. Gersabeck11 , M. Gersabeck51 , T. Gershon45,35 , Ph. Ghez4 ,
V. Gibson44 , V.V. Gligorov35 , C. Göbel54 , D. Golubkov28 , A. Golutvin50,28,35 , A. Gomes2 ,
H. Gordon52 , M. Grabalosa Gándara5 , R. Graciani Diaz33 , L.A. Granado Cardoso35 ,
E. Graugés33 , G. Graziani17 , A. Grecu26 , E. Greening52 , S. Gregson44 , O. Grünberg55 , B. Gui53 ,
E. Gushchin30 , Yu. Guz32 , T. Gys35 , C. Hadjivasiliou53 , G. Haefeli36 , C. Haen35 , S.C. Haines44 ,
S. Hall50 , T. Hampson43 , S. Hansmann-Menzemer11 , N. Harnew52 , S.T. Harnew43 , J. Harrison51 ,
P.F. Harrison45 , T. Hartmann55 , J. He7 , V. Heijne38 , K. Hennessy49 , P. Henrard5 ,
J.A. Hernando Morata34 , E. van Herwijnen35 , E. Hicks49 , D. Hill52 , M. Hoballah5 , C. Hombach51 ,
P. Hopchev4 , W. Hulsbergen38 , P. Hunt52 , T. Huse49 , N. Hussain52 , D. Hutchcroft49 , D. Hynds48 ,
V. Iakovenko41 , P. Ilten12 , J. Imong43 , R. Jacobsson35 , A. Jaeger11 , E. Jans38 , F. Jansen38 ,
P. Jaton36 , F. Jing3 , M. John52 , D. Johnson52 , C.R. Jones44 , B. Jost35 , M. Kaballo9 ,
S. Kandybei40 , M. Karacson35 , T.M. Karbach35 , I.R. Kenyon42 , U. Kerzel35 , T. Ketel39 ,
A. Keune36 , B. Khanji20 , O. Kochebina7 , I. Komarov36,29 , R.F. Koopman39 , P. Koppenburg38 ,
– 14 –
JHEP04(2013)001
M. Korolev29 , A. Kozlinskiy38 , L. Kravchuk30 , K. Kreplin11 , M. Kreps45 , G. Krocker11 ,
P. Krokovny31 , F. Kruse9 , M. Kucharczyk20,23,j , V. Kudryavtsev31 , T. Kvaratskheliya28,35 ,
V.N. La Thi36 , D. Lacarrere35 , G. Lafferty51 , A. Lai15 , D. Lambert47 , R.W. Lambert39 ,
E. Lanciotti35 , G. Lanfranchi18,35 , C. Langenbruch35 , T. Latham45 , C. Lazzeroni42 , R. Le Gac6 ,
J. van Leerdam38 , J.-P. Lees4 , R. Lefèvre5 , A. Leflat29,35 , J. Lefrançois7 , O. Leroy6 , Y. Li3 ,
L. Li Gioi5 , M. Liles49 , R. Lindner35 , C. Linn11 , B. Liu3 , G. Liu35 , J. von Loeben20 , J.H. Lopes2 ,
E. Lopez Asamar33 , N. Lopez-March36 , H. Lu3 , J. Luisier36 , H. Luo47 , A. Mac Raighne48 ,
F. Machefert7 , I.V. Machikhiliyan4,28 , F. Maciuc26 , O. Maev27,35 , S. Malde52 , G. Manca15,d ,
G. Mancinelli6 , N. Mangiafave44 , U. Marconi14 , R. Märki36 , J. Marks11 , G. Martellotti22 ,
A. Martens8 , L. Martin52 , A. Martı́n Sánchez7 , M. Martinelli38 , D. Martinez Santos39 ,
D. Martins Tostes2 , A. Massafferri1 , R. Matev35 , Z. Mathe35 , C. Matteuzzi20 , M. Matveev27 ,
E. Maurice6 , A. Mazurov16,30,35,e , J. McCarthy42 , R. McNulty12 , B. Meadows57,52 , F. Meier9 ,
M. Meissner11 , M. Merk38 , D.A. Milanes13 , M.-N. Minard4 , J. Molina Rodriguez54 , S. Monteil5 ,
D. Moran51 , P. Morawski23 , R. Mountain53 , I. Mous38 , F. Muheim47 , K. Müller37 , R. Muresan26 ,
B. Muryn24 , B. Muster36 , P. Naik43 , T. Nakada36 , R. Nandakumar46 , I. Nasteva1 , M. Needham47 ,
N. Neufeld35 , A.D. Nguyen36 , T.D. Nguyen36 , C. Nguyen-Mau36,o , M. Nicol7 , V. Niess5 , R. Niet9 ,
N. Nikitin29 , T. Nikodem11 , S. Nisar56 , A. Nomerotski52 , A. Novoselov32 ,
A. Oblakowska-Mucha24 , V. Obraztsov32 , S. Oggero38 , S. Ogilvy48 , O. Okhrimenko41 ,
R. Oldeman15,d,35 , M. Orlandea26 , J.M. Otalora Goicochea2 , P. Owen50 , B.K. Pal53 ,
A. Palano13,b , M. Palutan18 , J. Panman35 , A. Papanestis46 , M. Pappagallo48 , C. Parkes51 ,
C.J. Parkinson50 , G. Passaleva17 , G.D. Patel49 , M. Patel50 , G.N. Patrick46 , C. Patrignani19,i ,
C. Pavel-Nicorescu26 , A. Pazos Alvarez34 , A. Pellegrino38 , G. Penso22,l , M. Pepe Altarelli35 ,
S. Perazzini14,c , D.L. Perego20,j , E. Perez Trigo34 , A. Pérez-Calero Yzquierdo33 , P. Perret5 ,
M. Perrin-Terrin6 , G. Pessina20 , K. Petridis50 , A. Petrolini19,i , A. Phan53 , E. Picatoste Olloqui33 ,
B. Pietrzyk4 , T. Pilař45 , D. Pinci22 , S. Playfer47 , M. Plo Casasus34 , F. Polci8 , G. Polok23 ,
A. Poluektov45,31 , E. Polycarpo2 , D. Popov10 , B. Popovici26 , C. Potterat33 , A. Powell52 ,
J. Prisciandaro36 , V. Pugatch41 , A. Puig Navarro36 , W. Qian4 , J.H. Rademacker43 ,
B. Rakotomiaramanana36 , M.S. Rangel2 , I. Raniuk40 , N. Rauschmayr35 , G. Raven39 ,
S. Redford52 , M.M. Reid45 , A.C. dos Reis1 , S. Ricciardi46 , A. Richards50 , K. Rinnert49 ,
V. Rives Molina33 , D.A. Roa Romero5 , P. Robbe7 , E. Rodrigues51 , P. Rodriguez Perez34 ,
G.J. Rogers44 , S. Roiser35 , V. Romanovsky32 , A. Romero Vidal34 , J. Rouvinet36 , T. Ruf35 ,
H. Ruiz33 , G. Sabatino22,k , J.J. Saborido Silva34 , N. Sagidova27 , P. Sail48 , B. Saitta15,d ,
C. Salzmann37 , B. Sanmartin Sedes34 , M. Sannino19,i , R. Santacesaria22 , C. Santamarina Rios34 ,
E. Santovetti21,k , M. Sapunov6 , A. Sarti18,l , C. Satriano22,m , A. Satta21 , M. Savrie16,e ,
D. Savrina28,29 , P. Schaack50 , M. Schiller39 , H. Schindler35 , S. Schleich9 , M. Schlupp9 ,
M. Schmelling10 , B. Schmidt35 , O. Schneider36 , A. Schopper35 , M.-H. Schune7 , R. Schwemmer35 ,
B. Sciascia18 , A. Sciubba18,l , M. Seco34 , A. Semennikov28 , K. Senderowska24 , I. Sepp50 ,
N. Serra37 , J. Serrano6 , P. Seyfert11 , M. Shapkin32 , I. Shapoval40,35 , P. Shatalov28 ,
Y. Shcheglov27 , T. Shears49,35 , L. Shekhtman31 , O. Shevchenko40 , V. Shevchenko28 , A. Shires50 ,
R. Silva Coutinho45 , T. Skwarnicki53 , N.A. Smith49 , E. Smith52,46 , M. Smith51 , K. Sobczak5 ,
M.D. Sokoloff57 , F.J.P. Soler48 , F. Soomro18,35 , D. Souza43 , B. Souza De Paula2 , B. Spaan9 ,
A. Sparkes47 , P. Spradlin48 , F. Stagni35 , S. Stahl11 , O. Steinkamp37 , S. Stoica26 , S. Stone53 ,
B. Storaci37 , M. Straticiuc26 , U. Straumann37 , V.K. Subbiah35 , S. Swientek9 , V. Syropoulos39 ,
M. Szczekowski25 , P. Szczypka36,35 , T. Szumlak24 , S. T’Jampens4 , M. Teklishyn7 ,
E. Teodorescu26 , F. Teubert35 , C. Thomas52 , E. Thomas35 , J. van Tilburg11 , V. Tisserand4 ,
M. Tobin37 , S. Tolk39 , D. Tonelli35 , S. Topp-Joergensen52 , N. Torr52 , E. Tournefier4,50 ,
S. Tourneur36 , M.T. Tran36 , M. Tresch37 , A. Tsaregorodtsev6 , P. Tsopelas38 , N. Tuning38 ,
M. Ubeda Garcia35 , A. Ukleja25 , D. Urner51 , U. Uwer11 , V. Vagnoni14 , G. Valenti14 ,
R. Vazquez Gomez33 , P. Vazquez Regueiro34 , S. Vecchi16 , J.J. Velthuis43 , M. Veltri17,g ,
G. Veneziano36 , M. Vesterinen35 , B. Viaud7 , D. Vieira2 , X. Vilasis-Cardona33,n , A. Vollhardt37 ,
D. Volyanskyy10 , D. Voong43 , A. Vorobyev27 , V. Vorobyev31 , C. Voß55 , H. Voss10 , R. Waldi55 ,
R. Wallace12 , S. Wandernoth11 , J. Wang53 , D.R. Ward44 , N.K. Watson42 , A.D. Webber51 ,
D. Websdale50 , M. Whitehead45 , J. Wicht35 , D. Wiedner11 , L. Wiggers38 , G. Wilkinson52 ,
M.P. Williams45,46 , M. Williams50,p , F.F. Wilson46 , J. Wishahi9 , M. Witek23 , W. Witzeling35 ,
S.A. Wotton44 , S. Wright44 , S. Wu3 , K. Wyllie35 , Y. Xie47,35 , F. Xing52 , Z. Xing53 , Z. Yang3 ,
R. Young47 , X. Yuan3 , O. Yushchenko32 , M. Zangoli14 , M. Zavertyaev10,a , F. Zhang3 , L. Zhang53 ,
W.C. Zhang12 , Y. Zhang3 , A. Zhelezov11 , A. Zhokhov28 , L. Zhong3 , A. Zvyagin35
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
Centro Brasileiro de Pesquisas Fı́sicas (CBPF), Rio de Janeiro, Brazil
Universidade Federal do Rio de Janeiro (UFRJ), Rio de Janeiro, Brazil
Center for High Energy Physics, Tsinghua University, Beijing, China
LAPP, Université de Savoie, CNRS/IN2P3, Annecy-Le-Vieux, France
Clermont Université, Université Blaise Pascal, CNRS/IN2P3, LPC, Clermont-Ferrand, France
CPPM, Aix-Marseille Université, CNRS/IN2P3, Marseille, France
LAL, Université Paris-Sud, CNRS/IN2P3, Orsay, France
LPNHE, Université Pierre et Marie Curie, Université Paris Diderot, CNRS/IN2P3, Paris, France
Fakultät Physik, Technische Universität Dortmund, Dortmund, Germany
Max-Planck-Institut für Kernphysik (MPIK), Heidelberg, Germany
Physikalisches Institut, Ruprecht-Karls-Universität Heidelberg, Heidelberg, Germany
School of Physics, University College Dublin, Dublin, Ireland
Sezione INFN di Bari, Bari, Italy
Sezione INFN di Bologna, Bologna, Italy
Sezione INFN di Cagliari, Cagliari, Italy
Sezione INFN di Ferrara, Ferrara, Italy
Sezione INFN di Firenze, Firenze, Italy
Laboratori Nazionali dell’INFN di Frascati, Frascati, Italy
Sezione INFN di Genova, Genova, Italy
Sezione INFN di Milano Bicocca, Milano, Italy
Sezione INFN di Roma Tor Vergata, Roma, Italy
Sezione INFN di Roma La Sapienza, Roma, Italy
Henryk Niewodniczanski Institute of Nuclear Physics Polish Academy of Sciences, Kraków, Poland
AGH University of Science and Technology, Kraków, Poland
National Center for Nuclear Research (NCBJ), Warsaw, Poland
Horia Hulubei National Institute of Physics and Nuclear Engineering, Bucharest-Magurele,
Romania
Petersburg Nuclear Physics Institute (PNPI), Gatchina, Russia
Institute of Theoretical and Experimental Physics (ITEP), Moscow, Russia
Institute of Nuclear Physics, Moscow State University (SINP MSU), Moscow, Russia
Institute for Nuclear Research of the Russian Academy of Sciences (INR RAN), Moscow, Russia
Budker Institute of Nuclear Physics (SB RAS) and Novosibirsk State University, Novosibirsk,
Russia
Institute for High Energy Physics (IHEP), Protvino, Russia
Universitat de Barcelona, Barcelona, Spain
Universidad de Santiago de Compostela, Santiago de Compostela, Spain
European Organization for Nuclear Research (CERN), Geneva, Switzerland
Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
Physik-Institut, Universität Zürich, Zürich, Switzerland
Nikhef National Institute for Subatomic Physics, Amsterdam, The Netherlands
Nikhef National Institute for Subatomic Physics and VU University Amsterdam, Amsterdam, The
Netherlands
– 15 –
JHEP04(2013)001
1
40
41
42
43
44
45
46
47
48
49
50
52
53
54
55
56
57
a
b
c
d
e
f
g
h
i
j
k
l
m
n
o
p
P.N. Lebedev Physical Institute, Russian Academy of Science (LPI RAS), Moscow, Russia
Università di Bari, Bari, Italy
Università di Bologna, Bologna, Italy
Università di Cagliari, Cagliari, Italy
Università di Ferrara, Ferrara, Italy
Università di Firenze, Firenze, Italy
Università di Urbino, Urbino, Italy
Università di Modena e Reggio Emilia, Modena, Italy
Università di Genova, Genova, Italy
Università di Milano Bicocca, Milano, Italy
Università di Roma Tor Vergata, Roma, Italy
Università di Roma La Sapienza, Roma, Italy
Università della Basilicata, Potenza, Italy
LIFAELS, La Salle, Universitat Ramon Llull, Barcelona, Spain
Hanoi University of Science, Hanoi, Viet Nam
Massachusetts Institute of Technology, Cambridge, MA, United States
– 16 –
JHEP04(2013)001
51
NSC Kharkiv Institute of Physics and Technology (NSC KIPT), Kharkiv, Ukraine
Institute for Nuclear Research of the National Academy of Sciences (KINR), Kyiv, Ukraine
University of Birmingham, Birmingham, United Kingdom
H.H. Wills Physics Laboratory, University of Bristol, Bristol, United Kingdom
Cavendish Laboratory, University of Cambridge, Cambridge, United Kingdom
Department of Physics, University of Warwick, Coventry, United Kingdom
STFC Rutherford Appleton Laboratory, Didcot, United Kingdom
School of Physics and Astronomy, University of Edinburgh, Edinburgh, United Kingdom
School of Physics and Astronomy, University of Glasgow, Glasgow, United Kingdom
Oliver Lodge Laboratory, University of Liverpool, Liverpool, United Kingdom
Imperial College London, London, United Kingdom
School of Physics and Astronomy, University of Manchester, Manchester, United Kingdom
Department of Physics, University of Oxford, Oxford, United Kingdom
Syracuse University, Syracuse, NY, United States
Pontifı́cia Universidade Católica do Rio de Janeiro (PUC-Rio), Rio de Janeiro, Brazil, associated
to2
Institut für Physik, Universität Rostock, Rostock, Germany, associated to11
Institute of Information Technology, COMSATS, Lahore, Pakistan, associated to53
University of Cincinnati, Cincinnati, OH, United States, associated to53