source: trunk/MagicSoft/Mars/macros/calibration.C@ 3048

Last change on this file since 3048 was 3039, checked in by gaug, 22 years ago
*** empty log message ***
File size: 18.1 KB
Line 
1/* ======================================================================== *\
2!
3! *
4! * This file is part of MARS, the MAGIC Analysis and Reconstruction
5! * Software. It is distributed to you in the hope that it can be a useful
6! * and timesaving tool in analysing Data of imaging Cerenkov telescopes.
7! * It is distributed WITHOUT ANY WARRANTY.
8! *
9! * Permission to use, copy, modify and distribute this software and its
10! * documentation for any purpose is hereby granted without fee,
11! * provided that the above copyright notice appear in all copies and
12! * that both that copyright notice and this permission notice appear
13! * in supporting documentation. It is provided "as is" without express
14! * or implied warranty.
15! *
16!
17!
18! Author(s): Markus Gaug, 11/2003 <mailto:markus@ifae.es>
19!
20! Copyright: MAGIC Software Development, 2000-2003
21!
22!
23\* ======================================================================== */
24
25//const TString pedfile = "/mnt/users/mdoro/Mars/Data/20040201_14418_P_OffMrk421-1_E.root";
26//const TString calfile = "/mnt/users/mdoro/Mars/Data/20040201_1441*_C_OffMrk421-1_E.root";
27
28const TString pedfile = "/mnt/Data/rootdata/CrabNebula/2004_01_27/20040126_12386_P_Cab-On_E.root";
29const TString calfile = "/mnt/Data/rootdata/CrabNebula/2004_01_27/20040126_12525_C_Cab-On_E.root";
30
31
32//const TString pedfile = "/mnt/Data/rootdata/Miscellaneous/2003_12_19/20031218_03522_P_Park_E.root";
33//const TString calfile = "/mnt/Data/rootdata/Miscellaneous/2003_12_19/20031218_03527_C_Park_E.root";
34
35void calibration(TString pedname=pedfile,
36 TString calname=calfile)
37{
38
39 //
40 // Create a empty Parameter List and an empty Task List
41 // The tasklist is identified in the eventloop by its name
42 //
43 MParList plist;
44
45 MTaskList tlist;
46 plist.AddToList(&tlist);
47
48 //
49 // Now setup the tasks and tasklist for the pedestals:
50 // ---------------------------------------------------
51 //
52
53 MReadMarsFile read("Events", pedname);
54 read.DisableAutoScheme();
55
56 MGeomApply geomapl;
57 MExtractSignal sigcalc;
58 //
59 // We saw that the signal jumps between slices,
60 // thus set the extraction range as high as possible
61 //
62 sigcalc.SetRange(1,14,1,14);
63
64 MPedCalcPedRun pedcalc;
65 MFillH fill("MPedestalCam", "MExtractedSignalCam");
66
67 tlist.AddToList(&read);
68 tlist.AddToList(&geomapl);
69 tlist.AddToList(&sigcalc);
70 tlist.AddToList(&pedcalc);
71 tlist.AddToList(&fill);
72
73 MGeomCamMagic geomcam;
74 MPedestalCam pedcam;
75 plist.AddToList(&geomcam);
76 plist.AddToList(&pedcam);
77
78 //
79 // Create and setup the eventloop
80 //
81 MEvtLoop evtloop;
82 evtloop.SetParList(&plist);
83
84 //
85 // Execute first analysis
86 //
87 if (!evtloop.Eventloop())
88 return;
89
90 tlist.PrintStatistics();
91
92 // pedcam(17).DrawClone();
93
94 MStatusDisplay *d1 = new MStatusDisplay;
95 d1->SetUpdateTime(3000);
96 d1->Resize(850,700);
97
98 MHCamera dispped0 (geomcam, "MPedestalPix;Pedestal", "Mean per Slice");
99 MHCamera dispped1 (geomcam, "MPedestalPix;PedestalErr", "Mean Error per Slice");
100 MHCamera dispped2 (geomcam, "MPedestalPix;PedestalRms", "RMS per Slice");
101 MHCamera dispped3 (geomcam, "MPedestalPix;PedestalRmsErr", "RMS Error per Slice");
102 MHCamera dispped4 (geomcam, "MPedestalPix;Mean", "Fitted Mean per Slice");
103 MHCamera dispped5 (geomcam, "MPedestalPix;MeanErr", "Fitted Error of Mean per Slice");
104 MHCamera dispped6 (geomcam, "MPedestalPix;Sigma", "Fitted Sigma per Slice");
105 MHCamera dispped7 (geomcam, "MPedestalPix;SigmaErr", "Fitted Error of Sigma per Slice");
106 MHCamera dispped8 (geomcam, "MPedestalPix;Prob", "Probability of Fit");
107 MHCamera dispped9 (geomcam, "MPedestalPix;DeltaPedestalMean", "Rel. Diff. Mean per Slice (Calc.-Fitte)");
108 MHCamera dispped11 (geomcam, "MPedestalPix;DeltaPedestalMeanError", "Rel. Diff. Mean Error per Slice (Calc.-Fitted)");
109 MHCamera dispped12 (geomcam, "MPedestalPix;DeltaRmsSigma", "Rel. Diff. RMS per Slice (Calc.-Fitted)");
110 MHCamera dispped14 (geomcam, "MPedestalPix;DeltaRmsSigmaError", "Rel. Diff. RMS Error per Slice (Calc.-Fitted)");
111 dispped0.SetCamContent(pedcam, 0);
112 dispped0.SetCamError(pedcam, 1);
113 dispped1.SetCamContent(pedcam, 1);
114 dispped2.SetCamContent(pedcam, 2);
115 dispped2.SetCamError(pedcam,3);
116 dispped3.SetCamContent(pedcam, 3);
117 dispped4.SetCamContent(pedcam, 4);
118 dispped4.SetCamError(pedcam, 5);
119 dispped5.SetCamContent(pedcam, 5);
120 dispped6.SetCamContent(pedcam, 6);
121 dispped6.SetCamError(pedcam, 7);
122 dispped7.SetCamContent(pedcam, 7);
123 dispped8.SetCamContent(pedcam, 8);
124 dispped9.SetCamContent(pedcam, 9);
125 dispped9.SetCamError(pedcam, 10);
126 dispped11.SetCamContent(pedcam, 11);
127 dispped12.SetCamContent(pedcam, 12);
128 dispped12.SetCamError(pedcam, 13);
129 dispped14.SetCamContent(pedcam, 14);
130
131 dispped0.SetYTitle("Calc. Pedestal per slice [FADC counts]");
132 dispped1.SetYTitle("Calc. Pedestal Error per slice [FADC counts]");
133 dispped2.SetYTitle("Calc. Pedestal RMS per slice [FADC counts]");
134 dispped3.SetYTitle("Calc. Pedestal RMS Error per slice [FADC counts]");
135 dispped4.SetYTitle("Fitted Mean per slice [FADC counts]");
136 dispped5.SetYTitle("Error of Fitted Mean per slice [FADC counts]");
137 dispped6.SetYTitle("Fitted Sigma per slice [FADC counts]");
138 dispped7.SetYTitle("Error of Fitted Sigma per slice [FADC counts]");
139 dispped8.SetYTitle("Fit Probability [1]");
140 dispped9.SetYTitle("Rel. Diff. Pedestal Calc.-Fitted per slice [1]");
141 dispped11.SetYTitle("Rel. Diff. Pedestal Error Calc.-Fitted per slice [1]");
142 dispped12.SetYTitle("Rel. Diff. Pedestal RMS Calc.-Fitted per slice [1]");
143 dispped14.SetYTitle("Rel. Diff. Pedestal RMS Error Calc.-Fitted per slice [1]");
144
145 gStyle->SetOptStat(1111);
146 gStyle->SetOptFit();
147
148 // Histogram values
149 TCanvas &b1 = d1->AddTab("Direct Calculation");
150 b1.Divide(4,3);
151
152 CamDraw(b1,dispped0,pedcam,1,4,0);
153 CamDraw(b1,dispped1,pedcam,2,4,1);
154 CamDraw(b1,dispped2,pedcam,3,4,1);
155 CamDraw(b1,dispped3,pedcam,4,4,1);
156
157 // Fitted values
158 TCanvas &b2 = d1->AddTab("Fits");
159 b2.Divide(4,3);
160
161 CamDraw(b2,dispped4,pedcam,1,4,0);
162 CamDraw(b2,dispped5,pedcam,2,4,1);
163 CamDraw(b2,dispped6,pedcam,3,4,1);
164 CamDraw(b2,dispped7,pedcam,4,4,1);
165
166
167 // Fits Probability
168 TCanvas &b3 = d1->AddTab("Fit Probabilities");
169 b3.Divide(1,3);
170
171 CamDraw(b3,dispped8,pedcam,1,1,3);
172
173 // Differences
174 TCanvas &c4 = d1->AddTab("Relative Difference Calculation-Fits");
175 c4.Divide(4,3);
176
177 CamDraw(c4,dispped9,pedcam,1,4,1);
178 CamDraw(c4,dispped11,pedcam,2,4,1);
179 CamDraw(c4,dispped12,pedcam,3,4,1);
180 CamDraw(c4,dispped14,pedcam,4,4,1);
181
182 //
183 // Create a empty Parameter List and an empty Task List
184 //
185 MParList plist2;
186 MTaskList tlist2;
187 plist2.AddToList(&tlist2);
188
189 MExtractedSignalCam sigcam;
190 MCalibrationCam calcam;
191 //
192 // Get the previously created MPedestalCam into the new Parameter List
193 //
194 plist2.AddToList(&geomcam);
195 plist2.AddToList(&pedcam);
196 plist2.AddToList(&sigcam);
197 plist2.AddToList(&calcam);
198
199 //
200 // Get the MAGIC geometry
201 //
202 tlist2.AddToList(&geomapl);
203
204 //
205 // Now setup the new tasks and tasklist for the calibration
206 // ---------------------------------------------------
207 //
208
209 MReadMarsFile read2("Events", calname);
210 read2.DisableAutoScheme();
211
212 //
213 // We saw that the signal jumps between slices,
214 // thus take the sliding window
215 //
216 MExtractSignal2 sigcalc2;
217 MArrivalTimeCalc timecalc;
218 MCalibrationCalc calcalc;
219
220 //
221 // Skip the HiGain vs. LoGain calibration
222 //
223 calcalc.SkipHiLoGainCalibration();
224
225 //
226 // Making the step size a bit bigger, gives us
227 // faster results
228 //
229 // timecalc.SetStepSize(0.05);
230
231 //
232 // As long, as we don't have digital modules,
233 // we have to set the color of the pulser LED by hand
234 //
235 calcalc.SetPulserColor(MCalibrationCalc::kECT1);
236
237 //
238 // In case, we want to exclude a pre-defined list of bad pixels:
239 // (This is a preliminary feature)
240 //
241 //calcalc.ExcludePixelsFromAsciiFile("badpixels_all.dat");
242
243 //
244 // In case, you want to skip the blind pixel method:
245 // (NOT RECOMMENDED!!!)
246 //
247 // calcalc.SkipBlindPixelFit();
248
249 //
250 // In case, you want to skip the cosmics rejection
251 // (NOT RECOMMENDED!!!)
252 //
253 // calcalc.SkipCosmicsRejection();
254
255 //
256 // In case, you want to skip the quality checks
257 // (NOT RECOMMENDED!!!)
258 //
259 // calcalc.SkipQualityChecks();
260
261 //
262 // In case, we want to apply another fit function to the
263 // blind pixel
264 //
265 MCalibrationBlindPix *bp = calcam.GetBlindPixel();
266// bp->ChangeFitFunc(MHCalibrationBlindPixel::kEPolya);
267// bp->ChangeFitFunc(MHCalibrationBlindPixel::kEPoisson4);
268
269 tlist2.AddToList(&read2);
270 tlist2.AddToList(&sigcalc2);
271 //
272 // In case, you want to skip the somewhat lengthy calculation
273 // of the arrival times using a spline, uncomment the next line
274 //
275 tlist2.AddToList(&timecalc);
276 tlist2.AddToList(&calcalc);
277
278 //
279 // Create and setup the eventloop
280 //
281 MEvtLoop evtloop2;
282 evtloop2.SetParList(&plist2);
283
284 MProgressBar bar;
285 evtloop2.SetProgressBar(&bar);
286
287 //
288 // Execute second analysis
289 //
290 if (!evtloop2.Eventloop())
291 return;
292
293 tlist2.PrintStatistics();
294
295 //
296 // print the most important results of all pixels
297 //
298 calcam.Print();
299
300 //
301 // just one example how to get the plots of individual pixels
302 //
303 // calcam[18].DrawClone();
304
305 MHCamera disp1 (geomcam, "MCalibrationPix;Charge", "Fitted Mean Charges");
306 MHCamera disp3 (geomcam, "MCalibrationPix;SigmaCharge", "Sigma of Fitted Charges");
307 MHCamera disp5 (geomcam, "MCalibrationPix;ChargeProb", "Probability of Fit");
308 MHCamera disp6 (geomcam, "MCalibrationPix;Time", "Rel. Arrival Times");
309 MHCamera disp7 (geomcam, "MCalibrationPix;SigmaTime", "Sigma of Rel. Arrival Times");
310 MHCamera disp8 (geomcam, "MCalibrationPix;TimeProb", "Probability of Time Fit");
311 MHCamera disp9 (geomcam, "MCalibrationPix;Ped", "Pedestals");
312 MHCamera disp10 (geomcam, "MCalibrationPix;PedRms", "Pedestal RMS");
313 MHCamera disp11 (geomcam, "MCalibrationPix;RSigma", "Reduced Sigmas");
314 MHCamera disp12 (geomcam, "MCalibrationPix;PheFFactorMethod", "Nr. of Phe's (F-Factor Method)");
315 MHCamera disp13 (geomcam, "MCalibrationPix;MeanConversionFFactorMethod",
316 "Conversion Factor (F-Factor Method)");
317 MHCamera disp14 (geomcam, "MCalibrationPix;MeanPhotInsidePlexiglass",
318 "Nr. of Photons (Blind Pixel Method)");
319 MHCamera disp15 (geomcam, "MCalibrationPix;MeanConversionBlindPixelMethod",
320 "Conversion Factor (Blind Pixel Method)");
321 MHCamera disp16 (geomcam, "MCalibrationPix;RSigma/Charge", "Reduced Sigma per Charge");
322 MHCamera disp17 (geomcam, "MCalibrationPix;AbsTimeMean", "Abs. Arrival Times");
323 MHCamera disp18 (geomcam, "MCalibrationPix;AbsTimeRms", "RMS of Arrival Times");
324
325 disp1.SetCamContent(calcam, 0);
326 disp1.SetCamError(calcam,1);
327
328 disp3.SetCamContent(calcam, 2);
329 disp3.SetCamError(calcam,3);
330
331 disp5.SetCamContent(calcam, 4);
332
333 disp6.SetCamContent(calcam, 5);
334 disp6.SetCamError(calcam, 6);
335 disp7.SetCamContent(calcam, 6);
336 disp8.SetCamContent(calcam, 7);
337
338 disp9.SetCamContent(calcam, 8);
339 disp9.SetCamError(calcam, 9);
340
341 disp10.SetCamContent(calcam, 9);
342 disp11.SetCamContent(calcam, 10);
343
344 disp12.SetCamContent(calcam, 11);
345 disp12.SetCamError(calcam, 12);
346
347 disp13.SetCamContent(calcam, 13);
348 disp13.SetCamError(calcam, 14);
349
350 disp14.SetCamContent(calcam, 15);
351 disp15.SetCamContent(calcam, 16);
352 disp16.SetCamContent(calcam, 17);
353
354 disp17.SetCamContent(calcam, 18);
355 disp17.SetCamError(calcam, 19);
356 disp18.SetCamContent(calcam, 20);
357 disp18.SetCamError(calcam, 21);
358
359 disp1.SetYTitle("Charge [FADC counts]");
360 disp3.SetYTitle("\\sigma_{Charge} [FADC counts]");
361 disp5.SetYTitle("P_{Charge} [1]");
362 disp6.SetYTitle("Rel. Arr. Time [ns]");
363 disp7.SetYTitle("\\sigma_{Time} [ns]");
364 disp8.SetYTitle("P_{Time} [1]");
365 disp9.SetYTitle("Ped [FADC Counts ]");
366 disp10.SetYTitle("RMS_{Ped} [FADC Counts ]");
367 disp11.SetYTitle("\\sqrt{\\sigma^{2}_{Charge} - RMS^{2}_{Ped}} [FADC Counts]");
368 disp12.SetYTitle("Nr. Photo-Electrons [1]");
369 disp13.SetYTitle("Conversion Factor [PhE/FADC Count]");
370 disp14.SetYTitle("Nr. Photons [1]");
371 disp15.SetYTitle("Conversion Factor [Phot/FADC Count]");
372 disp16.SetYTitle("Reduced Sigma / Charge [1]");
373 disp17.SetYTitle("Mean Abs. Time [FADC slice]");
374 disp18.SetYTitle("RMS Abs. Time [FADC slices]");
375
376 MStatusDisplay *d3 = new MStatusDisplay;
377 d3->SetUpdateTime(3000);
378 d3->Resize(850,700);
379
380 gStyle->SetOptStat(1111);
381 gStyle->SetOptFit();
382
383 // Charges
384 TCanvas &c1 = d3->AddTab("Fitted Charges");
385 c1.Divide(2,3);
386
387 CamDraw(c1,disp1,calcam,1,2,1);
388 CamDraw(c1,disp3,calcam,2,2,1);
389
390 // Fit Probability
391 TCanvas &c2 = d3->AddTab("Fit Prob.");
392 c2.Divide(1,3);
393
394 CamDraw(c2,disp5,calcam,1,1,3);
395
396 // Rel. Times
397 TCanvas &c3 = d3->AddTab("Fitted Rel. Times");
398 c3.Divide(3,3);
399
400 CamDraw(c3,disp6,calcam,1,3,1);
401 CamDraw(c3,disp7,calcam,2,3,1);
402 CamDraw(c3,disp8,calcam,3,3,0);
403
404 // Abs. Times
405 TCanvas &c35 = d3->AddTab("Abs. Times");
406 c35.Divide(2,3);
407
408 CamDraw(c35,disp17,calcam,1,2,1);
409 CamDraw(c35,disp18,calcam,2,2,1);
410
411 // Pedestals
412 TCanvas &c4 = d3->AddTab("Pedestals");
413 c4.Divide(2,3);
414
415 CamDraw(c4,disp9,calcam,1,2,0);
416 CamDraw(c4,disp10,calcam,2,2,1);
417
418 // Reduced Sigmas
419 TCanvas &c5 = d3->AddTab("Reduced Sigmas");
420 c5.Divide(2,3);
421
422 // CamDraw(c5,disp11,calcam,1,2,1);
423 CamDraw(c5,disp11,calcam,1,2,2);
424 CamDraw(c5,disp16,calcam,2,2,1);
425
426 // F-Factor Method
427 TCanvas &c6 = d3->AddTab("F-Factor Method");
428 c6.Divide(2,3);
429
430 CamDraw(c6,disp12,calcam,1,2,1);
431 CamDraw(c6,disp13,calcam,2,2,1);
432
433 // Blind Pixel Method
434 TCanvas &c7 = d3->AddTab("Blind Pixel Method");
435 c7.Divide(2, 3);
436
437 CamDraw(c7,disp14,calcam,1,2,9);
438 CamDraw(c7,disp15,calcam,2,2,1);
439
440}
441
442void CamDraw(TCanvas &c, MHCamera &cam, MCamEvent &evt, Int_t i, Int_t j, Int_t fit)
443{
444
445 c.cd(i);
446 gPad->SetBorderMode(0);
447 MHCamera *obj1=(MHCamera*)cam.DrawCopy("hist");
448 obj1->AddNotify(evt);
449
450 c.cd(i+j);
451 gPad->SetBorderMode(0);
452 obj1->Draw();
453 ((MHCamera*)obj1)->SetPrettyPalette();
454
455 c.cd(i+2*j);
456 gPad->SetBorderMode(0);
457 TH1D *obj2 = (TH1D*)obj1->Projection();
458
459 obj2->Draw();
460 obj2->SetBit(kCanDelete);
461
462
463 const Double_t min = obj2->GetBinCenter(obj2->GetXaxis()->GetFirst());
464 const Double_t max = obj2->GetBinCenter(obj2->GetXaxis()->GetLast());
465 const Double_t integ = obj2->Integral("width")/2.5066283;
466 const Double_t mean = obj2->GetMean();
467 const Double_t rms = obj2->GetRMS();
468 const Double_t width = max-min;
469
470 if (rms == 0. || width == 0. )
471 return;
472
473 switch (fit)
474 {
475 case 0:
476 TF1 *sgaus = new TF1("sgaus","gaus(0)",min,max);
477 sgaus->SetBit(kCanDelete);
478 sgaus->SetParNames("Area","#mu","#sigma");
479 sgaus->SetParameters(integ/rms,mean,rms);
480 sgaus->SetParLimits(0,0.,integ);
481 sgaus->SetParLimits(1,min,max);
482 sgaus->SetParLimits(2,0,width/1.5);
483 obj2->Fit("sgaus","QLR");
484 obj2->GetFunction("sgaus")->SetLineColor(kYellow);
485 break;
486
487 case 1:
488 TString dgausform = "([0]-[3])/[2]*exp(-0.5*(x-[1])*(x-[1])/[2]/[2])";
489 dgausform += "+[3]/[5]*exp(-0.5*(x-[4])*(x-[4])/[5]/[5])";
490 TF1 *dgaus = new TF1("dgaus",dgausform.Data(),min,max);
491 dgaus->SetBit(kCanDelete);
492 dgaus->SetParNames("A_{tot}","#mu_{1}","#sigma_{1}","A_{2}","#mu_{2}","#sigma_{2}");
493 dgaus->SetParameters(integ,(min+mean)/2.,width/4.,
494 integ/width/2.,(max+mean)/2.,width/4.);
495 // The left-sided Gauss
496 dgaus->SetParLimits(0,integ-1.5,integ+1.5);
497 dgaus->SetParLimits(1,min+(width/10.),mean);
498 dgaus->SetParLimits(2,0,width/2.);
499 // The right-sided Gauss
500 dgaus->SetParLimits(3,0,integ);
501 dgaus->SetParLimits(4,mean,max-(width/10.));
502 dgaus->SetParLimits(5,0,width/2.);
503 obj2->Fit("dgaus","QLRM");
504 obj2->GetFunction("dgaus")->SetLineColor(kYellow);
505 break;
506
507 case 2:
508 TString tgausform = "([0]-[3]-[6])/[2]*exp(-0.5*(x-[1])*(x-[1])/[2]/[2])";
509 tgausform += "+[3]/[5]*exp(-0.5*(x-[4])*(x-[4])/[5]/[5])";
510 tgausform += "+[6]/[8]*exp(-0.5*(x-[7])*(x-[7])/[8]/[8])";
511 TF1 *tgaus = new TF1("tgaus",tgausform.Data(),min,max);
512 tgaus->SetBit(kCanDelete);
513 tgaus->SetParNames("A_{tot}","#mu_{1}","#sigma_{1}",
514 "A_{2}","#mu_{2}","#sigma_{2}",
515 "A_{3}","#mu_{3}","#sigma_{3}");
516 tgaus->SetParameters(integ,(min+mean)/2,width/4.,
517 integ/width/3.,(max+mean)/2.,width/4.,
518 integ/width/3.,mean,width/2.);
519 // The left-sided Gauss
520 tgaus->SetParLimits(0,integ-1.5,integ+1.5);
521 tgaus->SetParLimits(1,min+(width/10.),mean);
522 tgaus->SetParLimits(2,width/15.,width/2.);
523 // The right-sided Gauss
524 tgaus->SetParLimits(3,0.,integ);
525 tgaus->SetParLimits(4,mean,max-(width/10.));
526 tgaus->SetParLimits(5,width/15.,width/2.);
527 // The Gauss describing the outliers
528 tgaus->SetParLimits(6,0.,integ);
529 tgaus->SetParLimits(7,min,max);
530 tgaus->SetParLimits(8,width/4.,width/1.5);
531 obj2->Fit("tgaus","QLRM");
532 obj2->GetFunction("tgaus")->SetLineColor(kYellow);
533 break;
534 case 3:
535 obj2->Fit("pol0","Q");
536 obj2->GetFunction("pol0")->SetLineColor(kYellow);
537 break;
538 case 9:
539 break;
540 default:
541 obj2->Fit("gaus","Q");
542 obj2->GetFunction("gaus")->SetLineColor(kYellow);
543 break;
544 }
545
546 gPad->Modified();
547 gPad->Update();
548
549 /*
550 TPaveStats *st = (TPaveStats*)gPad->GetPrimitive("stats");
551 st->SetName("Projection");
552 st->SetX1NDC(0.7);
553 st->SetX2NDC(1.2);
554 st->SetY1NDC(0.4);
555 st->SetY2NDC(1.2);
556
557 gPad->Modified();
558 gPad->Update();
559
560 */
561}
Note: See TracBrowser for help on using the repository browser.