Chapter 6
Input file configuration
6.1 Runtime input file parameters
MCFM execution is performed in the Bin/ directory, with syntax:
mcfm input.ini
If no command line options are given, then MCFM will default to using the file input.ini in the
current directory for choosing options. The input.ini
file can be in any directory and then
the first argument to mcfm
should be the location of the file. Furthermore, one can
overwrite or append single configuration options with additional parameters like:
./mcfm benchmark/input.ini general%part=nlo lhapdf%dopdferrors=.true.
Here specifying a parameter uses a single dash, then the section name as in the input file (see below), followed by a percent sign, followed by the option name, followed by an equal sign and the actual value of the setting.
All default settings in the input file are explained below, as well as further optional
parameters. The top level setting mcfm_version
specifies the input file version number and it
must match the version of the code being used.
The general structure of a fixedorder calculation up to NNLO is as follows:
$$\sigma ={\sigma}_{0}+\mathrm{\Delta}{\sigma}_{1}+\mathrm{\Delta}{\sigma}_{2}\phantom{\rule{0.17em}{0ex}},$$  (6.1) 
where $\mathrm{\Delta}{\sigma}_{k}$ is of order ${\alpha}_{s}^{k}$ with respect to the leading order cross section ${\sigma}_{0}$, thus representing the N${}^{k}$LO contribution to the cross section. When performing the NLO calculation using dipole subtraction its contribution to the cross section can be decomposed as,
$$\mathrm{\Delta}{\sigma}_{1}=\mathrm{\Delta}{\sigma}_{1}^{v}+\mathrm{\Delta}{\sigma}_{1}^{r}\phantom{\rule{0.17em}{0ex}}.$$  (6.2) 
$\mathrm{\Delta}{\sigma}_{1}^{v}$ includes virtual (loop) contributions, as well as counterterms that render them finite. $\mathrm{\Delta}{\sigma}_{1}^{r}$ includes contributions from diagrams involving real parton emission, again with counterterms to make them finite. Only the sum of $\mathrm{\Delta}{\sigma}_{1}^{v}$ and $\mathrm{\Delta}{\sigma}_{1}^{v}$ is physical.
This contribution can also be computed using a slicing method with the corresponding decomposition,
$$\mathrm{\Delta}{\sigma}_{1}^{a}=\mathrm{\Delta}{\sigma}_{1}^{a,<}+\mathrm{\Delta}{\sigma}_{1}^{a,>}\phantom{\rule{0.17em}{0ex}}.$$  (6.3) 
$a$ labels the slicing resolution variable, which in MCFM can be either 0jettiness, ${q}_{T}$ (of a colorsinglet system) or ${p}_{T}^{{j}_{1}}$ (lead jet ${p}_{T}$) (thus corresponding to a jet veto). $\mathrm{\Delta}{\sigma}_{1}^{a,<}$ is termed the belowcut slicing contribution which is computed by the means of a factorization theorem and includes loop contributions. $\mathrm{\Delta}{\sigma}_{1}^{a,>}$ is the abovecut contribution containing radiation of an additional parton. Only the sum $\mathrm{\Delta}{\sigma}_{1}^{a}$ is physical and contains a dependence on the slicing resolution variable ${a}_{\text{cut}}$ that tends to zero as ${a}_{\text{cut}}\to 0$
At NNLO only slicing calculations are available. The decomposition is,
$$\mathrm{\Delta}{\sigma}_{2}^{a}=\mathrm{\Delta}{\sigma}_{2}^{a,<}+\mathrm{\Delta}{\sigma}_{2}^{a,v>}+\mathrm{\Delta}{\sigma}_{2}^{a,r>}\phantom{\rule{0.17em}{0ex}}.$$  (6.4) 
$\mathrm{\Delta}{\sigma}_{2}^{a,<}$ is the belowcut slicing contribution containing 2loop contributions. $\mathrm{\Delta}{\sigma}_{1}^{a,v>}$ is the abovecut contribution containing loop corrections to radiation of an additional parton. $\mathrm{\Delta}{\sigma}_{1}^{a,r>}$ is the abovecut contribution representing radiation of up to two additional partons. Only the sum $\mathrm{\Delta}{\sigma}_{2}^{a}$ is physical and contains a dependence on the slicing resolution variable ${a}_{\text{cut}}$ that tends to zero as ${a}_{\text{cut}}\to 0$
The type of computation that is performed depends on the parameter part
in the
general
section. The list of possible values, and the associated meaning, is shown in
Tables 6.1 and 6.2. They can also be listed by setting part
equal to help in the input
file
.
part that correspond
to performing a fixedorder calculation.  

