efficiencyTracker/functions.js at master · AdamGarnsworthy/efficiencyTracker · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
//functions//////////////////////////////////////////////////////////////////////////////////

//an eighth order polynomial for logEff in logE with coef. param = [0th order, 1st order, ..., 8th order].
//lo and hiParam are the 1-sigma extremes of the parameters
//returns string formatted for inclusion in dygraphs customGraph object.
function HPGeEfficiency(param, loParam, hiParam, logE){
	var i,
		logEff = 0,
		loDelta = 0,
		hiDelta = 0,
		eff;

	if(logE < Math.log(5)) return '0;0;0';

	for(i=0; i<9; i++){
		logEff += param[i]*Math.pow(logE,i);
		//loDelta += Math.pow((param[i] - loParam[i])*Math.pow(logE,i), 2);
		//hiDelta += Math.pow((hiParam[i] - param[i])*Math.pow(logE,i), 2);
	}
	//loDelta = Math.sqrt(loDelta);  //leave these errors out until we have a better grasp of what they should be.
	//hiDelta = Math.sqrt(hiDelta);

	eff = Math.exp(logEff);
	return (eff - eff*loDelta) + ';' + eff + ';' + (eff + eff*hiDelta);
}

function LaBrEfficiency(param, loParam, hiParam, logE){
	var i,
		logEff = 0,
		loDelta = 0,
		hiDelta = 0,
		eff;

	if(logE < Math.log(40)) return '0;0;0';

	for(i=0; i<9; i++){
		logEff += param[i]*Math.pow(logE,i);
		//loDelta += Math.pow((param[i] - loParam[i])*Math.pow(logE,i), 2);
		//hiDelta += Math.pow((hiParam[i] - param[i])*Math.pow(logE,i), 2);
	}
	//loDelta = Math.sqrt(loDelta);  //leave these errors out until we have a better grasp of what they should be.
	//hiDelta = Math.sqrt(hiDelta);

	eff = Math.exp(logEff);
	return (eff - eff*loDelta) + ';' + eff + ';' + (eff + eff*hiDelta);
}

function SiLiEfficiency(param, loParam, hiParam, logE){
	var i,
		loDelta = 0,
		hiDelta = 0,
		E = Math.exp(logE),
		eff = param[5]*param[0]*Math.exp(param[1]*Math.pow(E, param[2]))*(1 - Math.exp(param[3]*Math.pow(E, param[4])));

	if(logE < Math.log(15)) return '0;0;0';

	return (eff - eff*loDelta) + ';' + eff + ';' + (eff + eff*hiDelta);
}

function DESCANTefficiency(logE){
	if(logE >= Math.log(1000) && logE <= Math.log(5000))
		return '0.27;0.27;0.27';  //100% efficient * 27% geometric acceptance
	else{
		confirm('DESCANT Energy Out of Range', 'DESCANT neutron efficiencies are not reported below 1 MeV or above 5 MeV at this time.');
		return 'NaN;NaN;NaN'
	}
}

function SCEPTARefficiency(detector, logE){

	if(detector == 'SCEPTAR')
		return '0.8;0.8;0.8';
	else if(detector == 'SCEPTARZDS')
		return '0.65;0.65;0.65';
	else if(detector == 'SCEPTARPACES')
		return '0.4;0.4;0.4';
	else if(detector == 'PACESZDS')
		return '0.25;0.25;0.25';
}