| 1 | |
| 2 | |
| 3 | == EW Higgs production at NNLO in QCD == |
| 4 | |
| 5 | The aim of this project is to compute the total VBF cross section at NNLO using the structure function approach. This approach (see hep-ph/9206246) consists in parametrizing the proton-weak boson interaction vertex with the DIS structure functions, and is exact at LO and NLO, while at NNLO neglects terms suppressed by 1/Nc^2. Since 1/Nc^2 is of the same order of the strong coupling, these neglected terms can be considered NNNLO. Furthermore the terms which are neglected can be estimated to be numericaly small also because they do not contain any enhancement by Sudakov (soft) logarithms. |
| 6 | |
| 7 | As shown in (hep-ph/9206246), the NLO corrections are of order 5-10% with respect to the LO cross section. We can guess-estimate the NNLO corrections to be of order 1% or less. |
| 8 | === Structure functions === |
| 9 | |
| 10 | In order to compute the NNLO cross-section, we need NNLO DIS structure functions. |
| 11 | ==== New in v1.2 ==== |
| 12 | |
| 13 | The file libstruct_12.f contains code for the structure functions, checked until NLO. |
| 14 | |
| 15 | LO and NLO structure fucntions can esaily be found in literature (eg. the Ellis book), while for the NNLO structure functions we use some parametrizations by A.Vogt which can be found at http://www.liv.ac.uk/~avogt/coeff.html#c23l2 . |
| 16 | |
| 17 | The structure function library works as follow: |
| 18 | |
| 19 | Each structure function Fi (i=1,2,3,L) has the following declaration: |
| 20 | |
| 21 | double precision function Fi(double precision x,double precision q,integer ord,integer v,double precision zz) |
| 22 | |
| 23 | where: |
| 24 | * x is the Bjorken x for the DIS |
| 25 | * q is the square root of the virtuality Q^2 of the boson vector (in !GeV) |
| 26 | * ord is the perturbative order (1=LO, 2=NLO, 3=NNLO). Note that, eg. if ord=2 only the O(alpha_s) contribution is given and not the O(alpha_s^0)+O(alpha_s) one |
| 27 | * v is the index for the vector boson, considered, as in DIS, entering the proton (1=W-, 2=W+, 3=Z0) |
| 28 | * zz is the variable used for convolutions beyond the LO. It is supposed to be generated between 0 and 1. The convolution can be performed using any external integration routine. At LO it can be set to any value, since it is not used. |
| 29 | |
| 30 | Some paramter are passed through common blocks: |
| 31 | * /nflav/nf is the number of active flavours |
| 32 | * /prec/eps is a parameter used to avoid the singularities of the coefficient functions at x=1. Suitable values are 1d-10 or less. |
| 33 | |
| 34 | In the main progam there should be a call to the subroutine SetCouplings() in order to initialize the couplings to the various vectors. |
| 35 | |
| 36 | The libstruct library uses LHAPDF pdf sets, that have to be initialized in the main program. |
| 37 | |
| 38 | Until now, all the scales are set equal to Q^2. |
| 39 | ==== New in v1.3 ==== |
| 40 | |
| 41 | Due to some troubles with nf and setcouplings, in v1.3 no call to setcouplings is needed, neither nf setting in the main program. |
| 42 | |
| 43 | The structure functions are automatically computed with 5 active flavours for the Z0 and 4 active flavours for the W+/- (top quark decoupled). |
| 44 | ==== New in v1.4 ==== |
| 45 | |
| 46 | It has been fixed an issue concerning the singlet and non singlet distribution. The NNLO coefficient function has been checked. Should be reliable up to NNLO |
| 47 | |
| 48 | -- Main.MarcoZaro - 2009-12-17 |
| 49 | ==== New in v2.0 ==== |
| 50 | |
| 51 | The full (factorzization and renormalization) scale dependence has been implemented until NNLO. It has been checked at NLO with the code of Spira, while it is being checked at NNLO (the mellin moments of the coefficent functions are ok, the flavour combination is to be checked). |
| 52 | |
| 53 | -- Main.MarcoZaro - 2010-02-05 |
| 54 | ==== Scale dependence ==== |
| 55 | |
| 56 | The tar-ball conv.tgz contains code with closed form expression for the necessary convolutions P_ij'''P_ij and c_i'''P_ij at NNLO. Two check programs are attached as well. |
| 57 | |
| 58 | Sven - 2010-01-07 |
| 59 | ==== To Do ==== |
| 60 | * check scale dependence |
| 61 | === VBF cross section intgration === |
| 62 | |
| 63 | The file vbf_hadr_xx.f contains the code for the computation of the hadronic VBF cross-section. It has been checked up to NLO with the code by Michael Spira (spira.zip, see attachment table). |
| 64 | |
| 65 | In order to be properly compiled, it requires libstruct_xx.f (xx=13 or later) and the integration routine vegas.f . All these files and a Makefile are included in the archive vbf.tar |
| 66 | |
| 67 | Up to now the parameters have to be set inside the code. The most relevant are |
| 68 | * double precision sqrts, the total center of mass energy (in GeV) |
| 69 | * double precision mh, the higgs mass |
| 70 | * integer icomp, fixes the accuracy of the computation: 1 LO+intf, 2 NLO, 3 NLO "improved", 4 NNLO, 5 NNLO "improved", where "improved" means that the alpha_s power counting is performed in the structure functions and not for the whole cross section. (E.g. the NLO "improved" cross-section contains also the NLOxNLO contribution of the structure function expansion, which is ont inculded in the NLO) |
| 71 | ==== New in v2.0 ==== |
| 72 | |
| 73 | The parameters are set outside the program, using the cards param_card.dat for Standard Model parameters (masses, widths,...) and run_card.dat for the parameters of the specific run (energy, higgs mass, s-channel inclusion, pdfs,...). |
| 74 | |
| 75 | Furthermore different pdf sets can be used for LO, NLO and NNLO computations. |
| 76 | |
| 77 | The icomp variable now runs from 1 to 6. |
| 78 | 1. LO, PDF LO |
| 79 | 1. NLO, PDF NLO |
| 80 | 1. NNLO, pdf NNLO |
| 81 | 1. NLO imp, pdf NLO |
| 82 | 1. NLO imp, pdf NNLO |
| 83 | 1. NNLO imp, pdf NNLO |
| 84 | |
| 85 | The bash script run.sh allow to take several runs with different parameters, the run_card being made by the fortran program makeruncard.f at each run. The output is collected in the output file out.dat |
| 86 | |
| 87 | -- Main.MarcoZaro - 2010-01-22 |
| 88 | ==== New in v2.5 ==== |
| 89 | |
| 90 | Some speed optimization has been performed. Basically, in v2.0 the evolvePDF routine was called many times for the single structure function call, basically one time for each flavour combination (singlet, gluon,...). |
| 91 | |
| 92 | In order to speed up the code, the evolvePDF routine is called once in the structure function. This turned into a moderate gain in the program runtime (about 10-15%) with respect v2.0. |
| 93 | |
| 94 | Note that the interface of the structure function is not changed in this version with respect to v2.0. |
| 95 | ==== New in v3.0 ==== |
| 96 | |
| 97 | In this version the evolvePDF routine is called as few times as possible. Basically it is called in the function me(...), and then the result is pssed as an argument to the structure functions, whose interface has been modified. |
| 98 | |
| 99 | The runtime gain with respect v2.0 is of about 35%. |
| 100 | ==== New in v3.5 ==== |
| 101 | |
| 102 | We improved the running speed by defining a subroutine sethpl(x) which computes once the harmonic polylogs and stores them in a common block. In this version this subroutine is invoked inside the structure functions, so only the files libstruct.f and libscale.f have been modified. The running time dramatically reduced: RT[3.5]/RT[3.0]=35% |
| 103 | |
| 104 | Also the correct top-triangle in the t-channel is included in this version (even if not called in the integration). |
| 105 | ==== New in v4.0 ==== |
| 106 | |
| 107 | We still improved the running time by making the call to sethpl external to the structure functions, i.e. it is called in the me() function in the file vbf.f. RT[4.0]/RT[3.0]=30% |
| 108 | |
| 109 | -- Main.MarcoZaro - 2010-03-19 |
| 110 | |
| 111 | === S-channel === |
| 112 | |
| 113 | The VBF does not include the s-channel contribution to the Higgs production, which is basically Drell-Yan plus the decay of the virtual vector particle into higgs and real vector, and the subsequent hadronic decay of the latter. This contribution cannot be computed within the framework of the structure function approach, but, at least for the initial state, the strong corrections are just the Drell-Yan ones, and so it can be easily included in our code. |
| 114 | |
| 115 | (See for example [http://arxiv.org/abs/hep-ph/0307206 Brein, Djouadi, Harlander]). |
| 116 | |
| 117 | The s channel (up to NLO in v 2.0), is included in the file libschannel_xx, basically using the Drell-yan approach. |
| 118 | |
| 119 | The flag isch in the vbf_hadr main program, allows to include it in the computation (0 no, 1 yes). |
| 120 | ==== To do ==== |
| 121 | * NNLO s-channel |
| 122 | === Interference === |
| 123 | |
| 124 | In order to reach a precision of order 1% we should also consider the interference terms between WW and ZZ fusion, as quell as the ones coming for VBF-like diagrams with the s-channel diagram which are not taken into account in the structure function approach. |
| 125 | |
| 126 | The interference between VBF-like diagrams is known to be small, about some %s, and hence for us it is sufficient to include it at LO. |
| 127 | |
| 128 | The one coming from VBF x S-channel is zero in the zero W width limit, and it can be estimated as of the order of the ratio between the W-width and its mass. |
| 129 | ==== libintf_15.f ==== |
| 130 | |
| 131 | The file libintf_15.f contains the LO interferences between the VBF-like diagrams (eg. ud->udh, uu->uuh) as well as the interferences between VBF and S-channel (udbar->udbarh). It as been checked with Madgraph point by point in the phase space and at the integrated level. |
| 132 | |
| 133 | The interference impact is at the per-mille level or less at 14 TeV |
| 134 | |
| 135 | The files 14tev_LO_intf.dat and 14tev_LO_intf_sch.dat in the tarball intf_bench.tar contain some data, repsectively including and excluding the interference of the S-channel diagram with the VBF like ones. The number in the files are the higgs mass, the LO VBF cross section (no intf), the interference, the total LO cross section and the interference impact. |
| 136 | |
| 137 | -- Main.MarcoZaro - 2010-01-14 |
| 138 | === Things we know we neglect and we should have an estimate of: === |
| 139 | * alpha_s^2/Nc^2 terms: |
| 140 | * These come from "interference" diagrams. One way to try to estimate them is to consider the QED^2 corrections to VBF with a Z exchange. |
| 141 | * alpha_s^3 alpha_ew^2 contributions: |
| 142 | * See paper by Jeppe, Thomas, and Jennifer: ["http://arxiv.org/abs/0709.3513"] |
| 143 | * top quark effects: |
| 144 | * In particular we neglect diagrams including a top box coupled to two gluons, a Z and the higgs, and their interference with diagrams included in our calculation. Those diagrams are suppressed in the mt->infty limit. However, this contribution can be substantial at the LHC. (See for example [http://arxiv.org/abs/hep-ph/0307206 Brein, Djouadi, Harlander]), and have been computed in this work by [http://dx.doi.org/10.1103/PhysRevD.42.2253 Kniehl] (see also attachment kniehl.pdf) |
| 145 | * |
| 146 | * b quark initiated charged current processes (i.e. no b->c, b->u mixing). Since b->u << b->c, we just need to check the effect of the latter. Being totally inclusive, we already take into account CKM mixing for the other flavours |
| 147 | * interferences beyond LO (ok if interferences are already small) |
| 148 | * Non DY-like diagrams in the s-channel (related to top quark effects, ex. gg->VH via top box or triangle, (see above) |
| 149 | * ||s-ch||^2 -> branching ratio into hadrons: we should understand if take the experimental value, or the fixed order QCD computation |
| 150 | |
| 151 | -- Main.MarcoZaro - 2009-12-18 |
| 152 | === Attachments === |
| 153 | * [attachment:libstructf_12.f libstructf_12.f]: Structure functions library. V1.2 |
| 154 | * [attachment:vegas.f vegas.f]: The vegas integration routine used in vbf_hadr |
| 155 | * [attachment:libstructf_13.f libstructf_13.f]: Structure functions library. V1.3 |
| 156 | * [attachment:vbf.tar vbf.tar]: the vbf code |
| 157 | |
| 158 | * [attachment:spira.zip spira.zip]: M. Spira code for VBF |
| 159 | |
| 160 | * [attachment:libstructf_14.f libstructf_14.f]: The structure function library, with coefficient function checked up to NNLO. All previous version are obsolete |
| 161 | |
| 162 | * [attachment:vbf_14.tar vbf_14.tar]: the 1.4 version of the vbf_hadr code, including libstructf_14 |
| 163 | |
| 164 | * [attachment:conv.tgz conv.tgz]: All convolution for the NNLO scale dependence |
| 165 | |
| 166 | * [attachment:vbf_15.tar vbf_15.tar]: LO interference included |
| 167 | |
| 168 | * [attachment:intf_bench.tar intf_bench.tar]: interference data |
| 169 | |
| 170 | * [attachment:vbf.pdf vbf.pdf]: Structure functions for EW gauge boson exchange to NNLO |
| 171 | |
| 172 | * [attachment:vbffinal.pdf vbffinal.pdf]: Structure functions for EW gauge boson exchange to NNLO (final version) |
| 173 | |
| 174 | * [attachment:COEFFX.tar.gz COEFFX.tar.gz]: van Neerven code for DIS structure functions to NNLO |
| 175 | |
| 176 | * [attachment:coefmom.f coefmom.f]: my adaptation of the main van Nerven program to compute a few Mellin moments (compiles with gfortran) |
| 177 | |
| 178 | * [attachment:vbf_20.tar vbf_20.tar]: the 2.0 version of the vbf code |
| 179 | |
| 180 | * [attachment:kniehl.pdf kniehl.pdf]: knihel paper about pp->ZH @NNLO |
| 181 | |
| 182 | * [attachment:vbf_20.tar vbf_20.tar]: the 2.0 version of the vbf code |
| 183 | |
| 184 | * [attachment:vbf_20.tar vbf_20.tar]: the vbf_20 version of the code, with the bash scipts.. |
| 185 | |
| 186 | * [attachment:vbf_20.tar vbf_20.tar]: the vbf_20 version of the code, with the bash scipts and the run results in the file results.dat |
| 187 | |
| 188 | * [attachment:vbf_20.tar vbf_20.tar]: the vbf_20 version of the code, with the bash scipts and the run results (version 01) in the file results.dat |
| 189 | |
| 190 | * [attachment:vbf_20.tar vbf_20.tar]: STRUCTURE FUNCTIONS FULLY CHECKED AT NNLO |
| 191 | |
| 192 | * [attachment:vbf_25.tar vbf_25.tar]: the 2.5 version of the vbf code. Speed improvement of about 10-15% with respect to the 2.0 version |
| 193 | |
| 194 | * [attachment:vbf_30.tar vbf_30.tar]: the 3.0 version of the vbf code. Speed improvement of about 35% with respect to the 2.0 version |
| 195 | |
| 196 | * [attachment:res100224.tar res100224.tar]: the results obtained with 3.0 version, both scale and pdf dependence are included |
| 197 | |
| 198 | * [attachment:vbfcolor.pdf vbfcolor.pdf]: Notes about the color structure at NNLO and discussion of the color suppressed contributions |
| 199 | |
| 200 | * [attachment:vbf_35.tar vbf_35.tar]: 3.5 version of the vbf code. Top triangle in t-channel corrected and hplog speed improved |
| 201 | |
| 202 | * [attachment:vbf_40.tar vbf_40.tar]: 4.0 version of the vbf code. Hplog speed improved (call to sethplog external to the structure functions) |
| 203 | |
| 204 | |
| 205 | |
| 206 | |
| 207 | |
| 208 | |
| 209 | |
| 210 | |
| 211 | |
| 212 | |
| 213 | |
| 214 | |
| 215 | |
| 216 | |
| 217 | |
| 218 | |
| 219 | |
| 220 | |
| 221 | |
| 222 | |
| 223 | |
| 224 | |
| 225 | |
| 226 | |