HiggsSignals: Confronting arbitrary Higgs sectors with measurements at the Tevatron and the LHC

HiggsSignals is a Fortran90 computer code that allows to test the compatibility of Higgs sector predictions against Higgs rates and masses measured at the LHC or the Tevatron. Arbitrary models with any number of Higgs bosons can be investigated using a model-independent input scheme based on HiggsBounds. The test is based on the calculation of a chi-squared measure from the predictions and the measured Higgs rates and masses, with the ability of fully taking into account systematics and correlations for the signal rate predictions, luminosity and Higgs mass predictions. It features two complementary methods for the test. First, the peak-centered method, in which each observable is defined by a Higgs signal rate measured at a specific hypothetical Higgs mass, corresponding to a tentative Higgs signal. Second, the mass-centered method, where the test is evaluated by comparing the signal rate measurement to the theory prediction at the Higgs mass predicted by the model. The program allows for the simultaneous use of both methods, which is useful in testing models with multiple Higgs bosons. The code automatically combines the signal rates of multiple Higgs bosons if their signals cannot be resolved by the experimental analysis. We compare results obtained with HiggsSignals to official ATLAS and CMS results for various examples of Higgs property determinations and find very good agreement. A few examples of HiggsSignals applications are provided, going beyond the scenarios investigated by the LHC collaborations. For models with more than one Higgs boson we recommend to use HiggsSignals and HiggsBounds in parallel to exploit the full constraining power of Higgs search exclusion limits and the measurements of the signal seen at around 125.5 GeV.


Introduction
Searches for a Higgs boson [1][2][3][4][5][6] have been one of the driving factors behind experimental particle physics over many years. Until recently, results from these searches have always been in the form of exclusion limits, where different Higgs mass hypotheses are rejected at a certain confidence level (usually 95 %) by the non-observation of any signal. This has been the case for Standard Model (SM) Higgs searches at LEP [7], the Tevatron [8], and (until July 2012) also for the LHC experiments [9,10]. Limits have also been presented on extended Higgs sectors in theories beyond the SM, where one prominent example are the combined limits on the Higgs sector of the minimal supersymmetric standard model (MSSM) from the LEP experiments [11,12]. To test the predictions of models with arbitrary Higgs sectors consistently against all the available experimental data on Higgs exclusion limits, we have presented the public tool HiggsBounds [13,14], which recently appeared in version 4.0.0 [15,16].
With the recent discovery of a new state-compatible with a SM Higgs boson-by the LHC experiments ATLAS [17] and CMS [18], models with extended Higgs sectors are facing new constraints. It is no longer sufficient to test for nonexclusion, but the model predictions must be tested against the measured mass and rates of the observed state, which contains more information. Testing the model predictions of a Higgs sector with an arbitrary number of Higgs bosons against this Higgs signal 1 (and potentially against other signals of additional Higgs states discovered in the future) is the purpose of a new public computer program, Higgs-Signals, which we present here.
HiggsSignals is a Fortran90/2003 code, which evaluates a χ 2 measure to provide a quantitative answer to the statistical question of how compatible the Higgs search data (measured signal strengths and masses) is with the model predictions. This χ 2 value can be evaluated with two distinct methods, namely the peak-centered and the mass-centered χ 2 method. In the peak-centered χ 2 method, the (neutral) Higgs signal rates and masses predicted by the model are tested against the various signal rate measurements published by the experimental collaborations for a fixed Higgs mass hypothesis. This hypothetical Higgs mass is typically motivated by the signal "peak" observed in the channels with high mass resolution, i.e. the searches for H → γ γ and H → Z Z ( * ) → 4 . In this way, the model is tested at the mass position of the observed peak. In the mass-centered χ 2 method on the other hand, HiggsSignals tries to find for every neutral Higgs boson in the model the corresponding signal rate measurements, which are performed under the assumption of a Higgs boson mass equal to the predicted Higgs mass. Thus, the χ 2 is evaluated at the model-predicted mass position. For this method to be applicable, the experimental measurements therefore have to be given for a certain mass range.
The input from the user is given in the form of Higgs masses, production cross sections, and decay rates in a format similar to that used in HiggsBounds. The experimental data from Tevatron and LHC Higgs searches is provided with the program, so there is no need for the user to include these values manually. However, it is possible for the user to modify or add to the data at will. Like HiggsBounds, the aim is to always keep HiggsSignals updated with the latest experimental results.
The usefulness of a generic code such as Higgs-Signals has become apparent in the last year, given the intense work by theorists to use the new Higgs measurements as constraints on the SM and theories for new physics . With HiggsSignals, there now exists a public tool that can be used for both model-independent and modeldependent studies of Higgs masses, couplings, rates, etc. in a consistent framework. The χ 2 output of HiggsSignals also makes it convenient to use it as direct input to global fits, where a first example application can be found in Ref. [69].
This document serves both as an introduction to the physics and statistical methods used by HiggsSignals and as a technical manual for users of the code. It is organized as follows. Section 2 contains a very brief review of Higgs searches at hadron colliders, focusing on the published data which provides the key experimental input for HiggsSignals and the corresponding theory predictions. In Section 3 we present the HiggsSignals algorithms, including the precise definitions of the two χ 2 methods mentioned above. Section 4 provides the technical description (user manual) for how to use the code. We discuss the performance of HiggsSignals and validate with official fit results for Higgs coupling scaling factors from ATLAS and CMS in Sect. 5. Furthermore, we give some examples of fit results, which can be obtained by interpreting all presently available Higgs measurements. We conclude in Sect. 6. In the appendix, details are given on the implementation of theory mass uncertainties in the mass-centered χ 2 method.

Higgs signals in collider searches
The experimental data used in HiggsSignals is collected at hadron colliders, mainly the LHC, but there are also some complementary measurements from the Tevatron collider. This will remain the case for the foreseeable future, but the HiggsSignals methods can be easily extended to include data from, for instance, a future e + e − linear collider. In this section we give a very brief review of Higgs searches at hadron colliders, focussing the description on the experimental data that provides the basic input for HiggsSignals. For a more complete review see, e.g., Ref. [70][71][72].
Most searches for Higgs bosons at the LHC are performed under the assumption of the SM. This fixes completely the couplings of the Higgs state to fermions and vector bosons, and both the cross sections and branching ratios are fully specified as a function of the Higgs boson mass, m H . Most up-to-date predictions, including an extensive list of references, can be found in [73,74]. This allows experiments to measure one-parameter scalings of the total SM rate of a certain (ensemble of) signal channel(s), so-called signal strength modifiers, corresponding to the best fit to the data. These measurements are the basic experimental input used by Higgs-Signals. Two examples of this (from ATLAS and CMS) are shown in Fig. 1. The left plot (taken from [75]) shows the measured value of the signal strength modifier, which we denote byμ, in the inclusive pp → H → Z Z ( * ) → 4 process as a function of m H (black line). The cyan band gives a ±1 σ uncertainty on the measured rate. Since the signal strength modifier is measured relative to its SM value (μ = 1, displayed in Fig. 1 by a dashed line), this contains also the theory uncertainties on the SM Higgs cross section and branching ratios [73,74,76]. As can be seen from Fig. 1, the measured value ofμ is allowed to take on negative values. In the absence of sizable signal-background interferenceas is the case for the SM-the signal model would not givê μ < 0. This must therefore be understood as statistical downward fluctuations of the data w.r.t. the background expectation (the average background-only expectation isμ = 0). To keepμ as an unbiased estimator of the true signal strength, it is however essential that the full range of values is retained. As we shall see in more detail below, the applicability of HiggsSignals is limited to the mass range for which measurements ofμ are reported. It is therefore highly desirable that experiments publish this information even for mass regions where a SM Higgs signal has been excluded.
A second example of HiggsSignals input, this time from CMS, is shown in the right plot of Fig. 1 (from [77]). This figure summarizes the measured signal strength modifiers for all relevant Higgs decay channels at an interesting value of the Higgs mass, here m H = 125.7 GeV. This particular value is typically selected to correspond to the maximal significance for a signal seen in the data. It is important to note that, once a value for m H has been selected, this plot shows a compilation of information for the separate channels that is also available directly from the mass-dependent plots (as shown in Fig. 1a). Again, the error bars on the measuredμ values correspond to 1σ uncertainties that include both experimental (systematic and statistical) uncertainties, as well as SM theory uncertainties.
The idea of HiggsSignals is to compare the experimental measurements of signal strength modifiers to the Higgs sector predictions in arbitrary models. The model predictions must be provided by the user for each parameter point to be tested. To be able to do this consistently, we here describe the basic definitions that we apply. The production of Higgs bosons at hadron colliders can essentially proceed through five partonic subprocesses: gluon fusion (ggf), vector boson fusion (vbf), associated production with a gauge boson (H W /H Z), or associated production with top quarks (tt H ), see [73,74] for details. In models with an enhanced Higgs coupling to bottom quarks, the process bb → H is usually added. In this five-flavor scheme a b quark parton distribution describes the collinear gluon splitting to pairs of bottom quarks inside the proton. This contribution should be matched consistently, and in most cases, added to the gluon fusion subprocess (as prescribed by the Santander matching procedure [78]). We therefore sometimes refer to the sum of the gluon fusion and bb → H subprocesses as single Higgs production (singleH). Internally, HiggsSignals uses the same LHC cross sections for SM Higgs production at √ s = 7 and 8 TeV as HiggsBounds-4 [15]. The same holds for the reference SM branching ratios, which follow the prescription of the LHC Higgs Cross Section Working Group [73,74], see also [76] for more details. These branching ratios are the same as those used by the LHC experiments.
The theory prediction for the signal strength modifier of one specific analysis, from a single Higgs boson H , is computed in HiggsSignals as where the sum runs over all channels considered in this analysis. A channel is characterized by one specific production and one specific decay mode. The individual channel signal strength is given by and the SM channel weight is The SM weights contain the relative experimental efficiencies, i , for the different channels. Unfortunately, these are rarely quoted in experimental publications. If they are available, these numbers can be used by HiggsSignals, which leads to a more reliable comparison between theory predictions and the experimental data for these channels. In the case of unknown efficiencies, all channels considered by the analysis are treated equally, i.e. we set all i ≡ 1. Note, however, that for many observables approximate numbers for the channel efficiencies can be inferred by reproducing official fit results on scale factors for production cross sections or coupling strengths, which will be further discussed in Sect. 5.2.
One final word of caution should be added here: If the model features a non-standard tensor structure for the par-ticles, which should be confronted with the data, these interactions might lead to observable differences in the experimentally measured kinematic distributions and therefore to changes of the signal acceptance/efficiency of the Higgs analyses. In order to obtain reliable results from HiggsSignals for these types of models, one needs to check whether these effects are negligible. An interface for HiggsSignals, where the user can insert model signal efficiencies for each analysis, which are changed with respect to the SM signal efficiencies, is a planned feature for future development. However, it is impossible to completely unfold this model dependence using only the currently available public information.

