print "running vertigo/tests/ktest.py ..." from vertigo import LoopSettings, EventFactory, RaveConstantMagneticField, RaveVacuumPropagator, RaveKTransientTrackParticle, RaveVertexFactory, RaveKTreeFactory, RaveKParticleContainer, RaveKTreeContainer, RaveKConstraintBuilder LoopSettings.Instance().setVerbosity(100) muonmass = 0.1056583 muonsigma = 0.0000000001 kchi2 = 0. kndf = 0. magneticfield = RaveConstantMagneticField( 0., 0., 4. ) eventfactory = EventFactory ( "gun:kinematics" ) vertexfactory = RaveVertexFactory ( magneticfield, RaveVacuumPropagator() ) treefactory = RaveKTreeFactory( magneticfield ) muonmass = 0.1056583 muonsigma = 0.0000000001 Zmass = 91.187 Zsigma = 2.4952 #print "\nskipping the first shot" #eventfactory.skip(2) #print "using the following shot" event = eventfactory.next() tracks = event.tracks() print "\nfound ", len( tracks ), " generated tracks" eparticles = event.particles()
for child in children: printParticle( " ","- ", child ) # LoopSettings.Instance().setVerbosity(100) muonmass = 0.1056583 muonsigma = 0.0000000001 kchi2 = 0. kndf = 0. magneticfield = RaveConstantMagneticField( 0., 0., 4. ) beamspot=RaveEllipsoid3D( RavePoint3D(), RaveCovariance3D( .0015**2,0.,0.,.0015**2,0.,5.3**2 ) ) # eventfactory = EventFactory ( "gun:kinematics" ) eventfactory = EventFactory ( "file:input6.txt" ) vertexfactory = RaveVertexFactory ( magneticfield, RaveVacuumPropagator(), beamspot ) treefactory = RaveKTreeFactory( magneticfield, RaveVacuumPropagator(), beamspot ) muonmass = 0.1056583 muonsigma = 0.0000000001 Zmass = 91.187 Zsigma = 2.4952 event = eventfactory.next() tracks=event.tracks() particles=event.particles() print "\nreconstructing the vertex using method 'kalman' with", len(tracks),"tracks" vtree = treefactory.useVertexFitter( particles ) printTree ( vtree )
RaveVacuumPropagator, RaveKTransientTrackParticle, \ RaveVertexFactory, RaveKTreeFactory, RaveKParticleContainer, \ RaveKTreeContainer, RaveKConstraintBuilder LoopSettings.Instance().setVerbosity(100) muonmass = 0.1056583 muonsigma = 0.0000000001 kchi2 = 0. kndf = 0. magneticfield = RaveConstantMagneticField( 0., 0., 4. ) # eventfactory = EventFactory ( "file:input6.txt" ) eventfactory = EventFactory ( "gun:kinematics" ) vertexfactory = RaveVertexFactory ( magneticfield, RaveVacuumPropagator() ) treefactory = RaveKTreeFactory( magneticfield ) muonmass = 0.1056583 muonsigma = 0.0000000001 Zmass = 91.187 Zsigma = 2.4952 # print "\nskipping the first shot" # eventfactory.skip(1) # print "using the following shot" event = eventfactory.next() tracks = event.tracks() print "\nfound ", len( tracks ), " generated tracks" for e in tracks: