1.62, 1.83),
scale_correction_coeff=cms.vdouble(53.94, -27.15),
dr_bylayer=cms.vdouble([0.015] * 15))
phase2_hgcalV10.toModify(
energy_interpretations_em,
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),
scale_correction_coeff=cms.vdouble(53.92, -27.53),
)
phase2_hgcalV11.toModify(
energy_interpretations_em,
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),
scale_correction_coeff=cms.vdouble(52.99, -24.96),
)
energy_interpretations = cms.VPSet(energy_interpretations_em)
be_proc = cms.PSet(
ProcessorName=cms.string('HGCalBackendLayer2Processor3DClustering'),
C3d_parameters=histoMax_C3d_params.clone(),
energy_interpretations=energy_interpretations)
hgcalBackEndLayer2Producer = cms.EDProducer(
"HGCalBackendLayer2Producer",
InputCluster=cms.InputTag(
'hgcalBackEndLayer1Producer:HGCalBackendLayer1Processor2DClustering'),
kind=cms.string('calo'),
))
phase2_hgcalV10.toModify(
pfClustersFromCombinedCaloHCal,
hadCorrector="L1Trigger/Phase2L1ParticleFlow/data/hadcorr_barrel_106X.root",
resol=cms.PSet(
etaBins=cms.vdouble(0.700, 1.200, 1.600),
offset=cms.vdouble(3.084, 2.715, 0.107),
scale=cms.vdouble(0.118, 0.130, 0.442),
kind=cms.string('calo'),
))
phase2_hgcalV11.toModify(
pfClustersFromCombinedCaloHCal,
hadCorrector="L1Trigger/Phase2L1ParticleFlow/data/hadcorr_barrel_110X.root",
resol=cms.PSet(
etaBins=cms.vdouble(0.700, 1.200, 1.600),
offset=cms.vdouble(2.909, 2.864, 0.294),
scale=cms.vdouble(0.119, 0.127, 0.442),
kind=cms.string('calo'),
))
pfTracksFromL1TracksBarrel = pfTracksFromL1Tracks.clone(
resolCalo=pfClustersFromCombinedCaloHCal.resol.clone(), )
pfClustersFromCombinedCaloHF = pfClustersFromCombinedCalo.clone(
ecalCandidates=[],
hcalHGCTowers=[],
phase2barrelCaloTowers=[],
hadCorrector=cms.string("L1Trigger/Phase2L1ParticleFlow/data/hfcorr.root"),
resol=cms.PSet(
etaBins=cms.vdouble(3.500, 4.000, 4.500, 5.000),
_correctedJets = cms.EDProducer(
"L1TCorrectedPFJetProducer",
jets=cms.InputTag("_tag_"),
correctorFile=cms.string(
"L1Trigger/Phase2L1ParticleFlow/data/jecs/jecs.PU200.root"),
correctorDir=cms.string("_dir_"),
copyDaughters=cms.bool(False))
# Using phase2_hgcalV10 to customize the config for all 106X samples, since there's no other modifier for it
from Configuration.Eras.Modifier_phase2_hgcalV10_cff import phase2_hgcalV10
phase2_hgcalV10.toModify(
_correctedJets,
correctorFile=
"L1Trigger/Phase2L1ParticleFlow/data/jecs/jecs.PU200_106X.root")
from Configuration.Eras.Modifier_phase2_hgcalV11_cff import phase2_hgcalV11
phase2_hgcalV11.toModify(
_correctedJets,
correctorFile=
"L1Trigger/Phase2L1ParticleFlow/data/jecs/jecs.PU200_110X.root")
ak4PFL1CaloCorrected = _correctedJets.clone(jets='ak4PFL1Calo',
correctorDir='L1CaloJets')
ak4PFL1PFCorrected = _correctedJets.clone(jets='ak4PFL1PF',
correctorDir='L1PFJets')
ak4PFL1PuppiCorrected = _correctedJets.clone(jets='ak4PFL1Puppi',
correctorDir='L1PuppiJets')
l1PFJets = cms.Sequence(ak4PFL1Calo + ak4PFL1PF + ak4PFL1Puppi +
ak4PFL1CaloCorrected + ak4PFL1PFCorrected +
ak4PFL1PuppiCorrected)
etaBins = cms.vdouble( 0.700, 1.200, 1.600),
offset = cms.vdouble( 0.873, 1.081, 1.563),
scale = cms.vdouble( 0.011, 0.015, 0.012),
kind = cms.string('calo'),
)
)
# use phase2_hgcalV10 to customize for 106X L1TDR MC even in the barrel, since there's no other modifier for it
from Configuration.Eras.Modifier_phase2_hgcalV10_cff import phase2_hgcalV10
from Configuration.Eras.Modifier_phase2_hgcalV11_cff import phase2_hgcalV11
phase2_hgcalV10.toModify(pfClustersFromL1EGClusters,
corrector = "", # In this setup, TP's are already calibrated correctly :-)
# L1Trigger/Phase2L1ParticleFlow/data/emcorr_barrel_106X.root",
resol = cms.PSet(
etaBins = cms.vdouble( 0.700, 1.200, 1.600),
offset = cms.vdouble( 0.946, 0.948, 1.171),
scale = cms.vdouble( 0.011, 0.018, 0.019),
kind = cms.string('calo')
)
)
phase2_hgcalV11.toModify(pfClustersFromL1EGClusters,
corrector = "", # In this setup, TP's are already calibrated correctly :-)
# L1Trigger/Phase2L1ParticleFlow/data/emcorr_barrel_110X.root",
resol = cms.PSet(
etaBins = cms.vdouble( 0.700, 1.200, 1.600),
offset = cms.vdouble( 0.838, 0.924, 1.101),
scale = cms.vdouble( 0.012, 0.017, 0.018),
kind = cms.string('calo')
)
)
]),
)
# Era Modification for HGCal v9: use correct training and WPs
phase2_hgcalV9.toModify(
egamma_identification_histomax,
Inputs=cms.vstring(input_features_histomax['v9_394']),
Weights=cms.vstring(bdt_weights_histomax['v9_394']),
WorkingPoints=cms.vdouble([
wps[eff]
for wps, eff in zip(working_points_histomax['v9_394'], tight_wp)
]))
phase2_hgcalV10.toModify(
egamma_identification_histomax,
Inputs=cms.vstring(input_features_histomax['v10_3151']),
Weights=cms.vstring(bdt_weights_histomax['v10_3151']),
WorkingPoints=cms.vdouble([
wps[eff]
for wps, eff in zip(working_points_histomax['v10_3151'], tight_wp)
]))
phase2_hgcalV11.toModify(
egamma_identification_histomax,
Inputs=cms.vstring(input_features_histomax['v10_3151']),
Weights=cms.vstring(bdt_weights_histomax['v10_3151']),
WorkingPoints=cms.vdouble([
wps[eff]
for wps, eff in zip(working_points_histomax['v10_3151'], tight_wp)
]))
resol=cms.PSet(etaBins=cms.vdouble(1.900, 2.200, 2.500, 2.800, 2.950),
offset=cms.vdouble(2.593, 3.089, 2.879, 2.664, 2.947),
scale=cms.vdouble(0.120, 0.098, 0.099, 0.098, 0.124),
kind=cms.string('calo')),
)
from Configuration.Eras.Modifier_phase2_hgcalV10_cff import phase2_hgcalV10
from Configuration.Eras.Modifier_phase2_hgcalV11_cff import phase2_hgcalV11
phase2_hgcalV10.toModify(
pfClustersFromHGC3DClusters,
corrector=
"L1Trigger/Phase2L1ParticleFlow/data/hadcorr_HGCal3D_TC_106X.root",
resol=cms.PSet(etaBins=cms.vdouble(1.700, 1.900, 2.200, 2.500, 2.800,
2.900),
offset=cms.vdouble(-0.819, 0.900, 2.032, 2.841, 2.865,
1.237),
scale=cms.vdouble(0.320, 0.225, 0.156, 0.108, 0.119, 0.338),
kind=cms.string('calo')),
)
phase2_hgcalV11.toModify(
pfClustersFromHGC3DClusters,
corrector=
"L1Trigger/Phase2L1ParticleFlow/data/hadcorr_HGCal3D_TC_110X.root",
resol=cms.PSet(
etaBins=cms.vdouble(1.700, 1.900, 2.200, 2.500, 2.800, 2.900),
offset=cms.vdouble(1.793, 1.827, 2.363, 2.538, 2.812, 2.642),
scale=cms.vdouble(0.138, 0.137, 0.124, 0.115, 0.106, 0.121),
kind=cms.string('calo'),
),
)
preEmId=cms.string(""),
resol=cms.PSet(
etaBins=cms.vdouble(1.900, 2.200, 2.500, 2.800, 2.950),
offset=cms.vdouble(0.566, 0.557, 0.456, 0.470, 0.324),
scale=cms.vdouble(0.030, 0.024, 0.024, 0.023, 0.042),
kind=cms.string('calo'),
))
from Configuration.Eras.Modifier_phase2_hgcalV10_cff import phase2_hgcalV10
from Configuration.Eras.Modifier_phase2_hgcalV11_cff import phase2_hgcalV11
phase2_hgcalV10.toModify(
pfClustersFromHGC3DClustersEM,
corrector="L1Trigger/Phase2L1ParticleFlow/data/emcorr_hgc_106X.root",
resol=cms.PSet(etaBins=cms.vdouble(1.700, 1.900, 2.200, 2.500, 2.800,
2.900),
offset=cms.vdouble(2.579, 2.176, 1.678, 0.911, 0.672,
-2.292),
scale=cms.vdouble(0.048, 0.026, 0.012, 0.016, 0.022, 0.538),
kind=cms.string('calo')),
)
phase2_hgcalV11.toModify(
pfClustersFromHGC3DClustersEM,
corrector="L1Trigger/Phase2L1ParticleFlow/data/emcorr_hgc_110X.root",
resol=cms.PSet(etaBins=cms.vdouble(1.700, 1.900, 2.200, 2.500, 2.800,
2.900),
offset=cms.vdouble(2.581, 2.289, 1.674, 0.927, 0.604,
-2.377),
scale=cms.vdouble(0.046, 0.025, 0.016, 0.017, 0.023, 0.500),
kind=cms.string('calo')),
)
etMin = cms.double(1.0),
resol = cms.PSet(
etaBins = cms.vdouble( 1.900, 2.200, 2.500, 2.800, 2.950),
offset = cms.vdouble( 2.593, 3.089, 2.879, 2.664, 2.947),
scale = cms.vdouble( 0.120, 0.098, 0.099, 0.098, 0.124),
kind = cms.string('calo')
),
)
from Configuration.Eras.Modifier_phase2_hgcalV10_cff import phase2_hgcalV10
from Configuration.Eras.Modifier_phase2_hgcalV11_cff import phase2_hgcalV11
phase2_hgcalV10.toModify(pfClustersFromHGC3DClusters,
corrector = "L1Trigger/Phase2L1ParticleFlow/data/hadcorr_HGCal3D_TC_106X.root",
resol = cms.PSet(
etaBins = cms.vdouble( 1.700, 1.900, 2.200, 2.500, 2.800, 2.900),
offset = cms.vdouble(-0.819, 0.900, 2.032, 2.841, 2.865, 1.237),
scale = cms.vdouble( 0.320, 0.225, 0.156, 0.108, 0.119, 0.338),
kind = cms.string('calo')
),
)
phase2_hgcalV11.toModify(pfClustersFromHGC3DClusters,
corrector = "L1Trigger/Phase2L1ParticleFlow/data/hadcorr_HGCal3D_TC_110X.root",
resol = cms.PSet(
etaBins = cms.vdouble( 1.700, 1.900, 2.200, 2.500, 2.800, 2.900),
offset = cms.vdouble( 1.413, 1.631, 2.366, 2.768, 2.908, 1.434),
scale = cms.vdouble( 0.169, 0.155, 0.132, 0.114, 0.137, 0.338),
kind = cms.string('calo')
),
)