Statistical approach in HiggsSignals
As mentioned already in the introduction, HiggsSignals contains two different statistical methods to test models against the experimental data. These methods are complementary, and to provide a full model test it is advisable in many situations to use both simultaneously. Nevertheless, we leave the final choice of method to the user, and we therefore first describe both methods separately, before discussing their combination in Sect. 3.3.
As already touched upon in the previous section, the search results of ATLAS and CMS are reported in the form of the signal strength modifierμ, the ratio of the best-fit signal strength to the expected SM strength of a signal in a certain channel, and its uncertainty μ. In the profile likelihood approach [79] used by the experimental collaborations, μ is derived from the allowed variation of the signal strength multiplier μ around the best fit valueμ. This is calculated using the likelihood ratio λ(μ) = L(μ,θ)/L(μ,θ); the ratio of the likelihood function L for a given μ with nuisance parametersθ optimized at the given value of μ, divided by L forμ andθ optimized simultaneously (see [79] for more details).
The uncertainty ofμ is then calculated using a test statistics based on −2 ln λ(μ). According to [80,81], this can be expressed as where N is the data sample size. Generally, as shown in [79], this converges quite quickly to a central or non-central χ 2 distribution, depending on the nuisance parameters. If the test statistics follows a χ 2 distribution, the uncertainties of the measurement can generally be treated as Gaussian, hence we interpret all uncertainties μ as Gaussian, and neglect the O(1/ √ N ) term. Looking at the experimental results used in HiggsSignals and the available event sample sizes, this is justified in almost all analyses, apart from H → Z Z * , where visible differences from the Gaussian approximation are still possible due to the small event sample size. The largest remaining effects of non-Gaussian distributions are taken into account in HiggsSignals by using asymmetric uncertainties on the measured signal strength in the χ 2 calculation, if published as such by the collaborations.
While the χ 2 calculated in HiggsSignals can be expected to statistically approximate the true −2 ln λ distribution, cf. Eq. (4), there are three relevant experimental input quantities which can systematically affect the accuracy of the HiggsSignals output in case they are not presented in a complete form in the publicly disclosed information: Firstly, the relative efficiencies i of the various Higgs channels/processes considered in the (categories of a) Higgs analysis, as introduced in Eq. (3). Secondly, the correlations of the relevant experimental systematic uncertainties (e.g. of the jet energy scale (JES), e ± /γ identification and energy scale, tagging efficiencies, etc.) between different Higgs search analyses. Thirdly, the use of continuous variables for classification of channels/production processes (e.g. by using multivariate techniques), which cannot be mapped directly onto signal strengths measurements for distinct categories used as experimental input for the χ 2 fit in HiggsSignals. An example for this is the CMS H → Z Z * → 4 analysis [82]. The effects of such an approach and an approximate solution to this problem within HiggsSignals is discussed in Sect. 5.2.
While the signal efficiencies, i , could be provided straight-forwardly for every analysis as public information, the communication of the (correlated) systematics, both from experimental and theoretical sources, used in a given analysis is not common. However, within the Gaussian approximation these could in principle be taken into account in HiggsSignals. For the future it would be desirable if this information was provided in a model-independent way. Some ideas on how information on correlated systematic uncertainties in Higgs boson rate measurements could be communicated can be found in Ref. [83]. We discuss the possible impact of including this information in Sect. 5.2 for a few relevant cases.
The χ 2 based approach in HiggsSignals could in principle be replaced by the use of likelihood curves from the collaborations, which are currently available in (m H ,μ) grids for a few analyses [77,84], albeit not for the categories individually. Once they are available for the majority of analyses and for every single (category of an) analysis, the χ 2 could partly be replaced by the use of these likelihoods. However, significant modifications of the final likelihood by a tool like HiggsSignals would still be required to make it applicable to arbitrary Higgs sectors, due to potentially different signal compositions and hence changed theoretical rate uncertainties. Moreover, the necessity of incorporating cor-related systematics, as mentioned above, remains also in this approach. Already with the currently available statistics the ignorance of efficiencies and correlations of experimental systematics are often the dominant effects for the typically small deviations between the official results by the collaborations and the HiggsSignals results. The assumption on the parabolic shape of the likelihood, on the other hand, has typically a relatively small impact. More details will be given in Sect. 5.2.

The peak-centered χ 2 method
The objective of this method is to perform a χ 2 test for the hypothesis that a local excess, "signal" (or "peak observable"), in the observed data at a specified mass is generated by the model. In short, this test tries to minimize the total χ 2 by assigning, to each Higgs signal in the experimental dataset used, any number of Higgs bosons of the model. From each signal, both the predicted signal strength modifiers and the corresponding predicted Higgs masses (for channels with good mass resolution) enter the total χ 2 evaluation in a correlated way. Schematically, the total χ 2 is given by where N H is the number of (neutral) Higgs bosons of the model. The calculation of the individual contributions from the signal strength modifiers, χ 2 μ , and the Higgs masses, χ 2 m i , will be discussed below.
The input data used in this method is based on the prejudice that a Higgs signal has been observed at a particular Higgs mass value, which does not necessarily have to be the exact same value for all observables. Technically, each observable is defined by a single text file, which contains all relevant information needed by HiggsSignals. An experimental dataset 2 is then a collection of observables, whose text files are stored in a certain subdirectory of the HiggsSignals distribution. Users may add, modify or remove the experimental data for their own purposes, see Sect. 4.6 for more details.
Currently, an obvious and prominent application of the peak-centered χ 2 method would be the test of a single Higgs boson against the rate and mass measurements performed at around 125-126 GeV in all channels reported by the experimental collaborations at the LHC and Tevatron. This scenario will be discussed in detail in Sect. 5. However, Higgs-Signals is implemented in a way that is much more general: Firstly, contributions from other Higgs bosons in the model to the Higgs signals will be considered, and if relevant, included in the test automatically. Secondly, the extension of this test to more Higgs signals (in other mass regions) can simply be achieved by the inclusion of the proper experimental data, or for a phenomenological study, the desired pseudo-data.

Signal strength modifiers
For N defined signal observables, the total χ 2 contribution is given by where the observed and predicted signal strength modifiers are contained in the N -dimensional vectorsμ and µ, respectively. C μ is the signal strength covariance matrix. The signal strength covariance matrix C μ is constructed in the following way. The diagonal elements (C μ ) αα (corresponding to signal observable α) should first of all contain the intrinsic experimental (statistical and systematic) 1 σ uncertainties on the signal strengths squared, denoted by ( μ * α ) 2 . These will be treated as uncorrelated uncertainties, since there is no information publicly available on their correlations. We define these uncorrelated uncertainties by subtracting from the total uncertainty μ α (which is given directly from the 1 σ error band in the experimental data, cf. Fig. 1) the luminosity uncertainty as well as the theory uncertainties on the predicted signal rate (which we shall include later as correlated uncertainties). Hereby, we assume that these uncertainties can be treated as Gaussian errors. This gives Here, L is the relative uncertainty on the luminosity, and c SM a is the SM channel rate uncertainty (for a total of k channels contributing to the analysis with signal α) given by where σ SM a and BR SM a are the relative systematic uncertainties of the production cross section σ a and branching ratio BR a , respectively, of the channel a in the SM. Their values are taken from the LHC Higgs Cross Section Working Group [73,74], evaluated around m H ∼ 125 GeV: (9) The SM channel weights, ω a , have been defined in (3). The advantage of extracting ( μ * α ) 2 via Eq. (7) over using the experimental values ( μ α ) 2 directly is that it allows for the correlations in the theory uncertainties on the different channel rates to be taken into account. These are correlated to other signals which use the same channels, and since we want to investigate other models beyond the SM, the theory uncertainties on the channel rates are in general different. The same applies for the relative luminosity uncertainties, which can usually be taken equal for all analyses within one collaboration, thus leading to manageable correlations in the signal strength modifiers.
In the next step, we insert these correlated uncertainties into the covariance matrix. To each matrix element (C μ ) αβ , including the diagonal, we add a term ( L αμα )( L βμβ ) if the signals α and β are observed in analyses from the same collaboration (note that usually the further simplification L α = L β applies in this case). We then add the correlated theory uncertainties of the signal rates, given by ⎛ Here, k α and k β are the respective numbers of Higgs (production × decay) channels considered in the experimental analyses where the signals α and β are observed. We use the index notation p(a) and d(a), to map the channel a onto its production and decay processes, respectively. In other words, analyses where the signals share a common production and/or decay mode have correlated systematic uncertainties. These channel rate uncertainties are inserted in the covariance matrix according to their relative contributions to the total signal rate in the model, i.e. via the channel weight evaluated from the model predictions, If the theory uncertainties on the Higgs production and decay rates, as well as the channel weights of the model under investigation, are equal to those in the SM, and also the predicted signal strength matches with the observed signal strength, the uncertainties ( μ α ) 2 extracted from the experimental data are exactly restored for the diagonal elements (C μ ) αα , cf. Eq. (7). Finally, it is worth emphasizing again that this procedure only takes into account the correlations of the luminosity and theoretical signal rate uncertainties, whereas correlations between common experimental uncertainties (energy scale uncertainties, etc.) are neglected. Since this information is not publicly available so far, it could not be included in HiggsSignals.

Higgs mass observables
The other type of observables that give contributions to the total χ 2 in the peak-centered method is the measured masses corresponding to the observed signals. Not all signals come with a mass measurement; this is something which is specified explicitly in the experimental input data. In general, a Higgs boson in the model that is not assigned to a signal (see below for the precise definition), receives a zero χ 2 contribution from this signal. This would be the case, for example, for multiple Higgs bosons that are not close in mass to the observed signal.
HiggsSignals allows the probability density function (pdf) for the Higgs boson masses to be modeled either as a uniform distribution (box), as a Gaussian, or as a box with Gaussian tails. In the Gaussian case, a full correlation in the theory mass uncertainty is taken into account for a Higgs boson that is considered as an explanation for two (or more) signal observables (which include a mass measurement).
Assume that a signal α is observed at the massm α , and that a Higgs boson h i with a predicted mass m i (potentially with a theory uncertainty m i ), is assigned to this signal. Its χ 2 contribution is then simply given by for a uniform (box) mass pdf, and for a box-shaped pdf with Gaussian tails. Here, we denote the experimental uncertainty of the mass measurement of the analysis associated to signal α by m α . The use of a box-shaped mass pdf, Eq. (12), is not recommended in situations where the theory mass uncertainty is small compared to the experimental precision of the mass measurement (and in particular when m i = 0), since this can lead to overly restrictive results in the assignment of the Higgs boson(s) to high-resolution channels. Moreover, a box-shaped pdf is typically not a good description of the experimental uncertainty of a mass measurement in general. We included this option mostly for illustrational purposes.
In the case of a Gaussian mass pdf the χ 2 calculation is performed in a similar way as the calculation of χ 2 μ in Eq. (6). We define for each Higgs boson h i where the αth entry of the predicted mass vector m i is given by m i , if the Higgs boson h i is assigned to the signal α, or m α otherwise (thus leading to a zero χ 2 contribution from this observable and this Higgs boson). As can be seen from Eq. (14), we construct a mass covariance matrix C m i for each Higgs boson h i in the model. The diagonal elements (C m i ) αα contain the experimental mass resolution squared, ( m α ) 2 , of the analysis in which the signal α is observed. The squared theory mass uncertainty, ( m i ) 2 , enters all matrix elements (C m i ) αβ (including the diagonal) where the Higgs boson h i is assigned to both signal observables α and β. Thus, the theoretical mass uncertainty is treated as fully correlated. The sign of this correlation depends on the relative position of the predicted Higgs boson mass, m i , with respect to the two (different) observed mass values,m α,β (where we assumem α <m β for the following discussion): If the predicted mass lies outside the two measurements, i.e. m i <m α ,m β or m i >m α ,m β , then the correlation is assumed to be positive. If it lies in between the two mass measurements,m α < m i <m β , the correlation is negative (i.e. we have anti-correlated observables). The necessity of this sign dependence can be illustrated as follows: Let us assume the predicted Higgs mass is varied within its theoretical uncertainty. In the first case, the deviations of m i from the theoretical mass uncertaintiesm α,β both either increase or decrease (depending on the direction of the mass variation). Thus, the theoretical mass uncertaintines are positively correlated. However, in the latter case, a variation of m i towards one mass measurement always corresponds to a larger deviation of m i from the other mass measurements. Therefore, the theoretical mass uncertainties for these observables have to be anti-correlated.

