Search for exotic decays of a Higgs boson into undetectable particles and one or more photons

A search is presented for exotic decays of a Higgs boson into undetectable particles and one or two isolated photons in pp collisions at a center-of-mass energy of 8 TeV. The data correspond to an integrated luminosity of up to 19.4 fb − 1 collected with the CMS detector at the LHC. Higgs bosons produced in gluon–gluon fusion and in association with a Z boson are investigated, using models in which the Higgs boson decays into a gravitino and a neutralino or a pair of neutralinos, followed by the decay of the neutralino to a gravitino and a photon. The selected events are consistent with the background-only hypothesis, and limits are placed on the product of cross sections and branching fractions. Assuming a standard model Higgs boson production cross section, a 95% conﬁdence level upper limit is set on the branching fraction of a 125 GeV Higgs boson decaying into undetectable particles and one or two iso- lated photons as a function of the neutralino mass. For this class of models and neutralino masses from 1 to 120 GeV an upper limit in the range of 7 to 13% is obtained. Further results are given as a function of the neutralino lifetime, and also for a range of Higgs boson masses.


Introduction
The detailed studies of the properties of the observed Higgs boson [1][2][3] are key components of the LHC physics program. In the standard model (SM) and for a given mass of the Higgs boson, all properties of the Higgs boson are predicted. Physics beyond the SM (BSM) might lead to deviations from these predictions. Thus far, measurements of the Higgs bosons couplings to fermions and bosons and of the tensor structure of the Higgs boson interaction with electroweak gauge bosons show no significant deviations [4,5] with respect to SM expectations.
Measurements of Higgs boson couplings performed for visible decay modes provide constraints on partial decay widths of the Higgs boson to BSM particles. Assuming that the couplings of the Higgs boson to W and Z bosons are smaller than the SM values, this indirect method provides an upper limit on the branching fraction of the 125 GeV Higgs boson to BSM particles of 57% at a 95% confidence level (CL) [4,6]. An explicit search for BSM Higgs boson decays presents an alternative opportunity for the discovery of BSM physics. The observation of a sizable decay branching fraction of the Higgs boson to undetected (e.g. invisible or largely invisible) E-mail address: cms-publication-committee-chair@cern.ch. final states would be a clear sign of BSM physics and could provide a window on dark matter [7][8][9][10].
Several BSM models predict Higgs boson decays to undetectable particles and photons. In certain low-scale supersymmetry (SUSY) models, the Higgs bosons are allowed to decay into a gravitino ( G) and a neutralino ( χ 0 1 ) or a pair of neutralinos [11,12]. The neutralino then decays into a photon and a gravitino, the lightest supersymmetric particle and dark matter candidate. Fig. 1 shows Feynman diagrams for such decay chains of the Higgs boson (H) produced by gluon-gluon fusion (ggH) or in association with a Z boson decaying to charged leptons (ZH).
As the gravitino in these models has a negligible mass [11,12], the remaining parameter is the neutralino mass. If its mass is in  tion of 3 × 10 −4 . The decay H → Zγ has been studied in Z → e + e − and Z → μ + μ − final states. Upper limits on the product of the cross section and branching fraction of about a factor of ten larger than the SM expectation have been set at the 95% CL [13,14]. With the available data set the search presented is not sensitive to this decay, but it is sensitive to enhancements in the Higgs boson decay rates to undetectable particles and photons arising from BSM physics.
Various background processes lead to the signal signatures and are estimated from simulation or from control samples in data. The dominant background processes are from γ + jets events and diboson events in the ggH and ZH search, respectively. Details of the background estimation techniques are discussed in Section 5. The strength of the ZH channel analysis is an almost background-free selection leading to a larger sensitivity in the model-dependent interpretation. While both the ggH and the ZH channels provide sensitivity to BSM Higgs boson signatures, the ggH channel allows a model-independent interpretation of the results.
This analysis presents a first search for decays of a scalar boson to undetectable particles and one or two isolated photons. The scalar boson is produced in ggH or in ZH. The data used correspond to an integrated luminosity of up to 19.4 ± 0.5 fb −1 at a center-of-mass energy of

The CMS experiment
The CMS detector, definitions of angular and spatial coordinates, and its performance can be found in Ref. [15]. The central feature of the CMS apparatus is a superconducting solenoid, of 6 m internal diameter, providing a magnetic field of 3.8 T. The field volume contains a silicon pixel and strip tracker, a crystal electromagnetic calorimeter (ECAL), and a brass and scintillator hadron calorimeter. Muons are measured in gas-ionization detectors embedded in the steel flux-return yoke of the magnet. The first level of the CMS trigger system, composed of specialized hardware processors, is designed to select the most interesting events within 3 μs, using information from the calorimeters and muon detectors. A high-level trigger processor farm is used to reduce the rate to a few hundred events per second before data storage.
A particle-flow algorithm [16,17] is used to reconstruct all observable particles in the event. The algorithm combines all subdetector information to reconstruct individual particles and identify them as charged hadrons, neutral hadrons, photons, and leptons. The missing transverse energy vector E miss T is defined as the negative vector sum of the transverse momenta of all reconstructed particles (charged or neutral) in the event, with E miss T = | E miss T |. Jets are reconstructed using the anti-k T clustering algorithm [18] with a distance parameter of R = 0.5, as implemented in the FastJet package [19,20]. A multivariate selection is applied to separate jets from the primary interaction and those reconstructed due to energy deposits associated with pileup interactions [21]. The discrimination is based on the differences in the jet shapes, on the relative multiplicity of charged and neutral components, and on the different fraction of transverse momentum which is carried by the hardest components. Photon identification requirements and other procedures used in selecting events can be found in Section 4.

Data and simulation events
In the search for Higgs bosons produced in ggH, the trigger system requires the presence of one high transverse energy (E γ T ) photon candidate and significant E miss T . The presence of a photon candidate with E γ T > 30 GeV is required within the ECAL barrel region (|η γ | < 1.44). At the trigger level E miss T is calculated from calorimeter information, and is not corrected for muons. A selection requirement of E miss T > 25 GeV is applied. The efficiency of the trigger is monitored and measured with two control triggers for the photon and the E miss T trigger requirement. The data recorded with this trigger correspond to an integrated luminosity of 7.4 fb −1 and were part of the CMS "data parking" program implemented for the last part of the data taking at √ s = 8 TeV in 2012. In that program, CMS recorded additional data with relaxed trigger requirements, planning for a delayed offline reconstruction in 2013 after the completion of the LHC Run 1. For the search for Higgs bosons produced in ZH, collision events were collected using single-electron and single-muon triggers which require the presence of an isolated lepton with p T in excess of 27 GeV and 24 GeV, respectively. Also a dilepton trigger was used, requiring two leptons with p T thresholds of 17 GeV and 8 GeV. The luminosity integrated with these triggers at Several Monte Carlo (MC) event generators are used to simulate signal and background processes. The simulated samples are used to optimize the event selection, evaluate selection efficiencies and systematic uncertainties, and compute expected event yields. In all cases the MC samples are reweighted to match the trigger efficiency measured in data.
The Vγ , WZ, ZZ, VVV (where V represents W or Z bosons), Drell-Yan (DY) production of qq → Z/γ * , and qq → W + W − processes are generated with the MadGraph 5.1 event generator [22] at leading-order (LO), the gg → W + W − process is generated with the LO event generator gg2ww 3.1 [23], and the tt and tW processes are generated with powheg 1.0 at next-to-leading-order (NLO). The signal samples are also produced with MadGraph. The cross sections at NLO or higher orders if available are used for a given process to renormalize the MC event generators. All processes are interfaced to the pythia 6.4 generator [24] for parton shower and hadronization. The CTEQ6L set of parton distribution functions (PDF) [25] is used for LO generators, while the CT10 [26] PDF set is used for NLO generators. For all processes, the detector response is simulated with a detailed description of the CMS detector, based on the Geant4 package [27]. Additional pp interactions overlapping the event of interest in data, denoted as pileup events, are accounted for by simulating pp interactions with the pythia generator and adding them to each MC sample. The MC samples are tuned to reproduce the distribution in the number of pileup events in data. The average number of pileup events is about 26 for the collected data used in the ggH channel, and is about 21 for the collected data used in the ZH channel.

Event selection
Two strategies are followed to isolate the Higgs boson events produced by ggH and by ZH from the background processes. The signal cross sections are several orders of magnitude smaller than the major reducible background processes, whose contributions are greatly reduced using the event selections described in the following sub-sections.