description 
lo/lord 
${\sigma}_{0}$ 
virt 
$\mathrm{\Delta}{\sigma}_{1}^{v}$ 
real 
$\mathrm{\Delta}{\sigma}_{1}^{r}$ 
nlocoeff/totacoeff 
$\mathrm{\Delta}{\sigma}_{1}$ 
nlo/tota 
${\sigma}_{0}+\mathrm{\Delta}{\sigma}_{1}$. For photon processes that include fragmentation, nlo also includes the calculation of the fragmentation (frag) contributions. 
frag 
Processes 280, 285, 290, 295, 300302, 305307, 820823 only, see sections 13.67, 13.72 and 13.73 below. 
nlodk/todk 
Processes 114, 161, 166, 171, 176, 181, 186, 141, 146, 149, 233, 238, 501, 511 only, see sections 13.39 and 13.41 below. 
snloR 
$\mathrm{\Delta}{\sigma}_{1}^{a,>}$ 
snloV 
$\mathrm{\Delta}{\sigma}_{1}^{a,<}$ 
snlocoeff/scetnlocoeff 
$\mathrm{\Delta}{\sigma}_{1}^{a}$ 
snlo/scetnlo 
${\sigma}_{0}+\mathrm{\Delta}{\sigma}_{1}^{a}$ 
nnloVVcoeff 
$\mathrm{\Delta}{\sigma}_{2}^{a,<}$ 
nnloRVcoeff 
$\mathrm{\Delta}{\sigma}_{2}^{a,v>}$ 
nnloRRcoeff 
$\mathrm{\Delta}{\sigma}_{2}^{a,r>}$ 
nnloVV 
$\mathrm{\Delta}{\sigma}_{1}^{a,<}+\mathrm{\Delta}{\sigma}_{2}^{a,<}$ 
nnloRV 
$\mathrm{\Delta}{\sigma}_{1}^{a,>}+\mathrm{\Delta}{\sigma}_{2}^{a,v>}$ 
nnloRR 
$\mathrm{\Delta}{\sigma}_{2}^{a,r>}$ 
nnlocoeff 
$\mathrm{\Delta}{\sigma}_{2}^{a}$ 
nnlo 
${\sigma}_{0}+\mathrm{\Delta}{\sigma}_{1}+\mathrm{\Delta}{\sigma}_{2}^{a}$ 




part that correspond
to performing a calculation including largelog resummation.  

description 
resLO 
NLL resummed and matched 
resonlyLO 
NLL resummed only 
resonlyLOp 
NLLp resummed only 
resexpNLO 
NNLL resummed expanded to NLO 
resonlyNLO 
NNLL resummed 
resaboveNLO 
fixedorder matching to NLO 
resmatchcorrNLO 
matching corrections at NLO 
resonlyNLOp 
NNLLp resummed 
resexpNNLO 
N${}^{3}$LL resummed expanded to NNLO 
resonlyNNLO 
N${}^{3}$LL resummed 
resaboveNNLO 
fixedorder matching to NLO 
resmatchcorrNNLO 
matching corrections at NLO 
resLOp 
NLLp resummed and matched 
resNLO 
NNLL resummed, matched to NLO 
resNLOp 
N${}^{3}$LL resummed, matched to NLO 
resNNLO 
N${}^{3}$LL resummed, matched to NNLO 
resNNLOp 
N${}^{3}$LLp resummed, matched to NNLO 
resonlyNNLOp 
N${}^{3}$LLp resummed 




