.test_CRESST-II.cfg

//DaMaSCUS Crust - Configuration File

//Simulation input parameter
	simID		=	"test";		//Simulation ID
	samplesize	=	100;				//velocity sample size

//Light DM option (if true, the nuclear form factor is assumed to be 1 speeding up the simulations)
	LDM			=	true;				//Light DM Option: true or false. 

//Interaction type
	DM_FormFactor 	=	"Contact";		//Options: "General","Contact", "Electric-Dipole", "Long-Range"
	ZorA			=	"A";		//couple to charge or mass
	
	Screening		=	false;		//Coupling to charge. Always true for "Electric-Dipole" and "Long-Range"
	mMediator		=	0.0;		//mediator mass in MeV. Only relevant for "General".

//Analysis parameter
	cl 			=	0.9;				//Certainty level for constraints
	experiment	=	"CRESST-II";		//Options: "DAMIC", "XENON1T", "CRESST-II", "CRESST-surface","Semiconductor"

//Halo parameters
	rhoDM		=	0.3;				//in GeV/cm^3
	v0			=	220.0;				//in km/sec
	vEarth 		=	240.0;				//in km/sec
	vEscape		=	544.0;				//in km/sec

//Options for "Semiconductor": DM Electron scattering experiments with semiconductors
	target 		=	"Si";				//Options: "Si" or "Ge"
	Qthreshold 	=	4;					//in electron hole pairs
	exposure	=	1.0;				//in gram years
	efficiency	=	1.0;				//flat efficiency
	events		=	0;					//observed events

//Mass Parameter scan
	mMin		=	0.5;				//in GeV
	mMax		=	20.0;				//in GeV
	masses		=	1;
//Cross section scan
	dSigma		=	0.5;				//Initial CS stepsize (log10)
	
//Variation reduction
//Importance sampling parameters
	is_angle		=	0.6; 			//Scattering Angle (only for DM_FormFactor= "Contact")
	is_mfp			=	0.6;			//MFP
//Geometric Importance Splitting
	splitting 		=	false;		//(de-)activate importance splitting	
	splits 			=	1.0;			//number of splits at splitting surface
	kappa 			=	3.0;		//choose between 1.0 and 10.0. Determines how fast the number of importance boundaries increases

//Shielding
	//Atmosphere
		atmosphere  	=	false;
		atmo_layers		=	1;		//number of atmospheric layers
		altitude		=	0.0;	//in meter (0.0 for underground experiments)
	//User defined layers
		layers 		=	
		(
			{
				name		=	"Crust";
				thickness	=	1400.0;			//in m
				density		=	2.7;			//in g/cm^3
				composition =	
				(
					(0.466,	8.0,	16.0),		//Oxygen			O
					(0.277,	14.0,	28.0),		//Silicon			Si
					(0.081,	13.0,	27.0),		//Aluminum			Al
					(0.05,	26.0,	56.0),		//Iron				Fe
					(0.036,	20.0,	40.0),		//Calcium			Ca
					(0.028,	19.0,	39.0),		//Kalium/Potassium	Ka
					(0.026,	11.0,	23.0),		//Natrium/Sodium	Na
					(0.021,	12.0,	24.0)		//Magnesium			Mg
				);
			}
		);