Event selection in the ggH channel
In the ggH channel, each selected event is required to have at least one photon candidate with E γ T > 45 GeV and |η γ | < 1.44 using a cut-based selection [28,29]. To reduce the SM backgrounds arising from the leptonic decays of W and Z bosons, a lepton veto is applied. Events are rejected if they have one or more electrons fulfilling a loose identification requirement [30] and p e T > 10 GeV and |η e | < 2.5, excluding the transition region of 1.44 < |η e | ≤ 1.57 since the reconstruction of an electron object in this region is not optimal. Similarly, events containing muon candidates with p μ T > 10 GeV, |η μ | < 2.1, and well separated from the photon candidate requiring R(γ , μ) = ( η) 2 + ( φ) 2 > 0.3 (where the φ is azimuthal angle in radians), are rejected. In addition to the selection requirements described above, the E miss T is required to be greater than 40 GeV. This level of selection is referred to as the preselection. Additional selection criteria are applied to search for new physics in either a quasi model-independent way or optimized for a SUSY benchmark model. In this channel jets can arise from initial-state radiation. For both search strategies jets are required to have p j T > 30 GeV and |η j | < 2.4. These jets must not overlap with the photon candidate below R(γ , jet) < 0.5.
In the model-independent analysis, events with two or more jets are rejected. For events with one jet the azimuthal angle between the photon and the jet ( φ(γ , jet)) is required to be smaller than 2.5. This selection requirement rejects the dominant γ + jet background, where the photon and the jet tend to be back-to-back in the transverse plane.
In the model-dependent analysis developed for SUSY scenarios, no requirement is applied on jet multiplicity. In order to minimize the contribution from processes such as γ + jets and multijet events, two methods are used for identifying events with mismeasured E miss T . The E miss T significance method [31] takes account of reconstructed objects for each event and their known resolutions to compute an event-by-event estimate of the likelihood that the observed E miss T is consistent with zero. In addition, a minimization method [29] constructs a χ 2 function of the form where (p T reco ) i are the scalar transverse momenta of the reconstructed objects, such as jets and photons that pass the abovementioned identification criteria, the (σ p T ) i are the expected resolutions in each object, the σ E miss E miss x, y = E miss,reco In events with no genuine E miss T , the mismeasured quantities are re-distributed back into the particle momenta to minimize the ) is less than 45 GeV and the chi-square probability is larger than 10 −3 .
To further suppress multijet backgrounds, events are not considered if the scalar sum of the transverse momenta of the identified jets in the event (H T ) is greater than 100 GeV. An additional requirement is applied on the angle (α) between the beam direction and the major axis of the supercluster [28] in order to reject non-prompt photons that have showers elongated along the beam line.
Finally, the transverse mass, formed by the photon candidate, E miss T , and their opening angle, is required to be greater than 100 GeV. Photons from the continuum Zγ background have a harder spectrum than the photons resulting from the Higgs decay in the SUSY benchmark models considered. To further reduce the continuum Zγ background and for models with m H = 125 GeV a cut of E γ T < 60 GeV is applied. For higher masses the cut is optimized depending on each mass hypothesis going from 60 GeV up to 200 GeV for m H = 400 GeV.
The list of selection criteria used in the model-independent and the SUSY benchmark model analyses are given in Table 1, together with the cumulative efficiencies relative to the preselection for signal and background processes.

Event selection in the ZH channel
The leptonic decays of the Z boson, consisting of two oppositely charged same-flavor high-p T isolated leptons (e + e − , μ + μ − ), are used to tag the Higgs boson candidate events. Large missing transverse energy from the undetectable particles, at least one isolated high-E T photon, and little or moderate jet activity are required to select the signal events.
The details of the lepton candidate selection and missing transverse energy reconstruction are given in Ref. [32]. In addition, photon requirements based on a multivariate selection discussed in Refs. [28,33] have been used. The kinematic selection requires two leptons with p T > 20 GeV and one photon with E γ T > 20 GeV.
Furthermore, the dilepton mass must be compatible with that of a Z boson within 15 GeV of the pole mass.
To reduce the background from WZ events, events are removed if an additional loosely identified lepton is reconstructed with p T > 10 GeV. To reject most of the top-quark background, an event Table 1 Summary of ggH selection for both the quasi model-independent analysis and the analysis with the SUSY benchmark model with the cumulative efficiencies of the selection requirements relative to the preselection for Zγ → ννγ , γ + jet and for a signal in a SUSY benchmark model with ggH production of a Higgs boson with mass 125 GeV decaying into a neutralino of mass 120 GeV and a photon.

