TridHistoGLFrame.cxx

Go to the documentation of this file.

 
#include "TridHistoGLFrame.h"
#include "TridPage.h"
#include "TridSketches.h"
#include "TridOpenGLGlobal.h"
#include "TridControl.h"

#include "TROOT.h"
#include "TGX11.h"
#include "TVirtualX.h"
#include "KeySymbols.h"

#include <GL/gl.h>
#include <GL/glx.h>
#include <GL/glu.h>
#include "glf.h"
#include <cassert>
#include <iostream>
#include <cmath>
using namespace std;

TridHistoGLFrame::TridHistoGLFrame( TridPage* tp, 
                                    TGWindow& parent, 
                                    TH1* hist,
                                    TridControl* tc,
                                    HistoType_t type)
  : TridFlatGLFrame( tp, parent, tc ),
    fHist(hist),
    fHistoType(type)
{
  cout << "TridHistoGLFrame constructor." << endl;
}

TridHistoGLFrame::~TridHistoGLFrame()
{
  cout << "TridHistoGLFrame constructor." << endl;
}

void TridHistoGLFrame::SetupProjection( void )
{
  // Set up the projection for the 3d window.
  fAspectRatio = 1;
  if(fWidth>0) fAspectRatio = (double)fHeight/(double)fWidth;

  glOrtho( -0.2, 1.1,
           -0.3, 1.1,
          -1.,1.);
}



void TridHistoGLFrame::TransformToPOV( void )
{
  NormalizeAndClipPOV();
  
  // Need to manipulate the ModelView matrix to move our model around.
  glMatrixMode(GL_MODELVIEW);
  
  // Reset to 0,0,0; no rotation, no scaling.
  glLoadIdentity();   
}



Bool_t TridHistoGLFrame::HandleMotion(Event_t* ev)
{
  if(fButton != 0) {
   
    int delx, dely;
    delx = ev->fX - fMouseStart_x;
    dely = ev->fY - fMouseStart_y;
    fButtonTraveled+= abs(delx) + abs(dely);
    
    if(abs(delx) + abs(dely)>0){
      if( (fButton ==kButton1) || (fButton ==kButton3) ) {
        // Slide.
        
        GLfloat z;
        GLdouble x0, y0, z0;
        GLdouble x1, y1, z1;
        GLdouble x2, y2, z2;
        GLdouble  model[16];
        GLdouble  proj[16];
        GLint     view[4];
        glGetDoublev(GL_MODELVIEW_MATRIX, model);
        glGetDoublev(GL_PROJECTION_MATRIX, proj);
        glGetIntegerv(GL_VIEWPORT, view); 
        
        //glReadPixels(fMouseStart_x, fHeight-fMouseStart_y, 1, 1, GL_DEPTH_COMPONENT, GL_FLOAT, &z);
        z = 0;
        gluUnProject(fMouseStart_x, fHeight-fMouseStart_y,  z, model, proj, view, &x0, &y0, &z0);
        gluUnProject(fMouseLast_x, fHeight-fMouseLast_y,  z, model, proj, view, &x1, &y1, &z1);
        gluUnProject(ev->fX,        fHeight-ev->fY,       z, model, proj, view, &x2, &y2, &z2);
        
        // Click at x1,y1
        // Drag to x2,y2

        if(fTridControl==0) return true;

        // Determine where.
        if( (x0>0) && ( y0>0) ) {
          // Click was in histogram area. Use this to shift histogram without scaling.
          // Insist that minimum drag is ~ 1/2 bin width
          if( fabs(x1-x2)>0.01 ) {
            double xmin, xmax;
            if(fHistoType == kColor) {
              fTridControl->GetColorRange(xmin,xmax);
              double xshift = (x1-x2)*(xmax-xmin);
              fTridControl->ChangeColorRange(xmin+xshift,xmax+xshift);
            } else {
              fTridControl->GetTransRange(xmin,xmax);
              double xshift = (x2-x1)*(xmax-xmin);
              fTridControl->ChangeTransRange(xmin-xshift,xmax-xshift);
            }
            fMouseLast_x = ev->fX;
            fMouseLast_y = ev->fY;
          }
        }

        if( (x0<0) && ( y0>0) ) {
          // Click was in y-scale area. Scale drag.
          if( fabs(y1-y2)>0.01 ) {
            double yscale = (y1/y2);
            if(yscale<0.1) yscale = 0.1;
            if(yscale>10.0) yscale = 10.0;
            
            double newmax = fHist->GetMaximum()*yscale;
            if(newmax<0.1) newmax = 0.1;
            if(newmax>1e6) newmax = 1.e6;
            fHist->SetMaximum(newmax);      

            fMouseLast_x = ev->fX;
            fMouseLast_y = ev->fY;
          }
        }
        
        if( (x0>0) && ( y0<0) && (x1>0) && (y1<0) && (x2>0) && (y2<0) ) {
          // Click was in x-scale area. Scale drag.
          if( fabs(x1-x2)>0.01 ) {
            double xmin, xmax;
            if(fHistoType == kColor) {
              fTridControl->GetColorRange(xmin,xmax);
            } else {
              fTridControl->GetTransRange(xmin,xmax);
            }
            
            // Obsolete: allows dragging scale on left OR right.
            //if(x1<0.5) {
            //  xmin = (x1*(xmax-xmin)+xmin-xmax*x2)/(1-x2);
            //} else {
            //  xmax = xmin + x1*(xmax-xmin)/x2;
            //}
            
            xmax = xmin + x1*(xmax-xmin)/x2;

            if(fHistoType == kColor) {
              fTridControl->ChangeColorRange(xmin,xmax);
            } else {
              fTridControl->ChangeTransRange(xmin,xmax);
            }

            
            fMouseLast_x = ev->fX;
            fMouseLast_y = ev->fY;
          }
        }
                
        Update();
      }
    }
  }
  return true;
}


