Pre-approval comments v1

Questionnaire is filled in and submitted.

An svn repository (svn+ssh://svn.cern.ch/reps/cmshcg/trunk/cadi/HIG-17-028) is requested, and all datacards will be uploaded once the repository is created. Note that almost all datacards will be updated once the new hgg binning is obtained.

Since this analysis concerns a combination, filling in the entire Muon questionnaire would be an exercise in copying information from the Hgg and Hzz entries. It is preferred to simply refer to the Muon questionnaire entries for the individual analyses: HIG-16-041 , HIG-17-025 .

A line indicating the SM prediction is added to all the plots concerning combinations of differential cross sections. The 2D coupling scans already had SM predictions.

The normalization of the SM line is set to the YR4 N3LO cross section for gluon fusion: 48.52 fb.

pT, for smH and ggH+xH (split):

nJets, pT_jet and rapidity:

The uncertainty denoted by the vertical line corresponds to the full uncertainty (statistical and systematic), whereas the filled area denotes only the statistical uncertainty (the total uncertainty is dominated by the statistical uncertainty; see Question 5).

The uncertainties previously calculated concerned the combined uncertainty of statistical and systematic uncertainties. In order to calculate only the statistical uncertainties, the scans were repeated freezing all the nuisance parameters to their best-fit-value. The uncertainties split into stat and stat+syst are shown in the following plots, for the combined spectra.

The full set of scans from which the uncertainties were obtained is here: http://tklijnsm.web.cern.ch/tklijnsm/differentials2017/ptCombination_v3/preapp/statsyst_scans/ . In the majority of these scans, the systematic uncertainty is very small with respect to the statistical uncertainty.

The numerical values are also given in the tables in Question 6.

For all the observables, tables of their mu (ratio to SM) and cross section where computed. The resulting tables are displayed below, and will also be included in the paper draft.

The tables contain uncertainties in terms of ratio to SM (mu).

A separate scan was performed in which all the systematic nuisances were frozen to the best fit, yielding the only statistical uncertainties. This shape is shown in grey.

I performed this cross check for all the observables we are doing. The results agree within uncertainties, and are shown below.

Scanning the ratio of BRs per pT bin would only have been possible if both hgg and hzz had identical binning schemes. Currently however, floating a BR-modifier per hgg-bin and hzz-bin yields too many degrees of freedom (8 modifiers for hgg + 8 modifiers for hzz > 8 reco bins hgg + 5 reco bins hzz).

A 2-sigma line is included in the plot; note that 2-sigma (95.45%) differs slightly from the asked 95% CL. Also the numerical result is printed in a more clear way.

Observed:

Expected:

I added a point at (0,0). The contribution that is left is due to the kappa_top term, that is fixed to SM.

In black is shown the 'nominal' result: Branching ratios fixed to their SM value, and the normalization fixed to the ggH cross section calculated in Yellow Report 4. The result obtained by making the branching ratios depend on the couplings is shown in blue; the result is more constrained due to the constraining effect of the total width. In red a further modification is applied, where also the overall normalization is floated in the fit. As expected, this is less constraining than keeping it fixed.

The results shown here were computed using the Asimov (expected) dataset.

When the constraint from the normalization is dropped, the shape depends only on the ratio of kappa_t and kappa_g (this is a mathematical consequence of the parametrization). The plot below confirms this expectation: Different levels of exclusion are only obtained by varying the ratio.

In black is shown the 'nominal' result: Branching ratios fixed to their SM value, and the normalization fixed to the ggH cross section calculated in Yellow Report 4. The result obtained by making the branching ratios depend on the couplings is shown in blue; the result is more constrained due to the gamma-gamma loop effects (the gamma-gamma width goes to zero around kappa_t ~ 2.0).

In red a further modification is applied, where also the overall normalization is floated in the fit. As expected, this is less constraining than keeping it fixed. The discrimination due to the shape is mostly on kappa_t, whereas on kappa_g hardly a clear constraint can be made.

The results shown here were computed using the Asimov (expected) dataset.

I would suggest to use "c_g" here, which is what the authors behind the kappa_t/c_g originally called it. Once agreed upon I will change all plotting labels.

As the hgg cards are not yet migrated to the new binning scheme, I performed a combination of hzz and hbb:

The result is comparable to that of hgg alone. As soon as the new hgg datacards are in, I will restart all the relevant scans. The results shown here were computed using the Asimov (expected) dataset.

So far, no additional theory uncertainty for the mtop treatment above 350 GeV has been included.