OutputLevel = INFO)
HCalExtBcells = CellPositionsHCalBarrelNoSegTool("CellPositionsHCalExtBarrel",
readoutName = hcalExtBarrelReadoutName,
OutputLevel = INFO)
HECcells = CellPositionsCaloDiscsTool("CellPositionsHEC",
readoutName = hcalEndcapReadoutName,
OutputLevel = INFO)
HCalFwdcells = CellPositionsCaloDiscsTool("CellPositionsHCalFwd",
readoutName = hcalFwdReadoutName,
OutputLevel = INFO)
# cell positions
from Configurables import CreateCellPositions
positionsEcalBarrel = CreateCellPositions("positionsEcalBarrel",
positionsTool=ECalBcells,
hits = "ECalBarrelCells",
positionedHits = "ECalBarrelCellPositions",
OutputLevel = INFO)
positionsHcalBarrel = CreateCellPositions("positionsHcalBarrel",
positionsTool=HCalBcells,
hits = "HCalBarrelCells",
positionedHits = "HCalBarrelCellPositions",
OutputLevel = INFO)
positionsHcalExtBarrel = CreateCellPositions("positionsHcalExtBarrel",
positionsTool=HCalExtBcells,
hits = "HCalExtBarrelCells",
positionedHits = "HCalExtBarrelCellPositions",
OutputLevel = INFO)
positionsEcalEndcap = CreateCellPositions("positionsEcalEndcap",
positionsTool=EMECcells,
hits = "ECalEndcapCells",
HCalBcellVols = CellPositionsHCalBarrelNoSegTool("CellPositionsHCalBarrelVols",
readoutName = hcalBarrelReadoutName,
OutputLevel = INFO)
HCalBsegcells = CellPositionsHCalBarrelTool("CellPositionsHCalSegBarrel",
readoutName = hcalBarrelReadoutNamePhiEta,
radii = [2910.5, 3010.5, 3135.5, 3285.5, 3435.5, 3585.5, 3785.5, 4135.5, 4285.5, 4535.5],
OutputLevel = INFO)
############################################################################################################## ####### RESEGMENT HCAL ############# ##############################################################################################################
from Configurables import CreateVolumeCaloPositions,RedoSegmentation,CreateCaloCells,CreateCellPositions
# Create cells in HCal
# 2. step - rewrite the cellId using the Phi-Eta segmentation
# 3. step - merge new cells corresponding to eta-phi segmentation
# Hcal barrel cell positions
posHcalBarrel = CreateCellPositions("posHcalBarrel",
positionsTool=HCalBcellVols,
hits = prefix+"HCalBarrelCells",
positionedHits = "HCalBarrelPositions")
# Use Phi-Eta segmentation in Hcal barrel
resegmentHcalBarrel = RedoSegmentation("ReSegmentationHcal",
# old bitfield (readout)
oldReadoutName = hcalBarrelReadoutName,
# specify which fields are going to be altered (deleted/rewritten)
oldSegmentationIds = ["module","row"],
# new bitfield (readout), with new segmentation
newReadoutName = hcalBarrelReadoutNamePhiEta,
inhits = "HCalBarrelPositions",
outhits = "HCalBarrelCellsStep2")
createHcalBarrelCells = CreateCaloCells("CreateHCalBarrelCells",
doCellCalibration=False,recalibrateBaseline =False,
createHcalBarrelCellsNoise.cells.Path = "HCalBarrelCellsNoise"
else:
createHcalBarrelCellsNoise = CreateCaloCells("CreateHCalBarrelCellsNoise",
geometryTool = hcalgeo,
doCellCalibration = False, recalibrateBaseline =False,
addCellNoise = True, filterCellNoise = False,
noiseTool = noiseHcal)
createHcalBarrelCellsNoise.hits.Path = hcalCells
createHcalBarrelCellsNoise.cells.Path = "HCalBarrelCellsNoise"
# cell positions
from Configurables import CreateCellPositions
positionsEcalBarrel = CreateCellPositions("positionsEcalBarrel",
positionsTool=ECalBcells,
hits = ecalBarrelCellsForPositions,
positionedHits = "ECalBarrelCellPositions",
OutputLevel = INFO)
positionsHcalBarrel = CreateCellPositions("positionsHcalBarrel",
positionsTool=HCalBcells,
hits = hcalBarrelCellsForPositions,
positionedHits = "HCalBarrelCellPositions",
OutputLevel = INFO)
positionsHcalExtBarrel = CreateCellPositions("positionsHcalExtBarrel",
positionsTool=HCalExtBcells,
hits = prefix+"HCalExtBarrelCells",
positionedHits = "HCalExtBarrelCellPositions",
OutputLevel = INFO)
positionsHcalSegBarrel = CreateCellPositions("positionsSegHcalBarrel",
positionsTool=HCalBsegcells,
hits = hcalBarrelCellsForPositions,