A method and tool for combining differential or inclusive measurements obtained with simultaneously constrained uncertainties

A method is discussed that allows combining sets of differential or inclusive measurements. It is assumed that at least one measurement was obtained with simultaneously fitting a set of nuisance parameters, representing sources of systematic uncertainties. As a result of beneficial constraints from the data all such fitted parameters are correlated among each other. The best approach for a combination of these measurements would be the maximization of a combined likelihood, for which the full fit model of each measurement and the original data are required. However, only in rare cases this information is publicly available. In absence of this information most commonly used combination methods are not able to account for these correlations between uncertainties, which can lead to severe biases as shown in this article. The method discussed here provides a solution for this problem. It relies on the public result and its covariance or Hessian, only, and is validated against the combined-likelihood approach. A dedicated software package implementing this method is also presented. It provides a text-based user interface alongside a C++ interface. The latter also interfaces to ROOT classes for simple combination of binned measurements such as differential cross sections.


Introduction
A common technique to reduce the impact of systematic uncertainties, in particular in precision measurements in high energy physics, is to constrain the range of their variations by the data. A simultaneous fit of these variations and the parameters to be measured can be performed based on prior knowledge of the uncertainties and suitable distributions to constrain them. This technique can reduce the total uncertainty in many cases significantly, as in Refs. [1][2][3], and can be used to measure several parameters simultaneously (see e.g. Refs. [4,5]). However, it leads to non-negligible correlations between all fitted parameters. These are particularly important when at least one measurement that uses this technique contributes to a combination.
The most consistent approach to such combination would be to define a combined likelihood based on the original models, including all systematic variations, and the original data the models were fit to. Also other commonly used combination techniques and the corresponding software tools would require this information [6][7][8][9], which is publicly available only in very rare cases and often unrecoverable. This poses a serious problem for a consistent combination involving measurements obtained with simultaneous nuisance parameter fits.
The method described here provides a solution for this problem, since it is based on the central results and their covariance or Hessians, only. It allows separating constraints and correlations imposed by the previously fitted data from those that stem from prior knowledge of the systematic variations.
Therefore, the combination can be performed accounting for correlations between the measurements as well as for correlations and constraints within each individual measurement.
The dedicated software tool "Convino" is also presented in this article. It is specifically developed to perform combinations based on the method described here and provides a simple text-based user interface that can be used without knowledge of any programming language. Assumptions on correlations can be varied in an automated way. Moreover, partially correlated measurements of different quantities (e.g. bins of a differential distribution) can be combined simultaneously accounting for all correlations. In addition to the text-based interface, a C++ interface is provided to define the input to the combination. This interface can either read basic C++ standard library data types or ROOT [10] histogram and graph classes, which are commonly used in high-energyphysics analyses.
The combination method is described in Sect. 2. It is validated against the combined likelihood approach in Sect. 3. The effect of neglecting correlations within the same measurement is studied in Sect. 4. The installation and the user interface of the Convino program is described in Sect. 5.

