windP = 17
posE = 5
posP = 11
dupE = 7
dupP = 13
finE = 9
finP = 17
threshold = 12
createClusters = CreateCaloClustersSlidingWindow("CreateClusters",
towerTool=towers,
nEtaWindow=windE,
nPhiWindow=windP,
nEtaPosition=posE,
nPhiPosition=posP,
nEtaDuplicates=dupE,
nPhiDuplicates=dupP,
nEtaFinal=finE,
nPhiFinal=finP,
energyThreshold=threshold,
ellipse=True,
OutputLevel=DEBUG)
createClusters.clusters.Path = "CaloClusters"
createClusters.clusterCells.Path = "CaloClusterCells"
out = PodioOutput("out",
filename="output_allCalo_reco.root",
OutputLevel=DEBUG)
out.outputCommands = ["keep *"]
#CPU information
hcalFwdReadoutName=hcalFwdReadoutName)
towersNoise.ecalBarrelCells.Path = "ECalBarrelCellsNoise"
towersNoise.ecalEndcapCells.Path = "ECalEndcapCells"
towersNoise.ecalFwdCells.Path = "ECalFwdCells"
# towersNoise.hcalBarrelCells.Path = "HCalBarrelCells"
# towersNoise.hcalExtBarrelCells.Path = "HCalExtBarrelCells"
towersNoise.hcalBarrelCells.Path = "emptyCaloCells"
towersNoise.hcalExtBarrelCells.Path = "emptyCaloCells"
towersNoise.hcalEndcapCells.Path = "HCalEndcapCells"
towersNoise.hcalFwdCells.Path = "HCalFwdCells"
createclustersNoise = CreateCaloClustersSlidingWindow(
"CreateCaloClustersNoise",
towerTool=towersNoise,
nEtaWindow=7,
nPhiWindow=15,
nEtaPosition=3,
nPhiPosition=11,
nEtaDuplicates=5,
nPhiDuplicates=11,
nEtaFinal=winEta,
nPhiFinal=winPhi,
energyThreshold=enThreshold)
createclustersNoise.clusters.Path = "caloClustersNoise"
from Configurables import CorrectCluster
correctClusters = CorrectCluster("CorrectCluster",
energyAxis=energy,
numLayers=8,
etaValues=[0, 0.25],
presamplerShiftP0=[0.05938, 0.05938],
presamplerShiftP1=[0.0001833, 0.0001833],
presamplerScaleP0=[2.4, 2.4],
presamplerScaleP1=[-0.006838, -0.006838],
from Configurables import CreateCaloClustersSlidingWindow, SingleCaloTowerTool
from GaudiKernel.PhysicalConstants import pi
towers = SingleCaloTowerTool("towers",
deltaEtaTower=0.01,
deltaPhiTower=2 * pi / 629.,
readoutName=ecalReadoutName,
OutputLevel=DEBUG)
towers.cells.Path = "caloCells"
createclusters = CreateCaloClustersSlidingWindow("CreateCaloClusters",
towerTool=towers,
nEtaWindow=7,
nPhiWindow=15,
nEtaPosition=5,
nPhiPosition=11,
nEtaDuplicates=5,
nPhiDuplicates=11,
nEtaFinal=7,
nPhiFinal=15,
energyThreshold=8,
OutputLevel=DEBUG)
createclusters.clusters.Path = "caloClusters"
out = PodioOutput("output",
filename="output_ecalReco_flatNoise_test.root",
OutputLevel=DEBUG)
out.outputCommands = ["keep *"]
#CPU information
from Configurables import AuditorSvc, ChronoAuditor
OutputLevel=DEBUG)
createcells.DataInputs.hits.Path = "ECalHits"
createcells.DataOutputs.cells.Path = "caloCells"
#Create calo clusters
from Configurables import CreateCaloClustersSlidingWindow
from GaudiKernel.PhysicalConstants import pi
createclusters = CreateCaloClustersSlidingWindow("CreateCaloClusters",
readoutName=ecalReadoutName,
fieldNames=ecalFieldNames,
fieldValues=ecalFieldValues,
deltaEtaTower=0.01,
deltaPhiTower=2 * pi / 629.,
nEtaWindow=5,
nPhiWindow=15,
nEtaPosition=5,
nPhiPosition=5,
nEtaDuplicates=5,
nPhiDuplicates=15,
nEtaFinal=5,
nPhiFinal=15,
energyThreshold=7,
OutputLevel=DEBUG)
createclusters.DataInputs.cells.Path = "caloCells"
createclusters.DataOutputs.clusters.Path = "caloClusters"
out = PodioOutput("out",
filename="output_ecalReco_flatNoise_test.root",
OutputLevel=DEBUG)
out.outputCommands = ["keep *"]
towers.hcalEndcapCells.Path = hcalEndcapCellsName
towers.hcalFwdCells.Path = hcalFwdCellsName
# Cluster variables
windE = 9
windP = 17
posE = 5
posP = 11
dupE = 7
dupP = 13
finE = 9
finP = 17
threshold = 12
from Configurables import CreateCaloClustersSlidingWindow
createClusters = CreateCaloClustersSlidingWindow("CreateClusters")
createClusters.towerTool = towers
createClusters.nEtaWindow = windE
createClusters.nPhiWindow = windP
createClusters.nEtaPosition = posE
createClusters.nPhiPosition = posP
createClusters.nEtaDuplicates = dupE
createClusters.nPhiDuplicates = dupP
createClusters.nEtaFinal = finE
createClusters.nPhiFinal = finP
createClusters.energyThreshold = threshold
createClusters.clusters.Path = "CaloClusters"
createClusters.AuditExecute = True
ApplicationMgr().TopAlg +=[createClusters]
out = PodioOutput("out")
createcells.hits.Path = "ECalHits"
createcells.cells.Path = "caloCells"
#Create calo clusters
from Configurables import CreateCaloClustersSlidingWindow, SingleCaloTowerTool
from GaudiKernel.PhysicalConstants import pi
towers = SingleCaloTowerTool("towers",
deltaEtaTower = 0.01, deltaPhiTower = 2*pi/629.,
readoutName = ecalReadoutName,
OutputLevel = DEBUG)
towers.cells.Path = "caloCells"
createclusters = CreateCaloClustersSlidingWindow("CreateCaloClusters",
towerTool = towers,
nEtaWindow = 7, nPhiWindow = 15,
nEtaPosition = 5, nPhiPosition = 11,
nEtaDuplicates = 5, nPhiDuplicates = 11,
nEtaFinal = 7, nPhiFinal = 15,
energyThreshold = 8,
OutputLevel = DEBUG)
createclusters.clusters.Path = "caloClusters"
out = PodioOutput("output", filename = "output_ecalReco_flatNoise_test.root",
OutputLevel = DEBUG)
out.outputCommands = ["keep *"]
#CPU information
from Configurables import AuditorSvc, ChronoAuditor
chra = ChronoAuditor()
audsvc = AuditorSvc()
audsvc.Auditors = [chra]
podioinput.AuditExecute = True
activeVolumeName="LAr_sensitive")
positionhit.DataInputs.caloCells.Path = "caloCells"
positionhit.DataOutputs.caloPositionedHits.Path = "caloCellsPositions"
#Create calo clusters
from Configurables import CreateCaloClustersSlidingWindow
from GaudiKernel.PhysicalConstants import pi
createclusters = CreateCaloClustersSlidingWindow(
"CreateCaloClusters",
readoutName="ECalHitsPhiEta",
fieldNames=["system", "ECAL_Cryo", "bath", "EM_barrel"],
fieldValues=[5, 1, 1, 1],
deltaEtaTower=0.01,
deltaPhiTower=2 * pi / 629.,
nEtaWindow=9,
nPhiWindow=9,
nEtaPosition=7,
nPhiPosition=7,
nEtaDuplicates=9,
nPhiDuplicates=9,
energyThreshold=3,
OutputLevel=DEBUG)
createclusters.DataInputs.cells.Path = "caloCells"
createclusters.DataOutputs.clusters.Path = "caloClusters"
out = PodioOutput("out",
filename="output_ecalReco_test.root",
OutputLevel=DEBUG)
out.outputCommands = ["keep *"]
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"
towers.hcalFwdCells.Path = "emptyCaloCells"
from Configurables import CreateCaloClustersSlidingWindow
createClusters = CreateCaloClustersSlidingWindow("CreateClusters")
# Cluster variables
createClusters.towerTool = towers
createClusters.nEtaWindow = 9
createClusters.nPhiWindow = 17
createClusters.nEtaPosition = 5
createClusters.nPhiPosition = 11
createClusters.nEtaDuplicates = 7
createClusters.nPhiDuplicates = 13
createClusters.nEtaFinal = 9
createClusters.nPhiFinal = 17
# Minimal energy to create a cluster in GeV (FCC-ee detectors have to reconstruct low energy particles)
createClusters.energyThreshold = 0.1
createClusters.attachCells = True
createClusters.clusters.Path = "CaloClusters"
createClusters.clusterCells.Path = "CaloClusterCells"
