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  #With the new regional em factors there are 7 different factors used.
220  #Six of them are for silicon and one for scint. For silicon it is in the following order
221  # CE_E_120um, CE_E_200um, CE_E_300um, CE_H_120um, CE_H_200um, CE_H_300um
222  thicknessCorrection = cms.vdouble(1.132,1.092,1.084, 1.0, 1.0, 1.0), # 100, 200, 300 um
223  deltasi_index_regemfac = cms.int32(3),
224  #One factor for scint
225  sciThicknessCorrection = cms.double(1.0),
226  thicknessNoseCorrection = cms.vdouble(1.132,1.092,1.084), # 100, 200, 300 um
227 
228  HGCEE_noise_fC = hgceeDigitizer.digiCfg.noise_fC,
229  HGCEE_cce = hgceeDigitizer.digiCfg.chargeCollectionEfficiencies,
230  HGCHEF_noise_fC = hgchefrontDigitizer.digiCfg.noise_fC,
231  HGCHEF_cce = hgchefrontDigitizer.digiCfg.chargeCollectionEfficiencies,
232  HGCHEB_noise_MIP = hgchebackDigitizer.digiCfg.noise,
233  HGCHFNose_noise_fC = hfnoseDigitizer.digiCfg.noise_fC,
234  HGCHFNose_cce = hfnoseDigitizer.digiCfg.chargeCollectionEfficiencies,
235 
236  # expected resolution on time for recHits - ns units
237  minValSiPar = cms.double(10.),
238  maxValSiPar = cms.double(1.e4),
239  noiseSiPar = cms.double(5.5),
240  constSiPar = cms.double(0.02),
241 
242  # algo
243  algo = cms.string("HGCalRecHitWorkerSimple")
244 
245  )
246 
247 # For silicon the order is: CE_E_120um, CE_E_200um, CE_E_300um, CE_H_120um, CE_H_200um, CE_H_300um
248 phase2_hgcalV9.toModify( HGCalRecHit , thicknessCorrection = [0.759,0.760,0.773, 1.0, 1.0, 1.0] )
249 phase2_hgcalV10.toModify( HGCalRecHit , thicknessCorrection = [0.77, 0.77, 0.77, 0.84, 0.84, 0.84] , sciThicknessCorrection = cms.double(0.90) )
250 
251 phase2_hfnose.toModify( HGCalRecHit , thicknessNoseCorrection = [0.759,0.760,0.773])
hgcalDigitizer_cfi
HGCalUncalibRecHit_cfi
Generated for CMSSW Reference Manual by doxygen 1.8.16