Combination method
The combination is performed using a χ 2 minimisation. The χ 2 is defined as It is composed of three terms: the term χ 2 s,α represents the results of each measurement α and its statistical uncertainties. It follows a Neyman or Pearson χ 2 definition, with the statistical uncertainty being fixed for each measurement or being scaled with the combined value, respectively. A measurement can aim to determine a set of quantities, e.g. bins of a differential cross section, where the Pearson definition is more applicable. Alternatively a single quantity can be measured, e.g. the mass of a particle from a fit of an invariant mass peak position, where the Neyman definition is presumably better suited to describe the measurement. In both cases, the measured quantities are referred to as estimates in the following. The additional term χ 2 u,α describes the correlations between the systematic uncertainties and constraints on them from the data for each measurement α. The last term, χ 2 p , incorporates prior knowledge of the systematic uncertainties and correlation assumptions between uncertainties of the measurements to be combined.
In a typical measurement that exploits simultaneous constraints on the uncertainties from the data, the sources of uncertainties are uncorrelated prior to the fit to the data and the knowledge about their variations is modelled by independent penalty terms in the original likelihood. The goal of the method described here is to find an approximation for this likelihood that allows disentangling the independent penalty terms from the constraints and correlations between them, which are typically introduced by the data. 1 Therefore, the first central assumption is that the original likelihood for a measurement α can be approximated by: where (P α i (λ i )) 2 represents the penalty term for a systematic uncertainty i parametrised by a continuous parameter λ i , such that λ i = 1 corresponds to a 1σ variation. The terms χ 2 s,α and χ 2 u,α are defined as: Here, x α μ is the estimate μ obtained in measurement α andx μ the combined value to be determined. The relation between both is given by τ α μ = (X μ /x α μ ) 0.5 for the Pearson χ 2 definition. In case of the Neyman χ 2 , all τ α μ = 1. The matrix M represents the inverted statistical covariance of the estimates. The parameters k α μi and K α μi model the effect of the systematic variations on the estimates for absolute or relative uncertainties, respectively. Here, a relative uncertainty is an uncertainty that scales with the measured value such as e.g. the luminosity uncertainty in a cross-section measurement, while absolute uncertainties have a constant value irrespective of the central result. In principle, also other classes of dependencies can be incorporated through additional terms in ξ α μ . The correlations between the uncertainties and the constraints that stem from the fit to the data are described by the matrix D α . Throughout this paper, indices μ and ν are used for estimates, while systematic uncertainties are denoted with indices i and j.
The procedure to obtain the parameters of χ 2 α from the measurements to be combined are discussed in the following -firstly for results obtained through a simultaneous nuisance parameter fit and secondly for the specific case of orthogonal uncertainties.

Measurements obtained with simultaneous fits
The second central assumption of the method is that the parameters of χ 2 α can be determined from the Hessian of measurement α evaluated at the best-fit values,H α in . The entries of the Hessian can be ordered such that the matrix can be split in the following sub-matrices: whereM describes the relation between the estimates x α μ and x α ν . The relation between systematic variations and the estimates is described by κ. The matrixD quantifies the relation between the systematic variations.
All parameters of χ 2 α are determined by calculating analytically the Hessian of χ 2 α , H α (0), and identifying the resulting terms with their counterparts of the inputH α in . Here 0 means The components are calculated as follows: withk α νi = K α νi +k α νi . The matrix M can be directly identified withM. SinceM stems from a measurement of a physics quantity,M is positive definite and therefore invertible. Thus, the parametersk α νi can be determined as: Since a variation i is either relative or absolute,k α νi equals either K α νi or k α νi , with the other parameter being 0. The terms describing the analytic relations between the systematic uncertainties are calculated as: Only Gaussian penalty terms describing the prior knowledge of the uncertainties are considered 2 (P α i (λ i ) = λ i ), which simplify the equations, since in this case 2 In principle also other penalty terms (e.g. log-normal priors) can be accounted for. However, these lead to a non-constant D α .
In consequence D α becomes: such that all parameters of Eq. 2 are defined.

Measurements with orthogonal uncertainties
For orthogonal uncertainties, the same initial χ 2 described in Eq. 2 is used. However, the calculation of its parameters does not necessarily require the full Hessian. Instead, the calculation of the parameters simplifies to: with σ α μi being the contribution of uncertainty i to the total uncertainty σ α μ total of estimate x α μ . The matrix D α is 0, the terms of M are calculated as: Here, ρ μν is the statistical correlation between estimate μ and ν, and σ μ and σ ν are the corresponding statistical uncertainties.
For the orthogonal uncertainties as well as for measurements obtained by simultaneous fits, the constraints from the prior knowledge of the uncertainties are implemented in Eq. 1 through the term: with C being the matrix describing the correlation assumptions between the systematic uncertainties. In case no correlations are assumed, the term simplifies to: Only Gaussian penalty terms are considered in the following, such that P i (λ i ) = λ i . For a combination, C will be of the structure with matrices A describing the correlation assumptions, and 1 being the identity matrix.

