CMS 3D CMS Logo

HGCalRecHit_cfi.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
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 
61 dEdX_weights_v9 = cms.vdouble(0.0, # there is no layer zero
62  8.366557, # Mev
63  10.425456,
64  10.425456,
65  10.425456,
66  10.425456,
67  10.425456,
68  10.425456,
69  10.425456,
70  10.425456,
71  10.425456,
72  10.425456,
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  31.497849,
90  51.205434,
91  52.030486,
92  52.030486,
93  52.030486,
94  52.030486,
95  52.030486,
96  52.030486,
97  52.030486,
98  52.030486,
99  52.030486,
100  52.030486,
101  71.265149,
102  90.499812,
103  90.894274,
104  90.537470,
105  89.786205,
106  89.786205,
107  89.786205,
108  89.786205,
109  89.786205,
110  89.786205,
111  89.786205,
112  89.786205,
113  89.786205)
114 
115 
116 from Configuration.Eras.Modifier_phase2_hgcalV9_cff import phase2_hgcalV9
117 phase2_hgcalV9.toModify( dEdX, weights = dEdX_weights_v9 )
118 
119 dEdX_weights_v10 = cms.vdouble(0.0, # there is no layer zero
120  8.894541, # Mev
121  10.937907,
122  10.937907,
123  10.937907,
124  10.937907,
125  10.937907,
126  10.937907,
127  10.937907,
128  10.937907,
129  10.932882,
130  10.932882,
131  10.937907,
132  10.937907,
133  10.938169,
134  10.938169,
135  10.938169,
136  10.938169,
137  10.938169,
138  10.938169,
139  10.938169,
140  10.938169,
141  10.938169,
142  10.938169,
143  10.938169,
144  10.938169,
145  10.938169,
146  10.938169,
147  32.332097,
148  51.574301,
149  51.444192,
150  51.444192,
151  51.444192,
152  51.444192,
153  51.444192,
154  51.444192,
155  51.444192,
156  51.444192,
157  51.444192,
158  51.444192,
159  69.513118,
160  87.582044,
161  87.582044,
162  87.582044,
163  87.582044,
164  87.582044,
165  87.214571,
166  86.888309,
167  86.929520,
168  86.929520,
169  86.929520)
170 
171 
172 from Configuration.Eras.Modifier_phase2_hgcalV10_cff import phase2_hgcalV10
173 phase2_hgcalV10.toModify( dEdX, weights = dEdX_weights_v10 )
174 
175 # HGCAL rechit producer
176 HGCalRecHit = cms.EDProducer(
177  "HGCalRecHitProducer",
178  HGCEErechitCollection = cms.string('HGCEERecHits'),
179  HGCEEuncalibRecHitCollection = cms.InputTag('HGCalUncalibRecHit:HGCEEUncalibRecHits'),
180  HGCHEFrechitCollection = cms.string('HGCHEFRecHits'),
181  HGCHEFuncalibRecHitCollection = cms.InputTag('HGCalUncalibRecHit:HGCHEFUncalibRecHits'),
182  HGCHEBrechitCollection = cms.string('HGCHEBRecHits'),
183  HGCHEBuncalibRecHitCollection = cms.InputTag('HGCalUncalibRecHit:HGCHEBUncalibRecHits'),
184 
185  # digi constants
186  HGCEE_keV2DIGI = hgceeDigitizer.digiCfg.keV2fC,
187  HGCEE_isSiFE = HGCalUncalibRecHit.HGCEEConfig.isSiFE,
188  HGCEE_fCPerMIP = HGCalUncalibRecHit.HGCEEConfig.fCPerMIP,
189  HGCHEF_keV2DIGI = hgchefrontDigitizer.digiCfg.keV2fC,
190  HGCHEF_isSiFE = HGCalUncalibRecHit.HGCHEFConfig.isSiFE,
191  HGCHEF_fCPerMIP = HGCalUncalibRecHit.HGCHEFConfig.fCPerMIP,
192  HGCHEB_keV2DIGI = hgchebackDigitizer.digiCfg.keV2MIP,
193  HGCHEB_isSiFE = HGCalUncalibRecHit.HGCHEBConfig.isSiFE,
194  # don't produce rechit if detid is a ghost one
195  rangeMask = cms.uint32(4294442496),
196  rangeMatch = cms.uint32(1161838592),
197 
198 
199  # EM Scale calibrations
200  layerWeights = dEdX.weights,
201 
202  thicknessCorrection = cms.vdouble(1.132,1.092,1.084), # 100, 200, 300 um
203  HGCEE_noise_fC = hgceeDigitizer.digiCfg.noise_fC,
204  HGCEE_cce = hgceeDigitizer.digiCfg.chargeCollectionEfficiencies,
205  HGCHEF_noise_fC = hgchefrontDigitizer.digiCfg.noise_fC,
206  HGCHEF_cce = hgchefrontDigitizer.digiCfg.chargeCollectionEfficiencies,
207  HGCHEB_noise_MIP = hgchebackDigitizer.digiCfg.noise_MIP,
208  # algo
209  algo = cms.string("HGCalRecHitWorkerSimple")
210 
211  )
212 
213 phase2_hgcalV9.toModify( HGCalRecHit , thicknessCorrection = [0.759,0.760,0.773] ) #120um, 200um, 300um
214 phase2_hgcalV10.toModify( HGCalRecHit , thicknessCorrection = [0.781,0.775,0.769] ) #120um, 200um, 300um