Tools to scan the mSUGRA phasespace

This document shows how to use two simple macros, susygen.py and susymap.py, to scan any part of the 5D mSUGRA phasespace. These macros can only scan for information calculated by IsaSugra (mostly only the SUSY mass spectrum), as no events are generated. To automatically generate AtlFast SUSY events look at: * Generating AtlFast SUSY Events

IsaSugra

IsaSugra is a subprogram of IsaJet that calculates the mass-spectrum and decay-channels for a number of SUSY models, including mSUGRA. As input parameters, for mSUGRA, it requires the mSUGRA coordinate and the top mass. The output can take different forms of which two are important here. The first is a text file including all SUSY and Higgs masses and decay channels. This file can be used directly in Jimmy to produce events. The second is a more easily read text file also including all masses and decay channels, but in addition also containing some SM values that may be affected by SUSY, such as $Br(B->s\gamma )$, $Br(B->\mu \mu )$, $\Delta a_{\mu }$ and $\Omega _{CDM}h^{2}$ which is the cold dark matter relic density this SUSY point would produce.

If the susy point under consideration is theoretically not allowed (for example because at this point no electroweak symmetry breaking is possible), IsaSugra will give an error, indicating the reason why the susy point was rejected.

At Nikhef the most resent IsaSugra release (7.75) is installed at /project/atlas/users/nicoleru/Store/isajet/isasugra.x

When studying a susypoint with |A0|>0, keep in mind that some theoretical papers use a deffinition of A0 that differs by a minus sign from the one used by programs such as IsaSugra.

Scanning the mSUGRA phasespace

The macro susygen.py runs IsaSugra for a certain susypoint and checks if it is theoretically valid. If it is valid it retrieves certain values from the output files. The information gathered by susygen.py is printed to a text file. At present the following values are written (in this order on one line):

M0, M12, A0, tan(beta), sgn(mu), theory, lightest neutralino mass, lightest Higgs-boson mass, the LSP, BR(b->s+gamma), omega*h^2, Br(Bs->mu+mu) and delta a(mu).

The values are given in GeV (when applicable). The theory value is either 0 (ok) or 1 (not valid) and if 1 no further information is included. The LSP is indicated using the particle numbering of Herwig (lightest neutralino = 450, stau = 429).

susymap.py is a simple program that defines a grid (with evenly separated point in any variable) in the 5D mSUGRA phasespace and starts the text file in which all information will be printed (it also adds a line indicating which variables will be printed). susymap.py will call susygen.py for each point on the grid. The result of this program is a file indicating all relevant information about a large number of susypoints.

Examples of results

As an example of the use of these macros I have run them for a 2D grid, with tan(beta)=50, A0=0 and sgn(mu)=+1 fixed. Both m0 and m12 are varied in 25 steps from 100 to 2500 (so basically 625 point). On stoomboot (submitted in 1 job) this took approximately 3 hours. The output file is located at /project/atlas/users/nicoleru/Store/map_A0_0_tan_50_si_1.txt. The first few lines of that file are the following:

   m0   m12    A0   tan  sign  EWSB   neutalino   Higgsmass  LSP  BR(b->s+gamma)      omega*h^2    Br(Bs->mu+mu)    dela a(mu)
  100   100     0    50     1     1
  100   200     0    50     1     1
  100   300     0    50     1     1
  100   400     0    50     1     1
  100   500     0    50     1     1
  100   600     0    50     1     1
  100   700     0    50     1     0  292.791500  116.384600  429    2.747800e-04   1.000000e+04    1.181600e-08   2.434000e-09
  100   800     0    50     1     0  336.540300  116.977000  429    2.909800e-04   1.000000e+04    9.644800e-09   1.876200e-09
  100   900     0    50     1     0  380.830800  117.482900  429    3.016500e-04   1.000000e+04    8.222600e-09   1.488700e-09
  100  1000     0    50     1     0  424.931700  117.990400  429    3.092800e-04   1.000000e+04    7.275100e-09   1.209400e-09
  100  1100     0    50     1     0  469.779800  118.380700  429    3.149100e-04   1.000000e+04    6.617900e-09   1.001000e-09
  100  1200     0    50     1     0  514.302500  118.781100  429    3.194900e-04   1.000000e+04    6.133900e-09   8.418700e-10
  100  1300     0    50     1     0  559.054500  119.079700  429    3.213500e-04   1.000000e+04    5.769000e-09   7.175600e-10

Figure 1: The m0-m12 plane for tan(beta)=50, sgn(mu)=+1 and A0=0. The black area indicated the theoretically excluded region. In the shaded area the stau is the LSP. The colored lines show the excuded regions.

Using a root file to put all this information in a plot leads to figure 1. (/project/atlas/users/nicoleru/Store/rootfilefigure1.C)

Of course this program can also be used to plot a A0-tan(beta) plain, or to see how certain variables change close to the point you are studying.

As an example of other possible use of this program, I have extended the grid used in the first example to include different values of A0 (and changed tan(beta) to 10), namely A0={-3,-2,-1,0,1,2,3}*m12. Due to WMAP information, most of the m0-m12 plane is excluded. Only two very thin lines remain, one that borders the LSP=stau area, and one line at very low m12. Due to this, we can basically reduce the 5D mSUGRA phasespace to 4D, by relating m0 to m12. Using another root file, it is possible to plot the (lightest) Higgs-mass, or Br(B->s+gamma) or any other variable, as function of m12, for fixed A0, tan(beta) and sgn(mu). (/project/atlas/users/nicoleru/Store/rootfilefigure2.C)

Figure 2: Values on the WMAP allowed line (as function of m12) for tan(beta)=10, sgn(mu)=1 and different values of A0.

How to use

Both macros (susygen.py and susymap.py) should be easy to use. Just copy them to your own directory. They are both located at: /project/atlas/users/nicoleru/Store/.

As long as you are only interested in the information currently optained by susygen.py, you only need to change susymap.py to form your own 5D susy grid. If you need additional information, the mass of one of the susy particles for example, it will take slightly more work, as you need to know exactly where that information is printed in the IsaSugra output files. If you have any problems getting this additional information, just ask me (Nicole).

Making the plots is more difficult as I don't have one root file that, when you feed it a susymap file, simply gives you a plot similar to figure 1. I will work on that. In the mean time, the root file used to make figure 1 should easily be convertible to other grids.