1 import FWCore.ParameterSet.Config
as cms
8 from RecoEgamma.EgammaIsolationAlgos.egammaHBHERecHitThreshold_cff
import egammaHBHERecHit
10 from RecoEgamma.EgammaIsolationAlgos.egammaEcalPFClusterIsolationProducerRecoPhoton_cfi
import egammaEcalPFClusterIsolationProducerRecoPhoton
11 from RecoEgamma.EgammaIsolationAlgos.egammaHcalPFClusterIsolationProducerRecoPhoton_cfi
import egammaHcalPFClusterIsolationProducerRecoPhoton
16 photons = cms.EDProducer(
"GEDPhotonProducer",
17 photonProducer = cms.InputTag(
"photonCore"),
18 reconstructionStep = cms.string(
"tmp"),
19 outputPhotonCollection = cms.string(
""),
20 pfEgammaCandidates = cms.InputTag(
""),
21 valueMapPhotons = cms.string(
""),
23 regressionWeightsFromDB = cms.bool(
True),
24 energyRegressionWeightsFileLocation = cms.string(
'/afs/cern.ch/user/b/bendavid/cmspublic/regweights/gbrph.root'),
25 energyRegressionWeightsDBLocation = cms.string(
'wgbrph'),
26 superClusterEnergyCorrFunction = cms.string(
"EcalClusterEnergyCorrection"),
27 superClusterEnergyErrorFunction = cms.string(
"EcalClusterEnergyUncertainty"),
28 superClusterCrackEnergyCorrFunction = cms.string(
"EcalClusterCrackCorrection"),
29 photonEcalEnergyCorrFunction = cms.string(
"EcalClusterEnergyCorrectionObjectSpecific"),
31 candidateP4type = cms.string(
"fromEcalEnergy"),
32 isolationSumsCalculatorSet = cms.PSet(isolationSumsCalculator),
33 mipVariableSet = cms.PSet(mipVariable),
34 mvaBasedHaloVariableSet = cms.PSet(mvaHaloVariable),
35 usePrimaryVertex = cms.bool(
True),
36 primaryVertexProducer = cms.InputTag(
'offlinePrimaryVerticesWithBS'),
37 posCalc_t0_endcPresh = cms.double(3.6),
38 posCalc_logweight = cms.bool(
True),
39 posCalc_w0 = cms.double(4.2),
40 hbheInstance = cms.string(
''),
41 posCalc_t0_endc = cms.double(6.3),
42 barrelEcalHits = cms.InputTag(
"ecalRecHit",
"EcalRecHitsEB"),
43 hbheModule = cms.string(
'hbhereco'),
44 endcapEcalHits = cms.InputTag(
"ecalRecHit",
"EcalRecHitsEE"),
45 preshowerHits = cms.InputTag(
"ecalPreshowerRecHit",
"EcalRecHitsES"),
46 runMIPTagger = cms.bool(
True),
47 runMVABasedHaloTagger = cms.bool(
False),
48 highEt = cms.double(100.),
49 minR9Barrel = cms.double(0.94),
50 minR9Endcap = cms.double(0.95),
51 multThresEB = cms.double(1.00),
52 multThresEE = cms.double(1.25),
53 hOverEConeSize = cms.double(0.15),
54 hbheRecHits = egammaHBHERecHit.hbheRecHits,
55 recHitEThresholdHB = egammaHBHERecHit.recHitEThresholdHB,
56 recHitEThresholdHE = egammaHBHERecHit.recHitEThresholdHE,
57 maxHcalRecHitSeverity = egammaHBHERecHit.maxHcalRecHitSeverity,
58 hcalRun2EffDepth = cms.bool(
False),
59 posCalc_x0 = cms.double(0.89),
60 posCalc_t0_barl = cms.double(7.7),
61 minSCEtBarrel = cms.double(10.0),
62 minSCEtEndcap = cms.double(10.0),
63 maxHoverEBarrel = cms.double(0.5),
64 maxHoverEEndcap = cms.double(0.5),
65 ecalRecHitSumEtOffsetBarrel = cms.double(999999999),
66 ecalRecHitSumEtSlopeBarrel = cms.double(0.),
67 ecalRecHitSumEtOffsetEndcap = cms.double(999999999),
68 ecalRecHitSumEtSlopeEndcap = cms.double(0.),
69 hcalRecHitSumEtOffsetBarrel = cms.double(999999999),
70 hcalRecHitSumEtSlopeBarrel = cms.double(0.),
71 hcalRecHitSumEtOffsetEndcap = cms.double(999999999),
72 hcalRecHitSumEtSlopeEndcap = cms.double(0.),
73 nTrackSolidConeBarrel =cms.double(999999999),
74 nTrackSolidConeEndcap =cms.double(999999999),
75 nTrackHollowConeBarrel =cms.double(999999999),
76 nTrackHollowConeEndcap =cms.double(999999999),
77 trackPtSumSolidConeBarrel =cms.double(999999999),
78 trackPtSumSolidConeEndcap =cms.double(999999999),
79 trackPtSumHollowConeBarrel =cms.double(999999999),
80 trackPtSumHollowConeEndcap =cms.double(999999999),
81 sigmaIetaIetaCutBarrel=cms.double(999999999),
82 sigmaIetaIetaCutEndcap=cms.double(999999999),
83 posCalcParameters = cms.PSet( T0_barl = cms.double(7.4),
84 T0_endc = cms.double(6.3),
85 T0_endcPresh = cms.double(3.6),
86 LogWeighted = cms.bool(
True),
91 RecHitFlagToBeExcludedEB = cleanedHybridSuperClusters.RecHitFlagToBeExcluded,
92 RecHitSeverityToBeExcludedEB = cleanedHybridSuperClusters.RecHitSeverityToBeExcluded,
93 RecHitFlagToBeExcludedEE = multi5x5BasicClustersCleaned.RecHitFlagToBeExcluded,
94 RecHitSeverityToBeExcludedEE = cleanedHybridSuperClusters.RecHitSeverityToBeExcluded,
95 checkHcalStatus = cms.bool(
True),
96 PhotonDNNPFid = cms.PSet(
97 enabled = cms.bool(
False),
98 inputTensorName = cms.string(
"FirstLayer_input"),
99 outputTensorName = cms.string(
"sequential/FinalLayer/Sigmoid"),
100 modelsFiles = cms.vstring(
101 'RecoEgamma/PhotonIdentification/data/Photon_PFID_dnn/v1/EB/EB_modelDNN.pb',
102 'RecoEgamma/PhotonIdentification/data/Photon_PFID_dnn/v1/EE/EE_modelDNN.pb'),
103 scalersFiles = cms.vstring(
104 'RecoEgamma/PhotonIdentification/data/Photon_PFID_dnn/v1/EB/EB_scaler.txt',
105 'RecoEgamma/PhotonIdentification/data/Photon_PFID_dnn/v1/EE/EE_scaler.txt' 107 outputDim = cms.vuint32(1,1),
108 useEBModelInGap = cms.bool(
True)
110 pfECALClusIsolCfg = cms.PSet(
111 pfClusterProducer = egammaEcalPFClusterIsolationProducerRecoPhoton.pfClusterProducer,
112 drMax = egammaEcalPFClusterIsolationProducerRecoPhoton.drMax,
113 drVetoBarrel = egammaEcalPFClusterIsolationProducerRecoPhoton.drVetoBarrel,
114 drVetoEndcap = egammaEcalPFClusterIsolationProducerRecoPhoton.drVetoEndcap,
115 etaStripBarrel = egammaEcalPFClusterIsolationProducerRecoPhoton.etaStripBarrel,
116 etaStripEndcap = egammaEcalPFClusterIsolationProducerRecoPhoton.etaStripEndcap,
117 energyBarrel = egammaEcalPFClusterIsolationProducerRecoPhoton.energyBarrel,
118 energyEndcap = egammaEcalPFClusterIsolationProducerRecoPhoton.energyEndcap
121 pfHCALClusIsolCfg = cms.PSet(
123 pfClusterProducerHCAL = egammaHcalPFClusterIsolationProducerRecoPhoton.pfClusterProducerHCAL,
124 useHF = egammaHcalPFClusterIsolationProducerRecoPhoton.useHF,
125 pfClusterProducerHFEM = egammaHcalPFClusterIsolationProducerRecoPhoton.pfClusterProducerHFEM,
126 pfClusterProducerHFHAD = egammaHcalPFClusterIsolationProducerRecoPhoton.pfClusterProducerHFHAD,
127 drMax = egammaHcalPFClusterIsolationProducerRecoPhoton.drMax,
128 drVetoBarrel = egammaHcalPFClusterIsolationProducerRecoPhoton.drVetoBarrel,
129 drVetoEndcap = egammaHcalPFClusterIsolationProducerRecoPhoton.drVetoEndcap,
130 etaStripBarrel = egammaHcalPFClusterIsolationProducerRecoPhoton.etaStripBarrel,
131 etaStripEndcap = egammaHcalPFClusterIsolationProducerRecoPhoton.etaStripEndcap,
132 energyBarrel = egammaHcalPFClusterIsolationProducerRecoPhoton.energyBarrel,
133 energyEndcap = egammaHcalPFClusterIsolationProducerRecoPhoton.energyEndcap,
134 useEt = egammaHcalPFClusterIsolationProducerRecoPhoton.useEt,
140 photonsHGC = photons.clone(
141 photonProducer =
'photonCoreHGC' 144 islandPhotons = cms.EDProducer(
"PhotonProducer",
145 photonCoreProducer = cms.InputTag(
"islandPhotonCore"),
146 regressionWeightsFromDB = cms.bool(
True),
147 energyRegressionWeightsFileLocation = cms.string(
'/afs/cern.ch/user/b/bendavid/cmspublic/regweights/gbrph.root'),
148 energyRegressionWeightsDBLocation = cms.string(
'wgbrph'),
149 superClusterEnergyCorrFunction = cms.string(
"EcalClusterEnergyCorrection"),
150 superClusterEnergyErrorFunction = cms.string(
"EcalClusterEnergyUncertainty"),
151 superClusterCrackEnergyCorrFunction = cms.string(
"EcalClusterCrackCorrection"),
152 photonEcalEnergyCorrFunction = cms.string(
"EcalClusterEnergyCorrectionObjectSpecific"),
153 candidateP4type = cms.string(
"fromEcalEnergy"),
154 isolationSumsCalculatorSet = cms.PSet(isolationSumsCalculator),
155 mipVariableSet = cms.PSet(mipVariable),
156 usePrimaryVertex = cms.bool(
True),
157 primaryVertexProducer = cms.InputTag(
'offlinePrimaryVerticesWithBS'),
158 posCalc_t0_endcPresh = cms.double(3.6),
159 posCalc_logweight = cms.bool(
True),
160 posCalc_w0 = cms.double(4.2),
161 hbheInstance = cms.string(
''),
162 posCalc_t0_endc = cms.double(6.3),
163 barrelEcalHits = cms.InputTag(
"ecalRecHit",
"EcalRecHitsEB"),
164 hbheModule = cms.string(
'hbhereco'),
165 endcapEcalHits = cms.InputTag(
"ecalRecHit",
"EcalRecHitsEE"),
166 runMIPTagger = cms.bool(
True),
167 runMVABasedHaloTagger = cms.bool(
False),
168 highEt = cms.double(100.),
169 minR9Barrel = cms.double(10.0),
170 minR9Endcap = cms.double(10.0),
171 hOverEConeSize = cms.double(0.15),
172 hbheRecHits = egammaHBHERecHit.hbheRecHits,
173 recHitEThresholdHB = egammaHBHERecHit.recHitEThresholdHB,
174 recHitEThresholdHE = egammaHBHERecHit.recHitEThresholdHE,
175 maxHcalRecHitSeverity = egammaHBHERecHit.maxHcalRecHitSeverity,
176 hcalRun2EffDepth = cms.bool(
False),
177 posCalc_x0 = cms.double(0.89),
178 posCalc_t0_barl = cms.double(7.7),
179 minSCEtBarrel = cms.double(5.0),
180 minSCEtEndcap = cms.double(15.0),
181 maxHoverEBarrel = cms.double(0.99),
182 maxHoverEEndcap = cms.double(0.5),
183 ecalRecHitSumEtOffsetBarrel = cms.double(999999999),
184 ecalRecHitSumEtSlopeBarrel = cms.double(0.),
185 ecalRecHitSumEtOffsetEndcap = cms.double(999999999),
186 ecalRecHitSumEtSlopeEndcap = cms.double(0.),
187 hcalRecHitSumEtOffsetBarrel = cms.double(999999999),
188 hcalRecHitSumEtSlopeBarrel = cms.double(0.),
189 hcalRecHitSumEtOffsetEndcap = cms.double(999999999),
190 hcalRecHitSumEtSlopeEndcap = cms.double(0.),
191 nTrackSolidConeBarrel =cms.double(999999999),
192 nTrackSolidConeEndcap =cms.double(999999999),
193 nTrackHollowConeBarrel =cms.double(999999999),
194 nTrackHollowConeEndcap =cms.double(999999999),
195 trackPtSumSolidConeBarrel =cms.double(999999999),
196 trackPtSumSolidConeEndcap =cms.double(999999999),
197 trackPtSumHollowConeBarrel =cms.double(999999999),
198 trackPtSumHollowConeEndcap =cms.double(999999999),
199 sigmaIetaIetaCutBarrel=cms.double(999999999),
200 sigmaIetaIetaCutEndcap=cms.double(999999999),
201 posCalcParameters = cms.PSet( T0_barl = cms.double(7.4),
202 T0_endc = cms.double(6.3),
203 T0_endcPresh = cms.double(3.6),
204 LogWeighted = cms.bool(
True),
205 W0 = cms.double(4.2),
206 X0 = cms.double(0.89)
209 RecHitFlagToBeExcludedEB = cleanedHybridSuperClusters.RecHitFlagToBeExcluded,
210 RecHitSeverityToBeExcludedEB = cleanedHybridSuperClusters.RecHitSeverityToBeExcluded,
211 RecHitFlagToBeExcludedEE = multi5x5BasicClustersCleaned.RecHitFlagToBeExcluded,
212 RecHitSeverityToBeExcludedEE = cleanedHybridSuperClusters.RecHitSeverityToBeExcluded,
215 from Configuration.ProcessModifiers.egamma_lowPt_exclusive_cff
import egamma_lowPt_exclusive
216 egamma_lowPt_exclusive.toModify(photons,
1.8.14