6.1.1 General
Section general  Description


The process to be studied is given by choosing a process number, according to Tables in Section 4. $f({p}_{i})$ denotes a generic partonic jet. Processes denoted as “LO” may only be calculated in the Born approximation. For photon processes, “NLO+F” signifies that the calculation may be performed both at NLO and also including the effects of photon fragmentation and experimental isolation. In contrast, “NLO” for a process involving photons means that no fragmentation contributions are included and isolation is performed according to the procedure of Frixione [55]. 

The type of calculation to be performed. Possible values are given in Tables 6.1 and 6.2. 

When MCFM is run, it will write output to several files. The label runstring will be included in the names of these files. 

Directory for output and snapshot files 

Center of mass energy in GeV. 

The identities of the incoming hadrons may be set with these parameters, allowing simulations for both $p\overline{p}$ (such as the Tevatron) and $\mathit{pp}$ (such as the LHC). Setting ih1 equal to $+1$ corresponds to a proton, whilst $1$ corresponds to an antiproton. 

When set to .true. then all bosons are produced onshell. This is appropriate for calculations of total crosssections (such as when using removebr equal to .true., below). When interested in decay products of the bosons this should be set to .false.. 

When set to .true. the branching ratios are removed for unstable particles such as vector bosons or top quarks. See the process notes in Section 13, or the process webpages accessed via the list of processes for further details. 

Specifies whether or not to compute EW corrections for the process. Default is none. May be set to exact or sudakov for processes 31 (neutralcurrent DY), 157 (toppair production) and 190 (dijet production). For more details see section 5.3. 




6.1.2 Resummation
Section resummation  Description


If




Output
directory
for
LHAPDF
grid
files,
for
example


Input
directory
for
LHAPDF
grid
files,
for
example


Integration
range
of
purely
resummed
part,
for
example


Integration
range
of
fixedorder
expanded
resummed
part,
for
example


Lower
${q}_{T}$
cutoff
${q}_{0}$
for
the
fixedorder
part.
Typically
the
value
should
agree
with
the
lower
range
of


Parameter passed to the plotting routine to modify the transition function, see text. 




6.1.3 NNLO
Section nnlo  Description
 

Optional. This sets the value of the jettiness variable ${\tau}_{\text{cut}}$, as a multiple of the invariant mass of the Born system, i.e.
This variable separates the resolved and unresolved regions in NNLO calculations that use zerojettiness. The default value results in total inclusive cross sections with less than $1\%$ residual cutoff effects. 


Optional. Array that specifies multiple taucut values that should be sampled
on the fly in addition to the nominal taucut value. Both larger and smaller
values than the nominal one can be specified, although uncertainties for
smaller values will be large. We generally do not recommend smaller values
than the nominal one chosen with 


Optional. If 


Flag to use
${q}_{T}$
slicing, rather than 0jettiness, in the calculation of NNLO contributions.
Default is 


If 


If 


When using 0jettiness, perform the slicing cut in the centreofmass of the
color singlet system. Default is 


When using 0jettiness, include leading power corrections in the belowcut
calculation. Default is 


When using 0jettiness, only compute the power corrections to the belowcut
calculation. Default is 