double GetGoodTickWidth( double min, double max, double maxticks, char* format )
{
  double dumbTickWidth = (max-min)/maxticks;
  //cout << "Dumb width " << dumbTickWidth << endl;
  double thelog = log10(dumbTickWidth);
  double multiplier = pow(10,floor(thelog));
  double abcissa = pow(10,thelog)/multiplier; // Gives a number between 1 and 9.999
  double goodTickWidth;
  //if(abcissa<3.0) {
  //  goodTickWidth = 1.0*multiplier;
  //} else if(abcissa>3.0 && abcissa<7.0) {
  int sigfigs = 1;
  if(abcissa<3.0) {
    goodTickWidth = 2.5*multiplier;
    sigfigs = 2;
  }  else if(abcissa<7.0) {
    goodTickWidth = 5.0*multiplier;
  } else {
    goodTickWidth = 10.0 * multiplier;
  }

  double extreme = fabs(max);
  if(fabs(min)>fabs(max)) extreme = fabs(min);
  if(format) {
    if(extreme>10000) strcpy(format,"%1.1e");
    else if(extreme>10) strcpy(format,"%1.1f");
    else if(extreme>1) strcpy(format,"%1.1f");
    else if(extreme>0.001) strcpy(format,"%g");
    else strcpy(format,"%1.1e");
  }
  //cout << "Good width " << goodTickWidth << endl;
  return goodTickWidth;
}

