Minimum-bias and underlying-event tunes

For a detailled discussion of the models used in this tunes, we refer to Eur.Phys.J. C72 (2012) 2225. That paper also explains the procedure used to deduce the tunes.

H7-UE-MMHT (default tune)

It was shown in Ref.[JHEP 1310 (2013) 113] that a good description of both underlying event and double parton scattering data Ref.[JHEP 1303 (2013) 129] can be obtained if one includes the latter in the data being fit to with a sufficiently high weight. We followed the procedure described in Ref.[JHEP 1310 (2013) 113] using MMHT2014 LO parton distribution function [1] and obtained a tune consistent with double parton scattering data (\(\sigma_{\!\textit{eff}}\approx15\,\mathrm{mb}\)) which also gives a good description of the underlying event data from Tevatron’s lowest energy point \(\sqrt{s}=300\,\mathrm{GeV}\) to the LHC’s \(\sqrt{s}=7\,\mathrm{TeV}\).

Parameter H7-UE-MMHT
\(p_{\perp,0}^{\min}\) 4.39
\(b\) 0.366
\(\mu^2/{\rm GeV}^2\) 2.30
\(p_{\rm disrupt}\) 0.798
\(p_{\rm reco}\) 0.4276

Hw7.0 is released together with the tune H7-UE-MMHT, by default.

Please note that the tune is not suitable for lower energy runs (below \(\sqrt{s}=300\,\mathrm{GeV}\)). Please contact the authors to get more information.

Tunes suitable for older version of Herwig++

The improvements to both shower modules, as well as the inclusion of next-to-leading order cross sections in Herwig 7 have required a new tune to \(e^+ e^−\) data. Therefore, MPI tunes suitable for older version of Herwig can not be used with Herwig 7.

Herwig++ 2.7

Parameter UE-EE-5-CTEQ6L1 UE-EE-5 (MRST LO✶✶)
\(p_{\perp,0}^{\min}\) 3.91 4.620
\(b\) 0.33 0.314
\(\mu^2/{\rm GeV}^2\) 2.30 2.240
\(p_{\rm disrupt}\) 0.80 0.860
\(p_{\rm reco}\) 0.49 0.420

Herwig++ 2.6

From Herwig++ 2.6 we use an automatic energy evolution of the \(p_\perp^{\min}\) governed by the equation:

\(p_\perp^{\min}=p_{\perp,0}^{\min}\left(\frac{\sqrt{s}}{E_0}\right)^b\)

In the table below we present the parameters of the energy-extrapolating underlying-event tunes. Two different tunes for the plain colour reconnection (PCR) model using two different PDF sets MRST LO✶✶ (UE-EE-4) and CTEQ6L1 (UE-EE-4-CTEQ6L1) and one tune for the statistical colour reconnection (SCR) using CTEQ6L1 PDF set, UE-EE-SCR-CTEQ6L1. The first two parameters in the table below describe the running of \(p_\perp^{\min}\) according to the equation above.

Parameter UE-EE-4-CTEQ6L1 UE-EE-4 (MRST LO✶✶) UE-EE-SCR-CTEQ6L1
\(p_{\perp,0}^{\min}\) 2.81 3.11 2.64
\(b\) 0.24 0.21 0.21
\(\mu^2/{\rm GeV}^2\) 1.35 1.11 1.5
\(p_{\rm disrupt}\) 0.75 0.80 0.8
\(p_{\rm reco}\) 0.61 0.54 N/A
\(c\) N/A N/A 0.01
\(f\) N/A N/A 0.21
\(N_{\rm steps}\) N/A N/A 10
\(\alpha\) N/A N/A 0.66

These tunes are implemented in the following input files:

N.B.

  • The reference scale] \(E_0\) is set to 7000 GeV.
  • Please note that LHC-UE-EE-4 is identical to LHC-UE-EE-3, but uses the value for \(p_{\perp}^{\min}\) from the extrapolation formula also for 900, 1800 and 7000 GeV. There are no significant differences in the observables. The same is true for LHC-UE-EE-4-CTEQ6L1 and LHC-UE-EE-3-CTEQ6L1.

Herwig++ 2.5

UE-EE-3 for MRST LO✶✶ (default in Herwig++ 2.5) and CTEQ6L1 PDFs

The third generation of energy-extrapolation tunes for two PDF sets can be obtained here:

The tunes provide a better energy-extrapolation and can be used at 14 TeV.

Plots for various analyses are provided Herwig++ 2.5 EE3.

PDF CTEQ6L1 MRST LO✶✶ Plots UE-EE-2 vs UE-EE-3 (CTEQ6L1 and MRST LO✶✶)
\(\sqrt{s}\) \(p_\perp^{\min}\)  
900 1.55 1.86 ATLAS_2010_S8894728
1800 2.26 2.55 CDF_2001_S4751469
2760 2.33 2.62 No data
7000 2.752 3.06 ATLAS_2010_S8894728 vs UE-EE-2, vs UE7-2
8000 2.85 3.21 No data
14000 3.16 3.53 No data

