void tornado()
{
   gBenchmark->Start("tornado");

   double PI = 3.141592653;
   int d = 16;
   int numberOfPoints = 200;
   int numberOfCircles = 40;

   // create and open a canvas
   TCanvas *sky = new TCanvas("sky", "Tornado", 300, 10, 700, 500);
   sky->SetFillColor(14);

   // creating view
   TView *view = TView::CreateView(1, 0, 0);
   float range = numberOfCircles * d;
   view->SetRange(0, 0, 0, 4.0 * range, 2.0 * range, range);

   for (int j = d; j < numberOfCircles * d; j += d) {

      // create a PolyMarker3D
      TPolyMarker3D *pm3d = new TPolyMarker3D(numberOfPoints);

      float x, y, z;

      // set points
      for (int i = 1; i < numberOfPoints; i++) {
         float csin = sin(2 * PI / (double)numberOfPoints * (double)i) + 1;
         float ccos = cos(2 * PI / (double)numberOfPoints * (double)i) + 1;
         float esin = sin(2 * PI / (double)(numberOfCircles * d) * (double)j) + 1;
         x = j * (csin + esin);
         y = j * ccos;
         z = j;
         pm3d->SetPoint(i, x, y, z);
      }

      // set marker size, color & style
      pm3d->SetMarkerSize(1);
      pm3d->SetMarkerColor(2 + (d == (j & d)));
      pm3d->SetMarkerStyle(3);

      // draw
      pm3d->Draw();
   }

   gBenchmark->Show("tornado");

   Float_t ct = gBenchmark->GetCpuTime("tornado");

   TPaveText *text = new TPaveText(0.1, 0.81, 0.9, 0.97);
   text->SetFillColor(42);
   text->AddText("ROOT example: tornado.C");
   text->AddText(Form("Execution time: %g sec.", ct));
   text->Draw();
   sky->Update();
}