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 )
