Search for physics beyond the standard model in events with τ leptons, jets, and large transverse momentum imbalance in pp collisions at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\sqrt{s}=7\ \mathrm{TeV}$\end{document}

A search for physics beyond the standard model is performed with events having one or more hadronically decaying τ leptons, highly energetic jets, and large transverse momentum imbalance. The data sample corresponds to an integrated luminosity of 4.98 fb−1 of proton-proton collisions at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\sqrt{s}=7 ~\text {TeV} $\end{document} collected with the CMS detector at the LHC in 2011. The number of observed events is consistent with predictions for standard model processes. Lower limits on the mass of the gluino in supersymmetric models are determined.


Introduction
The standard model (SM) of particle physics has been successful in explaining a wide variety of data. In spite of this, the SM is incomplete. For example, it possesses a divergence in the Higgs sector [1] and has no cold dark matter (DM) candidate [2]. Many models of physics beyond the SM (BSM) have been proposed in order to address these problems.
DM particles, if produced in proton-proton collisions at the CERN Large Hadron Collider (LHC), would escape detection and result in a significant transverse momentum (p T ) imbalance in the detector. Additionally, cascade decays of heavy colored particles to final states with a high multiplicity of energetic jets and τ leptons appear very naturally in many BSM physics scenarios. Hence, events with multiple τ lepton candidates, large jet multiplicity, and significant transverse momentum imbalance, represent a distinct signature of new physics. In this paper, focus is placed on final states with hadronically decaying τ leptons. In what follows, the visible part of a hadronically decaying τ lepton will be referred to as τ h . * e-mail: cms-publication-committee-chair@cern.ch In certain models of supersymmetry (SUSY), the lightest supersymmetric particle (LSP) is a candidate for DM. It has been appreciated for some time that the DM relic density may be sensitive to coannihilation processes involving the LSP and the next-to-lightest supersymmetric particle (NLSP). Coannihilation is characterized by a mass difference ( M) between the NLSP and the LSP of approximately 5-15 GeV [3][4][5][6]. This small mass difference would be necessary to allow the NLSP to coannihilate with the LSP in the early universe, leading to the dark matter abundance that is currently observed [7]. If the supersymmetric partner of the τ lepton, the stau ( τ ), is the NLSP, and if the τ decays primarily to a τ lepton and the LSP, small values of M would lead to final states with low-energy τ leptons (p T ∼ M) [8]. Decays of colored SUSY particles can produce the τ via chargino ( χ ± ) or neutralino ( χ 0 ) intermediate states (e.g., χ 0 2 → τ τ → τ τ χ 0 1 ), resulting in final states with at least one τ h .
We present a search for BSM particles in events with exactly one τ h lepton and jets (single-τ h final state), and in events with jets and two or more τ h leptons (multiple-τ h final state). These two topologies provide complementary sensitivity to models with a wide range of M values. For example, in the case of very small values of M (∼5 GeV), the low-energy τ h cannot be effectively detected and the search for new physics in the single-τ h final state has better sensitivity. The analysis is performed using proton-proton collision data at √ s = 7 TeV collected with the Compact Muon Solenoid (CMS) detector [9] at the LHC in 2011. The data sample corresponds to an integrated luminosity of 4.98 ± 0.11 fb −1 . The search is characterized by methods that determine the backgrounds directly from data, to reduce the reliance on simulation. To illustrate the sensitivity of this search to BSM processes, the constrained minimal supersymmetric extension of the standard model, or minimal supergravity, is chosen as the benchmark [3,10,11]; we denote this benchmark as "CMSSM". An interpretation of the results in the context of simplified model spectra (SMS) [12,13] is also presented. The ATLAS collaboration has published a result on a search for one or more hadronically decaying tau leptons, highly energetic jets, and a large transverse momentum imbalance probing minimal Gauge Mediated Symmetry Breaking (GMSB) models [14].

