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 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 , , and 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/.
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.
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
Using a root file to put all this information in a plot leads to figure 1. Similar plots for a large variety of tan(beta), A0 and sgn(mu) can be found here.
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, 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 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).
How to use
The macros should be easy to use. Just copy them to the directory you want to use. They are both located at: /project/atlas/users/nicoleru/Store/.
As long as you are only interested in the information currently checked by susygen.py, you only need to change susymap.py to form you own 5D susy grid. 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, when you feed it a susymap file. I will work on that. In the mean time, the root file used to make figure 1 should easily be convertible to other grids.