Selection requirements
Model-independent SUSY benchmark model Table 2 Summary of ZH selection.

Variable Selection
Leptons 2 leptons, p T > 20 GeV is rejected if it passes the b-tagging selection (anti b-tagging) or if there is a selected jet with p T larger than 30 GeV (jet veto). The b-tagging selection is based on the presence of a muon in the event from the semileptonic decay of a bottom-quark, and on the impact parameters of the constituent tracks in jets containing decays of bottom-quarks [34]. The set of b-tagging veto criteria retain about 95% of the light-quark jets, while rejecting about 70% of the b-jets.
The signal topology is characterized by a Z( ) system with large transverse momentum balanced in the transverse plane by a E miss T + E smaller than 0.5, and the azimuthal angle between the two leptons φ is required to be smaller than 2.25 radians. Finally, p T is required to be larger than 60 GeV, and E miss T is required to be larger than 60 GeV. A summary of the selection for the analysis is shown in Table 2.
The signal-to-background fraction depends on the |η γ |, the pseudorapidity of the photon, with greater discrimination at lower values. To exploit this effect and improve sensitivity, the selected events are subdivided according to whether the photon is reconstructed in the barrel or endcap regions, as explained in Section 7.2.

Background estimation
The background estimation techniques and the composition of all backgrounds in the search with the ggH and ZH signatures are discussed below. The yield for the irreducible background from H → Zγ → ννγ is negligible and is therefore ignored in the analysis.

Background estimation in the ggH channel
The dominant background for the γ + E miss T signal in the ggH channel is the process γ + jet. Other SM backgrounds include Zγ → ννγ , Wγ , W → eν, W → μν, W → τ ν, multijet, and diphoton events. Background events that do not arise from pp collisions are also considered in the analysis. These backgrounds can be categorized broadly into three categories, as described below.

Background estimates from simulation
The γ + jet process surviving the various E miss T selection requirements is one of the most significant backgrounds in this analysis due to the presence of an isolated photon and its large production cross section. The MC normalization of this background is corrected using control samples in data for two event classes, events without jets and those with one or more jets. The control samples in data are obtained using events collected with a prescaled single-photon trigger and with the E miss T requirement reversed to ensure orthogonality to the signal phase space, where the contamination from other processes is minimal. Multiplicative correction factors (C ) are obtained after normalizing the event yield in the simulation to match the data in the control region, separately for events with no jets (C = 1.7) and one or more jets (C = 1.1). These correction factors are used to normalize the simulated γ + jet event yield in the signal region. An uncertainty of 16% is obtained for these correction factors based on the difference between the corrected and uncorrected simulation and the relative fraction of no jet events (about 10% of the events in the control region) and one or more jet events. The background processes Zγ → γ and W → μν contribute only a small fraction of the total background prediction, due to the lepton veto applied at the preselection stage, and are modeled using simulated samples.
To take into account differences between data and simulation due to imperfect MC modeling, various scale factors (SF) are applied to correct the estimates from simulation. These SFs are defined by the ratio of the efficiency in data to the efficiency in simulation for a given selection. The SF for photon reconstruction and identification is estimated from Z → e + e − decays [35] and is consistent with unity.

Background estimates from data
The contamination from jets misidentified as photons (jet → γ ) is estimated in a data control sample enriched with multijet events defined by E miss T < 40 GeV. This sample is used to measure the ratio of the number of candidates that pass the photon identification criteria to those failing the isolation requirements. The numerator of this ratio is further corrected for the photon contamination due to direct photon production using an isolation side band. The Table 3 Summary of all relative systematic uncertainties in percent for the signal and background estimates for the Higgs model (model-independent in parenthesis) selection in the ggH analysis.
Photon energy scale ±1% corrected ratio is applied to data events which pass the denominator selection and all other event requirements in the signal region.
The systematic uncertainty of this method is dominated by the choice of the isolation sideband, and is estimated to be 35% by changing the isolation criteria in the sideband region definition. The other sources of systematic uncertainty are determined by changing the E miss T selection for the control region, and the loose identification requirements on the photons, all of which are found to be comparatively small.
Events with single electrons misidentified as photons (electron → γ ) are another major source of background. This background is estimated with a tag-and-probe method using Z → e + e − events [36]. The efficiency to identify electrons ( γ e ) is estimated in the Z boson peak mass window of 60-120 GeV. The inefficiency (1 − γ e ) is found to be 2.31 ± 0.03%. The ratio (1 − γ e )/ γ e , which represents the electron misidentification rate, is applied to a sample where candidates are required to have hits in the pixel detector, and is used to estimate the contamination in the signal region. The misidentification rate is found to be dependent on the number of vertices reconstructed in the event and the number of tracks associated with the selected primary vertex. The difference in the final yields taking this dependence into account or neglecting it, using the inclusive measurement of γ e , is less than 5%.