void TridHistoGLFrame::DrawObjects()
{
  if(fHist==0) return;
  assert(fTridControl.get());

  // Draw everything.
  glMatrixMode(GL_MODELVIEW);
  
  if(!fHist) return;

  // Backing box.
  TridSketchPlane backplane ( TVector3(0.5, 0.5, -0.2),
                              kv_x, 10.0,
                              kv_y, 10.0
                              );
  backplane.SetColor(fTridControl->GetBackgroundColor());
  backplane.Draw();


  TridSketch::SetDefaultColor(fTridControl->GetForegroundColor());
  TridSketch::SetDefaultTrans(1);

  

  //int iminbin = fHist->GetXaxis()->GetFirst()-1;
  //int imaxbin = fHist->GetXaxis()->GetLast()+1;
  //int nbins   = imaxbin - iminbin;
  //double xmin = fHist->GetXaxis()->GetBinLowEdge(iminbin);
  //double xmax = fHist->GetXaxis()->GetBinLowEdge(imaxbin);
  //double ymin = fHist->GetYaxis()->GetBinLowEdge(0);
  //double ymax = fHist->GetMaximum();

  int nbins = fHist->GetNbinsX() + 2;
  double xmin,xmax;
  if(fHistoType == kColor) {
    fTridControl->GetColorRange(xmin,xmax);
  } else {
    fTridControl->GetTransRange(xmin,xmax);
  }

  //double ymin = fHist->GetYaxis()->GetBinLowEdge(0);
  double ymin = 0;
  double ymax = fHist->GetMaximum();
  if(ymax<=0) ymax = 1.0;

  double textScaleX = 1.0;
  double textScaleY = 1.0;
  if(fAspectRatio>0.4) {
    textScaleY = 1.0/fAspectRatio;
  } else {
    textScaleY = 1.0/0.4;
    textScaleX = fAspectRatio/0.4;
  }

  // Draw axes.
  
  // Draw lines.
  TridSketchLine xline(TVector3(0,0,0), TVector3(1,0,0));
  TridSketchLine yline(TVector3(0,0,0), TVector3(0,1,0));
  
  xline.Draw();
  yline.Draw();

  
  // Draw ticks and labels

  // XAxis: Max 20 ticks, each has a label.
  char format[10];

  double xTickWidth = GetGoodTickWidth(xmin,xmax,6.0,format);
  double xtick = floor(xmin/xTickWidth)*xTickWidth;
  while(xtick<xmin) xtick+= xTickWidth;

  while(xtick <= xmax) {
    double screenx = (xtick-xmin)/(xmax-xmin);
    TridSketchLine tick(TVector3(screenx,0,0),TVector3(screenx,-0.05,0));
    tick.Draw();
    char buff[10];
    sprintf(buff,format,xtick);
    TridSketchText ticklabel(TVector3(screenx,-0.05,0),kv_y,kv_z,buff);
    ticklabel.SetScale(0.02*textScaleX,0.02*textScaleY);
    ticklabel.SetJustify(TridSketchText::kJustifyTop);
    ticklabel.Draw();
    xtick+= xTickWidth;    
  }

  double yTickWidth = GetGoodTickWidth(ymin,ymax,6.0,format);
  double ytick = (yTickWidth - fmod(ymin,yTickWidth)) + ymin;
    
  while(ytick <= ymax) {
    double screeny = (ytick-ymin)/(ymax-ymin);
    TridSketchLine tick(TVector3(-0.02,screeny,-0.01),TVector3(1,screeny,-0.01));
    tick.Draw();
    char buff[10];
    sprintf(buff,format,ytick);
    TridSketchText ticklabel(TVector3(-0.02,screeny,0),kv_y,kv_z,buff);
    ticklabel.SetScale(0.02*textScaleX,0.02*textScaleY);
    ticklabel.SetJustify(TridSketchText::kJustifyRight);
    ticklabel.Draw();
    ytick+= yTickWidth;    
  }
  
 
  // Draw labels.

  // Draw titles
  const char* xtitle_string;
  if(fHistoType==kColor) {
    xtitle_string =  fTridControl->GetColorModeName();
  } else {
    xtitle_string =  fTridControl->GetTransModeName();
  };
  TridSketchText xtitle(TVector3(1.0,-0.15,0),kv_y,kv_z,
                          xtitle_string);
  xtitle.SetJustify(TridSketchText::kJustifyRight | TridSketchText::kJustifyTop );
  xtitle.SetScale(0.03*textScaleX,0.03*textScaleY);
  xtitle.Draw();

  TridSketchText ytitle(TVector3(-0.15,1.0,0),kv_x,kv_z,"Counts");
  ytitle.SetJustify(TridSketchText::kJustifyRight | TridSketchText::kJustifyBottom );
  ytitle.SetScale(0.03*textScaleY,0.03*textScaleX);
  ytitle.Draw();


  TridSketchText histtitle(TVector3(-0.2,-0.3,0),kv_y,kv_z,fHist->GetName());
  histtitle.SetColor(1,0.2,0);
  histtitle.SetJustify(TridSketchText::kJustifyLeft | TridSketchText::kJustifyBottom );
  histtitle.SetScale(0.02*textScaleX,0.02*textScaleY);
  histtitle.Draw();
 
  // Draw bins.

  double pitch = 1.0/(double)(nbins-1);

  glEnable(GL_BLEND);

  for(int i = -1; i<nbins; i++) {
    double lowedge = pitch*(double)(i);
    double highedge = lowedge + pitch *0.8;
 
    double top =  (fHist->GetBinContent(i) - ymin)/(ymax-ymin);
    if(top>1.05) top = 1.05;

    TVector3 x1(lowedge,  0, 0);
    TVector3 x2(highedge, 0, 0);
    TVector3 x3(highedge, top, 0);
    TVector3 x4(lowedge,  top, 0);
    
    TridSketchPlane graphbin(x1,x2,x3,x4);
    if(fHistoType==kColor) {
      graphbin.SetColor(fTridControl->GetColor((lowedge+highedge)*0.5));
    } else {
      graphbin.SetTransparency(lowedge);
    }
    graphbin.Draw();
  }

  glDisable(GL_BLEND);
}

void TridHistoGLFrame::DrawHUD(const char*)
{
}

---