CMS 3D CMS Logo
HGCalRecHit_cfi.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
2 from SimCalorimetry.HGCalSimProducers.hgcalDigitizer_cfi import *
3 from RecoLocalCalo.HGCalRecProducers.HGCalUncalibRecHit_cfi import *
4 
5 dEdX = cms.PSet(
6  weights = cms.vdouble(0.0, # there is no layer zero
7  8.603, # Mev
8  8.0675,
9  8.0675,
10  8.0675,
11  8.0675,
12  8.0675,
13  8.0675,
14  8.0675,
15  8.0675,
16  8.9515,
17  10.135,
18  10.135,
19  10.135,
20  10.135,
21  10.135,
22  10.135,
23  10.135,
24  10.135,
25  10.135,
26  11.682,
27  13.654,
28  13.654,
29  13.654,
30  13.654,
31  13.654,
32  13.654,
33  13.654,
34  38.2005,
35  55.0265,
36  49.871,
37  49.871,
38  49.871,
39  49.871,
40  49.871,
41  49.871,
42  49.871,
43  49.871,
44  49.871,
45  49.871,
46  62.005,
47  83.1675,
48  92.196,
49  92.196,
50  92.196,
51  92.196,
52  92.196,
53  92.196,
54  92.196,
55  92.196,
56  92.196,
57  92.196,
58  46.098),
59 
60  weightsNose = cms.vdouble(0.0, # there is no layer zero
61  39.500245, # Mev
62  39.756638,
63  39.756638,
64  39.756638,
65  39.756638,
66  66.020266,
67  92.283895,
68  92.283895)
69 )
70 
71 dEdX_weights_v9 = cms.vdouble(0.0, # there is no layer zero
72  8.366557, # Mev
73  10.425456,
74  10.425456,
75  10.425456,
76  10.425456,
77  10.425456,
78  10.425456,
79  10.425456,
80  10.425456,
81  10.425456,
82  10.425456,
83  10.425456,
84  10.425456,
85  10.425456,
86  10.425456,
87  10.425456,
88  10.425456,
89  10.425456,
90  10.425456,
91  10.425456,
92  10.425456,
93  10.425456,
94  10.425456,
95  10.425456,
96  10.425456,
97  10.425456,
98  10.425456,
99  31.497849,
100  51.205434,
101  52.030486,
102  52.030486,
103  52.030486,
104  52.030486,
105  52.030486,
106  52.030486,
107  52.030486,
108  52.030486,
109  52.030486,
110  52.030486,
111  71.265149,
112  90.499812,
113  90.894274,
114  90.537470,
115  89.786205,
116  89.786205,
117  89.786205,
118  89.786205,
119  89.786205,
120  89.786205,
121  89.786205,
122  89.786205,
123  89.786205)
124 
125 
126 from Configuration.Eras.Modifier_phase2_hgcalV9_cff import phase2_hgcalV9
127 phase2_hgcalV9.toModify( dEdX, weights = dEdX_weights_v9 )
128 
129 dEdX_weights_v10 = cms.vdouble(0.0, # there is no layer zero
130  8.894541, # Mev
131  10.937907,
132  10.937907,
133  10.937907,
134  10.937907,
135  10.937907,
136  10.937907,
137  10.937907,
138  10.937907,
139  10.932882,
140  10.932882,
141  10.937907,
142  10.937907,
143  10.938169,
144  10.938169,
145  10.938169,
146  10.938169,
147  10.938169,
148  10.938169,
149  10.938169,
150  10.938169,
151  10.938169,
152  10.938169,
153  10.938169,
154  10.938169,
155  10.938169,
156  10.938169,
157  32.332097,
158  51.574301,
159  51.444192,
160  51.444192,
161  51.444192,
162  51.444192,
163  51.444192,
164  51.444192,
165  51.444192,
166  51.444192,
167  51.444192,
168  51.444192,
169  69.513118,
170  87.582044,
171  87.582044,
172  87.582044,
173  87.582044,
174  87.582044,
175  87.214571,
176  86.888309,
177  86.929520,
178  86.929520,
179  86.929520)
180 
181 
182 from Configuration.Eras.Modifier_phase2_hgcalV10_cff import phase2_hgcalV10
183 phase2_hgcalV10.toModify( dEdX, weights = dEdX_weights_v10 )
184 
185 # HGCAL rechit producer
186 HGCalRecHit = cms.EDProducer(
187  "HGCalRecHitProducer",
188  HGCEErechitCollection = cms.string('HGCEERecHits'),
189  HGCEEuncalibRecHitCollection = cms.InputTag('HGCalUncalibRecHit:HGCEEUncalibRecHits'),
190  HGCHEFrechitCollection = cms.string('HGCHEFRecHits'),
191  HGCHEFuncalibRecHitCollection = cms.InputTag('HGCalUncalibRecHit:HGCHEFUncalibRecHits'),
192  HGCHEBrechitCollection = cms.string('HGCHEBRecHits'),
193  HGCHEBuncalibRecHitCollection = cms.InputTag('HGCalUncalibRecHit:HGCHEBUncalibRecHits'),
194  HGCHFNoserechitCollection = cms.string('HGCHFNoseRecHits'),
195  HGCHFNoseuncalibRecHitCollection = cms.InputTag('HGCalUncalibRecHit:HGCHFNoseUncalibRecHits'),
196 
197 
198  # digi constants
199  HGCEE_keV2DIGI = hgceeDigitizer.digiCfg.keV2fC,
200  HGCEE_isSiFE = HGCalUncalibRecHit.HGCEEConfig.isSiFE,
201  HGCEE_fCPerMIP = HGCalUncalibRecHit.HGCEEConfig.fCPerMIP,
202  HGCHEF_keV2DIGI = hgchefrontDigitizer.digiCfg.keV2fC,
203  HGCHEF_isSiFE = HGCalUncalibRecHit.HGCHEFConfig.isSiFE,
204  HGCHEF_fCPerMIP = HGCalUncalibRecHit.HGCHEFConfig.fCPerMIP,
205  HGCHEB_keV2DIGI = hgchebackDigitizer.digiCfg.keV2MIP,
206  HGCHEB_isSiFE = HGCalUncalibRecHit.HGCHEBConfig.isSiFE,
207  HGCHFNose_keV2DIGI = hfnoseDigitizer.digiCfg.keV2fC,
208  HGCHFNose_isSiFE = HGCalUncalibRecHit.HGCHFNoseConfig.isSiFE,
209  HGCHFNose_fCPerMIP = HGCalUncalibRecHit.HGCHFNoseConfig.fCPerMIP,
210  # don't produce rechit if detid is a ghost one
211  rangeMask = cms.uint32(4294442496),
212  rangeMatch = cms.uint32(1161838592),
213 
214 
215  # EM Scale calibrations
216  layerWeights = dEdX.weights,
217  layerNoseWeights = dEdX.weightsNose,
218 
219  thicknessCorrection = cms.vdouble(1.132,1.092,1.084), # 100, 200, 300 um
220  thicknessNoseCorrection = cms.vdouble(1.132,1.092,1.084), # 100, 200, 300 um
221  HGCEE_noise_fC = hgceeDigitizer.digiCfg.noise_fC,
222  HGCEE_cce = hgceeDigitizer.digiCfg.chargeCollectionEfficiencies,
223  HGCHEF_noise_fC = hgchefrontDigitizer.digiCfg.noise_fC,
224  HGCHEF_cce = hgchefrontDigitizer.digiCfg.chargeCollectionEfficiencies,
225  HGCHEB_noise_MIP = hgchebackDigitizer.digiCfg.noise,
226  HGCHFNose_noise_fC = hfnoseDigitizer.digiCfg.noise_fC,
227  HGCHFNose_cce = hfnoseDigitizer.digiCfg.chargeCollectionEfficiencies,
228 
229  # algo
230  algo = cms.string("HGCalRecHitWorkerSimple")
231 
232  )
233 
234 phase2_hgcalV9.toModify( HGCalRecHit , thicknessCorrection = [0.759,0.760,0.773] ) #120um, 200um, 300um
235 phase2_hgcalV10.toModify( HGCalRecHit , thicknessCorrection = [0.781,0.775,0.769] ) #120um, 200um, 300um
236 
237 phase2_hfnose.toModify( HGCalRecHit , thicknessNoseCorrection = [0.759,0.760,0.773])
hgcalDigitizer_cfi
HGCalUncalibRecHit_cfi
Generated for CMSSW Reference Manual by doxygen 1.8.11