Non-collision background estimates from data
The search is susceptible to contamination from non-collision backgrounds, which arise from cosmic ray interactions, spurious signals in the ECAL, and accelerator-induced secondary particles. The distribution of arrival-times of photons from these backgrounds is different from that of prompt photons produced in hard scattering. To quantify the contamination from these backgrounds a fit is performed to the candidate-time distribution using background distributions from the data [29]. The contamination due to out-of-time background contributions is found to be less than one percent of the total background and is therefore not included in the final event yield.

Background estimation in ZH channel
Processes that contribute significantly to the SM expectation in the ZH channel are listed below.

Non-resonant dilepton backgrounds
The contributions from W + W − , top-quark, W + jets, and Z/γ * → τ + τ − processes are estimated by exploiting the lepton flavor symmetry in the final states of these processes [37]. The branching fraction to the e ± μ ∓ final state is twice that of the e + e − or μ + μ − final states. Therefore, the e ± μ ∓ control region is used to extrapolate these backgrounds to the e + e − and μ + μ − channels. The method considers differences between the electron and muon identification efficiencies. The data driven estimates agree well with the number of background events expected when applying the same method to simulation. The small difference between the prediction and the obtained value using simulated events is taken as a systematic uncertainty.
The limited number of simulated events is also considered as part of the systematic uncertainty. In summary, the total systematic uncertainty is about 75%. Only two events were selected in the data control region.

Resonant background with three leptons in the final state
The WZ → ν process dominates the resonant backgrounds with three leptons in the final state. The electron → γ misidentification rate is measured in Z → e + e − events by comparing the ratios of electron-electron versus electron-photon pairs in data and in simulation, as described in Section 5.1.2. The average misidentification rate is 1-2% with the larger values at higher |η| γ .

Resonant background with two leptons in the final state
The WZ → ν process with failure to identify the lepton from W boson decays and the ZZ → 2 2ν process dominate these types of events. The jet → γ misidentification rate is measured in a sample containing a muon and a photon. This sample is expected to be dominated by jets misidentified as photons, with some contamination from W/Zγ events, which are subtracted in the study using the simulated prediction. The misidentification rate is similar to the obtained values in the ggH channel.

Resonant background with no genuine missing transverse energy
The background from Zγ or Z + jets events is predicted by the simulation to be about 15% of the total background. Several data regions are studied to verify that the background is estimated correctly. A good agreement between data and simulation is found in all cases. A 50% uncertainty, the statistical uncertainty of the control region, is taken for these backgrounds estimated from simulation.

