readNoisyCellsMap = TopoCaloNoisyCells(
"ReadNoisyCellsMap",
fileName=
"/eos/project/f/fccsw-web/testsamples/calo/cellNoise_map_segHcal_constNoiseLevel.root",
OutputLevel=DEBUG)
createTopoClusters = CaloTopoCluster(
"CreateTopoClusters",
TopoClusterInput=createTopoInput,
#expects neighbours map from cellid->vec < neighbourIds >
neigboursTool=readNeighboursMap,
#tool to get noise level per cellid
noiseTool=readNoisyCellsMap,
#cell positions tools for all sub - systems
positionsECalBarrelTool=ECalBcells,
positionsHCalBarrelTool=HCalBcells,
positionsHCalExtBarrelTool=HCalExtBcells,
positionsEMECTool=EMECcells,
positionsHECTool=HECcells,
positionsEMFwdTool=ECalFwdcells,
positionsHFwdTool=HCalFwdcells,
seedSigma=4,
neighbourSigma=0,
lastNeighbourSigma=0,
OutputLevel=DEBUG)
createTopoClusters.clusters.Path = "caloClustersBarrel"
createTopoClusters.clusterCells.Path = "caloClusterBarrelCells"
#Fill a collection of CaloHitPositions for detailed Cluster analysis
from Configurables import CreateCaloCellPositions
positionsClusterBarrel = CreateCaloCellPositions(
#Noise levels per cell
readNoisyCellsMap = TopoCaloNoisyCells("ReadNoisyCellsMap",
fileName = "/afs/cern.ch/work/c/cneubuse/public/FCChh/"
"cellNoise_map_segHcal_electronicsNoiseLevel.root",
OutputLevel = DEBUG)
createTopoClusters = CaloTopoCluster("CreateTopoClusters",
TopoClusterInput = createTopoInput,
#expects neighbours map from cellid->vec < neighbourIds >
neigboursTool = readNeighboursMap,
#tool to get noise level per cellid
noiseTool = readNoisyCellsMap,
#cell positions tools for all sub - systems
positionsECalBarrelTool = ECalBcells,
positionsHCalBarrelTool = HCalBcells,
positionsHCalExtBarrelTool = HCalExtBcells,
positionsEMECTool = EMECcells,
positionsHECTool = HECcells,
positionsEMFwdTool = ECalFwdcells,
positionsHFwdTool = HCalFwdcells,
seedSigma = 4,
neighbourSigma = 2,
lastNeighbourSigma = 0,
OutputLevel = DEBUG)
createTopoClusters.clusters.Path ="caloClustersBarrel"
createTopoClusters.clusterCells.Path = "caloClusterBarrelCells"
#Fill a collection of CaloHitPositions for detailed Cluster analysis
from Configurables import CreateCaloCellPositions
positionsClusterBarrel =CreateCaloCellPositions("positionsClusterBarrel",
positionsECalBarrelTool = ECalBcells,
readNoisyCellsMap = TopoCaloNoisyCells(
"ReadNoisyCellsMap",
fileName=
"/afs/cern.ch/work/c/cneubuse/public/FCChh/cellNoise_map_segHcal.root",
OutputLevel=DEBUG)
# Topo-Cluster Algorithm
# Seed and neighbour thresholds 4 - 2 - 0 w/noise
createTopoClustersNoise = CaloTopoCluster(
"CreateTopoClustersNoise",
TopoClusterInput=createTopoInputNoise,
# expects neighbours map from cellid->vec<neighbourIds>
neigboursTool=readNeighboursMap,
# tool to get noise level per cellid
noiseTool=readNoisyCellsMap,
# cell positions tools for all sub-systems
positionsECalBarrelTool=ECalBcells,
positionsHCalBarrelTool=HCalBcellVols,
seedSigma=4,
neighbourSigma=2,
lastNeighbourSigma=0,
OutputLevel=INFO)
createTopoClustersNoise.clusters.Path = "caloClustersBarrelNoise"
createTopoClustersNoise.clusterCells.Path = "caloClusterBarrelNoiseCells"
if (calib):
from Configurables import CreateCaloClusters
calibrateClustersNoise = CreateCaloClusters(
"CalibrateClustersNoise",
clusters="caloClustersBarrelNoise",
outClusters="calibCaloClustersBarrelNoise",
createTopoInputNoise.hcalExtBarrelCells.Path = "emptyCaloCells"
createTopoInputNoise.hcalEndcapCells.Path = "emptyCaloCells"
createTopoInputNoise.hcalFwdCells.Path = "emptyCaloCells"
readNoisyCellsMap = TopoCaloNoisyCells("ReadNoisyCellsMap",
fileName=inputNoisePerCell)
# Topo-Cluster Algorithm
# Seed and neighbour thresholds 4 - 2 - 0 w/noise
createTopoClustersNoise = CaloTopoCluster(
"CreateTopoClustersNoise",
TopoClusterInput=createTopoInputNoise,
# expects neighbours map from cellid->vec<neighbourIds>
neigboursTool=readNeighboursMap,
# tool to get noise level per cellid
noiseTool=readNoisyCellsMap,
# cell positions tools for all sub-systems
positionsECalBarrelTool=ECalBcells,
positionsHCalBarrelTool=HCalBsegcells,
positionsHCalBarrelNoSegTool=HCalBcellVols,
noSegmentationHCal=noSegmentationHCal,
seedSigma=sigma1,
neighbourSigma=sigma2,
lastNeighbourSigma=sigma3)
createTopoClustersNoise.clusters.Path = "caloClustersBarrelNoise"
createTopoClustersNoise.clusterCells.Path = "caloClusterBarrelNoiseCells"
from Configurables import CreateCaloCellPositions
if resegmentHCal:
positionsClusterBarrelNoise = CreateCaloCellPositions(
"positionsClusterBarrelNoise",
positionsECalBarrelTool=ECalBcells,