Exemplo n.º 1
0
        THistSvc().Output = [
            "rec DATAFILE='calibrateClusterNoise_histograms.root' TYP='ROOT' OPT='RECREATE'"
        ]
        THistSvc().PrintAll = True
        THistSvc().AutoSave = True
        THistSvc().AutoFlush = True
        THistSvc().OutputLevel = INFO

        from Configurables import CreateCaloCellPositions
        positionsCalibClusterBarrelNoise = CreateCaloCellPositions(
            "positionsCalibClusterBarrelNoise",
            positionsECalBarrelTool=ECalBcells,
            positionsHCalBarrelTool=HCalBcellVols,
            positionsHCalExtBarrelTool=HCalBcellVols,
            positionsEMECTool=EMECcells,
            positionsHECTool=HECcells,
            positionsEMFwdTool=ECalFwdcells,
            positionsHFwdTool=HCalFwdcells,
            hits="calibCaloClusterBarrelNoiseCells",
            positionedHits="calibCaloClusterBarrelNoiseCellPositions",
            OutputLevel=INFO)

    from Configurables import CreateCaloCellPositions
    positionsClusterBarrelNoise = CreateCaloCellPositions(
        "positionsClusterBarrelNoise",
        positionsECalBarrelTool=ECalBcells,
        positionsHCalBarrelTool=HCalBcellVols,
        positionsHCalExtBarrelTool=HCalBcellVols,
        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(
    "positionsClusterBarrel",
    positionsECalBarrelTool=ECalBcells,
    positionsHCalBarrelTool=HCalBcells,
    positionsHCalExtBarrelTool=HCalExtBcells,
    positionsEMECTool=EMECcells,
    positionsHECTool=HECcells,
    positionsEMFwdTool=ECalFwdcells,
    positionsHFwdTool=HCalFwdcells,
    hits="caloClusterBarrelCells",
    positionedHits="caloClusterBarrelCellPositions",
    OutputLevel=INFO)

out = PodioOutput("out",
                  filename="output_BarrelTopo_50GeVe_3ev.root",
                  OutputLevel=DEBUG)
out.outputCommands = [
    "drop *", "keep GenParticles", "keep GenVertices",
    "keep caloClustersBarrel", "keep caloClusterBarrelCells",
    "keep caloClusterBarrelCellPositions"
]
                                     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,
                                                positionsHCalBarrelTool = HCalBcells,
                                                positionsHCalExtBarrelTool = HCalExtBcells,
                                                positionsEMECTool = EMECcells,
                                                positionsHECTool = HECcells,
                                                positionsEMFwdTool = ECalFwdcells,
                                                positionsHFwdTool = HCalFwdcells,
                                                hits = "caloClusterBarrelCells",
                                                positionedHits = "caloClusterBarrelCellPositions",
                                                OutputLevel = INFO)

out = PodioOutput("out", filename = "output_BarrelTopo_electrNoise_50GeVe_3ev.root", OutputLevel = DEBUG)
out.outputCommands =["drop *", "keep GenParticles", "keep GenVertices", "keep caloClustersBarrel","keep caloClusterBarrelCells", "keep caloClusterBarrelCellPositions"]

#CPU information
from Configurables import AuditorSvc,ChronoAuditor 
chra = ChronoAuditor() 
audsvc = AuditorSvc() 
audsvc.Auditors =[chra] 
podioinput.AuditExecute = True