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,
예제 #3
0
        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,
예제 #5
0
                                                    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,
예제 #6
0
        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,
예제 #7
0
##############################################################################################################
#######                                       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,