This flag has two separate uses. In a fixedorder sliciing calculation, e.g. part
is equal to snlo or nnlo, the code uses
${p}_{T}^{\mathit{veto}}$
slicing, rather than 0jettiness, in the calculation of higherorder contributions.
In a resummed calculation, e.g. part is equal to resNLO or resNNLO, it enables
the use of
${p}_{T}^{\mathit{veto}}$
resummation rather that
${q}_{T}$
resummation. In this case the value of the jet veto ( 


Implements ${p}_{T}^{\mathit{veto}}$
formalism using the refactorized approach of Ref. [58]. Otherwise uses original
‘$B\times B\times S$’
factorization of belowcut crosssection into beam and soft functions (that
gives identical results). Default is 





6.1.4 PDFs
Section pdf  Description


This specifies the parton distributions used in the case when the code has been
built with 




6.1.5 LHAPDF
Section lhapdf  Description


Specifies the parton distributions used in the case when the code has been built
with 

Specifies the individual members of the parton distribution sets. A value of zero corresponds to the central value for Hessian sets. In the case when multiple sets have been specified above, each one needs a member number separated by space. 

When this is set to 




6.1.6 Scales
Section scales  Description
 

This parameter may be used to adjust the value of the renormalization scale. This is the scale at which ${\alpha}_{S}$ is evaluated and will typically be set to a mass scale appropriate to the process (${M}_{W}$, ${M}_{Z}$, ${m}_{t}$ for instance). 


This parameter may be used to adjust the value of the factorization scale and will typically be set to a mass scale appropriate to the process (${M}_{W}$, ${M}_{Z}$, ${m}_{t}$ for instance). 


This character string is used to specify whether the renormalization,
factorization and fragmentation scales are dynamic, i.e. recalculated on an
eventbyevent basis. If this string is set to ‘none’ then the scales are fixed
for all events at the values specified by renscale, facscale as well as
The type of dynamic scale to be used is selected by using a particular string for the variable dynamicscale, as indicated in Table 6.10. Not all scales are defined for each process, with program execution halted if an invalid selection is made in the input file. The selection chooses a reference scale, ${\mu}_{0}$. The actual scales used in the code are then,
Note that, for simplicity, the fragmentation scale (relevant only for processes involving photons) is set equal to the renormalization scale. In some cases it is possible for the dynamic scale to become very large. This can cause problems with the interpolation of data tables for the PDFs and fragmentation functions. As a result if a dynamic scale exceeds a maximum of $60$ TeV (PDF) or $990$ GeV (fragmentation) this value is set by default to the maximum. 


This additional option can be set to
The histograms corresponding to these different choices are included in the output file, from which an envelope of theoretical uncertainty may be constructed by the user. 


Number of additional scale variation points to choose, can be set to two or six. For two it just samples the first two variations as in Eq. 6.7. 





dynamic scale  ${\mu}_{0}^{2}$  comments 
m(34)  ${({p}_{3}+{p}_{4})}^{2}$  
m(345)  ${({p}_{3}+{p}_{4}+{p}_{5})}^{2}$  
m(3456)  ${({p}_{3}+{p}_{4}+{p}_{5}+{p}_{6})}^{2}$  
sqrt(M^ 2+pt34^ 2)  ${M}^{2}+{({\overrightarrow{{p}_{T}}}_{3}+{\overrightarrow{{p}_{T}}}_{4})}^{2}$  $M=$ mass of particle 3+4 
sqrt(M^ 2+pt345^ 2)  ${M}^{2}+{({\overrightarrow{{p}_{T}}}_{3}+{\overrightarrow{{p}_{T}}}_{4}+{\overrightarrow{{p}_{T}}}_{5})}^{2}$  $M=$ mass of particle 3+4+5 
sqrt(M^ 2+pt5^ 2)  ${M}^{2}+{\overrightarrow{{p}_{T}}}_{5}^{2}$  $M=$ mass of particle 3+4 
sqrt(M^ 2+ptj1^ 2)  ${M}^{2}+{\overrightarrow{{p}_{T}}}_{{j}_{1}}^{2}$  $M=$ mass(3+4), ${j}_{1}=$ leading ${p}_{T}$ jet 
pt(photon)  ${\overrightarrow{{p}_{T}}}_{\gamma}^{2}$  
pt(j1)  ${\overrightarrow{{p}_{T}}}_{{j}_{1}}^{2}$  
HT  ${\sum}_{i=1}^{n}{{p}_{T}}_{i}$  $n$ particles (partons, not jets) 
6.1.7 Masses
Section masses  Description


Higgs pole mass 

Topquark pole mass 

Bottomquark pole mass 

Charmquark pole mass 

Wboson pole mass 

Zboson pole mass 




6.1.8 Basic jets
Section basicjets  Description
 

This logical parameter chooses whether the calculated crosssection should be inclusive in the number of jets found at NLO. An exclusive crosssection contains the same number of jets at nexttoleading order as at leading order. An inclusive crosssection may instead contain an extra jet at NLO. 


This specifies the jetfinding algorithm that is used, and can take the values ktal (for the Run II ${k}_{T}$algorithm), ankt (for the “anti${k}_{T}$” algorithm [59]), cone (for a midpoint cone algorithm), hqrk (for a simplified cone algorithm designed for heavy quark processes) and none (to specify no jet clustering at all). The latter option is only a sensible choice when the leading order crosssection is welldefined without any jet definition: e.g. the single top process, $q{\overline{q}}^{\prime}\to t\overline{b}$, which is finite as ${p}_{T}(\overline{b})\to 0$. 


These specify the values of ${p}_{T,\mathrm{min}}$ and $\eta {\text{}}_{\mathrm{max}}$ for the jets that are found by the algorithm. 


Optional parameter for setting a minimum jet rapidity $\eta {\text{}}_{\mathrm{min}}$. 


Optional parameter for setting maximum jet ${p}_{T,\mathrm{max}}$ 


If the final state of the chosen process contains either quarks or gluons then for each event an attempt will be made to form them into jets. For this it is necessary to define the jet separation
so that after jet combination, all jet pairs are separated by $\mathrm{\Delta}R>$ Rcutjet. 


Optional parameter for using jet rapidity rather than pseudorapidity when performing jet cuts. Default is .true.. 





6.1.9 Mass cuts
Section masscuts  Description

m34min, m34max, m56min, m56max, m3456min, m3456max 
These parameters represent a basic set of mass cuts that are be applied to the calculated crosssection. The only events that contribute to the crosssection will have, for example, m34min $<$ m34$<$ m34max where m34 is the invariant mass of particles 3 and 4 that are specified by nproc. m34min$>0$ is obligatory for processes which can involve a virtual photon, such as nproc=31. By default, the maximum settings are set to $\sqrt{s}$. 




6.1.10 Cuts
Section cuts  Description


If this parameter is set to .false., then no additional cuts are applied to the events (apart from those in Sections 6.1.8 and 6.1.9) and the remaining parameters in this section are ignored. Otherwise, events will be rejected according to a set of cuts that is specified below. Further options may be implemented by editing src/User/gencuts_user.f90. 
ptleptmin, etaleptmax 
These specify the values of
${p}_{T,\mathrm{min}}$
and
$\eta {\text{}}_{\mathrm{max}}$
for one of the leptons produced in the process. One can also introduce optional
settings 
etaleptveto 
This should be specified as a pair of double precision numbers that indicate a rapidity range that should be excluded for the lepton that passes the above cuts. 
ptminmiss 
Specifies the minimum missing transverse momentum (coming from neutrinos). 
ptlept2min,

These specify the values of
${p}_{T,\mathrm{min}}$
and
$\eta {\text{}}_{\mathrm{max}}$
for the remaining leptons in the process. This allows for staggered cuts where,
for instance, only one lepton is required to be hard and central. One can also
introduce optional settings 
etalept2veto 
This should be specified as a pair of double precision numbers (separated by a space) that indicate a rapidity range that should be excluded for the remaining leptons. 




6.1.11 Cuts (continued)
Section cuts  Description
 
m34transmin 
For general processes, this specifies the minimum transverse mass of particles 3 and 4,
For the $W(\to \mathit{\ell \nu})\gamma $ process the role of this cut changes, to become instead a cut on the transverse cluster mass of the $(\mathit{\ell \gamma},\nu )$ system,
For the $\mathit{Z\gamma}$ process this parameter specifies a simple invariant mass cut,
A final mode of operation applies to the
$\mathit{W\gamma}$
process and is triggered by a negative value of
In each case the screen output indicates the cut that is applied. 

Rjlmin 
Using the definition of $\mathrm{\Delta}R$ given above in Eq.6.8), requires that all jetlepton pairs are separated by $\mathrm{\Delta}R>$ R(jet,lept)_min. 

