hcalBarrelReadoutName=hcalBarrelReadoutName,
hcalExtBarrelReadoutName="",
hcalEndcapReadoutName="",
hcalFwdReadoutName="",
OutputLevel=DEBUG)
createTopoInput.ecalBarrelCells.Path = "ECalBarrelCells"
createTopoInput.ecalEndcapCells.Path = "emptyCaloCells"
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,
ecalFwdReadoutName = "",
hcalBarrelReadoutName = hcalBarrelReadoutName, #"BarHCal_Readout_phieta",
hcalExtBarrelReadoutName = "",
hcalEndcapReadoutName = "",
hcalFwdReadoutName = "",
OutputLevel = DEBUG)
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
hcalBarrelReadoutName="BarHCal_Readout_phieta",
hcalExtBarrelReadoutName="",
hcalEndcapReadoutName="",
hcalFwdReadoutName="",
OutputLevel=DEBUG)
createTopoInput.ecalBarrelCells.Path = "ECalBarrelCellsNoise"
createTopoInput.ecalEndcapCells.Path = "emptyCaloCells"
createTopoInput.ecalFwdCells.Path = "emptyCaloCells"
createTopoInput.hcalBarrelCells.Path = "HCalBarrelCellsNoise"
createTopoInput.hcalExtBarrelCells.Path = "emptyCaloCells"
createTopoInput.hcalEndcapCells.Path = "emptyCaloCells"
createTopoInput.hcalFwdCells.Path = "emptyCaloCells"
readNeighboursMap = TopoCaloNeighbours(
"ReadNeighboursMap",
fileName=
"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",
OutputLevel=DEBUG)
createTopoClusters = CaloTopoCluster(
"CreateTopoClusters",
TopoClusterInput=createTopoInput,
#expects neighbours map from cellid->vec < neighbourIds >
neigboursTool=readNeighboursMap,
##############################################################################################################
####### PREPARE TOPO CLUSTERING : #############
####### PREPARE EMPTY CELL COLLECTIONS #############
####### READ NEIGHBOURS MAP #############
##############################################################################################################
# Fill empty collections for un-used Calo systems
from Configurables import CreateEmptyCaloCellsCollection
createemptycells = CreateEmptyCaloCellsCollection("CreateEmptyCaloCells")
createemptycells.cells.Path = "emptyCaloCells"
from Configurables import TopoCaloNeighbours, TopoCaloNoisyCells
# read the neighbours map
readNeighboursMap = TopoCaloNeighbours(
"ReadNeighboursMap",
fileName=
"/afs/cern.ch/work/c/cneubuse/public/FCChh/neighbours_map_segHcal.root",
OutputLevel=DEBUG)
##############################################################################################################
####### NOISE/NO NOISE TOOL FOR CLUSTER THRESHOLDS #############
##############################################################################################################
if elNoise:
from Configurables import CreateCaloCells, NoiseCaloCellsFromFileTool, TubeLayerPhiEtaCaloTool, CalibrateCaloHitsTool, NoiseCaloCellsFlatTool, NestedVolumesCaloTool
# ECal Barrel noise
noiseBarrel = NoiseCaloCellsFromFileTool(
"NoiseBarrel",
readoutName=ecalBarrelReadoutName,
noiseFileName=ecalBarrelNoisePath,
print 'Use postion tool w/o segmentation : ', noSegmentationHCal
print 'input topo cells : ', inputTopoCellCollectionHCalBarrel
print 'readout name: ', hcalBarrelReadoutName
if puNoise:
print 'pileup noise file name: ', inputPileupNoisePerCell
print 'Full HCal granularity used: ', noSegmentationHCal
##############################################################################################################
####### READ NEIGHBOURS MAP #############
##############################################################################################################
from Configurables import TopoCaloNeighbours, TopoCaloNoisyCells
# read the neighbours map
readNeighboursMap = TopoCaloNeighbours("ReadNeighboursMap",
fileName=inputNeighboursMap,
OutputLevel=DEBUG)
##############################################################################################################
####### NOISE/NO NOISE TOOL FOR CLUSTER THRESHOLDS #############
##############################################################################################################
if elNoise and not puNoise:
# Apply cell thresholds for electronics noise only if no pileup events have been merged
if not resegmentHCal:
inputNoisePerCell = "/afs/cern.ch/work/c/cneubuse/public/FCChh/cellNoise_map_segHcal_electronicsNoiseLevel.root"
if addedPU > 0:
# no offset, but rms from MinBias/simuPU200/
# input cells already resegmented HCal segmentation
inputNoisePerCell = "/afs/cern.ch/work/c/cneubuse/public/FCChh/inBfield/cellNoise_map_electronicsNoiseLevel_forPU" + str(
addedPU) + "_recTopo.root"