The CMS detector
The central feature of the CMS apparatus is a superconducting solenoid, of 6 m inner diameter, providing a magnetic field of 3.8 T. Within the field volume are a silicon pixel and strip tracker, a crystal electromagnetic calorimeter (ECAL), which includes a silicon sensor preshower detector in front of the ECAL endcaps, and a brass-scintillator hadron calorimeter. Muons are measured in gas-ionization detectors embedded in the steel return yoke. In addition to the barrel and endcap detectors, CMS has extensive forward calorimetry.
The inner tracker measures charged particles within |η| < 2.5 and provides an impact parameter resolution of about 15 µm and a p T resolution of about 1.5 % for 100 GeV particles. Collision events are selected with a first-level trigger based on fast electronics, and a higher-level trigger that runs a version of the offline reconstruction program optimized for speed.
The CMS experiment uses a right-handed coordinate system, with the origin at the nominal interaction point, the x axis pointing to the center of the LHC ring, the y axis pointing up (perpendicular to the plane of the LHC ring), and the z axis along the counterclockwise beam direction. The polar angle θ is measured from the positive z axis and the azimuthal angle in the x-y plane. The pseudorapidity is given by η = − ln[tan(θ/2)].

Object reconstruction and identification
Jets in the detector are reconstructed using particle-flow (PF) objects [15]. In the PF approach, information from all subdetectors is combined to reconstruct and identify final-state particles (muons, electrons, photons, and charged and neutral hadrons) produced in the collision. The anti-k T clustering algorithm [16] with a distance parameter R = 0.5 is used for jet clustering. Jets are required to satisfy criteria designed to identify anomalous behavior in the calorimeters, and to be well separated from any identified τ lepton.
Validation and efficiency studies are performed utilizing events with a τ h lepton and a light-lepton , with representing an electron (e) or muon (μ). Muons are reconstructed using the tracker and muon chambers. Selection requirements based on the minimum number of hits in the silicon tracker, pixel detector, and muon chambers are applied to suppress muon backgrounds from decays-in-flight of pions or kaons [17]. Electrons are reconstructed by combining tracks with ECAL clusters. Requirements are imposed to distinguish between prompt and non-prompt electrons, where the latter can arise from charged pion decay or photon conversion [18]. The light-lepton candidates are required to satisfy both track and ECAL isolation requirements. The track isolation variable is defined as the sum of the p T of the tracks, as measured by the tracking system, within an isolation cone of radius R = ( η) 2 + ( φ) 2 = 0.4 centered on the light-lepton track. The ECAL isolation variable is based on the amount of energy deposited in the ECAL within the same isolation cone. In both cases the contribution from the light-lepton candidate is removed from the sum.
Reconstruction of hadronically decaying τ leptons is performed using the hadron-plus-strips (HPS) algorithm [19], designed to optimize the performance of τ h reconstruction by considering specific τ h decay modes. To suppress backgrounds in which light-quark or gluon jets mimic hadronic τ decays, a τ h candidate is required to be spatially isolated from other energy deposits in the calorimeter. Charged hadrons and photons not considered in the reconstruction of the τ h decay mode are used to calculate the isolation. Additionally, τ h candidates are required to be distinguished from electrons and muons in the event. In this analysis, two HPS isolation definitions are used. The τ h isolation definition used for single-τ h final states rejects a τ h candidate if one or more charged hadrons with p T > 1.0 GeV or one or more photons with transverse energy E T > 1.5 GeV is found within an isolation cone of radius R = 0.5. The τ h isolation definition used for multiple-τ h final states rejects a τ h candidate if one or more charged hadrons with p T > 1.5 GeV or one or more photons with transverse energy E T > 2.0 GeV is found within an isolation cone of radius R = 0.3. The isolation criteria used for the multiple-τ h final state increases the signal-to-background ratio while reducing the rate of τ h misidentification. This affects the yield of events with light-quark or gluon jets that are misidentified as τ h leptons, which depends on the square of the misidentification rate. Here a final state with exactly two τ h candidates is considered since events with more than two τ h candidates are only a small fraction (<1 %) of events.
The missing transverse momentum / H T is defined as: where the sum runs over all the jets with p jet T > 30 GeV inside the fiducial detector volume of |η| < 5. The vector / H T is the negative of the vector sum in Eq. (1). The observable H T = p jet T is used to estimate the overall energy scale of the event. For the single-τ h final state, H T is calculated using jets with p T > 50 GeV and will be referred to as H 50 T . For the multiple-τ h final state, H T is calculated using jets with p T > 30 GeV and will be referred to as H 30 T . In both instances of the H T calculation, we consider all jets in |η| < 5 (the fiducial detector limit). The use of a lower p T threshold for the jets in the multiple-τ h final state increases the efficiency of signal events without significantly increasing the background.