Assignment of multiple Higgs bosons
If a model contains an extended (neutral) Higgs sector, it is a priori not clear which Higgs boson(s) give the best explanation of the experimental observations. Moreover, possible superpositions of the signal strengths of the Higgs bosons have to be taken into account. Another (yet hypothetical) complication arises if more than one Higgs signal has been discovered in the same Higgs search, indicating the discovery of another Higgs boson. In this case, care has to be taken that a Higgs boson of the model is only considered as an explanation of one of these signals.
In the peak-centered χ 2 method, these complications are taken into account by the automatic assignment of the Higgs bosons in the model to the signal observables. In this procedure, HiggsSignals tests whether the combined signal strength of several Higgs bosons might yield a better fit than the assignment of a single Higgs boson to one signal in an analysis. Moreover, based on the predicted and observed Higgs mass values, as well as their uncertainties, the program decides whether a comparison of the predicted and observed signal rates is valid for the considered Higgs boson. A priori, all possible Higgs combinations which can be assigned to the observed signal(s) of an analysis are considered. If more than one signal exists in one analysis, it is taken care of that each Higgs boson is assigned to at most one signal to avoid double-counting. A signal to which no Higgs boson is assigned contributes a χ 2 penalty given by Eq. (6) with the corresponding model prediction μ α = 0. This corresponds to the case where an observed signal cannot be explained by any of the Higgs bosons in the model.
For each Higgs search analysis the best Higgs boson assignment is found in the following way: For every possible assignment η of a Higgs boson combination to the signal α observed in the analysis, its corresponding tentative χ 2 contribution, χ 2 α,η , based on both the signal strength and potentially the Higgs mass measurement, is evaluated. In order to be considered for the assignment, the Higgs combination has to fulfill the following requirements: • Higgs bosons which have a mass m i close enough to the signal massm α , i.e.
are required to be assigned to the signal α. Here, denotes the assignment range, which can be modified by the user, see Sect. 4.4 (the default setting is = 1). • If the χ 2 contribution from the measured Higgs mass is deactivated for this signal, combinations with a Higgs boson that fulfills Eq. (15) are taken into account for a possible assignment, and not taken into account otherwise. • If the χ 2 contribution from the measured Higgs mass is activated, combinations with a Higgs boson mass which does not fulfill Eq. (15) are still considered. Here, the difference of the measured and predicted Higgs mass is automatically taken into account by the χ 2 contribution from the Higgs mass, χ 2 m .
In the case where multiple Higgs bosons are assigned to the same signal, the combined signal strength modifier μ is taken as the sum over their predicted signal strength modifiers (corresponding to incoherently adding their rates). The best Higgs-to-signals assignment η 0 in an analysis is that which minimizes the overall χ 2 contribution, i.e.
Here, the sum runs over all signals observed within this particular analysis. In this procedure, HiggsSignals only considers assignments η where each Higgs boson is not assigned to more than one signal within the same analysis in order to avoid double counting. There is also the possibility to enforce that a collection of peak observables is either assigned or not assigned in parallel. This can be useful if certain peak observables stem from the same Higgs analysis but correspond to measurements performed for specific tags or categories (e.g. as presently used in H → γ γ analyses). See Sect. 4.6 for a description of these assignment groups.
A final remark should be made on the experimental resolution, m α , which enters Eq. (15). In case the analysis has an actual mass measurement that enters the χ 2 contribution from the Higgs mass, m α gives the uncertainty of the mass measurement. If this is not the case, m α is an estimate of the mass range in which two Higgs boson signals cannot be resolved. This is taken to be the mass resolution quoted by the experimental analysis. Typical values are, for instance, 10 % (for V H → V (bb) [85]) and 20 % (for H → τ τ [86] and H → W W ( * ) → ν ν [87]) of the assumed Higgs mass. It should be kept in mind that the HiggsSignals procedure to automatically assign (possibly several) Higgs bosons to the signals potentially introduces sharp transitions from assigned to unassigned signals at certain mass values, see Sect. 5.1.1 for a further discussion. More detailed studies of overlapping signals from multiple Higgs bosons, where possible interference effects are taken into account, are desirable in case evidence for such a scenario emerges in the future data.

The mass-centered χ 2 method
The mass-centered χ 2 method is complementary to the peakcentered χ 2 method, since it allows for a more general test of the model against the experimental data without reference to particular signals. This method uses the data where the measured best-fit signal strength modifiers are published as a function of the Higgs mass over the (full) investigated mass range, as shown in Fig. 1a. 3 A χ 2 test can then be performed directly at the predicted Higgs mass(es), m i , of the model if these fall within the experimentally investigated mass range of an analysis a (denoted by G a ). For Higgs bosons that are outside this mass range, HiggsSignals provides no information. Also in this method, like in the peak-centered case, it can be necessary to consider the combined rates of several Higgs bosons which are close in mass compared to the experimental resolution. We begin with a general discussion of the single Higgs (non-mass-degenerate) case, and outline the combination scheme below.

Theory mass uncertainties
In theμ plot the experimental mass uncertainty is already taken into account in the experimental analysis. However, we also want to take into account a possible theoretical uncertainty on the predicted Higgs mass, m i . Higgs-Signals provides two different methods to include theoretical Higgs mass uncertainties in the mass-centered χ 2 evaluation: (i) (default setting) In the first method the predicted Higgs mass is varied around m i within its uncertainties. We denote this varied mass by m in the following. For a uniform (box) parametrization of the theoretical mass uncertainty, we have the allowed mass range A tentative χ 2 distribution is evaluated as a function of m , which, in the uniform (box) parametrization, takes the form For the Gaussian parametrization, we have In these expressions, n denotes the total number of considered analyses. Note that the predicted signal strengths, μ a , are always calculated at the predicted central values for the Higgs mass, m i , (from the user input), and the signal strength is held fixed in the mass variation. This is clearly an approximation, but for small theory mass uncertainties m i it is reasonable to treat resulting variations in μ as a second-order effect. 4 From a practical viewpoint, it also reduces significantly the amount of model information that has to be supplied by the user.
The final values forμ and μ are chosen for each Higgs boson h i at the mass value m 0 i = m , where χ 2 i (m ) is minimized (i.e. for each Higgs boson separately, but combining all channels). In this way, the most conservative value of the predicted Higgs mass, within its theory uncertainty, is used to define the measured signal strength modifiers for the final χ 2 evaluation. (ii) In the second approach to include theory mass uncertainties, HiggsSignals convolves the experimentally measured signal strength modifier,μ a (m), with a theory mass pdf, g(m , m), resulting in The theory mass pdf g(m , m) can again be chosen to be either a uniform (box) distribution or a Gaussian, both centered around the predicted mass value, m, and with a box width of ± m or a Gaussian width m, respectively. The pdf is normalized to unity over the mass range G a in order to preserve probability. In the case of zero theoretical Higgs mass uncertainty, 5 The model prediction is therefore tested directly against the measured valueμ(m) at the predicted (exact) value for the mass m.
The observed signal strength modifier after convolution, μ conv a , now includes contributions to the measured signal strength modifier from the mass region close to the predicted Higgs mass (weighted by g(m , m)). Similarly, the upper and lower experimental 1σ uncertainty (cyan) band values, μ a , are smeared In this case it is the smeared quantities, evaluated from Eqs. (20) and (21), that enter the χ 2 test.

The Stockholm clustering scheme
If more than one neutral Higgs boson of the model has a mass in the relevant region of an analysis, m i ∈ G a , possible superpositions of their signal rates have to be taken into account without double-counting. In order to determine the relevant combinations (out of the potentially many options), we use a prescription inspired by jet clustering. In a similar spirit, we call this the Stockholm clustering scheme: 1. Determine the nearest neighboring Higgs bosons h i and h j by their mass difference m i j = |m i − m j |. If min( m i j ) is larger than the experimental mass resolution of the analysis, the clustering is finished, and we proceed to step 4. If it is smaller, the two Higgs bosons h i and h j will be clustered (combined).

The combination of two adjacent Higgs bosons h i and
h j defines a new Higgs cluster h k with the following properties: • If both Higgs bosons h i and h j have non-zero theoretical mass uncertainties ( m i = 0 and m j = 0) the combined mass is obtained from a Gaussian average (regardless of the choice for Higgs mass pdf), with the combined theoretical mass uncertainty • If either m i or m j is known exactly, for instance m i = 0, the mass of the new Higgs cluster is chosen equal to this mass, m k = m i , with zero combined theory mass uncertainty, m k = m i = 0.
• If both m i and m j are known exactly, m i = m j = 0, the Higgs cluster is assigned an averaged mass 3. The procedure is repeated from step 1. The entities considered for further clustering include both the unclustered (initial) Higgs bosons, as well as the already combined Higgs clusters. The single Higgs bosons which form part of a cluster are no longer present. 4. Each single Higgs boson or Higgs cluster h k that remains after the clustering according to steps 1-3 enters the mass-centered χ 2 test. Their predicted signal strength modifiers are formed from the incoherent sum (again, neglecting interference effects) of the individual signal strength modifiers for the combined Higgs bosons, In this way, the predictions that are compared to one implemented analysis are determined. HiggsSignals repeats this procedure for all implemented experimental analyses.
Since the experimental mass resolution can vary significantly between different analyses, the resulting clustering in each case may differ.
The two different treatments of the theoretical mass uncertainties, as discussed above, have to be slightly extended for the case of Higgs clusters: (i) If the Higgs boson h i is contained within a Higgs cluster h k for one analysis a, the considered mass region for the variation of m in (18) in the case of a uniform (box) Higgs mass pdf. 6 We denote the resulting tentative total χ 2 from the variation of the mass of Higgs boson h i by χ 2 i . The variation is done for every Higgs boson contained in the cluster h k . When the cluster h k is evaluated against the observed results for analysis a, the observed valuesμ a and μ a are defined at the value of m where the global χ 2 , composed of all χ 2 i distributions, is minimal. 7 (ii) In the second approach, the convolution of the experimentalμ values with theory uncertainties is performed separately for each Higgs boson, or Higgs cluster k, with the combined Higgs mass pdf The normalization factor N = M k dm g k (m , m) to preserve probability. The sum runs over all Higgs bosons which have been combined for this cluster.
Once all model predictions and mass-centered observables have been defined, when necessary using Stockholm clustering as discussed above, the total mass-centered χ 2 is evaluated with a signal strength vector 8 and covariance matrix constructed analogously as in the peak-centered χ 2 method, cf. (6). The uncertainties of production cross sections, decay rates, and the luminosity are again treated as fully correlated Gaussian errors. Note that, in this method, there is no contribution from Higgs mass measurements to the total χ 2 , since the evaluation is done directly against the experimental data at the predicted Higgs mass values (within their uncertainties).
As a final remark, we would like to point out that theμ plots necessary for this method are so far only published for a 6 If M i ∩ M k = ∅, we increase M k until there is a (minimal) overlap. This will effectively lead to an evaluation of the tentative χ 2 at the boundary of M i which is closest to the mass m k of the Higgs cluster. 7 The global χ 2 is defined in the mass region few selected analyses. 9 Thus, there is not (yet) a full coverage of the various Higgs signal topologies with the mass-centered χ 2 method. Furthermore, the published results cover only a limited range in the Higgs mass, which is a further limit to its applicability.

Simultaneous use of both methods
Since the two methods presented here are complementarythey test inherently different statistical hypotheses-HiggsSignals allows for the possibility to apply the peak-centered and mass-centered χ 2 methods simultaneously. We present here one approach, which attempts to make maximal use of the available experimental information when testing models with multiple Higgs bosons. The user of HiggsSignals is of course free to use other combinations of the two results, which can be derived completely independently.
In the provided combined approach, HiggsSignals first runs the peak-centered χ 2 method and assigns the Higgs bosons to the observed signals, tracing the assigned combination for each analysis. In the second step, all remaining Higgs bosons (which have not been assigned) are considered with the mass-centered χ 2 method; their respective (masscentered) χ 2 contributions are constructed. In this way, a possible double-counting, where a Higgs boson is tested with both the peak-and mass-centered χ 2 method against the same data, is avoided. In the last step, the total χ 2 is evaluated. Here, the Higgs mass χ 2 from the (relevant) signals, as well as the χ 2 from combined signal strength vectors from both the peak-centered and the mass-centered approach, are evaluated with a full covariance matrix. This method thus tests the model predictions against the data in the maximal possible way, while ensuring that no Higgs boson is tested more than once against the same experimental data.
As a final recommendation, it should be noted that the mass ranges for the measuredμ values are still much smaller than the mass ranges for (SM) Higgs exclusion limits. To constrain theories with Higgs bosons outside this smaller range (or below the lower limit of the range currently considered by LHC searches), it is still highly recommended to run HiggsBounds [13][14][15] in parallel to HiggsSignals.

Installation
The latest version of HiggsSignals can be downloaded from the webpage http://higgsbounds.hepforge.org which is 9 Currently, theμ plots are published only for the H → γ γ , H → Z Z ( * ) and H → W W ( * ) searches. also the home of HiggsBounds. Since HiggsSignals depends on the HiggsBounds libraries, this code (version 4.0.0 or newer) should be downloaded and installed as well. For further detail on how to do this, we refer to the Higgs-Bounds manual [13][14][15]. Like HiggsBounds, Higgs-Signals is written in Fortran 90/2003. Both codes can be compiled, for example, using gfortran (version 4.2 or higher). After unpacking the downloaded source files, which should create a new directory for HiggsSignals, the user possibly needs to set the correct path to the Higgs-Bounds installation in the configure file. Optionally, the path to a FeynHiggs installation (version 2.9.4 or higher recommended) [88][89][90][91][92] can be set in order to use some of the example programs which use FeynHiggs subroutines (see below). Furthermore, compiler flags necessary for specific platforms can be placed here. Configuration and installation starts with running ./configure which will generate a makefile from the initial file makefile.in. Once this is done, run make to produce the HiggsSignals Fortran library (called libHS.a) and the command line executable. In addition, the user may conveniently use a bash script,

./run_tests.bat
to build the HiggsSignals library and executable as well as the provided example programs (described in Sect. 4.5). The script will then perform a few test runs.