Summary of systematic uncertainties
Systematic uncertainties in the background estimates from control samples in data are described in Section 5. A summary of the systematic uncertainties considered in each channel are listed in Tables 3 and 4. A common source of systematic uncertainty is associated with the measurement of the integrated luminosity, determined to 2.6% [38]. The uncertainties in the normalization of signal and  simulation-based backgrounds are obtained by varying the renormalization and factorization scales, and the parton distribution functions [26,[39][40][41][42][43].
Because the model-independent and model-specific selections differ significantly in the ggH channel, the systematic uncertainties are evaluated separately for each selection. The photon energy scale uncertainty [28] of about 1% affects the signal and background predictions by 4% for the model-specific selection and by 0.5% for the model-independent selection. Similarly, the jet energy scale uncertainty affects the signal and background predictions by 2-5% depending on the process and selection. After changing the photon or jet energy scales, the E miss T is also recomputed. In addition, the systematic uncertainty associated with the jet energy resolution and unclustered energy scale are propagated to the E miss T computation, and affect the signal and background predictions by 2-4%. As described in the previous section, a 16% uncertainty is applied to the γ + jet normalization due to the difference in the jet multiplicity distribution between the data and background predictions in the γ + jet control region. The uncertainty due to the pileup modeling is found to be 1%, and is estimated by shifting the central value of the total inelastic cross section within its uncertainty.
In the ZH channel, lepton-reconstruction and identification scale factors are measured using a control sample of Z/γ * → + − events in the Z peak region [36]. The associated uncertainty is about 2% per lepton. The photon identification uncertainty is taken to be 3% [33]. The effect of uncertainties in jet energy scale and E miss T on the analysis is also considered. The uncertainty in the b-tagging efficiency is estimated to be about 0.7% comparing inclusive Z/γ * → + − samples in data and simulation. The total uncertainty in the background estimates in the signal region is 36%, which is dominated by the statistical uncertainties in the data control samples from which they are derived.
The impact of the systematic uncertainties in the ggH channel is relatively important: the sensitivity is increased by about 50% if all the systematic uncertainties are removed, where the normalization uncertainties on the γ + jet and jet → γ background processes dominate. The ZH channel is limited by the statistical uncertainty, and the effect of the systematic uncertainties reduce the sensitivity by less than 10%.
Correlations between systematic uncertainties in the two channels are taken into account. In particular, the main sources of correlated systematic uncertainties are those in the experimental measurements such as the integrated luminosity, photon identification, the jet energy scale, and missing transverse energy resolution. All other systematic uncertainties are uncorrelated between them given the different signal and background processes.

Results
The results from the two searches and their combination are reported in this section. In the absence of deviations from the standard model predictions, the modified frequentist method, CL s [44][45][46], is used to define the exclusion limits.

Model-independent results in the ggH channel
Because of the variety of possible BSM signals that could contribute to this final state, the results are presented for a signal with the model-independent selection described in Section 4. The total number of observed events and the estimated SM backgrounds are summarized in Table 5, and found to be compatible within their uncertainties. Fig. 2

Model-specific results in the ggH channel
Imposing the model-specific selection described in Section 4 for the ggH channel, 1296 events are selected in data with a total estimated background of 1232 ± 188. The yields for this selection are shown in Table 6 and estimated for Higgs boson decays (H → G χ 0 1 , χ 0 1 → Gγ ) assuming the ggH production rate for SM Higgs bosons and a 100% branching fraction for this decay. Fig. 4 shows the transverse energy distribution of photons for data, the background estimates, and signal after the model-dependent selection, except the upper selection on the photon, for the ggH channel.

Results in the ZH channel and combinations
A total of four events are selected with the search in ZH. The background yield is estimated to 4.1 ± 1.8. The numbers of observed and expected events are shown in Table 7. The signal model assumes a SM ZH production rate and a 100% branching fraction to undetectable particles and one or two photons. The      selected ranges. A 95% CL upper limit on the branching fraction of 10% is set for a neutralino mass of 95 GeV. Expected and observed limits are also shown for the decay of possible heavier scalar Higgs bosons as a function of the Higgs boson mass in Fig. 7. The requirement on E γ T used in the ggH channel is removed. A lower threshold on E γ T is added, optimized to max- imize the sensitivity for each mass hypothesis. A combination of the two channels is not performed because the assumption of a common SM Higgs boson cross section is not justified.
As discussed in the introduction, some BSM models predict χ 0 1 neutralinos with sizable lifetimes. The performance of the searches has been evaluated for finite lifetimes without modifying the analysis strategy. The expected and observed limits are shown in    can be caused by a kink in the trajectory or by a lower velocity of the neutralino.

Summary
A search is presented for exotic decays of a Higgs boson into undetectable particles and one or two isolated photons in pp collisions at a center-of-mass energy of 8 TeV. The data correspond to an integrated luminosity of up to 19.4 fb −1 collected with the CMS detector at the LHC. Higgs bosons produced in gluon-gluon fusion or in association with a Z boson are investigated. Models including Higgs boson decays into a gravitino and a neutralino or a pair of neutralinos, followed by the neutralino decay to a gravitino and a photon, are tested. The measurements for the selected events in data are consistent with the background-only hypothesis, and the results are interpreted as limits on the product of cross sections and branching fractions. Assuming a standard model Higgs production cross section, a 95% CL upper limit is set on the branching fraction of a 125 GeV Higgs boson decaying into undetectable particles and one or two isolated photons as a function of the neutralino mass. For neutralino masses from 1 to 120 GeV an upper limit in the range of 7 to 13% is obtained. Further results are given as a function of the neutralino lifetime, and also for a range of Higgs boson masses.

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 centers 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: BMWFW and FWF (Austria); F.R.S.