Given an input model, the first task of SModelS is to decompose the full model into a sum of elements. Based on the input format, which can be
as explained in the Basic Input section, two distinct (but similar) decomposition methods are applied: the SLHA-based or the LHE-based decomposition.
The SLHA file describing the input model is required to contain the masses of all the BSM states as well as their production cross-sections and decay branching ratios. All the above information must follow the guidelines of the SLHA format. In particular, the cross-sections also have to be included as SLHA blocks according to the SLHA cross-section format.
Once the production cross-sections are read from the input file, all the cross-sections for production of two Z2-odd states are stored and serve as the initial step for the decomposition. (All the other cross-sections with a different number of Z2-odd states are ignored.) Starting from these primary mothers, all the possible decays are generated according to the information contained in the DECAY blocks. This procedure is represented in the figure below:
Each of the possible cascade decays for each mother corresponds to a branch. In order to finally generate elements, all the branches are combined in pairs according to the production cross-sections, as shown below:
For instance, assume [b1,b2,b3] and [B1,B2] represent all possible branches (or cascade decays) for the primary mothers A and B, respectively. Then, if a production cross-section for \(pp \rightarrow A+B\) is given in the input file, the following elements will be generated:
[b1,B1], [b1,B2], [b2,B1], [b2,B2], [b3,B1] and [b3,B2]
Each of the elements generated according to the procedure just described will also store its weight, which equals its production cross-section times all the branching ratios appearing in it. In order to avoid a too large number of elements, only those satisfying a minimum weight requirement are kept. Furthermore, the elements are grouped according to their topologies. The final output of the SLHA decomposition is a list of such topologies, where each topology contains a list of the elements generated during the decomposition.
Some models may contain a large number of new states and each may have a large number of possible decays. As a result, long cascade decays are possible and the number of elements generated by the decomposition process may become too large, and the computing time too long. For most practical purposes, however, elements with extremely small weights (cross-section times BRs) can be discarded, since they will fall well below the experimental limits. Therefore, during the SLHA decomposition, whenever an element is generated with a weight below some minimum value, this element (and all elements derived from it) is ignored. The minimum weight to be considered is given by the sigcut parameter and is easily adjustable (see slhaDecomposer.decompose)
Note that, when computing the theory predictions, the weight of several elements can be combined together. Hence it is recommended to set the value of sigcut approximately one order of magnitude below the minimum signal cross-sections the experimental data can constrain.
More general models can be input through an LHE event file containing parton-level events, including the production of the primary mothers and their cascade decays. Each event can then be directly mapped to an element with the element weight corresponding to the event weight. Finally, identical elements can be combined together (adding their weights). The procedure is represented in the example below:
Notice that, for the LHE decomposition, the elements generated are restricted to the events in the input file. Hence, the uncertainties on the elements weights (and which elements are actually generated by the model) are fully dependent on the Monte Carlo statistics used to generate the LHE file. Also, when generating the events it is important to ensure that the no mass smearing is applied, so the events always contain the mass value for a given particle.
During the decomposition process it is possible to perform several simplifications on the elements generated. In both the LHE and SLHA-based decompositions, two useful simplifications are possible: Mass Compression and Invisible Compression. The main advantage of performing these compressions is that the simplified element is always shorter (has fewer cascade decay steps), which makes it more likely to be constrained by experimental analyses. The details behind the compression methods are as follows:
In case of small mass differences, the decay of an intermediate state to a nearly degenerate one will in most cases produce soft final states, which can not be experimentally detected. Consequently, it is a good approximation to neglect the soft final states and compress the respective decay, as shown below:
After the compression, only the lightest of the two near-degenerate masses are kept in the element, as shown above. The main parameter which controls the compression is minmassgap, which corresponds to the maximum value of \(\epsilon\) in the figure above to which the compression is performed:
Note that the compression is an approximation since the final states, depending on the boost of the parent state, may not always be soft. It is recommended to choose values of minmassgap between 1-10 GeV; the default value is 5 GeV.
and can be easily turned on/off by the flag doCompress in the SLHA or LHE decompositions
Another class of compression is possible when the final states in the last decay as well as the last intermediate state are both invisible experimentally. The most common example is
as the last step of the decay chain, where \(B\) is an insivible particle leading to a MET signature. Since both the neutrino and \(B\) are invisible, for all experimental purposes the effective MET object is \(B + \nu = A\). Hence it is possible to omit the last step in the cascade decay, resulting in a compressed element. Note that this compression can be applied consecutively to several steps of the cascade decay if all of them contain only invisible final states:
and can be easily turned on/off by the flag doInvisible in the SLHA or LHE decompositions