Combined search for the quarks of a sequential fourth generation S. Chatrchyan et al.* (CMS Collaboration) (Received 5 September 2012; published 12 December 2012) Results are presented from a search for a fourth generation of quarks produced singly or in pairs in a data set corresponding to an integrated luminosity of 5 fb�1 recorded by the CMS experiment at the LHC in 2011. A novel strategy has been developed for a combined search for quarks of the up and down type in decay channels with at least one isolated muon or electron. Limits on the mass of the fourth-generation quarks and the relevant Cabibbo-Kobayashi-Maskawa matrix elements are derived in the context of a simple extension of the standard model with a sequential fourth generation of fermions. The existence of mass-degenerate fourth-generation quarks with masses below 685 GeV is excluded at 95% confidence level for minimal off-diagonal mixing between the third- and the fourth-generation quarks. With a mass difference of 25 GeV between the quark masses, the obtained limit on the masses of the fourth-generation quarks shifts by about �20 GeV. These results significantly reduce the allowed parameter space for a fourth generation of fermions. DOI: 10.1103/PhysRevD.86.112003 PACS numbers: 13.85.Rm, 12.60.�i, 14.65.Jk I. INTRODUCTION The existence of three generations of fermions has been firmly established experimentally [1]. The possibility of a fourth generation of fermions has not been excluded, although it is strongly constrained by precision measure- ments of electroweak observables. These observables are mainly influenced by the mass differences between the fourth-generation leptons or quarks. In particular, scenar- ios with a mass difference between the fourth-generation quarks smaller than the mass of theW boson are preferred, and even fourth-generation quarks with degenerate masses are allowed [2,3]. A new generation of fermions requires not only the existence of two additional quarks and two additional leptons, but also an extension of the Cabibbo-Kobayashi- Maskawa (CKM) [4,5] and Pontecorvo-Maki-Nakagawa- Sakata [6,7] matrices. New CKM (quark mixing) and Pontecorvo-Maki-Nakagawa-Sakata (lepton mixing) matrix elements are constrained by the requirement of consistency with electroweak precision measurements [8]. Previous searches at hadron colliders have considered either pair production or single production of one of the fourth-generation quarks [9–15]. The most stringent limits exclude the existence of a down-type (up-type) fourth- generation quark with a mass below 611 (570) GeV [14,15]. These limits on the quark mass values enter a region where the coupling of fourth-generation quarks to the Higgs field becomes large and perturbative calculations for the weak interaction start to fail, assuming the absence of other phenomena beyond the standard model [16]. To increase the sensitivity and to use a consistent approach while searching for a new generation of quarks, we have developed a simultaneous search for the up-type and down- type fourth-generation quarks, based on both the electro- weak and strong production mechanisms. If a fourth generation of quarks exists, their production cross sections and decay branching fractions will be gov- erned by an extended 4� 4 CKM matrix, V4�4 CKM, in which we denote the up- and down-type fourth-generation quarks as t0, and b0, respectively. For simplicity, we assume a model with one free parameter, A, where 0 � A � 1: V4�4 CKM ¼ Vud Vus Vub Vub0 Vcd Vcs Vcb Vcb0 Vtd Vts Vtb Vtb0 Vt0d Vt0s Vt0b Vt0b0 0 BBBBB@ 1 CCCCCA ¼ Oð1Þ Oð0Þ Oð0Þ 0 Oð0Þ Oð1Þ Oð0Þ 0 Oð0Þ Oð0Þ ffiffiffiffi A p ffiffiffiffiffiffiffiffiffiffiffiffiffi 1� A p 0 0 � ffiffiffiffiffiffiffiffiffiffiffiffiffi 1� A p ffiffiffiffi A p 0 BBBBB@ 1 CCCCCA : The complex phases are not shown for clarity. Within this model, mixing is allowed only between the third and the fourth generations. This is a reasonable assumption since the mixing between the third and the first two generations is observed to be small [17]. However, the limits presented in this paper would be too stringent if there is a fourth generation that mixes only with the first two generations, or the size of the mixing with the third generation is about the same as the mixing with the first two generations. With this search,we set limits on themasses of the fourth- generation quarks as a function of A. Since ffiffiffiffi A p ¼jVtbj, the lower limit of jVtbj> 0:81 from the single-top production *Full author list given at the end of the article. Published by the American Physical Society under the terms of the Creative Commons Attribution 3.0 License. Further distri- bution of this work must maintain attribution to the author(s) and the published article’s title, journal citation, and DOI. PHYSICAL REVIEW D 86, 112003 (2012) 1550-7998=2012=86(11)=112003(20) 112003-1 � 2012 CERN, for the CMS Collaboration http://dx.doi.org/10.1103/PhysRevD.86.112003 http://creativecommons.org/licenses/by/3.0/ cross section measurements [18] translates into a lower limit on the mixing between the third- and fourth- generation quarks in our model of A > 0:66. Using the data collected from ffiffiffi s p ¼ 7 TeV proton- proton collisions at the Large Hadron Collider (LHC), we search for fourth-generation quarks that are produced in pairs, namely b0 �b0 and t0 �t0, or through electroweak produc- tion, in particular tb0, t0b, and t0b0, where the charges are omitted in the notation. While the cross sections of the pair production processes do not depend on the value of A, the production cross sections of the tb0 and t0b processes depend linearly on ð1� AÞ, and the single-top and t0b0 cross sections on A. We assume the t0 and b0 masses to be degenerate within 25 GeV. In the case they are degenerate, they will decay in 100% of the cases to the third-generation quarks, since the decay of one fourth-generation quark to the other is kine- matically not allowed. However, even for nonzero mass differences, the branching fractions of the t0 ! bW and the b0 ! tW ! ðbWÞW decays are close to 100%, provided that the mass difference is small [19]. For instance, for a mass splitting of 25 GeV, and for Vt0b0 ¼ 0:005 (which would correspond to A ¼ 0:99975 in our model), less than 5% of the decays will be b0 ! t0W� (in the case mt0 mb0). For larger values of Vt0b, the branching fractions of b0 ! t0W� (or t0 ! b0W�) decrease even further. Therefore, the decay chains remain unchanged as long as the mass splitting is relatively small. We expect the following final states: (i) t0b ! bWb; (ii) t0 �t0 ! bWbW; (iii) b0t ! tWbW ! bWWbW; (iv) b0t0 ! tWbW ! bWWbW; (v) b0 �b0 ! tWtW ! bWWbWW. These decay chains imply that two jets from b quarks and one to four W bosons are expected in the final state for fourth-generation quarks produced both singly and in pairs. The W bosons decay to either hadronic or leptonic final states. Events with either one isolated lepton (muon or electron) or two same-sign dileptons or three leptons are selected. The different production processes are classified according to the number of observed W bosons. II. THE COMPACT MUON SOLENOID DETECTOR The central feature of the Compact Muon Solenoid (CMS) detector is a superconducting solenoid, 13 m in length and 6 m in internal diameter, providing an axial magnetic field of 3.8 T. The inside of the solenoid is equipped with various particle detection systems. Charged particle trajectories are measured by a silicon pixel and strip tracker, covering 0<�< 2� in azimuth and j�j< 2:5, where the pseudorapidity � is defined as � ln½tanð�=2Þ�, and � is the polar angle of the trajectory with respect to the anticlockwise-beam direction. A crystal electromagnetic calorimeter and a brass/scintillator hadron calorimeter surround the tracking volume and provide high-resolution energy and direction measurements of electrons, photons, and hadronic jets. Muons are measured in gas-ionization detectors embedded in the steel return yoke outside the solenoid. The CMS detector also has extensive forward calorimetry covering up to j�j< 5. The detector is nearly hermetic, allowing for energy bal- ance measurements in the plane transverse to the beam directions. A two-tier trigger system selects the most inter- esting proton collision events for use in physics analysis. A more detailed description of the CMS detector can be found elsewhere [20]. III. EVENT SELECTION AND SIMULATION The search for the fourth-generation quarks is performed using the ffiffiffi s p ¼ 7 TeV proton-proton collisions recorded by the CMS experiment at the LHC. We have analyzed the full data set collected in 2011 corresponding to an inte- grated luminosity of ð5:0� 0:1Þ fb�1. Events are selected with a trigger requiring an isolated muon or electron, where the latter is accompanied by at least one jet identi- fied as a b jet. The muon system, the calorimetry, and the tracker are used for the particle-flow event reconstruction [21]. Jets are reconstructed using the anti-kT algorithm [22] with a size parameter of 0.5. Events are further selected with at least one high-quality isolated muon or elec- tron with a transverse momentum (pT) exceeding 40 GeV in the acceptance range j�j< 2:1 for muons and j�j<2:5 for electrons. The relative isolation, Irel, is calculated from the other particle-flow particles within a cone of �R ¼ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffið��Þ2 þ ð��Þ2p < 0:4 around the axis of the lepton. It is defined as Irel ¼ ðEcharged T þ Ephoton T þ Eneutral T Þ=pT, where E charged T and E photon T are the transverse energies deposited by charged hadrons and photons, respectively, and Eneutral T is the transverse energy deposited by neutral particles other than photons. We identify muons and elec- trons as isolated when Irel < 0:125 and Irel < 0:1, respec- tively. The requirement on the relative isolation for electrons is tighter than for muons because the back- grounds for electrons are higher than for muons. Electron candidates in the transition region between electromag- netic calorimeter barrel and end cap (1:44< j�j< 1:57) are excluded because the reconstruction of an electron object in this region is not optimal. We require a missing transverse momentum 6ET of at least 40 GeV. The 6ET is calculated as the absolute value of the vector sum of the pT of all reconstructed objects. Jets are required to have a pT > 30 GeV. The jet energies are corrected to establish a uniform response of the calorimeter in � and a calibrated absolute response in pT. Furthermore, a correction is applied to take into account the energy clustered in jets due to additional proton interactions in the same bunch crossing. S. CHATRCHYAN et al. PHYSICAL REVIEW D 86, 112003 (2012) 112003-2 The observed data are compared to simulated data gen- erated with POWHEG 301 [23,24] for the single-top process, PYTHIA 6.4.22 [25] for the diboson processes, and MADGRAPH 5.1.1 [26] for the signal and other standard model processes. The POWHEG and MADGRAPH generators are interfaced with PYTHIA for the decay of the particles as well as the hadronization and the implementation of a CMS custom underlying event tuning (tune Z2) [27]. The match- ing of the matrix-element partons to the parton showers is obtained using the MLM matching algorithm [28]. The CTEQ6L1 leading-order (LO) parton distributions are used in the event generation [29]. The generated events are passed through the CMS detector simulation based on GEANT4 [30], and then processed by the same reconstruc- tion software as the collision data. The simulated events are reweighted to match the observed distribution of the num- ber of simultaneous proton interactions. For the full data set collected in 2011, we observe on average about nine interactions in each event. We smear the jet energies in the simulation to match the resolutions measured with data [31]. At least one of the jets within the tracker acceptance (j�j< 2:4) needs to be identified as a b jet. For the b-jet identification, we require the signed impact parameter significance of the third track in the jet (sorted by decreas- ing significance) to be larger than a value chosen such that the probability for a light quark jet to be misidentified as a b jet is about 1%. We apply scale factors measured from data to the simulated events to take into account the different b-jet efficiency and the different probability that a light quark or gluon is identified as a b jet in data and simulation [32]. The top-quark pair as well as the W and Z production cross section values used in the analysis correspond to the measured values from CMS [33,34]. We use the predicted cross section values for the single-top, t�tþW, t�tþ Z, and same-signWW processes [35–38]. The cross section values for the diboson production are obtained with the MCFM next-to-leading-order parton-level integrator [39,40]. For the pair-production of the fourth-generation quarks, we use the approximate next-to-next-to-leading-order cross section values from Ref. [41]. For the electroweak production processes mentioned above, we rescale the next-to-leading-order cross sections at 14 TeV [42] to 7 TeV using a scale factor defined as the ratio of the LO cross section at 7 TeVand the LO cross section at 14 TeVas obtained by the MADGRAPH event generator. The resulting production cross sections are maximal, hence assuming jVtb0 j ¼ jVt0bj ¼ jVt0b0 j ¼ 1, and are rescaled according to the value of A. IV. EVENT CLASSIFICATION Different channels are defined according to the number ofW bosons in the final state. Given that the t0 decay mode is the same as the top-quark decay mode, the t0b and t0 �t0 processes will yield signatures that are very similar to, respectively, the single-top and t�t processes in the standard model. We select these processes through the single-lepton decay channel. In the signal final states that contain a b0 quark, we expect three or fourW bosons. If two or more of theseW bosons decay to leptons, we may have events with two leptons of the same charge or with three charged leptons. Although the branching fraction of these decays is small compared to that of other decay channels, these final states are very interesting because of the low back- ground that is expected from standard model processes. A. The single-electron and single-muon decay channels On top of the aforementioned event selection criteria, we veto events with additional electrons or muons with Irel < 0:2 and pT > 10 GeV for muons and pT > 15 GeV for electrons. We divide the selected single-lepton events into different subsamples according to the signal final states. Therefore, we define a procedure to count the num- ber of W-boson candidates. Each event has at least one W boson that decays to leptons, consistent with the require- ments of an isolated lepton and a large missing transverse momentum from the neutrino, which escapes detection. The decays ofW bosons to q �q final states are reconstructed with the following procedure. For each event, we have a collection of selected jets used as input for the reconstruc- tion of the W-boson candidates. The one or two jets that are identified as b jets are removed from the collection. W-boson candidates are constructed from all possible pairs of the remaining jets in the collection. We use both the expected mass, mfit W ¼ 84:3 GeV, and the width, �fit mW ¼ 9:6 GeV, from a Gaussian fit to the reconstructed mass distribution of jet pairs from the decay of a W boson in simulated t�t events. The W-boson candidate with a mass that matches the value ofmfit W best is chosen as aW boson if its mass is within a �1�fit mW window around mfit W . The jet pair that provided the hadronically decaying W boson is removed from the collection, and the procedure is repeated until no more candidates are found for W bosons decaying to jets. Different exclusive subsamples are defined accord- ing to the number of b jets (exactly one or at least two) and the number of W-boson candidates (one, two, three, and at least four). There are seven subsamples, because we do not consider the subsample with only one b jet and one W boson. The subsample with two b jets and one W boson is dominated by singly produced t0 events. In this subsample, we apply a veto for additional jets with a transverse momentum exceeding 30 GeV. Furthermore, since b �b background tends to have jets which are produced back- to-back with balanced pT , we remove this background by requiring ��ðj1; j2Þ< � 2 þ �ðpj1 T � pj2 T Þ=ðpj1 T þ pj2 T Þ. Table I summarizes the requirements that define the different single-lepton decay subsamples, after the criteria on the 6ET, and the lepton and jet pT and � are applied. Table II shows the observed and predicted event yields. After the selection criteria, the dominant background COMBINED SEARCH FOR THE QUARKS OF A . . . PHYSICAL REVIEW D 86, 112003 (2012) 112003-3 contributions result from the production of top-quark pairs, W þ jets, and single top. Other processes with very small contributions to the total background are Zþ jets and diboson production, and also top-quark pairs produced in association with aW or Z boson. The combined event yield of these processes is about 1% of the total standard-model contribution. The multijet background is found to be negli- gible in each of the subsamples. The reason is the require- ments of an isolated muon or electron with pT > 40 GeV, a missing transversemomentum of 40GeV, and at least one jet identified as a b jet. Data and simulation are found to agree within the combined statistic and systematic uncertainties. B. The same-sign dilepton and trilepton decay channels The transverse momentum of at least one of the leptons in the multilepton channel is required to be larger than 40 GeV, while the threshold is reduced to 20 GeV for additional leptons. Events with two muons or electrons with a mass within 10 GeV of the Z-boson mass are rejected to reduce the standard model background with Z bosons in the final state. We require at least four jets for the same-sign dilepton events. In the case of the trilepton events, the minimum number of required jets is reduced to two. Table III summarizes the event selection require- ments defining the same-sign dilepton and trilepton decay channels that are applied on top of the other requirements on the 6ET and lepton and jet pT and �. There are several contributions to the total standard- model background for the same-sign dilepton events. One of these contributions comes from events for which the charge of one of the leptons is misreconstructed, for instance in t�t events with two W bosons decaying into leptons. Second, there are events with one prompt lepton and one nonprompt lepton passing the isolation and iden- tification criteria. Finally, there is an irreducible contribu- tion from standard-model processes with two prompt leptons of the same sign; e.g. W�W�, WZ, ZZ, t�tþW, and t�tþ Z. Except forW�W�, these processes are also the main contributions to the total background for the trilepton subsample. The event yields for the irreducible component of the background for the same-sign dilepton channel and the total background in the case of the trilepton subsample are taken from the simulation. We obtain from the data the predicted number of background events for the first two contributions to the total background in the same-sign dilepton subsample. For the same-sign dilepton events with at least one electron, the background is estimated from control samples. We determine the charge misidentification rate for electrons using a double-isolated-electron trigger. We require two isolated electrons with the dielectron in- variant mass within 10 GeVof the Z-boson mass. We select TABLE II. Event yields in the single lepton channel. Uncertainties reflect the combined statistical and systematic uncertainties. The prediction for the signal is shown for two different values of A and for a fourth-generation-quark mass mq0 ¼ 550 GeV. 1b 2W 1b 3W 1b 4W 2b 1W 2b 2W 2b 3W 2b 4W t�tþ jets 5630� 410 230þ29 �26 3:0þ1:9 �1:3 819þ59 �62 2810� 240 85þ12�10 0:6þ0:8 �0:5 W þ jets 490� 180 8:0þ3:1 �3:0 0:3þ0:9 �0:3 150þ47 �46 37� 12 1:1þ1:0 �0:4 0:0þ0:8 �0:0 Zþ jets 36þ5 �6 1:0þ0:2 �0:1 0 7:1þ1:0 �0:6 2:8þ1:0 �0:3 0 0 Single top 346� 64 6:5þ1:6 �1:5 0:2þ0:3 �0:2 200� 34 110� 19 2:5þ0:7 �0:5 0:0þ0:1 �0:0 VV 15� 2 0:4þ0:3 �0:1 0:0þ0:1 �0:0 15� 2 1:8� 0:3 0:0þ0:1 �0:0 0:0þ0:1 �0:0 t�tV 28� 3 3:4� 0:5 0:1� 0:0 0:7� 0:2 15� 5 1:5þ0:3 �0:2 0 Total background 6550� 450 249þ29 �26 3:6þ2:1�1:3 1190þ83 �85 2970� 240 91þ12�10 0:6þ1:2 �0:5 Observed 7003 242 8 1357 3043 91 4 Signal (A ¼ 1) 55� 1 12� 1 0:9� 0:2 1:0þ0:2 �0:3 49� 2 8:1� 0:4 0:5� 0:2 Signal (A ¼ 0:8) 85� 2 14� 1 1:0� 0:2 69� 3 66� 2 9:2� 0:4 0:5� 0:2 TABLE III. Overview of the event selection requirements spe- cific to the same-sign dilepton and trilepton decay channels. Same-sign dilepton Trilepton ¼ 2 isolated leptons with same sign ¼ 3 isolated leptons �4 jets ðpT>30GeV;j�j<2:4Þ �2jets ðpT>30GeV; j�j<2:4Þ � 1b jet � 1b jet TABLE I. Overview of the event selection requirements defin- ing the different subsamples in the single-lepton decay channel. The single-lepton decay channel is divided in seven different subsamples according to the number of b jets and the number of W-boson candidates. Single-lepton decay channel 1W 2W 3W 4W ¼ 2 jets � 4 jets � 6 jets � 8 jets ¼ 2b jets either ¼ 1 or � 2b jets ��ðj1; j2Þ requirement 1W ! q �q 2W ! q �q 3W ! q �q S. CHATRCHYAN et al. PHYSICAL REVIEW D 86, 112003 (2012) 112003-4 events with 6ET < 20 GeV and a transverse mass MT ¼ffiffiffiffiffiffiffiffiffiffiffiffiffiffi 2p‘ T 6ET q ½1� cosð��ð‘; 6ETÞÞ� less than 25 GeV to sup- press background from top-quark and W þ jets events. We define the charge misidentification ratio R as the num- ber of events with two electrons of the same sign divided by twice the number of events with two electrons of opposite sign, i.e. R ¼ NSS=2NOS. We obtain 0.14% and 1.4% for barrel and end-cap electron candidates, respec- tively. After the full event selection is applied, with the exception of the electron sign requirement, we obtain a number of selected data events with two electrons and with an electron and a muon in the final state. The background with two electrons or with an electron and a muon with the same sign is obtained by taking the number of opposite- sign events and scaling it with R. The pT spectrum of the electrons in the control sample and the signal region is similar. Therefore, no correction is applied for the pT dependency of the charge misidentification ratio. Another important background contribution to the same- sign dilepton channel originates from jets being misidenti- fied as an electron or a muon (‘‘fake’’ leptons). Two collections of leptons, ‘‘loose’’ and ‘‘tight’’, are defined based on the isolation and identification criteria. Loose leptons are required to fulfill Irel < 0:2, in contrast with Irel < 0:125 ð0:1Þ for tight muons (electrons). Moreover, we require j�j< 2:5 and pT < 10 ð15Þ for loose muons (electrons). Additionally, several identification criteria, intended to ensure the consistency of the lepton track with the primary vertex, are relaxed. We require at least one loose electron or muon. Additionally, we require 6ET < 20 GeV and MT < 25 GeV to suppress background from top-quark and W þ jets events. Moreover, we veto events with leptons of the same flavor which have a dilepton mass within 20 GeVof the Z-boson mass. We count the number of loose and tight leptons with a pT below 35 GeV. The threshold on the pT is required to suppress contamination from W þ jets events, which would bias the estimation, because leptons produced in jets have typically a soft pT spectrum. The probability that a loose (L) lepton passes the tight (T) selection criteria is then given by the ratio �TL ¼ NT=NL. To estimate the number of events from the background source with a nonprompt lepton, we count the number of events in data that pass the event selection criteria with one lepton passing the tight selection criteria and a second lepton passing the loose, but not the tight, criteria. This yield is multiplied by �TLð1� �TLÞ to deter- mine the number of events with a nonprompt lepton in the analysis. The statistical uncertainty on the estimated num- ber of events is large because only a few events are selected with one tight and one loose, but not tight, lepton. The total number of expected background events for the same-sign dilepton and trilepton channels is given in Table IV. V. SETTING LOWER LIMITS ON THE FOURTH-GENERATION QUARK MASSES We have defined different subsamples according to the reconstructed final state. In each of the different subsam- ples, we reconstruct observables that are sensitive to the presence of the fourth-generation quarks. These observ- ables are used as input to a fit of the combined distributions for the standard-model (background-only) hypothesis and the signal-plus-background hypothesis. With the profile likelihood ratio as a test statistic, we calculate the 95% confidence level (CL) upper limits on the combined input cross section of the signal as a function of the V4�4 CKM parameter A and the mass of the fourth-generation quarks. A. Observables sensitive to the fourth-generation quark production The expected number of events is small in the subsamples with two leptons of the same sign, the trilepton subsample, and the two single-lepton subsamples with four W-boson candidates. As a consequence, the event counts in each of these subsamples are used as the observable. Table IV sum- marizes the event counts for the subsamples with two leptons of the same sign and the trilepton subsample. In the single-lepton subsamples with one or three W bosons, we use ST as the observable to discriminate between the standard model background and the fourth- generation signal, where ST is defined as the scalar sum of the transverse momenta of the reconstructed objects in the final state, namely: TABLE IV. The prediction for the total number of background events compared with the number of observed events in the same-sign dilepton and the trilepton subsamples. The numbers of expected signal events are also shown for two possible scenarios. Type 2 muons 2 electrons Electronþmuon Trilepton Irreducible background 0:77� 0:08 0:59� 0:08 1:10� 0:11 0:96� 0:12 Background from charge misid � � � 0:47� 0:08 0:71� 0:06 � � � Background from fake leptons 0:06� 0:06 0:30� 0:15 0:46� 0:17 � � � Total background 0:83� 0:11 1:36� 0:19 2:27� 0:22 0:96� 0:12 Observed 2 2 2 1 Signal (A ¼ 1, mq0 ¼ 550 GeV) 3:31� 0:15 2:03� 0:36 5:29� 0:19 3:37� 0:16 Signal (A ¼ 0:8, mq0 ¼ 550 GeV) 3:79� 0:15 2:29� 0:36 6:00� 0:19 3:65� 0:16 COMBINED SEARCH FOR THE QUARKS OF A . . . PHYSICAL REVIEW D 86, 112003 (2012) 112003-5 ST ¼ 6ET þ pl T þ pb T þ pj T þ XN i¼0 p Wi q �q T ; (1) where the sum runs over the number of reconstructed hadronically decaying W bosons; pl T is the pT of the lepton, pb T the pT of the b jet, pj T the pT of the second b jet or, if there is no additional jet identified as a b jet, the pT of the jet with the highest transverse momentum in the event that is not used in the W-boson reconstruction, and p Wi q �q T the pT of the ith reconstructed W boson decaying to jets. In general, the decay products of the fourth-generation quarks are expected to have higher transverse momenta compared to the standard-model background. This is shown in Fig. 1 for three of the subsamples. The dominant contribution to the selected signal events in the subsample with two b jets and one W boson would come from the t0b process. Almost no signal events are selected for A ¼ 1, because in that case, the production cross section of t0b is equal to zero. The subsamples with two W bosons are dominated by t�t events. In this case, we use two sensitive observables: ST and the mass of the hadronic bW system, (GeV)TS 200 300 400 500 600 700 800 900 E ve nt s / 5 0 G eV 100 200 300 400 500 = 7 TeVs at -1CMS, 5 fb lepton+jets, 2b 1W = 550 GeVq'm tt ν l→W Single top Other Signal A=1 (X 8) Signal A=0.8 (X 8) Observed Systematic Uncertainty (a) (GeV)TS 200 300 400 500 600 700 800 900 1000 E ve nt s / 1 00 G eV 20 40 60 80 100 = 7 TeVs at -1CMS, 5 fb lepton+jets, 1b 3W = 550 GeVq'm tt ν l→W Single top Other Signal A=1 (X 8) Signal A=0.8 (X 8) Observed Systematic Uncertainty (b) (GeV)TS 200 300 400 500 600 700 800 900 1000 E ve nt s / 5 0 G eV 100 200 300 400 500 600 700 = 7 TeVs at -1CMS, 5 fb lepton+jets, 2b 2W = 550 GeVq'm tt ν l→W Single top Other Signal A=1 (X 8) Signal A=0.8 (X 8) Observed Systematic Uncertainty (c) (GeV)bWm 0 100 200 300 400 500 600 700 E ve nt s / 1 0 G eV 1 10 210 310 = 7 TeVs at -1CMS, 5 fb lepton+jets, 2b 2W = 550 GeVq'm tt ν l→W Other Signal A=1 (X 8) Signal A=0.8 (X 8) Observed Systematic Uncertainty (d) FIG. 1 (color online). The ST distribution for the subsamples with two b jets and oneW boson (a), one b jet and threeW bosons (b), two b jets and two W bosons (c), and the mbW distribution for the subsample with two b jets and two W bosons (d). The data distributions of these observables are compared to their expectation from the simulation assuming the fitted nuisance parameters. The fitted values of the nuisance parameters represent the systematic shifts that are applied on the simulation to fit the data in the background-only hypothesis. As an illustration, the total uncertainty band is shown around the simulated expected distribution before taking into account the fitted values of the nuisance parameters. The expected distribution for a signal is shown for two different values of the V4�4 CKM parameter A and for b0 and t0 masses of 550 GeV. The cross section of the signal in the plots is scaled by a factor of eight for visibility. The last bin in all the histograms includes the overflow. We do not expect much signal for A ¼ 1 in (a), because the subsample with two b jets and one W boson is mainly sensitive to single t0-quark production. S. CHATRCHYAN et al. PHYSICAL REVIEW D 86, 112003 (2012) 112003-6 mbW . The latter observable is sensitive to the fourth- generation physics, because of the higher mass of a hypothetical fourth-generation t0 quark compared to the top-quark mass. To obtain a higher sensitivity with thembW observable, four jets need to be assigned to the quarks to reconstruct the final state t0 �t0 ! WbWb ! q �qb‘�‘b. Therefore, six observables with discriminating power between correct and wrong jet/quark assignments are com- bined with a likelihood ratio method. These observables are angles between the decay products, theW-boson mass, the transverse momentum of the top quark decaying to hadrons, and an observable related to the values of the b-jet identification variable for the jets. The jet/quark assignment with the largest value of the likelihood ratio is chosen. The mass of the bW system is then reconstructed from this chosen jet/quark assignment. The lower plots in Fig. 1 show the projections of the two-dimensional ST versus mbW distribution. An overview of the observables used in the fit for the presence of the fourth-generation quarks is presented in Table V. B. Fitting for the presence of fourth-generation quarks We construct a single histogram ‘‘template’’ that con- tains the information of the sensitive observables from all the subsamples. Different template distributions are made for the signal corresponding to the different values of A and the fourth-generation quark masses mq0 . The binning of the two-dimensional observable distribution in the single- lepton subsamples with two W bosons is defined using the following procedure. We use a binning in the dimension of mbW such that the top-quark pair background events are uniformly distributed over the bins. Second, the binning in the dimension of ST in each of the mbW bins is chosen to obtain uniformly distributed top-quark pair events also in this dimension. The templates of the sensitive observables are used as input to obtain the likelihoods for the background-only and the signal-plus-background hypotheses. Systematic uncer- tainties are taken into account by introducing nuisance parameters, which may affect the shape and the normal- ization of the templates. In a case where the systematic uncertainty alters the shape of the templates, template morphing [43,44] is used to interpolate linearly on a bin-by-bin basis between the nominal templates and sys- tematically shifted ones. The normalization of the templates is affected by the uncertainty in the integrated luminosity, the lepton effi- ciency, and the normalization of the background processes. The integrated luminosity is measured with a precision of 2.2% [45] and has the same normalization effect on all the templates. The uncertainty in the lepton efficiency is a combination of the uncertainties in the trigger, selection, and identification efficiencies, which amounts to 3% and 5% for muon and electron, respectively. For the uncertainty in the normalization of the background processes, we use the uncertainties in the production cross section of the various standard-model processes. The most important contributions that affect the normalization of the templates are the 12% [33] (30%) uncertainty for the top-quark pair (single-top) production cross section and a 50% uncer- tainty for the W production cross section because of the large fraction of selected events with jets from heavy-flavor quarks. For the multilepton channel, we take into account the uncertainties in the background estimation obtained from the data. We also include the uncertainties in the production cross sections of Z (5% [34]), WW (35%), WZ (42%), ZZ (27%), t�tþW (19%), t�tþ Z (28%), and W�W� (49%). The uncertainties in the normalization of diboson and top-quark pair production in association with a boson are taken from a comparison of the next-to-leading- order and the LO predictions. The largest systematic effects on the shape of the tem- plates originate from the jet energy corrections [31] and the scale factors between data and simulation for the b-jet efficiency and the probability that a light quark or gluon is identified as a b jet [32]. These effects are estimated by varying the nominal value by �1 standard deviation. The uncertainty in the jet energy resolution of about 10% has a relatively small effect on the expected limits. The same is true for the uncertainty in the modeling of multiple inter- actions in the same beam crossing. The latter effect is evaluated by varying the average number of interactions in the simulation by 8%. The probability density functions of the background- only and the signal-plus-background hypotheses are fitted to the data to fix the nuisance parameters in both models. In the signal-plus-background model, an additional variable, defined as the cross section for the fourth-generation signal obtained by combining the separate search channels, is included. In the combined cross section variable, the rela- tive fraction of each fourth-generation signal process is fixed according to the probed model parameters ðA;mq0 Þ. Using a Gaussian approximation for the probability density function of the test statistic, we determine the 95% CL expected and observed limits on the combined cross sec- tion variable using the CLs criterion [46–48]. We exclude the point ðA;mq0 Þ at the 95% CL if the upper limit on the combined cross section variable is smaller than its TABLE V. Overview of the observables used in the limit calculation. Subsample Observable Single-lepton 1W ST Single-lepton 2W ST and mbW Single-lepton 3W ST Single-lepton 4W Event yield Same-sign dilepton Event yield Trilepton Event yield COMBINED SEARCH FOR THE QUARKS OF A . . . PHYSICAL REVIEW D 86, 112003 (2012) 112003-7 predicted value within the fourth-generation model. The procedure is repeated for each value of A and mq0 . C. Results and discussion We use the CLs procedure to calculate the combined limit for the single-muon, single-electron, same-sign di- lepton, and trilepton channels. When the value of the V4�4 CKM parameter A approaches unity, the standard model single-top and the t0b0 processes reach their maximal values for the production cross section. When the value of A decreases, the cross section of these processes decreases linearly with A. At the same time, the expected cross section of the t0b and tb0 processes increases with (1� A) and is equal to zero for A ¼ 1. Therefore, the t0b and tb0 processes are expected to enhance the sensitivity for fourth-generation quarks when the parameter A decreases. This is visible in the upper part of Fig. 2 where both the expected and observed limits on mq0 are more stringent for smaller values of A. For instance, the limit on the fourth-generation quark masses increases by 70 GeV for A ¼ 0:9 compared to the value of the limit for A 1. While the t0b and tb0 processes do not contribute for A 1, the inclusion of the t0b0 process results in a more stringent limit (a difference of about 30 GeV) compared to when this process is not taken into account. The existence of fourth-generation quarks with degen- erate masses is excluded for all parameter values below the line using the assumed model of the V4�4 CKM matrix. In particular, fourth-generation quarks with a degenerate mass below 685 GeV are excluded at the 95% CL for a parameter value of A 1. It is worth noting that no limits can be set for A exactly equal to unity (A ¼ 1), because in this special case, the fourth-generation quarks would be stable in the assumed model. The analysis is, however, valid for values of A extremely close to unity. The distance between the primary vertex and the decay vertex of the fourth-generation quarks is less than 1 mm for 1� A > 2� 10�14, a number obtained using the LO formula for the decay width of the top quark in which the top-quark mass is replaced with a fourth-generation-quark mass of 600 GeV. Up to now, the masses of the fourth-generation quarks were assumed to be degenerate. However, if a fourth generation of chiral quarks exists, this is not necessarily the case. Therefore, it is interesting to study how the limit would change for nondegenerate quark masses. If we assume nondegenerate masses, another decay channel for the fourth-generation quarks is possible. Namely, the branching fraction for the decay of t0 ðb0Þ into b0 ðt0Þ, and an off-shell W boson becomes nonzero. For values of the mass splitting up to about 25 GeV, this branching fraction is small as noted in the introduction. We assume a mass splitting of 25 GeVand unchanged branching fractions for the t0 and b0 decays. The sensitivity of the analysis increases or decreases depending on the specific values of the masses and hence the production cross sections of the fourth-generation quarks. The effect of the mass dif- ference between the fourth-generation quarks on the ex- clusion limit is shown in the bottom plot of Fig. 2 for a V4�4 CKM parameter A 1. For instance, in case mt0 ¼ mb0 þ 25 GeV (mt0 ¼ mb0 � 25 GeV), the limit on mt0 increases aboutþ20ð�20Þ GeV with respect to the degenerate-mass 2 tb' = 1 - V2 t'b = 1 - V2 t'b' = V2 tbA = V 0.8 0.85 0.9 0.95 1 ( G eV ) b' = m t' m 600 650 700 750 = 7 TeVs at -1CMS, 5 fb obs. limit 95% CL exp. limit 95% CL σ 1 ± σ 2 ± (GeV)b' - mt'm -20 -10 0 10 20 ( G eV ) t' m 500 550 600 650 700 750 = 7 TeVs at -1CMS, 5fb obs. limit 95% CL (A=1) exp. limit 95% CL (A=1) σ 1 ± σ 2 ± FIG. 2 (color online). Top: Exclusion limit on mt0 ¼ mb0 as a function of the V4�4 CKM parameter A. The parameter values below the solid line are excluded at 95% CL. The inner (outer) band indicates the 68% (95%) confidence interval around the expected limit. The slope indicates the sensitivity of the analysis to the t0b and tb0 processes. Bottom: For a V4�4 CKM parameter value A 1, the exclusion limit on mt0 versus mt0 �mb0 is shown. The exclusion limit is calculated for mass differences up to 25 GeV. The existence of up-type fourth-generation quarks with mass values below the observed limit are excluded at the 95% CL. S. CHATRCHYAN et al. PHYSICAL REVIEW D 86, 112003 (2012) 112003-8 case. To obtain this limit, we do not take into account the electroweak t0b0 process, which results in more conserva- tive exclusion limits. In particular, one observes that quarks with degenerate masses below about 655 GeVare excluded at the 95% CL compared to 685 GeV when the t0b0 process is included. VI. SUMMARY Results from a search for a fourth generation of quarks have been presented. A simple model for a unitary CKM matrix has been defined based on a single parameter A ¼ jVtbj2 ¼ jVt0b0 j2. Degenerate masses have been assumed for the fourth-generation quarks, hence mt0 ¼ mb0 . The information is combined from different subsam- ples corresponding to different final states with at least one electron or muon. Observables have been constructed in each of the subsamples and used to differentiate between the standard-model background and the processes with fourth-generation quarks. With this strategy, the search for singly and pair-produced t0 and b0 quarks has been combined in a coherent way into a single analysis. Model-dependent limits are derived on the mass of the quarks and the V4�4 CKM matrix element A. The existence of fourth-generation quarks with masses below 685 GeV is excluded at 95% confidence level for minimal off-diagonal mixing between the third- and the fourth-generation quarks. A nonzero cross section for the single fourth- generation quark production processes, corresponding to a value of the V4�4 CKM parameter A < 1, gives rise to a more stringent limit. When a mass difference of 25 GeV is assumed between t0 and b0 quarks, the limit on mt0 shifts by about þ20ð�20Þ GeV for mt0 ¼ mb0 þ 25 GeV (mt0 ¼ mb0 � 25 GeV). These results significantly reduce the allowed parameter space for a fourth generation of fermi- ons and raise the lower limits on the masses of the fourth generation quarks to the region where nonperturbative effects of the weak interactions are important. ACKNOWLEDGMENTS We congratulate our colleagues in the CERN accelerator departments for the excellent performance of the LHC machine. We thank the technical and administrative staff at CERN and other CMS institutes, and acknowledge sup- port from BMWF and FWF (Austria); FNRS and FWO (Belgium); CNPq, CAPES, FAPERJ, and FAPESP (Brazil); MES (Bulgaria); CERN; CAS, MoST, and NSFC (China); COLCIENCIAS (Colombia); MSES (Croatia); RPF (Cyprus); MoER, SF0690030s09 and ERDF (Estonia); Academy of Finland, MEC, and HIP (Finland); CEA and CNRS/IN2P3 (France); BMBF,DFG, andHGF (Germany); GSRT (Greece); OTKA and NKTH (Hungary); DAE and DST (India); IPM (Iran); SFI (Ireland); INFN (Italy); NRF and WCU (Korea); LAS (Lithuania); CINVESTAV, CONACYT, SEP, and UASLP-FAI (Mexico); MSI (New Zealand); PAEC (Pakistan); MSHE and NSC (Poland); FCT (Portugal); JINR (Armenia, Belarus, Georgia, Ukraine, Uzbekistan); MON, RosAtom, RAS and RFBR (Russia); MSTD (Serbia); SEIDI and CPAN (Spain); Swiss Funding Agencies (Switzerland); NSC (Taipei); ThEP, IPST and NECTEC (Thailand); TUBITAK and TAEK (Turkey); NASU (Ukraine); STFC (United Kingdom); DOE and NSF (USA). Individuals have received support from the Marie-Curie program and the European Research Council (European Union); the Leventis Foundation; the A. P. Sloan Foundation; the Alexander von Humboldt Foundation; the Austrian Science Fund (FWF); the Belgian Federal Science Policy Office; the Fonds pour la Formation à la Recherche dans l’Industrie et dans l’Agriculture (FRIA-Belgium); the Agentschap voor Innovatie door Wetenschap en Technologie (IWT- Belgium); the Ministry of Education, Youth and Sports (MEYS) of Czech Republic; the Council of Science and Industrial Research, India; the Compagnia di San Paolo (Torino); and the HOMING PLUS program of Foundation for Polish Science, cofinanced from European Union, Regional Development Fund. [1] D. E. Groom et al. (Particle Data Group), Eur. Phys. J. C 15, 357 (2000). [2] H. Flaecher, M. Goebel, J. Haller, A. Hoecker, K. Moenig, and J. Stelzer, Eur. Phys. J. C 60, 543 (2009). [3] M. Buchkremer, J.-M. Gerard, and F. Maltoni, J. High Energy Phys. 06 (2012) 135. [4] N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963). [5] M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 652 (1973). [6] B. Pontecorvo, Sov. Phys. JETP 6, 429 (1957). [7] Z. Maki, M. Nakagawa, and S. Sakata, Prog. Theor. Phys. 28, 870 (1962). [8] E. Asilar, E. Cavlan, O. Dogangun, S. Kefeli, E. Ozcan, M. Sahin, and G. Unel, Eur. Phys. J. C 72, 1966 (2012). [9] G. Aad et al. (ATLAS), Phys. Rev. Lett. 109, 032001 (2012). [10] G. Aad et al. (ATLAS), Phys. Rev. D 86, 012007 (2012). [11] G. Aad et al. (ATLAS), J. High Energy Phys. 04 (2012) 069. [12] G. Aad et al. (ATLAS), Phys. Rev. Lett. 108, 261802 (2012). [13] S. Chatrchyan et al. (CMS), Phys. Lett. B 716, 103 (2012). [14] S. Chatrchyan et al. (CMS), J. High Energy Phys. 05 (2012) 123. [15] S. Chatrchyan et al. (CMS), Phys. Lett. B 718, 307 (2012). COMBINED SEARCH FOR THE QUARKS OF A . . . PHYSICAL REVIEW D 86, 112003 (2012) 112003-9 http://dx.doi.org/10.1007/BF02683445 http://dx.doi.org/10.1007/BF02683445 http://dx.doi.org/10.1140/epjc/s10052-009-0966-6 http://dx.doi.org/10.1140/epjc/s10052-009-0966-6 http://dx.doi.org/10.1007/JHEP06(2012)135 http://dx.doi.org/10.1007/JHEP06(2012)135 http://dx.doi.org/10.1103/PhysRevLett.10.531 http://dx.doi.org/10.1143/PTP.49.652 http://dx.doi.org/10.1143/PTP.49.652 http://dx.doi.org/10.1143/PTP.28.870 http://dx.doi.org/10.1143/PTP.28.870 http://dx.doi.org/10.1140/epjc/s10052-012-1966-5 http://dx.doi.org/10.1103/PhysRevLett.109.032001 http://dx.doi.org/10.1103/PhysRevLett.109.032001 http://dx.doi.org/10.1103/PhysRevD.86.012007 http://dx.doi.org/10.1007/JHEP04(2012)069 http://dx.doi.org/10.1007/JHEP04(2012)069 http://dx.doi.org/10.1103/PhysRevLett.108.261802 http://dx.doi.org/10.1103/PhysRevLett.108.261802 http://dx.doi.org/10.1016/j.physletb.2012.07.059 http://dx.doi.org/10.1007/JHEP05(2012)123 http://dx.doi.org/10.1007/JHEP05(2012)123 http://dx.doi.org/10.1016/j.physletb.2012.10.038 [16] M. S. Chanowitz, M.A. Furman, and I. Hinchliffe, Nucl. Phys. B153, 402 (1979). [17] K. Nakamura et al. (Particle Data Group), Phys. Rev. D 86, 010001 (2012). [18] V.M. Abazov et al. (D0), Phys. Rev. D 85, 091104 (2012). [19] Y. Chao, K.-F. Chen, S.-K. Chen, W.-S. Hou, B.-Y. Huang, and Y.-J. Lei, Phys. Rev. D 84, 014029 (2011). [20] S. Chatrchyan et al. (CMS), JINST 3, S08004 (2008). [21] CMS Collaboration (CMS), Commissioning of the Particle-Flow Reconstruction in Minimum-Bias and Jet Events from Collisions at 7 TeV, CMS Physics Analysis Summary CMS-PAS-PFT-10-002 (2010), http://cdsweb .cern.ch/record/1279341. [22] M. Cacciari, G. P. Salam, and G. Soyez, J. High Energy Phys. 04 (2008) 063. [23] S. Alioli, P. Nason, C. Oleari, and E. Re, J. High Energy Phys. 09 (2009) 111. [24] E. Re, Eur. Phys. J. C 71, 1547 (2011). [25] T. Sjöstrand, S. Mrenna, and P. Z. Skands, J. High Energy Phys. 05 (2006) 026. [26] J. Alwall, M. Herquet, F. Maltoni, O. Mattelaer, and T. Stelzer, J. High Energy Phys. 06 (2011) 128. [27] R. Field, arXiv:1010.3558. [28] M. L. Mangano, M. Moretti, F. Piccinini, and M. Treccani, J. High Energy Phys. 01 (2007) 013. [29] J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P. Nadolsky, and W.K. Tung, J. High Energy Phys. 07 (2002) 012. [30] J. Allison et al., IEEE Trans. Nucl. Sci. 53, 270 (2006). [31] S. Chatrchyan et al. (CMS), JINST 6, P11002 (2011). [32] CMS Collaboration (CMS), b-Jet Identification in the CMS Experiment, CMS Physics Analysis Summary CMS-PAS-BTV-11-004 (2011), http://cdsweb.cern.ch/ record/1427247. [33] S. Chatrchyan et al. (CMS), Phys. Rev. D 84, 092004 (2011). [34] S. Chatrchyan et al. (CMS), J. High Energy Phys. 10 (2011) 132. [35] N. Kidonakis, Phys. Rev. D 83, 091503 (2011). [36] N. Kidonakis, Phys. Rev. D 81, 054028 (2010). [37] N. Kidonakis, Phys. Rev. D 82, 054018 (2010). [38] V. Hirschi, R. Frederix, S. Frixione, M.V. Garzelli, F. Maltoni, and R. Pittau, J. High Energy Phys. 05 (2011) 044. [39] J.M. Campbell and R.K. Ellis, Nucl. Phys. B, Proc. Suppl. 205-206, 10 (2010). [40] J.M. Campbell and R.K. Ellis, Phys. Rev. D 60, 113006 (1999). [41] M. Aliev, H. Lacker, U. Langenfeld, S. Moch, P. Uwer, and M. Wiedermann, Comput. Phys. Commun. 182, 1034 (2011). [42] J.M. Campbell, R. Frederix, F. Maltoni, and F. Tramontano, J. High Energy Phys. 10 (2009) 042. [43] A. L. Read, Nucl. Instrum. Methods Phys. Res., Sect. A 425, 357 (1999). [44] J. S. Conway, in Proceedings of PHYSTAT 2011 Workshop on Statistical Issues Related to Discovery Claims in Search Experiments and Unfolding, edited by H. B. Propser and L. Lyons (CERN, 2001), p. 115 [http:// cdsweb.cern.ch/record/1306523]. [45] CMS Collaboration, CMS Physics Analysis Summary Report No. CMS-PAS-SMP-12-008, 2012 [http://cdsweb .cern.ch/record/1434360]. [46] G. Cowan, K. Cranmer, E. Gross, and O. Vitells, Eur. Phys. J. C 71, 1554 (2011). [47] T. Junk, Nucl. Instrum. Methods Phys. Res., Sect. A 434, 435 (1999). [48] A. L. Read, J. Phys. G 28, 2693 (2002). S. Chatrchyan,1 V. Khachatryan,1 A.M. Sirunyan,1 A. Tumasyan,1 W. Adam,2 E. Aguilo,2 T. Bergauer,2 M. Dragicevic,2 J. Erö,2 C. Fabjan,2,b M. Friedl,2 R. Frühwirth,2,b V.M. Ghete,2 J. Hammer,2 N. Hörmann,2 J. Hrubec,2 M. Jeitler,2,b W. Kiesenhofer,2 V. Knünz,2 M. Krammer,2,b I. Krätschmer,2 D. Liko,2 I. Mikulec,2 M. Pernicka,2,a B. Rahbaran,2 C. Rohringer,2 H. Rohringer,2 R. Schöfbeck,2 J. Strauss,2 A. Taurok,2 W. Waltenberger,2 G. Walzel,2 E. Widl,2 C.-E. Wulz,2,b V. Mossolov,3 N. Shumeiko,3 J. Suarez Gonzalez,3 M. Bansal,4 S. Bansal,4 T. Cornelis,4 E. A. De Wolf,4 X. Janssen,4 S. Luyckx,4 L. Mucibello,4 S. Ochesanu,4 B. Roland,4 R. Rougny,4 M. Selvaggi,4 Z. Staykova,4 H. Van Haevermaet,4 P. Van Mechelen,4 N. Van Remortel,4 A. Van Spilbeeck,4 F. Blekman,5 S. Blyweert,5 J. D’Hondt,5 R. Gonzalez Suarez,5 A. Kalogeropoulos,5 M. Maes,5 A. Olbrechts,5 W. Van Doninck,5 P. Van Mulders,5 G. P. Van Onsem,5 I. Villella,5 B. Clerbaux,6 G. De Lentdecker,6 V. Dero,6 A. P. R. Gay,6 T. Hreus,6 A. Léonard,6 P. E. Marage,6 A. Mohammadi,6 T. Reis,6 L. Thomas,6 G. Vander Marcken,6 C. Vander Velde,6 P. Vanlaer,6 J. Wang,6 V. Adler,7 K. Beernaert,7 A. Cimmino,7 S. Costantini,7 G. Garcia,7 M. Grunewald,7 B. Klein,7 J. Lellouch,7 A. Marinov,7 J. Mccartin,7 A.A. Ocampo Rios,7 D. Ryckbosch,7 N. Strobbe,7 F. Thyssen,7 M. Tytgat,7 P. Verwilligen,7 S. Walsh,7 E. Yazgan,7 N. Zaganidis,7 S. Basegmez,8 G. Bruno,8 R. Castello,8 L. Ceard,8 C. Delaere,8 T. du Pree,8 D. Favart,8 L. Forthomme,8 A. Giammanco,8,c J. Hollar,8 V. Lemaitre,8 J. Liao,8 O. Militaru,8 C. Nuttens,8 D. Pagano,8 A. Pin,8 K. Piotrzkowski,8 N. Schul,8 J.M. Vizan Garcia,8 N. Beliy,9 T. Caebergs,9 E. Daubie,9 G. H. Hammad,9 G.A. Alves,10 M. Correa Martins Junior,10 D. De Jesus Damiao,10 T. Martins,10 M. E. Pol,10 M.H.G. Souza,10 W. L. Aldá Júnior,11 W. Carvalho,11 A. Custódio,11 E.M. Da Costa,11 C. De Oliveira Martins,11 S. Fonseca De Souza,11 D. Matos Figueiredo,11 L. Mundim,11 H. Nogima,11 V. Oguri,11 W. L. Prado Da Silva,11 A. Santoro,11 L. Soares Jorge,11 A. Sznajder,11 T. S. Anjos,12,d C. A. Bernardes,12,d F. A. Dias,12,e T.R. Fernandez Perez Tomei,12 E.M. Gregores,12,d C. Lagana,12 S. CHATRCHYAN et al. PHYSICAL REVIEW D 86, 112003 (2012) 112003-10 http://dx.doi.org/10.1016/0550-3213(79)90606-0 http://dx.doi.org/10.1016/0550-3213(79)90606-0 http://dx.doi.org/10.1103/PhysRevD.86.010001 http://dx.doi.org/10.1103/PhysRevD.86.010001 http://dx.doi.org/10.1103/PhysRevD.85.091104 http://dx.doi.org/10.1103/PhysRevD.84.014029 http://dx.doi.org/10.1088/1748-0221/3/08/S08004 http://cdsweb.cern.ch/record/1279341 http://cdsweb.cern.ch/record/1279341 http://dx.doi.org/10.1088/1126-6708/2008/04/063 http://dx.doi.org/10.1088/1126-6708/2008/04/063 http://dx.doi.org/10.1088/1126-6708/2009/09/111 http://dx.doi.org/10.1088/1126-6708/2009/09/111 http://dx.doi.org/10.1140/epjc/s10052-011-1547-z http://dx.doi.org/10.1088/1126-6708/2006/05/026 http://dx.doi.org/10.1088/1126-6708/2006/05/026 http://dx.doi.org/10.1007/JHEP06(2011)128 http://arXiv.org/abs/1010.3558 http://dx.doi.org/10.1088/1126-6708/2007/01/013 http://dx.doi.org/10.1088/1126-6708/2002/07/012 http://dx.doi.org/10.1109/TNS.2006.869826 http://dx.doi.org/10.1109/TNS.2006.869826 http://dx.doi.org/10.1088/1748-0221/6/11/P11002 http://cdsweb.cern.ch/record/1427247 http://cdsweb.cern.ch/record/1427247 http://dx.doi.org/10.1103/PhysRevD.84.092004 http://dx.doi.org/10.1103/PhysRevD.84.092004 http://dx.doi.org/10.1007/JHEP10(2011)132 http://dx.doi.org/10.1007/JHEP10(2011)132 http://dx.doi.org/10.1103/PhysRevD.83.091503 http://dx.doi.org/10.1103/PhysRevD.81.054028 http://dx.doi.org/10.1103/PhysRevD.82.054018 http://dx.doi.org/10.1007/JHEP05(2011)044 http://dx.doi.org/10.1007/JHEP05(2011)044 http://dx.doi.org/10.1016/j.nuclphysbps.2010.08.011 http://dx.doi.org/10.1016/j.nuclphysbps.2010.08.011 http://dx.doi.org/10.1103/PhysRevD.60.113006 http://dx.doi.org/10.1103/PhysRevD.60.113006 http://dx.doi.org/10.1016/j.cpc.2010.12.040 http://dx.doi.org/10.1016/j.cpc.2010.12.040 http://dx.doi.org/10.1088/1126-6708/2009/10/042 http://dx.doi.org/10.1016/S0168-9002(98)01347-3 http://dx.doi.org/10.1016/S0168-9002(98)01347-3 http://cdsweb.cern.ch/record/1306523 http://cdsweb.cern.ch/record/1306523