In [1]:

    

# Load in the module that does the heavy lifting
import CalcDiscrimination as CD


    

    




 
 **********************************
 *      CalcDiscrimination.py     *
 **********************************
 Code for calculating significance of discriminating between Dirac and Majorana DM...
 Version 1.0.2 - BJK 06/07/2017
 

In [2]:

    

# Pick some parameter values - see the paper for definitions
ensemble = "A"
mDM = 50.0    #DM mass
f = -0.995   #Interference parameter
r_np = 0.75   #Ratio of neutron to proton couplings


    

In [3]:

    

# Calculate the discrimination significances for 10 samps
# This will probably take about a minute or so depending on the parameters...
#CD.CalcDiscrim("A", mDM, f, r_np)
CD.CalcDiscrim(ensemble, mDM, f, r_np)


    

    




 Loading experiments for ensemble A ...
 
 DM mass [GeV]: 50.0
 lambda [GeV^-2]: [  1.35282138e-08  -1.00954295e-08   0.00000000e+00   1.01334697e-09]
 f = -0.995
 c_n/c_p = 0.75
 
 Ne(Xenon):  314.80595829 ; sig_p = 1e-46
 Ne(Argon):  177.115707225 ; sig_p = 1.2373669816e-46
 Ne(Silicon):  32.4760197221 ; sig_p = 6.63721752583e-46
 
 Generating 10 samples...
   Sample 1  - Discrimination significance: 5.00935142916 sigma
   Sample 2  - Discrimination significance: 4.58410708657 sigma
   Sample 3  - Discrimination significance: 3.8925486155 sigma
   Sample 4  - Discrimination significance: 4.63984025919 sigma
   Sample 5  - Discrimination significance: 5.17557832285 sigma
   Sample 6  - Discrimination significance: 4.66555903788 sigma
   Sample 7  - Discrimination significance: 5.42110285602 sigma
   Sample 8  - Discrimination significance: 4.52742153712 sigma
   Sample 9  - Discrimination significance: 5.16371596128 sigma
   Sample 10  - Discrimination significance: 4.30115625824 sigma
 Median significance: 4.65269964853

In [ ]: