DarkSusy

DarkSusy is a fortran 77 package developed mainly for dark matter computations in the MSSM framework. Although it can be used as a full spectrum generator it works, in general, at tree level. Thus the whole spectrum is passed from SoftSusy, where a full 1-loop computation is performed.

The package computes the relic density of the lightest neutralino as the lightest supersymmetric particle (LSP) including all posible annihilations and coannihilations channels. Also it computes dark matter direct and indirect detection observables. Direct detection observables computed by SuperBayes are the spin-independent (dependent) nucleon-neutralino cross sections. Those are not included in the likelihood analysis due to the large uncertainties arising from astrophysics and nuclear physics in their calculation. We define the nucleon-neutralino spin-independent cross section as the largest among the proton and neutron ones. For spin-dependent both proton and neutron-neutralino rates are output. Currently the direct detection observables are computed as implemented independently for us. For details on our implementation see hep-ph/0208069. For plotting purposes current experimental data have been added in SuperBayeS/data.
Indirect detection quantities are also implemented as it is described in astro-ph/0406204, and as for the direct detection ones, those are not included in the likelihood analysis. Details about settings for their computations can be found on the interface file described below.

Although the main goal of the code is to be a dark matter calculator, it computes collider observables as the SUSY contributions to &delta a_&mu at one-loop level and the branching ratio of the rare process B _s -> X_s &gamma too.
Furthermore we have added the leading two-loop SUSY corrections to &delta a_&mu and made an independent implementation of B _s -> X_s &gamma where the new Standard Model result by Misiak et al. ( hep-ph/0609241) is added. Beyond this, we have implemented the computation of the branching ratio of the process B _s -> &mu⁺ &mu^- and the sign of the total amplitude for the decay B _s -> X_s &gamma and that of the Standard Model. Excepting in highly non-standard scenarios they have to be the same from experimental constraints on the branching ratio of the process B -> X_s l⁺ l^- (see hep-ph/0410155). Finally electroweak observables as the W gauge boson mass M_W and the sine squared of the effective weak mixing angle sin²&theta_eff are computed. The details on how all these quantities are computed are stated in hep-ph/0602028 and hep-ph/0611173.

Moreover the package is interfaced to the fortran 77 code Hdecay which calculates the decay widths and branching ratios of the Standard Model Higgs boson and on the neutral and charged Higgses of the MSSM as an alternative way to FeynHiggs.

The comunication with the MonteCarlo is controlled by an interface routine which is set in the source directory of the package. A full debugging is available from there. The current output from the package can be seen in the Outputs part of the dictionary provided in the documentation of the SuperBayeS package.