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,