mcremone / decaf

0 stars 12 forks source link

NNLO uncertainties #23

Closed mcremone closed 1 year ago

mcremone commented 4 years ago

In principle, composed by three pieces: epsilonQCD, epsilonEW, and epsilonMIX.

epsilonQCD

It describes uncertainties related to variations of the renormalization and factorization scales which are performed to estimate the uncertainty of the theoretical prediction due to missing higher-order contributions.

The uncertainty is split into three parts:

All three nuisance parameters need to be treated uncorrelated, however each parameter is correlated for all V+jets processes and all bins of boson pT.

Questions:

epsilonEW

It describes uncertainties related to missing even higher-order(NNLO) contributions:

Questions:

epsilonMIX

It is an additional uncertainty regarding mixed contributions that cannot be described with the multiplicative or factorized approach. This nuisance parameter is chosen to be completely uncorrelated between the different processes.

Questions

So far I see at least three nuisances that for sure need to be accounted for in the fit, that are epsilon2, epsilon3, and epsilonMIX. It needs to be clarified if they are shape or normalization effects. We are going to have a different nuisance per V+jets process, therefore we will have: epsilon2^{V}, epsilon3^{V}, and epsilonMIX^{V}, where V can be W, Z, or gamma. The remaining nuisances, epsilonQCD1, epsilonQCD2, epsilonQCD3, and epsilon1, if they not only move in correlated fashion for all V+jets processes, but they also move by the same quantity, they are going to cancel out in the transfer factor ratio and therefore they can be excluded from the pool of uncertainties. If they do not, we are going to have one nuisance for all processes.

mcremone commented 4 years ago

Answers from Michael:

It turns out (see slides) that epsilonEW1, epsilonEW2, epsilonEW3, epsilonMIX, and epsilonQCD3 are actually very small effects, of the order of a fraction of percent. Therefore the can be neglected.

The only two to consider are epsilonQCD1 and epsilonQCD2. Looking at the slides, it appears that:

This should be confirmed looking at the following ratios:

since QCD1 and QCD2 are correlated between all processes, we need to check only the ratios of up variations for Z with up variations for W and gamma, and similarly and for down variations.

Three things can happen:

mcremone commented 4 years ago

Renormalization and Factorization Scale Uncertainties

Need to take into account these too. We will label renormalization as muR and factorization as mF. In the files, alternate shapes are defined as:

Obtained by multiplying by 0.5 muR and muF scales for the down variations, and multiplying by 2 muR and muF scales for the up variations.

First to check would be the absolute scale of the variations with respect to the nominal:

where muR_nominal=muF_nominal=*_NNLO_NLO_nnn_nnn_n. These ratios should be computed per process, that means that the same process should appear in both the numerator and the denominator. For example, muR_down/muR_nominal for W+jets will be:

evj_NNLO_NLO_nnn_nnn_n_Weight_scale_variation_muR_0p5_muF_1p0 / evj_NNLO_NLO_nnn_nnn_n

Second, we need to compute ratios as done for the other uncertainties:

and same for muF.