ROOT   Reference Guide
Loading...
Searching...
No Matches
TGraphPolargram.cxx
Go to the documentation of this file.
1// @(#)root/graf:$Id$
2// Author: Sebastian Boser, Mathieu Demaret 02/02/06
3
4/*************************************************************************
5 * Copyright (C) 1995-2000, Rene Brun and Fons Rademakers. *
6 * All rights reserved. *
7 * *
8 * For the licensing terms see $ROOTSYS/LICENSE. * 9 * For the list of contributors see$ROOTSYS/README/CREDITS. *
10 *************************************************************************/
11
12/** \class TGraphPolargram
13\ingroup BasicGraphics
14
15To draw polar axis
16
17TGraphPolargram draw the polar axis of the TGraphPolar.
18
19Example:
20
21Begin_Macro(source)
22{
23 TCanvas * CPol = new TCanvas("CPol","TGraphPolar Examples",500,500);
24
25 Double_t rmin=0;
26 Double_t rmax=TMath::Pi()*2;
27 Double_t r[1000];
28 Double_t theta[1000];
29
30 TF1 * fp1 = new TF1("fplot","cos(x)",rmin,rmax);
31 for (Int_t ipt = 0; ipt < 1000; ipt++) {
32 r[ipt] = ipt*(rmax-rmin)/1000+rmin;
33 theta[ipt] = fp1->Eval(r[ipt]);
34 }
35 TGraphPolar * grP1 = new TGraphPolar(1000,r,theta);
36 grP1->SetTitle("");
37 grP1->SetLineColor(2);
38 grP1->Draw("AOL");
39}
40End_Macro
41*/
42
43#include "TGraphPolar.h"
44#include "TGraphPolargram.h"
45#include "TGaxis.h"
46#include "THLimitsFinder.h"
47#include "TVirtualPad.h"
48#include "TLatex.h"
49#include "TEllipse.h"
50#include "TMath.h"
51
53
54////////////////////////////////////////////////////////////////////////////////
55/// TGraphPolargram Constructor.
56
58 Double_t tmin, Double_t tmax):
59 TNamed(name,"Polargram")
60{
61 Init();
62 fNdivRad = 508;
63 fNdivPol = 508;
64 fPolarLabels = nullptr;
65 fRwrmax = rmax;
66 fRwrmin = rmin;
67 fRwtmin = tmin;
68 fRwtmax = tmax;
69}
70
71////////////////////////////////////////////////////////////////////////////////
72/// Short constructor used in the case of a spider plot.
73
75 TNamed(name,"Polargram")
76{
77 Init();
78 fNdivRad = 0;
79 fNdivPol = 0;
80 fPolarLabels = nullptr;
81 fRwrmax = 1;
82 fRwrmin = 0;
83 fRwtmax = 0;
84 fRwtmin = 0;
85}
86
87////////////////////////////////////////////////////////////////////////////////
88/// TGraphPolargram destructor.
89
91{
92 if (fPolarLabels) delete [] fPolarLabels;
93}
94
95////////////////////////////////////////////////////////////////////////////////
96/// Set the Polar range.
97/// \param[in] tmin the start number.
98/// \param[in] tmax the end number.
99
101{
102 if (tmin < tmax) {
103 fRwtmin = tmin;
104 fRwtmax = tmax;
105 }
106 if (gPad) gPad->Modified();
107}
108
109////////////////////////////////////////////////////////////////////////////////
110/// Everything within the circle belongs to the TGraphPolargram.
111
113{
114 if (!gPad) return 9999;
115 Int_t i;
116 Double_t x = gPad->AbsPixeltoX(px);
117 Double_t y = gPad->AbsPixeltoY(py);
118
119 // Check if close to a (major) radial line.
120 Double_t rad = TMath::Sqrt(x*x+y*y);
121 Int_t div = (Int_t)rad*(fNdivRad%100);
122 Double_t dr = TMath::Min(TMath::Abs(rad-div*1./(fNdivRad%100)),
123 TMath::Abs(rad-(div+1)*1./(fNdivRad%100)));
124 Int_t drad = gPad->XtoPixel(dr)-gPad->XtoPixel(0);
125
126 // Check if close to a (major) Polar line.
127 // This is not a proper calculation, but rather fast.
128 Int_t dt = kMaxPixel;
129 for (i=0; i<(fNdivPol%100); i++) {
130 Double_t theta = i*2*TMath::Pi()/(fNdivPol%100);
131
132 // Attention: px,py in pixel units, line given in user coordinates.
133 Int_t dthis = DistancetoLine(px,py,0.,0.,TMath::Cos(theta),
134 TMath::Sin(theta));
135
136 // Fails if we are outside box described by the line.
137 // (i.e for all hor/vert lines)
138 if (dthis==9999) {
139
140 // Outside -> Get distance to endpoint of line.
141 if (rad>1) {
142 dthis = (Int_t)TMath::Sqrt(
143 TMath::Power(px-gPad->XtoPixel(TMath::Cos(theta)),2)+
144 TMath::Power(py-gPad->YtoPixel(TMath::Sin(theta)),2));
145 } else {
146
147 // Check for horizontal line.
148 if (((TMath::Abs(theta-TMath::Pi())<0.1) &&
149 ((px-gPad->XtoPixel(0))<0)) ||
150 ((TMath::Abs(theta)<0.1) &&
151 ((px-gPad->XtoPixel(0))>0))) {
152 dthis = TMath::Abs(py-gPad->YtoPixel(0.));
153 }
154
155 //Check for vertical line.
156 if (((TMath::Abs(theta-TMath::PiOver2())<0.1) &&
157 ((py-gPad->YtoPixel(0))>0)) ||
158 ((TMath::Abs(theta-3*TMath::PiOver2())<0.1) &&
159 (py-gPad->YtoPixel(0))<0)) {
160 dthis = TMath::Abs(px-gPad->XtoPixel(0.));
161 }
162 if (dthis==9999) {
163
164 // Inside, but out of box for nonorthogonal line ->
165 // get distance to start point.
166 dthis = (Int_t)TMath::Sqrt(
167 TMath::Power(px-gPad->XtoPixel(0.),2)+
168 TMath::Power(py-gPad->YtoPixel(0.),2));
169 }
170 }
171 }
172
173 // Take distance to closes line.
174 dt = TMath::Min(dthis,dt);
175 }
176 return TMath::Min(drad, dt);
177}
178
179////////////////////////////////////////////////////////////////////////////////
180/// Draw Polargram.
181
183{
184 Paint(options);
185 AppendPad(options);
186}
187
188////////////////////////////////////////////////////////////////////////////////
189/// Indicate that there is something to click here.
190
192{
193 if (!gPad) return;
194
195 Int_t kMaxDiff = 20;
196 static Int_t d1, d2, d3, px1, py1, px3, py3;
197 static Bool_t p1, p2, p3, p4, p5, p6, p7, p8;
198 Double_t px2, py2;
199 p2 = p3 = p4 = p5 = p6 = p7 = p8 = kFALSE;
200 if (!gPad->IsEditable()) return;
201 switch (event) {
202 case kMouseMotion:
203 px1 = gPad->XtoAbsPixel(TMath::Cos(GetAngle()));
204 py1 = gPad->YtoAbsPixel(TMath::Sin(GetAngle()));
205 d1 = TMath::Abs(px1 - px) + TMath::Abs(py1-py); //simply take sum of pixels differences
206 p1 = kFALSE;
207 px2 = gPad->XtoAbsPixel(-1);
208 py2 = gPad->YtoAbsPixel(1);
209 d2 = (Int_t)(TMath::Abs(px2 - px) + TMath::Abs(py2 - py)) ;
210 px3 = gPad->XtoAbsPixel(-1);
211 py3 = gPad->YtoAbsPixel(-1);
212 d3 = TMath::Abs(px3 - px) + TMath::Abs(py3 - py) ; //simply take sum of pixels differences
213 // check if point is close to the radial axis
214 if (d1 < kMaxDiff) {
215 gPad->SetCursor(kMove);
216 p1 = kTRUE;
217 }
218 // check if point is close to the left high axis
219 if ( d2 < kMaxDiff) {
220 gPad->SetCursor(kHand);
221 p7 = kTRUE;
222 }
223 // check if point is close to the left down axis
224 if ( d3 < kMaxDiff) {
225 gPad->SetCursor(kHand);
226 p8 = kTRUE;
227 }
228 // check if point is close to a main circle
229 if (!p1 && !p7 ) {
230 p6 = kTRUE;
231 gPad->SetCursor(kHand);
232 }
233 break;
234
235 case kButton1Down:
236 // Record initial coordinates
237 //px4 = px;
238 //py4 = py;
239
240 case kButton1Motion:
241 if (p1) {
242 px2 = gPad->AbsPixeltoX(px);
243 py2 = gPad->AbsPixeltoY(py);
244 if ( px2 < 0 && py2 < 0) {p2 = kTRUE;};
245 if ( px2 < 0 && py2 > 0 ) {p3 = kTRUE;};
246 if ( px2 > 0 && py2 > 0 ) {p4 = kTRUE;};
247 if ( px2 > 0 && py2 < 0 ) {p5 = kTRUE;};
248 px2 = TMath::ACos(TMath::Abs(px2));
249 py2 = TMath::ASin(TMath::Abs(py2));
250 if (p2) {
251 fAxisAngle = TMath::Pi()+(px2+py2)/2;
252 p2 = kFALSE;
253 };
254 if (p3) {
255 fAxisAngle = TMath::Pi()-(px2+py2)/2;
256 p3 = kFALSE;
257 };
258 if (p4) {
259 fAxisAngle = (px2+py2)/2;
260 p4 = kFALSE;
261 };
262 if (p5) {
263 fAxisAngle = -(px2+py2)/2;
264 p5 = kFALSE;
265 };
266 }
267 break;
268
269 case kButton1Up:
270 Paint();
271 }
272}
273
274////////////////////////////////////////////////////////////////////////////////
275/// Find the alignement rule to apply for TText::SetTextAlign(Short_t).
276
278{
279 Double_t pi = TMath::Pi();
280
281 while(angle < 0 || angle > 2*pi){
282 if(angle < 0) angle+=2*pi;
283 if(angle > 2*pi) angle-=2*pi;
284 }
286 if(angle > 0 && angle < pi/2) return 11;
287 else if(angle > pi/2 && angle < pi) return 31;
288 else if(angle > pi && angle < 3*pi/2) return 33;
289 else if(angle > 3*pi/2 && angle < 2*pi) return 13;
290 else if(angle == 0 || angle == 2*pi) return 12;
291 else if(angle == pi/2) return 21;
292 else if(angle == pi) return 32;
293 else if(angle == 3*pi/2) return 23;
294 else return 0;
295 }
296 else{
297 if(angle >= 0 && angle <= pi/2) return 12;
298 else if((angle > pi/2 && angle <= pi) || (angle > pi && angle <= 3*pi/2)) return 32;
299 else if(angle > 3*pi/2 && angle <= 2*pi) return 12;
300 else return 0;
301 }
302}
303
304////////////////////////////////////////////////////////////////////////////////
305/// Determine the orientation of the polar labels according to their angle.
306
308{
309 Double_t pi = TMath::Pi();
310 Double_t convraddeg = 180.0/pi;
311
312 while(angle < 0 || angle > 2*pi){
313 if(angle < 0) angle+=2*pi;
314 if(angle > 2*pi) angle-=2*pi;
315 }
316
317 if(angle >= 0 && angle <= pi/2) return angle*convraddeg;
318 else if(angle > pi/2 && angle <= pi) return (angle + pi)*convraddeg;
319 else if(angle > pi && angle <= 3*pi/2) return (angle - pi)*convraddeg;
320 else if(angle > 3*pi/2 && angle <= 2*pi) return angle*convraddeg;
321 else return 0;
322}
323
324////////////////////////////////////////////////////////////////////////////////
325/// Initialize some of the fields of TGraphPolargram.
326
328{
329 fAxisAngle = 0;
330 fCutRadial = 0;
331 fDegree = kFALSE;
332 fGrad = kFALSE;
333 fLineStyle = 3;
335 fPolarLabelFont = 62;
336 fPolarOffset = 0.04;
337 fPolarTextSize = 0.04;
338 fRadialOffset = 0.025;
339 fRadian = kTRUE;
341 fRadialLabelFont = 62;
342 fRadialTextSize = 0.035;
343 fTickpolarSize = 0.02;
344}
345
346////////////////////////////////////////////////////////////////////////////////
347/// Paint TGraphPolargram.
348
350{
351 Int_t optionpoldiv, optionraddiv;
352 Bool_t optionLabels = kTRUE;
353
354 TString opt = chopt;
355 opt.ToUpper();
356
357 if(opt.Contains('P')) optionpoldiv=1; else optionpoldiv=0;
358 if(opt.Contains('R')) optionraddiv=1; else optionraddiv=0;
361 if(!opt.Contains('P') && !opt.Contains('R')) optionpoldiv=optionraddiv=1;
362 if(opt.Contains('N')) optionLabels = kFALSE;
363
364 if(optionraddiv) PaintRadialDivisions(kTRUE);
366 if(optionpoldiv) PaintPolarDivisions(optionLabels);
367}
368
369////////////////////////////////////////////////////////////////////////////////
370/// This is simplified from TEllipse::PaintEllipse.
371/// Draw this ellipse with new coordinates.
372
374 Double_t phimin, Double_t phimax, Double_t theta)
375{
376 if (!gPad) return;
377
378 Int_t i;
379 const Int_t np = 200; // Number of point to draw circle
380 static Double_t x[np+3], y[np+3];
381
382 // Set number of points approximatively proportional to the ellipse
383 // circumference.
384
385 Double_t circ = TMath::Pi()*2*r*(phimax-phimin)/36;
386 Int_t n = (Int_t)(np*circ/((gPad->GetX2()-gPad->GetX1())+
387 (gPad->GetY2()-gPad->GetY1())));
388 if (n < 8) n = 8;
389 if (n > np) n = np;
390 Double_t angle,dx,dy;
391 Double_t dphi = (phimax-phimin)*TMath::Pi()/(180*n);
392 Double_t ct = TMath::Cos(TMath::Pi()*theta/180);
393 Double_t st = TMath::Sin(TMath::Pi()*theta/180);
394 for (i=0; i<=n; i++) {
395 angle = phimin*TMath::Pi()/180 + Double_t(i)*dphi;
396 dx = r*TMath::Cos(angle);
397 dy = r*TMath::Sin(angle);
398 x[i] = x1 + dx*ct - dy*st;
399 y[i] = y1 + dx*st + dy*ct;
400 }
401 gPad->PaintPolyLine(n+1,x,y);
402}
403
404////////////////////////////////////////////////////////////////////////////////
405/// Draw Polar divisions.
406/// Check for editable pad or create default.
407
409{
410 if (!gPad) return ;
411 Int_t i, j, rnum, rden, first, last;
412
413 gPad->RangeAxis(-1,-1,1,1);
414 gPad->Range(-1.25,-1.25,1.25,1.25);
415 Int_t ndivMajor = fNdivPol%100;
416 Int_t ndivMinor = fNdivPol/100;
417
418 if (!gPad->GetLogy()) {
419 for (i=0; i<ndivMajor; i++) {
420 Double_t txtval = fRwtmin + i*(fRwtmax-fRwtmin)/ndivMajor;
421 Double_t theta = i*2*TMath::Pi()/ndivMajor;
422 Double_t costheta = TMath::Cos(theta);
423 Double_t sintheta = TMath::Sin(theta);
424 Double_t tantheta = TMath::Tan(theta);
425 Double_t costhetas = (1+fPolarOffset)*costheta;
426 Double_t sinthetas = (1+fPolarOffset)*sintheta;
427 Double_t corr = 0.01;
428
429 TLatex textangular;
430 textangular.SetTextColor(GetPolarColorLabel());
431 textangular.SetTextFont(GetPolarLabelFont());
432
433 TString form = " ";
434 TGaxis axis;
436 // Polar numbers are aligned with their axis.
437 if(!fPolarLabels && optionLabels){;
438 if (fRadian) {
439 // Radian case.
440 ReduceFraction(2*i, ndivMajor, rnum, rden); // Reduces the fraction.
441 if (rnum == 0) form.Form("%d",rnum);
442 if (rnum == 1 && rden == 1) form = "#pi";
443 if (rnum == 1 && rden != 1) form.Form("#frac{#pi}{%d}",rden);
444 if (rnum != 1 && rden == 1 && i !=0) form.Form("%d#pi",rnum);
445 if (rnum != 1 && rden != 1) form.Form("#frac{%d#pi}{%d}",rnum,rden);
446 textangular.SetTextAlign(FindAlign(theta));
447 textangular.PaintLatex(costhetas,
448 sinthetas, FindTextAngle(theta),
449 GetPolarLabelSize(), form.Data());
450 } else {
451 // Any other cases: numbers are aligned with their axis.
452 form.Form("%5.3g",txtval);
453 axis.LabelsLimits(form.Data(),first,last);
454 TString s = form;
455 if (first != 0) s.Remove(0, first);
456 textangular.SetTextAlign(FindAlign(theta));
457 textangular.PaintLatex(costhetas,
458 sinthetas, FindTextAngle(theta),
459 GetPolarLabelSize(), s);
460 }
461 } else if (fPolarLabels){
462 // print the specified polar labels
463 textangular.SetTextAlign(FindAlign(theta));
464 textangular.PaintLatex(costhetas,sinthetas,FindTextAngle(theta),
466 }
467 } else {
468 // Polar numbers are shown horizontally.
469 if(!fPolarLabels && optionLabels){
470 if (fRadian) {
471 // Radian case
472 ReduceFraction(2*i, ndivMajor, rnum, rden);
473 if (rnum == 0) form.Form("%d",rnum);
474 if (rnum == 1 && rden == 1) form = "#pi";
475 if (rnum == 1 && rden != 1) form.Form("#frac{#pi}{%d}",rden);
476 if (rnum != 1 && rden == 1 && i !=0) form.Form("%d#pi",rnum);
477 if (rnum != 1 && rden != 1) form.Form("#frac{%d#pi}{%d}",rnum,rden);
478 if(theta >= 3*TMath::Pi()/12.0 && theta < 2*TMath::Pi()/3.0) corr=0.04;
479 textangular.SetTextAlign(FindAlign(theta));
480 textangular.PaintLatex(costhetas,corr+sinthetas,0,
481 GetPolarLabelSize(),form.Data());
482 } else {
483 // Any other cases where numbers are shown horizontally.
484 form.Form("%5.3g",txtval);
485 axis.LabelsLimits(form.Data(),first,last);
486 TString s = form;
487 if (first != 0) s.Remove(0, first);
488 if(theta >= 3*TMath::Pi()/12.0 && theta < 2*TMath::Pi()/3.0) corr=0.04;
489 textangular.SetTextAlign(FindAlign(theta));
490 textangular.PaintLatex(costhetas, //j'ai efface des offset la
491 corr+sinthetas,0,GetPolarLabelSize(),s);
492 }
493 } else if (fPolarLabels) {
494 // print the specified polar labels
495 textangular.SetTextAlign(FindAlign(theta));
496 textangular.PaintText(costhetas,sinthetas,fPolarLabels[i]);
497 }
498 }
500 //Check if SetTickPolar is activated, and draw tick marks
501 Bool_t issettickpolar = gPad->GetTicky();
502
503 if (issettickpolar) {
504 if (theta != 0 && theta !=TMath::Pi()) {
505 gPad->PaintLine((sintheta-GetTickpolarSize())/tantheta,sintheta-GetTickpolarSize(),
506 (sintheta+GetTickpolarSize())/tantheta,sintheta+GetTickpolarSize());
507 }
508 if (theta == 0 || theta ==TMath::Pi()) {
509 gPad->PaintLine(1-GetTickpolarSize(),0,1+GetTickpolarSize(),0);
510 gPad->PaintLine(-1+GetTickpolarSize(),0,-1-GetTickpolarSize(),0);
511 }
512 }
515 gPad->PaintLine(0.,0.,costheta,sintheta);
516 // Add minor lines w/o text.
517 Int_t oldLineStyle = GetLineStyle();
518 TAttLine::SetLineStyle(2); //Minor lines always in this style.
519 TAttLine::Modify(); //Changes line attributes apart from style.
520 for (j=1; j<ndivMinor; j++) {
521 Double_t thetamin = theta+j*2*TMath::Pi()/(ndivMajor*ndivMinor);
522 gPad->PaintLine(0.,0.,TMath::Cos(thetamin),TMath::Sin(thetamin));
523 }
524 TAttLine::SetLineStyle(oldLineStyle);
526 }
527 } else {
528 Int_t big = (Int_t)fRwtmax;
529 Int_t test= 1;
530 while (big >= 10) {
531 big = big/10;
532 test++;
533 }
534 for (i=1; i<=test; i++) {
535 Double_t txtval = pow((double)10,(double)(i-1));
536 Double_t theta = (i-1)*2*TMath::Pi()/(double)(test);
537 Double_t costheta = TMath::Cos(theta);
538 Double_t sintheta = TMath::Sin(theta);
539 Double_t tantheta = TMath::Tan(theta);
540 Double_t costhetas = (1+fPolarOffset)*costheta;
541 Double_t sinthetas = (1+fPolarOffset)*sintheta;
542 Double_t corr = 0.01;
543
544 TLatex textangular;
545 textangular.SetTextColor(GetPolarColorLabel());
546 textangular.SetTextFont(GetPolarLabelFont());
547
548 TString form = " ";
549 TGaxis axis;
550
552 if(!fPolarLabels && optionLabels){
553 // Polar numbers are aligned with their axis.
554 form.Form("%5.3g",txtval);
555 axis.LabelsLimits(form.Data(),first,last);
556 TString s = form;
557 if (first != 0) s.Remove(0, first);
558 textangular.SetTextAlign(FindAlign(theta));
559 textangular.PaintLatex(costhetas,
560 sinthetas, FindTextAngle(theta), GetPolarLabelSize(), s);
561 }
562 else if (fPolarLabels){
563 // print the specified polar labels
564 textangular.SetTextAlign(FindAlign(theta));
565 textangular.PaintText(costhetas,sinthetas,fPolarLabels[i]);
566 }
567
568 } else {
569 if(!fPolarLabels && optionLabels){
570 // Polar numbers are shown horizontally.
571 form.Form("%5.3g",txtval);
572 axis.LabelsLimits(form.Data(),first,last);
573 TString s = form;
574 if (first != 0) s.Remove(0, first);
575 if(theta >= 3*TMath::Pi()/12.0 && theta < 2*TMath::Pi()/3.0) corr=0.04;
576 textangular.SetTextAlign(FindAlign(theta));
577 textangular.PaintLatex(costhetas,
578 corr+sinthetas,0,GetPolarLabelSize(),s);
579 } else if (fPolarLabels){
580 // print the specified polar labels
581 textangular.SetTextAlign(FindAlign(theta));
582 textangular.PaintText(costhetas,sinthetas,fPolarLabels[i]);
583 }
584 }
585
587 //Check if SetTickPolar is activated, and draw tick marks
588 Bool_t issettickpolar = gPad->GetTicky();
589 if (issettickpolar) {
590 if (theta != 0 && theta !=TMath::Pi()) {
591 gPad->PaintLine((sintheta-GetTickpolarSize())/tantheta,sintheta-GetTickpolarSize(),
592 (sintheta+GetTickpolarSize())/tantheta,sintheta+GetTickpolarSize());
593 }
594 if (theta == 0 || theta ==TMath::Pi()) {
595 gPad->PaintLine(1-GetTickpolarSize(),0,1+GetTickpolarSize(),0);
596 gPad->PaintLine(-1+GetTickpolarSize(),0,-1-GetTickpolarSize(),0);
597 }
598 }
601 gPad->PaintLine(0.,0.,costheta,sintheta);
602 // Add minor lines w/o text.
603 Int_t oldLineStyle = GetLineStyle();
604 TAttLine::SetLineStyle(2); //Minor lines always in this style.
605 TAttLine::Modify(); //Changes line attributes apart from style.
606 Double_t a=0;
607 Double_t b,c,d;
608 b = TMath::Log(10)*test;
609 d= 2*TMath::Pi()/(double)test;
610 for (j=1; j<9; j++) {
611 a=TMath::Log(j+1)-TMath::Log(j)+a;
612 c=a/b*6.28+d*(i-1);
613 gPad->PaintLine(0.,0.,TMath::Cos(c),TMath::Sin(c));
614 }
615 TAttLine::SetLineStyle(oldLineStyle);
617 }
618 }
619}
620
621////////////////////////////////////////////////////////////////////////////////
622/// Paint radial divisions.
623/// Check for editable pad or create default.
624
626{
627 if (!gPad) return ;
628
629 static char chopt[8] = "";
630 Int_t i,j;
631 Int_t ndiv = TMath::Abs(fNdivRad);
632 Int_t ndivMajor = ndiv%100;
633 Int_t ndivMinor = ndiv/100;
634 Int_t ndivmajor = 0;
635 Double_t frwrmin = 0., frwrmax = 0., binWidth = 0;
636
637 THLimitsFinder::Optimize(fRwrmin,fRwrmax,ndivMajor,frwrmin,
638 frwrmax, ndivmajor,binWidth,"");
639
640 if (!gPad->GetLogx()) {
641 gPad->RangeAxis(-1,-1,1,1);
642 gPad->Range(-1.25,-1.25,1.25,1.25);
643 Double_t umin = fRwrmin;
644 Double_t umax = fRwrmax;
645 Double_t rmajmin = (frwrmin-fRwrmin)/(fRwrmax-fRwrmin);
646 Double_t rmajmax = (frwrmax-fRwrmin)/(fRwrmax-fRwrmin);
647 Double_t dist = (rmajmax-rmajmin)/ndivmajor;
648 Int_t ndivminor = 0;
649
650 chopt[0] = 0;
651 strncat(chopt, "SDH", 4);
652 if (fNdivRad < 0) strncat(chopt, "N",2);
653 if(drawaxis){
654 // Paint axis.
655 TGaxis axis;
661 umin, umax, ndiv, chopt, 0., kFALSE);
662 }
663
664 // Paint Circles.
665 // First paint main circle.
666 PaintCircle(0.,0.,1,0.,360,0);
667 // Optimised case.
668 if (fNdivRad>0 ) {
669 Double_t frwrmini = 0., frwrmaxi = 0., binWidth2 =0;
670 THLimitsFinder::Optimize(frwrmin,frwrmin+binWidth,ndivMinor,frwrmini,
671 frwrmaxi, ndivminor,binWidth2,"");
672 Double_t dist2 = dist/(ndivminor);
673 // Paint major circles.
674 for (i=1; i<=ndivmajor+2; i++) {
677 PaintCircle(0.,0.,rmajmin,0.,360,0);
678
679 //Paint minor circles.
682 for (j=1; j<ndivminor+1; j++) {
683 if (rmajmin+j*dist2<=1) PaintCircle(0.,0.,rmajmin+j*dist2,0.,360,0);
684 }
685 rmajmin = (frwrmin-fRwrmin)/(fRwrmax-fRwrmin)+(i-1)*dist;
686 }
687 // Non-optimized case.
688 } else {
689
690 // Paint major circles.
691 for (i=1; i<=ndivMajor; i++) {
694 Double_t rmaj = i*1./ndivMajor;
695 PaintCircle(0.,0.,rmaj,0.,360,0);
696
697 // Paint minor circles.
698 for (j=1; j<ndivMinor; j++) {
701 PaintCircle(0.,0.,rmaj- j*1./(ndivMajor*ndivMinor),0.,360,0);
702 }
703 }
704 }
705 } else {
706 // Draw Log scale on radial axis if option activated.
707 Int_t big = (Int_t)fRwrmax;
708 Int_t test= 0;
709 while (big >= 10) {
710 big = big/10;
711 test++;
712 }
713 for (i=1; i<=test; i++) {
716 Double_t ecart;
717 ecart = ((double) i)/ ((double) test);
718 PaintCircle(0.,0.,ecart,0,360,0);
721 Double_t a=0;
722 Double_t b,c,d;
723 b = TMath::Log(10)*test;
724 d = 1/(double)test;
725 for (j=1; j<9; j++) {
726 a = TMath::Log(j+1)-TMath::Log(j)+a;
727 c = a/b+d*(i-1);
728 PaintCircle(0,0.,c,0.,360,0);
729 }
730 }
731 }
734}
735
736////////////////////////////////////////////////////////////////////////////////
737/// Reduce fractions.
738
740{
741 Int_t a = 0;
742 Int_t b = 0;
743 Int_t i = 0;
744 Int_t j = 0;
745 a = den;
746 b = num;
747 if (b > a) {
748 j = b;
749 } else {
750 j = a;
751 }
752 for (i=j; i > 1; i--) {
753 if ((a % i == 0) && (b % i == 0)) {
754 a = a/i;
755 b = b/i;
756 }
757 }
758 rden = a;
759 rnum = b;
760}
761
762////////////////////////////////////////////////////////////////////////////////
763/// Set axis angle.
764
766{
767 fAxisAngle = angle/180*TMath::Pi();
768}
769
770////////////////////////////////////////////////////////////////////////////////
771/// Set the number of Polar divisions: enter a number ij with 0<i<99 and 0<j<99
772/// - i sets the major Polar divisions.
773/// - j sets the minor Polar divisions.
774
776{
777 if (ndiv > 0)
778 fNdivPol = ndiv;
779 if (gPad) gPad->Modified();
780}
781
782////////////////////////////////////////////////////////////////////////////////
783/// Set the number of radial divisions: enter a number ij with 0<i<99 and 0<j<99
784/// - i sets the major radial divisions.
785/// - j sets the minor radial divisions.
786
788{
789 fNdivRad = ndiv;
790 if (gPad) gPad->Modified();
791}
792
793////////////////////////////////////////////////////////////////////////////////
794/// Set some specified polar labels, used in the case of a spider plot.
795
797{
798 if(!fPolarLabels)
800 fPolarLabels[div] = label;
801 if (gPad) gPad->Modified();
802}
803
804////////////////////////////////////////////////////////////////////////////////
805/// Set Polar labels color.
806
808{
809 fPolarLabelColor = tcolorangular;
810}
811
812////////////////////////////////////////////////////////////////////////////////
813
815{
816 // Set Polar label font.
817
818 fPolarLabelFont = tfontangular;
819}
820
821////////////////////////////////////////////////////////////////////////////////
822/// Set angular labels size.
823
825{
826 fPolarTextSize = angularsize;
827}
828
829////////////////////////////////////////////////////////////////////////////////
830/// Set the labels offset.
831
833{
834 fPolarOffset = angularOffset;
835 if (gPad) gPad->Modified();
836}
837
838////////////////////////////////////////////////////////////////////////////////
839/// Set radial labels color.
840
842{
843 fRadialLabelColor = tcolorradial;
844}
845
846////////////////////////////////////////////////////////////////////////////////
847/// Set radial label font.
848
850{
851 fRadialLabelFont = tfontradial;
852}
853
854////////////////////////////////////////////////////////////////////////////////
855/// Set radial labels size.
856
858{
859 fRadialTextSize = radialsize;
860}
861
862////////////////////////////////////////////////////////////////////////////////
863/// Set the labels offset.
864
866{
867 fRadialOffset = radialOffset;
868 if (gPad) gPad->Modified();
869}
870
871////////////////////////////////////////////////////////////////////////////////
872/// Allows to change range Polar.
873/// \param[in] tmin the start number.
874/// \param[in] tmax the end number.
875
877{
878 fDegree = kFALSE;
879 fGrad = kFALSE;
880 fRadian = kFALSE;
881
882 if (tmin < tmax) {
883 fRwtmin = tmin;
884 fRwtmax = tmax;
885 }
886 if (gPad) gPad->Modified();
887}
888
889////////////////////////////////////////////////////////////////////////////////
890/// Set the radial range.
891/// \param[in] rmin radius at center of the circle.
892/// \param[in] rmax radius at the intersection of the right X axis part and the circle.
893
895{
896 if (rmin < rmax) {
897 fRwrmin = rmin;
898 fRwrmax = rmax;
899 }
900 if (gPad) gPad->Modified();
901}
902
903////////////////////////////////////////////////////////////////////////////////
904/// Set polar ticks size.
905
907{
908 fTickpolarSize = tickpolarsize;
909}
910
911////////////////////////////////////////////////////////////////////////////////
912/// The Polar circle is labelled using degree.
913
915{
916 fDegree = kTRUE;
917 fGrad = kFALSE;
918 fRadian = kFALSE;
919 ChangeRangePolar(0,360);
920}
921
922////////////////////////////////////////////////////////////////////////////////
923/// The Polar circle is labelled using gradian.
924
926{
927 fGrad = kTRUE;
928 fRadian = kFALSE;
929 fDegree = kFALSE;
930 ChangeRangePolar(0,200);
931}
932
933////////////////////////////////////////////////////////////////////////////////
934/// The Polar circle is labelled using radian.
935
937{
938 fRadian = kTRUE;
939 fGrad = kFALSE;
940 fDegree = kFALSE;
942}
943
944////////////////////////////////////////////////////////////////////////////////
945/// Set range from 0 to 2*pi
946
948{
950}
@ kMouseMotion
Definition Buttons.h:23
@ kButton1Motion
Definition Buttons.h:20
@ kButton1Up
Definition Buttons.h:19
@ kButton1Down
Definition Buttons.h:17
@ kMove
Definition GuiTypes.h:374
@ kHand
Definition GuiTypes.h:374
const Int_t kMaxPixel
Max value for an int.
Definition GuiTypes.h:369
#define d(i)
Definition RSha256.hxx:102
#define b(i)
Definition RSha256.hxx:100
#define c(i)
Definition RSha256.hxx:101
#define a(i)
Definition RSha256.hxx:99
int Int_t
Definition RtypesCore.h:45
short Color_t
Definition RtypesCore.h:92
short Font_t
Definition RtypesCore.h:88
constexpr Bool_t kFALSE
Definition RtypesCore.h:101
double Double_t
Definition RtypesCore.h:59
constexpr Bool_t kTRUE
Definition RtypesCore.h:100
const char Option_t
Definition RtypesCore.h:66
#define ClassImp(name)
Definition Rtypes.h:377
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void char Point_t Rectangle_t WindowAttributes_t Float_t Float_t Float_t Int_t Int_t UInt_t UInt_t Rectangle_t Int_t Int_t Window_t TString Int_t GCValues_t GetPrimarySelectionOwner GetDisplay GetScreen GetColormap GetNativeEvent const char const char dpyName wid window const char font_name cursor keysym reg const char only_if_exist regb h Point_t np
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void char Point_t Rectangle_t WindowAttributes_t Float_t r
Option_t Option_t TPoint TPoint const char x1
Option_t Option_t TPoint TPoint angle
Option_t Option_t TPoint TPoint const char y1
char name[80]
Definition TGX11.cxx:110
#define gPad
virtual void SetLineStyle(Style_t lstyle)
Set the line style.
Definition TAttLine.h:42
virtual Style_t GetLineStyle() const
Return the line style.
Definition TAttLine.h:34
virtual void Modify()
Change current line attributes if necessary.
Definition TAttLine.cxx:245
Style_t fLineStyle
Line style.
Definition TAttLine.h:22
Int_t DistancetoLine(Int_t px, Int_t py, Double_t xp1, Double_t yp1, Double_t xp2, Double_t yp2)
Compute distance from point px,py to a line.
Definition TAttLine.cxx:209
virtual void SetTextAlign(Short_t align=11)
Set the text alignment.
Definition TAttText.h:42
virtual void SetTextColor(Color_t tcolor=1)
Set the text color.
Definition TAttText.h:44
virtual void SetTextFont(Font_t tfont=62)
Set the text font.
Definition TAttText.h:46
The axis painter class.
Definition TGaxis.h:24
virtual void PaintAxis(Double_t xmin, Double_t ymin, Double_t xmax, Double_t ymax, Double_t &wmin, Double_t &wmax, Int_t &ndiv, Option_t *chopt="", Double_t gridlength=0, Bool_t drawGridOnly=kFALSE)
Control function to draw an axis.
Definition TGaxis.cxx:1008
void SetLabelFont(Int_t labelfont)
Definition TGaxis.h:105
void SetLabelOffset(Float_t labeloffset)
Definition TGaxis.h:106
void SetLabelColor(Int_t labelcolor)
Definition TGaxis.h:104
void SetLabelSize(Float_t labelsize)
Definition TGaxis.h:107
void LabelsLimits(const char *label, Int_t &first, Int_t &last)
Internal method to find first and last character of a label.
Definition TGaxis.cxx:2533
To draw polar axis.
Double_t fRwrmin
Minimal radial value (real world)
void ReduceFraction(Int_t Num, Int_t Denom, Int_t &rnum, Int_t &rden)
Reduce fractions.
Double_t GetPolarLabelSize()
Color_t fRadialLabelColor
Set color of the radial labels.
void SetPolarLabelColor(Color_t tcolorangular=1)
Set Polar labels color.
void PaintPolarDivisions(Bool_t noLabels)
Draw Polar divisions.
Int_t DistancetoPrimitive(Int_t px, Int_t py) override
Everything within the circle belongs to the TGraphPolargram.
void SetRangePolar(Double_t tmin, Double_t tmax)
Allows to change range Polar.
Double_t fPolarTextSize
Set Polar text size.
Double_t fRwtmin
Minimal angular value (real world)
Double_t GetRadialOffset()
void ExecuteEvent(Int_t event, Int_t px, Int_t py) override
Indicate that there is something to click here.
void Init()
Initialize some of the fields of TGraphPolargram.
Double_t fRadialTextSize
Double_t GetTickpolarSize()
void SetPolarOffset(Double_t PolarOffset=0.04)
Set the labels offset.
void SetTwoPi()
Set range from 0 to 2*pi.
Double_t FindTextAngle(Double_t theta)
Determine the orientation of the polar labels according to their angle.
void SetRadialLabelColor(Color_t tcolorradial=1)
Set radial labels color.
void PaintCircle(Double_t x, Double_t y, Double_t r, Double_t phimin, Double_t phimax, Double_t theta)
This is simplified from TEllipse::PaintEllipse.
void SetRangeRadial(Double_t rmin, Double_t rmax)
Set the radial range.
void SetTickpolarSize(Double_t tickpolarsize=0.02)
Set polar ticks size.
void SetAxisAngle(Double_t angle=0)
Set axis angle.
void SetToDegree()
The Polar circle is labelled using degree.
void SetNdivPolar(Int_t Ndiv=508)
Set the number of Polar divisions: enter a number ij with 0<i<99 and 0<j<99.
Int_t fNdivRad
Number of radial divisions.
Font_t GetPolarLabelFont()
void Draw(Option_t *options="") override
Draw Polargram.
Double_t fTickpolarSize
Set size of Tickmarks.
void SetRadialLabelSize(Double_t radialsize=0.035)
Set radial labels size.
Font_t fPolarLabelFont
Set font of angular labels.
TString * fPolarLabels
! [fNdivPol] Specified polar labels
void PaintRadialDivisions(Bool_t drawaxis)
Paint radial divisions.
void SetToGrad()
The Polar circle is labelled using gradian.
Double_t GetAngle()
~TGraphPolargram() override
TGraphPolargram destructor.
void ChangeRangePolar(Double_t tmin, Double_t tmax)
Set the Polar range.
void SetPolarLabelSize(Double_t angularsize=0.04)
Set angular labels size.
void SetRadialLabelFont(Font_t tfontradial=62)
Set radial label font.
Double_t fAxisAngle
Set angle of the radial axis.
void SetToRadian()
The Polar circle is labelled using radian.
Double_t fRwtmax
Minimal angular value (real world)
Double_t fRwrmax
Maximal radial value (real world)
void SetPolarLabel(Int_t div, const TString &label)
Set some specified polar labels, used in the case of a spider plot.
TGraphPolargram(const char *name, Double_t rmin, Double_t rmax, Double_t tmin, Double_t tmax)
TGraphPolargram Constructor.
Font_t fRadialLabelFont
Set font of radial labels.
Int_t fNdivPol
Number of polar divisions.
Font_t GetRadialLabelFont()
Int_t fCutRadial
if fCutRadial = 0, circles are cut by radial axis if fCutRadial = 1, circles are not cut
Color_t fPolarLabelColor
Set color of the angular labels.
Double_t fPolarOffset
Offset for Polar labels.
void Paint(Option_t *options="") override
Paint TGraphPolargram.
void SetPolarLabelFont(Font_t tfontangular=62)
void SetRadialOffset(Double_t RadialOffset=0.025)
Set the labels offset.
Color_t GetRadialColorLabel()
Int_t FindAlign(Double_t angle)
Find the alignement rule to apply for TText::SetTextAlign(Short_t).
void SetNdivRadial(Int_t Ndiv=508)
Set the number of radial divisions: enter a number ij with 0<i<99 and 0<j<99.
Double_t GetRadialLabelSize()
Color_t GetPolarColorLabel()
Double_t fRadialOffset
Offset for radial labels.
static void Optimize(Double_t A1, Double_t A2, Int_t nold, Double_t &BinLow, Double_t &BinHigh, Int_t &nbins, Double_t &BWID, Option_t *option="")
Static function to compute reasonable axis limits.
To draw Mathematical Formula.
Definition TLatex.h:18
virtual void PaintLatex(Double_t x, Double_t y, Double_t angle, Double_t size, const char *text)
Main drawing function.
Definition TLatex.cxx:2057
The TNamed class is the base class for all named ROOT classes.
Definition TNamed.h:29
R__ALWAYS_INLINE Bool_t TestBit(UInt_t f) const
Definition TObject.h:199
virtual void AppendPad(Option_t *option="")
Append graphics object to current pad.
Definition TObject.cxx:184
void SetBit(UInt_t f, Bool_t set)
Set or unset the user status bits as specified in f.
Definition TObject.cxx:780
void ResetBit(UInt_t f)
Definition TObject.h:198
Basic string class.
Definition TString.h:139
const char * Data() const
Definition TString.h:380
void ToUpper()
Change string to upper case.
Definition TString.cxx:1183
TString & Remove(Ssiz_t pos)
Definition TString.h:685
void Form(const char *fmt,...)
Formats a string using a printf style format descriptor.
Definition TString.cxx:2334
Bool_t Contains(const char *pat, ECaseCompare cmp=kExact) const
Definition TString.h:636
virtual void PaintText(Double_t x, Double_t y, const char *text)
Draw this text with new coordinates.
Definition TText.cxx:752
Double_t y[n]
Definition legend1.C:17
Double_t x[n]
Definition legend1.C:17
const Int_t n
Definition legend1.C:16
Double_t ACos(Double_t)
Returns the principal value of the arc cosine of x, expressed in radians.
Definition TMath.h:632
Double_t ASin(Double_t)
Returns the principal value of the arc sine of x, expressed in radians.
Definition TMath.h:624
constexpr Double_t PiOver2()
Definition TMath.h:51
Double_t Log(Double_t x)
Returns the natural logarithm of x.
Definition TMath.h:756
Double_t Sqrt(Double_t x)
Returns the square root of x.
Definition TMath.h:662
LongDouble_t Power(LongDouble_t x, LongDouble_t y)
Returns x raised to the power y.
Definition TMath.h:721
Short_t Min(Short_t a, Short_t b)
Returns the smallest of a and b.
Definition TMathBase.h:198
Double_t Cos(Double_t)
Returns the cosine of an angle of x radians.
Definition TMath.h:594
constexpr Double_t Pi()
Definition TMath.h:37
Double_t Sin(Double_t)
Returns the sine of an angle of x radians.
Definition TMath.h:588
Double_t Tan(Double_t)
Returns the tangent of an angle of x radians.
Definition TMath.h:600
Short_t Abs(Short_t d)
Returns the absolute value of parameter Short_t d.
Definition TMathBase.h:123
Definition first.py:1
Definition test.py:1