HiEgammaAlgos/python/HiEgamma_cff.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 # clustering sequence
0004 from RecoHI.HiEgammaAlgos.HiIslandClusteringSequence_cff import *
0005 from RecoEcal.EgammaClusterProducers.hybridClusteringSequence_cff import *
0006 from RecoEcal.EgammaClusterProducers.multi5x5ClusteringSequence_cff import *
0007 from RecoEcal.EgammaClusterProducers.multi5x5PreshowerClusteringSequence_cff import *
0008 from RecoEcal.EgammaClusterProducers.preshowerClusteringSequence_cff import *
0009 from RecoHI.HiEgammaAlgos.HiIsolationCommonParameters_cff import *
0010 from RecoEcal.EgammaClusterProducers.particleFlowSuperClusterECAL_cfi import *
0011
0012 particleFlowSuperClusterECAL.regressionConfig.vertexCollection = 'hiSelectedVertex'
0013
0014 hiEcalClusteringTask = cms.Task(islandClusteringTask,hybridClusteringTask,multi5x5ClusteringTask,multi5x5PreshowerClusteringTask,preshowerClusteringTask,particleFlowSuperClusterECAL)
0015 hiEcalClusteringSequence = cms.Sequence(hiEcalClusteringTask)
0016
0017
0018 # reco photon producer
0019 from RecoEgamma.EgammaPhotonProducers.photonSequence_cff import *
0020
0021 # use island for the moment
0022 photonCore.scHybridBarrelProducer = "correctedIslandBarrelSuperClusters"
0023 photonCore.scIslandEndcapProducer = "correctedIslandEndcapSuperClusters"
0024 photonCore.minSCEt    = 8.0
0025 photons.minSCEtBarrel = 5.0
0026 photons.minSCEtEndcap = 15.0
0027 photons.minR9Barrel   = 10.  #0.94
0028 photons.minR9Endcap   = 10.  #0.95
0029 photons.maxHoverEEndcap = 0.5   #0.5
0030 photons.maxHoverEBarrel = 0.99  #0.5
0031 photons.primaryVertexProducer = 'hiSelectedVertex' # replace the primary vertex
0032 photons.isolationSumsCalculatorSet.trackProducer = isolationInputParameters.track    # cms.InputTag("highPurityTracks")
0033
0034 from RecoHI.HiEgammaAlgos.photonIsolationHIProducer_cfi import photonIsolationHIProducer
0035 hiPhotonTask = cms.Task(photonTask,photonIsolationHIProducer)
0036 hiPhotonSequence = cms.Sequence(hiPhotonTask)
0037
0038 # HI Ecal reconstruction
0039 hiEcalClustersTask = cms.Task(hiEcalClusteringTask)
0040 hiEcalClusters = cms.Sequence(hiEcalClustersTask)
0041 hiEgammaTask = cms.Task(hiPhotonTask)
0042 hiEgammaSequence = cms.Sequence(hiEgammaTask)
0043
0044 # HI Spike Clean Sequence
0045 import RecoHI.HiEgammaAlgos.hiSpikeCleaner_cfi
0046 hiSpikeCleanedSC = RecoHI.HiEgammaAlgos.hiSpikeCleaner_cfi.hiSpikeCleaner.clone()
0047 cleanPhotonCore = photonCore.clone(
0048     scHybridBarrelProducer = "hiSpikeCleanedSC"
0049 )
0050 cleanPhotons = photons.clone(
0051     photonCoreProducer = cms.InputTag("cleanPhotonCore")
0052 )
0053
0054 hiPhotonCleaningTask = cms.Task(hiSpikeCleanedSC,
0055                                 cleanPhotonCore,
0056                                 cleanPhotons)
0057 hiPhotonCleaningSequence = cms.Sequence(hiPhotonCleaningTask)