# active material volume name
hcalVolumeName = ["moduleVolume", "wedgeVolume", "layerVolume"]
# ECAL bitfield names & values
hcalFieldNames=["system"]
hcalFieldValues=[8]
############################################################################################################## ####### CELL POSITIONS TOOLS ############# ############################################################################################################## # Configure tools for calo cell positions
from Configurables import CellPositionsECalBarrelTool, CellPositionsHCalBarrelNoSegTool, CellPositionsHCalBarrelTool
ECalBcells = CellPositionsECalBarrelTool("CellPositionsECalBarrel",
readoutName = ecalBarrelReadoutNamePhiEta,
OutputLevel = INFO)
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",
smear = SimG4SmearGenParticles("smear",
inParticles = "GenParticles",
smearedParticles = "SmearedParticles",
smearTool = smeartool,
OutputLevel=DEBUG)
#Configure tools for calo reconstruction
from Configurables import ConstNoiseTool
noiseTool = ConstNoiseTool("ConstNoiseTool")
#Configure tools for calo cell positions
from Configurables import CellPositionsECalBarrelTool, CellPositionsHCalBarrelTool, CellPositionsHCalBarrelNoSegTool, CellPositionsCaloDiscsTool, CellPositionsTailCatcherTool
ECalBcells = CellPositionsECalBarrelTool("CellPositionsECalBarrel",
readoutName = ecalBarrelReadoutName)
HCalBcells = CellPositionsHCalBarrelTool("CellPositionsHCalBarrel",
readoutName = hcalBarrelReadoutName, #"BarHCal_Readout_phieta",
radii = [291.05, 301.05, 313.55, 328.55, 343.55, 358.55, 378.55, 413.55, 428.55, 453.55],
OutputLevel = INFO)
HCalExtBcells = CellPositionsHCalBarrelNoSegTool("CellPositionsHCalExtBarrel",
readoutName = hcalExtBarrelReadoutName)
EMECcells = CellPositionsCaloDiscsTool("CellPositionsEMEC",
readoutName = ecalEndcapReadoutName,
OutputLevel = INFO)
ECalFwdcells = CellPositionsCaloDiscsTool("CellPositionsECalFwd",
readoutName = ecalFwdReadoutName,
OutputLevel = INFO)
HECcells = CellPositionsCaloDiscsTool("CellPositionsHEC",
readoutName = hcalEndcapReadoutName,
OutputLevel = INFO)
HCalFwdcells = CellPositionsCaloDiscsTool("CellPositionsHCalFwd",
readoutName = hcalFwdReadoutName,
OutputLevel = INFO)
from Configurables import CellPositionsCaloDiscsTool
EMECcells = CellPositionsCaloDiscsTool("CellPositionsEMEC")
EMECcells.readoutName = ecalBarrelReadoutName
from Configurables import CellPositionsCaloDiscsTool
ECalFwdcells = CellPositionsCaloDiscsTool("CellPositionsECalFwd")
ECalFwdcells.readoutName = ecalFwdReadoutName
ECalFwdcells.OutputLevel = INFO
from Configurables import CellPositionsHCalBarrelNoSegTool
HCalBcellVols = CellPositionsHCalBarrelNoSegTool("CellPositionsHCalBarrelVols")
HCalBcellVols.readoutName = "HCalBarrelReadout"
HCalBcellVols.OutputLevel = INFO
from Configurables import CellPositionsHCalBarrelTool
HCalBsegcells = CellPositionsHCalBarrelTool("CellPositionsHCalBarrel")
HCalBsegcells.readoutName = "BarHCal_Readout_phieta"
HCalBsegcells.radii = [291.05, 301.05, 313.55, 328.55, 343.55, 358.55, 378.55, 413.55, 428.55, 453.55]
HCalBsegcells.OutputLevel = INFO
from Configurables import CellPositionsCaloDiscsTool
HECcells = CellPositionsCaloDiscsTool("CellPositionsHEC")
HECcells.readoutName = hcalEndcapReadoutName
from Configurables import CellPositionsCaloDiscsTool
HCalFwdcells = CellPositionsCaloDiscsTool("CellPositionsHCalFwd")
HCalFwdcells.readoutName = hcalFwdReadoutName
# additionally for HCal
from Configurables import RewriteBitfield
#Configure tools for calo cell positions
from Configurables import CellPositionsECalBarrelTool, CellPositionsHCalBarrelTool, CellPositionsCaloDiscsTool, CellPositionsTailCatcherTool
ECalBcells = CellPositionsECalBarrelTool("CellPositionsECalBarrel",
readoutName=ecalBarrelReadoutName,
OutputLevel=INFO)
EMECcells = CellPositionsCaloDiscsTool("CellPositionsEMEC",
readoutName=ecalEndcapReadoutName,
OutputLevel=INFO)
ECalFwdcells = CellPositionsCaloDiscsTool("CellPositionsECalFwd",
readoutName=ecalFwdReadoutName,
OutputLevel=INFO)
HCalBcells = CellPositionsHCalBarrelTool("CellPositionsHCalBarrel",
readoutName=hcalBarrelReadoutName,
radii=[
291.05, 301.05, 313.55, 328.55,
343.55, 358.55, 378.55, 403.55,
428.55, 453.55
],
OutputLevel=INFO)
HCalExtBcells = CellPositionsHCalBarrelTool(
"CellPositionsHCalExtBarrel",
readoutName=hcalExtBarrelReadoutName,
radii=[
356.05, 373.55, 398.55, 423.55, 291.05, 301.05, 313.55, 328.55, 348.55,
373.55, 398.55, 423.55
],
OutputLevel=INFO)
HECcells = CellPositionsCaloDiscsTool("CellPositionsHEC",
readoutName=hcalEndcapReadoutName,
OutputLevel=INFO)
HCalFwdcells = CellPositionsCaloDiscsTool("CellPositionsHCalFwd",
##############################################################################################################
#Configure tools for calo cell positions
from Configurables import CellPositionsECalBarrelTool, CellPositionsHCalBarrelTool, CellPositionsCaloDiscsTool, CellPositionsCaloDiscsTool, CellPositionsTailCatcherTool
ECalBcells = CellPositionsECalBarrelTool(
"CellPositionsECalBarrel",
readoutName=ecalBarrelReadoutNamePhiEta,
OutputLevel=INFO)
EMECcells = CellPositionsCaloDiscsTool("CellPositionsEMEC",
readoutName=ecalEndcapReadoutName,
OutputLevel=INFO)
ECalFwdcells = CellPositionsCaloDiscsTool("CellPositionsECalFwd",
readoutName=ecalFwdReadoutName,
OutputLevel=INFO)
HCalBcells = CellPositionsHCalBarrelTool("CellPositionsHCalBarrel",
readoutName=hcalBarrelReadoutName,
OutputLevel=INFO)
HCalExtBcells = CellPositionsHCalBarrelTool(
"CellPositionsHCalExtBarrel",
readoutName=hcalExtBarrelReadoutName,
OutputLevel=INFO)
HECcells = CellPositionsCaloDiscsTool("CellPositionsHEC",
readoutName=hcalEndcapReadoutName,
OutputLevel=INFO)
HCalFwdcells = CellPositionsCaloDiscsTool("CellPositionsHCalFwd",
readoutName=hcalFwdReadoutName,
OutputLevel=INFO)
TailCatchercells = CellPositionsTailCatcherTool(
"CellPositionsTailCatcher",
readoutName=tailCatcherReadoutName,
centralRadius=901.5,
readoutName = ecalBarrelReadoutName,
OutputLevel = INFO)
EMECcells = CellPositionsCaloDiscsTool("CellPositionsEMEC",
readoutName = ecalEndcapReadoutName,
OutputLevel = INFO)
ECalFwdcells = CellPositionsCaloDiscsTool("CellPositionsECalFwd",
readoutName = ecalFwdReadoutName,
OutputLevel = INFO)
HCalBcells = CellPositionsHCalBarrelNoSegTool("CellPositionsHCalBarrel",
readoutName = hcalBarrelReadoutName,
OutputLevel = INFO)
HCalExtBcells = CellPositionsHCalBarrelNoSegTool("CellPositionsHCalExtBarrel",
readoutName = hcalExtBarrelReadoutName,
OutputLevel = INFO)
HCalBsegcells = CellPositionsHCalBarrelTool("CellPositionsHCalSegBarrel",
readoutName = hcalBarrelReadoutNamePhiEta,
radii = [2910.5, 3010.5, 3135.5, 3285.5, 3435.5, 3585.5, 3785.5, 4035.5, 4285.5, 4535.5],
OutputLevel = INFO)
HCalExtBsegcells = CellPositionsHCalBarrelTool("CellPositionsHCalExtSegBarrel",
readoutName = hcalExtBarrelReadoutNamePhiEta,
radii = [ 356.05
, 373.55
, 398.55
, 423.55
, 291.05
, 301.05
, 313.55
, 328.55
, 348.55
, 373.55
, 398.55
, 423.55