Python RaveKTreeFactoryの例

コード例 #1

ファイル: ktest1.py プロジェクト: newtrino/vertigo

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()

コード例 #2

ファイル: test22.py プロジェクト: newtrino/vertigo

  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 )

コード例 #3

ファイル: test6.py プロジェクト: newtrino/vertigo

                    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: