neighbour_weights=neighbour_weights_1stOrder,
seed_position=cms.string("TCWeighted"), #BinCentre, TCWeighted
seeding_space=cms.string("RPhi"), # RPhi, XY
seed_smoothing_ecal=seed_smoothing_ecal,
seed_smoothing_hcal=seed_smoothing_hcal,
)
histoMax_C3d_clustering_params = cms.PSet(
dR_multicluster=cms.double(0.03),
dR_multicluster_byLayer_coefficientA=cms.vdouble(),
dR_multicluster_byLayer_coefficientB=cms.vdouble(),
shape_threshold=cms.double(1.),
shape_distance=cms.double(0.015),
minPt_multicluster=cms.double(0.5), # minimum pt of the multicluster (GeV)
cluster_association=cms.string("NearestNeighbour"),
EGIdentification=egamma_identification_histomax.clone(),
)
# V9 samples have a different definition of the dEdx calibrations. To account for it
# we rescale the thresholds of the clustering seeds
# (see https://indico.cern.ch/event/806845/contributions/3359859/attachments/1815187/2966402/19-03-20_EGPerf_HGCBE.pdf
# for more details)
phase2_hgcalV10.toModify(
histoMax_C3d_seeding_params,
threshold_histo_multicluster=7.5, # MipT
)
histoMaxVariableDR_C3d_params = histoMax_C3d_clustering_params.clone(
dR_multicluster=cms.double(0.),
dR_multicluster_byLayer_coefficientA=cms.vdouble(dr_layerbylayer),
dR_multicluster_byLayer_coefficientB=cms.vdouble([0] * (MAX_LAYERS + 1)))
dist_dbscan_multicluster=cms.double(0.005),
minN_dbscan_multicluster=cms.uint32(3),
EGIdentification=egamma_identification_drnn_dbscan.clone())
histoMax_C3d_params = cms.PSet(type_multicluster=cms.string('HistoMaxC3d'),
dR_multicluster=cms.double(0.03),
dR_multicluster_byLayer_coefficientA=cms.vdouble(),
dR_multicluster_byLayer_coefficientB=cms.vdouble(),
minPt_multicluster=cms.double(0.5), # minimum pt of the multicluster (GeV)
nBins_R_histo_multicluster=cms.uint32(36),
nBins_Phi_histo_multicluster=cms.uint32(216),
binSumsHisto=binSums,
threshold_histo_multicluster=cms.double(10.),
cluster_association=cms.string("NearestNeighbour"),
EGIdentification=egamma_identification_histomax.clone(),
neighbour_weights=neighbour_weights_1stOrder
)
# V9 samples have a different defintiion of the dEdx calibrations. To account for it
# we reascale the thresholds of the clustering seeds
# (see https://indico.cern.ch/event/806845/contributions/3359859/attachments/1815187/2966402/19-03-20_EGPerf_HGCBE.pdf
# for more details)
phase2_hgcalV9.toModify(histoMax_C3d_params,
threshold_histo_multicluster=7.5, # MipT
)
histoMaxVariableDR_C3d_params = histoMax_C3d_params.clone(
dR_multicluster = cms.double(0.),
dR_multicluster_byLayer_coefficientA = cms.vdouble(dr_layerbylayer),
dR_multicluster_byLayer_coefficientB = cms.vdouble([0]*(MAX_LAYERS+1))