Rllmin 
When nonzero, all leptonlepton pairs must be separated by $\mathrm{\Delta}R>$ R(lept,lept)_min. 

delyjjmin 
This enforces a rapidity gap between the two hardest jets ${j}_{1}$ and ${j}_{2}$, so that: ${\eta}^{{j}_{1}}{\eta}^{{j}_{2}}>$ delyjjmin. 

jetsopphem 
If this parameter is set to .true., then the two hardest jets are required to lie in opposite hemispheres, ${\eta}^{{j}_{1}}\cdot {\eta}^{{j}_{2}}<0$. 

lbjscheme 
This integer parameter provides no additional cuts when it takes the value 0. When equal to 1 or 2, leptons are required to lie between the two hardest jets. With the ordering ${\eta}^{{j}_{\text{\u2212}}}<{\eta}^{{j}_{\text{+}}}$ for the rapidities of jets ${j}_{1}$ and ${j}_{2}$: lbjscheme = 1 : ${\eta}^{{j}_{\text{\u2212}}}<{\eta}^{\mathrm{leptons}}<{\eta}^{{j}_{\text{+}}}$; lbjscheme = 2 : ${\eta}^{{j}_{\text{\u2212}}}+$ Rcutjet$<{\eta}^{\mathrm{leptons}}<{\eta}^{{j}_{\text{+}}}$Rcutjet. 

