Logo ROOT   6.14/05
Reference Guide
TGLBoxPainter.cxx
Go to the documentation of this file.
1 // @(#)root/gl:$Id$
2 // Author: Timur Pocheptsov 31/08/2006
3 
4 /*************************************************************************
5  * Copyright (C) 1995-2006, 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 #include <cctype>
12 
13 #include "KeySymbols.h"
14 #include "TVirtualX.h"
15 #include "Buttons.h"
16 #include "TString.h"
17 #include "TROOT.h"
18 #include "TMath.h"
19 #include "TClass.h"
20 #include "TColor.h"
21 #include "TStyle.h"
22 #include "TH3.h"
23 #include "TVirtualMutex.h"
24 
25 #include "TPolyMarker3D.h"
26 #include "TGLPlotCamera.h"
27 #include "TGLBoxPainter.h"
28 #include "TGLIncludes.h"
29 
30 /** \class TGLBoxPainter
31 \ingroup opengl
32 Paints TH3 histograms by rendering variable-sized boxes matching the
33 bin contents.
34 */
35 
37 
38 ////////////////////////////////////////////////////////////////////////////////
39 /// Normal constructor.
40 
42  : TGLPlotPainter(hist, cam, coord, kTRUE, kTRUE, kTRUE),
43  fXOZSlice("XOZ", (TH3 *)hist, coord, &fBackBox, TGLTH3Slice::kXOZ),
44  fYOZSlice("YOZ", (TH3 *)hist, coord, &fBackBox, TGLTH3Slice::kYOZ),
45  fXOYSlice("XOY", (TH3 *)hist, coord, &fBackBox, TGLTH3Slice::kXOY),
46  fType(kBox),
47  fPolymarker(0)
48 {
49 }
50 
51 
52 ////////////////////////////////////////////////////////////////////////////////
53 /// Normal constructor.
54 
56  TGLPlotCamera *cam, TGLPlotCoordinates *coord)
57  : TGLPlotPainter(hist, cam, coord, kFALSE, kFALSE, kFALSE),
58  fXOZSlice("XOZ", (TH3 *)hist, coord, &fBackBox, TGLTH3Slice::kXOZ),
59  fYOZSlice("YOZ", (TH3 *)hist, coord, &fBackBox, TGLTH3Slice::kYOZ),
60  fXOYSlice("XOY", (TH3 *)hist, coord, &fBackBox, TGLTH3Slice::kXOY),
61  fType(kBox),
62  fPolymarker(pm)
63 {
64 }
65 
66 ////////////////////////////////////////////////////////////////////////////////
67 ///Show box info (i, j, k, binContent).
68 
70 {
71  fPlotInfo = "";
72 
73  if (fSelectedPart) {
75  if (fHist->Class())
76  fPlotInfo += fHist->Class()->GetName();
77  fPlotInfo += "::";
78  fPlotInfo += fHist->GetName();
79  } else if (!fHighColor){
80  const Int_t arr2Dsize = fCoord->GetNYBins() * fCoord->GetNZBins();
81  const Int_t binI = (fSelectedPart - fSelectionBase) / arr2Dsize + fCoord->GetFirstXBin();
82  const Int_t binJ = (fSelectedPart - fSelectionBase) % arr2Dsize / fCoord->GetNZBins() + fCoord->GetFirstYBin();
83  const Int_t binK = (fSelectedPart - fSelectionBase) % arr2Dsize % fCoord->GetNZBins() + fCoord->GetFirstZBin();
84 
85  fPlotInfo.Form("(binx = %d; biny = %d; binz = %d; binc = %f)", binI, binJ, binK,
86  fHist->GetBinContent(binI, binJ, binK));
87  } else
88  fPlotInfo = "Switch to true color mode to get correct info";
89  }
90 
91  return (Char_t *)fPlotInfo.Data();
92 }
93 
94 ////////////////////////////////////////////////////////////////////////////////
95 ///Set ranges, find min and max bin content.
96 
98 {
102 
103  if (!fCoord->SetRanges(fHist, kFALSE, kTRUE))//kFALSE == drawErrors, kTRUE == zAsBins
104  return kFALSE;
105 
108 
110  fMinMaxVal.first = fMinMaxVal.second;
111  //Bad. You can up-date some bin value and get wrong picture.
112  for (Int_t ir = fCoord->GetFirstXBin(); ir <= fCoord->GetLastXBin(); ++ir) {
113  for (Int_t jr = fCoord->GetFirstYBin(); jr <= fCoord->GetLastYBin(); ++jr) {
114  for (Int_t kr = fCoord->GetFirstZBin(); kr <= fCoord->GetLastZBin(); ++kr) {
115  fMinMaxVal.second = TMath::Max(fMinMaxVal.second, fHist->GetBinContent(ir, jr, kr));
116  fMinMaxVal.first = TMath::Min(fMinMaxVal.first, fHist->GetBinContent(ir, jr, kr));
117  }
118  }
119  }
120 
124 
125  if (fPolymarker) {
126  const Double_t xScale = fCoord->GetXScale();
127  const Double_t yScale = fCoord->GetYScale();
128  const Double_t zScale = fCoord->GetZScale();
129 
130  fPMPoints.assign(fPolymarker->GetP(), fPolymarker->GetP() + fPolymarker->GetN() * 3);
131  for (unsigned i = 0; i < fPMPoints.size(); i += 3) {
132  fPMPoints[i] *= xScale;
133  fPMPoints[i + 1] *= yScale;
134  fPMPoints[i + 2] *= zScale;
135  }
136  }
137 
138  if (fCoord->Modified()) {
144  }
145 
146  return kTRUE;
147 }
148 
149 ////////////////////////////////////////////////////////////////////////////////
150 /// User clicks right mouse button (in a pad).
151 
153 {
154  fMousePosition.fX = px;
156  fCamera->StartPan(px, py);
157  fBoxCut.StartMovement(px, fCamera->GetHeight() - py);
158 }
159 
160 ////////////////////////////////////////////////////////////////////////////////
161 /// User's moving mouse cursor, with middle mouse button pressed (for pad).
162 /// Calculate 3d shift related to 2d mouse movement.
163 
165 {
166  if (fSelectedPart >= fSelectionBase) {//Pan camera.
169 
170  fCamera->SetCamera();
172  fCamera->Pan(px, py);
173 
176  } else if (fSelectedPart > 0) {
177  //Convert py into bottom-top orientation.
178  //Possibly, move box here
179  py = fCamera->GetHeight() - py;
182 
183  fCamera->SetCamera();
185 
186 
187  if (!fHighColor) {
189  fBoxCut.MoveBox(px, py, fSelectedPart);
190  else
191  MoveSection(px, py);
192  } else {
193  MoveSection(px, py);
194  }
195 
198  }
199 
200  fMousePosition.fX = px, fMousePosition.fY = py;
202 }
203 
204 ////////////////////////////////////////////////////////////////////////////////
205 /// Box1 == spheres.
206 
208 {
209  using namespace std;//isdigit must be in std. But ...
210 
211  const Ssiz_t boxPos = option.Index("box");//"box" _already_ _exists_ in a string.
212  if (boxPos + 3 < option.Length() && isdigit(option[boxPos + 3]))
213  option[boxPos + 3] - '0' == 1 ? fType = kBox1 : fType = kBox;
214  else
215  fType = kBox;
216  option.Index("z") == kNPOS ? fDrawPalette = kFALSE : fDrawPalette = kTRUE;
217 }
218 
219 ////////////////////////////////////////////////////////////////////////////////
220 /// Remove sections.
221 
223 {
224  if (event == kButton1Double && (HasSections() || fBoxCut.IsActive())) {
228  if (fBoxCut.IsActive())
229  fBoxCut.TurnOnOff();
230  if (!gVirtualX->IsCmdThread())
231  gROOT->ProcessLineFast(Form("((TGLPlotPainter *)0x%lx)->Paint()", (ULong_t)this));
232  else
233  Paint();
234  } else if (event == kKeyPress && (py == kKey_c || py == kKey_C)) {
235  if (fHighColor)
236  Info("ProcessEvent", "Switch to true color mode to use box cut");
237  else {
238  fBoxCut.TurnOnOff();
240  }
241  }
242 }
243 
244 ////////////////////////////////////////////////////////////////////////////////
245 /// Initialize some gl state variables.
246 
248 {
249  glEnable(GL_DEPTH_TEST);
250  glEnable(GL_LIGHTING);
251  glEnable(GL_LIGHT0);
252  //For box option back polygons are culled (but not for dynamic profiles).
253  glEnable(GL_CULL_FACE);
254  glCullFace(GL_BACK);
255 
256  glLightModeli(GL_LIGHT_MODEL_TWO_SIDE, GL_TRUE);
257 }
258 
259 ////////////////////////////////////////////////////////////////////////////////
260 ///Return back some gl state variables.
261 
263 {
264  glDisable(GL_DEPTH_TEST);
265  glDisable(GL_LIGHTING);
266  glDisable(GL_LIGHT0);
267  glDisable(GL_CULL_FACE);
268  glLightModeli(GL_LIGHT_MODEL_TWO_SIDE, GL_FALSE);
269 }
270 
271 namespace {
272 
273  /////////////////////////////////////////////////////////////////////////////
274  ///
275 
276  void DrawMinusSigns(Double_t xMin, Double_t xMax, Double_t yMin, Double_t yMax,
277  Double_t zMin, Double_t zMax, Int_t fp, Bool_t onSphere, Bool_t transp)
278  {
279  const TGLDisableGuard depthTest(GL_DEPTH_TEST);
280  const TGLDisableGuard cullFace(GL_CULL_FACE);
281 
282  const Double_t ratio = onSphere ? 0.4 : 0.15;
283  const Double_t leftX = xMin + ratio * (xMax - xMin), rightX = xMax - ratio * (xMax - xMin);
284  const Double_t leftY = yMin + ratio * (yMax - yMin), rightY = yMax - ratio * (yMax - yMin);
285  const Double_t lowZ = zMin / 2. + zMax / 2. - 0.1 * (zMax - zMin);
286  const Double_t upZ = zMin / 2. + zMax / 2. + 0.1 * (zMax - zMin);
287 
288 
289  const Double_t minusVerts[][3] = {{xMin, leftY, lowZ}, {xMin, leftY, upZ}, {xMin, rightY, upZ}, {xMin, rightY, lowZ},
290  {leftX, yMin, lowZ}, {rightX, yMin, lowZ}, {rightX, yMin, upZ}, {leftX, yMin, upZ},
291  {xMax, leftY, lowZ}, {xMax, rightY, lowZ}, {xMax, rightY, upZ}, {xMax, leftY, upZ},
292  {rightX, yMax, lowZ}, {leftX, yMax, lowZ}, {leftX, yMax, upZ}, {rightX, yMax, upZ}};
293  const Int_t minusQuads[][4] = {{0, 1, 2, 3}, {4, 5, 6, 7}, {8, 9, 10, 11}, {12, 13, 14, 15}};
294 
295 
296  TGLDisableGuard light(GL_LIGHTING);
297  glColor3d(1., 0., 0.);
298 
299  const Int_t frontPlanes[][2] = {{0, 1}, {1, 2}, {2, 3}, {3, 0}};//Code duplication again :(
300  const Int_t *verts = minusQuads[frontPlanes[fp][0]];
301 
302  glBegin(GL_POLYGON);
303  glVertex3dv(minusVerts[verts[0]]);
304  glVertex3dv(minusVerts[verts[1]]);
305  glVertex3dv(minusVerts[verts[2]]);
306  glVertex3dv(minusVerts[verts[3]]);
307  glEnd();
308 
309  verts = minusQuads[frontPlanes[fp][1]];
310 
311  glBegin(GL_POLYGON);
312  glVertex3dv(minusVerts[verts[0]]);
313  glVertex3dv(minusVerts[verts[1]]);
314  glVertex3dv(minusVerts[verts[2]]);
315  glVertex3dv(minusVerts[verts[3]]);
316  glEnd();
317 
318  const Float_t nullEmission[] = {0.f, 0.f, 0.f, 1.f};
319  glMaterialfv(GL_FRONT_AND_BACK, GL_EMISSION, nullEmission);
320  glMaterialfv(GL_FRONT_AND_BACK, GL_SPECULAR, nullEmission);
321 
322  glColor4d(0., 0., 0., 0.25);
323  glPolygonMode(GL_FRONT, GL_LINE);
324 
325  if (!transp) {
326  glEnable(GL_BLEND);
327  glBlendFunc(GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA);
328  }
329 
330  glEnable(GL_LINE_SMOOTH);
331  glHint(GL_LINE_SMOOTH_HINT, GL_NICEST);
332 
333  verts = minusQuads[frontPlanes[fp][0]];
334 
335  glBegin(GL_POLYGON);
336  glVertex3dv(minusVerts[verts[0]]);
337  glVertex3dv(minusVerts[verts[1]]);
338  glVertex3dv(minusVerts[verts[2]]);
339  glVertex3dv(minusVerts[verts[3]]);
340  glEnd();
341 
342  verts = minusQuads[frontPlanes[fp][1]];
343 
344  glBegin(GL_POLYGON);
345  glVertex3dv(minusVerts[verts[0]]);
346  glVertex3dv(minusVerts[verts[1]]);
347  glVertex3dv(minusVerts[verts[2]]);
348  glVertex3dv(minusVerts[verts[3]]);
349  glEnd();
350 
351  glPolygonMode(GL_FRONT, GL_FILL);
352 
353  if (!transp)
354  glDisable(GL_BLEND);
355  }
356 
357 }
358 
359 ////////////////////////////////////////////////////////////////////////////////
360 
362 {
363  if (fPolymarker)
364  return DrawCloud();
365 
366  // Draw set of boxes (spheres)
367 
368  //Shift plot to point of origin.
369  const Rgl::PlotTranslation trGuard(this);
370 
372  glDisable(GL_CULL_FACE);
373  DrawSections();
374  glEnable(GL_CULL_FACE);
375 
376  if (!fSelectionPass) {
377  glEnable(GL_POLYGON_OFFSET_FILL);//[0
378  glPolygonOffset(1.f, 1.f);
379  SetPlotColor();
380  if (HasSections()) {
381  //Boxes are semi-transparent if we have any sections.
382  glEnable(GL_BLEND);//[1
383  glBlendFunc(GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA);
384  }
385  }
386 
387  //Using front point, find the correct order to draw boxes from
388  //back to front/from bottom to top (it's important only for semi-transparent boxes).
389  const Int_t frontPoint = fBackBox.GetFrontPoint();
390  Int_t irInit = fCoord->GetFirstXBin(), iInit = 0;
391  const Int_t nX = fCoord->GetNXBins();
392  Int_t jrInit = fCoord->GetFirstYBin(), jInit = 0;
393  const Int_t nY = fCoord->GetNYBins();
394  Int_t krInit = fCoord->GetFirstZBin(), kInit = 0;
395  const Int_t nZ = fCoord->GetNZBins();
396 
397  const Int_t addI = frontPoint == 2 || frontPoint == 1 ? 1 : (iInit = nX - 1, irInit = fCoord->GetLastXBin(), -1);
398  const Int_t addJ = frontPoint == 2 || frontPoint == 3 ? 1 : (jInit = nY - 1, jrInit = fCoord->GetLastYBin(), -1);
399  const Int_t addK = fBackBox.Get2DBox()[frontPoint + 4].Y() < fBackBox.Get2DBox()[frontPoint].Y() ? 1
400  : (kInit = nZ - 1, krInit = fCoord->GetLastZBin(),-1);
401  const Double_t xScale = fCoord->GetXScale();
402  const Double_t yScale = fCoord->GetYScale();
403  const Double_t zScale = fCoord->GetZScale();
404  const TAxis *xA = fXAxis;
405  const TAxis *yA = fYAxis;
406  const TAxis *zA = fZAxis;
407 
408  if (fSelectionPass && fHighColor)
410 
411  Double_t maxContent = TMath::Max(TMath::Abs(fMinMaxVal.first), TMath::Abs(fMinMaxVal.second));
412  if(!maxContent)//bad, find better way to check zero.
413  maxContent = 1.;
414 
415  for(Int_t ir = irInit, i = iInit; addI > 0 ? i < nX : i >= 0; ir += addI, i += addI) {
416  for(Int_t jr = jrInit, j = jInit; addJ > 0 ? j < nY : j >= 0; jr += addJ, j += addJ) {
417  for(Int_t kr = krInit, k = kInit; addK > 0 ? k < nZ : k >= 0; kr += addK, k += addK) {
418  const Double_t binContent = fHist->GetBinContent(ir, jr, kr);
419  const Double_t w = TMath::Power(TMath::Abs(binContent) / maxContent,1./3.);
420  if (!w)
421  continue;
422 
423  const Double_t xMin = xScale * (xA->GetBinLowEdge(ir) / 2 + xA->GetBinUpEdge(ir) / 2 - w * xA->GetBinWidth(ir) / 2);
424  const Double_t xMax = xScale * (xA->GetBinLowEdge(ir) / 2 + xA->GetBinUpEdge(ir) / 2 + w * xA->GetBinWidth(ir) / 2);
425  const Double_t yMin = yScale * (yA->GetBinLowEdge(jr) / 2 + yA->GetBinUpEdge(jr) / 2 - w * yA->GetBinWidth(jr) / 2);
426  const Double_t yMax = yScale * (yA->GetBinLowEdge(jr) / 2 + yA->GetBinUpEdge(jr) / 2 + w * yA->GetBinWidth(jr) / 2);
427  const Double_t zMin = zScale * (zA->GetBinLowEdge(kr) / 2 + zA->GetBinUpEdge(kr) / 2 - w * zA->GetBinWidth(kr) / 2);
428  const Double_t zMax = zScale * (zA->GetBinLowEdge(kr) / 2 + zA->GetBinUpEdge(kr) / 2 + w * zA->GetBinWidth(kr) / 2);
429 
430  if (fBoxCut.IsActive() && fBoxCut.IsInCut(xMin, xMax, yMin, yMax, zMin, zMax))
431  continue;
432 
433  const Int_t binID = fSelectionBase + i * fCoord->GetNZBins() * fCoord->GetNYBins() + j * fCoord->GetNZBins() + k;
434 
435  if (fSelectionPass && !fHighColor)
437  else if(!fHighColor && fSelectedPart == binID)
438  glMaterialfv(GL_FRONT, GL_EMISSION, Rgl::gOrangeEmission);
439 
440  if (fType == kBox) {
441  Rgl::DrawBoxFront(xMin, xMax, yMin, yMax, zMin, zMax, frontPoint);
442  } else {
443  Rgl::DrawSphere(&fQuadric, xMin, xMax, yMin, yMax, zMin, zMax);
444  }
445 
446  if (binContent < 0. && !fSelectionPass)
447  DrawMinusSigns(xMin, xMax, yMin, yMax, zMin, zMax, frontPoint, fType != kBox, HasSections());
448 
449  if (!fSelectionPass && !fHighColor && fSelectedPart == binID)
450  glMaterialfv(GL_FRONT, GL_EMISSION, Rgl::gNullEmission);
451  }
452  }
453  }
454 
455  if (fBoxCut.IsActive())
457 
458  if (!fSelectionPass && fType != kBox1) {
459  glDisable(GL_POLYGON_OFFSET_FILL);//0]
460  TGLDisableGuard lightGuard(GL_LIGHTING);//[2 - 2]
461  glColor4d(0., 0., 0., 0.25);
462  glPolygonMode(GL_FRONT, GL_LINE);//[3
463 
464  const TGLEnableGuard blendGuard(GL_BLEND);//[4-4] + 1]
465  glBlendFunc(GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA);
466  const TGLEnableGuard smoothGuard(GL_LINE_SMOOTH);//[5-5]
467  glHint(GL_LINE_SMOOTH_HINT, GL_NICEST);
468 
469  for(Int_t ir = irInit, i = iInit; addI > 0 ? i < nX : i >= 0; ir += addI, i += addI) {
470  for(Int_t jr = jrInit, j = jInit; addJ > 0 ? j < nY : j >= 0; jr += addJ, j += addJ) {
471  for(Int_t kr = krInit, k = kInit; addK > 0 ? k < nZ : k >= 0; kr += addK, k += addK) {
472  const Double_t w = TMath::Power(TMath::Abs(fHist->GetBinContent(ir, jr, kr)) / maxContent,1./3.);
473 
474  if (!w)
475  continue;
476 
477  const Double_t xMin = xScale * (xA->GetBinLowEdge(ir) / 2 + xA->GetBinUpEdge(ir) / 2 - w * xA->GetBinWidth(ir) / 2);
478  const Double_t xMax = xScale * (xA->GetBinLowEdge(ir) / 2 + xA->GetBinUpEdge(ir) / 2 + w * xA->GetBinWidth(ir) / 2);
479  const Double_t yMin = yScale * (yA->GetBinLowEdge(jr) / 2 + yA->GetBinUpEdge(jr) / 2 - w * yA->GetBinWidth(jr) / 2);
480  const Double_t yMax = yScale * (yA->GetBinLowEdge(jr) / 2 + yA->GetBinUpEdge(jr) / 2 + w * yA->GetBinWidth(jr) / 2);
481  const Double_t zMin = zScale * (zA->GetBinLowEdge(kr) / 2 + zA->GetBinUpEdge(kr) / 2 - w * zA->GetBinWidth(kr) / 2);
482  const Double_t zMax = zScale * (zA->GetBinLowEdge(kr) / 2 + zA->GetBinUpEdge(kr) / 2 + w * zA->GetBinWidth(kr) / 2);
483 
484  if (fBoxCut.IsActive() && fBoxCut.IsInCut(xMin, xMax, yMin, yMax, zMin, zMax))
485  continue;
486 
487  Rgl::DrawBoxFront(xMin, xMax, yMin, yMax, zMin, zMax, frontPoint);
488  }
489  }
490  }
491 
492  glPolygonMode(GL_FRONT, GL_FILL);//3]
493  }
494 
496  DrawPalette();
497 }
498 
499 ////////////////////////////////////////////////////////////////////////////////
500 ///Draw a frame and a polymarker inside.
501 
503 {
504  //Shift plot to the point of origin.
505  const Rgl::PlotTranslation trGuard(this);
506 
507  //Frame.
509 
511  glColor3fv(fPhysicalShapeColor);
512 
513  glDisable(GL_LIGHTING);
514 
515  const TGLVertex3 *bb = fBackBox.Get3DBox();
516  const Double_t dX = (bb[1].X() - bb[0].X()) / 40.;
517  const Double_t dY = (bb[3].Y() - bb[0].Y()) / 40.;
518  const Double_t dZ = (bb[4].Z() - bb[0].Z()) / 40.;
519  //Now, draw the cloud of points (polymarker) inside the frame.
521 
522  glEnable(GL_LIGHTING);
523 }
524 
525 ////////////////////////////////////////////////////////////////////////////////
526 /// Set boxes color.
527 
529 {
530  Float_t diffColor[] = {0.8f, 0.8f, 0.8f, 0.05f};
531 
532  if (fPhysicalShapeColor) {
533  diffColor[0] = fPhysicalShapeColor[0];
534  diffColor[1] = fPhysicalShapeColor[1];
535  diffColor[2] = fPhysicalShapeColor[2];
536  } else {
537  if (fHist->GetFillColor() != kWhite)
538  if (const TColor *c = gROOT->GetColor(fHist->GetFillColor()))
539  c->GetRGB(diffColor[0], diffColor[1], diffColor[2]);
540  }
541 
542  glMaterialfv(GL_FRONT_AND_BACK, GL_DIFFUSE, diffColor);
543  const Float_t specColor[] = {1.f, 1.f, 1.f, 1.f};
544  glMaterialfv(GL_FRONT_AND_BACK, GL_SPECULAR, specColor);
545  glMaterialf(GL_FRONT_AND_BACK, GL_SHININESS, 70.f);
546 }
547 
548 ////////////////////////////////////////////////////////////////////////////////
549 /// Draw XOZ parallel section.
550 
552 {
553  if (fSelectionPass)
554  return;
556 }
557 
558 ////////////////////////////////////////////////////////////////////////////////
559 /// Draw YOZ parallel section.
560 
562 {
563  if (fSelectionPass)
564  return;
566 }
567 
568 ////////////////////////////////////////////////////////////////////////////////
569 /// Draw XOY parallel section.
570 
572 {
573  if (fSelectionPass)
574  return;
576 }
577 
578 ////////////////////////////////////////////////////////////////////////////////
579 /// Check, if any section exists.
580 
582 {
583  return fXOZSectionPos > fBackBox.Get3DBox()[0].Y() || fYOZSectionPos> fBackBox.Get3DBox()[0].X() ||
585 }
586 
587 ////////////////////////////////////////////////////////////////////////////////
588 ///Draw. Palette.
589 ///Originally, fCamera was never null.
590 ///It can be a null now because of gl-viewer.
591 
593 {
594  if (!fCamera) {
595  //Thank you, gl-viewer!
596  return;
597  }
598 
599  const TGLLevelPalette * palette = 0;
600  const TGLVertex3 *frame = fBackBox.Get3DBox();
601 
602  if (fXOZSectionPos > frame[0].Y())
603  palette = &fXOZSlice.GetPalette();
604  else if (fYOZSectionPos > frame[0].X())
605  palette = &fYOZSlice.GetPalette();
606  else if (fXOYSectionPos > frame[0].Z())
607  palette = &fXOYSlice.GetPalette();
608 
609  if (!palette || !palette->GetPaletteSize()) {
610  return;
611  }
612 
613  Rgl::DrawPalette(fCamera, *palette);
614 
615  glFinish();
616 
617  fCamera->SetCamera();
619 }
620 
621 ////////////////////////////////////////////////////////////////////////////////
622 ///Draw. Palette. Axis.
623 
625 {
626  if (HasSections()) {
627  gVirtualX->SetDrawMode(TVirtualX::kCopy);//TCanvas by default sets in kInverse
629  }
630 }
virtual const char * GetName() const
Returns name of object.
Definition: TNamed.h:47
void DrawBox(Bool_t selectionPass, Int_t selected) const
Draw cut as a semi-transparent box.
Int_t GetFrontPoint() const
The nearest point.
Definition: TGLPlotBox.cxx:294
void RestoreModelviewMatrix() const
EBoxType fType
Definition: TGLBoxPainter.h:43
Camera for TGLPlotPainter and sub-classes.
Definition: TGLPlotCamera.h:21
Int_t GetFirstXBin() const
Definition: Buttons.h:29
Int_t GetNYBins() const
Number of Y bins.
void MoveBox(Int_t px, Int_t py, Int_t axisID)
Move box cut along selected direction.
const Float_t gNullEmission[]
Definition: TGLUtil.cxx:2814
float Float_t
Definition: RtypesCore.h:53
Bool_t GetZLog() const
Get Z log.
TGLBoxPainter(const TGLBoxPainter &)
void StartPan(Int_t px, Int_t py)
User clicks right mouse button (in a pad).
virtual Double_t GetBinLowEdge(Int_t bin) const
Return low edge of bin.
Definition: TAxis.cxx:504
const Ssiz_t kNPOS
Definition: RtypesCore.h:111
Double_t fYOZSectionPos
EGLCoordType GetCoordType() const
Get coordinates type.
void DrawSections() const
Draw sections (if any).
void StartPan(Int_t px, Int_t py)
User clicks somewhere (px, py).
virtual Double_t GetBinContent(Int_t bin) const
Return content of bin number bin.
Definition: TH1.cxx:4770
Int_t GetPaletteSize() const
Get. Palette. Size.
Definition: TGLUtil.cxx:4221
const TGLVertex3 * Get3DBox() const
Get 3D box.
Definition: TGLPlotBox.cxx:303
SCoord_t fX
Definition: TPoint.h:35
#define gROOT
Definition: TROOT.h:410
Ssiz_t Index(const char *pat, Ssiz_t i=0, ECaseCompare cmp=kExact) const
Definition: TString.h:634
SCoord_t fY
Definition: TPoint.h:36
Basic string class.
Definition: TString.h:131
Int_t GetNXBins() const
Number of X bins.
void DrawCloud() const
Draw a frame and a polymarker inside.
#define f(i)
Definition: RSha256.hxx:104
Short_t Min(Short_t a, Short_t b)
Definition: TMathBase.h:168
int Int_t
Definition: RtypesCore.h:41
bool Bool_t
Definition: RtypesCore.h:59
void TurnOnOff()
Turn the box cut on/off.
void DrawSectionXOZ() const
Draw XOZ parallel section.
TGLTH3Slice fXOZSlice
Definition: TGLBoxPainter.h:34
void DrawPalette(const TGLPlotCamera *camera, const TGLLevelPalette &palette)
Draw. Palette.
void DrawPaletteAxis(const TGLPlotCamera *camera, const Range_t &minMax, Bool_t logZ)
static void RenderPolyMarkers(const TAttMarker &marker, Char_t transp, Float_t *p, Int_t n, Int_t pick_radius=0, Bool_t selection=kFALSE, Bool_t sec_selection=kFALSE)
Render polymarkers at points specified by p-array.
Definition: TGLUtil.cxx:1941
STL namespace.
void SetMinMax(const Rgl::Range_t &newRange)
std::vector< Double_t > fPMPoints
Definition: TGLBoxPainter.h:51
void DrawPalette() const
Draw.
A slice of a TH3.
Short_t Abs(Short_t d)
Definition: TMathBase.h:108
void DrawPlot() const
Paints TH3 histograms by rendering variable-sized boxes matching the bin contents.
Definition: TGLBoxPainter.h:32
LongDouble_t Power(LongDouble_t x, LongDouble_t y)
Definition: TMath.h:734
virtual Double_t GetBinUpEdge(Int_t bin) const
Return up edge of bin.
Definition: TAxis.cxx:514
Int_t GetNZBins() const
Number of Z bins.
void Apply(Double_t phi, Double_t theta) const
Applies rotations and translations before drawing.
void ProcessEvent(Int_t event, Int_t px, Int_t py)
Remove sections.
3 component (x/y/z) vertex class.
Definition: TGLUtil.h:82
Double_t fXOZSectionPos
Bool_t IsActive() const
Bool_t InitGeometry()
Set ranges, find min and max bin content.
void Info(const char *location, const char *msgfmt,...)
void SetCamera() const
Viewport and projection.
Int_t GetFirstYBin() const
void SetYLog(Bool_t yLog)
If log changed, sections must be reset, set fModified.
TGLPlotBox fBackBox
void AddOption(const TString &stringOption)
Box1 == spheres.
void SetViewVolume(const TGLVertex3 *box)
&#39;box&#39; is the TGLPlotPainter&#39;s back box&#39;s coordinates.
The 3-D histogram classes derived from the 1-D histogram classes.
Definition: TH3.h:31
void DrawSectionXOY() const
Draw XOY parallel section.
const TPolyMarker3D * fPolymarker
Definition: TGLBoxPainter.h:50
Definition: Rtypes.h:58
Double_t fPadTheta
void DrawPaletteAxis() const
Draw. Palette. Axis.
void ObjectIDToColor(Int_t objectID, Bool_t highColor)
Object id encoded as rgb triplet.
Definition: TGLUtil.cxx:2858
const Rgl::Range_t & GetZRangeScaled() const
Scaled range.
TGLPlotCamera * fCamera
void RestoreProjectionMatrix() const
void DrawBoxFront(Double_t xMin, Double_t xMax, Double_t yMin, Double_t yMax, Double_t zMin, Double_t zMax, Int_t fp)
Draws lego&#39;s bar as a 3d box.
Definition: TGLUtil.cxx:2974
Int_t GetLastYBin() const
Class to manage histogram axis.
Definition: TAxis.h:30
Helper class for plot-painters holding information about axis ranges, numbers of bins and flags if ce...
void ResetModified()
Reset modified.
void Pan(Int_t px, Int_t py)
Pan camera.
void Form(const char *fmt,...)
Formats a string using a printf style format descriptor.
Definition: TString.cxx:2264
void DrawSectionYOZ() const
Draw YOZ parallel section.
char * Form(const char *fmt,...)
Ssiz_t Length() const
Definition: TString.h:405
void DrawSlice(Double_t pos) const
Draw slice.
void SetXLog(Bool_t xLog)
If log changed, sections must be reset, set fModified.
Base class for plot-painters that provide GL rendering of various 2D and 3D histograms, functions and parametric surfaces.
char * GetPlotInfo(Int_t px, Int_t py)
Show box info (i, j, k, binContent).
std::vector< Double_t > fZLevels
#define gVirtualX
Definition: TVirtualX.h:350
TString fPlotInfo
Definition: TGLBoxPainter.h:45
const Bool_t kFALSE
Definition: RtypesCore.h:88
void DrawBox(Int_t selectedPart, Bool_t selectionPass, const std::vector< Double_t > &zLevels, Bool_t highColor) const
Draw back box for a plot.
Definition: TGLPlotBox.cxx:184
PyObject * fType
Rgl::Range_t fMinMaxVal
Definition: TGLBoxPainter.h:46
int Ssiz_t
Definition: RtypesCore.h:63
void DeInitGL() const
Return back some gl state variables.
Double_t Z() const
Definition: TGLUtil.h:122
Double_t GetYScale() const
void Pan(Int_t px, Int_t py)
User&#39;s moving mouse cursor, with middle mouse button pressed (for pad).
#define ClassImp(name)
Definition: Rtypes.h:359
TGLTH3Slice fYOZSlice
Definition: TGLBoxPainter.h:35
double Double_t
Definition: RtypesCore.h:55
void InitGL() const
Initialize some gl state variables.
TGLQuadric fQuadric
Definition: TGLBoxPainter.h:48
Bool_t Modified() const
Modified.
unsigned long ULong_t
Definition: RtypesCore.h:51
void SetPlotColor() const
Set boxes color.
virtual Color_t GetFillColor() const
Return the fill area color.
Definition: TAttFill.h:30
void SaveProjectionMatrix() const
The TH1 histogram class.
Definition: TH1.h:56
The color creation and management class.
Definition: TColor.h:19
Int_t GetFirstZBin() const
Bool_t SetRanges(const TH1 *hist, Bool_t errors=kFALSE, Bool_t zBins=kFALSE)
Set bin ranges, ranges.
const Rgl::Range_t & GetYRangeScaled() const
Scaled range.
Int_t GetLastXBin() const
void DrawSphere(TGLQuadric *quadric, Double_t xMin, Double_t xMax, Double_t yMin, Double_t yMax, Double_t zMin, Double_t zMax)
Cylinder for lego3.
Definition: TGLUtil.cxx:3242
char Char_t
Definition: RtypesCore.h:29
Double_t fXOYSectionPos
virtual Double_t GetBinWidth(Int_t bin) const
Return bin width.
Definition: TAxis.cxx:526
A 3D polymarker.
Definition: TPolyMarker3D.h:32
const Float_t * fPhysicalShapeColor
void SetZLog(Bool_t zLog)
If log changed, sections must be reset, set fModified.
virtual Int_t GetN() const
Definition: TPolyMarker3D.h:58
void StartMovement(Int_t px, Int_t py)
Start cut&#39;s movement.
void MoveSection(Int_t px, Int_t py)
Create dynamic profile using selected plane.
Short_t Max(Short_t a, Short_t b)
Definition: TMathBase.h:200
void SetPlotBox(const Rgl::Range_t &xRange, const Rgl::Range_t &yRange, const Rgl::Range_t &zRange)
Set up a frame box.
Definition: TGLPlotBox.cxx:198
virtual Float_t * GetP() const
Definition: TPolyMarker3D.h:59
const Rgl::Range_t & GetXRangeScaled() const
Scaled range.
#define c(i)
Definition: RSha256.hxx:101
Bool_t HasSections() const
Check, if any section exists.
Bool_t fUpdateSelection
Int_t GetLastZBin() const
TGLTH3Slice fXOYSlice
Definition: TGLBoxPainter.h:36
Int_t GetHeight() const
viewport[3]
void SaveModelviewMatrix() const
TGLPlotCoordinates * fCoord
Double_t Y() const
Definition: TGLUtil.h:120
virtual void Paint()
Draw lego/surf/whatever you can.
const Bool_t kTRUE
Definition: RtypesCore.h:87
const Float_t gOrangeEmission[]
Definition: TGLUtil.cxx:2811
TGLBoxCut fBoxCut
Double_t GetZScale() const
Double_t GetXScale() const
const TGLLevelPalette & GetPalette() const
const char * Data() const
Definition: TString.h:364
const TGLVertex3 * Get2DBox() const
Get 2D box.
Definition: TGLPlotBox.cxx:312
Bool_t IsInCut(Double_t xMin, Double_t xMax, Double_t yMin, Double_t yMax, Double_t zMin, Double_t zMax) const
Check, if box defined by xmin/xmax etc. is in cut.
Double_t X() const
Definition: TGLUtil.h:118