Documentation for version 2.2.CPC (December 2008)

For more information on micrOMEGAs please consult the following references:

For micromegas_2.2 :
G. Bélanger, F. Boudjema, A. Pukhov, A. Semenov, arXiv:0803.2360 [hep-ph]
and
 G. Bélanger, F. Boudjema, A. Pukhov, A. Semenov, Comput.Phys.Commun.176:367-382,2007 hep-ph/0607059

1. For the MSSM model see micromegas_1.3, if you are using this module please cite
G. Bélanger, F. Boudjema, A. Pukhov and A. Semenov, Comput. Phys. Commun. 174 (2006) 577; hep-ph/0405253.
and
G. Bélanger, F. Boudjema, A. Pukhov and A. Semenov, Comput. Phys. Commun. 149 (2002) 103; hep-ph/0112278.

2. For a dexcription and physics studies in the NMSSM see
G. Bélanger, F. Boudjema,  C. hugonie, A. Pukhov, A. Semenov, JCAP0509:001(2005); hep-ph/0505142

 3. For a description and physics studies in the CPVMSSM see
G. Bélanger, F. Boudjema,  S. Kraml, A. Pukhov, A. Semenov, Phys. Rev. D73:115007 (2006); hep-ph/0604150
4. The RHN model is a generic model with a Dirac right-handed neutrino, for a description and physics studies see
G. Belanger, A.Pukhov, G.Servant, JCAP 0801:009 (2008); arXiv:0706.0526 [hep-ph].

5. If you are using the  little Higgs model, please cite
 A. Belyaev, C.R.Chen, K. Tobe, C.P. Yuan,  Phys.Rev. D74:115020 (2006)

 For more information see the README files


Changes from micromegas 2.2
 

1. Using PYTHIA simulation we recomputed the spectra of gamma, positron, antiproton, neutrinos in the channels
  DM,DM-> e,E
       -> u,U
       -> d,D
       -> s,S
In the previous version for light leptons(quarks) we were using instead the spectra for 
    DM,DM-> m,M
    DM,DM-> c,C
With the new spectra we find large corrections in models where DM is not a Majorana particle. On the other hand 
in MSSM-like models the light fermion annihilation channels are suppressed so do not need to be compute precisely. 


2.The version of CalCHEP has been updated.

3. When implementing a new model there is no need to put a  '*' before the masses of  Odd paricles in 
the Variable and Constraints tables. micrOMEGAs now knows about all  particle masses before process generation.


4. The function findParam has been removed.  Use findVal(findValW) instead.
With findVal one can access  not only the free parameters of the model  but also all constrained
parameters that are listed above the last  mass parameter in the list of Constraints.
To use findVal for other constrained parameters add into Constraints model list  a comment line with
   %Local!
Then all parameters above this line will be accessible by findVal
 

5.The version of CalCHEP has been updated.

6. In  MSSM-EWSB the soft masses for fermions of first generation are now model parameters. In the function 
    readVarMSSM
the mass of the first generation are still taken automatically to be equal to those of the second generation.On the other hand
  readVar
reads all parameters. So if you use assignVal to redefine the soft mass for a second generation fermion remember to use
also assignVal to redefine that of the  first generation, otherwise it will be set to some default value NOT 
to the value of the second generation fermion.

7. In the function  nucleusRecoil0, the Fermi distribution used to  describe the spin-dependent nucleus  form factors
has been replaced by a Gauss distribution (see manual).

8. We have added 2 new  routines 

extern double nucleusRecoilAux(rho, vfv, A, Z, J, S00, S01, S11,Mwimp, cs_SI_P, cs_SI_N, cs_SD_P, cs_SD_N,dNdE);
extern double nucleusRecoil0Aux(rho, vfv,A, Z, J, Sp,Sn,Mwimp, cs_SI_P, cs_SI_N, cs_SD_P, cs_SD_N,dNdE);

which are similar to nucleusRecoil and nucleusRecoil0 but use as input
Mwimp   - LSP mass in [GeV]
cs_SI_P - proton  spin-independent cross section in [pb]
cs_SI_N - neutron spin-independent cross section in [pb]
cs_SD_P - proton  spin-dependent   cross section in [pb]
cs_SD_N - neutron spin-dependent   cross section in [pb]
These routines depend only on the parameters specified as their arguments.  In particular the couplings of WIMP-nucleon interactions are extracted from 
the cross section. In order  to simulate a sign difference between  proton and neutron couplings one has to enter a negative cross section. 

 
9. A new function 

setRecoilEnergyGrid(step, nStep)
 
redefines the recoil energy grid. The nucleusRecoil function retuns recoil energy distribution
dNdE[i] =  dN/dE_i
where  E_i= i*step*KeV  i=0,...,nStep-1