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 -> Xs
&gamma too.
Furthermore we have added the leading two-loop SUSY
corrections to &delta a&mu and made an
independent implementation
of B
s -> Xs &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 ->
Xs &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 -> Xs l+
l- (see
hep-ph/0410155).
Finally electroweak observables as the W gauge boson mass
MW
and the sine squared of the effective weak mixing angle
sin2&thetaeff
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.
|