The SUSY routines and common-block variables are described in section . To illustrate the usage of the switches and parameters, we give six simple examples.
Example 1: Light Stop
The first example is an MSSM model with a light neutralino and a light stop , so that can occur. The input parameters are
IMSS(1) = 1, RMSS(1) = 70., RMSS(2) = 70., RMSS(3) = 225., RMSS(4) = -40.,
RMSS(5) = 1.5, RMSS(6) = 100., RMSS(7) = 125., RMSS(8) = 250.,
RMSS(9) = 250., RMSS(10) = 1500., RMSS(11) = 1500., RMSS(12) = -128.,
RMSS(13) = 100., RMSS(14) = 125., RMSS(15) = 800., RMSS(16) = 800.,
RMSS(17) = 0., and RMSS(19) = 400.0.
The top mass is fixed at 175 GeV, PMAS(6,1) = 175.0. The resulting model has GeV and GeV. IMSS(1) = 1 turns on the MSSM simulation. By default, there are no intrinsic relations between the gaugino masses, so GeV, GeV, and GeV. The pole mass of the gluino is slightly higher than the parameter , and the decay occurs almost 100% of the time.
Example 2: SUSY Les Houches Accord spectrum
The second example shows how to input a spectrum file in the SUSY Les Houches Accord format [Ska03] to PYTHIA. First, you should set IMSS(1) = 11 and open the spectrum file you want to use on some unused Logical Unit Number. Then, set IMSS(21) equal to that number, to tell PYTHIA where to read the spectrum file from. This should be done somewhere in your main program before calling PYINIT. During the call to PYINIT, PYTHIA will read the spectrum file, perform a number of consistency checks and issue warning messages if it finds something it does not understand or which seems inconsistent. E.g. BLOCK GAUGE will normally be present in the spectrum file, but since PYTHIA currently cannot use the information in that block, it will issue a warning that the block will be ignored. In case a decay table is also desired to be read in, the Logical Unit Number on which the decay table is opened should be put in IMSS(22). To avoid inconsistencies, the spectrum and the decay table should normally go together, so IMSS(22) should normally be equal to IMSS(21).
Example 3: Calling ISASUSY 7.71 at runtime using
IMSS(1) = 12
The third example shows how to use the built-in run-time interface to ISASUSY with the IMSS(1) = 12 option. First, the PYTHIA source code needs to be changed. Rename the function VISAJE to, for example, FDUMMY, rename the subroutines SUGRA and SSMSSM to e.g. SDUMM1 and SDUMM2, and recompile. In the calling program, set IMSS(1) = 12 and the RMSS input parameters exactly as in example 5, and compile the executable while linked to both ISAJET and the modified PYTHIA. The resulting mass and mixing spectrum is printed in the PYTHIA output.
Example 4: Calling ISASUSY 7.71 at runtime using
IMSS(1) = 13
The fourth example shows how to use the built-in run-time interface to ISASUSY with the IMSS(1) = 13 option. First, the PYTHIA source code needs to be changed, cf. the previous example. In the calling program, set IMSS(1) = 13 and open an ISAJET SUSY model input file on any available Logical Unit Number. The contents of the file should be exactly identical to what would normally be typed when using the ISAJET RGE executable stand-alone (normally isasugra.x). Then, store that Unit Number in IMSS(20), that will enable PYTHIA to access the correct file during initialization. Compile the executable while linked to both ISAJET and the modified PYTHIA. The resulting mass and mixing spectrum is printed in the PYTHIA output.
Example 5: Approximate SUGRA
This example shows you how to get a (very) approximate SUGRA model. Note that this way of obtaining the SUSY spectrum should never be used for serious studies. The input parameters are
IMSS(1) = 2, RMSS(1) = 200., RMSS(4) = 1., RMSS(5) = 10., RMSS(8) = 800., and RMSS(16) = 0.0.
The resulting model has GeV, GeV, GeV, GeV, GeV, GeV, GeV, GeV, and GeV. It corresponds to the choice =800 GeV, 200 GeV, , , and sign(). The output is similar to an ISASUSY run, but there is not exact agreement.
Example 6: ISASUSY 7.71 Model
The final example demonstrates how to convert the output of an ISASUSY run directly into the PYTHIA format, i.e. if SLHA output is not available. This assumes that you already made an ISASUSY run, e.g. with the equivalents of the input parameters above. From the output of this run you can now extract those physical parameters that need to be handed to PYTHIA, in the above example
IMSS(1) = 1, IMSS(3) = 1, IMSS(8) = 0, IMSS(9) = 1, RMSS(1) = 79.61,
RMSS(2) = 155.51, RMSS(3) = 533.1, RMSS(4) = 241.30, RMSS(5) = 10.,
RMSS(6) = 808.0, RMSS(7) = 802.8, RMSS(8) = 878.4, RMSS(9) = 877.1,
RMSS(10) = 743.81, RMSS(11) = 871.26, RMSS(12) = 569.87, RMSS(13) = 803.20,
RMSS(14) = 794.71, RMSS(15) = -554.96, RMSS(16) = -383.23, RMSS(17) = -126.11,
RMSS(19) = 829.94 and RMSS(22) = 878.5.