glsurfaces.C

Go to the documentation of this file.

/// \file
/// \ingroup tutorial_gl
/// Various surfaces rendered with GL.
/// This example draws 6 surfaces using OpenGL in pad (one is remake
/// of a classic surfaces.C, another contains 4 surfaces).
///
/// The commands used are exactly the same as with a normal pad.
/// The only command to add is: `gStyle->SetCanvasPreferGL(true);`
///
/// \macro_image(nobatch)
/// \macro_code
///
/// \authors Rene Brun, Timur Pocheptsov
 
void glsurfaces()
{
   gStyle->SetPalette(0);
 
   // after this command all legos surfaces (surf/surf1/surf2/surf4/tf3
   // options) are automatically rendered with OpenGL.
   gStyle->SetCanvasPreferGL(kTRUE);
 
   TCanvas *c1 = new TCanvas("glc1", "Surfaces Drawing Options", 200, 10, 700, 900);
   c1->SetFillColor(42);
   gStyle->SetFrameFillColor(42);
   auto title = new TPaveText(0.2, 0.96, 0.8, 0.995);
   title->SetFillColor(33);
   title->AddText("Examples of Surface options");
   title->Draw();
 
   TPad *pad1 = new TPad("pad1", "Gouraud shading", 0.03, 0.50, 0.98, 0.95, 21);
   TPad *pad2 = new TPad("pad2", "Color mesh", 0.03, 0.02, 0.98, 0.48, 21);
   pad1->Draw();
   pad2->Draw();
   // We generate a 2-D function
   TF2 *f2 = new TF2("f2", "x**2 + y**2 - x**3 -8*x*y**4", -1., 1.2, -1.5, 1.5);
   // Draw this function in pad1 with Gouraud shading option
   pad1->cd();
   pad1->SetLogz();
   f2->SetFillColor(45);
   f2->Draw("glsurf4");
 
   TF2 *f2clone = new TF2("f2clone", "x**2 + y**2 - x**3 -8*x*y**4", -1., 1.2, -1.5, 1.5);
   // Draw this function in pad2 with color mesh option
   pad2->cd();
   pad2->SetLogz();
   f2clone->Draw("glsurf1");
 
   // add axis titles. The titles are set on the intermediate
   // histogram used for visualisation. We must force this histogram
   // to be created, then force the redrawing of the two pads
   pad2->Update();
   f2->GetHistogram()->GetXaxis()->SetTitle("x title");
   f2->GetHistogram()->GetYaxis()->SetTitle("y title");
   f2->GetHistogram()->GetXaxis()->SetTitleOffset(1.4);
   f2->GetHistogram()->GetYaxis()->SetTitleOffset(1.4);
   f2clone->GetHistogram()->GetXaxis()->SetTitle("x title");
   f2clone->GetHistogram()->GetYaxis()->SetTitle("y title");
   f2clone->GetHistogram()->GetXaxis()->SetTitleOffset(1.4);
   f2clone->GetHistogram()->GetYaxis()->SetTitleOffset(1.4);
   pad1->Modified();
   pad2->Modified();
 
   TCanvas *c2 = new TCanvas("glc2", "Surfaces Drawing Options with gl", 700, 10, 700, 700);
   c2->SetFillColor(42);
   gStyle->SetFrameFillColor(42);
 
   c2->Divide(2, 2);
 
   c2->cd(1);
   TF2 *fun1 = new TF2("fun1", "1000*((sin(x)/x)*(sin(y)/y))+200", -6., 6., -6., 6.);
   fun1->SetNpx(30);
   fun1->SetNpy(30);
   fun1->SetFillColor(kGreen);
   fun1->Draw("glsurf3");
 
   c2->cd(2);
   TF2 *fun2 = new TF2("fun2", "cos(y)*sin(x)+cos(x)*sin(y)", -6., 6., -6., 6.);
   fun2->Draw("glsurf1cyl");
 
   c2->cd(3);
   TF2 *fun3 = new TF2("fun3", "sin(x) / x * cos(y) * y", -6., 6., -6., 6.);
   fun3->Draw("glsurfpol");
 
   c2->cd(4);
   TF3 *fun4 = new TF3("fun4", "sin(x * x + y * y + z * z - 4)", -2.5, 2.5, -2.5, 2.5, -2.5, 2.5);
   Int_t colInd = TColor::GetColor(1.f, 0.5f, 0.f);
   fun4->SetFillColor(colInd);
   fun4->Draw("gl"); // tf3 option
}