ptbjetmin, etabjetmax 
If a process involving $b$quarks
is being calculated, then these can be used to specify stricter values of
${p}_{T}^{\mathrm{min}}$
and $\eta {\text{}}^{\mathrm{max}}$
for $b$jets.
Similarly, values for 





6.1.12 Photon
Note that all the photon cuts specified in this section of the input file, are applied even if makecuts is set to .false..
Section photon  Description
 
fragmentation 
This parameter is a logical variable that determines whether the production of photons by a parton fragmentation process is included. If fragmentation is set to .true., the code uses a standard cone isolation procedure (that includes LO fragmentation contributions in the NLO calculation). If fragmentation is set to .false., the code implements a Frixionestyle photon cut [55],
In this equation, ${R}_{0}$, ${\mathit{\epsilon}}_{h}$ and $n$ are defined by cone_ang, epsilon_h and n_pow respectively (see below). ${E}_{T,i}^{j}$ is the transverse energy of a parton, ${E}_{T}^{\gamma}$ is the transverse energy of the photon and $R$ is the separation between the photon and the parton using the usual definition
$n$ is an integer parameter which by default is set to 1. 

fragmentation_set 
A length eight character variable that is used to choose the particular photon fragmentation set. Currently implemented fragmentation functions can be called with ‘BFGSet_I’, ‘BFGSetII’ [60] or ‘GdRG__LO’ [61]. 

fragmentation_scale 
A double precision variable that will be used to choose the scale at which the photon fragmentation is evaluated. 

gammptmin 
This specifies the value of
${p}_{T}^{\mathrm{min}}$
for the photon with the largest transverse momentum. Note that this cut,
together with all the photon cuts specified in this section of the input file, are
applied even if makecuts is set to .false.. One can also add an entry for


gammrapmax 
This specifies the value of
$y{\text{}}^{\mathrm{max}}$
for any photons produced in the process. One can also add an entry for


gammpt2, gammpt3 
The values of ${p}_{T}^{\mathrm{min}}$ for the second and third photons, ordered by ${p}_{T}$. 

Rgalmin 
Using the usual definition of $R$ above, this requires that all photonlepton pairs are separated by $R>$ Rgalmin. This parameter must be nonzero for processes in which photon radiation from leptons is included. 

Rgagamin 
Using the usual definition of $R$ above, this requires that all photon pairs are separated by $R>$ Rgagamin. 

Rgajetmin 
Using the usual definition of $R$ above, this requires that all photonjet pairs are separated by $R>$ Rgajetmin. 

cone_ang 
Fixes the cone size (${R}_{0}$) for photon isolation. This cone is used in both forms of isolation. 

epsilon_h 
This cut controls the amount of radiation allowed in cone when fragmentation is set to .true.. If epsilon_h $<1$ then the photon is isolated using ${\sum}_{\in {R}_{0}}{E}_{T}(\mathrm{had})<{\mathit{\epsilon}}_{h}\phantom{\rule{0.17em}{0ex}}{p}_{T}^{\gamma}.$ Otherwise epsilon_h$>1$ sets ${E}_{T}(\mathit{max})$ in ${\sum}_{\in {R}_{0}}{E}_{T}(\mathrm{had})<{E}_{T}(\mathit{max})$. 

n_pow 
When using the Frixione isolation prescription, the exponent $n$ in Eq. (6.13). 

fixed_coneenergy 
This is only operational when using the Frixione isolation prescription. If fixed_coneenergy is .false. then ${\mathit{\epsilon}}_{h}$ controls the amount of hadronic energy allowed inside the cone using the Frixione isolation prescription (see above, Eq. (6.13)) If fixed_coneenergy is .true. then this formula is replaced by one where ${\mathit{\epsilon}}_{h}{E}_{T}^{\gamma}\to {\mathit{\epsilon}}_{h}$. 

hybrid, R_inner 
If hybrid is set to .true. use a hybrid isolation scheme with Frixione isolation on an inner cone of radius R_inner. 





6.1.13 Histograms
Section histogram  Description


Write output histograms suitable as input for topdrawer. 

Write output histograms as whitespaceseparated columns. 

Use the new plotting infrastructure introduced in MCFM10.0 




6.1.14 Imtegration
Section integration  Description


When 

Initialization seed for MT19937 pseudo random number generator. 

Relative precision goal for the integration. 

When 

When 

Sets the relative precision goal for the warmup run. Unless this precision is reached, the number of calls for the warmup is increased. 

Sets the ${\chi}^{2}$ per iteration goal for the warmup run. Unless the ${\chi}^{2}\u2215\text{it.}$ of the warmup is below this target, the number of calls for the warmup is increased. 




6.2 Process specific options
6.2.1 Single Top
Section singletop  Description


Sets real Wilson coefficient of $\mathcal{\mathcal{Q}}(3,33)\mathit{\phi q}$ for processes 164 and 169. See 13.40 and ref. [33]. 

Sets real and imaginary part of the $\mathcal{\mathcal{Q}}33\mathit{\phi ud}$ Wilson coefficient. 

Sets real and imaginary part of the $\mathcal{\mathcal{Q}}33\mathit{uW}$ Wilson coefficient. 

Sets real and imaginary part of the $\mathcal{\mathcal{Q}}33\mathit{dW}$ Wilson coefficient. 

Sets real and imaginary part of the $\mathcal{\mathcal{Q}}33\mathit{uG}$ Wilson coefficient. 

Sets real and imaginary part of the $\mathcal{\mathcal{Q}}33\mathit{dG}$ Wilson coefficient. 
 

Enable contributions of order
$1\u2215{\Lambda}^{4}$
when set to 

When set to 

When set to 

At NNLO there are several different contributions from vertex
corrections on the lightquark line, heavyquark line in production,
and heavyquark line in the topquark decay. Additionally there
are oneloop times oneloop interference contributions between all
three contributions. For a fully inclusive calculation without decay





6.2.2 Anomalous $W\u2215Z$ couplings
Section anom_wz  Description

enable 
Boolean flag to enable anomalous Wboson and Zboson coupling contributions for certain processes. False has the same effect as setting all anomalous couplings to zero, but additionally skips computation of anomalous coupling code parts. 


delg1_z 
$\mathrm{\Delta}{g}_{1}^{Z}$ See section 13.24. 
delk_z 
$\mathrm{\Delta}{\kappa}^{Z}$ See section 13.24. 
delk_g 
$\mathrm{\Delta}{\kappa}^{\gamma}$ See sections 13.24 and 13.72. 
lambda_z 
${\Lambda}^{Z}$ See section 13.24. 
lambda_g  


