def __init__(self,
name='InDetTrackingGeometryBuilder',
namePrefix='',
setLayerAssociation=True,
buildTrtStrawLayers=False):
# get the ToolSvc
from AthenaCommon.AppMgr import ToolSvc, ServiceMgr
# the Detector flags to be imported
from TrkDetDescrSvc.TrkDetDescrJobProperties import TrkDetFlags
# beampipe
from InDetTrackingGeometry.InDetTrackingGeometryConf import InDet__BeamPipeBuilder
BeamPipeBuilder = InDet__BeamPipeBuilder(name=namePrefix +
'BeamPipeBuilder')
BeamPipeBuilder.BeamPipeRadius = TrkDetFlags.BeamPipeRadius()
BeamPipeBuilder.BeamPipeThickness = TrkDetFlags.BeamPipeThickness()
BeamPipeBuilder.BeamPipeMaterialBinsZ = TrkDetFlags.BeamPipeLayerMaterialBinsZ(
)
BeamPipeBuilder.OutputLevel = TrkDetFlags.InDetBuildingOutputLevel()
ToolSvc += BeamPipeBuilder
# the layer builders
layerbuilders = []
binnings = []
colors = []
# @TODO put DetFlags isOn here, but make sure that this is really when the detector is not built
# PIXEL building
from InDetTrackingGeometry.InDetTrackingGeometryConf import InDet__SiLayerBuilder
PixelLayerBuilder = InDet__SiLayerBuilder(name=namePrefix +
'PixelLayerBuilder')
PixelLayerBuilder.PixelCase = True
PixelLayerBuilder.Identification = 'Pixel'
PixelLayerBuilder.SiDetManagerLocation = 'Pixel'
# additionall layers - handle with care !
PixelLayerBuilder.BarrelAdditionalLayerRadii = [130] # The PST
PixelLayerBuilder.BarrelAdditionalLayerType = [
0
] # -- will shift volume boundary to PST
PixelLayerBuilder.EndcapAdditionalLayerPositionsZ = [-1900.,
1900.] # DBM
PixelLayerBuilder.EndcapAdditionalLayerType = [1, 1] # DBM
# Pixel barrel specifications
PixelLayerBuilder.BarrelLayerBinsZ = TrkDetFlags.PixelBarrelLayerMaterialBinsZ(
)
PixelLayerBuilder.BarrelLayerBinsPhi = TrkDetFlags.PixelBarrelLayerMaterialBinsPhi(
)
PixelLayerBuilder.EndcapLayerBinsR = TrkDetFlags.PixelEndcapLayerMaterialBinsR(
)
PixelLayerBuilder.EndcapLayerBinsPhi = TrkDetFlags.PixelEndcapLayerMaterialBinsPhi(
)
# set the layer association
PixelLayerBuilder.SetLayerAssociation = setLayerAssociation
# output level
PixelLayerBuilder.OutputLevel = TrkDetFlags.PixelBuildingOutputLevel()
# the binning type of the layers
PixelLayerBinning = 2
# add it to tool service
ToolSvc += PixelLayerBuilder
# put them to the caches
layerbuilders += [PixelLayerBuilder]
binnings += [PixelLayerBinning]
colors += [3]
# SCT building
SCT_LayerBuilder = InDet__SiLayerBuilder(name=namePrefix +
'SCT_LayerBuilder')
SCT_LayerBuilder.PixelCase = False
SCT_LayerBuilder.Identification = 'SCT'
SCT_LayerBuilder.SiDetManagerLocation = 'SCT'
# additionall layers - handle with care !
SCT_LayerBuilder.EndcapAdditionalLayerPositionsZ = [-2850, 2850]
SCT_LayerBuilder.EndcapAdditionalLayerType = [0, 0]
# SCT barrel specifications
SCT_LayerBuilder.BarrelLayerBinsZ = TrkDetFlags.SCT_BarrelLayerMaterialBinsZ(
)
SCT_LayerBuilder.BarrelLayerBinsPhi = TrkDetFlags.SCT_BarrelLayerMaterialBinsPhi(
)
# SCT endcap specifications
SCT_LayerBuilder.EndcapLayerBinsR = TrkDetFlags.SCT_EndcapLayerMaterialBinsR(
)
SCT_LayerBuilder.EndcapLayerBinsPhi = TrkDetFlags.SCT_EndcapLayerMaterialBinsPhi(
)
SCT_LayerBuilder.EndcapComplexRingBinning = TrkDetFlags.SCT_EndcapLayerDynamicRings(
)
# set the layer association
SCT_LayerBuilder.SetLayerAssociation = setLayerAssociation
# output level
SCT_LayerBuilder.OutputLevel = TrkDetFlags.SCT_BuildingOutputLevel()
# the binning type of the layer
SCT_LayerBinning = 2
# SCT -> ToolSvc
ToolSvc += SCT_LayerBuilder
# put them to the caches
layerbuilders += [SCT_LayerBuilder]
binnings += [SCT_LayerBinning]
colors += [4]
from InDetTrackingGeometry.InDetTrackingGeometryConf import InDet__TRT_LayerBuilder
TRT_LayerBuilder = InDet__TRT_LayerBuilder(name=namePrefix +
'TRT_LayerBuilder')
# TRT barrel specifications
TRT_LayerBuilder.ModelBarrelLayers = TrkDetFlags.TRT_BarrelModelLayers(
)
TRT_LayerBuilder.BarrelLayerBinsZ = TrkDetFlags.TRT_BarrelLayerMaterialBinsZ(
)
TRT_LayerBuilder.BarrelLayerBinsPhi = TrkDetFlags.TRT_BarrelLayerMaterialBinsPhi(
)
# SCT endcap specifications
TRT_LayerBuilder.ModelEndcapLayers = TrkDetFlags.TRT_EndcapModelLayers(
)
TRT_LayerBuilder.EndcapLayerBinsR = TrkDetFlags.TRT_EndcapLayerMaterialBinsR(
)
TRT_LayerBuilder.EndcapLayerBinsPhi = TrkDetFlags.TRT_EndcapLayerMaterialBinsPhi(
)
# set the binning from bi-aequidistant to arbitrary for complex TRT volumes
TRT_LayerBinning = 1
if buildTrtStrawLayers or TrkDetFlags.TRT_BuildStrawLayers():
TRT_LayerBinning = 2
TRT_LayerBuilder.ModelLayersOnly = False
# output level
TRT_LayerBuilder.OutputLevel = TrkDetFlags.TRT_BuildingOutputLevel()
# TRT -> ToolSvc
ToolSvc += TRT_LayerBuilder
# put them to the caches
layerbuilders += [TRT_LayerBuilder]
binnings += [TRT_LayerBinning]
colors += [5]
# helpers for the InDetTrackingGeometry Builder : layer array creator
from TrkDetDescrTools.TrkDetDescrToolsConf import Trk__LayerArrayCreator
InDetLayerArrayCreator = Trk__LayerArrayCreator(
name='InDetLayerArrayCreator')
InDetLayerArrayCreator.EmptyLayerMode = 2 # deletes empty material layers from arrays
InDetLayerArrayCreator.OutputLevel = TrkDetFlags.InDetBuildingOutputLevel(
)
# add to ToolSvc
ToolSvc += InDetLayerArrayCreator
# helpers for the InDetTrackingGeometry Builder : volume array creator
from TrkDetDescrTools.TrkDetDescrToolsConf import Trk__TrackingVolumeArrayCreator
InDetTrackingVolumeArrayCreator = Trk__TrackingVolumeArrayCreator(
name='InDetTrackingVolumeArrayCreator')
InDetTrackingVolumeArrayCreator.OutputLevel = TrkDetFlags.InDetBuildingOutputLevel(
)
# add to ToolSvc
ToolSvc += InDetTrackingVolumeArrayCreator
# helpers for the InDetTrackingGeometry Builder : tracking voluem helper for glueing
from TrkDetDescrTools.TrkDetDescrToolsConf import Trk__TrackingVolumeHelper
InDetTrackingVolumeHelper = Trk__TrackingVolumeHelper(
name='InDetTrackingVolumeHelper')
InDetTrackingVolumeHelper.OutputLevel = TrkDetFlags.InDetBuildingOutputLevel(
)
# the material bins
InDetTrackingVolumeHelper.BarrelLayerBinsZ = TrkDetFlags.InDetPassiveLayerMaterialBinsRz(
)
InDetTrackingVolumeHelper.BarrelLayerBinsPhi = TrkDetFlags.InDetPassiveLayerMaterialBinsPhi(
)
InDetTrackingVolumeHelper.EndcapLayerBinsR = TrkDetFlags.InDetPassiveLayerMaterialBinsRz(
)
InDetTrackingVolumeHelper.EndcapLayerBinsPhi = TrkDetFlags.InDetPassiveLayerMaterialBinsPhi(
)
# add to ToolSvc
ToolSvc += InDetTrackingVolumeHelper
# helpers for the InDetTrackingGeometry Builder : cylinder volume creator
from TrkDetDescrTools.TrkDetDescrToolsConf import Trk__CylinderVolumeCreator
InDetCylinderVolumeCreator = Trk__CylinderVolumeCreator(
name='InDetCylinderVolumeCreator')
# give him the layer array creator
InDetCylinderVolumeCreator.LayerArrayCreator = InDetLayerArrayCreator
InDetCylinderVolumeCreator.TrackingVolumeArrayCreator = InDetTrackingVolumeArrayCreator
InDetCylinderVolumeCreator.TrackingVolumeHelper = InDetTrackingVolumeHelper
# specifiy the binning, passive layers, entry layers
InDetCylinderVolumeCreator.PassiveLayerThickness = TrkDetFlags.InDetPassiveLayerThickness(
)
InDetCylinderVolumeCreator.PassiveLayerBinsPhi = TrkDetFlags.InDetPassiveLayerMaterialBinsPhi(
)
InDetCylinderVolumeCreator.PassiveLayerBinsRZ = TrkDetFlags.InDetPassiveLayerMaterialBinsRz(
)
# output level
InDetCylinderVolumeCreator.OutputLevel = TrkDetFlags.InDetBuildingOutputLevel(
)
# add to ToolSvc
ToolSvc += InDetCylinderVolumeCreator
# the envelope definition service
from SubDetectorEnvelopes.SubDetectorEnvelopesConf import DetDescrDBEnvelopeSvc
AtlasEnvelopeSvc = DetDescrDBEnvelopeSvc(
'AtlasEnvelopeSvcDefinitionSvc')
# set the output level for the Envelope service
AtlasEnvelopeSvc.OutputLevel = TrkDetFlags.InDetBuildingOutputLevel()
# add to SvcMgr
ServiceMgr += AtlasEnvelopeSvc
# the tracking geometry builder
InDet__RobustTrackingGeometryBuilder.__init__(self,namePrefix+name,\
BeamPipeBuilder = BeamPipeBuilder,\
LayerBuilders = layerbuilders,
LayerBinningType = binnings,
ColorCodes = colors,
BarrelEntryLayers = [ 2, 2, 2 ],
EndcapEntryLayers = [ 1, 0, 1 ],
EnvelopeDefinitionSvc = AtlasEnvelopeSvc,
VolumeEnclosureDiscPositions = [ 3000., 3450. ],
TrackingVolumeCreator = InDetCylinderVolumeCreator,
LayerArrayCreator = InDetLayerArrayCreator,
BuildBoundaryLayers = TrkDetFlags.InDetBuildMaterialBoundaries(),
ReplaceAllJointBoundaries = TrkDetFlags.InDetBuildJointBoundaries(),
OutputLevel = TrkDetFlags.InDetBuildingOutputLevel(),
ExitVolumeName = TrkDetFlags.InDetContainerName(),
MagneticFieldMode = TrkDetFlags.MagneticFieldMode())
def __init__(self,
name='InDetTrackingGeometryBuilder',
namePrefix='',
setLayerAssociation=True,
buildTrtStrawLayers=False,
buildCustomPixel=False,
buildCustomSCT=False,
buildCustomTRT=False,
createXML=True,
dictionaryFileName='InDetTrackingGeometry.xml'):
# get the ToolSvc
from AthenaCommon.AppMgr import ToolSvc, ServiceMgr
# the Detector flags to be imported
from TrkDetDescrSvc.TrkDetDescrJobProperties import TrkDetFlags
from ISF_FatrasDetDescrTools.ISF_FatrasDetDescrToolsConf import iFatras__InputLayerMaterialProvider
ISF_InputLayerMaterialProvider = iFatras__InputLayerMaterialProvider(
name=namePrefix + 'ISF_InputLayerMaterialProvider')
#ISF_InputLayerMaterialProvider.OutputLevel = 1
ToolSvc += ISF_InputLayerMaterialProvider
print ISF_InputLayerMaterialProvider
# beampipe
from InDetTrackingGeometry.InDetTrackingGeometryConf import InDet__BeamPipeBuilder
BeamPipeBuilder = InDet__BeamPipeBuilder(name=namePrefix +
'BeamPipeBuilder')
BeamPipeBuilder.BeamPipeRadius = TrkDetFlags.BeamPipeRadius()
BeamPipeBuilder.BeamPipeThickness = TrkDetFlags.BeamPipeThickness()
BeamPipeBuilder.BeamPipeMaterialBinsZ = TrkDetFlags.BeamPipeLayerMaterialBinsZ(
)
BeamPipeBuilder.OutputLevel = TrkDetFlags.InDetBuildingOutputLevel()
ToolSvc += BeamPipeBuilder
print BeamPipeBuilder
# the layer providers
from TrkDetDescrTools.TrkDetDescrToolsConf import Trk__LayerProvider
BeamPipeProvider = Trk__LayerProvider(name=namePrefix +
'BeamPipeProvider')
BeamPipeProvider.LayerBuilder = BeamPipeBuilder
ToolSvc += BeamPipeProvider
# binning of the beam pipe
BeamPipeBinning = 2
# the layer builders
layerProviders = [BeamPipeProvider]
binningsCenter = [BeamPipeBinning]
binningsEndcap = [BeamPipeBinning]
colors = [2]
IBLLayerBuilder = None
PixelLayerBuilder = None
if buildCustomPixel:
from ISF_FatrasDetDescrTools.ISF_FatrasDetDescrToolsConf import iFatras__PlanarDetLayerBuilder
## EXT IBL building
#IBLLayerBuilder = iFatras__PlanarDetLayerBuilder(name=namePrefix+'IBLLayerBuilder')
#IBLLayerBuilder.PixelCase = True
#IBLLayerBuilder.Identification = 'Pixel'
#IBLLayerBuilder.InputLayerMaterialProvider = ISF_InputLayerMaterialProvider
## Assign custom material
#IBLLayerBuilder.BarrelLayerBinsZ = TrkDetFlags.PixelBarrelLayerMaterialBinsZ()
#IBLLayerBuilder.BarrelLayerBinsPhi = TrkDetFlags.PixelBarrelLayerMaterialBinsPhi()
#IBLLayerBuilder.EndcapLayerBinsR = TrkDetFlags.PixelEndcapLayerMaterialBinsR()
#IBLLayerBuilder.EndcapLayerBinsPhi = TrkDetFlags.PixelEndcapLayerMaterialBinsPhi()
#IBLLayerBuilder.CustomMaterial = False
#IBLLayerBuilder.CustomMaterialThickness = 0.250
#IBLLayerBuilder.CustomMaterialX0 = 8.333
#IBLLayerBuilder.CustomMaterialL0 = 100
#IBLLayerBuilder.CustomMaterialA = 14
#IBLLayerBuilder.CustomMaterialZ = 28.0855
#IBLLayerBuilder.CustomMaterialRho = 0.00233
## INNER PIXEL BARREL
#INPB_NR = 1
#IBLLayerBuilder.BarrelLayers = INPB_NR
#IBLLayerBuilder.LayersZsectors = [40] # [58, 58]
#IBLLayerBuilder.LayerPhiSectors = [16]
#IBLLayerBuilder.LayerTilt = [-14.0 ] * INPB_NR #degrees
#IBLLayerBuilder.LayerMinPhi = [-180.0 ] * INPB_NR #degrees
#IBLLayerBuilder.LayerMaxPhi = [ 180.0 ] * INPB_NR #degrees
#IBLLayerBuilder.LayerMinZ = [-40.5 * z/2. for z in IBLLayerBuilder.LayersZsectors] #mm
#IBLLayerBuilder.LayerMaxZ = [40.5 * z/2. for z in IBLLayerBuilder.LayersZsectors] #mm
#IBLLayerBuilder.LayerRadius = [39.]
#IBLLayerBuilder.LayerThickness = [0.150 + 0.150 + 0.182] * INPB_NR # chip + sensor + hybrid # CHECK HYBRID MATERIAL
#IBLLayerBuilder.LayerLengthY = [ 40.5 ]
#IBLLayerBuilder.LayerLengthXmin = [ 16.81 ]
#IBLLayerBuilder.LayerPitchX = [0.050] * INPB_NR
#IBLLayerBuilder.LayerPitchY = [0.050] * INPB_NR
#IBLLayerBuilder.LayerRotation = [0.] * INPB_NR
#IBLLayerBuilder.AdditionalLayerRadius = []
#IBLLayerBuilder.EndcapDiscs = 0
## set the layer association
#IBLLayerBuilder.SetLayerAssociation = setLayerAssociation
#IBLLayerBuilder.OutputLevel = TrkDetFlags.PixelBuildingOutputLevel()
## the binning type of the layers
#IBLLayerBinning = 2
## add it to tool service
#ToolSvc += IBLLayerBuilder
#print IBLLayerBuilder
## the layer providers
#IBLLayerProvider = Trk__LayerProvider(name=namePrefix+'IBLLayerProvider')
#IBLLayerProvider.LayerBuilder = IBLLayerBuilder
#ToolSvc += IBLLayerProvider
## put them to the caches
#layerProviders += [ IBLLayerProvider ]
#binningsCenter += [ IBLLayerBinning ]
#binningsEndcap += [ IBLLayerBinning ]
#colors += [ 3 ]
# PIXEL building
PixelLayerBuilder = iFatras__PlanarDetLayerBuilder(
name=namePrefix + 'PixelLayerBuilder')
PixelLayerBuilder.PixelCase = True
PixelLayerBuilder.Identification = 'Pixel'
PixelLayerBuilder.InputLayerMaterialProvider = ISF_InputLayerMaterialProvider
# Pixel barrel specifications
# Noemi's configuration
# Assign custom material
PixelLayerBuilder.BarrelLayerBinsZ = TrkDetFlags.PixelBarrelLayerMaterialBinsZ(
)
PixelLayerBuilder.BarrelLayerBinsPhi = 100 #TrkDetFlags.PixelBarrelLayerMaterialBinsPhi()
PixelLayerBuilder.EndcapLayerBinsR = TrkDetFlags.PixelEndcapLayerMaterialBinsR(
)
PixelLayerBuilder.EndcapLayerBinsPhi = TrkDetFlags.PixelEndcapLayerMaterialBinsPhi(
)
PixelLayerBuilder.CustomMaterial = True
PixelLayerBuilder.CustomMaterialThickness = 0.250
PixelLayerBuilder.CustomMaterialX0 = 8.333
PixelLayerBuilder.CustomMaterialL0 = 100
PixelLayerBuilder.CustomMaterialA = 14
PixelLayerBuilder.CustomMaterialZ = 28.0855
PixelLayerBuilder.CustomMaterialRho = 0.00233
## BARREL
PixelLayerBuilder.BarrelLayers = 3
if IBLLayerBuilder:
PixelLayerBuilder.LayerOffset = IBLLayerBuilder.BarrelLayers
PixelLayerBuilder.LayersZsectors = [13, 13, 13]
PixelLayerBuilder.LayerPhiSectors = [22, 38, 52]
PixelLayerBuilder.LayerTilt = [-20.0, -20.0, -20.0] #degree
PixelLayerBuilder.LayerPhiOffset = [0., 0.,
360. / 52. * 0.4] #degree
PixelLayerBuilder.LayerMinPhi = [-180.0, -180.0, -180.0] #degree
PixelLayerBuilder.LayerMaxPhi = [180.0, 180.0, 180.0] #degree
PixelLayerBuilder.LayerMinZ = [-400.5, -400.5, -400.5]
PixelLayerBuilder.LayerMaxZ = [400.5, 400.5, 400.5]
PixelLayerBuilder.LayerRadius = [50.5, 88.5, 122.5]
PixelLayerBuilder.LayerThickness = [0.250, 0.250, 0.250]
PixelLayerBuilder.LayerLengthY = [60.8, 60.8, 60.8]
PixelLayerBuilder.LayerLengthXmin = [16.4, 16.4, 16.4]
PixelLayerBuilder.LayerPitchX = [0.010, 0.010, 0.010]
PixelLayerBuilder.LayerPitchY = [0.055, 0.055, 0.055]
#PixelLayerBuilder.LayerPitchY = [ 6.08, 6.08, 6.08 ]
#PixelLayerBuilder.LayerPitchX = [ 1.64, 1.64, 1.64 ]
PixelLayerBuilder.LayerRotation = [1., 1., 1.] #degree
PixelLayerBuilder.LayerSCTlike = False
PixelLayerBuilder.AdditionalLayerRadius = [] #45.5, 83.5, 117.5 ]
# ENDCAPS
PixelLayerBuilder.EndcapDiscs = 3
if IBLLayerBuilder:
PixelLayerBuilder.DiscOffset = IBLLayerBuilder.EndcapDiscs
PixelLayerBuilder.DiscPhiSectors = [[48], [48], [48]]
PixelLayerBuilder.DiscZpos = [
-649.0, -579.0, -493.5, 496.5, 581.0, 651.0
]
PixelLayerBuilder.DiscRingMinR = [[88.8], [88.8], [88.8]]
PixelLayerBuilder.DiscRingMaxR = [[149.6], [149.6], [149.6]]
PixelLayerBuilder.DiscMinPhi = [[-180.0], [-180.0], [-180.0]]
PixelLayerBuilder.DiscMaxPhi = [[180.0], [180.0], [180.0]]
PixelLayerBuilder.DiscThickness = [0.250, 0.250, 0.250]
PixelLayerBuilder.DiscLengthY = [[60.8], [60.8], [60.8]]
PixelLayerBuilder.DiscLengthXmin = [[16.4], [16.4], [16.4]]
PixelLayerBuilder.DiscPitchX = [[0.010], [0.010], [0.010]]
PixelLayerBuilder.DiscPitchY = [[0.055], [0.055], [0.055]]
#PixelLayerBuilder.DiscPitchY = [ [6.08], [6.08], [6.08] ]
#PixelLayerBuilder.DiscPitchX = [ [1.64], [1.64], [1.64] ]
PixelLayerBuilder.DiscSeparation = [[4.], [4.], [4.]]
PixelLayerBuilder.AdditionalDiscZpos = [
-1900., 1900.
] #-645.0, -575.0, -445.0, 445.0, 575.0, 645.0]
#ENDCAPS
#PixelLayerBuilder.EndcapDiscs = 7
#PixelLayerBuilder.DiscPhiSectors = [[48], [48], [48], [48], [48], [48], [48]]
#PixelLayerBuilder.DiscZpos = [-850, -800, -750, -700, -650.0, -580.0, -450.0, 450.0, 580.0, 650.0, 700, 750, 800, 850]
#PixelLayerBuilder.DiscRingMinR = [[88.8], [88.8], [88.8], [88.8], [88.8], [88.8], [88.8] ]
#PixelLayerBuilder.DiscRingMaxR = [[149.6], [149.6], [149.6], [149.6], [149.6], [149.6], [149.6]]
#PixelLayerBuilder.DiscMinPhi = [[-180.0], [-180.0], [-180.0], [-180.0], [-180.0], [-180.0], [-180.0]]
#PixelLayerBuilder.DiscMaxPhi = [[180.0], [180.0], [180.0], [180.0], [180.0], [180.0], [180.0]]
#PixelLayerBuilder.DiscThickness = [0.250, 0.250, 0.250, 0.250, 0.250, 0.250, 0.250]
#PixelLayerBuilder.DiscLengthY = [[60.8], [60.8], [60.8], [60.8], [60.8], [60.8], [60.8]]
#PixelLayerBuilder.DiscLengthXmin = [[16.4], [16.4], [16.4], [16.4], [16.4], [16.4], [16.4]]
#PixelLayerBuilder.DiscPitchX = [[0.8], [0.8], [0.8], [0.8], [0.8], [0.8], [0.8]]
#PixelLayerBuilder.DiscPitchY = [[0.8], [0.8], [0.8], [0.8], [0.8], [0.8], [0.8]]
##PixelLayerBuilder.DiscStereo = [[0], [0], [0]]
##PixelLayerBuilder.DiscStereoSeparation = [[0], [0], [0]]
#PixelLayerBuilder.AdditionalDiscZpos = [] #-645.0, -575.0, -445.0, 445.0, 575.0, 645.0]
# BARREL
#PixelLayerBuilder.BarrelLayers = 1
#PixelLayerBuilder.LayersZsectors = [ 13 ]
#PixelLayerBuilder.LayerPhiSectors = [ 22 ]
#PixelLayerBuilder.LayerTilt = [ -20.0 ] #degree
#PixelLayerBuilder.LayerMinPhi = [ -180.0 ] #degree
#PixelLayerBuilder.LayerMaxPhi = [ 180.0 ] #degree
#PixelLayerBuilder.LayerMinZ = [ -400.5 ]
#PixelLayerBuilder.LayerMaxZ = [ 400.5 ]
#PixelLayerBuilder.LayerRadius = [ 50.5 ]
#PixelLayerBuilder.LayerThickness = [ 0.250 ]
#PixelLayerBuilder.LayerLengthY = [ 60.8 ]
#PixelLayerBuilder.LayerLengthXmin = [ 16.4]
#PixelLayerBuilder.LayerPitchX = [ 0.010 ]
#PixelLayerBuilder.LayerPitchY = [ 0.055 ]
#PixelLayerBuilder.LayerRotation = [ 1 ] #degree
#PixelLayerBuilder.LayerStereo = [ 0 ]
#PixelLayerBuilder.AdditionalLayerRadius = [ 45.5 ] #, 83.5, 117.5 ]
#PixelLayerBuilder.BarrelLayers = 0
#PixelLayerBuilder.LayersZsectors = [ 13 , 13]
#PixelLayerBuilder.LayerPhiSectors = [ 22, 38 ]
#PixelLayerBuilder.LayerTilt = [ -20.0 , -20.0]
#PixelLayerBuilder.LayerMinPhi = [ -180.0, -180.0]
#PixelLayerBuilder.LayerMaxPhi = [ 180.0, 180.0]
#PixelLayerBuilder.LayerMinZ = [ -400.5, -400.5 ]
#PixelLayerBuilder.LayerMaxZ = [ 400.5, 400.5]
#PixelLayerBuilder.LayerRadius = [ 50.5, 88.5]
#PixelLayerBuilder.LayerThickness = [0.250, 0.250]
#PixelLayerBuilder.LayerPitchX = [ 0.4, 0.4 ]
#PixelLayerBuilder.LayerPitchY = [ 0.4, 0.4 ]
#PixelLayerBuilder.LayerLengthY = [ 60.8, 60.8 ]
#PixelLayerBuilder.LayerLengthXmin = [ 16.4, 16.4 ]
#PixelLayerBuilder.LayerRotation = [1, 1 ]
#PixelLayerBuilder.LayerStereo = [0, 0]
##PixelLayerBuilder.LayerStereoSeparation = [0]
#PixelLayerBuilder.AdditionalLayerRadius = []
#PixelLayerBuilder.BarrelLayers = 4
#PixelLayerBuilder.LayersZsectors = [ 13, 13, 13, 13]
#PixelLayerBuilder.LayerPhiSectors = [ 22, 38, 52 , 80]
#PixelLayerBuilder.LayerTilt = [ -20.0, -20.0, -20.0 , -20.0] #degree
#PixelLayerBuilder.LayerMinPhi = [ -180.0, -180.0, -180.0 , -180.0] #degree
#PixelLayerBuilder.LayerMaxPhi = [ 180.0, 180.0, 180.0 , 180.0] #degree
#PixelLayerBuilder.LayerMinZ = [ -400.5, -400.5, -400.5 , -400.5]
#PixelLayerBuilder.LayerMaxZ = [ 400.5, 400.5, 400.5, 400.5]
#PixelLayerBuilder.LayerRadius = [ 50.5, 85.5, 110.5, 135.5]
#PixelLayerBuilder.LayerThickness = [0.250, 0.250, 0.250, 0.250]
#PixelLayerBuilder.LayerLengthY = [ 60.8, 60.8, 60.8, 60.8]
#PixelLayerBuilder.LayerLengthXmin = [ 16.4 , 16.4 , 16.4, 16.4]
#PixelLayerBuilder.LayerPitchX = [ 0.4, 0.4, 0.4, 0.4 ]
#PixelLayerBuilder.LayerPitchY = [ 0.4, 0.4, 0.4, 0.4 ]
#PixelLayerBuilder.LayerRotation = [1, 1, 1, 1] #degree
#PixelLayerBuilder.LayerStereo = [0, 0, 0, 0]
##PixelLayerBuilder.LayerStereoSeparation = [0, 0, 0]
#PixelLayerBuilder.AdditionalLayerRadius = [] #45.5, 83.5, 117.5 ]
# set the layer association
PixelLayerBuilder.SetLayerAssociation = setLayerAssociation
# output level
#TrkDetFlags.InDetBuildingOutputLevel = 0
PixelLayerBuilder.OutputLevel = TrkDetFlags.PixelBuildingOutputLevel(
)
# the binning type of the layers
PixelLayerBinning = 2
# add it to tool service
ToolSvc += PixelLayerBuilder
print PixelLayerBuilder
# the layer providers
PixelLayerProvider = Trk__LayerProvider(name=namePrefix +
'PixelLayerProvider')
PixelLayerProvider.LayerBuilder = PixelLayerBuilder
ToolSvc += PixelLayerProvider
# put them to the caches
layerProviders += [PixelLayerProvider]
binningsCenter += [PixelLayerBinning]
binningsEndcap += [PixelLayerBinning]
colors += [3]
SCT_LayerBuilder = None
if buildCustomSCT:
# SCT building
from ISF_FatrasDetDescrTools.ISF_FatrasDetDescrToolsConf import iFatras__PlanarDetLayerBuilder
SCT_LayerBuilder = iFatras__PlanarDetLayerBuilder(
name=namePrefix + 'SCT_LayerBuilder')
SCT_LayerBuilder.PixelCase = False
SCT_LayerBuilder.Identification = 'SCT'
SCT_LayerBuilder.InputLayerMaterialProvider = ISF_InputLayerMaterialProvider
SCT_LayerBuilder.SiDetManagerLocation = 'SCT'
# SCT barrel specifications
SCT_LayerBuilder.BarrelLayerBinsZ = TrkDetFlags.SCT_BarrelLayerMaterialBinsZ(
)
SCT_LayerBuilder.BarrelLayerBinsPhi = TrkDetFlags.SCT_BarrelLayerMaterialBinsPhi(
)
# SCT endcap specifications
SCT_LayerBuilder.EndcapLayerBinsR = TrkDetFlags.SCT_EndcapLayerMaterialBinsR(
)
SCT_LayerBuilder.EndcapLayerBinsPhi = TrkDetFlags.SCT_EndcapLayerMaterialBinsPhi(
)
SCT_LayerBuilder.CustomMaterial = False
SCT_LayerBuilder.CustomMaterialThickness = 0.250
SCT_LayerBuilder.CustomMaterialX0 = 8.333
SCT_LayerBuilder.CustomMaterialL0 = 100
SCT_LayerBuilder.CustomMaterialA = 14
SCT_LayerBuilder.CustomMaterialZ = 28.0855
SCT_LayerBuilder.CustomMaterialRho = 0.00233
SCT_LayerBuilder.DetElementMapName = "SCT_IdHashDetElementMap"
SCT_LayerBuilder.BarrelLayers = 4
SCT_LayerBuilder.LayerSCTlike = True
#SCT_LayerBuilder.LayerInvertOtherSurface = True
SCT_LayerBuilder.LayersZsectors = [12, 12, 12, 12]
SCT_LayerBuilder.LayerPhiSectors = [32, 40, 48, 56]
SCT_LayerBuilder.LayerPhiOffset = [
(-360. / phi * 0.2 + 180. / phi)
for phi in SCT_LayerBuilder.LayerPhiSectors
] #degrees
SCT_LayerBuilder.LayerTilt = [11.0, 11.0, 11.25, 11.25] #degree
SCT_LayerBuilder.LayerMinPhi = [-180.0, -180.0, -180.0,
-180.0] #degree
SCT_LayerBuilder.LayerMaxPhi = [180.0, 180.0, 180.0,
180.0] #degree
SCT_LayerBuilder.LayerMinZ = [
-742.095, -742.095, -742.095, -742.095
]
SCT_LayerBuilder.LayerMaxZ = [742.095, 742.095, 742.095, 742.095]
SCT_LayerBuilder.LayerRadius = [299., 371., 443., 514.]
SCT_LayerBuilder.LayerRSeparation = [-3., -3., -3., -3.]
SCT_LayerBuilder.LayerThickness = [0.2850, 0.2850, 0.2850, 0.2850]
SCT_LayerBuilder.LayerLengthY = [126.09, 126.09, 126.09, 126.09]
SCT_LayerBuilder.LayerLengthXmin = [61.44, 61.44, 61.44, 61.44]
SCT_LayerBuilder.LayerPitchX = [0.080, 0.080, 0.080, 0.080]
#SCT_LayerBuilder.LayerPitchX = [ 6.144 , 6.144 , 6.144 , 6.144]
SCT_LayerBuilder.LayerPitchY = [126.09, 126.09, 126.09, 126.09]
SCT_LayerBuilder.LayerRotation = [0., 0., 0., 0.] #degree
SCT_LayerBuilder.LayerStereo = [-1.15, -1.15, -1.15, -1.15]
SCT_LayerBuilder.LayerStereoSeparation = [1., 1., 1., 1.]
SCT_LayerBuilder.AdditionalLayerRadius = []
SCT_LayerBuilder.EndcapDiscs = 9
SCT_LayerBuilder.DiscSCTlike = True
#SCT_LayerBuilder.DiscSCTlike = False
# discs are ordered from the smaller z to the bigger one
# from the smallest ring to the biggest
SCT_LayerBuilder.DiscPhiSectors = [[40, 52], [40, 40, 52],
[40, 40, 52], [40, 40, 52],
[40, 40, 52], [40, 40, 52],
[40, 52], [40, 52], [48]]
#SCT_LayerBuilder.DiscZpos = [-2735.45, -2520.25, -2130.45, -1786.65, -1415.95, -1315.15, -1106.75, -949.25, -869.05,
# 869.05, 949.25, 1106.75, 1315.15, 1415.95, 1786.65, 2130.45, 2520.25, 2735.45]
SCT_LayerBuilder.DiscZpos = [
-2742.45, -2511.25, -2123.45, -1778.55, -1406.95, -1307.15,
-1098.75, -941.25, -861.05, 860.05, 940.25, 1097.75, 1306.15,
1405.95, 1777.65, 2121.45, 2511.25, 2741.45
]
SCT_LayerBuilder.DiscRingMinR = [[337.6, 438.77],
[275.0, 337.6, 438.77],
[275.0, 337.6, 438.77],
[275.0, 337.6, 438.77],
[275.0, 337.6, 438.77],
[275.0, 337.6, 438.77],
[337.6, 438.77], [402.8, 438.77],
[438.77]]
SCT_LayerBuilder.DiscRingMaxR = [[455.3, 560.0],
[334.1, 455.3, 560.0],
[334.1, 455.3, 560.0],
[334.1, 455.3, 560.0],
[334.1, 455.3, 560.0],
[334.1, 455.3, 560.0],
[455.3, 560.0], [455.3, 560.0],
[560.0]]
SCT_LayerBuilder.DiscMinPhi = [[-180.0, -180.0],
[-180.0, -180.0, -180],
[-180.0, -180.0, -180],
[-180.0, -180.0, -180],
[-180.0, -180.0, -180],
[-180.0, -180.0, -180],
[-180.0, -180.0], [-180.0, -180.0],
[-180.0]]
SCT_LayerBuilder.DiscMaxPhi = [[180.0, 180.0], [180.0, 180.0, 180],
[180.0, 180.0, 180],
[180.0, 180.0, 180],
[180.0, 180.0, 180],
[180.0, 180.0, 180], [180.0, 180.0],
[180.0, 180.0], [180.0]]
SCT_LayerBuilder.DiscThickness = [
0.2850, 0.2850, 0.2850, 0.2850, 0.2850, 0.2850, 0.2850, 0.2850,
0.2850
]
#SCT_LayerBuilder.DiscLengthY = [[117.7, 117.7], [117.7, 117.7, 117.7], [117.7, 117.7, 117.7], [117.7, 117.7, 117.7], [117.7, 117.7, 117.7], [117.7, 117.7, 117.7], [117.7, 117.7], [117.7, 117.7], [117.7]]
#SCT_LayerBuilder.DiscLengthXmin = [[53.78, 53.78], [53.78, 53.78, 53.78], [53.78, 53.78, 53.78], [53.78, 53.78, 53.78], [53.78, 53.78, 53.78], [53.78, 53.78, 53.78], [53.78, 53.78], [53.78, 53.78], [53.78]]
#SCT_LayerBuilder.DiscLengthXmax = [[72.53, 72.53], [72.53, 72.53, 72.53], [72.53, 72.53, 72.53], [72.53, 72.53, 72.53], [72.53, 72.53, 72.53], [72.53, 72.53, 72.53], [72.53, 72.53], [72.53, 72.53], [72.53]]
#SCT_LayerBuilder.DiscLengthXmin = [[83.78, 83.78], [83.78, 83.78, 83.78], [83.78, 83.78, 83.78], [83.78, 83.78, 83.78], [83.78, 83.78, 83.78], [83.78, 83.78, 83.78], [83.78, 83.78], [83.78, 83.78], [83.78]]
#SCT_LayerBuilder.DiscLengthXmax = [[92.53, 92.53], [92.53, 92.53, 92.53], [92.53, 92.53, 92.53], [92.53, 92.53, 92.53], [92.53, 92.53, 92.53], [92.53, 92.53, 92.53], [92.53, 92.53], [92.53, 92.53], [92.53]]
SCT_LayerBuilder.DiscLengthXmin = [[53.784, 54.490],
[43.810, 53.784, 54.490],
[43.810, 53.784, 54.490],
[43.810, 53.784, 54.490],
[43.810, 53.784, 54.490],
[43.810, 53.784, 54.490],
[53.784, 54.490],
[64.174, 54.490], [54.490]]
SCT_LayerBuilder.DiscLengthXmax = [[72.534, 69.548],
[53.226, 72.534, 69.548],
[53.226, 72.534, 69.548],
[53.226, 72.534, 69.548],
[53.226, 72.534, 69.548],
[53.226, 72.534, 69.548],
[72.534, 69.548],
[72.534, 69.548], [69.548]]
SCT_LayerBuilder.DiscLengthY = [[117.700, 121.236],
[59.100, 117.700, 121.236],
[59.100, 117.700, 121.236],
[59.100, 117.7, 121.236],
[59.100, 117.7, 121.236],
[59.100, 117.7, 121.236],
[117.7, 121.236],
[52.474, 121.236], [121.236]]
#SCT_LayerBuilder.DiscPitchX = [[0.080, 0.080], [0.080, 0.080, 0.080], [0.080, 0.080, 0.080], [0.080, 0.080, 0.080], [0.080, 0.080, 0.080], [0.080, 0.080, 0.080], [0.080, 0.080], [0.080, 0.080], [0.080]]
#SCT_LayerBuilder.DiscPitchY = [[117.700, 121.236], [59.100, 117.700, 121.236], [59.100, 117.700, 121.236], [59.100, 117.7, 121.236], [59.100, 117.7, 121.236], [59.100, 117.7, 121.236], [117.7, 121.236], [52.474, 121.236], [121.236]]
SCT_LayerBuilder.DiscPitchX = [[207.0e-06, 161.5e-06],
[207.0e-06, 207.0e-06, 161.5e-06],
[207.0e-06, 207.0e-06, 161.5e-06],
[207.0e-06, 207.0e-06, 161.5e-06],
[207.0e-06, 207.0e-06, 161.5e-06],
[207.0e-06, 207.0e-06, 161.5e-06],
[207.0e-06, 161.5e-06],
[207.0e-06, 161.5e-06], [161.5e-06]
] #in rad for DiscSurfaces
SCT_LayerBuilder.DiscPitchY = [[117.700, 121.236],
[59.100, 117.700, 121.236],
[59.100, 117.700, 121.236],
[59.100, 117.7, 121.236],
[59.100, 117.7, 121.236],
[59.100, 117.7, 121.236],
[117.7, 121.236], [52.474, 121.236],
[121.236]]
#SCT_LayerBuilder.DiscPitchX = [[5.3784, 5.4490], [4.3810, 5.3784, 5.4490], [4.3810, 5.3784, 5.4490], [4.3810, 5.3784, 5.4490], [4.3810, 5.3784, 5.4490], [4.3810, 5.3784, 5.4490], [5.3784, 5.4490], [6.4174, 5.4490], [5.4490]]
#SCT_LayerBuilder.DiscPitchY = [[117.700, 121.236], [59.100, 117.700, 121.236], [59.100, 117.700, 121.236], [59.100, 117.7, 121.236], [59.100, 117.7, 121.236], [59.100, 117.7, 121.236], [117.7, 121.236], [52.474, 121.236], [121.236]]
SCT_LayerBuilder.DiscSeparation = [[5.3, 5.3], [5.3, 5.3, 5.3],
[5.1, 5.1,
5.1], [5.3, 5.3, 5.3],
[5.3, 5.3, 5.3],
[5.2, 5.2, 5.2], [5.2, 5.2],
[5.2, 5.2], [5.1]]
SCT_LayerBuilder.RingDisplacement = [[15.1, 15.1],
[15.1, 15.1, 15.1],
[15.2, 15.2, 15.2],
[15.1, 15.1, 15.1],
[15.1, 15.1, 15.1],
[15.15, 15.15, 15.15],
[15.15, 15.15],
[15.15, 15.15], [0.]]
SCT_LayerBuilder.DiscStereo = [[2.3, 2.3], [2.3, 2.3, 2.3],
[2.3, 2.3, 2.3], [2.3, 2.3, 2.3],
[2.3, 2.3, 2.3], [2.3, 2.3, 2.3],
[2.3, 2.3, 2.3], [2.3, 2.3],
[2.3, 2.3], [2.3]]
#SCT_LayerBuilder.DiscStereo = [[10.0, 10.0], [10., 10., 10.], [10., 10., 10.], [10., 10., 10.], [10., 10., 10.], [10., 10., 10.], [10., 10., 10.], [10., 10.], [10., 10.], [10.]]
SCT_LayerBuilder.DiscStereoSeparation = [[1., 1.], [1., 1., 1.],
[1., 1.,
1.], [1., 1., 1.],
[1., 1., 1.],
[1., 1., 1.], [1., 1.],
[1., 1.], [1.]]
#SCT_LayerBuilder.EndcapDiscs = 1
#SCT_LayerBuilder.DiscSCTlike = True
##SCT_LayerBuilder.DiscSCTlike = False
## discs are ordered from the smaller z to the bigger one
## from the smallest ring to the biggest
#SCT_LayerBuilder.DiscPhiSectors = [[16]]
#SCT_LayerBuilder.DiscZpos = [-1000, 1000]
#SCT_LayerBuilder.DiscRingMinR = [[337.6]]
#SCT_LayerBuilder.DiscRingMaxR = [[455.3]]
#SCT_LayerBuilder.DiscMinPhi = [[-180.0]]
#SCT_LayerBuilder.DiscMaxPhi = [[180.0]]
#SCT_LayerBuilder.DiscThickness = [0.2850]
#SCT_LayerBuilder.DiscLengthXmin = [[40]]
#SCT_LayerBuilder.DiscLengthXmax = [[80]]
#SCT_LayerBuilder.DiscLengthY = [[117.70]]
#SCT_LayerBuilder.DiscPitchX = [[0.080]]
#SCT_LayerBuilder.DiscPitchY = [[117.70]]
#SCT_LayerBuilder.DiscStereo = [[90.]]
#SCT_LayerBuilder.DiscStereoSeparation = [[1.]]
#SCT_LayerBuilder.EndcapDiscs = 1
#SCT_LayerBuilder.DiscSCTlike = True
##SCT_LayerBuilder.DiscSCTlike = False
## discs are ordered from the smaller z to the bigger one
## from the smallest ring to the biggest
#SCT_LayerBuilder.DiscPhiSectors = [[16]]
#SCT_LayerBuilder.DiscZpos = [-900., 900.]
#SCT_LayerBuilder.DiscRingMinR = [[337.6]]
#SCT_LayerBuilder.DiscRingMaxR = [[455.3]]
#SCT_LayerBuilder.DiscMinPhi = [[-180.0]]
#SCT_LayerBuilder.DiscMaxPhi = [[180.0]]
#SCT_LayerBuilder.DiscThickness = [0.2850]
#SCT_LayerBuilder.DiscLengthXmin = [[53.784]]
#SCT_LayerBuilder.DiscLengthXmax = [[72.534]]
#SCT_LayerBuilder.DiscLengthY = [[117.700]]
#SCT_LayerBuilder.DiscPitchX = [[207.0e-06]]
#SCT_LayerBuilder.DiscPitchX = [[0.0079739*2.]]
#SCT_LayerBuilder.DiscPitchY = [[117.700]]
#SCT_LayerBuilder.DiscStereo = [[10.]]
#SCT_LayerBuilder.DiscStereoSeparation = [[1.]]
SCT_LayerBuilder.AdditionalDiscZpos = [-2850, 2850]
# set the layer association
SCT_LayerBuilder.SetLayerAssociation = setLayerAssociation
# output level
SCT_LayerBuilder.OutputLevel = TrkDetFlags.SCT_BuildingOutputLevel(
)
# the binning type of the layer
SCT_LayerBinning = 2
# SCT -> ToolSvc
ToolSvc += SCT_LayerBuilder
print SCT_LayerBuilder
# the layer providers
SCT_LayerProvider = Trk__LayerProvider(name=namePrefix +
'SCT_LayerProvider')
SCT_LayerProvider.LayerBuilder = SCT_LayerBuilder
ToolSvc += SCT_LayerProvider
# put them to the caches
layerProviders += [SCT_LayerProvider]
binningsCenter += [SCT_LayerBinning]
binningsEndcap += [SCT_LayerBinning]
colors += [4]
TRT_LayerBuilder = None
if buildCustomTRT:
print "No TRT for the moment"
# TRT building
#from InDetTrackingGeometry.InDetTrackingGeometryConf import InDet__TRT_LayerBuilder
#TRT_LayerBuilder = InDet__TRT_LayerBuilder(name=namePrefix+'TRT_LayerBuilder')
## TRT barrel specifications
#TRT_LayerBuilder.ModelBarrelLayers = TrkDetFlags.TRT_BarrelModelLayers()
#TRT_LayerBuilder.BarrelLayerBinsZ = TrkDetFlags.TRT_BarrelLayerMaterialBinsZ()
#TRT_LayerBuilder.BarrelLayerBinsPhi = TrkDetFlags.TRT_BarrelLayerMaterialBinsPhi()
## SCT endcap specifications
#TRT_LayerBuilder.ModelEndcapLayers = TrkDetFlags.TRT_EndcapModelLayers()
#TRT_LayerBuilder.EndcapLayerBinsR = TrkDetFlags.TRT_EndcapLayerMaterialBinsR()
#TRT_LayerBuilder.EndcapLayerBinsPhi = TrkDetFlags.TRT_EndcapLayerMaterialBinsPhi()
## set the binning from bi-aequidistant to arbitrary for complex TRT volumes
#TRT_LayerBinning = 1
#if buildTrtStrawLayers or TrkDetFlags.TRT_BuildStrawLayers() :
#TRT_LayerBinning = 2
#TRT_LayerBuilder.ModelLayersOnly = False
## output level
#TRT_LayerBuilder.OutputLevel = TrkDetFlags.TRT_BuildingOutputLevel()
## TRT -> ToolSvc
#ToolSvc += TRT_LayerBuilder
## put them to the caches
#layerbuilders += [ TRT_LayerBuilder ]
#binnings += [ TRT_LayerBinning ]
#colors += [ 5 ]
# helpers for the InDetTrackingGeometry Builder : layer array creator
if createXML:
from ISF_FatrasDetDescrTools import XMLDumper
InDetXMLDumper = XMLDumper.XMLDumper(name='InDetXMLDumper',
fileName=dictionaryFileName)
if IBLLayerBuilder:
InDetXMLDumper.addPixelBuilder(IBLLayerBuilder)
if PixelLayerBuilder:
InDetXMLDumper.addPixelBuilder(PixelLayerBuilder)
if SCT_LayerBuilder:
InDetXMLDumper.addSCTBuilder(SCT_LayerBuilder)
if TRT_LayerBuilder:
InDetXMLDumper.addTRTBuilder(TRT_LayerBuilder)
InDetXMLDumper.dump()
InDetXMLDumper.closeFile()
# layer && tracking volume array creator
from TrkDetDescrTools.TrkDetDescrToolsConf import Trk__LayerArrayCreator
InDetLayerArrayCreator = Trk__LayerArrayCreator(
name='InDetLayerArrayCreator')
InDetLayerArrayCreator.EmptyLayerMode = 2 # deletes empty material layers from arrays
InDetLayerArrayCreator.OutputLevel = TrkDetFlags.InDetBuildingOutputLevel(
)
# add to ToolSvc
ToolSvc += InDetLayerArrayCreator
# helpers for the InDetTrackingGeometry Builder : volume array creator
from TrkDetDescrTools.TrkDetDescrToolsConf import Trk__TrackingVolumeArrayCreator
InDetTrackingVolumeArrayCreator = Trk__TrackingVolumeArrayCreator(
name='InDetTrackingVolumeArrayCreator')
InDetTrackingVolumeArrayCreator.OutputLevel = TrkDetFlags.InDetBuildingOutputLevel(
)
# add to ToolSvc
ToolSvc += InDetTrackingVolumeArrayCreator
# helpers for the InDetTrackingGeometry Builder : tracking voluem helper for glueing
from TrkDetDescrTools.TrkDetDescrToolsConf import Trk__TrackingVolumeHelper
InDetTrackingVolumeHelper = Trk__TrackingVolumeHelper(
name='InDetTrackingVolumeHelper')
InDetTrackingVolumeHelper.OutputLevel = TrkDetFlags.InDetBuildingOutputLevel(
)
# the material bins
InDetTrackingVolumeHelper.BarrelLayerBinsZ = TrkDetFlags.InDetPassiveLayerMaterialBinsRz(
)
InDetTrackingVolumeHelper.BarrelLayerBinsPhi = TrkDetFlags.InDetPassiveLayerMaterialBinsPhi(
)
InDetTrackingVolumeHelper.EndcapLayerBinsR = TrkDetFlags.InDetPassiveLayerMaterialBinsRz(
)
InDetTrackingVolumeHelper.EndcapLayerBinsPhi = TrkDetFlags.InDetPassiveLayerMaterialBinsPhi(
)
# add to ToolSvc
ToolSvc += InDetTrackingVolumeHelper
# helpers for the InDetTrackingGeometry Builder : cylinder volume creator
from TrkDetDescrTools.TrkDetDescrToolsConf import Trk__CylinderVolumeCreator
InDetCylinderVolumeCreator = Trk__CylinderVolumeCreator(
name='InDetCylinderVolumeCreator')
# give him the layer array creator
InDetCylinderVolumeCreator.LayerArrayCreator = InDetLayerArrayCreator
InDetCylinderVolumeCreator.TrackingVolumeArrayCreator = InDetTrackingVolumeArrayCreator
InDetCylinderVolumeCreator.TrackingVolumeHelper = InDetTrackingVolumeHelper
# specifiy the binning, passive layers, entry layers
InDetCylinderVolumeCreator.PassiveLayerThickness = TrkDetFlags.InDetPassiveLayerThickness(
)
InDetCylinderVolumeCreator.PassiveLayerBinsPhi = TrkDetFlags.InDetPassiveLayerMaterialBinsPhi(
)
InDetCylinderVolumeCreator.PassiveLayerBinsRZ = TrkDetFlags.InDetPassiveLayerMaterialBinsRz(
)
# output level
InDetCylinderVolumeCreator.OutputLevel = TrkDetFlags.InDetBuildingOutputLevel(
)
# add to ToolSvc
ToolSvc += InDetCylinderVolumeCreator
# the envelope definition service
from SubDetectorEnvelopes.SubDetectorEnvelopesConf import DetDescrDBEnvelopeSvc
AtlasEnvelopeSvc = DetDescrDBEnvelopeSvc(
'AtlasEnvelopeSvcDefinitionSvc')
# set the output level for the Envelope service
AtlasEnvelopeSvc.OutputLevel = TrkDetFlags.InDetBuildingOutputLevel()
# add to SvcMgr
ServiceMgr += AtlasEnvelopeSvc
# the tracking geometry builder
InDet__StagedTrackingGeometryBuilder.__init__(
self,
namePrefix + name,
LayerBuilders=layerProviders,
LayerBinningTypeCenter=binningsCenter,
LayerBinningTypeEndcap=binningsEndcap,
ColorCodes=colors,
EnvelopeDefinitionSvc=AtlasEnvelopeSvc,
VolumeEnclosureDiscPositions=[3000., 3450.],
TrackingVolumeCreator=InDetCylinderVolumeCreator,
LayerArrayCreator=InDetLayerArrayCreator,
CheckForRingLayout=False,
BuildBoundaryLayers=TrkDetFlags.InDetBuildMaterialBoundaries(),
ReplaceAllJointBoundaries=TrkDetFlags.InDetBuildJointBoundaries(),
OutputLevel=TrkDetFlags.InDetBuildingOutputLevel(),
ExitVolumeName=TrkDetFlags.InDetContainerName(),
MagneticFieldMode=TrkDetFlags.MagneticFieldMode())
def __init__(self,name = 'InDetTrackingGeometryBuilder',
namePrefix = '',
setLayerAssociation = True):
# get the ToolSvc
from AthenaCommon.AppMgr import ToolSvc, ServiceMgr
from AthenaCommon.DetFlags import DetFlags
# the Detector flags to be imported
from TrkDetDescrSvc.TrkDetDescrJobProperties import TrkDetFlags
# beampipe
from InDetTrackingGeometry.InDetTrackingGeometryConf import InDet__BeamPipeBuilder
BeamPipeBuilder = InDet__BeamPipeBuilder(name=namePrefix+'BeamPipeBuilder')
BeamPipeBuilder.BeamPipeRadius = TrkDetFlags.BeamPipeRadius()
BeamPipeBuilder.BeamPipeThickness = TrkDetFlags.BeamPipeThickness()
BeamPipeBuilder.BeamPipeMaterialBinsZ = TrkDetFlags.BeamPipeLayerMaterialBinsZ()
BeamPipeBuilder.OutputLevel = TrkDetFlags.InDetBuildingOutputLevel()
# the binning type of the layers
BeamPipeBinning = 2
# where to build it from
BeamPipeBuilder.BeamPipeFromGeoModel = TrkDetFlags.BeamPipeFromGeoModel()
BeamPipeBuilder.BeamPipeEnvelope = TrkDetFlags.BeamPipeEnvelope()
if not TrkDetFlags.BeamPipeFromGeoModel():
BeamPipeBuilder.BeamPipeOffsetX = TrkDetFlags.BeamPipeOffsetX()
BeamPipeBuilder.BeamPipeOffsetY = TrkDetFlags.BeamPipeOffsetY()
BeamPipeBuilder.BeamPipeRadius = TrkDetFlags.BeamPipeRadius()
ToolSvc += BeamPipeBuilder
# the layer providers
from TrkDetDescrTools.TrkDetDescrToolsConf import Trk__LayerProvider
BeamPipeProvider = Trk__LayerProvider(name=namePrefix+'BeamPipeProvider')
BeamPipeProvider.LayerBuilder = BeamPipeBuilder
ToolSvc += BeamPipeProvider
# the layer builders
layerProviders = [ BeamPipeProvider ]
binningsCenter = [ BeamPipeBinning ]
binningsEndcap = [ BeamPipeBinning ]
colors = [ 2 ]
# PIXEL building
from InDetTrackingGeometry.InDetTrackingGeometryConf import InDet__SiLayerBuilder
PixelLayerBuilder = InDet__SiLayerBuilder(name=namePrefix+'PixelLayerBuilder')
PixelLayerBuilder.PixelCase = True
PixelLayerBuilder.Identification = 'Pixel1'
PixelLayerBuilder.SiDetManagerLocation = 'Pixel'
# The pixel split mode
PixelLayerBuilder.SplitMode = TrkDetFlags.PixelSplitMode()
# Pixel barrel specifications
PixelLayerBuilder.BarrelLayerBinsZ = TrkDetFlags.PixelBarrelLayerMaterialBinsZ()
PixelLayerBuilder.BarrelLayerBinsPhi = TrkDetFlags.PixelBarrelLayerMaterialBinsPhi()
PixelLayerBuilder.EndcapLayerBinsR = TrkDetFlags.PixelEndcapLayerMaterialBinsR()
PixelLayerBuilder.EndcapLayerBinsPhi = TrkDetFlags.PixelEndcapLayerMaterialBinsPhi()
# set the layer association
PixelLayerBuilder.SetLayerAssociation = setLayerAssociation
# output level
PixelLayerBuilder.OutputLevel = TrkDetFlags.PixelBuildingOutputLevel()
# the binning type of the layers
PixelLayerBinning = 2
# add it to tool service
ToolSvc += PixelLayerBuilder
# put them to the caches
PixelLayerProvider = Trk__LayerProvider(name=namePrefix+'PixelLayerProvider')
PixelLayerProvider.LayerBuilder = PixelLayerBuilder
ToolSvc += PixelLayerProvider
layerProviders += [ PixelLayerProvider ]
binningsCenter += [ PixelLayerBinning ]
binningsEndcap += [ PixelLayerBinning ]
colors += [ 3 ]
# check if the pixel split mode is one
if TrkDetFlags.PixelSplitMode() != 0:
# set split mode to pixel layer builder
print ("[SLHC_GeometryBuilder] Building Pixel layers in split mode")
# change the identificaiton of the first one
PixelLayerBuilder.Identification = 'Pixel1'
# and now configure the second one
PixelLayerBuilder2 = InDet__SiLayerBuilder(name=namePrefix+'PixelLayerBuilder2')
PixelLayerBuilder2.PixelCase = True
PixelLayerBuilder2.SplitMode = TrkDetFlags.PixelSplitMode()
PixelLayerBuilder2.Identification = 'Pixel2'
PixelLayerBuilder2.SiDetManagerLocation = 'Pixel'
# Pixel barrel specifications
PixelLayerBuilder2.BarrelLayerBinsZ = TrkDetFlags.PixelBarrelLayerMaterialBinsZ()
PixelLayerBuilder2.BarrelLayerBinsPhi = TrkDetFlags.PixelBarrelLayerMaterialBinsPhi()
PixelLayerBuilder2.EndcapLayerBinsR = TrkDetFlags.PixelEndcapLayerMaterialBinsR()
PixelLayerBuilder2.EndcapLayerBinsPhi = TrkDetFlags.PixelEndcapLayerMaterialBinsPhi()
# set the layer association
PixelLayerBuilder2.SetLayerAssociation = setLayerAssociation
# output level
PixelLayerBuilder2.OutputLevel = TrkDetFlags.PixelBuildingOutputLevel()
# add it to tool service
ToolSvc += PixelLayerBuilder2
# the layer providers
PixelLayerProvider2 = Trk__LayerProvider(name=namePrefix+'PixelLayerProvider2')
PixelLayerProvider2.LayerBuilder = PixelLayerBuilder2
ToolSvc += PixelLayerProvider2
# put them to the caches
layerProviders += [ PixelLayerProvider2 ]
binningsCenter += [ PixelLayerBinning ]
binningsEndcap += [ PixelLayerBinning ]
colors += [ 3 ]
if DetFlags.SCT_on() :
# SCT building
SCT_LayerBuilder = InDet__SiLayerBuilder(name=namePrefix+'SCT_LayerBuilder')
SCT_LayerBuilder.PixelCase = False
SCT_LayerBuilder.Identification = 'SCT1'
SCT_LayerBuilder.SiDetManagerLocation = 'SCT'
# general steering
SCT_LayerBuilder.SplitMode = TrkDetFlags.SCT_SplitMode()
# SCT barrel specifications
SCT_LayerBuilder.BarrelLayerBinsZ = TrkDetFlags.SCT_BarrelLayerMaterialBinsZ()
SCT_LayerBuilder.BarrelLayerBinsPhi = TrkDetFlags.SCT_BarrelLayerMaterialBinsPhi()
# SCT endcap specifications
SCT_LayerBuilder.EndcapLayerBinsR = TrkDetFlags.SCT_EndcapLayerMaterialBinsR()
SCT_LayerBuilder.EndcapLayerBinsPhi = TrkDetFlags.SCT_EndcapLayerMaterialBinsPhi()
# SCT_LayerBuilder.EndcapComplexRingBinning = TrkDetFlags.SCT_EndcapLayerDynamicRings()
# set the layer association
SCT_LayerBuilder.SetLayerAssociation = setLayerAssociation
# output level
SCT_LayerBuilder.OutputLevel = TrkDetFlags.SCT_BuildingOutputLevel()
# the binning type of the layer
SCT_LayerBinning = 2
# SCT -> ToolSvc
ToolSvc += SCT_LayerBuilder
# the layer providers
SCT_LayerProvider = Trk__LayerProvider(name=namePrefix+'SCT_LayerProvider')
SCT_LayerProvider.LayerBuilder = SCT_LayerBuilder
ToolSvc += SCT_LayerProvider
# put them to the caches
layerProviders += [ SCT_LayerProvider ]
binningsCenter += [ SCT_LayerBinning ]
binningsEndcap += [ SCT_LayerBinning ]
colors += [ 4 ]
# the split and the second SCT layer builder is needed if the SCT outer barrel
# spans over the full width
if TrkDetFlags.SCT_SplitMode() != 0 :
print ("[SLHC_GeometryBuilder] Building SCT layers in split mode")
from InDetTrackingGeometry.InDetTrackingGeometryConf import InDet__SiLayerBuilder
SCT_LayerBuilder2 = InDet__SiLayerBuilder(name=namePrefix+'SCT_LayerBuilder2')
# switch of the pixel case
SCT_LayerBuilder2.PixelCase = False
SCT_LayerBuilder2.SiDetManagerLocation = 'SCT'
# set the general parameters
SCT_LayerBuilder2.OutputLevel = 2
SCT_LayerBuilder2.SplitMode = TrkDetFlags.SCT_SplitMode()
SCT_LayerBuilder2.Identification = "SCT2"
# SCT barrel specifications
SCT_LayerBuilder2.BarrelLayerBinsZ = TrkDetFlags.SCT_BarrelLayerMaterialBinsZ()
SCT_LayerBuilder2.BarrelLayerBinsPhi = TrkDetFlags.SCT_BarrelLayerMaterialBinsPhi()
# SCT endcap specifications
SCT_LayerBuilder2.EndcapLayerBinsR = TrkDetFlags.SCT_EndcapLayerMaterialBinsR()
SCT_LayerBuilder2.EndcapLayerBinsPhi = TrkDetFlags.SCT_EndcapLayerMaterialBinsPhi()
# set the layer association
SCT_LayerBuilder2.SetLayerAssociation = setLayerAssociation
# output level
SCT_LayerBuilder2.OutputLevel = TrkDetFlags.SCT_BuildingOutputLevel()
ToolSvc += SCT_LayerBuilder2
# the layer providers
SCT_LayerProvider2 = Trk__LayerProvider(name=namePrefix+'SCT_LayerProvider2')
SCT_LayerProvider2.LayerBuilder = SCT_LayerBuilder2
ToolSvc += SCT_LayerProvider2
# put them to the caches
layerProviders += [ SCT_LayerProvider2 ]
binningsCenter += [ SCT_LayerBinning ]
binningsEndcap += [ SCT_LayerBinning ]
colors += [ 5 ]
# helpers for the InDetTrackingGeometry Builder : layer array creator
from TrkDetDescrTools.TrkDetDescrToolsConf import Trk__LayerArrayCreator
InDetLayerArrayCreator = Trk__LayerArrayCreator(name = 'InDetLayerArrayCreator')
InDetLayerArrayCreator.EmptyLayerMode = 2 # deletes empty material layers from arrays
InDetLayerArrayCreator.OutputLevel = TrkDetFlags.InDetBuildingOutputLevel()
# add to ToolSvc
ToolSvc += InDetLayerArrayCreator
# helpers for the InDetTrackingGeometry Builder : volume array creator
from TrkDetDescrTools.TrkDetDescrToolsConf import Trk__TrackingVolumeArrayCreator
InDetTrackingVolumeArrayCreator = Trk__TrackingVolumeArrayCreator(name = 'InDetTrackingVolumeArrayCreator')
InDetTrackingVolumeArrayCreator.OutputLevel = TrkDetFlags.InDetBuildingOutputLevel()
# add to ToolSvc
ToolSvc += InDetTrackingVolumeArrayCreator
# helpers for the InDetTrackingGeometry Builder : tracking voluem helper for glueing
from TrkDetDescrTools.TrkDetDescrToolsConf import Trk__TrackingVolumeHelper
InDetTrackingVolumeHelper = Trk__TrackingVolumeHelper(name ='InDetTrackingVolumeHelper')
InDetTrackingVolumeHelper.OutputLevel = TrkDetFlags.InDetBuildingOutputLevel()
# the material bins
InDetTrackingVolumeHelper.BarrelLayerBinsZ = TrkDetFlags.InDetPassiveLayerMaterialBinsRz()
InDetTrackingVolumeHelper.BarrelLayerBinsPhi = TrkDetFlags.InDetPassiveLayerMaterialBinsPhi()
InDetTrackingVolumeHelper.EndcapLayerBinsR = TrkDetFlags.InDetPassiveLayerMaterialBinsRz()
InDetTrackingVolumeHelper.EndcapLayerBinsPhi = TrkDetFlags.InDetPassiveLayerMaterialBinsPhi()
# add to ToolSvc
ToolSvc += InDetTrackingVolumeHelper
# helpers for the InDetTrackingGeometry Builder : cylinder volume creator
from TrkDetDescrTools.TrkDetDescrToolsConf import Trk__CylinderVolumeCreator
InDetCylinderVolumeCreator = Trk__CylinderVolumeCreator(name = 'InDetCylinderVolumeCreator')
# give him the layer array creator
InDetCylinderVolumeCreator.LayerArrayCreator = InDetLayerArrayCreator
InDetCylinderVolumeCreator.TrackingVolumeArrayCreator = InDetTrackingVolumeArrayCreator
InDetCylinderVolumeCreator.TrackingVolumeHelper = InDetTrackingVolumeHelper
# specifiy the binning, passive layers, entry layers
InDetCylinderVolumeCreator.PassiveLayerThickness = TrkDetFlags.InDetPassiveLayerThickness()
InDetCylinderVolumeCreator.PassiveLayerBinsPhi = TrkDetFlags.InDetPassiveLayerMaterialBinsPhi()
InDetCylinderVolumeCreator.PassiveLayerBinsRZ = TrkDetFlags.InDetPassiveLayerMaterialBinsRz()
# output level
InDetCylinderVolumeCreator.OutputLevel = TrkDetFlags.InDetBuildingOutputLevel()
# add to ToolSvc
ToolSvc += InDetCylinderVolumeCreator
# the envelope definition service
from SubDetectorEnvelopes.SubDetectorEnvelopesConf import DetDescrDBEnvelopeSvc
AtlasEnvelopeSvc = DetDescrDBEnvelopeSvc('AtlasEnvelopeSvcDefinitionSvc')
# set the output level for the Envelope service
AtlasEnvelopeSvc.OutputLevel = TrkDetFlags.InDetBuildingOutputLevel()
# add to SvcMgr
ServiceMgr += AtlasEnvelopeSvc
# the tracking geometry builder
InDet__StagedTrackingGeometryBuilder.__init__(self,namePrefix+name,\
BeamPipeBuilder = BeamPipeBuilder,\
LayerBuilders = layerProviders,
LayerBinningTypeCenter = binningsCenter,
LayerBinningTypeEndcap = binningsEndcap,
ColorCodes = colors,
EnvelopeDefinitionSvc = AtlasEnvelopeSvc,
# VolumeEnclosureDiscPositions = [ 3000., 3450. ],
TrackingVolumeCreator = InDetCylinderVolumeCreator,
LayerArrayCreator = InDetLayerArrayCreator,
BuildBoundaryLayers = True,
# ReplaceAllJointBoundaries = TrkDetFlags.InDetBuildMaterialBoundaries(),
OutputLevel = TrkDetFlags.InDetBuildingOutputLevel(),
ExitVolumeName = TrkDetFlags.InDetContainerName(),
MagneticFieldMode = TrkDetFlags.MagneticFieldMode(),
isSLHC = True)
