MappingmSUGRA

This document shows how to use two simple macros, susygen.py and susymap.py, to scan any part of the 5D mSUGRA phasespace. For now these macros can only scan for information calculated by IsaSugra (mostly only the SUSY mass spectrum), but they will soon be augmented with the possibility of producing AtlFast events and thereby scan the phasespace for information like cross-section, the fraction of events expected to have N leptons, ext.

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 requiers 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 easaly 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 dencity 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.

Scanning the mSUGRA phasespace

The macro susygen.py runs IsaSugra for a certain susypoint and cheacks 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 prensent 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 dela 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.

susymap.py is a simple program that defines a grid (with evenly seperated point in any variable) in the 5D mSUGRA phasespace and starts the textfile 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 macro's I have run them for a 2 dimentional grid, with tan(beta)=50, A0=0 and sgn(mu)=+1 fixed. Both m0 and m12 are varied in 25 steps form 100 to 2500 (so basicly 625 point). On stoomboot (submitted in 1 job) this took approximatly 3 hours. The output file can be found here. Using a root file to put all this information in a plot leads to figure 1. Similar plots for a large variaty of tan(beta), A0 and sgn(mu) can be found here.

Of course this program can aslo be used to plot a A0-tan(beta) plain, or to see how certain variables change in close proximation 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 diffent values of A0, namely A0={-3,-2,-1,0,1,2,3}. 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 basicly 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).

How to use

The macro's should be easy to use. Just copy them to the directory you want to use. As long as you are only interseted in the information currently cheacked by susygen.py, you only need to change susymap.py to form you own 5D susygrid. If you need additional information, the mass of one of the particles for example, it will take slightly more work, as you need to know exactly where that information is printed in de 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 simply gives you a plot similar to figure 1, whene you feed it an susymap file. I will work on that. In the mean time, the root file used to make figure 1 should easily be convertable to other grids.