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 )
def printTree( tree ):
tree.movePointerToTheTop()
printParticle( " ", "top particle", tree.currentParticle() )
print " vertex: ", tree.currentDecayVertex().position()
print " children:"
children = tree.finalStateParticles()
for child in children:
printParticle( " ","- ", child )
# --------------------------------------------------- Unconstrained Vertex Fitting
print "\nreconstructing the decay chain"
vtree = treefactory.useVertexFitter( particles )
printTree( vtree )
ZParticle = vtree.topParticle()
print "\nguessed value for Zmass = ", ZParticle.fullstate().m()
# --------------------------------------------------- TwoTrackMass Constraint
print "\napplying the mass constraint Zmass = ", Zmass
ttmconstraint = RaveKConstraintBuilder().createTwoTrackMassKinematicConstraint( Zmass )
ctree = treefactory.useVertexFitter( particles, ttmconstraint )
printTree( ctree )
# --------------------------------------------------- MassKinematicConstraint