with these fixed parameters

Parameter CTEQ6L1 MRST LO✶✶
\(p_{\rm disrupt}\) 0.75 0.80
\(\mu^2/{\rm GeV}^2\) 1.35 1.11
\(p_{\rm reco}\) 0.61 0.54

N.B You have to choose the desired center-of-mass energy in the input files! Please read the instructions inside the files carefully.

UE-EE-2

\(\mu^2 = 1.1/{\rm GeV}^2\), \(p_{\rm disrupt}=0.15\), \(p_{\rm reco}= 0.55\), (input file: LHC-UE-EE-2.in )
\(\sqrt{s}\) Name \(p_\perp^{\min}\) Plots UE-EE-2 vs UE-EE-1 vs (UE7 and MU900)
900 UE-EE-2 2.34 ATLAS_2010_S8894728 900
1800 UE-EE-2 3.09 plots
7000 UE-EE-2 4.02 ATLAS_2010_S8894728 7000

Based on this results we made interpolation to \(\sqrt(s) = 2.76 {\rm TeV}\) and suggest to use ‘’‘ptMin = 3.31’‘’ for UE studies at this energy (input file: LHC-UE-EE-2.in).

N.B. You have to choose the desired center-of-mass energy at two places in the LHC-UE-EE-2.in file! Please read the instructions inside the file carefully.

Minimum Bias and Underlying Event, 900 GeV and 7 TeV

Data Name Main parameters Input file Plots
7000 ATLAS UE LHC-UE7-2 \(p_{\perp}^{\min}=3.36\) \(\mu^2 =0.81 /{\rm GeV}^2\), \(p_{\rm disrupt}=0.35\), \(p_{\rm reco}=0.616\) LHC-UE7-2.in LHC-MU900-2.in
900 ATLAS UE&MB(Nch >=6) MU900-2 \(p_{\perp}^{\min}=2.60\) \(\mu^2 =1.15 /{\rm GeV}^2\), \(p_{\rm disrupt}=0.298\), \(p_{\rm reco}=0.41\) UE7-2 vs UE7-1, UE plots MU900-2 vs MU900-1, MB and UE plots

Note about deprecated tunes

The first generation of MB and UE tunes (LHC-MU900-1, LHC-UE7-1, LHC-UE-EE-*-1) overestimated the production of strange hadrons and the b-fragmentation function was also not described properly. Therefore, we discourage the use of the first generation of MB and UE tunes for production runs. This problem was fixed in the second generation of tunes, which you can obtain here.

How to use these tunes

To generate minimum-bias events, the provided run cards can be used as is. To generate hard processes with one of the tunes provided above, follow these steps in your run card.

  1. Set the colour reconnection, colour disrupt and MPI parameters as found in the respective input file (e.g. LHC-UE7-1.in). The following example uses the parameters from the LHC-UE7-1 tune

    # Colour reconnection settings
    set /Herwig/Hadronization/ColourReconnector:ColourReconnection Yes
    set /Herwig/Hadronization/ColourReconnector:ReconnectionProbability 0.6165547
    
    # Colour Disrupt settings
    set /Herwig/Partons/RemnantDecayer:colourDisrupt 0.3493643
    
    # inverse hadron radius
    set /Herwig/UnderlyingEvent/MPIHandler:InvRadius 0.81
    
    # Min KT parameter
    set /Herwig/UnderlyingEvent/KtCut:MinKT 3.36
    
    # This should always be 2*MinKT!!
    set  /Herwig/UnderlyingEvent/UECuts:MHatMin 6.72
    
  2. Adopt these MPI settings:

    set /Herwig/UnderlyingEvent/MPIHandler:softInt Yes
    set /Herwig/UnderlyingEvent/MPIHandler:twoComp Yes
    set /Herwig/UnderlyingEvent/MPIHandler:DLmode 3
    
  3. Select the matrix element for the hard subprocess, e.g.

    cd /Herwig/MatrixElements/
    insert SimpleQCD:MatrixElements[0] MEqq2gZ2ff
    
  4. Optional: select kinematic cuts to be applied for the hard subprocess, e.g.

    cd /Herwig/Cuts
    set LeptonKtCut:MinKt 10.0*GeV
    
  5. Specify which of the hard processes (selected above in SimpleQCD:MatrixElements) is identical to the UE one, i.e. QCD dijets.

    # -1 means, none of the matrix elements is QCD dijet production
    set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE -1
    

The main difference of these instructions to the ready-to-use minimum bias run cards is that the minimum-bias-specific settings are dropped:

set JetKtCut:MinKT  0.0*GeV
set QCDCuts:MHatMin 0.0*GeV
set QCDCuts:X1Min   0.055
set QCDCuts:X2Min   0.055
[1]In the near future we also plan to provide tunes using CT14 and NNPDF3.0 parton distribution functions.