Signal and background samples
The major sources of SM background are top-quark pair (tt) events and events with a W or Z boson accompanied by jets. Both tt and W + jets events can have genuine τ h leptons, large genuine / H T from W boson decays, and jets that can be misidentified as a τ h . Similarly, Z + jets events with Z(→ νν) and with one or more jets misidentified as a τ h lepton provide a source of background. Z + jets events with Z(→ νν) present a background because of the genuine τ h leptons and the genuine / H T from the neutrinos in the τ h decay. QCD multijet events can become a background when a mismeasured jet gives rise to large / H T and jets are misidentified as τ h leptons.
Data are compared with predictions obtained from samples of Monte Carlo (MC) simulated events. Signal and background MC samples are produced with the PYTHIA 6.4.22 [20] and MADGRAPH [21] generators using the Z2 tune [22] and the NLO CTEQ6L1 parton distribution function (PDF) set [23]. The τ lepton decays are simulated with the TAUOLA [24] program. The generated events are processed with a detailed simulation of the CMS apparatus using the GEANT4 package [25]. The MC yields are normalized to the integrated luminosity of the data using next-toleading order (NLO) cross sections [26][27][28][29][30][31]. For the 2011 LHC running conditions, the mean number of interactions in a single beam crossing is ∼10. The effect of multiple interactions per bunch crossing (pileup) is taken into account by superimposing MC minimum-bias events so that the probability distribution for overlapping pp collisions in the simulation matches the measured distribution.

Event selection
Events for both the single-and multiple-τ h final states are selected using a trigger that requires / H T > 150 GeV. This trigger allows us to maintain sensitivity in regions where the p T value of the τ h is small (p T ∼ 15 GeV). This trigger efficiency, for an offline selection requirement of / H T > 250 GeV, is 98.9 %. For the τ h efficiency and validation studies, samples are chosen using triggers that require the presence of both a τ h candidate and a muon.

Background estimate
The background contributions are categorized differently for the single-and multiple-τ h final states. For the single-τ h final state, the background contributions are divided into events containing a genuine τ h and events where a jet is misidentified as a τ h . For the multiple-τ h final state, the main background contribution arises from misidentified τ h leptons. We identify the different sources of background individually using dedicated data control regions (CR).

Estimate of backgrounds in the single-τ h final state
In the single-τ h final state, the largest background contribution comes from W + jets events that contain a genuine τ h lepton. The other significant contribution arises from QCD multijet events in which a jet is misidentified as a τ h . The W + jets background contribution is estimated using a sample of W + jets events with W → μν. The QCD multijet background is determined by selecting a QCD-dominated CR and evaluating the τ h misidentification rate.

