MadDM

Absence of the differential recoil energy rate in MadDM

Asked by Adil JUEID

Dear MadDM developers,

I have tried to compute the differential rates in dcostheta, dE, and in dE dcostheta for DM direct detection experiment but I got 0.000 meaning that MadDM doesn't compute the spectra.

The syntax I used was
import model MajDM
generate direct_detection
launch

I did not change the default option for the direct_detection (option is directional which I believe is the one that necessary for the computation of differential recoil spectra).

Did I miss something in this process?

Thanks,
Cheers,
Adil

Question information

Language:
English Edit question
Status:
Solved
For:
MadDM Edit question
Assignee:
No assignee Edit question
Solved by:
Olivier Mattelaer
Solved:
Last query:
Last reply:
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said : 		#1

Your model is not known by our database of model.
Therefore I'm not able to comment.

Cheers,

Olivier

> On 20 Nov 2019, at 17:17, Adil JUEID <email address hidden> wrote:
>
> import model MajDM
> generate direct_detection
> launch

Revision history for this message
Adil JUEID (adil-hep) said : 		#2

Hi Olivier,

Thanks for your reply!

Here is a link to the model file
https://www.dropbox.com/s/ug0l0b7t85c8r9m/MajDM_SM.tar.xz?dl=0

Your help is much appreciated!

Cheers,
Adil

Revision history for this message
Best Olivier Mattelaer (olivier-mattelaer) said : 		#3

Can you try the following patch? (it does the trick for me)

=== modified file 'Templates/src/directional_detection.f'
--- Templates/src/directional_detection.f 2019-02-15 13:58:45 +0000
+++ Templates/src/directional_detection.f 2019-11-27 11:26:49 +0000
@@ -928,10 +928,11 @@
       Allocate (sigmaWN(n))

 ! Some variables
+ c = 3.d+5 ! Speed of light in km/s

       MD0 = MD*(GeV/cs**2)
       E0 = (1.d0/2.d0)*MD0*((v0/c)**2)*GeVtoKeV ! most probable kinetic energy of WIMPs
-
+
       sigmawnSI = sigmawn0SI*picobarn !WIMP-proton spin independent cross-section in cm^2
       sigmawpSI = sigmawp0SI*picobarn !WIMP-neutron spin independent cross-section in cm^2

@@ -947,7 +948,7 @@
       AMU = 0.932d0*(GeV/cs**2) ! Atomic Mass Units
       RhoD = rhoDM*GeV*cs**(-2)*cm**(-3) !Density of Dark Matter in our local part of the Galaxy.

- c = 3.d+5 ! Speed of light in km/s
+

 c form factor contribution terms for the spin independent calculation.
       fn = sqrt((pi/4.d0)*sigmawnSI*(1.d0/mu_n)**2.d0) ! cm/GeV

Thanks for the bug report,

Olivier

Revision history for this message
Adil JUEID (adil-hep) said : 		#4

Thanks! That solved my problem!

Cheers,
Adil

Revision history for this message
Adil JUEID (adil-hep) said : 		#5

Thanks Olivier Mattelaer, that solved my question.