===============================================================================================================================================================================================
variant2: for each run define phi-averaged A for normalization channel (Dref,16) and then, divide Rijk on it, i.e. get RRijk
41 printf(
"reco: gROOT Reset \n");
43 gROOT->SetStyle(
"Plain");
50 gStyle->SetOptStat(101110);
55 gStyle->SetOptFit(0010);
64 gStyle->SetStatX(0.91);
65 gStyle->SetStatY(0.75);
66 gStyle->SetStatW(0.20);
67 gStyle->SetStatH(0.10);
73 Float_t LeftOffset = 0.12;
74 Float_t TopOffset = 0.22;
76 gStyle->SetLineWidth(1);
80 gStyle->SetOptTitle(0);
82 gStyle->SetTitleColor(1);
84 gStyle->SetTitleX(0.15);
85 gStyle->SetTitleH(0.15);
87 gStyle->SetTitleW(0.60);
88 gStyle->SetTitleFont(42);
89 gStyle->SetTitleFontSize(0.07);
91 gStyle->SetPalette(1);
104 gStyle->SetPadTopMargin(TopOffset);
105 gStyle->SetPadBottomMargin(LeftOffset);
106 gStyle->SetPadRightMargin(TopOffset);
107 gStyle->SetPadLeftMargin(LeftOffset);
155 TFile *hfile1 =
new TFile(
"Global_346300.root",
"READ");
167 TPostScript psfile(
"zadcampltest.ps", 111);
171 TCanvas *
c1 =
new TCanvas(
"c1",
"Hcal4test", 200, 10, 700, 900);
174 TDirectory *
dir = (TDirectory *)hfile1->FindObjectAny(dirnm.c_str());
186 TH2F *Zzzdepth1hef1 = (TH2F *)dir->FindObjectAny(
"h_mapDepth1_HE");
190 Zzzdepth1hef1->SetXTitle(
"#eta \b");
191 Zzzdepth1hef1->SetYTitle(
"#phi \b");
192 Zzzdepth1hef1->SetZTitle(
"h_mapDepth1_HE \b");
193 Zzzdepth1hef1->Draw(
"COLZ");
196 TH2F *Zzzdepth1hef2 = (TH2F *)dir->FindObjectAny(
"h_mapDepth2_HE");
200 Zzzdepth1hef2->SetXTitle(
"#eta \b");
201 Zzzdepth1hef2->SetYTitle(
"#phi \b");
202 Zzzdepth1hef2->SetZTitle(
"h_mapDepth2_HE \b");
203 Zzzdepth1hef2->Draw(
"COLZ");
206 TH2F *Zzzdepth1hef3 = (TH2F *)dir->FindObjectAny(
"h_mapDepth3_HE");
210 Zzzdepth1hef3->SetXTitle(
"#eta \b");
211 Zzzdepth1hef3->SetYTitle(
"#phi \b");
212 Zzzdepth1hef3->SetZTitle(
"h_mapDepth3_HE \b");
213 Zzzdepth1hef3->Draw(
"COLZ");
216 TH2F *Zzzdepth1hef4 = (TH2F *)dir->FindObjectAny(
"h_mapDepth4_HE");
220 Zzzdepth1hef4->SetXTitle(
"#eta \b");
221 Zzzdepth1hef4->SetYTitle(
"#phi \b");
222 Zzzdepth1hef4->SetZTitle(
"h_mapDepth4_HE \b");
223 Zzzdepth1hef4->Draw(
"COLZ");
226 TH2F *Zzzdepth1hef5 = (TH2F *)dir->FindObjectAny(
"h_mapDepth5_HE");
230 Zzzdepth1hef5->SetXTitle(
"#eta \b");
231 Zzzdepth1hef5->SetYTitle(
"#phi \b");
232 Zzzdepth1hef5->SetZTitle(
"h_mapDepth5_HE \b");
233 Zzzdepth1hef5->Draw(
"COLZ");
236 TH2F *Zzzdepth1hef6 = (TH2F *)dir->FindObjectAny(
"h_mapDepth6_HE");
240 Zzzdepth1hef6->SetXTitle(
"#eta \b");
241 Zzzdepth1hef6->SetYTitle(
"#phi \b");
242 Zzzdepth1hef6->SetZTitle(
"h_mapDepth6_HE \b");
243 Zzzdepth1hef6->Draw(
"COLZ");
246 TH2F *Zzzdepth1hef7 = (TH2F *)dir->FindObjectAny(
"h_mapDepth7_HE");
250 Zzzdepth1hef7->SetXTitle(
"#eta \b");
251 Zzzdepth1hef7->SetYTitle(
"#phi \b");
252 Zzzdepth1hef7->SetZTitle(
"h_mapDepth7_HE \b");
253 Zzzdepth1hef7->Draw(
"COLZ");
256 TH1F *jdjswi8 = (TH1F *)dir->FindObjectAny(
"h_ADCAmpl_HE");
257 gPad->SetLogy(kFALSE);
259 jdjswi8->SetMarkerStyle(20);
260 jdjswi8->SetMarkerSize(0.8);
261 jdjswi8->GetYaxis()->SetLabelSize(0.04);
262 jdjswi8->SetXTitle(
"h_ADCAmpl_HE \b");
263 jdjswi8->SetMarkerColor(2);
264 jdjswi8->SetLineColor(2);
278 TH1F *wqwhef1 = (TH1F *)dir->FindObjectAny(
"h_mapDepth1ADCAmpl_HE");
282 wqwhef1->SetZTitle(
"h_mapDepth1ADCAmpl_HE \b");
283 wqwhef1->Draw(
"COLZ");
286 TH1F *wqwhef2 = (TH1F *)dir->FindObjectAny(
"h_mapDepth2ADCAmpl_HE");
290 wqwhef2->SetZTitle(
"h_mapDepth2ADCAmpl_HE \b");
291 wqwhef2->Draw(
"COLZ");
294 TH1F *wqwhef3 = (TH1F *)dir->FindObjectAny(
"h_mapDepth3ADCAmpl_HE");
298 wqwhef3->SetZTitle(
"h_mapDepth3ADCAmpl_HE \b");
299 wqwhef3->Draw(
"COLZ");
302 TH1F *wqwhef4 = (TH1F *)dir->FindObjectAny(
"h_mapDepth4ADCAmpl_HE");
306 wqwhef4->SetZTitle(
"h_mapDepth4ADCAmpl_HE \b");
307 wqwhef4->Draw(
"COLZ");
310 TH1F *wqwhef5 = (TH1F *)dir->FindObjectAny(
"h_mapDepth5ADCAmpl_HE");
314 wqwhef5->SetZTitle(
"h_mapDepth5ADCAmpl_HE \b");
315 wqwhef5->Draw(
"COLZ");
318 TH1F *wqwhef6 = (TH1F *)dir->FindObjectAny(
"h_mapDepth6ADCAmpl_HE");
322 wqwhef6->SetZTitle(
"h_mapDepth6ADCAmpl_HE \b");
323 wqwhef6->Draw(
"COLZ");
326 TH1F *wqwhef7 = (TH1F *)dir->FindObjectAny(
"h_mapDepth7ADCAmpl_HE");
330 wqwhef7->SetZTitle(
"h_mapDepth7ADCAmpl_HE \b");
331 wqwhef7->Draw(
"COLZ");
344 TH1F *azzdepth1hef1 = (TH1F *)dir->FindObjectAny(
"h_sumADCAmplperLS3");
347 azzdepth1hef1->SetZTitle(
"h_sumADCAmplperLS3 \b");
348 azzdepth1hef1->Draw(
"");
351 TH1F *azzdepth1hef2 = (TH1F *)dir->FindObjectAny(
"h_sumADCAmplperLS4");
354 azzdepth1hef2->SetZTitle(
"h_sumADCAmplperLS3 \b");
355 azzdepth1hef2->Draw(
"");
358 TH1F *azzdepth1hef3 = (TH1F *)dir->FindObjectAny(
"h_sumADCAmplperLS5");
361 azzdepth1hef3->SetZTitle(
"h_sumADCAmplperLS3 \b");
362 azzdepth1hef3->Draw(
"");
365 TH1F *azzdepth1hef4 = (TH1F *)dir->FindObjectAny(
"h_sumADCAmplperLSdepth4HEu");
368 azzdepth1hef4->SetZTitle(
"h_sumADCAmplperLS3 \b");
369 azzdepth1hef4->Draw(
"");
372 TH1F *azzdepth1hef5 = (TH1F *)dir->FindObjectAny(
"h_sumADCAmplperLSdepth5HEu");
375 azzdepth1hef5->SetZTitle(
"h_sumADCAmplperLS3 \b");
376 azzdepth1hef5->Draw(
"");
379 TH1F *azzdepth1hef6 = (TH1F *)dir->FindObjectAny(
"h_sumADCAmplperLSdepth6HEu");
382 azzdepth1hef6->SetZTitle(
"h_sumADCAmplperLS3 \b");
383 azzdepth1hef6->Draw(
"");
386 TH1F *azzdepth1hef7 = (TH1F *)dir->FindObjectAny(
"h_sumADCAmplperLSdepth7HEu");
389 azzdepth1hef7->SetZTitle(
"h_sumADCAmplperLS3 \b");
390 azzdepth1hef7->Draw(
"");
403 TH1F *poqw1 = (TH1F *)dir->FindObjectAny(
"h_sumCutADCAmplperLS3");
406 poqw1->SetZTitle(
"h_sumCutADCAmplperLS3 \b");
410 TH1F *poqw2 = (TH1F *)dir->FindObjectAny(
"h_sumCutADCAmplperLS4");
413 poqw2->SetZTitle(
"h_sumCutADCAmplperLS3 \b");
417 TH1F *poqw3 = (TH1F *)dir->FindObjectAny(
"h_sumCutADCAmplperLS5");
420 poqw3->SetZTitle(
"h_sumCutADCAmplperLS3 \b");
424 TH1F *poqw4 = (TH1F *)dir->FindObjectAny(
"h_sumCutADCAmplperLSdepth4HEu");
427 poqw4->SetZTitle(
"h_sumCutADCAmplperLS3 \b");
431 TH1F *poqw5 = (TH1F *)dir->FindObjectAny(
"h_sumCutADCAmplperLSdepth5HEu");
434 poqw5->SetZTitle(
"h_sumCutADCAmplperLS3 \b");
438 TH1F *poqw6 = (TH1F *)dir->FindObjectAny(
"h_sumCutADCAmplperLSdepth6HEu");
441 poqw6->SetZTitle(
"h_sumCutADCAmplperLS3 \b");
445 TH1F *poqw7 = (TH1F *)dir->FindObjectAny(
"h_sumCutADCAmplperLSdepth7HEu");
448 poqw7->SetZTitle(
"h_sumCutADCAmplperLS3 \b");
462 TH1F *cdew1 = (TH1F *)dir->FindObjectAny(
"h_sum0ADCAmplperLS3");
465 cdew1->SetZTitle(
"h_sum0ADCAmplperLS3 \b");
469 TH1F *cdew2 = (TH1F *)dir->FindObjectAny(
"h_sum0ADCAmplperLS4");
472 cdew2->SetZTitle(
"h_sum0ADCAmplperLS3 \b");
476 TH1F *cdew3 = (TH1F *)dir->FindObjectAny(
"h_sum0ADCAmplperLS5");
479 cdew3->SetZTitle(
"h_sum0ADCAmplperLS3 \b");
483 TH1F *cdew4 = (TH1F *)dir->FindObjectAny(
"h_sum0ADCAmplperLSdepth4HEu");
486 cdew4->SetZTitle(
"h_sum0ADCAmplperLS3 \b");
490 TH1F *cdew5 = (TH1F *)dir->FindObjectAny(
"h_sum0ADCAmplperLSdepth5HEu");
493 cdew5->SetZTitle(
"h_sum0ADCAmplperLS3 \b");
497 TH1F *cdew6 = (TH1F *)dir->FindObjectAny(
"h_sum0ADCAmplperLSdepth6HEu");
500 cdew6->SetZTitle(
"h_sum0ADCAmplperLS3 \b");
504 TH1F *cdew7 = (TH1F *)dir->FindObjectAny(
"h_sum0ADCAmplperLSdepth7HEu");
507 cdew7->SetZTitle(
"h_sum0ADCAmplperLS3 \b");
printf("params %d %f %f %f\n", minT, eps, errmax, chi2max)