Estimate of the W + jets background in the single-τ h final state
To evaluate the W + jets background, we exploit the similarity between W decays to a muon and to a tau lepton and select a sample of W + jets events with W(→ μν). This sample will be referred to as the muon control sample. To select the muon control sample, events are required to contain exactly one muon and no reconstructed τ h or electron.
To emulate the τ h acceptance, the muon is required to satisfy |η| < 2.1. The yields in the muon control sample are corrected for muon reconstruction (ε reco μ ) and isolation efficiency (ε iso μ ). The muon reconstruction efficiency is derived from data using a sample of Z + jets events and parameterized as a function of p T and η. The muon isolation criteria help to distinguish between muons from the decay of the W boson and muons from semileptonic decays of c and b quarks. The isolation efficiency is parameterized as a function of the separation from the nearest jet and the momentum of the jet. A correction factor (P W μ ) is applied to the muons in the muon control sample to account for muons that do not come from a τ -lepton decay. This correction factor depends on the p T of the muon and the / H T value in the event and is derived from a simulated sample of W + jets events.
As the muons in the muon control sample are selected to mimic a τ h , a correction is applied to emulate the probability to reconstruct and identify a τ h lepton. The reconstruction and identification efficiency ε reco τ is parameterized as a function of the p T of the τ h candidate and as a function of the total number N of charged particles and photons in the isolation cone [ Fig. 1(a)]. Corrections are also applied to account for the hadronic branching fraction (f bf(hadr) τ ) of a τ lepton. Except for the f bf(hadr) τ the values of the correction factors differ in each event. The corrections are combined to define an overall event weight, defined as: A τ h response template is derived from simulated events. The response template is given by the ratio of the reconstructed energy of the τ h to the true generator-level energy. The τ h response depends on the transverse momentum of the generated τ lepton [ Fig. 1(b)] and on the number of reconstructed primary vertices in the event. The muon p T spectrum is smeared as a function of p T and the number of primary vertices to mimic the p T distribution of the τ h .
Fully simulated W+jets events are used to verify the procedure. Figure 2 shows the H 50 T and / H T distributions from simulated W + jets events for the single-τ h final state. These events satisfy the baseline selection described in Sect. 5. The reconstructed τ h is required to match a hadronically decaying generated tau lepton, to ensure that only the genuine tau background is addressed in this check. The event yield and distributions are compared with the prediction from the simulated muon control sample and agree within statistical uncertainties, thus verifying the closure of the method in MC simulation. Hence, the predicted H 50 T and / H T distributions from the muon control sample can be taken to describe a τ h sample within statistical uncertainties.
The muon control sample consists primarily of W + jets events, but also contains tt events in which one W boson decays into a muon while the other W boson decays either into an unidentified τ lepton or into a light lepton that is not reconstructed. Any isolated muons produced through the decay of b or c quarks can also contribute to the muon control sample. SM processes containing a Z boson or two W bosons can also contribute to the muon control sample if one of the two decay muons is not reconstructed.
The true event yields of each process as determined from simulation are summarized in Table 1 for the baseline and SR selections. For both selections the number of predicted events with a genuine τ h lepton is seen to agree with the true number of events. The value of ε reco τ that is used to calculate the predicted rate is measured in a sample of W + jets and is different from the value that would be measured in a sample of tt events. This leads to an overestimation of the tt contribution. A systematic uncertainty is assigned to account for this overestimation.

