CMS 3D CMS Logo

HGCalRecHit_cfi.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
4 
5 dEdX_weights = cms.vdouble(0.0, # there is no layer zero
6  8.603, # Mev
7  8.0675,
8  8.0675,
9  8.0675,
10  8.0675,
11  8.0675,
12  8.0675,
13  8.0675,
14  8.0675,
15  8.9515,
16  10.135,
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  11.682,
26  13.654,
27  13.654,
28  13.654,
29  13.654,
30  13.654,
31  13.654,
32  13.654,
33  38.2005,
34  55.0265,
35  49.871,
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  62.005,
46  83.1675,
47  92.196,
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  46.098)
58 
59 dEdX_weights_v9 = cms.vdouble(0.0, # there is no layer zero
60  8.366557, # Mev
61  10.425456,
62  10.425456,
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  31.497849,
88  51.205434,
89  52.030486,
90  52.030486,
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  71.265149,
100  90.499812,
101  90.894274,
102  90.537470,
103  89.786205,
104  89.786205,
105  89.786205,
106  89.786205,
107  89.786205,
108  89.786205,
109  89.786205,
110  89.786205,
111  89.786205)
112 
113 
114 # HGCAL rechit producer
115 HGCalRecHit = cms.EDProducer(
116  "HGCalRecHitProducer",
117  HGCEErechitCollection = cms.string('HGCEERecHits'),
118  HGCEEuncalibRecHitCollection = cms.InputTag('HGCalUncalibRecHit:HGCEEUncalibRecHits'),
119  HGCHEFrechitCollection = cms.string('HGCHEFRecHits'),
120  HGCHEFuncalibRecHitCollection = cms.InputTag('HGCalUncalibRecHit:HGCHEFUncalibRecHits'),
121  HGCHEBrechitCollection = cms.string('HGCHEBRecHits'),
122  HGCHEBuncalibRecHitCollection = cms.InputTag('HGCalUncalibRecHit:HGCHEBUncalibRecHits'),
123 
124  # digi constants
125  HGCEE_keV2DIGI = hgceeDigitizer.digiCfg.keV2fC,
126  HGCEE_isSiFE = HGCalUncalibRecHit.HGCEEConfig.isSiFE,
127  HGCEE_fCPerMIP = HGCalUncalibRecHit.HGCEEConfig.fCPerMIP,
128  HGCHEF_keV2DIGI = hgchefrontDigitizer.digiCfg.keV2fC,
129  HGCHEF_isSiFE = HGCalUncalibRecHit.HGCHEFConfig.isSiFE,
130  HGCHEF_fCPerMIP = HGCalUncalibRecHit.HGCHEFConfig.fCPerMIP,
131  HGCHEB_keV2DIGI = hgchebackDigitizer.digiCfg.keV2MIP,
132  HGCHEB_isSiFE = HGCalUncalibRecHit.HGCHEBConfig.isSiFE,
133  # don't produce rechit if detid is a ghost one
134  rangeMask = cms.uint32(4294442496),
135  rangeMatch = cms.uint32(1161838592),
136 
137 
138  # EM Scale calibrations
139  layerWeights = dEdX_weights,
140 
141  thicknessCorrection = cms.vdouble(1.132,1.092,1.084), # 100, 200, 300 um
142  HGCEE_noise_fC = hgceeDigitizer.digiCfg.noise_fC,
143  HGCEE_cce = hgceeDigitizer.digiCfg.chargeCollectionEfficiencies,
144  HGCHEF_noise_fC = hgchefrontDigitizer.digiCfg.noise_fC,
145  HGCHEF_cce = hgchefrontDigitizer.digiCfg.chargeCollectionEfficiencies,
146  HGCHEB_noise_MIP = hgchebackDigitizer.digiCfg.noise_MIP,
147  # algo
148  algo = cms.string("HGCalRecHitWorkerSimple")
149 
150  )
151 
152 from Configuration.Eras.Modifier_phase2_hgcalV9_cff import phase2_hgcalV9
153 phase2_hgcalV9.toModify( HGCalRecHit , layerWeights = dEdX_weights_v9 )
154 phase2_hgcalV9.toModify( HGCalRecHit , thicknessCorrection = cms.vdouble(0.759,0.760,0.773) ) #120um, 200um, 300um