Technical implementation
The final minimisation of Eq. 1 is performed using the Minuit algorithms [11]. The total uncertainty on each combined value is determined by scanning χ 2 = χ 2 min + 1 using the Minos algorithm. These algorithms as implemented in ROOT 6 as "TMinuit2" are employed.
The correlations that are assumed between systematic uncertainties can vary between −1 and 1. These extremes are special cases for which the correlation matrix C becomes non-invertible. In practice, a correlation of C i j = ±1 means that parameters i and j describe the same variation. In such cases, an entry C i j = ±1 is replaced by C i j = ±(1 − 10 −3 ). The difference to ±1 is almost negligible. For illustration, 2 × 2 parameters are chosen with C i j ≈ ±1. The affected part of the χ 2 , χ 2 F , can be simplified to Given that a variation of λ = ±1 corresponds to only a fraction of the total uncertainty on each estimate, the effect of the approximation C i j = ±(1 − 10 −3 ) is negligible.

Validation
The validation is based on pseudo-measurements. Each pseudo-measurement is a binned likelihood fit with steerable central value and bin-wise uncertainties. This has the advantage, that the full likelihood of each pseudo-measurement is known and can be adjusted to different scenarios. Therefore, it is possible to compare the results obtained with the method proposed here to the ones obtained using the combined likelihood as reference. Since the latter in principle contains arbitrarily more parameters, small deviations are expected.
The validation is first performed with respect to the statistical bias, only. Secondly, the modelling of systematic uncertainties is tested with respect to correlations between the uncertainties of the pseudo-measurements, and the modelling of relative uncertainties.
For each pseudo-measurement a Poisson likelihood is chosen to determine the central resultx μ . The bin-wise uncertainties are randomly generated and modelled by the parameters λ i . In the case that an uncertainty corresponds to an absolute variation, its effect on each bin is generated independently. For more than one bin (N bins > 1) this results in correlations between the uncertainties after the fit, as well as in constraints on their variations.
The likelihood for pseudo-measurement α is defined as: with P being the Poisson likelihood andP α i (λ i ) the Gaussian penalty terms modelling the prior knowledge of each uncertainty. The parameters X α μ andX α μi are given as: where K α ν j describes the magnitude of global relative variations and k α νi j absolute shape variations, different for each bin i. The value of x α μ is the input to each pseudo-measurement and corresponds to the number of events that would be observed in a real measurement. The elements of the matrices K α and k α i are chosen to describe different validation scenarios. Finally, the fit to determinex α μ is performed and the resulting Hessian is recorded.
The combined likelihood for several pseudo-measurements is given by: where φ models the prior knowledge of N systematic uncertainties and the correlation assumptions between them, analogue to χ 2 p in Eq. 16. For every validation step, the difference Δx between the result obtained with the method proposed in this document and using the combined likelihood is recorded. The compatibility of both approaches is is quantified by Δx/σx , which is the difference between their central results normalised to the total uncertainty of the combinedlikelihood combination.

Statistical bias
To evaluate the statistical bias, the impact of systematic uncertainties on each pseudo-measurement is set to 0, corresponding to K = 0 and k = 0. Only one quantity,x, is determined from two estimates x a and x b , chosen as: with s being a scaling factor and γ describing the compatibility between the estimates. The latter is chosen to be either γ = 10 to describe two very incompatible measurements or γ = 3 for a more realistic scenario where both estimates still agree well enough to enter a combination. Two pseudo-measurements are generated for each choice of s and γ and are . 1 Difference between the combined values using a direct Poissonlikelihood combination and the method proposed here with Neyman and Pearson χ 2 definition relative to the total uncertainty. The estimates to be combined differ by about γ σ and are displayed as a function of the first estimate's relative statistical uncertainty. The second estimate's statistical uncertainty scales accordingly combined either using a Pearson or Neyman χ 2 definition. For both choices, the uncertainties on the combined results agree very well with the ones obtained using the direct combination based on L comb . The bias of the central value is shown in Fig. 1 relative to the uncertainty of the combined value. It behaves as expected: it is smaller but of opposite sign for the Pearson χ 2 definition and is reduced with smaller statistical uncertainties and better compatibility between the results.

