This shows you the differences between two versions of the page.
| Both sides previous revision Previous revision Next revision | Previous revision | ||
|
idm_at_1loop [2021/07/31 16:51] shankha.banerjee |
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 2 NLO data files (NLO_KL2_2_bb_P58.txt and NLO_KL2_2_mumu_P58.txt, where KL2 is for \lambda_2 and the value is 2) | + | * 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 198: | Line 198: | ||
| 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 WW* (5%) and f\bar{f} (95%) contributing. We, have chosen just two channels here for the example, b\bar{b} and \mu^+\mu^-. We are also giving the LO tables for b\bar{b} and \mu^+\mu^- **generated with micrOMEGAs version 5.0.11**. We are also providing the \nu\nu tables such that you can include these in the code. Including the other processes can be done similarly. | + | 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 | ||
| - | 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 then come to the directory **microMEGAs_YourVersion/IDM/**. Then prepare a data file called data.par, which, for P58 | ||
| ===== Help and Contact===== | ===== Help and Contact===== | ||