.test_config_Dark_Photon.cfg

//Darphene - Configuration File

//ID
	ID		    =	".test_Dark_Photon";

// Run modus
    run_modus   =	"Custom"; 	//Options: "Energy-Spectrum", "Directional-Spectrum", "Daily-Modulation","All"
    MC_points   =	1000;
    grid_points =	100;

// Graphene Options
	carbon_wave_functions	= "RHF"; 		// Options: "Hydrogenic" and "Roothaan-Hartree-Fock"/"RHF"
											// Recommended: "RHF"
    work_function 			= 4.3; 			// in eV

//Dark matter particle
	DM_mass		  				=	100.0;		// in MeV
	DM_spin		  				=	0.5;
	DM_fraction					=	1.0;		// the DM particle's fractional abundance (set to 1.0 for 100%)

	DM_interaction				=	"Dark-Photon";		// Options: "Dark-Photon", "NREFT", "Electric-Dipole", "Magnetic-Dipole", "Anapole"

// Options for "Dark-Photon"
	DM_cross_section_electron	=	1.0e-37;	//in cm^2 (only relevant for dark photon model)
	DM_form_factor				=	"Contact";	// Options: "Contact", "Long-Range", "General"
	DM_mediator_mass			=	0.0;		// in MeV (only relevant if 'DM_form_factor' is "General")

// Options for "NREFT" (operator, effective coupling, form factor, parameter e.g. mediator mass [MeV])
    NREFT_couplings =   (
                        (1, 1.0, "Contact", 3.0),   
                        # (5, 1.0, "Contact", 3.0)        
    );
// Options for "Electric-Dipole", "Magnetic-Dipole", "Anapole"
    DM_coupling     =  1.0;		// either in GeV^-1 or GeV^-2 (depending on the interaction)

//Dark matter distribution
	DM_distribution 	=	"SHM";		//Options: "SHM", "SHM++"
	DM_local_density	=	0.4;		//in GeV / cm^3
	time 				=	6.0; 		// in hr.  t = 0 corresponds to the DM wind coming from the top.
	
	//Options for "SHM" and "SHM++"
		SHM_v0			=	238.0;					//in km/sec
		SHM_vObserver	=	(0.0, 0.0, 250.552);	//in km/sec
		SHM_vEscape		=	544.0;					//in km/sec
	//Options for "SHM++"
		SHMpp_eta		=	0.2;
		SHMpp_beta		=	0.9;

//Computation of exclusion limits
	exposure				=	1.0;	//in gram*year
	constraints_certainty	=	0.95;	//Certainty level
	constraints_mass_min	=	0.02;	//in GeV										
	constraints_mass_max	=	1.0;	//in GeV
	constraints_masses		=	10;