1 import FWCore.ParameterSet.Config
as cms
7 from RecoEgamma.EgammaIsolationAlgos.egammaHBHERecHitThreshold_cff
import egammaHBHERecHit
11 photons = cms.EDProducer(
"GEDPhotonProducer",
12 photonProducer = cms.InputTag(
"photonCore"),
13 reconstructionStep = cms.string(
"tmp"),
14 outputPhotonCollection = cms.string(
""),
15 pfEgammaCandidates = cms.InputTag(
""),
16 valueMapPhotons = cms.string(
""),
18 regressionWeightsFromDB = cms.bool(
True),
19 energyRegressionWeightsFileLocation = cms.string(
'/afs/cern.ch/user/b/bendavid/cmspublic/regweights/gbrph.root'),
20 energyRegressionWeightsDBLocation = cms.string(
'wgbrph'),
21 superClusterEnergyCorrFunction = cms.string(
"EcalClusterEnergyCorrection"),
22 superClusterEnergyErrorFunction = cms.string(
"EcalClusterEnergyUncertainty"),
23 superClusterCrackEnergyCorrFunction = cms.string(
"EcalClusterCrackCorrection"),
24 photonEcalEnergyCorrFunction = cms.string(
"EcalClusterEnergyCorrectionObjectSpecific"),
26 candidateP4type = cms.string(
"fromEcalEnergy"),
27 isolationSumsCalculatorSet = cms.PSet(isolationSumsCalculator),
28 mipVariableSet = cms.PSet(mipVariable),
29 usePrimaryVertex = cms.bool(
True),
30 primaryVertexProducer = cms.InputTag(
'offlinePrimaryVerticesWithBS'),
31 posCalc_t0_endcPresh = cms.double(3.6),
32 posCalc_logweight = cms.bool(
True),
33 posCalc_w0 = cms.double(4.2),
34 hbheInstance = cms.string(
''),
35 posCalc_t0_endc = cms.double(6.3),
36 barrelEcalHits = cms.InputTag(
"ecalRecHit",
"EcalRecHitsEB"),
37 hbheModule = cms.string(
'hbhereco'),
38 endcapEcalHits = cms.InputTag(
"ecalRecHit",
"EcalRecHitsEE"),
39 preshowerHits = cms.InputTag(
"ecalPreshowerRecHit",
"EcalRecHitsES"),
40 runMIPTagger = cms.bool(
True),
41 highEt = cms.double(100.),
42 minR9Barrel = cms.double(0.94),
43 minR9Endcap = cms.double(0.95),
44 multThresEB = cms.double(1.00),
45 multThresEE = cms.double(1.25),
46 hOverEConeSize = cms.double(0.15),
47 hbheRecHits = egammaHBHERecHit.hbheRecHits,
48 recHitEThresholdHB = egammaHBHERecHit.recHitEThresholdHB,
49 recHitEThresholdHE = egammaHBHERecHit.recHitEThresholdHE,
50 maxHcalRecHitSeverity = egammaHBHERecHit.maxHcalRecHitSeverity,
51 hcalRun2EffDepth = cms.bool(
False),
52 posCalc_x0 = cms.double(0.89),
53 posCalc_t0_barl = cms.double(7.7),
54 minSCEtBarrel = cms.double(10.0),
55 minSCEtEndcap = cms.double(10.0),
56 maxHoverEBarrel = cms.double(0.5),
57 maxHoverEEndcap = cms.double(0.5),
58 ecalRecHitSumEtOffsetBarrel = cms.double(999999999),
59 ecalRecHitSumEtSlopeBarrel = cms.double(0.),
60 ecalRecHitSumEtOffsetEndcap = cms.double(999999999),
61 ecalRecHitSumEtSlopeEndcap = cms.double(0.),
62 hcalRecHitSumEtOffsetBarrel = cms.double(999999999),
63 hcalRecHitSumEtSlopeBarrel = cms.double(0.),
64 hcalRecHitSumEtOffsetEndcap = cms.double(999999999),
65 hcalRecHitSumEtSlopeEndcap = cms.double(0.),
66 nTrackSolidConeBarrel =cms.double(999999999),
67 nTrackSolidConeEndcap =cms.double(999999999),
68 nTrackHollowConeBarrel =cms.double(999999999),
69 nTrackHollowConeEndcap =cms.double(999999999),
70 trackPtSumSolidConeBarrel =cms.double(999999999),
71 trackPtSumSolidConeEndcap =cms.double(999999999),
72 trackPtSumHollowConeBarrel =cms.double(999999999),
73 trackPtSumHollowConeEndcap =cms.double(999999999),
74 sigmaIetaIetaCutBarrel=cms.double(999999999),
75 sigmaIetaIetaCutEndcap=cms.double(999999999),
76 posCalcParameters = cms.PSet( T0_barl = cms.double(7.4),
77 T0_endc = cms.double(6.3),
78 T0_endcPresh = cms.double(3.6),
79 LogWeighted = cms.bool(
True),
84 RecHitFlagToBeExcludedEB = cleanedHybridSuperClusters.RecHitFlagToBeExcluded,
85 RecHitSeverityToBeExcludedEB = cleanedHybridSuperClusters.RecHitSeverityToBeExcluded,
86 RecHitFlagToBeExcludedEE = multi5x5BasicClustersCleaned.RecHitFlagToBeExcluded,
87 RecHitSeverityToBeExcludedEE = cleanedHybridSuperClusters.RecHitSeverityToBeExcluded,
88 checkHcalStatus = cms.bool(
True)
91 photonsFromMultiCl = photons.clone(
92 photonProducer =
'photonCoreFromMultiCl'
95 islandPhotons = cms.EDProducer(
"PhotonProducer",
96 photonCoreProducer = cms.InputTag(
"islandPhotonCore"),
97 regressionWeightsFromDB = cms.bool(
True),
98 energyRegressionWeightsFileLocation = cms.string(
'/afs/cern.ch/user/b/bendavid/cmspublic/regweights/gbrph.root'),
99 energyRegressionWeightsDBLocation = cms.string(
'wgbrph'),
100 superClusterEnergyCorrFunction = cms.string(
"EcalClusterEnergyCorrection"),
101 superClusterEnergyErrorFunction = cms.string(
"EcalClusterEnergyUncertainty"),
102 superClusterCrackEnergyCorrFunction = cms.string(
"EcalClusterCrackCorrection"),
103 photonEcalEnergyCorrFunction = cms.string(
"EcalClusterEnergyCorrectionObjectSpecific"),
104 candidateP4type = cms.string(
"fromEcalEnergy"),
105 isolationSumsCalculatorSet = cms.PSet(isolationSumsCalculator),
106 mipVariableSet = cms.PSet(mipVariable),
107 usePrimaryVertex = cms.bool(
True),
108 primaryVertexProducer = cms.InputTag(
'offlinePrimaryVerticesWithBS'),
109 posCalc_t0_endcPresh = cms.double(3.6),
110 posCalc_logweight = cms.bool(
True),
111 posCalc_w0 = cms.double(4.2),
112 hbheInstance = cms.string(
''),
113 posCalc_t0_endc = cms.double(6.3),
114 barrelEcalHits = cms.InputTag(
"ecalRecHit",
"EcalRecHitsEB"),
115 hbheModule = cms.string(
'hbhereco'),
116 endcapEcalHits = cms.InputTag(
"ecalRecHit",
"EcalRecHitsEE"),
117 runMIPTagger = cms.bool(
True),
118 highEt = cms.double(100.),
119 minR9Barrel = cms.double(10.0),
120 minR9Endcap = cms.double(10.0),
121 hOverEConeSize = cms.double(0.15),
122 hbheRecHits = egammaHBHERecHit.hbheRecHits,
123 recHitEThresholdHB = egammaHBHERecHit.recHitEThresholdHB,
124 recHitEThresholdHE = egammaHBHERecHit.recHitEThresholdHE,
125 maxHcalRecHitSeverity = egammaHBHERecHit.maxHcalRecHitSeverity,
126 hcalRun2EffDepth = cms.bool(
False),
127 posCalc_x0 = cms.double(0.89),
128 posCalc_t0_barl = cms.double(7.7),
129 minSCEtBarrel = cms.double(5.0),
130 minSCEtEndcap = cms.double(15.0),
131 maxHoverEBarrel = cms.double(0.99),
132 maxHoverEEndcap = cms.double(0.5),
133 ecalRecHitSumEtOffsetBarrel = cms.double(999999999),
134 ecalRecHitSumEtSlopeBarrel = cms.double(0.),
135 ecalRecHitSumEtOffsetEndcap = cms.double(999999999),
136 ecalRecHitSumEtSlopeEndcap = cms.double(0.),
137 hcalRecHitSumEtOffsetBarrel = cms.double(999999999),
138 hcalRecHitSumEtSlopeBarrel = cms.double(0.),
139 hcalRecHitSumEtOffsetEndcap = cms.double(999999999),
140 hcalRecHitSumEtSlopeEndcap = cms.double(0.),
141 nTrackSolidConeBarrel =cms.double(999999999),
142 nTrackSolidConeEndcap =cms.double(999999999),
143 nTrackHollowConeBarrel =cms.double(999999999),
144 nTrackHollowConeEndcap =cms.double(999999999),
145 trackPtSumSolidConeBarrel =cms.double(999999999),
146 trackPtSumSolidConeEndcap =cms.double(999999999),
147 trackPtSumHollowConeBarrel =cms.double(999999999),
148 trackPtSumHollowConeEndcap =cms.double(999999999),
149 sigmaIetaIetaCutBarrel=cms.double(999999999),
150 sigmaIetaIetaCutEndcap=cms.double(999999999),
151 posCalcParameters = cms.PSet( T0_barl = cms.double(7.4),
152 T0_endc = cms.double(6.3),
153 T0_endcPresh = cms.double(3.6),
154 LogWeighted = cms.bool(
True),
155 W0 = cms.double(4.2),
156 X0 = cms.double(0.89)
159 RecHitFlagToBeExcludedEB = cleanedHybridSuperClusters.RecHitFlagToBeExcluded,
160 RecHitSeverityToBeExcludedEB = cleanedHybridSuperClusters.RecHitSeverityToBeExcluded,
161 RecHitFlagToBeExcludedEE = multi5x5BasicClustersCleaned.RecHitFlagToBeExcluded,
162 RecHitSeverityToBeExcludedEE = cleanedHybridSuperClusters.RecHitSeverityToBeExcluded,
165 from Configuration.ProcessModifiers.egamma_lowPt_exclusive_cff
import egamma_lowPt_exclusive
166 egamma_lowPt_exclusive.toModify(photons,
1.8.16