Input and output
HiggsSignals is designed to require mostly the same input as HiggsBounds, so that users already familiar with this code should be able to transfer their existing analyses to also use HiggsSignals with a minimal amount of extra work. There are two ways to run HiggsSignals: either from the command line, or via the subroutines contained in the HiggsSignals library libHS.a. For the command line version, the model predictions (Higgs masses, their theory uncertainties, total widths, production and decay rates) have to be specified in data files using the same format as HiggsBounds-4, see Ref. [15]. The command line version of HiggsSignals is presented in more detail in Sect. 4.3.
In the subroutine version, the model predictions (which can be given as effective couplings, or as cross sections either at partonic or hadronic level) have to be provided via subroutines. Most of these subroutines are shared with the HiggsBounds library (for details we refer again to [15]). In addition to the HiggsBounds input, HiggsSignals requires input of the theoretical uncertainties on both the Higgs masses and the rate predictions. Therefore, Higgs-Signals contains two additional input subroutines to set these quantities, see Sect. 4.4 for more details. An accessible demonstration of how to use the HiggsSignals subroutines is provided by the example programs, discussed further in Sect. 4.5.
As already mentioned, the required input of Higgs production and decay rates can be given either as effective couplings, or as cross sections at partonic or hadronic level. For supersymmetric models there is an option of using the SUSY Les Houches Accord (SLHA) [93,94] for input (either using data files or subroutines). In this case, the production rates are always approximated using the effective couplings specified in the two HiggsBounds specific input SLHA blocks (as specified in Ref. [15]), whereas the Higgs branching ratios are taken directly from the corresponding decay blocks. If present, the theoretical mass uncertainties are read in from the SLHA block DMASS (as available e.g. from FeynHiggs). Otherwise, since there is no consensus yet on how to encode the theoretical rate uncertainties in the SLHA format, these have to be given to HiggsSignals explicitly by hand. 10 The main results from HiggsSignals are reported in the form of a χ 2 value and the number of considered observables. For reference, the code also calculates the p value associated to the total χ 2 and the number of degrees of freedom N . The user may specify the number of free model parameters N p (see below). Then, the number of degrees 10 This can be done by either calling the subroutine setup_ rate_uncertainties (see below) or by including the rate uncertainties directly in the file usefulbits_HS.f90 in case the subroutine cannot be used (i.e. if HiggsSignals is run on the command line). If the user does not specify the rate uncertainties (in either case), they are assumed to be identical to the SM rate uncertainties, Eq. (9). of freedom is given by N = N obs − N p , where N obs is the total number of the included observables. Note that if the user does not specify N p , the p value is evaluated assuming N p = 0.
In the case of running with input data files, the Higgs-Signals output is written into new files as described in Sect. 4.3. There also exist subroutines, see Sect. 4.4, to specify the extent of screen output and to retrieve many quantities of interest for further analysis.
If HiggsSignals is run in the SLHA mode, the results can be appended to the SLHA file in the form of new SLHAinspired 11 blocks. The main results are then collected in BLOCK HiggsSignalsResults, as shown for a specific example in Table 1. The first entries of this BLOCK contain general information on the global settings of the HiggsSignals run, i.e. the version number, the experimental data set, the χ 2 method and the Higgs mass parametrization used. Moreover, it lists the number of analyzed observables of the different types (BLOCK entries 4-6), as well as the total number (BLOCK entry 7). Next, it gives the corresponding χ 2 values separately from the signal strength peak observables (BLOCK entry 8), the Higgs mass peak observables (BLOCK entry 9), and the mass-centered observables (BLOCK entry 10). The total signal strength χ 2 for both methods (the sum of BLOCK entries 8 and 10) is provided (BLOCK entry 11), as is the total χ 2 sum (BLOCK entry 12). The final element (BLOCK entry 13) gives the reference p value, as discussed above.
Additional output specific to the peak-centered χ 2 method is collected in BLOCK HiggsSignalsPeakObservables.
We show an excerpt from this extensive BLOCK for an example (MSSM) parameter point in Table 2. The first identifier, OBS, in the BLOCK enumerates the peak observables, whereas the second number, FLAG, labels the specific quantity (for this peak observable). For every peak observable, the first entries (FLAG=1-11) give general information about the experimental data defining the observable. This is followed by model-specific information and the results from the HiggsSignals run. FLAG=12 displays a binary code representing the Higgs boson combination which has been assigned to the signal. It has the same length as the number of Higgs bosons, 12 such that an assigned Higgs boson with index k corresponds to the binary value 2 k−1 . A code of only zeroes means that no Higgs boson has been assigned to this peak observable. In the specific example shown in Table 2, the lightest of the three neutral Higgs bosons in the MSSM (with k = 1) has been assigned.
This BLOCK also contains additional information (index i, Particle data group (PDG) number, mass, and signal strength contribution under FLAG=13-16) about the assigned Higgs boson that gives the largest contribution to the total predicted signal strength. The total predicted signal strength is given by FLAG=17. The HiggsSignals results (FLAG=18-20) contain the χ 2 contribution from the signal strength and Higgs mass test from this observable, as well as the total χ 2 contribution obtained for the assigned Higgs boson combination. Finally, the χ 2 obtained for the case with no predicted signal, μ = 0, is given for FLAG=21. It should be noted that the quoted χ 2 values correspond to intermediate results in the total χ 2 evaluation, where correlated uncertainties are taken into account by the covariance matrix. For instance, the signal strength χ 2 (FLAG=18) corresponds to χ 2 μ,α in Eq. (6), where α is the index of the peak observable given in the first column of the BLOCK. Thus, this quantity differs from the naïvely calculated χ 2 = (μ −μ) 2 /( μ) 2 , and might in the extreme case even be negative due to the impact of correlated uncertainties.
The results from the mass-centered χ 2 method are summarized in

BLOCK HiggsSignalsMassCenteredObservables
in a similar way as in BLOCK HiggsSignalsPeak Observables. An example is given in Table 3. The modelindependent information about the observable (FLAG=1-7) is identical to the corresponding information in BLOCK HiggsSignalsPeakObservables. However, since the evaluated experimental quantities of the mass-centered observable depend on the model prediction, cf. Sect. 3.2, we give the information of the tested Higgs boson (cluster) at first (FLAG=8-10), corresponding to Eqs. (22)- (24). The number and binary code of the combined Higgs bosons, which form a Stockholm Higgs cluster, is given by FLAG=11 and 12, respectively. From the experimental data is given the mass position (FLAG=13), and the measured signal strength with its lower and upper uncertainties (FLAG=14-16). Finally, the resulting χ 2 con-tribution from this mass-centered observable is given at FLAG=17.
Note that there is also the possibility to create a new SLHA file with the HiggsSignals output blocks even if the input was not provided in SLHA format. Moreover, HiggsSignals can give an extensive screen output with similar information as encoded in the three SLHA output blocks. The level of information that is desired should then be specified before the HiggsSignals run via the subroutine setup_output_level. See Sect. 4.4 for more details.

Running HiggsSignals on the command line
HiggsSignals can be run on the command line as follows: This command line call is very similar to the one of Higgs-Bounds and the last four arguments have been directly taken over from HiggsBounds. The user may consult the HiggsBounds manual [15] for more details on these arguments. The number of neutral and charged Higgs bosons of the model are specified by nHzero and nHplus, respectively. As in HiggsBounds, the model predictions are read in from the data files specified by prefix. Which data files are required as input depends on the argument whichinput, which can take the string values effC, part, hadr and SLHA for the various input formats. The theory mass uncertainties are read in from the data file <prefix>MHall_uncertainties.dat for both the neutral and charged Higgs bosons. If this file is absent these uncertainties are set to zero. For more information of the data file structure we refer to the HiggsBounds-4 manual [15]. Note that for whichinput=SLHA, all the input is read in from the SLHA input file which, like the ordinary data files, should be specified by <prefix>.
The first three arguments are intrinsic HiggsSignals options. The string < expdata > specifies which experimental data set should be used. HiggsSignals will read in the observables found in the directory Expt_tables/ < expdata >. The second argument, < mode >, specifies which χ 2 method should be used; it can take the string values peak (for the peak-centered χ 2 method, described in Sect. 3.1), mass (for the mass-centered χ 2 method, see Sect. 3.2), or both (for the simultaneous use of both methods, as described in Sect. 3.3). Finally, the < pdf > argument takes an integer selecting the parametrization for the Higgs mass uncertainty as either 1 (box), 2 (Gaussian), or 3 (box+Gaussian) pdf.
As an example, the user may run ./HiggsSignals latestresults peak 2 effC 3 1 example_data/mhmax/mhmax_ which runs the peak-centered χ 2 method on the provided parameter points in the (M A , tan β) plane of the m max h benchmark scenario [95] of the MSSM, using the most recent Higgs data contained in the directory Expt_tables/ latestresults/.
The HiggsSignals output from a successful command line run is collected in the data file < prefix > HiggsSignals_results.dat, except for the case whichinput=SLHA, where the results are attached as SLHA output blocks to the SLHA file, cf. Sect. 4.2. The SUSY spectrum generator SPheno [96,97], used in conjunction with the model building tool SARAH [98][99][100], can write directly the HiggsBounds (and thus Higgs-Signals) data files for input in the effective couplings format.

HiggsSignals subroutines
In this section we present the subroutines needed for the use of HiggsSignals. First, we go step-by-step through the user subroutines encountered during a normal run of HiggsSignals. Then, we list additional (optional) subroutines for specific applications of HiggsSignals, and for a convenient handling of the output.

Main user subroutines
The subroutine that is usually called first is which sets up the HiggsSignals framework: It allocates internal arrays according to the number of neutral (nHzero) and charged 13 (nHplus) Higgs bosons in the model and reads in the tables for the SM branching ratios in the same way as done in HiggsBounds. Furthermore, it calls the subroutine setup_observables, which reads in the experimental data contained in the directory Expt_tables/(expdata). The user may create a new directory in Expt_tables/ containing the relevant observables for his study, see Sect. 4.6 for more details. For convenience, we also provide a wrapper subroutine which does not require the third argument but uses the experimental data from the folder Expt_tables/ latestresults/. The next step is to specify the probability density function (pdf) for the Higgs masses, which is done using setup_pdf. Available settings are pdf = 1 for a uniform (box-shaped) distribution, pdf = 2 for a Gaussian, and pdf = 3 for a box-shaped pdf with Gaussian tails. The impact of this choice has been discussed in detail in Sect. 3 and will furthermore be demonstrated in Sect. 5. With the subroutine values for the theory mass uncertainties m i can be specified. This subroutine sets the theoretical uncertainties of the neutral Higgs boson masses (in GeV) of the model via the array dMh. The default values (in case this subroutine is not invoked) is for all uncertainties to be zero. Note that HiggsBounds-4 also contains a similar subroutine (set_mass_uncertainties) to set theoretical Table 4 Ordering of the elements of the input arrays dCS and dBR for the relative uncertainties of the hadronic production cross sections and branching ratios, respectively Array Element Recall that the hadronic production mode "singleH" usually contains both the partonic processes gg → H and bb → H , currently assuming equal experimental efficiencies. The latter can change in the future once search categories with b-tags are included. This table will possibly be extended once measurements in new channels (e.g. H → Z γ ) are performed mass uncertainties of the neutral and charged Higgs bosons. These uncertainties are taken into account via mass variation in the HiggsBounds run. Since the treatment of these uncertainties is intrinsically different between the two codes, we allow the user to set the theoretical mass uncertainties for HiggsSignals independently using this subroutine. 14 For models with different uncertainties on the Higgs production cross sections and branching ratios than those for a SM Higgs boson, these should be specified using this subroutine, which sets the theoretical uncertainties of the production and decay rates (in %) in the considered model. In the current implementation, LHC and Tevatron channels are considered to have the same relative rate uncertainties, and the rate uncertainties are assumed to be the same for all neutral Higgs bosons, independent of their masses. The input arrays should follow the structure of Table 4.
The remaining required input (Higgs boson masses, total widths, branching ratios, cross sections) is identical to the HiggsBounds input and should be set via the Higgs-Bounds input subroutines, cf. Ref. [15].
In order to evaluate a meaningful p value during the HiggsSignals run, the program has to know the number of free model parameters, N p , cf. Sect. 4.2. This number is specified by the subroutine setup_nparam. If this subroutine is not called before the main HiggsSignals run, the code assumes no free model parameters, N p = 0.
14 The use of different theoretical mass uncertainties in Higgs-Bounds and HiggsSignals is restricted to the subroutine version. In the command line version of both programs, the theoretical uncertainties will be read in from the same data file, namely < prefix >MHall_uncertainties.dat.
Once all the input has been specified, the main Higgs-Signals evaluation can be run by calling the run _HiggsSignals subroutine to start the χ 2 evaluation. The mode flag specifies the χ 2 method which is used in the following evaluation process. Possible values are mode = 1 (peak-centered method, cf. Sect. 3.1), mode = 2 (masscentered method, cf. Sect. 3.2), or mode = 3 (simultaneous use of both methods, cf. Sect. 3.3). After a successful run, this subroutine returns the χ 2 contribution from the signal strength measurements (csqmu), 15 the χ 2 contribution from the Higgs mass measurements (csqmh), and the total χ 2 value (csqtot). It also returns the number of observables involved in the χ 2 evaluation (nobs). If the mass-centered χ 2 method is employed, it is important to realize that nobs can depend on many parameters, such as the Higgs boson masses of the model (which may be inside or outside the range of an analysis). The Stockholm clustering can also affect the number of observables that are evaluated in the final χ 2 calculation. Finally, the associated p value (Pvalue) for the total χ 2 with nobs−N p degrees of freedom is calculated.
At the end of a HiggsSignals run, the user should call this routine to deallocate all internal arrays.