h1Z 
${h}_{1}^{Z}$ Anomalous couplings for $\mathit{Z\gamma}$ process at NNLO. See section 13.73. 
h1gam 
${h}_{1}^{\gamma}$ See section 13.73. 
h2Z 
${h}_{2}^{Z}$ See section 13.73. 
h2gam 
${h}_{2}^{\gamma}$ See section 13.73. 
h3Z 
${h}_{3}^{Z}$ See section 13.73. 
h3gam 
${h}_{3}^{\gamma}$ See section 13.73. 
h4Z 
${h}_{4}^{Z}$ See section 13.73. 
h4gam 
${h}_{4}^{\gamma}$ See section 13.73. 


tevscale 
Formfactor scale, in TeV See section 13.24. No formfactors are applied to the anomalous couplings if this value is negative. 




6.2.3 $W\u2215Z$+2 jets
Section wz2jet  Description


This only has an effect when running a $W+2$ jets or $Z+2$ jets process. When .true., it includes the effect of fourquark processes. Please see section 13.10 below. 

This only has an effect when running a $W+2$ jets or $Z+2$ jets process. When .true., it includes the effect of twoquark, twogluon processes. Please see section 13.10 below. 




6.2.4 H jetmass
Section hjetmass  Description

mtex 
Sets the order $k=0,2,4$ of the $1\u2215{m}_{t}^{k}$ expansion for virtual corrections in the $H+$jet process 200. See section 13.43. 




6.2.5 Anomalous $H$ couplings
Section anom_higgs  Description

hwidth_ratio 
For processes 123–126, 128–133 only, this variable provides a rescaling of the width of the Higgs boson. Couplings are rescaled such that the corresponding cross section close to the Higgs boson peak is unchanged. Further details of this procedure are given in arXiv:1311.3589. 

See arXiv:1311.3589. 




6.2.6 Extra
Section extra  Description

debug 
A logical variable which can be used during a debugging phase to mandate special behaviours. Passed by common block common/debug/debug. 
verbose 
A logical variable which can be used during a debugging phase to write special information. Passed in common block common/verbose/verbose. 
new_pspace 
A logical variable which can be used during a debugging phase to test alternative versions of the phase space. Passed in common block common/new_pspace/new_pspace. 
spira 
A logical variable. If spira is 
noglue 
A logical variable. The default value is 
ggonly 
A logical variable. The default value is 
gqonly 
A logical variable. The default value is 
omitgg 
A logical variable. The default value is 
clustering 
This logical parameter determines whether clustering is performed to yield jets.
Only during a debugging phase should this variable be set to 
colourchoice 
If colourchoice=0, all colour structures are included ($W,Z+2$ jets). If colourchoice=1, only the leading colour structure is included ($W,Z+2$ jets). 
rtsmin 
A minimum value of $\sqrt{{s}_{12}}$, which ensures that the invariant mass of the incoming partons can never be less than rtsmin. 

Flag to set the use of the userimplemented reweighting procedure





6.2.7 Dipoles
Section dipoles  Description

aii 
A double precision variable which can be used to limit the kinematic range for the subtraction of initialinitial dipoles as suggested by Trocsanyi and Nagy [63]. The value aii=1 corresponds to standard CataniSeymour subtraction. 
aif 
A double precision variable which can be used to limit the kinematic range for the subtraction of initialfinal dipoles as suggested by Trocsanyi and Nagy [63]. The value afi=1 corresponds to standard CataniSeymour subtraction. 
afi 
A double precision variable which can be used to limit the kinematic range for the subtraction of finalinitial dipoles as suggested by Trocsanyi and Nagy [63]. The value afi=1 corresponds to standard CataniSeymour subtraction. 
aff 
A double precision variable which can be used to limit the kinematic range for the subtraction of finalfinal dipoles as suggested by Trocsanyi and Nagy [63]. The value aff=1 corresponds to standard CataniSeymour subtraction. 
bfi 
A double precision variable which can be used to limit the kinematic range for the subtraction of finalinitial dipoles in the photon fragmentation case. 
bff 
A double precision variable which can be used to limit the kinematic range for the subtraction of finalfinal dipoles in the photon fragmentation case. 



