VINCIA Merging Settings

## Sector Merging (MESS) -- VINCIA only

The VINCIA sector shower employs its own CKKW-L merging scheme, which differs from the one implemented for the simple showers. The biggest difference is that (for gluon emissions) the VINCIA sector shower only possesses a single history (or branching tree), i.e., any given configuration produced by it can be uniquely traced back to every intermediate state the shower has produced on the way. While for gluon splittings, all possible quark permutations need to be taken into account, the VINCIA sector shower is still "maximally bijective", i.e., it has the lowest possible number of histories. MESS merging can therefore be viewed as a midpoint between the CKKW-L and the CKKW scheme, as only a single, deterministic history needs to be generated, which nevertheless exactly reflects the shower history and Sudakov factors are generated dynamically using trial showers. As such, it is specifically designed for merging with high-multiplicity matrix elements.

VINCIA's merging may be enabled by using the VINCIA sector shower and switching merging on:

    PartonShowers:model = 2
    Vincia:sectorShower = on
    Merging:doMerging = on

In addition, the user should set

    Vincia:kineMapFFsplit = 1

since the inverse kinematic map for other splitting maps are not currently available. We also advise running with

    Check:abortIfVeto = on

such that any errors which occur during merging will be flagged as an aborted event (rather than a zero weight event for normally vetoed events).

We note that, different to merging with the simple showers, `Merging:doMerging` should always be set to `on` to enable merging, no matter which merging scale definition is used.

MESS merging can generally replace CKKW-L merging with the simple showers and is illustrated in the command file `main162mess.cmnd` for `main162.cc` and `main164mess.cmnd` for `main164.cc`. However, note that a few modifications to the command file are needed compared to merging with the default shower, as outlined below.

### Matrix Elements

As for merging with the simple showers, the user has to provide LHE files with pre-generated events with up to `N` additional jets. The maximal number of additional jets relative to the Born that do not arise from resonance decays is specified by setting

    Merging:nJetMax = N

### Process Specification

The hard process should be specified through `Merging:Process`, using the following rules:

- The whole string should be encased in curly brackets, `{ ... }`.
- Specify particles one at a time, separated by whitespace, using Pythia's naming conventions in ParticleData.
- The intial and final state should be separated by `>`.
- Exactly two initial state particles should be specified
- Resonance decays may be specifed also using `>`, but if there is more than one particle in the intermediate final state it should be encased in round brackets, `( ... )`. In principle as many nested decays may be specified as required, with additional brackets for each decay. However, we make no claims on how physically sensible merging in such systems is.
- Resonances must decay to exactly 2 particles.
- In addition to Pythia's naming scheme, the following "multiparticles" are defined:
  - `p=p+=pbar=p-=n=nbar=j={1,2,3,4,5,-1,-2,-3,-4,-5,21}`
  - `q=Q=QUARK={1,2,3,4,5}`
  - `qbar=QBAR=ANTIQUARK={-1,-2,-3,-4,-5}`
  - `LEPTONS={11,-11,13,-13,15,-15}`
  - `l+={-11,-13,-15}`
  - `l-={11,13,15}`
  - `NEUTRINOS={12,-12,14,-14,16,-16}`
  - `nu={12,14,16}`
  - `nubar={-12,-14,-16}`
  - `gammaZ={22,23}`

We note that the option `Merging:Process = guess` is not supported, as we believe that the user should be fully aware and in control of what the hard process in the event is.

#### VBF and HEFT

Following the idea that the user needs to specify the event topology exactly, the following two switches can be used to further refine the process:

`flag` Vincia:MergeVBF (`default = off`)
Experimental switch to enable merging in VBF processes.

`flag` Vincia:MergeHEFT (`default = off`)
Experimental switch to consider HEFT couplings in the history construction.

### Merging Scale

MESS merging presently supports three different merging scale definitions. If only `Merging:doMerging = on` is set, the shower evolution variable (a generalised `Ariadne` pT) is used to define the merging scale. Additionally, the merging scale may be defined in terms of a kT cut by setting

    Merging:doKTMerging = on

and (optionally) specifying `Merging:ktType`, see CKKW-L Merging for details. In both cases, the merging scale should be specified through `Merging:TMS` while for merging with a kT cut, also the D parameter should be specified by `Merging:Dparameter`.

It is also possible to define the merging scale by a set of cuts imposed on the events upon generation by setting `Merging:doCutBasedMerging = on`. The user then has to provide the values of the three cuts: `Merging:pTiMS`, `Merging:QijMS`, and `Merging:dRijMS` representing the minimal transverse momentum pTi of a single jet and the minimal invariant mass Qij and minimal separation ΔRij of every pair of jets, as used to generate the events.

### Merging in Resonance Systems (Unvalidated)

The Vincia sector shower supports dedicated merging in resonance systems. At the current stage, however, this should be considered an experimental feature as it has yet to be validated. The following settings are therefore mainly intended for expert users, for instance in the context of getting started on performing such validations.

`flag` Vincia:MergeInResSystems (`default = off`)
Switch to enable merging of additional jets in resonance decay systems. Currenly handles colour-singlet resonances only.

`mode` Vincia:MergeNJetMaxRes (`default = 0`; `minimum = 0`)
Analogue to `Merging:nJetMax`, to communicate the maximum number of additional jets that can be produced from a given resonance decay system by the matrix-element generator. Only used if `Vincia:MergeInResSystems = on`.

`mode` Vincia:MergeNResSys (`default = 0`; `minimum = 0`)
The number of resonance systems allowed to produce jets if `Vincia:MergeInResSystems = on` and `Vincia:MergeNJetMaxRes > 0`.

Additionally, for simple topologies, resonances can be explicitly inserted into the event record if they have not been written to the event file by the matrix element generator. We do note that the sector merging heavily relies on this information and that attaching resonances to leptons is ambiguous for general processes.

`flag` Vincia:InsertResInMerging (`default = off`)
If set to `on`, Vincia tries to attach resonances to final-state leptons and insert them into the event record before constructing the shower history. If turned on, the hard process must explicitly contain a resonance, e.g. `Merging:Process = { p p > Z0 }`.

### Treatment of QED (and other non-QCD) Corrections

Vincia's sector merging algorithm itself is so far restricted to pure QCD corrections. Thus, it is not possible to include any matrix elements with higher-order QED (or other non-QCD) corrections among the merged samples.

Vincia's QED showers can still be enabled, and will be on by default. While no fixed-order QED corrections can be included in the merging at the moment, QED showers can still be used to dress accepted events with logarithmically enhanced radiation after the merging has taken place. This means that the merging algorithm does not include QED clusterings when constructing shower histories, and QED showers are switched off during trial showers. In particular, QED branchings above the (QCD) merging scale in showers off accepted events will always be allowed, regardless of the jet multiplicity.

Vincia's electroweak showers, on the other hand, are currently not supported for merging and should not be switched on. The reason is that the emission of electroweak bosons alters the definition of the hard process, rendering the number of emissions ambiguous. An example is the process `p p > j j` , which may evolve to `p p > j j W+` when electroweak showers are turned on. From the point of view of the QCD showers, this process may now equally well look like a Drell-Yan process `p p > W+` with two QCD branchings. It is therefore only sensible to allow for electroweak showers in the merging when these are also taken into account in the construction of shower histories. As this is currently not the case, electroweak showers are not supported in the sector merging at the moment.