Specific user subroutines
This section provides a list (alphabetically ordered) of subroutines handling more special features of HiggsSignals.
If the user wants to perform a dedicated statistical study using pseudo-measurements (also called toy-measurements) for the Higgs signal rates and mass measurements, they can be set via this subroutine for the peak observable with the identification number obsID. This observable ID is unique to the peak observable and is encoded in the experimental data, see Sect. 4.6 for more details. After a (dummy) run of HiggsSignals the observable ID can also be read out with the subroutine get_ID_of_peakobservable (see below). The arguments mu_obs and mh_obs are the pseudo-measured values for the signal strength modifierμ and the Higgs massm. Note that the uncertainties are kept at their original values.
If the user wants to scale the uncertainties of the Higgs signal rate and mass measurements, this can be done via this subroutine in an analogous way as setting the toy measurements (using assign_toyvalues_to_peak).
Here, scale_mu is the scale factor for the experimental uncertainty on the signal strength of the peak with identification number obsID. The theoretical rate uncertainties, which can be set independently via the subroutine setup_rate_uncertainties (see above), are unaffected by this scale factor. In this way, Higgs-Signals allows the user to scale the experimental and theoretical rate uncertainties independently. This is useful if the user is interested in a future projection of the compatibility between the model and the experimental data, assuming that a certain improvement in the precision of the measurements and/or theoretical predictions can be achieved.
After the HiggsSignals run the user can employ the following "get_" subroutines to obtain useful information from the HiggsSignals output. The following three subroutines are contained in the Fortran module io.
If the peak-centered χ 2 method is used, the peak observables are internally enumerated in HiggsSignals based on their alphabetical appearance in the directory Expt_ tables/(expdata) of the used experimental dataset. This ordering is reflected e.g. in the screen output and the SLHA output. However, a safer way to access the peak observables (for instance to set toy observables) is to use the unique observable ID of the peak observable. For this, the user may call this subroutine which returns the observable ID obsID internally structured at the position i.
This subroutine returns the total number of various observables: ntotal is the total number of observables, npeakmu and npeakmh are the number of signal strength and Higgs mass observables entering the peak-centered χ 2 method, respectively, nmpred is the number of observables consid- More information about the HiggsSignals result can be obtained by calling this subroutine. It returns the total predicted signal strength modifier, the index of the dominantly contributing Higgs boson and the number of combined Higgs bosons for the peak observable with observable identifier obsID as mupred, npeak and nHcomb, respectively.
The user may apply the subroutine get_Pvalue to evaluate the p value again after run_HiggsSignals, with the possibility to vary N p . The result is based on the total χ 2 and the total number of observables from the last Higgs-Signals run as well as the number of free parameters, Np, which are passed as input to this subroutine.
This subroutine allows the user to read out the predicted signal rate for an arbitrary channel combination. This channel combination is specified by the number of combined channels, Nchannels, and the array IDchannels, which contains the two-digit IDs of these channels as specified in cf. Table 5. The output (rate) is the combined rate. It is more general than get_Rvalues (see below). This returns the model-predicted signal rates (normalized to the SM signal rates) of Higgs boson i for the six different processes listed in Table 6. These signal rates are calculated via Eq. (1), assuming that all channels have the same relative If the user does not use the SLHA input format of Higgs-Signals, or rather wants to write the output into a different file, this subroutine can be used to create a new file as specified by the argument filename. If this file already exists, HiggsSignals will not overwrite this file but give a warning. The integer argument detailed takes values of 0 or 1, determining whether only the block HiggsSignalsResults or all possible output blocks (i.e. also the block HiggsSignalsPeakObservables and/or HiggsSignalsMassCenteredObservables), respectively, are written to the file. The wrapper subroutine does the same but for the default filename called HSoutput.slha.
If HiggsSignals is run on an SLHA input file, the subroutine HiggsSignals_SLHA_output appends the HiggsSignals results as blocks to the SLHA input file.
The following "setup_" subroutines can be used to change the default settings of the HiggsSignals run. Thus, they should be called before the subroutine run_ HiggsSignals.
This subroutine can be used to change the mass range, in which a Higgs boson is forced to be assigned to a peak observable, see Sect. 3.1.3. The value Lambda corresponds to in Eq. (15).
The subroutine can be used to switch off (on) the correlations among the systematic uncertainties in the χ 2 evaluation of the signal strength [Higgs mass] part by setting corr_mu [corr_mh] = 0 (1). If this subroutine is not called, the default is to evaluate the χ 2 with correlated uncertainties (corr_mu = corr_mh = 1).
If the mass-centered χ 2 method is used, the treatment of the Higgs mass theory uncertainty can be set by calling this subroutine with mode=1 to use the mass variation (default), or mode=2 for convolving the theory mass uncertainty with theμ plot. See Sect. 3.2 for more details of these methods.
The user may control the screen output from the Higgs-Signals run with the subroutine, where level takes values from 0 to 3, corresponding to the following output: 0 Silent mode (suitable for model parameter scans, etc.) (default), 1 Screen output for each analysis with its peak and/or masscentered observables. The channel signal strength modifiers and SM channel weights, cf. Eqs. (2) and (3), respectively, are given for all channels considered by the analysis. 2 Screen output of the essential experimental data of the peak observables and/or implementedμ plots (as used for the mass-centered χ 2 method). For each observable, the signal channels are listed with the implemented efficiencies. 3 Creates text files holding essential information about the experimental data and the model predictions for each observable. In the peak-centered χ 2 run mode, the files peak_information.txt and peak_ massesandrates.txt are created. The first file lists all peak observables, including a description and references to the publications, whereas the second file gives the observed and model-predicted values for the Higgs mass 16 and signal rates and their corresponding pull values, which we define as: pull value = predicted value − observed value (Gaussian combined) uncertainty (26) Note that in this expression the effect of correlated uncertainties is not taken into account. In the mass-centered χ 2 run mode, the files mctables_information.txt and mcobservables_information.txt are created. The first file gives general information about the analyses with an implementedμ-plot. The second file lists all mass-centered observables, which have been constructed during the HiggsSignals run, including the mass position, the observed and predicted signal strength values as well as their pull values. The third example program, HSwithSLHA, uses the SLHA input of HiggsBounds, i.e. an SLHA file which contains the two special input blocks for HiggsBounds. It can be executed with ./HSwithSLHA < number of SLHA files > < SLHAfilename > The program can test several SLHA files in one call. The total number of SLHA files must therefore be given as the first argument. The SLHA files must all have the same name, and should be enumerated by SLHA_filename.x, where x is a number. Running, for example, ./HSwithSLHA 2 SLHAexample.fh would require the two SLHA files SLHAexample.fh.1 and SLHAexample.fh.2 to be present. The output is written as SLHA blocks, cf. Sect. 4.2, which are appended to each input SLHA file. The example program HBandHSwithSLHA can be run in an analogous way. It employs both Higgs-Bounds and HiggsSignals on the provided SLHA file(s), demonstrating how these two codes can be run together efficiently.
The example program HSwithToys demonstrates how to set new values (corresponding to pseudo-measurements) forμ andm for each signal. In the code, HiggsSignals is first run on the SM with a Higgs mass around 126 GeV using the effective couplings input format. Then, the predicted signal strengths are read out from the HiggsSignals output and set as pseudo-measurements. A second Higgs-Signals run on these modified observables then results in a total χ 2 of zero.
The example program HS_scale_uncertainties also runs on the SM with a Higgs mass around 126 GeV. It scans over a universal scale factor for (i) the experimental uncertainty of the signal strengthμ only, (ii) the theoretical uncertainties of the production cross sections and branching ratios only, and finally (iii) both experimental and theoretical uncertainties. The output of each scan is saved in text files. In this way, rough projections of the model compatibility to a more accurate measurement in the future (with the same central values) can be made.
The last example, HBandHSwithFH, demonstrates how to run HiggsBounds and HiggsSignals simultaneously on a realistic model, in this case the MSSM. Here, FeynHiggs [88][89][90][91] is used to calculated the MSSM pre- . For a detailed description of each line in the file, see Table 9 dictions needed as input for HiggsBounds and Higgs-Signals.

Input of new experimental data into HiggsSignals
The ambition with HiggsSignals is to always keep the code updated with the latest experimental results. Nevertheless, there are several situations when a user may want to manually add new data (or pseudo-data) to the program, for example to assess the impact of a hypothetical future measurement. For advanced users, we therefore provide a full description of the data file format used by HiggsSignals.
For each observable that should be considered by Higgs-Signals, there must exist a textfile (file suffix: .txt). This file should be placed in a directory

