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 |
|
|
)
|