import ROOT from ROOT import ertool from ROOT import geoalgo from basictool.geoviewer import GeoViewer import matplotlib.pyplot as plt import random as rnd import time import datetime plt.ion() tottime = 0 # track: viewer = GeoViewer() viewer._use_box = False geoAlgo = geoalgo.GeoAlgo() line1 = geoalgo.LineSegment(105.727, -66.7871, 814.841, 105.438, -67.4666, 815.464) line2 = geoalgo.LineSegment(105.229, -67.5519, 815.806, 108.149, -62.1392, 810.246) pt1 = geoalgo.Vector() pt2 = geoalgo.Vector() d = geoAlgo.SqDist(line1,line2,pt1,pt2) line3 = geoalgo.LineSegment(pt1,pt2) viewer.add(line1,'ln1','b') viewer.add(line2,'ln2','r') viewer.add(line3,'ln3','k') print 'distance is %.02f'%d viewer.show() try: counter = input('Hit ENTER when you are ready to close the viewer') except:
import ROOT from ROOT import ertool from ROOT import geoalgo from basictool.geoviewer import GeoViewer import matplotlib.pyplot as plt import random as rnd import time import datetime plt.ion() tottime = 0 # track: viewer = GeoViewer() viewer._use_box = False geoAlgo = geoalgo.GeoAlgo() line1 = geoalgo.LineSegment(0, 0, 0, 0, 0, 1) line2 = geoalgo.LineSegment(1, 0, 0, 1, 0, 1) pt1 = geoalgo.Vector() pt2 = geoalgo.Vector() d = geoAlgo.SqDist(line1, line2, pt1, pt2) line3 = geoalgo.LineSegment(pt1, pt2) viewer.add(line1, 'ln1', 'b') viewer.add(line2, 'ln2', 'r') viewer.add(line3, 'ln3', 'k') print 'distance is %.04f' % d viewer.show() try: counter = input('Hit ENTER when you are ready to close the viewer') except:
import ROOT from ROOT import ertool from ROOT import geoalgo from basictool.geoviewer import GeoViewer import matplotlib.pyplot as plt import random as rnd import time import datetime plt.ion() tottime = 0 # track: viewer2 = GeoViewer() viewer2._use_box = False geoAlgo = geoalgo.GeoAlgo() ntries = 1000 npoints = 800 #for x in xrange(ntries): track = ertool.Track() #viewer2._fig.clear() #for n in xrange(npoints): #time.sleep(0.01) point = geoalgo.Vector(rnd.random(),rnd.random(),rnd.random()) n = 0 viewer2.add(point,'pt %i'%n,'k') track.push_back(point) viewer2.add(point,'pt %i'%n,'k') track.push_back(point) viewer2.add(point,'pt %i'%n,'k')
import ROOT from ROOT import ertool from ROOT import geoalgo from basictool.geoviewer import GeoViewer import matplotlib.pyplot as plt import random as rnd import time import datetime plt.ion() tottime = 0 # track: viewer2 = GeoViewer() viewer2._use_box = False geoAlgo = geoalgo.GeoAlgo() ntries = 1000 npoints = 800 #for x in xrange(ntries): track = ertool.Track() #viewer2._fig.clear() #for n in xrange(npoints): #time.sleep(0.01) point = geoalgo.Vector(rnd.random(), rnd.random(), rnd.random()) n = 0 viewer2.add(point, 'pt %i' % n, 'k') track.push_back(point) viewer2.add(point, 'pt %i' % n, 'k') track.push_back(point) viewer2.add(point, 'pt %i' % n, 'k')