createcells.hits.Path = "mergedECalHits"
createcells.cells.Path = "ecalBarrelCells"
from Configurables import CreateEmptyCaloCellsCollection
createemptycells = CreateEmptyCaloCellsCollection("CreateEmptyCaloCells")
createemptycells.cells.Path = "emptyCaloCells"
#Create calo clusters
from Configurables import CreateCaloClustersSlidingWindow, CaloTowerTool
from GaudiKernel.PhysicalConstants import pi
towers = CaloTowerTool("towers",
deltaEtaTower=0.01,
deltaPhiTower=2 * pi / 629.,
ecalBarrelReadoutName=ecalReadoutName,
ecalEndcapReadoutName="",
ecalFwdReadoutName="",
hcalBarrelReadoutName="",
hcalExtBarrelReadoutName="",
hcalEndcapReadoutName="",
hcalFwdReadoutName="",
OutputLevel=DEBUG)
towers.ecalBarrelCells.Path = "ecalBarrelCells"
towers.ecalEndcapCells.Path = "emptyCaloCells"
towers.ecalFwdCells.Path = "emptyCaloCells"
towers.hcalBarrelCells.Path = "emptyCaloCells"
towers.hcalExtBarrelCells.Path = "emptyCaloCells"
towers.hcalEndcapCells.Path = "emptyCaloCells"
towers.hcalFwdCells.Path = "emptyCaloCells"
createclusters = CreateCaloClustersSlidingWindow("CreateCaloClusters",
towerTool=towers,
fieldNames = ["system"],
fieldValues = [5],
activeVolumesNumber = 8)
from Configurables import CreateEmptyCaloCellsCollection
createemptycells = CreateEmptyCaloCellsCollection("CreateEmptyCaloCells")
createemptycells.cells.Path = "emptyCaloCells"
#Create calo clusters
from Configurables import CreateCaloClustersSlidingWindow, CaloTowerTool
from GaudiKernel.PhysicalConstants import pi
towers = CaloTowerTool("towers",
deltaEtaTower = 0.01, deltaPhiTower = 2*pi/704.,
ecalBarrelReadoutName = ecalBarrelReadoutNamePhiEta,
ecalEndcapReadoutName = "",
ecalFwdReadoutName = "",
hcalBarrelReadoutName = "",
hcalExtBarrelReadoutName = "",
hcalEndcapReadoutName = "",
hcalFwdReadoutName = "")
towers.ecalBarrelCells.Path = "ECalBarrelCells"
towers.ecalEndcapCells.Path = "emptyCaloCells"
towers.ecalFwdCells.Path = "emptyCaloCells"
towers.hcalBarrelCells.Path = "emptyCaloCells"
towers.hcalExtBarrelCells.Path = "emptyCaloCells"
towers.hcalEndcapCells.Path = "emptyCaloCells"
towers.hcalFwdCells.Path = "emptyCaloCells"
# call pileup tool
# prepare TH2 histogram with pileup per abs(eta)
from Configurables import PreparePileup
from Configurables import CreateCaloCells
createEcalEndcapCells = CreateCaloCells("CreateECalEndcapCells")
createEcalEndcapCells.geometryTool = endcapGeometry
createEcalEndcapCells.doCellCalibration=False # already calibrated
createEcalEndcapCells.addCellNoise=True
createEcalEndcapCells.filterCellNoise=False
createEcalEndcapCells.noiseTool = noiseEndcap
createEcalEndcapCells.hits="newECalEndcapCells"
createEcalEndcapCells.cells=ecalEndcapCellsName+"Noise"
ApplicationMgr().TopAlg +=[createEcalEndcapCells]
#Create calo clusters
from GaudiKernel.PhysicalConstants import pi
from Configurables import CaloTowerTool
towers = CaloTowerTool("towers")
towers.deltaEtaTower = 0.01
towers.deltaPhiTower = 2*pi/704.
towers.ecalBarrelReadoutName = ecalBarrelReadoutName
towers.ecalEndcapReadoutName = ecalEndcapReadoutName
towers.ecalFwdReadoutName = ecalFwdReadoutName
towers.hcalBarrelReadoutName = hcalBarrelReadoutPhiEtaName
towers.hcalExtBarrelReadoutName = hcalExtBarrelReadoutPhiEtaName
towers.hcalEndcapReadoutName = hcalEndcapReadoutName
towers.hcalFwdReadoutName = hcalFwdReadoutName
towers.ecalBarrelCells.Path = ecalBarrelCellsName + "Noise"
towers.ecalEndcapCells.Path = ecalEndcapCellsName + "Noise"
towers.ecalFwdCells.Path = ecalFwdCellsName
towers.hcalBarrelCells.Path = "newHCalBarrelCells"
towers.hcalExtBarrelCells.Path ="newHCalExtBarrelCells"
createemptycells = CreateEmptyCaloCellsCollection("CreateEmptyCaloCells")
createemptycells.cells.Path = "emptyCaloCells"
#Create calo clusters
from Configurables import CreateCaloClustersSlidingWindow
from GaudiKernel.PhysicalConstants import pi
from Configurables import CaloTowerTool
towers = CaloTowerTool("towers",
deltaEtaTower=0.01,
deltaPhiTower=2 * pi / 704.,
ecalBarrelReadoutName=ecalBarrelReadoutName,
ecalEndcapReadoutName=ecalEndcapReadoutName,
ecalFwdReadoutName=ecalFwdReadoutName,
hcalBarrelReadoutName=hcalBarrelReadoutPhiEtaName,
hcalExtBarrelReadoutName=hcalExtBarrelReadoutPhiEtaName,
hcalEndcapReadoutName=hcalEndcapReadoutName,
hcalFwdReadoutName=hcalFwdReadoutName,
OutputLevel=INFO)
towers.ecalBarrelCells.Path = ecalBarrelCellsName + "Noise"
towers.ecalEndcapCells.Path = "emptyCaloCells" #ecalEndcapCellsName + "Noise"
towers.ecalFwdCells.Path = "emptyCaloCells" #ecalFwdCellsName
towers.hcalBarrelCells.Path = "emptyCaloCells"
towers.hcalExtBarrelCells.Path = "emptyCaloCells" # "newHCalExtBarrelCells"
towers.hcalEndcapCells.Path = "emptyCaloCells" #hcalEndcapCellsName
towers.hcalFwdCells.Path = "emptyCaloCells" #hcalFwdCellsName
# Cluster variables
windE = 9
createHcalBarrelCells.addPosition = True
createHcalBarrelCells.addCellNoise = False
createHcalBarrelCells.filterCellNoise = False
createHcalBarrelCells.hits = "HCalBarrelPositionedHits"
createHcalBarrelCells.cells = "HCalBarrelCells"
ApplicationMgr().TopAlg += [createHcalBarrelCells]
# sliding window clustering #FIXME not yet ready for key4hep
#Empty cells for parts of calorimeter not implemented yet
from Configurables import CreateEmptyCaloCellsCollection
createemptycells = CreateEmptyCaloCellsCollection("CreateEmptyCaloCells")
createemptycells.cells.Path = "emptyCaloCells"
ApplicationMgr().TopAlg += [createemptycells]
from Configurables import CaloTowerTool
towers = CaloTowerTool("towers")
towers.deltaEtaTower = 0.01
towers.deltaPhiTower = 2 * constants.pi / 768.
towers.ecalBarrelReadoutName = ecalBarrelReadoutNamePhiEta
towers.ecalEndcapReadoutName = ""
towers.ecalFwdReadoutName = ""
towers.hcalBarrelReadoutName = ""
towers.hcalExtBarrelReadoutName = ""
towers.hcalEndcapReadoutName = ""
towers.hcalFwdReadoutName = ""
towers.ecalBarrelCells.Path = EcalBarrelCellsName
towers.ecalEndcapCells.Path = "emptyCaloCells"
towers.ecalFwdCells.Path = "emptyCaloCells"
towers.hcalBarrelCells.Path = "emptyCaloCells"
towers.hcalExtBarrelCells.Path = "emptyCaloCells"
towers.hcalEndcapCells.Path = "emptyCaloCells"