MappingmSUGRA

From Atlas Wiki
Revision as of 18:52, 12 March 2008 by Nicoleru (talk | contribs)
Jump to navigation Jump to search

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 , , and which is the cold dark matter relic dencity this SUSY point would produce.

If the susy point uncer 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 cheackes 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 pahsespace 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 susypoint.