Hall-D Software  alpha
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
ST_Tresolution/macros/SC_PTC_ResolvsZPlot.C
Go to the documentation of this file.
1 //// SC_PTC_ResolvsZPlot
2 /// Plot the time resolution Vs path length along the paddle.
3 
4 #include "Riostream.h"
5 #include "TCanvas.h"
6 #include "TCut.h"
7 TCanvas *PT_can[30];
8 void SC_PTC_ResolvsZPlot() {
9 
10  ifstream in;
11 
12  in.open("SCTR_ResVSz.txt");
13  Float_t sector, z, z_err, t, t_err,sigma,sigma_err;
14  Int_t nlines = 0;
15  std::ofstream SC_PTC;
16  SC_PTC.open ("Resol_1DFits.txt", std::ofstream::out);
17  auto f = TFile::Open("PTC_NS.root","RECREATE");
18  TNtuple ntuple("ntuple","data from ascii file","sector:z:z_err:sigma:sigma_err");
19  while (1) {
20  in >> sector >> z >> z_err >> sigma >> sigma_err;
21  if (!in.good()) break;
22  if (nlines < 5) printf("sector=%8f, z=%8f, z_err=%8f, sigma=%8f sigma_err=%8f\n",sector,z,z_err,sigma,sigma_err);
23  ntuple.Fill(sector,z,z_err,sigma,sigma_err);
24  nlines++;
25  }
26  auto c1 = new TCanvas("c1","c1",600,500);
27  c1->Divide(3,2);
28  c1->cd(1);
29  ntuple.Draw("z:z_err:sigma:sigma_err","sector >0 & sector <2 " );
30  TGraphErrors *gr = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
31  gr->Draw("ap"); //draw graph in current pad
32  gr->SetTitle("Sector 1");
33  gr->SetMaximum(0.4);
34  gr->SetMinimum(0.1);
35  TFitResultPtr sec1_r1 = gr->Fit("pol1","S","",10,30);
36  double sec1_r1_intcp = sec1_r1->Parameter(0);
37  double sec1_r1_slope = sec1_r1->Parameter(1);
38  TFitResultPtr sec1_r2 = gr->Fit("pol1","S+","",30,41);
39  double sec1_r2_intcp = sec1_r2->Parameter(0);
40  double sec1_r2_slope = sec1_r2->Parameter(1);
41  TFitResultPtr sec1_r3 = gr->Fit("pol1","S+","",41,54);
42  double sec1_r3_intcp = sec1_r3->Parameter(0);
43  double sec1_r3_slope = sec1_r3->Parameter(1);
44  TFitResultPtr sec1_r4 = gr->Fit("pol1","S+","",54,60);
45  double sec1_r4_intcp = sec1_r4->Parameter(0);
46  double sec1_r4_slope = sec1_r4->Parameter(1);
47  c1->cd(2);
48  ntuple.Draw("z:z_err:sigma:sigma_err","sector >1 & sector <3 " );
49  TGraphErrors *gr2 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
50  gr2->Draw("ap"); //draw graph in current pad
51  gr2->SetTitle("Sector 2");
52  gr2->SetMaximum(0.4);
53  gr2->SetMinimum(0.1);
54  TFitResultPtr sec2_r1 = gr2->Fit("pol1","S","",10,30);
55  double sec2_r1_intcp = sec2_r1->Parameter(0);
56  double sec2_r1_slope = sec2_r1->Parameter(1);
57  TFitResultPtr sec2_r2 = gr2->Fit("pol1","S+","",30,41);
58  double sec2_r2_intcp = sec2_r2->Parameter(0);
59  double sec2_r2_slope = sec2_r2->Parameter(1);
60  TFitResultPtr sec2_r3 = gr2->Fit("pol1","S+","",41,54);
61  double sec2_r3_intcp = sec2_r3->Parameter(0);
62  double sec2_r3_slope = sec2_r3->Parameter(1);
63  TFitResultPtr sec2_r4 = gr2->Fit("pol1","S+","",54,60);
64  double sec2_r4_intcp = sec2_r4->Parameter(0);
65  double sec2_r4_slope = sec2_r4->Parameter(1);
66 c1->cd(3);
67  ntuple.Draw("z:z_err:sigma:sigma_err","sector >2 & sector <4 " );
68  TGraphErrors *gr3 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
69  gr3->Draw("ap"); //draw graph in current pad
70  gr3->SetTitle("Sector 3");
71  gr3->SetMaximum(0.4);
72  gr3->SetMinimum(0.1);
73  TFitResultPtr sec3_r1 = gr3->Fit("pol1","S","",10,30);
74  double sec3_r1_intcp = sec3_r1->Parameter(0);
75  double sec3_r1_slope = sec3_r1->Parameter(1);
76  TFitResultPtr sec3_r2 = gr3->Fit("pol1","S+","",30,41);
77  double sec3_r2_intcp = sec3_r2->Parameter(0);
78  double sec3_r2_slope = sec3_r2->Parameter(1);
79  TFitResultPtr sec3_r3 = gr3->Fit("pol1","S+","",41,54);
80  double sec3_r3_intcp = sec3_r3->Parameter(0);
81  double sec3_r3_slope = sec3_r3->Parameter(1);
82  TFitResultPtr sec3_r4 = gr3->Fit("pol1","S+","",54,60);
83  double sec3_r4_intcp = sec3_r4->Parameter(0);
84  double sec3_r4_slope = sec3_r4->Parameter(1);
85  c1->cd(4);
86  ntuple.Draw("z:z_err:sigma:sigma_err","sector >3 & sector <5 " );
87  TGraphErrors *gr4 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
88  gr4->Draw("ap"); //draw graph in current pad
89  gr4->SetTitle("Sector 4");
90  gr4->SetMaximum(0.4);
91  gr4->SetMinimum(0.1);
92  TFitResultPtr sec4_r1 = gr4->Fit("pol1","S","",10,30);
93  double sec4_r1_intcp = sec4_r1->Parameter(0);
94  double sec4_r1_slope = sec4_r1->Parameter(1);
95  TFitResultPtr sec4_r2 = gr4->Fit("pol1","S+","",30,41);
96  double sec4_r2_intcp = sec4_r2->Parameter(0);
97  double sec4_r2_slope = sec4_r2->Parameter(1);
98  TFitResultPtr sec4_r3 = gr4->Fit("pol1","S+","",41,54);
99  double sec4_r3_intcp = sec4_r3->Parameter(0);
100  double sec4_r3_slope = sec4_r3->Parameter(1);
101  TFitResultPtr sec4_r4 = gr4->Fit("pol1","S+","",54,60);
102  double sec4_r4_intcp = sec4_r4->Parameter(0);
103  double sec4_r4_slope = sec4_r4->Parameter(1);
104  c1->cd(5);
105  ntuple.Draw("z:z_err:sigma:sigma_err","sector >4 & sector <6 " );
106  TGraphErrors *gr5 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
107  gr5->Draw("ap"); //draw graph in current pad
108  gr5->SetTitle("Sector 5");
109  gr5->SetMaximum(0.4);
110  gr5->SetMinimum(0.1);
111  TFitResultPtr sec5_r1 = gr5->Fit("pol1","S","",10,30);
112  double sec5_r1_intcp = sec5_r1->Parameter(0);
113  double sec5_r1_slope = sec5_r1->Parameter(1);
114  TFitResultPtr sec5_r2 = gr5->Fit("pol1","S+","",30,41);
115  double sec5_r2_intcp = sec5_r2->Parameter(0);
116  double sec5_r2_slope = sec5_r2->Parameter(1);
117  TFitResultPtr sec5_r3 = gr5->Fit("pol1","S+","",41,54);
118  double sec5_r3_intcp = sec5_r3->Parameter(0);
119  double sec5_r3_slope = sec5_r3->Parameter(1);
120  TFitResultPtr sec5_r4 = gr5->Fit("pol1","S+","",54,60);
121  double sec5_r4_intcp = sec5_r4->Parameter(0);
122  double sec5_r4_slope = sec5_r4->Parameter(1);
123  c1->cd(6);
124  ntuple.Draw("z:z_err:sigma:sigma_err","sector >5 & sector <7 " );
125  TGraphErrors *gr6 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
126  gr6->Draw("ap"); //draw graph in current pad
127  gr6->SetTitle("Sector 6");
128  gr6->SetMaximum(0.4);
129  gr6->SetMinimum(0.1);
130  TFitResultPtr sec6_r1 = gr6->Fit("pol1","S","",10,30);
131  double sec6_r1_intcp = sec6_r1->Parameter(0);
132  double sec6_r1_slope = sec6_r1->Parameter(1);
133  TFitResultPtr sec6_r2 = gr6->Fit("pol1","S+","",30,41);
134  double sec6_r2_intcp = sec6_r2->Parameter(0);
135  double sec6_r2_slope = sec6_r2->Parameter(1);
136  TFitResultPtr sec6_r3 = gr6->Fit("pol1","S+","",41,54);
137  double sec6_r3_intcp = sec6_r3->Parameter(0);
138  double sec6_r3_slope = sec6_r3->Parameter(1);
139  TFitResultPtr sec6_r4 = gr6->Fit("pol1","S+","",54,60);
140  double sec6_r4_intcp = sec6_r4->Parameter(0);
141  double sec6_r4_slope = sec6_r4->Parameter(1);
142  auto c2 = new TCanvas("c2","c2",600,500);
143  c2->Divide(3,2);
144  c2->cd(1);
145  ntuple.Draw("z:z_err:sigma:sigma_err","sector >6 & sector <8 " );
146  TGraphErrors *gr7 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
147  gr7->Draw("ap"); //draw graph in current pad
148  gr7->SetTitle("Sector 7");
149  gr7->SetMaximum(0.4);
150  gr7->SetMinimum(0.1);
151  TFitResultPtr sec7_r1 = gr7->Fit("pol1","S","",10,30);
152  double sec7_r1_intcp = sec7_r1->Parameter(0);
153  double sec7_r1_slope = sec7_r1->Parameter(1);
154  TFitResultPtr sec7_r2 = gr7->Fit("pol1","S+","",30,41);
155  double sec7_r2_intcp = sec7_r2->Parameter(0);
156  double sec7_r2_slope = sec7_r2->Parameter(1);
157  TFitResultPtr sec7_r3 = gr7->Fit("pol1","S+","",41,54);
158  double sec7_r3_intcp = sec7_r3->Parameter(0);
159  double sec7_r3_slope = sec7_r3->Parameter(1);
160  TFitResultPtr sec7_r4 = gr7->Fit("pol1","S+","",54,60);
161  double sec7_r4_intcp = sec7_r4->Parameter(0);
162  double sec7_r4_slope = sec7_r4->Parameter(1);
163  c2->cd(2);
164  ntuple.Draw("z:z_err:sigma:sigma_err","sector >7 & sector <9 " );
165  TGraphErrors *gr8 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
166  gr8->Draw("ap"); //draw graph in current pad
167  gr8->SetTitle("Sector 8");
168  gr8->SetMaximum(0.4);
169  gr8->SetMinimum(0.1);
170  TFitResultPtr sec8_r1 = gr8->Fit("pol1","S","",10,30);
171  double sec8_r1_intcp = sec8_r1->Parameter(0);
172  double sec8_r1_slope = sec8_r1->Parameter(1);
173  TFitResultPtr sec8_r2 = gr8->Fit("pol1","S+","",30,41);
174  double sec8_r2_intcp = sec8_r2->Parameter(0);
175  double sec8_r2_slope = sec8_r2->Parameter(1);
176  TFitResultPtr sec8_r3 = gr8->Fit("pol1","S+","",41,54);
177  double sec8_r3_intcp = sec8_r3->Parameter(0);
178  double sec8_r3_slope = sec8_r3->Parameter(1);
179  TFitResultPtr sec8_r4 = gr8->Fit("pol1","S+","",54,60);
180  double sec8_r4_intcp = sec8_r4->Parameter(0);
181  double sec8_r4_slope = sec8_r4->Parameter(1);
182 
183 c2->cd(3);
184  ntuple.Draw("z:z_err:sigma:sigma_err","sector >8 & sector <10 " );
185  TGraphErrors *gr9 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
186  gr9->Draw("ap"); //draw graph in current pad
187  gr9->SetTitle("Sector 9");
188  gr9->SetMaximum(0.4);
189  gr9->SetMinimum(0.1);
190  TFitResultPtr sec9_r1 = gr9->Fit("pol1","S","",10,30);
191  double sec9_r1_intcp = sec9_r1->Parameter(0);
192  double sec9_r1_slope = sec9_r1->Parameter(1);
193  TFitResultPtr sec9_r2 = gr9->Fit("pol1","S+","",30,41);
194  double sec9_r2_intcp = sec9_r2->Parameter(0);
195  double sec9_r2_slope = sec9_r2->Parameter(1);
196  TFitResultPtr sec9_r3 = gr9->Fit("pol1","S+","",41,54);
197  double sec9_r3_intcp = sec9_r3->Parameter(0);
198  double sec9_r3_slope = sec9_r3->Parameter(1);
199  TFitResultPtr sec9_r4 = gr9->Fit("pol1","S+","",54,60);
200  double sec9_r4_intcp = sec9_r4->Parameter(0);
201  double sec9_r4_slope = sec9_r4->Parameter(1);
202 
203  c2->cd(4);
204  ntuple.Draw("z:z_err:sigma:sigma_err","sector >9 & sector <11 " );
205  TGraphErrors *gr10 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
206  gr10->Draw("ap"); //draw graph in current pad
207  gr10->SetTitle("Sector 10");
208  gr10->SetMaximum(0.4);
209  gr10->SetMinimum(0.1);
210  TFitResultPtr sec10_r1 = gr10->Fit("pol1","S","",10,30);
211  double sec10_r1_intcp = sec10_r1->Parameter(0);
212  double sec10_r1_slope = sec10_r1->Parameter(1);
213  TFitResultPtr sec10_r2 = gr10->Fit("pol1","S+","",30,41);
214  double sec10_r2_intcp = sec10_r2->Parameter(0);
215  double sec10_r2_slope = sec10_r2->Parameter(1);
216  TFitResultPtr sec10_r3 = gr10->Fit("pol1","S+","",41,54);
217  double sec10_r3_intcp = sec10_r3->Parameter(0);
218  double sec10_r3_slope = sec10_r3->Parameter(1);
219  TFitResultPtr sec10_r4 = gr10->Fit("pol1","S+","",54,60);
220  double sec10_r4_intcp = sec10_r4->Parameter(0);
221  double sec10_r4_slope = sec10_r4->Parameter(1);
222 
223  c2->cd(5);
224  ntuple.Draw("z:z_err:sigma:sigma_err","sector >10 & sector <12 " );
225  TGraphErrors *gr11 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
226  gr11->Draw("ap"); //draw graph in current pad
227  gr11->SetTitle("Sector 11");
228  gr11->SetMaximum(0.4);
229  gr11->SetMinimum(0.1);
230  TFitResultPtr sec11_r1 = gr11->Fit("pol1","S","",10,30);
231  double sec11_r1_intcp = sec11_r1->Parameter(0);
232  double sec11_r1_slope = sec11_r1->Parameter(1);
233  TFitResultPtr sec11_r2 = gr11->Fit("pol1","S+","",30,41);
234  double sec11_r2_intcp = sec11_r2->Parameter(0);
235  double sec11_r2_slope = sec11_r2->Parameter(1);
236  TFitResultPtr sec11_r3 = gr11->Fit("pol1","S+","",41,54);
237  double sec11_r3_intcp = sec11_r3->Parameter(0);
238  double sec11_r3_slope = sec11_r3->Parameter(1);
239  TFitResultPtr sec11_r4 = gr11->Fit("pol1","S+","",54,60);
240  double sec11_r4_intcp = sec11_r4->Parameter(0);
241  double sec11_r4_slope = sec11_r4->Parameter(1);
242  c2->cd(6);
243  ntuple.Draw("z:z_err:sigma:sigma_err","sector >11 & sector <13 " );
244  TGraphErrors *gr12 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
245  gr12->Draw("ap"); //draw graph in current pad
246  gr12->SetTitle("Sector 12");
247  gr12->SetMaximum(0.4);
248  gr12->SetMinimum(0.1);
249  TFitResultPtr sec12_r1 = gr12->Fit("pol1","S","",10,30);
250  double sec12_r1_intcp = sec12_r1->Parameter(0);
251  double sec12_r1_slope = sec12_r1->Parameter(1);
252  TFitResultPtr sec12_r2 = gr12->Fit("pol1","S+","",30,41);
253  double sec12_r2_intcp = sec12_r2->Parameter(0);
254  double sec12_r2_slope = sec12_r2->Parameter(1);
255  TFitResultPtr sec12_r3 = gr12->Fit("pol1","S+","",41,54);
256  double sec12_r3_intcp = sec12_r3->Parameter(0);
257  double sec12_r3_slope = sec12_r3->Parameter(1);
258  TFitResultPtr sec12_r4 = gr12->Fit("pol1","S+","",54,60);
259  double sec12_r4_intcp = sec12_r4->Parameter(0);
260  double sec12_r4_slope = sec12_r4->Parameter(1);
261  auto c3 = new TCanvas("c3","c3",600,500);
262  c3->Divide(3,2);
263  c3->cd(1);
264  ntuple.Draw("z:z_err:sigma:sigma_err","sector >12 & sector <14 " );
265  TGraphErrors *gr13 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
266  gr13->Draw("ap"); //draw graph in current pad
267  gr13->SetTitle("Sector 13");
268  gr13->SetMaximum(0.4);
269  gr13->SetMinimum(0.1);
270  TFitResultPtr sec13_r1 = gr13->Fit("pol1","S","",10,30);
271  double sec13_r1_intcp = sec13_r1->Parameter(0);
272  double sec13_r1_slope = sec13_r1->Parameter(1);
273  TFitResultPtr sec13_r2 = gr13->Fit("pol1","S+","",30,41);
274  double sec13_r2_intcp = sec13_r2->Parameter(0);
275  double sec13_r2_slope = sec13_r2->Parameter(1);
276  TFitResultPtr sec13_r3 = gr13->Fit("pol1","S+","",41,54);
277  double sec13_r3_intcp = sec13_r3->Parameter(0);
278  double sec13_r3_slope = sec13_r3->Parameter(1);
279  TFitResultPtr sec13_r4 = gr13->Fit("pol1","S+","",54,60);
280  double sec13_r4_intcp = sec13_r4->Parameter(0);
281  double sec13_r4_slope = sec13_r4->Parameter(1);
282  c3->cd(2);
283  ntuple.Draw("z:z_err:sigma:sigma_err","sector >13 & sector <15 " );
284  TGraphErrors *gr14 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
285  gr14->Draw("ap"); //draw graph in current pad
286  gr14->SetTitle("Sector 14");
287  gr14->SetMaximum(0.4);
288  gr14->SetMinimum(0.1);
289  TFitResultPtr sec14_r1 = gr14->Fit("pol1","S","",10,30);
290  double sec14_r1_intcp = sec14_r1->Parameter(0);
291  double sec14_r1_slope = sec14_r1->Parameter(1);
292  TFitResultPtr sec14_r2 = gr14->Fit("pol1","S+","",30,41);
293  double sec14_r2_intcp = sec14_r2->Parameter(0);
294  double sec14_r2_slope = sec14_r2->Parameter(1);
295  TFitResultPtr sec14_r3 = gr14->Fit("pol1","S+","",41,54);
296  double sec14_r3_intcp = sec14_r3->Parameter(0);
297  double sec14_r3_slope = sec14_r3->Parameter(1);
298  TFitResultPtr sec14_r4 = gr14->Fit("pol1","S+","",54,60);
299  double sec14_r4_intcp = sec14_r4->Parameter(0);
300  double sec14_r4_slope = sec14_r4->Parameter(1);
301 c3->cd(3);
302  ntuple.Draw("z:z_err:sigma:sigma_err","sector >14 & sector <16 " );
303  TGraphErrors *gr15 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
304  gr15->Draw("ap"); //draw graph in current pad
305  gr15->SetTitle("Sector 15");
306  gr15->SetMaximum(0.4);
307  gr15->SetMinimum(0.1);
308  TFitResultPtr sec15_r1 = gr15->Fit("pol1","S","",10,30);
309  double sec15_r1_intcp = sec15_r1->Parameter(0);
310  double sec15_r1_slope = sec15_r1->Parameter(1);
311  TFitResultPtr sec15_r2 = gr15->Fit("pol1","S+","",30,41);
312  double sec15_r2_intcp = sec15_r2->Parameter(0);
313  double sec15_r2_slope = sec15_r2->Parameter(1);
314  TFitResultPtr sec15_r3 = gr15->Fit("pol1","S+","",41,54);
315  double sec15_r3_intcp = sec15_r3->Parameter(0);
316  double sec15_r3_slope = sec15_r3->Parameter(1);
317  TFitResultPtr sec15_r4 = gr15->Fit("pol1","S+","",54,60);
318  double sec15_r4_intcp = sec15_r4->Parameter(0);
319  double sec15_r4_slope = sec15_r4->Parameter(1);
320  c3->cd(4);
321  ntuple.Draw("z:z_err:sigma:sigma_err","sector >15 & sector <17 " );
322  TGraphErrors *gr16 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
323  gr16->Draw("ap"); //draw graph in current pad
324  gr16->SetTitle("Sector 16");
325  gr16->SetMaximum(0.4);
326  gr16->SetMinimum(0.1);
327  TFitResultPtr sec16_r1 = gr16->Fit("pol1","S","",10,30);
328  double sec16_r1_intcp = sec16_r1->Parameter(0);
329  double sec16_r1_slope = sec16_r1->Parameter(1);
330  TFitResultPtr sec16_r2 = gr16->Fit("pol1","S+","",30,41);
331  double sec16_r2_intcp = sec16_r2->Parameter(0);
332  double sec16_r2_slope = sec16_r2->Parameter(1);
333  TFitResultPtr sec16_r3 = gr16->Fit("pol1","S+","",41,54);
334  double sec16_r3_intcp = sec16_r3->Parameter(0);
335  double sec16_r3_slope = sec16_r3->Parameter(1);
336  TFitResultPtr sec16_r4 = gr16->Fit("pol1","S+","",54,60);
337  double sec16_r4_intcp = sec16_r4->Parameter(0);
338  double sec16_r4_slope = sec16_r4->Parameter(1);
339  c3->cd(5);
340  ntuple.Draw("z:z_err:sigma:sigma_err","sector >16 & sector <18 " );
341  TGraphErrors *gr17 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
342  gr17->Draw("ap"); //draw graph in current pad
343  gr17->SetTitle("Sector 17");
344  gr17->SetMaximum(0.4);
345  gr17->SetMinimum(0.1);
346  TFitResultPtr sec17_r1 = gr17->Fit("pol1","S","",10,30);
347  double sec17_r1_intcp = sec17_r1->Parameter(0);
348  double sec17_r1_slope = sec17_r1->Parameter(1);
349  TFitResultPtr sec17_r2 = gr17->Fit("pol1","S+","",30,41);
350  double sec17_r2_intcp = sec17_r2->Parameter(0);
351  double sec17_r2_slope = sec17_r2->Parameter(1);
352  TFitResultPtr sec17_r3 = gr17->Fit("pol1","S+","",41,54);
353  double sec17_r3_intcp = sec17_r3->Parameter(0);
354  double sec17_r3_slope = sec17_r3->Parameter(1);
355  TFitResultPtr sec17_r4 = gr17->Fit("pol1","S+","",54,60);
356  double sec17_r4_intcp = sec17_r4->Parameter(0);
357  double sec17_r4_slope = sec17_r4->Parameter(1);
358  c3->cd(6);
359  ntuple.Draw("z:z_err:sigma:sigma_err","sector >17 & sector <19 " );
360  TGraphErrors *gr18 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
361  gr18->Draw("ap"); //draw graph in current pad
362  gr18->SetTitle("Sector 18");
363  gr18->SetMaximum(0.4);
364  gr18->SetMinimum(0.1);
365  TFitResultPtr sec18_r1 = gr18->Fit("pol1","S","",10,30);
366  double sec18_r1_intcp = sec18_r1->Parameter(0);
367  double sec18_r1_slope = sec18_r1->Parameter(1);
368  TFitResultPtr sec18_r2 = gr18->Fit("pol1","S+","",30,41);
369  double sec18_r2_intcp = sec18_r2->Parameter(0);
370  double sec18_r2_slope = sec18_r2->Parameter(1);
371  TFitResultPtr sec18_r3 = gr18->Fit("pol1","S+","",41,54);
372  double sec18_r3_intcp = sec18_r3->Parameter(0);
373  double sec18_r3_slope = sec18_r3->Parameter(1);
374  TFitResultPtr sec18_r4 = gr18->Fit("pol1","S+","",54,60);
375  double sec18_r4_intcp = sec18_r4->Parameter(0);
376  double sec18_r4_slope = sec18_r4->Parameter(1);
377  auto c4 = new TCanvas("c4","c4",600,500);
378  c4->Divide(3,2);
379  c4->cd(1);
380  ntuple.Draw("z:z_err:sigma:sigma_err","sector >18 & sector <20 " );
381  TGraphErrors *gr19 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
382  gr19->Draw("ap"); //draw graph in current pad
383  gr19->SetTitle("Sector 19");
384  gr19->SetMaximum(0.4);
385  gr19->SetMinimum(0.1);
386  TFitResultPtr sec19_r1 = gr19->Fit("pol1","S","",10,30);
387  double sec19_r1_intcp = sec19_r1->Parameter(0);
388  double sec19_r1_slope = sec19_r1->Parameter(1);
389  TFitResultPtr sec19_r2 = gr19->Fit("pol1","S+","",30,41);
390  double sec19_r2_intcp = sec19_r2->Parameter(0);
391  double sec19_r2_slope = sec19_r2->Parameter(1);
392  TFitResultPtr sec19_r3 = gr19->Fit("pol1","S+","",41,54);
393  double sec19_r3_intcp = sec19_r3->Parameter(0);
394  double sec19_r3_slope = sec19_r3->Parameter(1);
395  TFitResultPtr sec19_r4 = gr19->Fit("pol1","S+","",54,60);
396  double sec19_r4_intcp = sec19_r4->Parameter(0);
397  double sec19_r4_slope = sec19_r4->Parameter(1);
398  c4->cd(2);
399  ntuple.Draw("z:z_err:sigma:sigma_err","sector >19 & sector <21 " );
400  TGraphErrors *gr20 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
401  gr20->Draw("ap"); //draw graph in current pad
402  gr20->SetTitle("Sector 20");
403  gr20->SetMaximum(0.4);
404  gr20->SetMinimum(0.1);
405  TFitResultPtr sec20_r1 = gr20->Fit("pol1","S","",10,30);
406  double sec20_r1_intcp = sec20_r1->Parameter(0);
407  double sec20_r1_slope = sec20_r1->Parameter(1);
408  TFitResultPtr sec20_r2 = gr20->Fit("pol1","S+","",30,41);
409  double sec20_r2_intcp = sec20_r2->Parameter(0);
410  double sec20_r2_slope = sec20_r2->Parameter(1);
411  TFitResultPtr sec20_r3 = gr20->Fit("pol1","S+","",41,54);
412  double sec20_r3_intcp = sec20_r3->Parameter(0);
413  double sec20_r3_slope = sec20_r3->Parameter(1);
414  TFitResultPtr sec20_r4 = gr20->Fit("pol1","S+","",54,60);
415  double sec20_r4_intcp = sec20_r4->Parameter(0);
416  double sec20_r4_slope = sec20_r4->Parameter(1);
417 c4->cd(3);
418  ntuple.Draw("z:z_err:sigma:sigma_err","sector >20 & sector <22 " );
419  TGraphErrors *gr21 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
420  gr21->Draw("ap"); //draw graph in current pad
421  gr21->SetTitle("Sector 21");
422  gr21->SetMaximum(0.4);
423  gr21->SetMinimum(0.1);
424  TFitResultPtr sec21_r1 = gr21->Fit("pol1","S","",10,30);
425  double sec21_r1_intcp = sec21_r1->Parameter(0);
426  double sec21_r1_slope = sec21_r1->Parameter(1);
427  TFitResultPtr sec21_r2 = gr21->Fit("pol1","S+","",30,41);
428  double sec21_r2_intcp = sec21_r2->Parameter(0);
429  double sec21_r2_slope = sec21_r2->Parameter(1);
430 TFitResultPtr sec21_r3 = gr21->Fit("pol1","S+","",41,54);
431  double sec21_r3_intcp = sec21_r3->Parameter(0);
432  double sec21_r3_slope = sec21_r3->Parameter(1);
433  TFitResultPtr sec21_r4 = gr21->Fit("pol1","S+","",54,60);
434  double sec21_r4_intcp = sec21_r4->Parameter(0);
435  double sec21_r4_slope = sec21_r4->Parameter(1);
436  c4->cd(4);
437  ntuple.Draw("z:z_err:sigma:sigma_err","sector >21 & sector <23 " );
438  TGraphErrors *gr22 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
439  gr22->Draw("ap"); //draw graph in current pad
440  gr22->SetTitle("Sector 22");
441  gr22->SetMaximum(0.4);
442  gr22->SetMinimum(0.1);
443  TFitResultPtr sec22_r1 = gr22->Fit("pol1","S","",10,30);
444  double sec22_r1_intcp = sec22_r1->Parameter(0);
445  double sec22_r1_slope = sec22_r1->Parameter(1);
446  TFitResultPtr sec22_r2 = gr22->Fit("pol1","S+","",30,41);
447  double sec22_r2_intcp = sec22_r2->Parameter(0);
448  double sec22_r2_slope = sec22_r2->Parameter(1);
449  TFitResultPtr sec22_r3 = gr22->Fit("pol1","S+","",41,54);
450  double sec22_r3_intcp = sec22_r3->Parameter(0);
451  double sec22_r3_slope = sec22_r3->Parameter(1);
452  TFitResultPtr sec22_r4 = gr22->Fit("pol1","S+","",54,60);
453  double sec22_r4_intcp = sec22_r4->Parameter(0);
454  double sec22_r4_slope = sec22_r4->Parameter(1);
455  c4->cd(5);
456  ntuple.Draw("z:z_err:sigma:sigma_err","sector >22 & sector <24 " );
457  TGraphErrors *gr23 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
458  gr23->Draw("ap"); //draw graph in current pad
459  gr23->SetTitle("Sector 23");
460  gr23->SetMaximum(0.4);
461  gr23->SetMinimum(0.1);
462  TFitResultPtr sec23_r1 = gr23->Fit("pol1","S","",10,30);
463  double sec23_r1_intcp = sec23_r1->Parameter(0);
464  double sec23_r1_slope = sec23_r1->Parameter(1);
465  TFitResultPtr sec23_r2 = gr23->Fit("pol1","S+","",30,41);
466  double sec23_r2_intcp = sec23_r2->Parameter(0);
467  double sec23_r2_slope = sec23_r2->Parameter(1);
468  TFitResultPtr sec23_r3 = gr23->Fit("pol1","S+","",41,54);
469  double sec23_r3_intcp = sec23_r3->Parameter(0);
470  double sec23_r3_slope = sec23_r3->Parameter(1);
471  TFitResultPtr sec23_r4 = gr23->Fit("pol1","S+","",54,60);
472  double sec23_r4_intcp = sec23_r4->Parameter(0);
473  double sec23_r4_slope = sec23_r4->Parameter(1);
474  c4->cd(6);
475  ntuple.Draw("z:z_err:sigma:sigma_err","sector >23 & sector <25 " );
476  TGraphErrors *gr24 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
477  gr24->Draw("ap"); //draw graph in current pad
478  gr24->SetTitle("Sector 24");
479  gr24->SetMaximum(0.4);
480  gr24->SetMinimum(0.1);
481  TFitResultPtr sec24_r1 = gr24->Fit("pol1","S","",10,30);
482  double sec24_r1_intcp = sec24_r1->Parameter(0);
483  double sec24_r1_slope = sec24_r1->Parameter(1);
484  TFitResultPtr sec24_r2 = gr24->Fit("pol1","S+","",30,41);
485  double sec24_r2_intcp = sec24_r2->Parameter(0);
486  double sec24_r2_slope = sec24_r2->Parameter(1);
487  TFitResultPtr sec24_r3 = gr24->Fit("pol1","S+","",41,54);
488  double sec24_r3_intcp = sec24_r3->Parameter(0);
489  double sec24_r3_slope = sec24_r3->Parameter(1);
490  TFitResultPtr sec24_r4 = gr24->Fit("pol1","S+","",54,60);
491  double sec24_r4_intcp = sec24_r4->Parameter(0);
492  double sec24_r4_slope = sec24_r4->Parameter(1);
493  auto c5 = new TCanvas("c5","c5",600,500);
494  c5->Divide(3,2);
495  c5->cd(1);
496  ntuple.Draw("z:z_err:sigma:sigma_err","sector >24 & sector <26 " );
497  TGraphErrors *gr25 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
498  gr25->Draw("ap"); //draw graph in current pad
499  gr25->SetTitle("Sector 25");
500  gr25->SetMaximum(0.4);
501  gr25->SetMinimum(0.1);
502  TFitResultPtr sec25_r1 = gr25->Fit("pol1","S","",10,30);
503  double sec25_r1_intcp = sec25_r1->Parameter(0);
504  double sec25_r1_slope = sec25_r1->Parameter(1);
505  TFitResultPtr sec25_r2 = gr25->Fit("pol1","S+","",30,41);
506  double sec25_r2_intcp = sec25_r2->Parameter(0);
507  double sec25_r2_slope = sec25_r2->Parameter(1);
508  TFitResultPtr sec25_r3 = gr25->Fit("pol1","S+","",41,54);
509  double sec25_r3_intcp = sec25_r3->Parameter(0);
510  double sec25_r3_slope = sec25_r3->Parameter(1);
511  TFitResultPtr sec25_r4 = gr25->Fit("pol1","S+","",54,60);
512  double sec25_r4_intcp = sec25_r4->Parameter(0);
513  double sec25_r4_slope = sec25_r4->Parameter(1);
514  c5->cd(2);
515  ntuple.Draw("z:z_err:sigma:sigma_err","sector >25 & sector <27 " );
516  TGraphErrors *gr26 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
517  gr26->Draw("ap"); //draw graph in current pad
518  gr26->SetTitle("Sector 26");
519  gr26->SetMaximum(0.4);
520  gr26->SetMinimum(0.1);
521  TFitResultPtr sec26_r1 = gr26->Fit("pol1","S","",10,30);
522  double sec26_r1_intcp = sec26_r1->Parameter(0);
523  double sec26_r1_slope = sec26_r1->Parameter(1);
524  TFitResultPtr sec26_r2 = gr26->Fit("pol1","S+","",30,41);
525  double sec26_r2_intcp = sec26_r2->Parameter(0);
526  double sec26_r2_slope = sec26_r2->Parameter(1);
527  TFitResultPtr sec26_r3 = gr26->Fit("pol1","S+","",41,54);
528  double sec26_r3_intcp = sec26_r3->Parameter(0);
529  double sec26_r3_slope = sec26_r3->Parameter(1);
530  TFitResultPtr sec26_r4 = gr26->Fit("pol1","S+","",54,60);
531  double sec26_r4_intcp = sec26_r4->Parameter(0);
532  double sec26_r4_slope = sec26_r4->Parameter(1);
533 c5->cd(3);
534  ntuple.Draw("z:z_err:sigma:sigma_err","sector >26 & sector <28 " );
535  TGraphErrors *gr27 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
536  gr27->Draw("ap"); //draw graph in current pad
537  gr27->SetTitle("Sector 27");
538  gr27->SetMaximum(0.4);
539  gr27->SetMinimum(0.1);
540  TFitResultPtr sec27_r1 = gr27->Fit("pol1","S","",10,30);
541  double sec27_r1_intcp = sec27_r1->Parameter(0);
542  double sec27_r1_slope = sec27_r1->Parameter(1);
543  TFitResultPtr sec27_r2 = gr27->Fit("pol1","S+","",30,41);
544  double sec27_r2_intcp = sec27_r2->Parameter(0);
545  double sec27_r2_slope = sec27_r2->Parameter(1);
546  TFitResultPtr sec27_r3 = gr27->Fit("pol1","S+","",41,54);
547  double sec27_r3_intcp = sec27_r3->Parameter(0);
548  double sec27_r3_slope = sec27_r3->Parameter(1);
549  TFitResultPtr sec27_r4 = gr27->Fit("pol1","S+","",54,60);
550  double sec27_r4_intcp = sec27_r4->Parameter(0);
551  double sec27_r4_slope = sec27_r4->Parameter(1);
552  c5->cd(4);
553  ntuple.Draw("z:z_err:sigma:sigma_err","sector >27 & sector <29 " );
554  TGraphErrors *gr28 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
555  gr28->Draw("ap"); //draw graph in current pad
556  gr28->SetTitle("Sector 28");
557  gr28->SetMaximum(0.4);
558  gr28->SetMinimum(0.1);
559  TFitResultPtr sec28_r1 = gr28->Fit("pol1","S","",10,30);
560  double sec28_r1_intcp = sec28_r1->Parameter(0);
561  double sec28_r1_slope = sec28_r1->Parameter(1);
562  TFitResultPtr sec28_r2 = gr28->Fit("pol1","S+","",30,41);
563  double sec28_r2_intcp = sec28_r2->Parameter(0);
564  double sec28_r2_slope = sec28_r2->Parameter(1);
565  TFitResultPtr sec28_r3 = gr28->Fit("pol1","S+","",41,54);
566  double sec28_r3_intcp = sec28_r3->Parameter(0);
567  double sec28_r3_slope = sec28_r3->Parameter(1);
568  TFitResultPtr sec28_r4 = gr28->Fit("pol1","S+","",54,60);
569  double sec28_r4_intcp = sec28_r4->Parameter(0);
570  double sec28_r4_slope = sec28_r4->Parameter(1);
571  c5->cd(5);
572  ntuple.Draw("z:z_err:sigma:sigma_err","sector >28 & sector <30 " );
573  TGraphErrors *gr29 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
574  gr29->Draw("ap"); //draw graph in current pad
575  gr29->SetTitle("Sector 29");
576  gr29->SetMaximum(0.4);
577  gr29->SetMinimum(0.1);
578  TFitResultPtr sec29_r1 = gr29->Fit("pol1","S","",10,30);
579  double sec29_r1_intcp = sec29_r1->Parameter(0);
580  double sec29_r1_slope = sec29_r1->Parameter(1);
581  TFitResultPtr sec29_r2 = gr29->Fit("pol1","S+","",30,41);
582  double sec29_r2_intcp = sec29_r2->Parameter(0);
583  double sec29_r2_slope = sec29_r2->Parameter(1);
584  TFitResultPtr sec29_r3 = gr29->Fit("pol1","S+","",41,54);
585  double sec29_r3_intcp = sec29_r3->Parameter(0);
586  double sec29_r3_slope = sec29_r3->Parameter(1);
587  TFitResultPtr sec29_r4 = gr29->Fit("pol1","S+","",54,60);
588  double sec29_r4_intcp = sec29_r4->Parameter(0);
589  double sec29_r4_slope = sec29_r4->Parameter(1);
590  c5->cd(6);
591  ntuple.Draw("z:z_err:sigma:sigma_err","sector >29 & sector <31 " );
592  TGraphErrors *gr30 = new TGraphErrors(ntuple.GetSelectedRows(), ntuple.GetV1(), ntuple.GetV3(),ntuple.GetV2(), ntuple.GetV4());
593  gr30->Draw("ap"); //draw graph in current pad
594  gr30->SetTitle("Sector 30");
595  gr30->SetMaximum(0.4);
596  gr30->SetMinimum(0.1);
597  TFitResultPtr sec30_r1 = gr30->Fit("pol1","S","",10,30);
598  double sec30_r1_intcp = sec30_r1->Parameter(0);
599  double sec30_r1_slope = sec30_r1->Parameter(1);
600  TFitResultPtr sec30_r2 = gr30->Fit("pol1","S+","",30,41);
601  double sec30_r2_intcp = sec30_r2->Parameter(0);
602  double sec30_r2_slope = sec30_r2->Parameter(1);
603  TFitResultPtr sec30_r3 = gr30->Fit("pol1","S+","",41,54);
604  double sec30_r3_intcp = sec30_r3->Parameter(0);
605  double sec30_r3_slope = sec30_r3->Parameter(1);
606  TFitResultPtr sec30_r4 = gr30->Fit("pol1","S+","",54,60);
607  double sec30_r4_intcp = sec30_r4->Parameter(0);
608  double sec30_r4_slope = sec30_r4->Parameter(1);
609  in.close();
610  f->Write();
611  SC_PTC << sec1_r1_intcp << "\t" << sec1_r1_slope << "\t" << sec1_r2_intcp << "\t" << sec1_r2_slope << "\t" <<sec1_r3_intcp << "\t" << sec1_r3_slope << sec1_r4_intcp << "\t" << sec1_r4_slope <<endl;
612  SC_PTC << sec2_r1_intcp << "\t" << sec2_r1_slope << "\t" << sec2_r2_intcp << "\t" << sec2_r2_slope << "\t" <<sec2_r3_intcp << "\t" << sec2_r3_slope <<sec2_r4_intcp << "\t" << sec2_r4_slope <<endl;
613  SC_PTC << sec3_r1_intcp << "\t" << sec3_r1_slope << "\t" << sec3_r2_intcp << "\t" << sec3_r2_slope << "\t" <<sec3_r3_intcp << "\t" << sec3_r3_slope <<sec3_r4_intcp << "\t" << sec3_r4_slope <<endl;
614  SC_PTC << sec4_r1_intcp << "\t" << sec4_r1_slope << "\t" << sec4_r2_intcp << "\t" << sec4_r2_slope << "\t" <<sec4_r3_intcp << "\t" << sec4_r3_slope <<sec4_r4_intcp << "\t" << sec4_r4_slope <<endl;
615  SC_PTC << sec5_r1_intcp << "\t" << sec5_r1_slope << "\t" << sec5_r2_intcp << "\t" << sec5_r2_slope << "\t" <<sec5_r3_intcp << "\t" << sec5_r3_slope <<sec5_r4_intcp << "\t" << sec5_r4_slope <<endl;
616  SC_PTC << sec6_r1_intcp << "\t" << sec6_r1_slope << "\t" << sec6_r2_intcp << "\t" << sec6_r2_slope << "\t" <<sec6_r3_intcp << "\t" << sec6_r3_slope <<sec6_r4_intcp << "\t" << sec6_r4_slope <<endl;
617  SC_PTC << sec7_r1_intcp << "\t" << sec7_r1_slope << "\t" << sec7_r2_intcp << "\t" << sec7_r2_slope << "\t" <<sec7_r3_intcp << "\t" << sec7_r3_slope <<sec7_r4_intcp << "\t" << sec7_r4_slope <<endl;
618  SC_PTC << sec8_r1_intcp << "\t" << sec8_r1_slope << "\t" << sec8_r2_intcp << "\t" << sec8_r2_slope << "\t" <<sec8_r3_intcp << "\t" << sec8_r3_slope <<sec8_r4_intcp << "\t" << sec8_r4_slope <<endl;
619  SC_PTC << sec9_r1_intcp << "\t" << sec9_r1_slope << "\t" << sec9_r2_intcp << "\t" << sec9_r2_slope << "\t" <<sec9_r3_intcp << "\t" << sec9_r3_slope <<sec9_r4_intcp << "\t" << sec9_r4_slope <<endl;
620 
621  SC_PTC << sec10_r1_intcp << "\t" << sec10_r1_slope << "\t" << sec10_r2_intcp << "\t" << sec10_r2_slope << "\t" <<sec10_r3_intcp << "\t" << sec10_r3_slope <<sec10_r4_intcp << "\t" << sec10_r4_slope <<endl;
622  SC_PTC << sec11_r1_intcp << "\t" << sec11_r1_slope << "\t" << sec11_r2_intcp << "\t" << sec11_r2_slope << "\t" <<sec11_r3_intcp << "\t" << sec11_r3_slope <<sec11_r4_intcp << "\t" << sec11_r4_slope <<endl;
623  SC_PTC << sec12_r1_intcp << "\t" << sec12_r1_slope << "\t" << sec12_r2_intcp << "\t" << sec12_r2_slope << "\t" <<sec12_r3_intcp << "\t" << sec12_r3_slope <<sec12_r4_intcp << "\t" << sec12_r4_slope <<endl;
624  SC_PTC << sec13_r1_intcp << "\t" << sec13_r1_slope << "\t" << sec13_r2_intcp << "\t" << sec13_r2_slope << "\t" <<sec13_r3_intcp << "\t" << sec13_r3_slope << sec13_r4_intcp << "\t" << sec13_r4_slope <<endl;
625  SC_PTC << sec14_r1_intcp << "\t" << sec14_r1_slope << "\t" << sec14_r2_intcp << "\t" << sec14_r2_slope << "\t" <<sec14_r3_intcp << "\t" << sec14_r3_slope <<sec14_r4_intcp << "\t" << sec14_r4_slope <<endl;
626  SC_PTC << sec15_r1_intcp << "\t" << sec15_r1_slope << "\t" << sec15_r2_intcp << "\t" << sec15_r2_slope << "\t" <<sec15_r3_intcp << "\t" << sec15_r3_slope<<sec15_r4_intcp << "\t" << sec15_r4_slope <<endl;
627  SC_PTC << sec16_r1_intcp << "\t" << sec16_r1_slope << "\t" << sec16_r2_intcp << "\t" << sec16_r2_slope << "\t" <<sec16_r3_intcp << "\t" << sec16_r3_slope<<sec16_r4_intcp << "\t" << sec16_r4_slope <<endl;
628  SC_PTC << sec17_r1_intcp << "\t" << sec17_r1_slope << "\t" << sec17_r2_intcp << "\t" << sec17_r2_slope << "\t" <<sec17_r3_intcp << "\t" << sec17_r3_slope<<sec17_r4_intcp << "\t" << sec17_r4_slope <<endl;
629  SC_PTC << sec18_r1_intcp << "\t" << sec18_r1_slope << "\t" << sec18_r2_intcp << "\t" << sec18_r2_slope << "\t" <<sec18_r3_intcp << "\t" << sec18_r3_slope<<sec18_r4_intcp << "\t" << sec18_r4_slope <<endl;
630  SC_PTC << sec19_r1_intcp << "\t" << sec19_r1_slope << "\t" << sec19_r2_intcp << "\t" << sec19_r2_slope << "\t" <<sec19_r3_intcp << "\t" << sec19_r3_slope<<sec19_r4_intcp << "\t" << sec19_r4_slope <<endl;
631  SC_PTC << sec20_r1_intcp << "\t" << sec20_r1_slope << "\t" << sec20_r2_intcp << "\t" << sec20_r2_slope << "\t" <<sec20_r3_intcp << "\t" << sec20_r3_slope<<sec20_r4_intcp << "\t" << sec20_r4_slope <<endl;
632  SC_PTC << sec21_r1_intcp << "\t" << sec21_r1_slope << "\t" << sec21_r2_intcp << "\t" << sec21_r2_slope << "\t" <<sec21_r3_intcp << "\t" << sec21_r3_slope<<sec21_r4_intcp << "\t" << sec21_r4_slope <<endl;
633  SC_PTC << sec22_r1_intcp << "\t" << sec22_r1_slope << "\t" << sec22_r2_intcp << "\t" << sec22_r2_slope << "\t" <<sec22_r3_intcp << "\t" << sec22_r3_slope<<sec22_r4_intcp << "\t" << sec22_r4_slope <<endl;
634  SC_PTC << sec23_r1_intcp << "\t" << sec23_r1_slope << "\t" << sec23_r2_intcp << "\t" << sec23_r2_slope << "\t" <<sec23_r3_intcp << "\t" << sec23_r3_slope<<sec23_r4_intcp << "\t" << sec23_r4_slope <<endl;
635  SC_PTC << sec24_r1_intcp << "\t" << sec24_r1_slope << "\t" << sec24_r2_intcp << "\t" << sec24_r2_slope << "\t" <<sec24_r3_intcp << "\t" << sec24_r3_slope<<sec24_r4_intcp << "\t" << sec24_r4_slope <<endl;
636  SC_PTC << sec25_r1_intcp << "\t" << sec25_r1_slope << "\t" << sec25_r2_intcp << "\t" << sec25_r2_slope << "\t" <<sec25_r3_intcp << "\t" << sec25_r3_slope<<sec25_r4_intcp << "\t" << sec25_r4_slope <<endl;
637  SC_PTC << sec26_r1_intcp << "\t" << sec26_r1_slope << "\t" << sec26_r2_intcp << "\t" << sec26_r2_slope << "\t" <<sec26_r3_intcp << "\t" << sec26_r3_slope<<sec26_r4_intcp << "\t" << sec26_r4_slope <<endl;
638  SC_PTC << sec27_r1_intcp << "\t" << sec27_r1_slope << "\t" << sec27_r2_intcp << "\t" << sec27_r2_slope << "\t" <<sec27_r3_intcp << "\t" << sec27_r3_slope<<sec27_r4_intcp << "\t" << sec27_r4_slope <<endl;
639  SC_PTC << sec28_r1_intcp << "\t" << sec28_r1_slope << "\t" << sec28_r2_intcp << "\t" << sec28_r2_slope << "\t" <<sec28_r3_intcp << "\t" << sec28_r3_slope<<sec28_r4_intcp << "\t" << sec28_r4_slope <<endl;
640  SC_PTC << sec29_r1_intcp << "\t" << sec29_r1_slope << "\t" << sec29_r2_intcp << "\t" << sec29_r2_slope << "\t" <<sec29_r3_intcp << "\t" << sec29_r3_slope<<sec29_r4_intcp << "\t" << sec29_r4_slope <<endl;
641  SC_PTC << sec30_r1_intcp << "\t" << sec30_r1_slope << "\t" << sec30_r2_intcp << "\t" << sec30_r2_slope << "\t" <<sec30_r3_intcp << "\t" << sec30_r3_slope<<sec30_r4_intcp << "\t" << sec30_r4_slope <<endl;
642  SC_PTC.close();
643 }
SC_PTC
void SC_PTC(char const *input_filename)
Definition: SC_PTC.C:89
gr3
static TGraph * gr3
Definition: JEventProcessor_BCAL_TimeCalibration.cc:107
gr
TGraph * gr[NCHANNELS]
Definition: st_tw_resols.C:43
c1
Double_t c1[2][NMODULES]
Definition: tw_corr.C:68
f
TF1 * f
Definition: FitGains.C:21
gr12
static TGraphErrors * gr12
Definition: JEventProcessor_BCAL_TimeCalibration.cc:103
c2
Double_t c2[2][NMODULES]
Definition: tw_corr.C:69
gr2
static TGraph * gr2
Definition: JEventProcessor_BCAL_TimeCalibration.cc:106
SC_PTC_ResolvsZPlot
void SC_PTC_ResolvsZPlot()
Definition: ST_Propagation_Time/macros/SC_PTC_ResolvsZPlot.C:8
sigma
Double_t sigma[NCHANNELS]
Definition: st_tw_resols.C:37
PT_can
TCanvas * PT_can[30]
Definition: ST_Propagation_Time/macros/Resolution.C:74
gr11
static TGraphErrors * gr11
Definition: JEventProcessor_BCAL_TimeCalibration.cc:101
c4
TCanvas * c4
Definition: plot_epics_dump.C:170
c5
TCanvas * c5
Definition: ST_Monitoring_Expert_Waveform.C:185
c3
TCanvas * c3
Definition: plot_epics_dump.C:151
printf
printf("string=%s", string)