RootFitLib/DrawResponse.cpp at main · terawatt93/RootFitLib · 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
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
#include "FitLib.hh"

void DrawResponse()
{
	vector<double> Peaks={263,278,305.5,355.5,385};

	ResponseFunction resp;


	resp.ReadResponse("LaBr");
	resp.UseExternal=true;
	resp.SubstrateFunction=TF1("sub","pol1",0,10000);
	resp.GenerateResponseFunction(100,500,Peaks);
	resp.SetNpx(2000);


	/*resp.SetParameter(0,0);
	resp.SetParameter(1,0);
	resp.SetParameter(2,0);
	resp.SetParameter(3,0);
	resp.SetParameter(4,0.2);
	resp.SetParameter(5,5000);
	resp.SetParameter(6,0.2);
	resp.SetParameter(7,5000);
	resp.SetParameter(8,0.2);
	resp.SetParameter(9,5000);
	resp.SetParameter(10,0.2);
	resp.SetParameter(11,5000);*/

	/*resp.SetParameter(0,0);
	resp.SetParameter(1,0.2);
	resp.SetParameter(2,5000);
	resp.SetParameter(3,0.2);
	resp.SetParameter(4,5000);
	resp.SetParameter(5,0.2);
	resp.SetParameter(6,5000);
	resp.SetParameter(7,0.2);
	resp.SetParameter(8,5000);*/


	/*par0 200000
par1 -0.1
par2 0.135902
par3 0.4878
par4 4.75842e+06
par5 9.99996e+06
par6 0.458399
par7 5.23115e+06
par8 9.99998e+06
par9 0.498138
par10 6.13842e+06
par11 9.99999e+06
par12 0.494848
par13 5.1607e+06
par14 9.99998e+06
par15 0.468554
par16 4.5638e+06
par17 9.99996e+06
*/

resp.SetParameter(0,200000);
resp.SetParameter(1,-0.1);
resp.SetParameter(2,0.135902);
resp.SetParameter(3,0.4878);
resp.SetParameter(4,4.75842e+06);
resp.SetParameter(5,9.99996e+06);
resp.SetParameter(6,0.458399);
resp.SetParameter(7,5.23115e+06);
resp.SetParameter(8,9.99998e+06);
resp.SetParameter(9,0.498138);
resp.SetParameter(10,6.13842e+06);
resp.SetParameter(11,9.99999e+06);
resp.SetParameter(12,0.494848);
resp.SetParameter(13,5.1607e+06);
resp.SetParameter(14,9.99998e+06);
resp.SetParameter(15,0.468554);
resp.SetParameter(16,4.5638e+06);
resp.SetParameter(17,9.99996e+06);

	/*resp.SetParameter(0,200000);
	resp.SetParameter(1,0.2);
	resp.SetParameter(2,5000);
	resp.SetParameter(3,5000);
	resp.SetParameter(4,0.2);
	resp.SetParameter(5,5000);
	resp.SetParameter(6,5000);
	resp.SetParameter(7,0.2);
	resp.SetParameter(8,5000);
	resp.SetParameter(9,5000);
	resp.SetParameter(10,0.2);
	resp.SetParameter(11,5000);
	resp.SetParameter(12,5000);*/

	resp.Draw();
	resp.PrintParameters();
	gPad->GetCanvas()->Print("func.pdf","pdf");
}