User Tools

Site Tools


idm_at_1loop

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
idm_at_1loop [2021/07/31 16:32]
shankha.banerjee [B. Interface for MicrOMEGAs for obtaining the relic]
idm_at_1loop [2021/08/01 18:24] (current)
boudjema [B. Interface for MicrOMEGAs for obtaining the relic]
Line 9: Line 9:
   -  Mathematica   -  Mathematica
   -  FeynArts, FormCalc and LoopTools [[http://​www.feynarts.de/​formcalc/​|(download here)]]   -  FeynArts, FormCalc and LoopTools [[http://​www.feynarts.de/​formcalc/​|(download here)]]
-  -  ​Micromegas ​[[http://​lapth.cnrs.fr/​micromegas/​|(download here)]]+  -  ​micrOMEGAs ​[[http://​lapth.cnrs.fr/​micromegas/​|(download here)]]
  
 (install these 3rd party codes/​programmes) (install these 3rd party codes/​programmes)
Line 25: Line 25:
  
  
-The upload consists of one single .zip file (IDM_modelfiles_masterfile.zip). {{ :​sloops_idm_public:​idm_modelfiles_masterfile.zip |Download it}}. +The upload consists of one single .zip file (IDM_modelfiles_masterfiles.zip). {{ :​sloops_idm_public:​idm_modelfiles_masterfiles.zip |Download it}}. 
  
 Unzip it. Upon unzip, you will have  Unzip it. Upon unzip, you will have 
Line 33: Line 33:
      ​* ​ SloopS-IDM_MSbar:​(model 1044)  is  for the Mixed OS- MSbar scheme for m_h< 2M_X      ​* ​ SloopS-IDM_MSbar:​(model 1044)  is  for the Mixed OS- MSbar scheme for m_h< 2M_X
      * SloopS-IDM_OS_CS:​(model103) is for the fully OS scheme for m_h > 2 M_ X and complex mass  relevant, for example, for resonant Higgs annihilation      * SloopS-IDM_OS_CS:​(model103) is for the fully OS scheme for m_h > 2 M_ X and complex mass  relevant, for example, for resonant Higgs annihilation
-     * ImprovedMO: contains the micrOMEGAs routine improveCS.c and other example ​files+     * ImprovedMO: contains the micrOMEGAs routine improveCS.c, 3 NLO data files (NLO_KL2_2_bb_P58.txt,​ NLO_KL2_2_mumu_P58.txt ​and NLO_KL2_2_ntnt_P58.txt,​ where KL2_2 stands for \lambda_2 = 2) and 3 LO data files (LO_bb_P58.txt,​ LO_mumu_P58.txt and LO_ntnt_P58.txt).
  
  
Line 55: Line 55:
       MH3 = 170D0       MH3 = 170D0
       MHC = 200D0 (mass of MH+)       MHC = 200D0 (mass of MH+)
-      mu2sq = 10000D0 (\mu2sq reconstructs \lambda ​see equation 2.8 in Paper 1)+      mu2sq = 10000D0 (\mu2sq reconstructs \lambda_L ​see equation 2.8 in Paper 1)
       la2 = 0.01D0 (\lambda2)       la2 = 0.01D0 (\lambda2)
  
Line 197: Line 197:
  
 You should have generated (for each process you want to correct) enough data for the corrected cross sections for well distributed values of the relative velocity v (Total energy/​Sqrt(S)). For each channel create a table of Pcm (the centre of mass momentum, in GeV) versus \sigma v (in GeV^{-2}). micrOMEGAs and improveCS.c requires us to convert the data in these units (Pcm in GeV and \sigma v in GeV^{-2}).This will replace the LO micrOMEGAs table through the MO routine improveCS.c (see the micrOMEGAs manual). See also the example improveCS.c along with the example files, in the main attachment. Keep all of these files in **micrOMEGAs_YourVersion/​IDM/​lib/​**. You should have generated (for each process you want to correct) enough data for the corrected cross sections for well distributed values of the relative velocity v (Total energy/​Sqrt(S)). For each channel create a table of Pcm (the centre of mass momentum, in GeV) versus \sigma v (in GeV^{-2}). micrOMEGAs and improveCS.c requires us to convert the data in these units (Pcm in GeV and \sigma v in GeV^{-2}).This will replace the LO micrOMEGAs table through the MO routine improveCS.c (see the micrOMEGAs manual). See also the example improveCS.c along with the example files, in the main attachment. Keep all of these files in **micrOMEGAs_YourVersion/​IDM/​lib/​**.
 +
 +To give you a bit more context, the example improveCS.c includes the paths to the corrected cross-sections. Please try to find this string "FILE *XX=". P58 (in paper1) has the processes XX -> WW* (5%) and AX -> f\bar{f} (95%) contributing. We, have chosen just two channels here for the example, b\bar{b} and \mu^+\mu^-. In addition, we are also providing the LO tables for b\bar{b} and \mu^+\mu^- **generated with micrOMEGAs version 5.0.11**. Such that the users can check, we are providing the \nu\nu LO and NLO tables, which are not included in the example version of improveCS.c. ​
 +The other processes can be included similarly.
 +
 +Now, once the improveCS.c and the data files are in the **micrOMEGAs_YourVersion/​IDM/​lib/​** and the **microMEGAs_YourVersion/​IDM/​work/​models/​func1.mdl** has been updated as explained above, come to the directory **microMEGAs_YourVersion/​IDM/​**. Then prepare a data file called data.par, which, for P58 reads
 +      MHX 58
 +      MH3 66 (H3 stands for A in our papers)
 +      MHC 110
 +      laL 0.0 (\lambda_L)
 +      la2 2.0 (\lambda_2)
 +      Mh 125
 +
 +Finally, the following steps need to be performed (from **microMEGAs_YourVersion/​IDM/​**) to extract the value of the relic density.
 +      make clean
 +      make main=main.c
 +      ./main data.par
 +
 +
 +
 ===== Help and Contact===== ===== Help and Contact=====
 If you need help contact If you need help contact
idm_at_1loop.1627741963.txt.gz · Last modified: 2021/07/31 16:32 by shankha.banerjee