createTopoInput.ecalBarrelCells.Path ="selectedECalBarrelCellsNoise"
createTopoInput.ecalEndcapCells.Path ="emptyCaloCells"
createTopoInput.ecalFwdCells.Path ="emptyCaloCells"
createTopoInput.hcalBarrelCells.Path ="selectedHCalBarrelCellsNoise"
createTopoInput.hcalExtBarrelCells.Path ="emptyCaloCells"
createTopoInput.hcalEndcapCells.Path ="emptyCaloCells"
createTopoInput.hcalFwdCells.Path = "emptyCaloCells"
readNeighboursMap =TopoCaloNeighbours("ReadNeighboursMap",
#fileName = "/afs/cern.ch/work/c/cneubuse/public/FCChh/neighbours_map_barrel.root",
fileName = "neighbours_map_fullGran_barrel.root", #http://fccsw.web.cern.ch/fccsw/testsamples/calo/neighbours_map_barrel.root",
OutputLevel = DEBUG)
#Noise levels per cell
readNoisyCellsMap = TopoCaloNoisyCells("ReadNoisyCellsMap",
fileName = "http://fccsw.web.cern.ch/fccsw/testsamples/calo/cellNoise_map_electronicsNoiseLevel.root",
# fileName = "http://fccsw.web.cern.ch/fccsw/testsamples/calo/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,
createTopoInput.ecalFwdCells.Path = "emptyCaloCells"
createTopoInput.hcalBarrelCells.Path = "HCalBarrelCells"
createTopoInput.hcalExtBarrelCells.Path = "emptyCaloCells"
createTopoInput.hcalEndcapCells.Path = "emptyCaloCells"
createTopoInput.hcalFwdCells.Path = "emptyCaloCells"
readNeighboursMap = TopoCaloNeighbours(
"ReadNeighboursMap",
fileName=
"/eos/project/f/fccsw-web/testsamples/calo/neighbours_map_segHcal.root",
OutputLevel=DEBUG)
#Thresholds estimated from atlas, without noise !!!
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,
hcalBarrelReadoutName=hcalBarrelReadoutVolume,
hcalExtBarrelReadoutName="",
hcalEndcapReadoutName="",
hcalFwdReadoutName="",
OutputLevel=DEBUG)
createTopoInputNoise.ecalBarrelCells.Path = "ECalBarrelCellsNoise"
createTopoInputNoise.ecalEndcapCells.Path = "emptyCaloCells"
createTopoInputNoise.ecalFwdCells.Path = "emptyCaloCells"
createTopoInputNoise.hcalBarrelCells.Path = "HCalBarrelCellsNoise"
createTopoInputNoise.hcalExtBarrelCells.Path = "emptyCaloCells"
createTopoInputNoise.hcalEndcapCells.Path = "emptyCaloCells"
createTopoInputNoise.hcalFwdCells.Path = "emptyCaloCells"
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,
createTopoInput.ecalFwdCells.Path = "emptyCaloCells"
createTopoInput.hcalBarrelCells.Path = "HCalBarrelCellsNoise"
createTopoInput.hcalExtBarrelCells.Path = "emptyCaloCells"
createTopoInput.hcalEndcapCells.Path = "emptyCaloCells"
createTopoInput.hcalFwdCells.Path = "emptyCaloCells"
readNeighboursMap = TopoCaloNeighbours(
"ReadNeighboursMap",
fileName=
"/afs/cern.ch/work/c/cneubuse/public/FCChh/neighbours_map_segHcal.root",
OutputLevel=DEBUG)
#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,
ecalFwdReadoutName = "",
hcalBarrelReadoutName = hcalBarrelReadoutName,
hcalExtBarrelReadoutName = "",
hcalEndcapReadoutName = "",
hcalFwdReadoutName = "",
OutputLevel = INFO)
createTopoInputNoise.ecalBarrelCells.Path = "ECalBarrelCellsNoise"
createTopoInputNoise.ecalEndcapCells.Path = "emptyCaloCells"
createTopoInputNoise.ecalFwdCells.Path = "emptyCaloCells"
createTopoInputNoise.hcalBarrelCells.Path = "HCalBarrelCellsNoise"
createTopoInputNoise.hcalExtBarrelCells.Path = "emptyCaloCells"
createTopoInputNoise.hcalEndcapCells.Path = "emptyCaloCells"
createTopoInputNoise.hcalFwdCells.Path = "emptyCaloCells"
readNoisyCellsMap = TopoCaloNoisyCells("ReadNoisyCellsMap",
fileName = "/afs/cern.ch/work/c/cneubuse/public/FCChh/cellNoise_map_segHcal_electronicsNoiseLevel.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 = sigma1,
neighbourSigma = sigma2,
ecalBarrelReadoutName=ecalBarrelReadoutName,
ecalEndcapReadoutName="",
ecalFwdReadoutName="",
hcalBarrelReadoutName=hcalBarrelReadoutName,
hcalExtBarrelReadoutName="",
hcalEndcapReadoutName="",
hcalFwdReadoutName="")
createTopoInputNoise.ecalBarrelCells.Path = cellInputECal
createTopoInputNoise.ecalEndcapCells.Path = "emptyCaloCells"
createTopoInputNoise.ecalFwdCells.Path = "emptyCaloCells"
createTopoInputNoise.hcalBarrelCells.Path = cellInputHCal
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,
##############################################################################################################
####### REBASE PU TO 0 #############
##############################################################################################################
# Need to resegment HCal cells
# pileup noise levels always in DeltaEta=0.025 cells
#Configure tools for calo cell positions
inputPileupNoisePerCell = "/afs/cern.ch/work/c/cneubuse/public/FCChh/inBfield/cellNoise_map_electronicsNoiseLevel_forPU" + str(
simargs.pileup) + ".root"
if resegmentHCal:
inputPileupNoisePerCell = "/afs/cern.ch/work/c/cneubuse/public/FCChh/inBfield/resegmentedHCal/cellNoise_map_electronicsNoiseLevel_forPU" + str(
simargs.pileup) + ".root"
from Configurables import TopoCaloNoisyCells
readNoisyCellsMap = TopoCaloNoisyCells("ReadNoisyCellsMap",
fileName=inputPileupNoisePerCell)
rebaseEcalBarrelCells = CreateCaloCells("RebaseECalBarrelCells",
doCellCalibration=False,
addCellNoise=False,
filterCellNoise=False,
recalibrateBaseline=True,
readCellNoiseTool=readNoisyCellsMap,
hits=ecalBarrelOutput1,
cells="mergedECalBarrelCells")
rebaseHcalBarrelCells = CreateCaloCells("RebaseHCalBarrelCells",
doCellCalibration=False,
addCellNoise=False,
filterCellNoise=False,
recalibrateBaseline=True,