# Integrators¶

Herwig 7 offer a new infrastructure for Monte Carlo integration, which has been developed along with the Matchbox module; a GeneralSampler class is responsible of exploiting a certain algorithm implemented in BinSampler classes to perform the integration of individual subprocesses. Currently, we offer three main algorithms: CellGridSampler, MonacoSampler and FlatBinSampler for validation purposes. All sampling objects are loacted in the /Herwig/Smaplers directory of the default repository.

## General options¶

We report some of the less-expert settings here; a full and detailed documentation can be found in the manual.

Switch to verbose mode to report results on individual sub-process contributions:

set Sampler:Verbose On


Minimum selection probability for a sub-process if its estimated cross section would give rise to a smaller probability:

set Sampler:MinSelection 0.01


Sample sub-processes flat (validation purposes only):

set Sampler:FlatSubprocesses On


Update grid information from statistics accumulated while running:

set Sampler:WriteGridsOnFinish Yes


## Common sampler options¶

Set the number of initial points to be used for estimating sub-process cross sections:

set Sampler:BinSampler:InitialPoints 10000


The following options make sense for the CellGridSampler and MonacoSampler:

Set the number of iterations to be performed (not applicable to flat sampling):

set Sampler:BinSampler:NIterations 3


Set the enhancement factor for the number of points to be used in the next iteration:

set Sampler:BinSampler:EnhancementFactor 1.2


Set a path for data on random number population versus cross section:

set Sampler:BinSampler:RandomNumbers ./population-data


The following options make sense for the CellGridSampler and FlatBinSampler:

Enable the remapping of phase space dimensions using a given number of points:

set Sampler:BinSampler:RemapperPoints 10000


Remap the dimension used for channel selection (only applicable for TreePhasespace):

set Sampler:BinSampler:RemapChannelDimension Yes


Remap the random numbers corresponding to the parton luminosity using a certain number of bins:

set Sampler:BinSampler:LuminsoityMapperBins 20


Remap all other dimensions using a certain number of bins:

set Sampler:BinSampler:GeneralMapperBins 20


Enforce a minimum selection probability for remapper bins even if the estimated cross section would have given rise to a smaller probability:

set Sampler:BinSampler:RemapperMinSelection 0.01


## CellGrid sampler¶

The CellGridSampler is extending the algorithm initially described in 1108.6182. It can be enabled using:

set Sampler:BinSampler CellGridSampler


Options accepted by the CellGridSampler include:

Set the number of exploration points to presample a cell:

set CellGridSampler:ExplorationPoints 1000


Set the number of exploration steps to be performed:

set CellGridSampler:ExplorationSteps 3


Set the gain parameter of the cell splitting algorithm:

set CellGridSampler:Gain 0.3


Set the epsilon parameter of the cell splitting algorithm:

set CellGridSampler:Epsilon 1.0


Enforce a minimum selection probability for cells even if the estimated cross section would have given rise to a smaller probability:

set CellGridSampler:MinimumSelection 0.01


## Monaco sampler¶

The MonacoSampler is implementing the algorithm described in 0811.4559. It can be enabled using:

set Sampler:BinSampler MonacoSampler


Options accepted by the MonacoSampler include:

Set the grid modification rate:

set MonacoSampler:Alpha 0.8


Set the number of grid divisions:

set MonacoSampler:GridDivisions 48


## Flat sampling¶

Flat sampling is offered as a cross check facitlity. It can be enabled via:

set Sampler:BinSampler FlatBinSampler


There are no parameters for flat sampling.

# NLO consistency checks¶

## Subtraction checks¶

The convergence of the unintegrated subtraction can be checked by indicating a directory where the relevant data should be stored:

set Factory:SubtractionData ./subtraction-checks


Keep in mind that phase space sampling will be highly restricted to the singular regions, so no physics analysis should be run when enabling this test.

Several types of plots are supported:

Perform enveloping of the ratio $$D/M$$:

set Factory:SubtractionPlotType LinRatio


Perform enveloping of the logarithm of the relative difference $$|D-M|/M$$:

set Factory:SubtractionPlotType LogRelDiff


## Pole cancellation¶

Similarly to the subtraction checks, histograms of pole cancellation measures can be performed by indicating a directory to store the data:

set Fatctory:PoleData ./pole-checks


## List of dipoles and mappings¶

The list of dipoles which have been found for a certain process and the mapping of real emission to Born legs can be accessed via:

set Factory:InitVerbose On