Logo ROOT   6.08/07
Reference Guide
tornado.C File Reference

Detailed Description

View in nbviewer Open in SWAN Show 3-d polymarker.

pict1_tornado.C.png
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();
}
char timeStr[60];
gBenchmark->Show("tornado");
Float_t ct = gBenchmark->GetCpuTime("tornado");
sprintf( timeStr, "Execution time: %g sec.", ct);
TPaveText *text = new TPaveText( 0.1, 0.81, 0.9, 0.97 );
text->SetFillColor( 42 );
text->AddText("ROOT example: tornado.C");
text->AddText(timeStr);
text->Draw();
sky->Update();
}
Author
Rene Brun

Definition in file tornado.C.