Expt_tables/(expdata)/
where (expdata) is the name identifying the new (or existing) experimental dataset. 17 All analysis files in this directory will then be read in automatically by HiggsSignals during the initialization.
As an example we show in Tables 7 and 8 the two data files for the inclusive measurement of the ATLAS H → Z Z ( * ) → 4 analysis [75], which define a peak observable and provide the fullμ plot as needed by the mass-centered χ 2 method, respectively. The first 11 rows of these files encode general information about the analysis and the observable (each row is required), as described in Table 9. Comments can be included in the top rows if they are starting with a # symbol. Note that the observable ID must be unique, 17 The identifier (expdata) is the argument which has to be passed to initialize_HiggsSignals at initialization, cf. Sect. 4.4. This file is located in Expt_tables/latestresults-1.0.0_ inclusive/ (with name ATL_H-ZZ-4l_7-8TeV_4.6fb-1_ 20.7fb-1_2013013201.txt). It is the same analysis for which we already defined a peak observable in Table 7. For a detailed description of each line in the file, see Table 9 whereas the analysis ID must be the same for (peak-or masscentered) observables, which correspond to the same analysis and where a multiple assignment of the same Higgs boson to the corresponding observables shall be avoided. In the (yet hypothetical) case that two distinct signals have been observed within the same analysis, their peak observables thus need to have the same analysis ID, otherwise a Higgs boson might be assigned to both signals. All integers should not have more than 10 digits. The channel codes in the 10th row are given as two-digit integers, where the first digit encodes the production mode, and the second digit the decay mode. The corresponding numbers are given in Table 5. For example, the channel code of ( pp) → H W → (bb)W is 35. In the example of Table 7, we thus consider all five production modes, but only a single decay mode, i.e. H → γ γ .
Channel efficiencies can be included in the 11th row. They correspond to the channels as defined by the channel codes on the previous row, and thus have to be given in the same order. If the experimental channel efficiencies are unknown (as in the given example of an inclusive measurement), the reference mass in the 9th row should be set equal to −1, in which case the 11th row will be ignored. Since it must still be present, it could be left blank for the sake of clarity. Note that the channel efficiencies are defined as the fraction of events passing the analysis cuts, and not the relative contribution of this channel to total signal yield. The latter would Number of search channels, reference mass for efficiencies (−1: no efficiencies given) 10 Search channel codes (see Table 5) (# entries must equal # channels)) 11 Channel efficiencies (# entries must equal # channels) use information about the channel cross section, which in our case is already taken care of by the channel weights ω, cf. Eq. (3). Furthermore, it is only the relative efficiencies among the channels that are important, and not their overall normalization (for the same reason). We therefore typically normalize the relative efficiencies such that the first element in the 11th row is equal to 1. As an example, the user may investigate one of the category measurements provided in the folder Expt_tables/latestresults-1.0.0/. From the 12th row onwards, the signal strength data is listed. Each row contains four values: the Higgs mass, the measured signal strength modifier at the lower edge of the 1σ uncertainty ("cyan") band,μ − μ, the central value (bestfit)μ, and finally the signal strength modifier at the upper edge of the 1 σ uncertainty band,μ + μ. In the case of a peak observable definition, as in Table 7, the data file ends after the 12th row, since the signal strength is only measured at a single Higgs mass value (corresponding to the signal). In contrast, for the construction of mass-centered observables, the data is listed here for the full investigated mass range, which is typically extracted from the correspondingμ plot using EasyNData [102].
As a further remark, we point out a general limitation in the implementation of experimental data: some results from the LHC experiments are given for the combination of data collected at different center-of-mass energies, e.g. at 7 and 8 TeV. These results cannot be disentangled by Higgs-Signals. Therefore, these observables are implemented as if the data was collected at the center-of-mass energy, which can be assumed to be dominating the experimental data. This approximation is valid, since both the observed and the predicted signal strengths are treated as SM normalized quantities. The only remaining inaccuracy lies in the SM channel weights, Eq. (3), which depend on the center-of-mass energy.
A complication arises in the assignment of Higgs observables if an analysis with one measured mass peak value is split up in several categories, each containing an individual signal rate measurement, see e.g. [84,103]. In this case, each category result defines a peak observable, however only one of these observables can be associated with the mass measurement from the analysis, which is going to contribute to the χ 2 . In all other categories this contribution has to be switched off. Nevertheless, this difference in the implementation can lead to inconsistent assignments of the Higgs boson(s) to the category observables. In order to enforce a consistent assignment, peak observables can build an assignment group. This enforces that the Higgs boson(s) are assigned to either all or none of the observables in this group, judged by the assignment status of the observable containing the mass measurement. For each peak observable, the assignment group can be specified in the experimental table, cf. Table 9. Note that the analysis IDs of the category peak observables have to be different from each other.

HiggsSignals applications
In this section we discuss a few example applications which demonstrate the performance of HiggsSignals. Most of the examples are chosen such that their results can be validated with official results from ATLAS and CMS. The quality of agreement of the reproduced HiggsSignals results with the official results justifies the Gaussian limit approximation in the statistical approach of HiggsSignals. Note that to a certain extent (which is difficult to estimate), the accuracy of the reproduced results suffers from the lack of publicly available information of the analysis efficiencies on the various production modes. At the end of this section, we briefly discuss a few HiggsSignals example applications, where the results incorporate all presently available Higgs data from the LHC and the Tevatron. Another example application of HiggsSignals within the context of the MSSM was presented in Ref. [16].

The peak-centered χ 2 method for a SM-like Higgs boson
As a first application we discuss the performance of the peakcentered χ 2 method on a SM-like Higgs boson. As already shown in Fig. 1b, a simple one parameter fit can be performed to the signal strength modifier μ, which scales the predicted signal rates of all investigated Higgs channels uniformly. In this fit the Higgs mass is held fixed at e.g. m H = 125.7 GeV. Using the signal strength measurements of the individual search channels obtained by the CMS collaboration [77], as given in Fig. 1b [117] and DØ [118], as summarized in Fig. 2. If possible, we implement results from the 7 and 8 TeV LHC runs as separate observables. However, if the only quoted result is a combination of both center-of-mass energies we implement it as an 8 TeV result. As mentioned in Sect. 3, we employ the quoted asymmetric uncertainties to account for the dominant effects of potentially remaining non-Gaussian behavior of the measurements. The H → γ γ and H → Z Z ( * ) → 4 analyses of ATLAS and CMS have a rather precise mass resolution, thus we treat the implemented mass value of their signal as a measurement which enters the Higgs mass part of the total χ 2 , cf. Sect. 3.1. Note however that the implemented mass value is not necessarily the most precise measurement of the Higgs mass but rather the mass value for which the signal strength was published by the experimental analysis. The Higgs mass can be determined more accurately from a simultaneous fit to the mass and the signal strength. This can be done with the mass-centered χ 2 method, as discussed in the next subsection. Note also that the Higgs mass values assumed in the signal strength measurements can differ by up to ∼2.5 GeV. It would be desirable if the experiments would present their best-fit signal strengths for all available channels (including specially tagged categories) also for a common Higgs mass (equal or close to the Higgs mass value preferred by the combined data) once a combination of different channels is performed. In the present case, global fits combining the signal strength measurements performed at different Higgs masses rely on the assumption that these measurements do not vary too much within these mass differences.
It can nevertheless be interesting to discuss the total χ 2 distribution obtained in the peak-centered χ 2 method as a function of the Higgs mass, m H . This serves as a demonstration of the three different Higgs mass uncertainty parametrizations (box, Gaussian, box+Gaussian pdfs), as well as the implications of taking into account the correlations among the systematic uncertainties in the χ 2 calculation. Furthermore, features of the automatic assignment of the Higgs boson to the peak observables can be studied. In the following example, we set the predicted signal strength for all Higgs channels to their SM values (μ i ≡ 1) and set the production and decay rate uncertainties to the values given in Eq. (9), as recommended by the LHC Higgs Cross Section Working Group for the SM Higgs boson around m H 125 GeV. We then evaluate the total peak-centered χ 2 for each Higgs boson mass m H ∈ [110, 140] GeV using the peak observables presented in Fig. 2. In the SM the Higgs mass is treated as a free parameter, which corresponds to setting the theory mass uncertainty to zero. In order to illustrate the effects of a non-zero theory mass uncertainty, we also consider a model with SM-like Higgs couplings, but which has a 2 GeV theory uncertainty on the predicted Higgs mass.
The total χ 2 mass distribution is shown in Fig. 3 for four different cases: In Fig. 3a,b the correlations among the systematic uncertainties of the signal rates, luminosity and Higgs mass predictions are neglected, whereas they are taken into account in Fig. 3c,d. In order to demonstrate the difference between the three parametrizations of the Higgs mass uncertainty we show the χ 2 distribution assuming a theoretical Higgs mass uncertainty of m = 0 GeV in Fig. 3a,c and m = 2 GeV in Fig. 3b,d, respectively. Furthermore, Fig. 3 includes the number of peak observables, which have been assigned with the Higgs boson, as a function of the Higgs mass. These are depicted by the faint graphs for each Higgs mass uncertainty parametrization.
The discontinuous shape of the χ 2 distribution is caused by changes in the Higgs boson assignment to the individual observables. Recall that, if the Higgs mass m H is too far away from the implemented mass position of the peak observable, the Higgs boson is not assigned to the signal. This yields a χ 2 contribution corresponding to no predicted signal, μ = 0, cf. Sect. 3.1. Most of the peak observables have different mass resolutions, therefore the χ 2 distribution has a staircase-like shape. At each step, the total number of peak observable assignments changes.
As can be seen in Fig. 3 all three parametrizations of the theoretical Higgs mass uncertainty yield the same total χ 2 values if the Higgs mass m H is far away from the implemented signal mass position, because typically observables which enter the Higgs mass part of the χ 2 in the Gaussian parametrization exhibit a decent mass resolution, and the  [117] and DØ [118], as they are implemented in Higgs-Signals-1.0.0 as peak observables. The left panel shows the Higgs mass value for which the signal strength was measured. A value with error bars indicates that the mass value is treated as a Higgs mass observable in the peak-centered χ 2 method, whereas a gray asterisk only serves as an indication of the Higgs mass value, which was assumed in the rate measurement. This value does not enter directly the total χ 2 .
For some LHC analyses, measurements for both the 7 and 8 TeV data exist, shown in blue and red, respectively. If the measurement is based on the combined 7/8 TeV dataset, we treat it as an 8 TeV measurement only. For the H → γ γ analyses from ATLAS and CMS, the special tagged categories were implemented as separate peak observables, including their efficiencies, but collected together in assignment groups. In total there are 4 Higgs mass observables and 63 Higgs signal rate observables. This data is used for the performance scans in Fig. 3 and the example applications in Sect. 5. 3 Higgs boson is only assigned if this χ 2 is low, i.e. m H ≈m. Conversely, at the χ 2 minimum at a Higgs mass m H ∼ 125-126 GeV, we obtain slightly different χ 2 values for the three parametrizations: Firstly, assuming that every observable is assigned with the Higgs boson, the minimal χ 2 is in general slightly higher in the Gaussian case than in the box and  A difference between the Gaussian and the theory box with experimental Gaussian (box+Gaussian) parametrization appears only for non-zero m. For m = 2 GeV the minimal χ 2 is obtained for a plateau m H ≈ (124.8-126.5) GeV in the box+Gaussian case, whereas in the Gaussian case we have a non-degenerate minimum at m H = 125.7 GeV. However, outside this plateau the χ 2 shape of the box+Gaussian increases faster than in the Gaussian case, since the uncertainty governing this Gaussian slope is smaller.
For the Gaussian parametrization of the theoretical Higgs mass uncertainty and no theoretical mass uncertainty the minimal χ 2 at m H = 125.7 GeV changes from 75.7 to 73.0 (for 63 signal strength observables and 4 mass observables) if we include the correlations among the systematic uncertainties in the χ 2 evaluation. In the case of a non-zero theoretical mass uncertainty, also the shape of the total χ 2 distribution can be affected when the correlations are taken into account. Recall that only in the Gaussian parametrization the correlations of the theoretical mass uncertainties enter the χ 2 evaluation, featuring a sign dependence on the relative position of the predicted Higgs mass value with respect to the two   [17] observed Higgs mass values, cf. Section 3.1.2. This results in a shallower slope of the χ 2 distribution at Higgs masses larger than all mass measurements, m H 126.8 GeV, since all mass observables are positively correlated in this case.
In conclusion we would like to emphasize that, although the direct χ 2 contribution from (the few) mass measurements to the total χ 2 might appear small in comparison to the χ 2 contribution from (many) signal strength measurements, the automatic assignment of Higgs boson(s) to the peak observables introduces a strong mass dependence, even for peak observables without an implemented mass measurement. Hereby, the procedure tries to ensure that a comparison of the predicted and observed signal strength is valid for each observable (depending on the mass resolution of the corresponding Higgs analysis), or otherwise considers the signal as not explainable by the model.

Combining search channels with the mass-centered
As a first demonstration of the mass-centered χ 2  searches. This is possible because the fullμ plot was published for these analyses for 7 and 8 TeV, except for the H → Z Z ( * ) → 4 search where only the combined 7/8 TeV result is available. 18 We scan the relevant Higgs mass range m H = (110 − 150) GeV, as well as the signal strength μ, and at each point (m H , μ) evaluate the mass-centered χ 2 using the correspondingμ plots as mass-centered observables. We then find the best-fit μ value (and the corresponding 1σ and 2σ regions) by minimizing the χ 2 (finding χ 2 = 1 and χ 2 = 4, respectively) for a fixed Higgs mass m H . This is shown in Fig. 4a and b for the H → γ γ channel and the combination of H → γ γ , H → W W ( * ) → ν ν and H → Z Z ( * ) → 4 , respectively. These results nicely agree with the corresponding official ATLAS results [17,105], which are shown in Fig. 4c,d for comparison. Especially at the signal around 126 GeV the Gaussian limit approxima- The obtained χ 2 values from this scan are summed and associated with the point (m H , μ). Thus we test the combined hypothesis of having a Higgs boson at m H with signal strength μ, and no signal elsewhere. The procedure is then repeated for all points in the two-dimensional (m H , μ) plane to obtain the 2D χ 2 likelihood map. The results are shown in Fig. 5 for each Higgs decay mode separately. For comparison, we also show the official ATLAS results [109, 110,122] as faintly colored contours. Qualitatively, the obtained 68 and 95 % C.L. regions (corresponding to χ 2 = 2.30 and χ 2 = 5.99, respectively) agree fairly well for H → Z Z and H → W W , whereas the H → γ γ result is shifted towards larger Higgs masses by around 0.8 GeV. A potential reason for this discrepancy is that effects of the mass scale systematic (MSS) uncertainty are only indirectly taken into account in HiggsSignals by simply using the corresponding plateau-shapedμ plot [122] instead of including the MSS uncertainty in the profile likelihood as a nuisance. Nevertheless, the 68 and 95 % C.L. regions still have a large overlap. Note also that the spiky structures of the contour ellipses in Fig. 5 are rather an artifact of our data extraction with EasyNData [102] than a physical effect. 20 A simultaneous fit to the ATLAS Higgs channels can also be performed. The best fit point of such a combination is found at where the uncertainties given refer to the 1D profiled 68 % confidence interval. We have verified that these results remain stable when varying the step sizes in the scan. The two discussed examples show the usefulness of the mass-centered χ 2 method. We focussed here on the validation of the method by comparing with official results from ATLAS. It is however easy to go beyond that and take all available data from ATLAS and CMS (and the Tevatron) into account for a simultaneous analysis. This we leave for a future study. However, we would like to emphasize again that the usefulness of this method strongly depends on the information (here in particular theμ plots for the individual channels) the experimental collaborations decide to publish.

Validation with official fit results for Higgs coupling scaling factors
A major task after the discovery of a Higgs-like state is the determination of its coupling properties and thus a thorough test of its compatibility with the SM. Both ATLAS [104,123] and CMS [77,115] have obtained results for Higgs coupling scaling factors in the framework of restricted benchmark models proposed by the LHC Higgs Cross Section Working Group [101]. Numerous other studies have been performed, both for Higgs coupling scaling factors  as well as for particular models, including composite Higgs scenarios [49][50][51][52], Two Higgs Doublet Models (2HDMs) [53][54][55][56][57], supersymmetric models [58][59][60][61][62][63][64][65][66] as well as other, more exotic extensions of the SM [67,68]. Here, we want to focus on the reproduction of the official ATLAS and CMS results using the Higgs coupling scaling factors as defined in the benchmark models of Ref. [101] in order to validate the Higgs-Signals implementation.
We validate with the ATLAS and CMS results, as presented at the Moriond 2013 conference [104,115]. The measurements from ATLAS and CMS, which are used as observ-   Tables 10  and 11 Table 11, have been performed for this assumed Higgs mass value, except for the H → γ γ categories being measured at m H = 125.0 GeV. Before we discuss the benchmark fits of Higgs coupling scale factors, we look at ATLAS and CMS fits that explicitly target the different production modes by combining channels with a particular decay mode. These fits allow to investigate sources of potential deviations between the official and the reproduced HiggsSignals results separately for each Higgs boson decay mode. Furthermore, unknown channel efficiencies can be adjusted within reasonable ranges, such that the agreement of the fit outcome is optimized. The signal composition of all included observables after this optimization is given in Tables 10 and 11. For the H → γ γ categories we use the published channel efficiencies. Two-parameter fits were performed for each decay mode to a signal strength modifier associated with the gluon fusion (ggf) and tt H production mechanisms, μ ggf+ttH , and a signal strength modifier for the VBF and V H production modes, μ qqH+VH . The results of the same fits performed with HiggsSignals are shown in Fig. 6 in direct comparison with the results from ATLAS [104,110] and CMS [77], which are faintly overlaid in the figure. Using the ATLAS results, Fig. 6a, the derived H → W W ellipse is in perfect agreement with the official result. Also the H → γ γ and H → Z Z ellipses agree reasonably well. The reproduced H → γ γ ellipse is slightly shifted towards larger values of μ ggf+ttH . A potential source of this discrepancy may be the different mass positions at which the measurements are performed. Moreover, the inclusion of correlations among the experimental systematic uncertainties becomes more important, the more the measurements are divided into smaller subsets/categories. These correlations are not publicly known and hence not taken into account by HiggsSignals. In the H → Z Z result, a significant difference between the approximations in HiggsSignals and the full profile likelihood (PLL) treatment can be observed. The PLL has a longer tail at large signal strengths, thus leading to extended 68 and 95 % C.L. regions at large values of μ qqH+VH . This is partly due to the Gaussian approximation, which is more constraining at large values than a Poisson distribution with the same central value, as is used in the PLL. This is especially relevant for the very small event count for VBF H → Z Z candidates. In addition, missing information about correlations of experimental systematics might contribute to the observed difference at large μ qqH+VH . Note also that one of the two H → Z Z category measurements that are publicly available [110], cf. Table 10, is a combination of the VBF and V H production channels, whereas the ATLAS analysis internally treats these channels as separate categories. The requirement of a positive probability density function (pdf) leads to the edge at negative μ qqH+VH in the official ATLAS result. We checked that adding the requirement of a positive signal strength modifier in HiggsSignalsthis edge is reproduced quite well.
Using the CMS results, Fig. 6b, we find reasonably good agreement between HiggsSignals and the official results for H → W W, bb, and τ τ . The H → γ γ ellipses roughly agree in the μ ggf+ttH range as well as in the correlations of the fit parameters (seen in the tilt of the ellipses). However, our reproduced ellipse is shifted towards lower values of μ qqH+VH . In order to investigate the influence of correlated experimental systematic uncertainties, we introduced a tunable degree of correlation among the VBF-tagged H → γ γ categories. A much better agreement between Higgs-Signals and the official result is obtained when around 30 % of the measured relative signal strength uncertainty of the VBF-tagged categories is treated as a fully correlated uncertainty. This indicates that including this type of (not public) information could potentially lead to an improvement of the HiggsSignals methodology in certain channels. A similar effect from correlations of experimental systematics may lead to the differences observed in the H → τ τ ellipses. The H → Z Z ellipse can only be roughly reproduced using the publicly available data for the two H → Z Z observables. Even after adjusting their production mode efficiencies, cf. Table 11, differences remain due to the Gaussian approximation and possibly further (publicly unavailable) information on the VBF-likeness of the observed signal events [82].
Using the results in Fig. 6, we can estimate the typical differences between the official results from ATLAS and CMS and the HiggsSignals implementation. We classify the difference in two ways: first, the χ 2 in our fit between the official best fit point from the collaboration and the best fit point from HiggsSignals, and second, the distance between the two best fit points in the parameter space relative to the 1σ uncertainty in the direction spanned by these two best-fit points.
For the comparison with the official ATLAS result, cf. Fig. 6a, the χ 2 is 0.158, 3.5 × 10 −4 and 3.6 × 10 −3 for H → γ γ , H → W W and H → Z Z, respectively. For H → γ γ the difference is small but non-negligible, as pointed out before. The latter two can be regarded as insignificant. The difference between the best fit points of ATLAS and HiggsSignals, relative to the corresponding 1σ uncertainty is 24, 6.6 and 7.7 %, respectively. Also here, a reasonable agreement well within 1σ is observed.
For the comparison with the official CMS result, cf. Fig. 6b, the differences in χ 2 between the best fit points are 0.51 for H → γ γ , 0.34 for H → Z Z, and less than 0.05 for the other channels. Plausible reasons for the differences in H → γ γ and H → Z Z are discussed above. For the remainder of channels there is very good agreement. The same picture arises for the relative distance of best fit points in parameter space with respect to the 1σ uncertainty measured in the same direction, where the largest deviation is observed for H → γ γ with 44 %. Still, this is well within 1σ and should be sufficient for exploratory studies of new physics models. All other channels agree significantly better.
We now turn to the discussion of global fits in the Higgs coupling scale factor benchmark scenarios. Regarding the interpretation of the following benchmark fits, it should be kept in mind that only two parameters are allowed to deviate from their SM values, while all other Higgs couplings and partial decay widths have been fixed to their SM values. The way an observed deviation from the SM manifests itself in the parameter space of coupling strength modifiers κ i will sensitively depend on how general the basis of the κ i is that one has chosen. Furthermore the framework of the coupling strength modifiers κ i as defined in Ref. [101] is designed for the analysis of relatively small deviations from the SM. In case a firm preference should be established in a parameter region that is very different from the SM case (e.g. a different relative sign of Higgs couplings), the framework of the coupling strength modifiers κ i would have to be replaced by a more general parametrization.
The first benchmark model we want to investigate is a two-dimensional fit to universal scale factors for the Higgs coupling to the massive SM vector bosons, κ V , and to SM fermions, κ F . In this fit it is assumed that no other modifications to the total width than those induced by the coupling scale factors κ F and κ V are present, allowing for a fit to the coupling strength modifiers individually rather than to ratios of the scale factors [101]. Note that the loop-induced effective H γ γ coupling is derived in this approximation from the (scaled) tree-level couplings Htt and H W + W − and thus exhibits a non-trivial scaling behavior. In particular the interference between the t and W boson loops introduces a dependence on the relative sign of the scale factors κ F and κ V . In the case of a relative minus sign this interference term gives a positive contribution to the H γ γ coupling.
The reconstructed ATLAS and CMS fits obtained with HiggsSignals are shown in Figs. 7a and 8a, The (2D) compatibility with the SM hypothesis of these points is 11.1 and 28.4 % for the reproduced ATLAS and CMS fit, respectively. In order to probe the presence of BSM physics in the Higgs boson phenomenology a fit to the loop-induced Higgs couplings to gluons, κ g , and photons, κ γ , can be performed. In this fit it is assumed that all other (tree-level) Higgs couplings are as in the SM and no new Higgs boson decay modes exist. Figures 9a and 10a show the 2D likelihood maps in the (κ γ , κ g ) parameter plane for the HiggsSignals result using the ATLAS and CMS observables, respectively. The corresponding official ATLAS and CMS results are given in Figs. 9b and 10b. Again, we observe reasonably good agreement with the official results. We find the best fit points at These are (2D) compatible with the SM at the level of 7.6 and 17.1 %, respectively. In the ATLAS fit, the best-fit region obtained by HiggsSignals is slightly shifted with respect to the official result towards lower values of κ g by roughly κ g ∼ 0.05-0.10, whereas the agreement in κ γ direction is very good. In the CMS fit, the agreement is better. Here, to photons, κ γ , and gluons, κ g . In these fits, the Higgs boson mass is assumed to be 126 GeV. The full available data from the Tevatron and LHC experiments as presented at the Moriond 2013 conference (and shortly after) is used. This data is summarized in Fig. 2 the HiggsSignals χ 2 distribution is slightly shallower than the official CMS likelihood at low values of κ γ , leading to slightly larger C.L. contours.
We conclude this section by pointing out that, despite some discrepancies that are observed in fits to single decay modes using subsets of the available measurements, Fig. 6, the combination of all available channels from each experiment reproduces the official results quite well. We are therefore confident that the accuracy of the HiggsSignals method is sufficient for surveys of new physics parameter spaces compatible with the Higgs measurements, and for simple coupling scale factor fits. For a more precise determination of the Higgs boson coupling structure with HiggsSignals, however, it would be desirable if the experimental collabo-rations made information on efficiencies, correlated experimental uncertainties and all category measurements publicly available in a more complete way. We would expect a significant reduction of the observed remaining discrepancies if this information was included in HiggsSignals.

Example applications of HiggsSignals
We now go beyond validation and repeat the two discussed Higgs coupling scaling factor fits including the full presently available data from the LHC and Tevatron experiments, as listed in Fig. 2. This includes data presented up until shortly after the Moriond 2013 conference. We assume a Higgs boson mass of 126 GeV. The fit results for the Higgs coupling scale factors (κ V , κ F ), defined in Sect. 5.2 and [101], are shown in Fig. 11a. The best-fit point is found at κ V = 0.99 +0.06 −0.06 , κ F = 0.86 +0.14 −0.10 , with χ 2 /ndf = 68.7/61, where the profiled one-dimensional 68 % C.L. uncertainties are given. For this fit the SM point is found to be located well within the 68 % C.L. contour, with a (2D) χ 2 compatibility with the best fit point of 59 A similar improvement is seen for the (κ γ , κ g ) fit, shown in Fig. 11b, where the best fit point is found at κ γ = 1.12 +0.10 −0.08 , κ g = 0.90 +0.09 −0.08 , with χ 2 /ndf = 67.6/61, which can be compared with Figs. 9 and 10. Here, the SM is compatible with the fit result at the level of 31.8 %. The fit shows a weak tendency towards slightly reduced κ g and slightly enhanced κ γ . The discrimination power on κ g will increase only slowly with more data, since the large uncertainty of the rate prediction for single Higgs production is already the dominant limitation of the precision of the combined fit [45]. As a further example application we performed fits in three of the MSSM benchmark scenarios recently proposed for the interpretation of the SUSY Higgs search results at the LHC [126]. These scenarios are defined in terms of two free parameters, tan β = v 2 /v 1 (the ratio of the vacuum expectation values of the two Higgs doublets), and either M A (the CP-odd Higgs boson mass) or μ (the Higgsino mass parameter). The other parameters are fixed to their default values as specified in [126] to exhibit certain features of the MSSM Higgs phenomenology. For each parameter point in these two-dimensional planes we calculated the model predictions with FeynHiggs-2.9.4 and evaluated the total χ 2 , comprised of the LEP Higgs exclusion χ 2 value [7,11] obtained from HiggsBounds-4 [15,16], as well as the total χ 2 from HiggsSignals using the peak-centered χ 2 method. The theoretical mass uncertainty of the lightest Higgs boson is set to 2 GeV when treated as a Gaussian uncertainty (i.e. in the LEP exclusion χ 2 from HiggsBounds and in Higgs-Signals), and to 3 GeV in the evaluation of 95 % C.L. LHC exclusions with HiggsBounds.
The first scenario is an updated version of the well-known m max h benchmark scenario [126,129], where the masses of the gluino and the squarks of the first and second generation were set to higher values in view of the latest bounds from SUSY searches at the LHC, see [126] for details. The  Fig. 12 in the (M A , tan β) plane. Besides the colors indicating the χ 2 = χ 2 − χ 2 best−fit distribution relative to the best-fit point (shown as a green star) we also show the parameter regions that are excluded at 95 % C.L. by LHC searches for a light charged Higgs boson (darkgreen, coarsely striped) [128], neutral Higgs boson(s) in the τ τ final state (orange, checkered) [127] and the combination of SM search channels (red, striped) [115], as obtained using HiggsBounds. As an indication for the parameter regions that are 95 % C.L. excluded by neutral Higgs searches at LEP [7,11] we include a corresponding contour (black, dashed) for the value χ 2 LEP,HB = 4.0. Conversely, the parameter regions favored by the fit are shown as 68 and 95 % C.L. regions (based on the 2D χ 2 probability w.r.t. the best fit point) by the solid and dashed gray lines, respectively.
As can be seen in the figure, the best fit regions are obtained in a strip at relatively small values of tan β ≈ 4.5-7, where in this scenario M h ∼ 125.5 GeV is found. At larger tan β values the light Higgs mass in this benchmark scenario (which was designed to maximise M h for a given tan β in the region of large M A ) turns out to be higher than the measured mass of the observed signal, resulting in a corresponding χ 2 penalty. At very low tan β values the light Higgs mass is found to be below the preferred mass region, again resulting in a χ 2 penalty. Here, the χ 2 steeply rises (for M h 122 GeV), because the mass-sensitive observables (H → γ γ, Z Z ( * ) ) cannot be explained by the light Higgs boson anymore, cf. Sect. 5    mainly SM-like couplings. Consequently, the χ 2 contribution from the rate measurements is similar to the one for a SM Higgs boson. In this regime, the Higgs mass dependence of the total χ 2 (from HiggsSignals) is comparable to the results shown in Fig. 3d. We find the best fit point at (M A , tan β ) = (674 GeV, 5.0) with χ 2 /ndf = 70.2/66. The number of degrees of freedom (ndf) comprises 63 signal strengths and 4 mass measurements presented in Fig. 2, as well as one LEP exclusion observable from HiggsBounds.
The second scenario that we discuss here is a modification of the m max h scenario with a lower value of X t , leading to M h ∼ 125.5 GeV over nearly the whole (M A , tan β) plane [126]. This so-called m mod+ h scenario is shown in Fig. 13 (with the same colors and meaning of thecontours as for the m max h scenario, Fig. 12). The best fit point is found at (M A , tan β) = (674 GeV, 9.3) with χ 2 = 70.7/66. Only slightly larger χ 2 values are found over the rest of the plane, except for the lowest M A and tan β values, where M h is found to be below the preferred mass region. As in the preferred region for the m max h scenario the lightest Higgs boson is mostly SM-like here, and the χ 2 from the rates is close to the one found in the m max h scenario. As a final example, we performed a fit in the low−M H benchmark scenario of the MSSM [126]. This scenario is based on the assumption that the Higgs observed at ∼ 125.5 GeV is the heavy CP-even Higgs boson of the MSSM. In this case the light CP-even Higgs has a mass below the LEP limit for a SM Higgs boson of 114.4 GeV [7], but is effectively decoupled from the SM gauge bosons. The other states of the Higgs spectrum are also rather light, with masses around ∼130 GeV, so that this scenario offers good prospects for the searches for additional Higgs bosons [19,65,66]. Since M A must be relatively small in this case the (μ, tan β) plane is scanned [126], where only tan β 10 is considered. The CP-odd Higgs boson mass is fixed to M A = 110 GeV.
Our results are shown in Fig. 14. The 95 % C.L. excluded regions are obtained from the same Higgs searches as in Fig. 12, except for the red patterned region, which results from applying the limit from the CMS SM Higgs search H → Z Z ( * ) → 4 [82] to the SM-like, heavy CP-even Higgs boson (see below). Two distinct best-fit regions are found [126]: The parameter space with μ ∼ (1.6-2.0) TeV and tan β ∼ 4-6 predicts a heavy CP-even Higgs boson with a well compatible mass value M H ≈ 126 GeV and SM-like couplings. However, large parts (at low tan β 4.9) of this region favored by the rate and mass measurements are severely constrained by charged Higgs searches [128]. The best-fit point is found at the edge of the excluded region at (μ, tan β) = (1850 GeV, 4.9). The second region favored by the fit is located at large values of μ ∼ (2.4-2.9) TeV and tan β ∼ 6-7. Here, the masses of the CP-even Higgs bosons are generally lower. For instance, at (μ, tan β) ∼ (3070 GeV, 6.0), we have M h ≈ 76.1 GeV and M H ≈ 122.8 GeV. For slightly larger (lower) values of μ (tan β) we find a steep edge in the HiggsSignals χ 2 distribution, because M H becomes too low to allow for an assignment of the heavy CP-even Higgs boson to all mass-sensitive peak observables, cf. the results shown in Fig. 3d, Sect. 5.1.1. Due to the low mass of the light CP-even Higgs boson in this region, the LEP channel e + e − → h A [11] is kinematically accessible and contributes a non-negligible χ 2 which increases with μ. The parameter space between the two preferred regions suffers a rather large χ 2 penalty, since in particular the predicted rates for the H → Z Z ( * ) , W W * channels are above the rates measured at the LHC, as can also be seen from the 95 % C.L. exclusion by HiggsBounds in this region.
At the best-fit point we find a χ 2 /ndf = 80.3/66. Compared with the light CP-even Higgs interpretation of the observed signal, as discussed in the m max h and m mod+ h scenarios, the fit quality is only slightly worse.

Conclusions
We have presented HiggsSignals, a public Fortran code to test the predictions of models with arbitrary Higgs sectors against measurements obtained from Higgs searches at the LHC, the Tevatron, and any potential future experiment. The code is publicly available at http://higgsbounds.hepforge.org The code features two statistical tests, one which determines the compatibility of the model with experimentally observed Higgs signals, and a second which tests for general compatibility with the observed Higgs data at the predicted mass(es) of the Higgs boson(s) in the theory. Since the two tests are complementary, we also provide a method to perform both simultaneously and use the combined results for models with multiple Higgs bosons.
The main experimental results used by HiggsSignals are the signal strength modifiers,μ, as a function of the Higgs mass in the various search channels. These results have to be supplemented by their respective experimental uncertainties, μ, and (preferably, if this information is available) with the experimental efficiencies and correlations. The information onμ and μ channel by channel constitutes the most general and robust experimental input for testing the theoretical predictions of different models, and we strongly encourage the experimental collaborations to continue to make them public with as much details provided as possible.
The default implementation of HiggsSignals uses theμ results available from the LHC and the Tevatron, and it is planned to continuously update these results in forthcoming versions of HiggsSignals. However, it is easily possible for the user to include additional experimental data. For assessing possible future projections it is also possible to implement hypothetical future experimental results.
The input that has to be provided by the user (and which is similar to the HiggsBounds input) consists of the Higgs boson masses, preferably the corresponding theory uncertainties, the Higgs production cross sections and decay branching ratios, where several levels of approximation are possible. In case of the MSSM also the SLHA [93,94] can be used as input/output format.
We presented in detail the two statistical methods provided by HiggsSignals: the peak-centered χ 2 method, in which each observable is defined by a Higgs signal rate measured at a specific hypothetical Higgs mass, corresponding to a tentative Higgs signal. In the second, the mass-centered χ 2 method, the χ 2 is evaluated by comparing the signal rate measurement to the theory prediction at the Higgs mass predicted by the model. It was described how these two methods can be combined, as it is an option of Higgs-Signals, to yield the most reliable consistency test. In this combination, the mass-centered χ 2 method is applied only to those Higgs bosons which have not yet been tested with the peak-centered χ 2 method against the same data. Similarly, in order to include a more complete set of constraints on the Higgs sector, it is recommended to use Higgs-Signals together with HiggsBounds to test the model under consideration also against the existing Higgs exclusion bounds.
The installation, usage and subroutines of Higgs-Signals were explained in detail, together with the various input and output formats. It was explained how the user can add new (hypothetical) experimental data. Several pre-defined example codes were presented that permit the user to get familiar with HiggsSignals and, by modifying the example codes, analyze own models of interest. As an example, by linking HiggsSignals to FeynHiggs, the consistency of any MSSM parameter point with the observed LHC signal can be analyzed in a simple way. Furthermore, some example codes demonstrate how to use HiggsBounds and HiggsSignals simultaneously in an efficient way.
We have presented several examples of the use of Higgs-Signals. As a first example the combined best-fit signal strength has been determined. For the peak-centered χ 2 method the mass dependence, the impact of correlations between the systematic uncertainties and the treatment of theoretical uncertainties has been discussed in detail. For the case of a SM-like Higgs boson, we demonstrated how the mass can be determined from a fit to the signal rate measurements as a function of the mass using the mass-centered χ 2 method. Moreover, we employed this method for a combination of different search channels over the full investigated mass range. Various fits for coupling strength modifiers have been carried out using the peak-centered χ 2 method. Their results have been compared for validation purposes with official results from the ATLAS and CMS collaborations, and very good agreement has been found.
It is expected that the agreement with the official results published by ATLAS and CMS could be improved even further if relative signal efficiencies of different production modes in all search channels would be publicly provided by the experimental collaborations. The same applies to a more complete description of the impact of individual experimental systematic uncertainties and their correlations amongst search channels. In particular, it would be useful if systematic uncertainties were given as a relative error on the quoted signal strength. We would furthermore welcome the publi-cation of the fullμ plot for every analysis to allow a χ 2 test at various Higgs masses.
Going beyond just a validation of HiggsSignals results, we have also given a few examples of Higgs-Signals applications. In particular, we have performed fits of Higgs coupling scaling factors including the full presently available data from both the LHC and the Tevatron. In the first step of the mass-centered χ 2 method, Higgs-Signals constructs possible Higgs boson combinations following the Stockholm clustering scheme. In our example, h 2 and h 3 are combined in a Higgs cluster, denoted by h 23 , for the H → W W analysis since their mass difference is lower than the experimental mass resolution. In all other cases, the Higgs bosons are tested singly, thus we have in total 11 observables. The mass and its uncertainty associated with the Higgs cluster h 23 are derived from Eqs. (22) and (23) to m 23 = 137.5 GeV and m 23 = 1.4 GeV. Its predicted signal strength is μ 23 = 0.7.
In the second step, the observed quantitiesμ α and μ α have to be determined from theμ plots for each observable α. In order to take into account the theoretical mass uncertainties, the relevant mass region is scanned to construct the tentative total χ 2 i (m ) distribution for each Higgs boson h i , as described in Sect. 3.2. For this example, the χ 2 i (m ) distributions for the box-shaped and Gaussian parametrization of the theoretical mass uncertainty are shown in Fig. 15a, b, respectively. At the mass positionm i , where χ 2 i (m ) is minimal, the observed quantitiesμ α and μ α are extracted from theμ plots for those observables α, which test the Higgs boson i.
In the box-shaped parametrization, the measured signal strengths of all mass-centered observables which test h 1 are defined atm 1 = 124.7 GeV, where χ 2 1 is minimal. In contrast, the Higgs bosons h 2 and h 3 form the Higgs cluster h 23 in the H → W W analysis, therefore their allowed mass variations are restricted to the overlap regions M 2 ∩ M 23 and M 3 ∩ M 23 , cf. Eq. (17), respectively. In those observables, where h 2 (h 3 ) is tested singly, the measured quantities are defined atm = 136.1 GeV (138.9 GeV). For the observable testing the Higgs cluster h 23 the observable is defined by the minimum of the joint χ 2 distribution, which is located at m = 138.9 GeV.
In the Gaussian parametrization the mass variation is less restricted. In contrast to the box-shaped parametrization, each mass variation is allowed over the full available mass range of the analyses, however, the additional contribution of the Higgs mass to the tentative χ 2 , cf. Eq. (19), tries to keep the varied mass close the its original predicted value. From the minimum of each tentative χ 2 distribution, the observed quantities of analyses, which test either h 1 , h 2 or h 3 singly, are defined atm = 124.8, 133. We want to illustrate how the experimental data changes, if we choose to fold the theoretical Higgs mass uncertainty, m, into the originalμ plot, as discussed in Sect. 3.2. For this, we look at theμ plot published by ATLAS for the H → Z Z ( * ) → 4 search [121] and convolve it with a uniform (box) or Gaussian Higgs mass pdf, centered at m H , for various theoretical mass uncertainties m = (0, 2, 5) GeV, following Eqs. (20) and (21). This is done over the full mass range, m H ∈ [112, 160] GeV, to obtain the results shown in Fig. 16. For m = 0 GeV, theμ plot is unchanged, whereas for increasing m it becomes smoother and fluctuations tend to vanish. This happens faster for the Gaussian pdf.