Systematic uncertainties
The effect of absolute systematic uncertainties is evaluated by combining two pseudo-measurements, with randomly chosen elements of the matrices k α ν . An upper threshold t is defined, such that for each element i, j: limiting the contribution of systematic uncertainties. Two bins, two systematic uncertainties, and one x α μ per pseudo-measurement are considered. The sign of k α νi j is chosen to be constant for each systematic uncertainty. The estimates x a and x b for measurement a and b are set to: x a = 30000 and (28) to reduce the effect of statistical uncertainties. The resulting statistical uncertainty of 0.6% does not account for the difference of 2% between both values, such that the modelling of the systematic uncertainties will affect the combination significantly.
For large systematic variations the maximisation of Eq. 21 with Minuit can become numerically unstable. This is the case when the variation becomes as large as the nominal entry, X α μ , in at least one of the bins. Therefore, the Poisson likelihood is approximated with a Gaussian form, which is valid for low statistical uncertainties such as in this test. Thus, L α becomes: The matrix S α allows modelling direct statistical correlations betweenX α μi andX α νi . Here, S is set to 1. In total 2 × 20,000 pseudo-measurements are generated, each with a different random choice of the uncertainties. The total relative uncertainty, σ x /x, on the estimate of pseudo-measurement a is shown in Fig. 2 for different values of the threshold t. Depending on t, the uncertainty varies from moderate values to more than 100%. The same applies to pseudomeasurement b (not displayed). The average constraints on the systematic uncertainties reach from about 90% (t = 0.01) to 50% (t = 1.0) with respect to their initial 1σ variation.
In a first validation step, each uncertainty of one pseudomeasurement is assumed to be highly correlated with exactly one uncertainty of the other pseudo-measurement by assigning a correlation factor c = 0.99. A total of 20,000 combinations are performed. The ratio Δ r σ between the uncertainty on the combined value obtained with the method proposed here and by maximising L comb is shown in Fig. 3 as a function of t. Asymmetric uncertainties on the combined value are accounted for and are equally well described. With an increasing contribution of the systematic uncertainties, the Δ r σ -distribution becomes slightly broader, but does not indicate any numerically relevant mismodelling. The resulting values for Δx/σx are illustrated in Fig. 4. For t = 0.01, the statistical uncertainties are non negligible. This leads to a small bias towards lower values introduced by the choice of Eq. 30 and discussed in the Sect. 3.1. However, for all choices of t and all pseudo experiments, the differences between the direct likelihood approach and the method described here are well below 5% of the total uncertainty and can therefore be considered negligible. Moreover, the dependence on the assumed correlation between the uncertainties of both pseudo-measurements is studied, as well as possible biases with respect to the number of bins in each pseudo-measurement. Figure 5 shows the dependence of Δx/σx on the choice for the correlation coefficients c for t = 1. The modelling worsens slightly when |c| decreases, but is below about 3% with respect to the total uncertainty on the combined value for all 20,000 pseudo experiments. Also, the total uncertainty remains well modeled with only a very moderate increase of combinations with |Δ r σ | slightly different from 1, as shown in Fig. 6. The same conclusion can be drawn when the procedure described here is repeated for a different number of bins in each pseudo-measurement (not shown here). All results for 2, 4, 20, and 100 bins show a good modelling of the combined likelihood approach with respect to the central values and the total uncertainties. In general, the central result and its uncertainty are very well modelled for a large range of relative contributions from systematic uncertainties, correlations among them, and the chosen number of bins in each pseudo-measurement. Very small deviations of the order of a few percent with respect to the total uncertainty are observed. These are expected as a result of reducing the binned information of the initial measurement likelihood to one or many estimates and a corresponding Hessian. The method described here shows a similar stability for different choices of x a and x b , and of the number of uncertainties. Moreover, it is also valid for multiple estimates within one pseudo-measurement without statistical correlations between them. The case of statistical correlations is discussed separately in Sect. 3.3.

