def visualizeMagFields(): checkMagFields.run()
def visualizeMagFields(): checkMagFields.run()
rtdb = run.GetRuntimeDb() # -----Create geometry---------------------------------------------- modules = shipDet_conf.configure(run, ShipGeo) import geomGeant4 if hasattr(ShipGeo.Bfield, "fieldMap"): fieldMaker = geomGeant4.addVMCFields(ShipGeo, '', True, withVirtualMC=False) else: print "no fieldmap given, geofile too old, not anymore support" exit(-1) sGeo = fgeo.FAIRGeom geoMat = ROOT.genfit.TGeoMaterialInterface() ROOT.genfit.MaterialEffects.getInstance().init(geoMat) bfield = ROOT.genfit.FairShipFields() bfield.setField(fieldMaker.getGlobalField()) fM = ROOT.genfit.FieldManager.getInstance() fM.init(bfield) volDict = {} i = 0 for x in ROOT.gGeoManager.GetListOfVolumes(): volDict[i] = x.GetName() i += 1 import checkMagFields checkMagFields.run(-6815) # checkMagFields.run(-5894)