3 from ROOT
import TCanvas, TView, TPolyMarker3D, TPaveText
4 from ROOT
import gROOT, gBenchmark
5 from math
import cos, sin, pi
9 gBenchmark.Start(
'tornado' )
16 sky = TCanvas(
'sky',
'Tornado', 300, 10, 700, 500 )
17 sky.SetFillColor( 14 )
21 rng = numberOfCircles * d
22 view.SetRange( 0, 0, 0, 4.0*rng, 2.0*rng, rng )
25 for j
in range( d, numberOfCircles * d, d ):
28 pm3d = TPolyMarker3D( numberOfPoints )
31 for i
in xrange( 1, numberOfPoints ) :
32 csin =
sin( 2*pi / numberOfPoints * i ) + 1
33 ccos =
cos( 2*pi / numberOfPoints * i ) + 1
34 esin =
sin( 2*pi / (numberOfCircles*d) * j ) + 1
35 x = j * ( csin + esin );
38 pm3d.SetPoint( i, x, y, z );
41 pm3d.SetMarkerSize( 1 )
42 pm3d.SetMarkerColor( 2 + ( d == ( j & d ) ) )
43 pm3d.SetMarkerStyle( 3 )
49 polymarkers.append( pm3d )
51 gBenchmark.Show(
'tornado' )
53 ct = gBenchmark.GetCpuTime(
'tornado' )
54 timeStr =
'Execution time: %g sec.' % ct
56 text = TPaveText( 0.1, 0.81, 0.9, 0.97 )
57 text.SetFillColor( 42 )
58 text.AddText(
'PyROOT example: tornado.py' )
59 text.AddText( timeStr )
static TView * CreateView(Int_t system=1, const Double_t *rmin=0, const Double_t *rmax=0)
Create a concrete default 3-d view via the plug-in manager.