Modeling of statistical correlations
The correct modelling of statistical correlations between the estimates within a measurement is tested by generating two Relative total uncertainty of the pseudo-measurement a for different values of the threshold t for systematic uncertainties. All pseudo-measurements comprise one relative and two absolute uncertainties depicted either. Different choices for x α μ were tested and confirm a good modelling with respect to d.

Relative uncertainties
The modelling of relative uncertainties is studied by generating two pseudo-measurements, each of them with one parameter to be combined, one relative uncertainty, and two absolute uncertainties. The relative uncertainty applies to all bins in the same way and will therefore not receive constraints. In consequence, it will be dominant. Thus, the total uncertainty of each pseudo-measurement will differ from the dependence on t previously illustrated in Fig. 4. A total of 2 × 5000 pseudo-measurements are generated. Figure 9 shows the relative uncertainty of pseudo-measurement a, including one relative uncertainty, as a function of t. For t larger than 0.15, the direct likelihood combination shows instabilities in some cases, likely related to the Gaussian penalty terms, while log-normal terms would be more suitable for large relative uncertainties. As shown in Figs. 10 and 11, also when combining pseudo-measurements with contributions from relative uncertainties, central values and uncertainties are well modelled, assuming the uncertainties of one pseudo-measurement to be uncorrelated with the uncertainties of the other. The same holds true for high correlations between the pseudo-measurements.
Additionally, the method is validated using exactly one estimate per pseudo-measurement and one large relative uncertainty of +15%. The input estimates are set to 30000 + β, where β is a randomly generated value between 0 and 750. The uncertainty is assumed to be fully correlated between the pseudo-measurements. This results in asymmetric uncertainties on each pseudo-measurement and the combined value. Moreover, for this particular choice of uncertainties, the combined value can be larger than the highest input estimate. When comparing the direct likelihood combination to the  . 11 Ratio of the uncertainties onx obtained with the method described in this document and using a direct likelihood combination, shown for different values of the threshold t for systematic uncertainties. All pseudo-measurements comprise one relative and two absolute uncertainties method proposed here, also in this case no bias with respect to the central value or its uncertainties can be observed. In summary, the combination method described here does not require the original data and the full fit model, but sufficiently describes the initial measurement for a large variety of possible central values, binning choices and uncertainties. In consequence, the combination results are numerically equivalent to a using the full likelihood information, in particular in case of dominant systematic uncertainties.

Neglecting correlations
For comparison, two pseudo-measurements a and b with two bins, two systematic uncertainties, and the same parameters described in Sect. 3.2 are combined neglecting correlations between uncertainties within the same pseudo-measurement, but still considering strong correlations between pseudo- measurements a and b. This approximates the situation in which the BLUE method [6][7][8] can be used for the combination. The correlations within one pseudo-measurement are removed by inverting (D α + 1) in Eq. 2, removing the off-diagonal elements of the resulting covariance matrix, and replacing D α by the inverse of this covariance matrix minus 1. By choosing the Neyman χ 2 definition in addition, this makes this test equivalent to the BLUE method. As shown in Fig. 12, this approximation can lead to wrong individual combination results with respect to the central value when the contribution of systematic uncertainties becomes nonnegligible. Also, the uncertainty on the combined value can be severely mismodelled if the correlations within one measurement are neglected, as displayed in Fig. 13. The total uncertainty can be underestimated or strongly overestimated, in particular if it is dominated by systematic uncertainties. Therefore, it is crucial to model these correlations consis-tently when performing a combination of results obtained in simultaneous fits of systematic uncertainties and the quantity to be determined.

