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)
