Changes between Initial Version and Version 1 of VBF


Ignore:
Timestamp:
03/20/12 16:17:47 (8 years ago)
Author:
trac
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • VBF

    v1 v1  
     1
     2
     3== EW Higgs production at NNLO in QCD ==
     4
     5The 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
     7As 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
     10In order to compute the NNLO cross-section, we need NNLO DIS structure functions.
     11==== New in v1.2 ====
     12
     13The file libstruct_12.f contains code for the structure functions, checked until NLO.
     14
     15LO 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
     17The structure function library works as follow:
     18
     19Each structure function Fi (i=1,2,3,L) has the following declaration:
     20
     21double precision function Fi(double precision x,double precision q,integer ord,integer v,double precision zz)
     22
     23where:
     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
     30Some 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
     34In the main progam there should be a call to the subroutine SetCouplings() in order to initialize the couplings to the various vectors.
     35
     36The libstruct library uses LHAPDF pdf sets, that have to be initialized in the main program.
     37
     38Until now, all the scales are set equal to Q^2.
     39==== New in v1.3 ====
     40
     41Due to some troubles with nf and setcouplings, in v1.3 no call to setcouplings is needed, neither nf setting in the main program.
     42
     43The 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
     46It 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
     51The 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
     56The 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
     58Sven - 2010-01-07
     59==== To Do ====
     60   * check scale dependence
     61=== VBF cross section intgration ===
     62
     63The 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
     65In 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
     67Up 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
     73The 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
     75Furthermore different pdf sets can be used for LO, NLO and NNLO computations.
     76
     77The 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
     85The 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
     90Some 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
     92In 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
     94Note that the interface of the structure function is not changed in this version with respect to v2.0.
     95==== New in v3.0 ====
     96
     97In 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
     99The runtime gain with respect v2.0 is of about 35%.
     100==== New in v3.5 ====
     101
     102We 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
     104Also 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
     107We 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
     113The 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
     117The s channel (up to NLO in v 2.0), is included in the file libschannel_xx, basically using the Drell-yan approach.
     118
     119The 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
     124In 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
     126The 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
     128The 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
     131The 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
     133The interference impact is at the per-mille level or less at 14 TeV
     134
     135The 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