c=RaveCovariance3D( width**2, 0,0,width**2,0,length**2 )
e=RaveEllipsoid3D( RavePoint3D(),c )
LoopSettings.Instance().setVerbosity(0)
ravefactory=RaveVertexFactory ( RaveConstantMagneticField(0.,0.,4.), RaveVacuumPropagator(), e, "avr-primcut:3.0" )
eventfactory=EventFactory ( "bjets.170.1.txt.gz" )
event=eventfactory.next()
for simvtx in event.simVertices():
print simvtx
print len(event.jets()),"jets in event."
print len(event.tracks()),"tracks in event."
primaries=RaveTrackContainer()
first=True
secondaries=None
for jet in event.jets():
tracks=jet.tracks()
print len(tracks), "tracks in jet."
if not first: # put all tracks but of "first jet" in "primaries"
for track in tracks:
primaries.append ( track )
else:
secondaries=tracks
first=False
vertices=ravefactory.create ( event.tracks(), True )
LoopSettings.Instance().setVerbosity(0)
ravefactory=RaveVertexFactory ( RaveConstantMagneticField(0.,0.,4.), RaveVacuumPropagator() )
eventfactory=EventFactory ( "input8.txt" )
event=eventfactory.next()
if doAvr:
vertices=ravefactory.create ( event.tracks(), "avr-primcut:2.0-seccut:3.0-smoothing:1-minweight:0.5" )
print "[avr]"
for vertex in vertices:
pprint ( vertex )
print
if doAvf:
print "[avf]"
vertex=ravefactory.create ( event.tracks(), "avf-sigmacut:2.0-smoothing:1-ratio:0.0078125" )[0]
newtracks=RaveTrackContainer()
pprint ( vertex )
wtracks = vertex.weightedTracks()
for t in wtracks:
if t[0] < 0.5: newtracks.push_back ( t[1] )
vertex=ravefactory.create ( newtracks, "avf-sigmacut:3.0-smoothing:1-ratio:0.25" )[0]
pprint ( vertex )
