MSEL = 19
In PYTHIA, the particle content of a two-Higgs-doublet scenario is included: two neutral scalar particles, 25 and 35, one pseudoscalar one, 36, and a charged doublet, . (Of course, these particles may also be associated with corresponding Higgs states in larger multiplets.) By convention, we choose to call the lighter scalar Higgs and the heavier . The pseudoscalar is called and the charged . Charged-Higgs production is covered in section .
A number of processes have been duplicated for and . The correspondence between ISUB numbers is shown in the table above: the first column of ISUB numbers corresponds to , the second to , and the third to . Note that several of these processes are not expected to take place at all, owing to vanishing Born term couplings. We have still included them for flexibility in simulating arbitrary couplings at the Born or loop level, or for the case of mixing between the scalar and pseudoscalar sectors.
A few Standard Model Higgs processes have no correspondence in the scheme above. These include
The processes 111-113, 183-185 and 188-190 have only been worked out in full detail for the Standard Model Higgs case, and not when e.g. squark loop contributions need be considered. The approximate procedure outlined in subsection , based on combining the kinematics shape from simple expressions in the limit with a normalization derived from the cross section, should therefore be viewed as a first ansatz only. In particular, it is not recommended to try the non-default MSTP(38)=0 option, which is incorrect beyond the Standard Model.
In processes 121, 122, 181, 182, 186 and 187 the recoiling heavy flavour is assumed to be top, which is the only one of interest in the Standard Model, and the one where the parton-distribution-function approach invoked in processes 3, 151 and 156 is least reliable. However, it is possible to change the quark flavour in 121 etc.; for each process ISUB this flavour is given by KFPR(ISUB,2). This may become relevant if couplings to states are enhanced, e.g. if in the MSSM. The matrix elements in this group are based on scalar Higgs couplings; differences for a pseudoscalar Higgs remains to be worked out.
By default, the has the couplings of the Standard Model Higgs, while the and have couplings set in PARU(171) - PARU(178) and PARU(181) - PARU(190), respectively. The default values for the and have no deep physics motivation, but are set just so that the program will not crash due to the absence of any couplings whatsoever. You should therefore set the above couplings to your desired values if you want to simulate either or . Also the couplings of the particle can be modified, in PARU(161) - PARU(165), provided that MSTP(4) is set to 1.
For MSTP(4)=2, the mass of the (in PMAS(25,1)) and the value (in PARU(141)) are used to derive the masses of the other Higgses, as well as all Higgs couplings. PMAS(35,1) - PMAS(37,1) and PARU(161) - PARU(195) are overwritten accordingly. The relations used are the ones of the Born-level MSSM [Gun90]. Loop corrections to those expressions have been calculated within specific supersymmetric scenarios, and are known to have a non-negligible effects on the resulting phenomenology. By switching on supersymmetry simulation and setting parameters appropriately, one will gain access to these mass formulae, see section .
Note that not all combinations of
allowed; for MSTP(4)=2 the requirement of a finite mass
imposes the constraint
A more realistic approach to the Higgs mass spectrum is to include radiative corrections to the Higgs potential. Such a machinery has never been implemented as such in PYTHIA, but appears as part of the Supersymmetry framework described in subsection . At tree level, the minimal set of inputs would be IMSS(1)=1 to switch on SUSY, RMSS(5) to set the value (this overwrites the PARU(141) value when SUSY is switched on) and RMSS(19) to set mass. However, the significant radiative corrections depend on the properties of all particles that couple to the Higgs boson, and the user may want to change the default values of the relevant RMSS inputs. In practice, the most important are those related indirectly to the physical masses of the third generation supersymmetric quarks and the Higgsino: RMSS(10) to set the left-handed doublet SUSY mass parameter, RMSS(11) to set the right stop mass parameter, RMSS(12) to set the right sbottom mass parameter, RMSS(4) to set the Higgsino mass and a portion of the squark mixing, and RMSS(16) and RMSS(17) to set the stop and bottom trilinear couplings, respectively, which specifies the remainder of the squark mixing. From these inputs, the Higgs masses and couplings would be derived. Note that switching on SUSY also implies that Supersymmetric decays of the Higgs particles become possible if kinematically allowed. If you do not want this to happen, you may want to increase the SUSY mass parameters. (Use CALL PYSTAT(2) after initialization to see the list of branching ratios.)
Pair production of Higgs states may be a relevant source, see section below.
Finally, heavier Higgses may decay into lighter ones, if kinematically allowed, in processes like or . Such modes are included as part of the general mixture of decay channels, but they can be enhanced if the uninteresting channels are switched off.