Estimate of the QCD multijet background in the single-τ h final state
To estimate the background where a jet is misidentified as a τ h lepton, a QCD-dominated control sample is obtained by selecting events with H 50 T > 350 GeV and 40 < / H T < 60 GeV. The control sample is selected using a prescaled H T trigger with criteria that lead to a sample where about 99 % of the events arise from QCD multijet production. The probability for a jet to be misidentified as a τ h lepton is measured by determining the fraction of jets from the single-τ h control sample that pass the τ h identification criteria. Jets considered in the calculation of the misidentification rate satisfy the requirements p T > 5 GeV and |η| < 2.5. The misidentification rates f i for each jet i depend on η and p T and are used to determine an overall weight, which is applied to each event.
The event weights are defined as: where n is the number of jets. The measured misidentification rates shown in Fig. 3(a) are applied to data events in the region with H 50 T > 350 GeV and with H 50 T > 600 for two regions of / H T : 60 < / H T < 80 GeV and 80 < / H T < 100 GeV. These four regions are dominated by QCD multijet events. The results for data and simulation, as well as the predicted fraction of QCD multijet events, are shown in Table 2. The ratio of selected events over predicted events is statistically compatible with one and stable over the range of / H T . Figure 3(  with similar selection requirements to those used in the SR, but are enriched with events from the background process in question. Correction factors and selection efficiencies are measured in those CRs and used to extrapolate to the SR. We use the observed jet multiplicity in each CR along with the measured rate at which a jet is misidentified as a τ h to calculate the yield in the SR. The following equation is used to estimate each background contribution B: where N SR B is the predicted rate in the SR, N CR B is the observed number of events in the CR, and α xy is the correction factor for acceptance and efficiency for events in the selected/predicted (data) 1.14 ± 0.26 0.97 ± 0.37 -CR with true physics objects "x" and "y". Here the physics object can be a τ h or a quark or gluon jet. Since the dominant SM backgrounds contribute to the SR when jets are misidentified as τ h lepton, the background estimation strategy outlined in Eq. (4) relies on the determination of the event probability P(m) for at least "m" jets to be misidentified as a τ h , where P(m) is the product of three factors: (i) the probability P (N) for an event to contain N jets, (ii) the number of possible ways for exactly n jets to pass the τ h identification criteria given N possible jets C(N, n) = N !/n!(N − n)!, and (iii) the probability f for a single jet to be misidentified as a τ h . The P(m) terms are given by: Equation (5) would be identical to Eq. (3) if used in the case of the single-τ h final state. Equation (4) is used to estimate the tt, W + jets, and Z + jets background contributions to the SR. The P (N) terms are determined from data using the jet multiplicity distribution in each CR, while the f terms are measured for each background process by determining the fraction of jets in each CR that pass the τ h identification criteria. Since the QCD multijet contribution to the SR for the multiple-τ h final state is negligible according to simulation, a data-to-MC scale factor is used to correct the QCD multijet prediction from simulation. In the sections that follow, the selections used to define high purity CRs are outlined and the correction factors α xy used in Eq. (4) are defined. The fraction of events with two τ h leptons is denoted A τ τ , the fraction with one τ h lepton and one jet misidentified as a τ h lepton is denoted A τj , and the fraction with two jets misidentified as τ h leptons is denoted A jj .

Estimate of the tt event background to the multiple-τ h final state
To estimate the contribution of tt events to the multipleτ h SR, a CR is selected by removing the τ h isolation re-quirement and by requiring the presence of at least two b-quark jets (b jets), identified using the track-countinghigh-efficiency (TCHE) algorithm at the medium working point [32]. Because QCD multijet, W + jets, Z(→ τ τ ) + jets and Z(→ νν) + jets events are unlikely to contain two b jets, this requirement provides a sample in which about 99 % of the events are tt events, according to simulation. Figure 4(a) shows the p T distribution of τ h leptons in the tt CR for data and simulation. According to simulation, the fraction of events in the tt control sample that contains one genuine τ h is A τj = 0.166 ± 0.011, while the fraction without a genuine τ h is A jj = 0.834 ± 0.025. The genuine τ h τ h contribution is negligible (A τ τ ∼ 0) according to simulation. Incomplete knowledge of the genuine τ h τ h contribution is included as a source of systematic uncertainty in the tt background prediction. Therefore, α τj in Eq. (4) is given by A τj ε iso τ /P (2 b jets), where ε iso τ is the probability for a τ h lepton to pass the isolation requirement, while α jj is given by A jj /P (2 b jets). The probability P (2 b jets) to identify two or more b jets is determined by the b jet identification efficiency factor [32]. The number of tt events in the SR is calculated as: The probability for a jet in a tt event to be misidentified as a τ h lepton has an average measured value of f = 0.022 ± 0.004. Cross checks are made to validate the use of the bjet identification efficiency as measured in Ref. [32] for this analysis. The estimated tt contribution in the SR is determined to be N SR tt = 2.03 ± 0.36.

Estimate of the Z(→ νν) + jets event background to the multiple-τ h final state
The contribution of Z(→ νν) + jets events to the multipleτ h SR is evaluated by selecting a sample of Z(→ μμ) + jets events and treating the muons as neutrinos. The sample is collected using a trigger designed to select a muon and a τ h . Jet selection criteria similar to those used for the SR sample are imposed. In addition, we require two muons passing the criteria outlined in Sect. 3. The control sample has a purity of about 99 % as estimated from simulation. The / H T distribution for events in this CR is shown in Fig. 4(b). The Z(→ νν) + jets background is estimated by interpreting the p T of the pair of muons as / H T . In order to predict the Z(→ νν) + jets rate in the SR, the Z(→ μμ) + jets sample is corrected for the ratio of the branching fractions R = B(Z → νν)/B(Z → μμ), for trigger efficiencies, for the geometric acceptance A μ as measured from simulation, and for the reconstruction efficiency ε reco μ as measured from data. Therefore, α jj in Eq. (4) is given by: Since there is no prompt production of a genuine τ h in the Z(→ μμ) + jets sample, α τj = 0 and α τ τ = 0. The Z(→ νν) + jets contribution to the SR is calculated as: H T > 250 GeV. The muon identification efficiency ε μ is measured using a "tag-and-probe" method. The probability for a jet to be misidentified as a τ h lepton has a measured value of f = 0.016 ± 0.002. The estimated Z(→ νν) + jets contribution to the SR is determined to be N SR Z(→νν) = 0.03 ± 0.02.

Estimate of the Z(→ τ τ ) + jets event background to the multiple-τ h final state
The contribution from Z → τ τ events is determined with the Z(→ μμ) + jets CR sample used to estimate the background from Z(→ νν) + jets, with the muons treated as τ h leptons. The α xy factors are more difficult to estimate for Z → τ τ events since there are several ways in which Z → τ τ events can contribute to the SR: (i) both τ h leptons pass the kinematic acceptance and identification criteria; (ii) both τ h leptons pass the kinematic acceptance criteria, but only one passes the identification criteria; (iii) one τ h fails the kinematic acceptance criteria, while the other τ h passes both the kinematic acceptance and identification criteria; or (iv) both τ h leptons fail the kinematic acceptance criteria. The Z(→ τ τ ) + jets contribution to the SR is calculated as: where R is given by: A τ is the τ h acceptance, ε τ is the τ h identification efficiency in this control sample, and f = 0.016 ± 0.002. The estimated Z (→ τ τ ) + jets contribution to the SR is determined to be N SR Z(→τ τ ) = 0.21 ± 0.13.

Estimate of the W + jets event background to the multiple-τ h final state
To select the W + jets CR, the τ h isolation requirement, which discriminates between a τ h lepton and other jets, is removed from the SR selection requirements. However, the lack of the τ h isolation requirement increases the contribution from other backgrounds as most of the backgrounds arise because jets are misidentified as a τ h lepton. To minimize the contribution from tt production, events are required to have no jets identified as a b jet. This requirement reduces the contamination from tt events to around 5 %. The purity of the W+jets CR is approximately 65 %. Figure 4(c) shows the / E T distribution, defined as the magnitude of the negative of the vector sum of the transverse momentum of all PF objects in the event, for events in the W + jets CR. The contributions of QCD multijet, tt, and Z(→ νν) + jets events are subtracted in order to determine the number of W + jets events in the CR. The predicted rates for QCD multijet, tt, and Z(→ νν) + jets events are determined by extrapolating from their corresponding CRs. Since there is no genuine multiple-τ h production in W + jets, α τ τ = 0. According to simulation, the fraction of events in the CR with one genuine τ h is A τj = 0.149 ± 0.016, while the fraction of events without a genuine τ h is A jj = 0.851 ± 0.038. Therefore, α τj in Eq. (4) is given by A τj ε iso τ /P (0 b jets), where ε iso τ is the probability for a τ h to pass the isolation requirement and P (0 b jets) is the probability to not have any light-quark or gluon jet misidentified as a b jet. Similarly, α jj is given by A jj /P (0 b jets). The contribution of W + jets events to the SR is then calculated as: (2) . (11) The average rate at which jets are misidentified as a τ h lepton is measured to be 0.019 ± 0.001. The rate f b at which light-quark jets or gluon jets are misidentified as a b jet is used to determine P (0 b jets). The estimated W + jets contribution to the SR is determined to be N SR W+jets = 5.20 ± 0.63.

Estimate of the QCD multijet event background
to the multiple-τ h final state QCD multijet events contribute to the multiple-τ h SR when mismeasurements of jet energies lead to large values of / H T and when jets are misidentified as τ h candidates. By removing the τ h isolation criteria and inverting the | φ(j 2 , / H T )| requirement, a QCD CR sample with about 99 % purity is obtained. Figure 4(d) shows the expected and observed / H T distributions for this sample. A scale factor is obtained from this CR and used to correct the signal prediction for QCD multijet events in simulation. The estimated contribution to the SR from QCD multijet events is determined to be N SR QCD = 0.02 ± 0.02.

Systematic uncertainties
Systematic uncertainties are taken into account for both signal and background events and are described separately. Both the signal and background are affected by the systematic uncertainty in the identification of the τ h candidate. The systematic uncertainty for τ h identification is obtained using a Z → τ τ enhanced region and by correcting this cross section by that measured for Z → ee and Z → μμ events. This uncertainty is validated on a control sample of Z → τ τ events. The level of agreement between data and simulation is found to be at the level of 7 %. Further validation of the performance of τ h identification in a SUSY-like environment is performed by selecting a W(→ τ ν → τ h νν) + jets CR with large hadronic activity (H T ) and large transverse momentum imbalance (/ H T ). The level of agreement between the predicted rate for W(→ τ ν → τ h νν) events and the observed number of events is within 7 % and is determined as a function of H T and / H T .

Systematic uncertainties on background events
The principal sources of systematic uncertainty on the background predictions arise from the correction factors, the finite number of events in the CRs, the measured rates at which jets are misidentified as a τ h lepton, and the level of agreement between the observed and predicted numbers of events in CRs. The contributions to the uncertainties on the correction factors are different for each background category. The dominant effect is due to the uncertainty in the τ h identification efficiency. In the multiple-τ h final state, uncertainties in the jet-energy scale (JES) [33] and the τ h -energy scale (TES) [34] are used to evaluate how changes in H T , / H T , and jet kinematics affect the correction factors. The systematic uncertainty on the correction factors due to the JES and TES is at most ∼3 % for all backgrounds. Smaller contributions to the uncertainties in the correction factors arise from the muon reconstruction and isolation efficiency (<1 %), the uncertainty in the branching fractions ( 1 %), and the uncertainties in trigger efficiency (1 %).
The systematic uncertainties on the measured rates for jet misidentification as a τ h lepton are dominated by the size of the jet sample used to measure these rates and range from 2 % for the single-τ h final state to 5.6-10 % for the multipleτ h final state. The level of agreement between the observed and predicted number of events in MC studies of the CRs is used to assign an additional systematic uncertainty and ranges from 2 % for the single-τ h final state to 3 % for the multiple-τ h final state. Finally, the systematic uncertainty arising from statistical uncertainties on the number of events in the CRs ranges from 2-5 % for the multiple-τ h final state to 3-10 % for the single-τ h final state.

Systematic uncertainties on signal events
The main sources of systematic uncertainties in the SR are due to trigger efficiencies, identification efficiencies, the energy, and momentum scales, the luminosity measurement and PDFs. The uncertainty on the luminosity measurement is 2.2 % [35]. Systematic uncertainties on the / H T triggers (2.5 %) are measured using a sample in which around 99 % of the events are tt events, which have a similar topology to events in the SR samples. The systematic uncertainties on the TES and JES (3.0 %) yield an uncertainty on the signal acceptance of 2.3 %. The uncertainty on the / H T scale depends on the uncertainty of the JES (2-5 % depending on the η and p T values of the jet) and on the unclustered energy scale (10 %). Unclustered energy is defined as the energy found "outside" any reconstructed lepton or jet with p T > 10 GeV. The unclustered energy scale uncertainty has a negligible systematic uncertainty on the signal acceptance. The systematic uncertainty due to imprecise knowledge of the PDFs (11 %) is determined by comparing the CTEQ6.6L [36], MSTW 2008 NLO [37], and NNPDF2.1 [38] PDFs with the default PDF [39]. The systematic uncertainty due to the imprecise modeling of the initial-state and final-state radiation [40] is negligible ( 1 %). The systematic uncertainties associated with event pileup are also  H T for the single-τ h final state. The points with errors represent data that satisfy the baseline selection while the filled green (light) and filled blue (dark) areas shows the predicted backgrounds due to events containing a genuine τ h and a misidentified τ h , respectively. The hatched area shows the total uncertainty on the prediction (Color figure online) negligible. Uncertainties on the theoretical cross sections are evaluated by varying the PDFs and by changing the renormalization and factorization scales by a factor of two [26][27][28][29][30][31].

Results
For the single-τ h final state, the number of background events containing a genuine τ h , as well as the number of background events containing a misidentified τ h , are estimated with data. The results for the baseline and the full selection are listed in Table 3. Figure 5 shows the H 50 T and / H T distributions of data and the different background predictions. The observed number of events in data is in agreement with the SM predictions. The largest sources of background for the multiple-τ h final state are from tt and W + jets events. A counting experiment is performed and the background predictions from data are compared with the observed number of events. Table 4 lists these background predictions and the observed number of events in the SR. Figure 6 shows the H 30 T as well as the M eff distributions in the SR, where M eff is the sum T . The background distributions in Fig. 6 are taken from simulation and normalized over the full spectrum. The estimated number of events due to the SM background processes is in agreement with the number of observed events in the SR.

Limits on new physics
The observed numbers of events in the single-τ h and multiple-τ h final states do not reveal any evidence of physics beyond the standard model. Exclusion limits are set using the CL s [41] criterion in the context of the CMSSM [42]. The CMSSM parameter space with tan β = 40, A 0 = −500 GeV, μ > 0, and M t = 173.2 GeV is chosen as a possible scenario with a light τ and a value of M ≤ 20 GeV. The excluded regions are shown for the single-τ h and multiple-τ h final states in Figs. 7(a) and 7(b), respectively. The limits are set using a simple counting experiment. Systematic uncertainties are treated as nuisance parameters and marginalized, and contamination from signal events in the control samples is taken into account. In the CLs method, both the background-only as well as the signal + background hypothesis are used to derive the confidence levels CL s and the resulting limits and the uncertainty bands on the exclusion contours. In the case of very small values of M(∼5 GeV), the lower-energy τ h cannot be effectively detected and only the energetic τ h from the decay of the neutralino can be observed. The search for new physics with a single τ lepton has a better sensitivity in this case. The single-τ h and multiple-τ h topologies thus have complementary sensitivity and together provide coverage for models with a wide range of M values.
Using the limits set by the single-τ h analysis, a common gaugino mass m 1/2 of <495 GeV is excluded at 95 % Confidence Level (CL) for a common scalar mass m 0 of The results for the multiple-τ h final states are also interpreted in the context of SMS [13]. The τ τ SMS scenario (T3tauh) is studied where gluinos are produced in pairs and subsequently decay to τ lepton pairs and an LSP via a neu- Fig. 8 Diagram for the T3tauh SMS model tralino ( g → qq χ 0 2 ; χ 0 2 → τ τ → τ τ χ 0 1 ). The diagram for the T3tauh model is given in Fig. 8. A gluino mass of <740 GeV is excluded at 95 % CL for LSP masses up to 205 GeV (here, the mass of χ 0 2 is the average of the masses of the gluino and the LSP). Figure 9(a) shows the 95 % CL exclusion region obtained for T3tauh. The limits on the mass of the gluino and LSP are shown with a solid red line.
In the simplified GMSB scenario, the τ is the NLSP and decays to a τ lepton and a gravitino G, with a mass of the order of ∼ keV [43][44][45] ( χ 0 2 → τ τ → τ τ G). The topology for this simplified GMSB scenario is similar to that of T3tauh except for the assumption that both the gluinos decay to τ -lepton pairs with a branching fraction of 100 %. Therefore, the results are also interpreted in the simplified GMSB scenario using the T3tauh scenario. The signal acceptance is corrected to account for the final state containing up to four τ leptons. A gluino with mass <860 GeV is excluded at 95 % CL. Figure 9(b) shows the exclusion limits for the simplified GMSB scenario as a function of the gluino mass.
Since the SMS topologies considered in this paper are characterized by two τ leptons in the final state, we do not present SMS limits for the single-τ h final state.

Summary
A search for physics beyond the standard model with one or more hadronically decaying τ leptons, highly energetic jets, and large transverse momentum imbalance in the final state is presented. The data sample corresponds to an integrated luminosity of 4.98 ± 0.11 fb −1 of pp collisions at √ s = 7 TeV collected with the CMS detector. The final number of events selected in data is consistent with the predictions for standard model processes. We set upper limits on the cross sections for the CMSSM, GMSB, and SMS scenarios. Within the CMSSM framework at tan β = 40, a gaugino mass m 1/2 < 495 GeV is excluded at 95 % CL for scalar masses m 0 < 440 GeV. This result sets a lower limit on the mass of the gluino at 1.15 TeV with 95 % CL in this region. In the multiple-τ h final state, a gluino with a mass less than 740 GeV is excluded for the T3tauh simplified model while a gluino with a mass less than 860 GeV is excluded for the simplified GMSB scenario at 95 % CL. Acknowledgements We congratulate our colleagues in the CERN accelerator departments for the excellent performance of the LHC and thank the technical and administrative staffs at CERN and at other CMS institutes for their contributions to the success of the CMS effort. In addition, we gratefully acknowledge the computing centres and personnel of the Worldwide LHC Computing Grid for delivering so effectively the computing infrastructure essential to our analyses. Finally, we acknowledge the enduring support for the construction and operation of the LHC and the CMS detector provided by the following funding agencies: the Austrian Federal Ministry of Science and Research; the Belgian Fonds de la Recherche Scientifique, and Fonds voor Wetenschappelijk Onderzoek; the Brazilian Funding Agencies