CMS 3D CMS Logo
hgcalBackEndLayer2Producer_cfi.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
2 
3 from L1Trigger.L1THGCal.egammaIdentification import egamma_identification_drnn_cone, \
4  egamma_identification_drnn_dbscan, \
5  egamma_identification_histomax
6 
7 from Configuration.Eras.Modifier_phase2_hgcalV9_cff import phase2_hgcalV9
8 from Configuration.Eras.Modifier_phase2_hgcalV10_cff import phase2_hgcalV10
9 from Configuration.Eras.Modifier_phase2_hgcalV11_cff import phase2_hgcalV11
10 from Configuration.Eras.Modifier_phase2_hfnose_cff import phase2_hfnose
11 
12 
13 binSums = cms.vuint32(13, # 0
14  11, 11, 11, # 1 - 3
15  9, 9, 9, # 4 - 6
16  7, 7, 7, 7, 7, 7, # 7 - 12
17  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, # 13 - 27
18  3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3 # 28 - 41
19  )
20 
21 EE_DR_GROUP = 7
22 FH_DR_GROUP = 6
23 BH_DR_GROUP = 12
24 MAX_LAYERS = 52
25 
26 
27 dr_layerbylayer = ([0] + # no layer 0
28  [0.015]*EE_DR_GROUP + [0.020]*EE_DR_GROUP + [0.030]*EE_DR_GROUP + [0.040]*EE_DR_GROUP + # EM
29  [0.040]*FH_DR_GROUP + [0.050]*FH_DR_GROUP + # FH
30  [0.050]*BH_DR_GROUP) # BH
31 
32 
33 dr_layerbylayer_Bcoefficient = ([0] + # no layer 0
34  [0.020]*EE_DR_GROUP + [0.020]*EE_DR_GROUP + [0.02]*EE_DR_GROUP + [0.020]*EE_DR_GROUP + # EM
35  [0.020]*FH_DR_GROUP + [0.020]*FH_DR_GROUP + # FH
36  [0.020]*BH_DR_GROUP) # BH
37 
38 
39 neighbour_weights_1stOrder = cms.vdouble(0, 0.25, 0,
40  0.25, 0, 0.25,
41  0, 0.25, 0)
42 
43 neighbour_weights_2ndOrder = cms.vdouble(-0.25, 0.5, -0.25,
44  0.5, 0, 0.5,
45  -0.25, 0.5, -0.25)
46 
47 
48 seed_smoothing_ecal = cms.vdouble(
49  1., 1., 1.,
50  1., 1.1, 1.,
51  1., 1., 1.,
52  )
53 seed_smoothing_hcal = cms.vdouble(
54  1., 1., 1., 1., 1.,
55  1., 1., 1., 1., 1.,
56  1., 1., 2., 1., 1.,
57  1., 1., 1., 1., 1.,
58  1., 1., 1., 1., 1.,
59  )
60 
61 
62 distance_C3d_params = cms.PSet(type_multicluster=cms.string('dRC3d'),
63  dR_multicluster=cms.double(0.01),
64  minPt_multicluster=cms.double(0.5), # minimum pt of the multicluster (GeV)
65  dist_dbscan_multicluster=cms.double(0.),
66  minN_dbscan_multicluster=cms.uint32(0),
67  EGIdentification=egamma_identification_drnn_cone.clone(),
68  )
69 
70 
71 dbscan_C3d_params = cms.PSet(type_multicluster=cms.string('DBSCANC3d'),
72  dR_multicluster=cms.double(0.),
73  minPt_multicluster=cms.double(0.5), # minimum pt of the multicluster (GeV)
74  dist_dbscan_multicluster=cms.double(0.005),
75  minN_dbscan_multicluster=cms.uint32(3),
76  EGIdentification=egamma_identification_drnn_dbscan.clone())
77 
78 
79 histoMax_C3d_seeding_params = cms.PSet(type_histoalgo=cms.string('HistoMaxC3d'),
80  nBins_X1_histo_multicluster=cms.uint32(42), # bin size of about 0.012
81  nBins_X2_histo_multicluster=cms.uint32(216), # bin size of about 0.029
82  binSumsHisto=binSums,
83  kROverZMin=cms.double(0.076),
84  kROverZMax=cms.double(0.58),
85  threshold_histo_multicluster=cms.double(10.),
86  neighbour_weights=neighbour_weights_1stOrder,
87  seed_position=cms.string("TCWeighted"),#BinCentre, TCWeighted
88  seeding_space=cms.string("RPhi"),# RPhi, XY
89  seed_smoothing_ecal=seed_smoothing_ecal,
90  seed_smoothing_hcal=seed_smoothing_hcal,
91  )
92 
93 histoMax_C3d_clustering_params = cms.PSet(dR_multicluster=cms.double(0.03),
94  dR_multicluster_byLayer_coefficientA=cms.vdouble(),
95  dR_multicluster_byLayer_coefficientB=cms.vdouble(),
96  shape_threshold=cms.double(1.),
97  shape_distance=cms.double(0.015),
98  minPt_multicluster=cms.double(0.5), # minimum pt of the multicluster (GeV)
99  cluster_association=cms.string("NearestNeighbour"),
100  EGIdentification=egamma_identification_histomax.clone(),
101  )
102 
103 
104 # V9 samples have a different defintion of the dEdx calibrations. To account for it
105 # we reascale the thresholds of the clustering seeds
106 # (see https://indico.cern.ch/event/806845/contributions/3359859/attachments/1815187/2966402/19-03-20_EGPerf_HGCBE.pdf
107 # for more details)
108 phase2_hgcalV9.toModify(histoMax_C3d_seeding_params,
109  threshold_histo_multicluster=7.5, # MipT
110  )
111 
112 
113 histoMaxVariableDR_C3d_params = histoMax_C3d_clustering_params.clone(
114  dR_multicluster = cms.double(0.),
115  dR_multicluster_byLayer_coefficientA = cms.vdouble(dr_layerbylayer),
116  dR_multicluster_byLayer_coefficientB = cms.vdouble([0]*(MAX_LAYERS+1))
117  )
118 
119 
120 histoSecondaryMax_C3d_params = histoMax_C3d_seeding_params.clone(
121  type_histoalgo = cms.string('HistoSecondaryMaxC3d')
122  )
123 
124 histoMaxXYVariableDR_C3d_params = histoMax_C3d_seeding_params.clone(
125  seeding_space=cms.string("XY"),
126  nBins_X1_histo_multicluster=cms.uint32(192),
127  nBins_X2_histo_multicluster=cms.uint32(192)
128  )
129 
130 histoInterpolatedMax_C3d_params = histoMax_C3d_seeding_params.clone(
131  type_histoalgo = cms.string('HistoInterpolatedMaxC3d')
132  )
133 
134 
135 histoThreshold_C3d_params = histoMax_C3d_seeding_params.clone(
136  type_histoalgo = cms.string('HistoThresholdC3d')
137  )
138 
139 
140 histoMax_C3d_params = cms.PSet(
141  type_multicluster=cms.string('Histo'),
142  histoMax_C3d_clustering_parameters = histoMaxVariableDR_C3d_params.clone(),
143  histoMax_C3d_seeding_parameters = histoMax_C3d_seeding_params.clone(),
144  )
145 
146 
147 energy_interpretations_em = cms.PSet(type = cms.string('HGCalTriggerClusterInterpretationEM'),
148  layer_containment_corrs = cms.vdouble(0., 0.0, 1.38, 0.97, 1.11, 0.92, 1.06, 1.01, 1.06, 0.89, 1.0, 1.06, 0.89, 1.62, 1.83),
149  scale_correction_coeff = cms.vdouble(-27.15, 53.94),
150  dr_bylayer = cms.vdouble([0.015]*15)
151  )
152 
153 phase2_hgcalV10.toModify(
154  energy_interpretations_em,
155  layer_containment_corrs=cms.vdouble(0., 0.0, 1.73, 0.97, 1.08, 1.1, 1.01, 0.96, 1.18, 0.98, 1.05, 0.99, 0.89, 1.75, 2.0),
156  scale_correction_coeff=cms.vdouble(-27.53, 53.92),
157  )
158 
159 phase2_hgcalV11.toModify(
160  energy_interpretations_em,
161  layer_containment_corrs=cms.vdouble(0., 0.0, 1.28, 1.09, 1.0, 1.07, 1.09, 1.04, 1.0, 1.09, 1.07, 1.03, 0.93, 1.4, 1.89),
162  scale_correction_coeff=cms.vdouble(-24.96, 52.99),
163  )
164 
165 
166 energy_interpretations = cms.VPSet(energy_interpretations_em)
167 
168 be_proc = cms.PSet(ProcessorName = cms.string('HGCalBackendLayer2Processor3DClustering'),
169  C3d_parameters = histoMax_C3d_params.clone(),
170  energy_interpretations = energy_interpretations
171  )
172 
173 hgcalBackEndLayer2Producer = cms.EDProducer(
174  "HGCalBackendLayer2Producer",
175  InputCluster = cms.InputTag('hgcalBackEndLayer1Producer:HGCalBackendLayer1Processor2DClustering'),
176  ProcessorParameters = be_proc.clone()
177  )
178 
179 
180 hgcalBackEndLayer2ProducerHFNose = hgcalBackEndLayer2Producer.clone(
181  InputCluster = cms.InputTag('hgcalBackEndLayer1ProducerHFNose:HGCalBackendLayer1Processor2DClustering'),
182  ProcessorParameters = dict(
183  C3d_parameters = dict(
184  histoMax_C3d_seeding_parameters = dict(
185 
186  nBins_X1_histo_multicluster = 4, # R bin size: 5 FullModules * 8 TP
187  binSumsHisto = cms.vuint32(13,11,9,9),
188  kROverZMin = 0.025,
189  kROverZMax = 0.1
190  )
191  )
192  )
193 )
egammaIdentification
Generated for CMSSW Reference Manual by doxygen 1.8.16