ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/CmsHi/JetAnalysis/python/ExtraEGammaReco_cff.py
Revision: 1.1.2.1
Committed: Tue Aug 9 14:27:12 2011 UTC (13 years, 9 months ago) by frankma
Content type: text/x-python
Branch: cmssw39x_branch
Changes since 1.1: +46 -0 lines
Log Message:
add photon sequences

File Contents

# User Rev Content
1 frankma 1.1.2.1 import FWCore.ParameterSet.Config as cms
2    
3     # reco
4     from RecoHI.HiEgammaAlgos.HiEgamma_cff import *
5    
6     # detector responce
7     from CmsHi.PhotonAnalysis.isoConeInspector_cfi import *
8     from CmsHi.PhotonAnalysis.ecalHistProducer_cfi import *
9     from CmsHi.PhotonAnalysis.SpikeInspector_cfi import *
10    
11     photonObj = "cleanPhotons"
12    
13     # spike cleaner for island superclsters
14     cleanPhotons.maxHoverEBarrel = cms.double(100)
15     PhotonIDProd.photonProducer = cms.string(photonObj)
16     gamIsoDepositTk.src = cms.InputTag(photonObj)
17     gamIsoDepositEcalFromHits.src = cms.InputTag(photonObj)
18     gamIsoDepositHcalFromTowers.src = cms.InputTag(photonObj)
19     gamIsoDepositHcalDepth1FromTowers.src = cms.InputTag(photonObj)
20     gamIsoDepositHcalDepth2FromTowers.src = cms.InputTag(photonObj)
21    
22     # clean photon filter
23     goodPhotons = cms.EDFilter("PhotonSelector",
24     src = cms.InputTag("cleanPhotons"),
25     cut = cms.string("pt > 18 & hadronicOverEm < 0.2 & abs(eta) < 1.45"),
26     )
27     filterGoodPhotons = cms.EDFilter("PhotonCountFilter",
28     src = cms.InputTag("goodPhotons"),
29     minNumber = cms.uint32(0), # tmp solution for both gamma and jets
30     maxNumber = cms.uint32(999999),
31     )
32    
33     barrelPhotonFilter = cms.Sequence ( goodPhotons * filterGoodPhotons )
34    
35     # pat
36     from PhysicsTools.PatAlgos.producersHeavyIons.heavyIonPhotons_cff import *
37    
38     patPhotons.photonSource = cms.InputTag(photonObj)
39     photonMatch.src = cms.InputTag(photonObj)
40    
41     # final seq
42     photon_extra_reco = cms.Sequence(
43     hiEcalClusteringSequence *
44     hiPhotonCleaningSequence *
45     barrelPhotonFilter
46     )