print "Simple btagging test, but now with beamspot constraint!" from vertigo import RaveVertexFactory, RaveFlavorTagFactory, RaveConstantMagneticField,\ RaveVacuumPropagator, EventFactory, VertigoToCmsObjects,\ RavePoint3D, LoopSettings, RaveTrackContainer, VertexMetaFactory_Instance, \ RaveEllipsoid3D, RaveCovariance3D, LoopSettings_Instance LoopSettings_Instance().setVerbosity(255) ravefactory=VertexMetaFactory_Instance().get("dummy") method="avr-primcut:2.5-primT:4096-primr:0.125-secT:1024-secr:0.25-smoothing:0-weightthreshold:0.0-seccut:2.3" ravefactory.setDefaultMethod ( method ) # beamspot=RaveEllipsoid3D() beamspot=RaveEllipsoid3D( RavePoint3D(), RaveCovariance3D( .0015**2,0.,0.,.0015**2,0.,5.3**2 ) ) ravefactory.setBeamSpot ( beamspot ) if False: print "bs=",beamspot.isValid() print "Has Beamspot?",ravefactory.hasBeamspot() print "dzz",ravefactory.getBeamspot().covariance().dzz() print "mz=",ravefactory.getMagneticField().inTesla ( RavePoint3D() ).z() btagfactory=RaveFlavorTagFactory ( ravefactory ) set=LoopSettings.Instance() set.setVerbosity(0) eventfactory=EventFactory ( "bb.1.txt" )