file frontends/DDCalc_2_0_0.hpp
[No description available] More…
Defines
| Name | |
|---|---|
| BACKENDNAME | |
| BACKENDLANG | |
| VERSION | |
| SAFE_VERSION | |
| REFERENCE |
Detailed Description
Author:
- Lars A. Dal (l.a.dal@fys.uio.no)
- Christopher Savage (chris@savage.name)
- Pat Scott (p.scott@imperial.ac.uk)
- Sebastian Wild (sebastian.wild@ph.tum.de)
- Felix Kahlhoefer (felix.kahlhoefer@desy.de)
Date:
- 2014 Jul
- 2014 Sept
- 2015 Jan,Feb,June
- 2016 Apr
- 2016 Aug
- 2016 Aug
Frontend header for the DDCalc backend.
Compile-time registration of available functions and variables from this backend.
Authors (add name and date if you modify):
Macros Documentation
define BACKENDNAME
#define BACKENDNAME DDCalc
define BACKENDLANG
#define BACKENDLANG Fortran
define VERSION
#define VERSION 2.0.0
define SAFE_VERSION
#define SAFE_VERSION 2_0_0
define REFERENCE
#define REFERENCE GAMBIT:2017fax,GAMBIT:2018eea
Source code
// GAMBIT: Global and Modular BSM Inference Tool
// *********************************************
/// \file
///
/// Frontend header for the DDCalc backend.
///
/// Compile-time registration of available
/// functions and variables from this backend.
///
/// *********************************************
///
/// Authors (add name and date if you modify):
///
/// \author Lars A. Dal
/// (l.a.dal@fys.uio.no)
/// \date 2014 Jul
///
/// \author Christopher Savage
/// (chris@savage.name)
/// \date 2014 Sept
/// \date 2015 Jan,Feb,June
///
/// \author Pat Scott
/// (p.scott@imperial.ac.uk)
/// \date 2016 Apr
///
/// \author Sebastian Wild
/// (sebastian.wild@ph.tum.de)
/// \date 2016 Aug
///
/// \author Felix Kahlhoefer
/// (felix.kahlhoefer@desy.de)
/// \date 2016 Aug
///
/// *********************************************
// Identify backend
#define BACKENDNAME DDCalc
#define BACKENDLANG Fortran
#define VERSION 2.0.0
#define SAFE_VERSION 2_0_0
#define REFERENCE GAMBIT:2017fax,GAMBIT:2018eea
// Load it
LOAD_LIBRARY
// Allow this version of DDCalc for all models but fermionic HP
BE_ALLOW_MODELS(MSSM63atQ, MSSM63atMGUT)
BE_ALLOW_MODELS(ScalarSingletDM_Z2,ScalarSingletDM_Z2_running)
BE_ALLOW_MODELS(VectorSingletDM_Z2)
BE_ALLOW_MODELS(ScalarSingletDM_Z3,ScalarSingletDM_Z3_running)
// BACKEND FUNCTIONS =======================================
/* Import functions.
BE_FUNCTION arguments:
* Function name used within GAMBIT.
* Function return type (void for fortran subroutines).
* Argument type list (pointers for fortran routines).
* Symbol name in compiled object file (see below).
* Capability name.
Naming conventions for the object/library symbols of Fortran module
routines are typically:
__<modulename>_MOD_<routinename> [gfortran]
<modulename>_mp_<routinename>_ [ifort]
where the module and routine names are in lower case. To avoid
compiler-dependendent symbol names, BIND() statements are used in
the Fortran source code to explicitly specify the symbol names.
[n.b.: For non-module routines, '<routinename>_' (again in lower
case) is the convention for both compilers.] We take as our naming
convention for the externally-accessible DDCalc routines:
C_DDCalc_<modulename(kinda)>_routinename>
The 'C' is to signify routines intended for calling from C/C++
(argument and return types are expressly declared to be type-
compatible with the C bool, int, and double types).
*/
// Default initialisation of three main classes via factory functions
BE_FUNCTION(DDCalc_InitWIMP, int, (), "C_DDWIMP_ddcalc_initwimp", "InitWIMP")
BE_FUNCTION(DDCalc_InitHalo, int, (), "C_DDHalo_ddcalc_inithalo", "InitHalo")
BE_FUNCTION(DDCalc_InitDetector, int, (const bool&), "C_DDExperiments_ddcalc_initdetector", "InitDetector")
// Initialization (specific experimental factory functions).
BE_FUNCTION(XENON100_2012_Init, int, (), "C_DDCalc_xenon100_2012_init", "XENON100_2012_Init")
BE_FUNCTION(XENON1T_2017_Init, int, (), "C_DDCalc_xenon1t_2017_init", "XENON1T_2017_Init")
BE_FUNCTION(XENON1T_2018_Init, int, (), "C_DDCalc_xenon1t_2018_init", "XENON1T_2018_Init")
BE_FUNCTION(LUX_2013_Init, int, (), "C_DDCalc_lux_2013_init", "LUX_2013_Init")
BE_FUNCTION(LUX_2016_Init, int, (), "C_DDCalc_lux_2016_init", "LUX_2016_Init")
BE_FUNCTION(PandaX_2016_Init, int, (), "C_DDCalc_pandax_2016_init", "PandaX_2016_Init")
BE_FUNCTION(PandaX_2017_Init, int, (), "C_DDCalc_pandax_2017_init", "PandaX_2017_Init")
BE_FUNCTION(LUX_2015_Init, int, (), "C_DDCalc_lux_2015_init", "LUX_2015_Init")
BE_FUNCTION(PICO_2L_Init, int, (), "C_DDCalc_pico_2l_init", "PICO_2L_Init")
BE_FUNCTION(PICO_60_Init, int, (), "C_DDCalc_pico_60_init", "PICO_60_Init")
BE_FUNCTION(PICO_60_2017_Init, int, (), "C_DDCalc_pico_60_2017_init", "PICO_60_2017_Init")
BE_FUNCTION(SuperCDMS_2014_Init, int, (), "C_DDCalc_supercdms_2014_init", "SuperCDMS_2014_Init")
BE_FUNCTION(CDMSlite_Init, int, (), "C_DDCalc_cdmslite_init", "CDMSlite_Init")
BE_FUNCTION(SIMPLE_2014_Init, int, (), "C_DDCalc_simple_2014_init", "SIMPLE_2014_Init")
BE_FUNCTION(CRESST_II_Init, int, (), "C_DDCalc_cresst_ii_init", "CRESST_II_Init")
BE_FUNCTION(LZ_Init, int, (), "C_DDCalc_lz_init", "LZ_Init")
BE_FUNCTION(PICO_500_Init, int, (), "C_DDCalc_pico_500_init", "PICO_500_Init")
BE_FUNCTION(DarkSide_50_Init, int, (), "C_DDCalc_darkside_50_init", "DarkSide_50_Init")
BE_FUNCTION(DARWIN_Init, int, (), "C_DDCalc_darwin_init", "DARWIN_Init")
//BE_FUNCTION(DARWIN_Ar_Init, int, (), "C_DDCalc_darwin_ar_init", "DARWIN_Ar_Init")
//BE_FUNCTION(DARWIN_Xe_Init, int, (), "C_DDCalc_darwin_xe_init", "DARWIN_Xe_Init")
// Set halo parameters (Standard Halo Model):
// rho [GeV/cm^3], vrot [km/s], v0 [km/s], vesc [km/s]
// This need only be called once at the beginning if the halo parameters will not be modified during a scan. Default
// values are already set via DDCalc_InitHalo routine, so it need not be called at all if the default values are to be used.
BE_FUNCTION(DDCalc_SetSHM, void, (const int&, const double&, const double&, const double&, const double&), "C_DDCalc_ddcalc_setshm", "SetSHM")
// Set the WIMP mass and couplings for the Higgs portal DM models.
// *higgsportal: mass, fsp, fsn, app, apn
// Units: mass [GeV]; f [GeV^-2] = pure scalar coupling; a [GeV^-2] = pure pseudoscalar coupling
// Convention: f = G/2 where G is the effective 4 vertex DM-nucleon coupling.
BE_FUNCTION(DDCalc_SetWIMP_higgsportal, void, (const int&, const double&, const double&, const double&, const double&, const double&), "C_DDCalc_ddcalc_setwimp_higgsportal", "SetWIMP_higgsportal")
// Get the WIMP mass and couplings for the Higgs portal DM models.
BE_FUNCTION(DDCalc_GetWIMP_higgsportal, void, (const int&, double&, double&, double&, double&, double&), "C_DDCalc_ddcalc_getwimp_higgsportal", "GetWIMP_higgsportal")
// Set the WIMP mass and couplings/cross-sections for standard SI/SD scattering.
// There are three versions, depending on how the couplings are specified:
// * mfa: mass, fp, fn, ap, an
// * mG: mass, Gp_SI, Gn_SI, Gp_SD, Gn_SD
// * msigma: mass, sigmapSI, sigmanSI, sigmapSD, sigmanSD
// Units: mass [GeV]; f, G [GeV^-2]; a [unitless]; sigma [pb]
// Here, f & a are the typical WIMP-nucleon couplings for spin-independent (SI) and spin-dependent (SD) interactions.
// The G's are the effective 4 fermion vertex couplings, related to f & a by a normalization factor. The sigmas are WIMP-
// nucleon scattering cross-sections; a negative value can be used to indicated the corresponding coupling should be taken
// to be negative.
BE_FUNCTION(DDCalc_SetWIMP_mfa, void, (const int&, const double&, const double&, const double&, const double&, const double&), "C_DDCalc_ddcalc_setwimp_mfa", "SetWIMP_mfa")
BE_FUNCTION(DDCalc_SetWIMP_mG, void, (const int&, const double&, const double&, const double&, const double&, const double&), "C_DDCalc_ddcalc_setwimp_mg", "SetWIMP_mG")
BE_FUNCTION(DDCalc_SetWIMP_msigma, void, (const int&, const double&, const double&, const double&, const double&, const double&), "C_DDCalc_ddcalc_setwimp_msigma", "SetWIMP_msigma")
// Get the WIMP mass and couplings/cross-sections. Same signature and units as above for setters. The only difference is
// that the WIMP-nucleon cross-sections are always positive (physical) values.
BE_FUNCTION(DDCalc_GetWIMP_mfa, void, (const int&,double&,double&,double&,double&,double&), "C_DDCalc_ddcalc_getwimp_mfa", "GetWIMP_mfa")
BE_FUNCTION(DDCalc_GetWIMP_mG, void, (const int&,double&,double&,double&,double&,double&), "C_DDCalc_ddcalc_getwimp_mg", "GetWIMP_mG")
BE_FUNCTION(DDCalc_GetWIMP_msigma, void, (const int&,double&,double&,double&,double&,double&), "C_DDCalc_ddcalc_getwimp_msigma", "GetWIMP_msigma")
// Set the WIMP mass, spin, and coupling structure within the non-relativistic effective theory of DM-nucleon interactions.
// - SetWIMP_NREffectiveTheory initializes a WIMP within the non-relativistic effective theory setup, setting all coefficients to zero.
// Arguments are the WIMP index, the mass of the WIMP in GeV, and the spin of the WIMP.
// - SetNRCoefficient sets the coefficient of a single operator to a given value.
// Arguments are:
// (1) the WIMP index
// (2) The operator index, i.e. an integer specifying the non-relativistic operator, e.g. 6 for O_6.
// For the specific cases of O_1 and O_4 one can also use the operators (q^2/mp^2) * O_1 and (q^2/mp^2) * O_4,
// by passing -1 and -4, respectively.
// (3) The isospin index: 0 for the isoscalar and 1 for the isovector component of the operator.
// (4) The desired value of the operator coefficient in units GeV^(-2).
BE_FUNCTION(DDCalc_SetWIMP_NREffectiveTheory, void, (const int&,const double&,const double&), "C_DDCalc_ddcalc_setwimp_nreffectivetheory", "SetWIMP_NREffectiveTheory")
BE_FUNCTION(DDCalc_SetNRCoefficient, void, (const int&,const int&,const int&,const double&), "C_DDCalc_ddcalc_setnrcoefficient", "SetNRCoefficient")
// Get the values of the isoscalar and isovector part of a given non-relativistic operator.
// Arguments are:
// (1) the WIMP index
// (2) the operator index (see description fir SetNRCoefficient above)
// (3) gives the value of the isoscalar component of the operator, in units GeV^(-2)
// (4) gives the value of the isovector component of the operator, in units GeV^(-2)
BE_FUNCTION(DDCalc_GetNRCoefficient, void, (const int&,const int&,double&,double&), "C_DDCalc_ddcalc_getnrcoefficient", "GetNRCoefficient")
// Specify the minimum recoil energy to be included in the rate calculations [keV]. Note the efficiency curves already account for
// detector and analysis thresholds regardless of this setting, so setting this to 0 keV (the default behavior when initialization is
// performed) does not imply that very low energy recoils actually contribute to the signal.
BE_FUNCTION(DDCalc_SetDetectorEmin, void, (const int&, const double&), "C_DDCalc_ddcalc_setdetectoremin", "DD_SetDetectorEmin")
// Calculation function; should be run once for each experiment and model prior to using event and likelihood routines below.
BE_FUNCTION(DDCalc_CalcRates, void, (const int&, const int&, const int&), "C_DDRates_ddcalc_calcrates", "DD_CalcRates")
// Observed and expected events, likelihoods and p values.
// Events: observed events
// Background: average background expectation
// Signal: average signal expectation
// LogLikelihood: log of the likelihood (not -2lnL)
// LogPValue: log of the p value
// Factor x by which sigma -> x*sigma would yield given p-value (given as log(p))
BE_FUNCTION(DDCalc_Events, int, (const int&), "C_DDRates_ddcalc_events", "DD_Events")
BE_FUNCTION(DDCalc_Background, double, (const int&), "C_DDRates_ddcalc_background", "DD_Background")
BE_FUNCTION(DDCalc_Signal, double, (const int&), "C_DDRates_ddcalc_signal", "DD_Signal")
BE_FUNCTION(DDCalc_Bins, int, (const int&), "C_DDRates_ddcalc_bins", "DD_Bins")
BE_FUNCTION(DDCalc_BinEvents, int, (const int&, const int&), "C_DDRates_ddcalc_binevents", "DD_BinEvents")
BE_FUNCTION(DDCalc_BinBackground, double, (const int&, const int&), "C_DDRates_ddcalc_binbackground", "DD_BinBackground")
BE_FUNCTION(DDCalc_BinSignal, double, (const int&, const int&), "C_DDRates_ddcalc_binsignal", "DD_BinSignal")
BE_FUNCTION(DDCalc_LogLikelihood, double, (const int&), "C_DDStats_ddcalc_loglikelihood", "DD_LogLikelihood")
// Do memory cleanup (nowhere to actually use these in GAMBIT proper, but they could be useful in standalones that hammer DDCalc).
BE_FUNCTION(DDCalc_FreeWIMPs, void, (), "C_DDUtils_ddcalc_freewimps", "FreeWIMPs")
BE_FUNCTION(DDCalc_FreeHalos, void, (), "C_DDUtils_ddcalc_freehalos", "FreeHalos")
BE_FUNCTION(DDCalc_FreeDetectors, void, (), "C_DDUtils_ddcalc_freedetectors", "FreeDetectorss")
BE_FUNCTION(DDCalc_FreeAll, void, (), "C_DDUtils_ddcalc_freeall", "FreeAll")
// DM mass, spin & conjugate nature, plus couplings and fraction of cosmological DM that is accounted for by model.
BE_INI_DEPENDENCY(mwimp, double)
BE_INI_DEPENDENCY(DD_couplings, DM_nucleon_couplings)
BE_INI_DEPENDENCY(spinwimpx2, unsigned int)
BE_INI_DEPENDENCY(wimp_sc, bool)
BE_INI_DEPENDENCY(RD_fraction, double)
BE_INI_DEPENDENCY(LocalHalo, LocalMaxwellianHalo)
// Convenience function for returning detector index given an analysis name.
BE_CONV_FUNCTION(DDCalc_Experiment, int, (const str&), "DD_Experiment")
// Convenience function for calling CalcRates with internally-initialised WIMP and halo objects.
BE_CONV_FUNCTION(DDCalc_CalcRates_simple, void, (const int&), "DD_CalcRates")
// Undefine macros to avoid conflict with other backends
#include "gambit/Backends/backend_undefs.hpp"
Updated on 2024-07-18 at 13:53:35 +0000