Program installation and user interface
The method described in Sect. 2 is implemented in the dedicated Convino program for the combination of experimental results. The source code can be found at https://github.com/ jkiesele/Convino/releases. It can be compiled using make with gcc version 4.9 or newer, or clang 8.0.0 or newer (OSX) and ROOT 6 installed on the system. Other versions might be sufficient but are not tested. The measurements and the configuration for the combination are contained in human-readable text files. Alternatively, a C++ library is provided with the software package, providing an interface to C++ standard-library or ROOT classes, the latter commonly used in high-energy physics. Both interfaces are described in the following, starting with the text-based interface. The discussion of the text-based interface serves as reference for the description of the C++ interface.

Text-based interface
The "convino" executable can be found in the base directory after compiling. It prints usage information and a list of options if the -h option is specified. Other options are: s perform correlation scan p save scan plots as .pdf in addition to a .root file d switch on debug printout --neyman uses a Neyman χ 2 instead of the Pearson χ 2 --prefix defines a prefix for all output files and directories In addition to the options, a text file is passed to the executable. It is referred to as base file in the following and is described in Sect. 5.1.2. Each measurement comprising one or a set of estimates is described in a measurement file. Well documented examples for both types of files are provided in the examples directory and should be consulted alongside this manual.

Measurement file
Each measurement file consists of blocks. Each block describes estimates or uncertainties. They are defined by a Hessian, a correlation matrix together with constraints, or a set of orthogonal uncertainties. The latter should be provided in the following format: The uncertainties sys_XX on the estimates esti-mate_XX are given in absolute values. The keyword stat is reserved for the statistical uncertainty. The uncertainties and their effect on the estimates in a measurement using a simultaneous nuisance parameter fit technique are described either by a Hessian or a correlation matrix. The Hessian must be written in the following form: while the correlation matrix has to include additional information about the constraints on the parameters. These constraints are given in units of 1σ variations for the systematic uncertainties, such that a value of 1 corresponds to no reduction and lower values indicate a constraint from the fit to the data. For estimates, the constraints are given in absolute units and correspond to the total uncertainties. In both cases, they are defined in parentheses, such that the correlation matrix is of the format:   Here, n_estimates gives the number of estimates.

Base file
The first block of the base file defines the number of measurement files (nFiles) to be considered for the combination and the corresponding file names. An example is given below: The files must be in the same directory as the base file. The second block defines the observables, the estimates should be combined to: [observables] combined_a = estimate_a1 + estimate_a2 combined_b = estimate_b1 + estimate_b2 [end observables] Here, estimate_a1 and estimate_a2 should be combined to combined_a, and similarly for esti-mate_b1 and estimate_b2. The number of estimates that should be combined to a single quantity is not limited, as well as the number of combined values. This makes it possible to combine simultaneously e.g. a large amount of bins from differential cross sections from various channels and experiments. However, in this case, the C++ interface is probably more practical.
The last block describes the correlations that should be assigned using the following syntax.
Here, a correlation coefficient of 0.2 is assigned between sys_b1 and sys_c2 and -0.3 between sys_c1 and s-ys_d2. The correlation assumptions between the parameters can be scanned in an automated way. In this case, the following syntax is used to define the scan ranges: In this case, the scan range for a single coefficient can start from positive values to negative values to allow accounting for anti-correlations between the parameters that are scanned simultaneously.
Correlation matrices are positive definite by definition, a correlation matrix C with large off-diagonal entries might lose this property if ill-posed assumptions are made, such as: In this case, the program exits and it is strongly advised to revise the plausibility of the correlation assumptions. The results of the combination are saved in the output file result.txt, or <prefix>_result.txt in case a prefix is specified. The output file contains the original input correlations, the combination results, the minimum χ 2 , and pulls and constraints on all parameters. The output of the scan, including all correlation matrices, is saved in the file scan_result.txt.

C++ interface
The C++ interface is optimized for the combination of differential distributions and provides three basic classes which will be described in the following: the class measurement, which is analogous to a measurement file discussed in the previous Section, the class combiner to perform the combination, and a class combinationResult that collects the output of the combination. The measurement class and the combinationResult class provide interfaces to C++ standard library std::vector<double> or alternatively to ROOT histograms and graphs. An example of the usage is provided in bin/differentialExample.cpp. Any cpp file that will be placed in the bin directory will be compiled automatically when running make. Alternatively, the compilation of the Convino package will create the library libconvino.so that can be linked against. The header files can be found in the include directory. Each class is documented in the corresponding header file. Therefore, the documentation here is limited to the general usage.

Measurement class
The measurement class provides the possibility to define a set of estimates, their statistical correlations and systematic uncertainties. Each object can only contain one set of estimates at once. In case the information is read from a ROOT TH1 histogram, each measurement class object can contain only one nominal histogram.
For a measurement with orthogonal uncertainties, the following procedure should be applied: the nominal values are set using the function setMeasured. Systematic uncertainties can be added in a second step to the measurement object with addSystematics. The type of each uncertainty is defined using the function using setParameter-Type after all uncertainties have been added. Here, it is recommended to use the parameter name to identify the correct uncertainty. In a last step, statistical correlations between the estimates can be set using the method setEstimateCorrelation.
If a measurement comprises correlated uncertainties, the corresponding measurement object should be configured using the function setHessian, which defines the uncertainties and estimates at once. Additional orthogonal uncertainties can be added in a subsequent step using addSystematics.

Combiner class
Once the individual measurement objects are defined, they are added to a combiner object using the function addMeasurement. For the following combination, it is assumed that the entries of each measurement in the same bin or with the same vector index should be combined. It is not possible to combine a number of estimates from one measurement object with a different number of estimates from another. The correlation assumptions are defined with set-SystCorrelation. It is advised to use the uncertainties names as input for unambiguous identification.
The combination is initiated by calling the method combine, which returns a combinationResult class object.

CombinationResult class
The combinationResult class is a container for all information regarding the inputs to the combination, the correlation matrices, and the combined values as well as the post-combination correlation matrices, pulls and constraints. If differential distributions are combined, the result can be fed back to a ROOT TH1 object or a TGraphAsymmErrors using the functions fillTH1 or fillTGraphAsymm-Errors.

Summary
The combination method presented in this document allows combining measurements obtained with simultaneous nuisance parameter fits consistently, taking into account the constraints from the data as well as correlations between systematic uncertainties within each measurement. In contrast to the optimal case of a direct likelihood combination, based on the product of the individual likelihoods of each measurement, the method does not require the full knowledge of the original data and the fit models. This information would also be required by other commonly used combination methods, however, it is publicly available only in rare cases. It is shown that not accounting for correlations between uncertainties within the same measurement can lead to nonnegligible deviations from the combined likelihood approach with respect to the combined value and its uncertainty.
The method described here does not introduce such deviations and relies on the central results and their covariances or Hessians, only, which makes it applicable to a significantly larger variety of combinations. An extensive validation is performed using pseudo-measurement with varying contributions of statistical and systematic uncertainties, correlation assumptions, binning choices, and pri-or statistical correlations. All obtained results and uncertainties are numerically equivalent to a direct likelihood combination. Only for measurements strongly limited by statistical precision, the same known caveats as in other χ 2 or least-squares-based approaches (e.g. the BLUE method) apply. In addition, the Convino program is presented. It is developed to perform combinations using the method described here and provides a text-based and a C++ user interface. The text-based user interface provides an automatic scan of correlation assumptions and creates the corresponding figures for graphical representation.