import FWCore.ParameterSet.Config as cms
#
# sequence to make photons from clusters in ECAL
#
# photon producer
from RecoEgamma.EgammaPhotonProducers.photonCore_cfi import *
from RecoEgamma.EgammaPhotonProducers.photons_cfi import *
photonSequence = cms.Sequence(photonCore + photons)
_photonSequenceFromMultiCl = photonSequence.copy()
_photonSequenceFromMultiCl += (photonCoreFromMultiCl + photonsFromMultiCl)
_photonSequenceWithIsland = photonSequence.copy()
_photonSequenceWithIsland += (islandPhotonCore + islandPhotons)
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toReplaceWith(photonSequence, _photonSequenceFromMultiCl)
from Configuration.Eras.Modifier_pA_2016_cff import pA_2016
from Configuration.Eras.Modifier_peripheralPbPb_cff import peripheralPbPb
from Configuration.Eras.Modifier_pp_on_XeXe_2017_cff import pp_on_XeXe_2017
for e in [pA_2016, peripheralPbPb, pp_on_XeXe_2017]:
e.toReplaceWith(photonSequence, _photonSequenceWithIsland)
zmumugammaOldAnalysis = DQMOffline.EGamma.zmumugammaAnalyzer_cfi.zmumugammaAnalysis.clone(
)
zmumugammaOldAnalysis.ComponentName = cms.string('zmumugammaOldAnalysis')
zmumugammaOldAnalysis.analyzerName = cms.string('zmumugammaOldValidation')
zmumugammaOldAnalysis.phoProducer = cms.InputTag('photons')
# HGCal customizations
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
stdPhotonAnalysisHGCal = stdPhotonAnalysis.clone()
stdPhotonAnalysisHGCal.ComponentName = 'stdPhotonAnalyzerHGCalFromMultiCl'
stdPhotonAnalysisHGCal.analyzerName = 'stdPhotonAnalyzerHGCalFromMultiCl'
stdPhotonAnalysisHGCal.phoProducer = 'photonsFromMultiCl'
stdPhotonAnalysisHGCal.isolationStrength = 2
stdPhotonAnalysisHGCal.etaMin = -3.0
stdPhotonAnalysisHGCal.etaMax = 3.0
stdPhotonAnalysisHGCal.maxPhoEta = 3.0
egammaDQMOffline = cms.Sequence(photonAnalysis * stdPhotonAnalysis *
zmumugammaOldAnalysis * zmumugammaAnalysis *
piZeroAnalysis * electronAnalyzerSequence)
_egammaDQMOfflineHGCal = egammaDQMOffline.copy()
_egammaDQMOfflineHGCal += stdPhotonAnalysisHGCal
phase2_hgcal.toReplaceWith(egammaDQMOffline, _egammaDQMOfflineHGCal)
from Configuration.Eras.Modifier_peripheralPbPb_cff import peripheralPbPb
from Configuration.ProcessModifiers.pp_on_AA_cff import pp_on_AA
from Configuration.Eras.Modifier_pp_on_XeXe_2017_cff import pp_on_XeXe_2017
for e in [peripheralPbPb, pp_on_AA, pp_on_XeXe_2017]:
e.toModify(stdPhotonAnalysis, phoProducer=cms.InputTag('islandPhotons'))
# calo geometry
#
# changed by tommaso. now the calibrations are read from Configuration/StaqndardSequences/data/*Conditions.cff
#
#HCAL reconstruction
from RecoLocalCalo.Configuration.hcalLocalReco_cff import *
from RecoLocalCalo.Configuration.hcalGlobalReco_cff import *
#
# sequence CaloLocalReco and CaloGlobalReco
#
calolocalrecoTask = cms.Task(ecalLocalRecoTask, hcalLocalRecoTask)
calolocalreco = cms.Sequence(calolocalrecoTask)
from RecoLocalCalo.HcalRecProducers.HcalHitSelection_cfi import *
reducedHcalRecHitsTask = cms.Task(reducedHcalRecHits)
reducedHcalRecHitsSequence = cms.Sequence(reducedHcalRecHitsTask)
#
# R.Ofierzynski (29.Oct.2009): add NZS sequence
#
from RecoLocalCalo.Configuration.hcalLocalRecoNZS_cff import *
calolocalrecoTaskNZS = cms.Task(ecalLocalRecoTask, hcalLocalRecoTask,
hcalLocalRecoTaskNZS)
calolocalrecoNZS = cms.Sequence(calolocalrecoTaskNZS)
from RecoLocalCalo.Configuration.hgcalLocalReco_cff import *
_phase2_calolocalrecoTask = calolocalrecoTask.copy()
_phase2_calolocalrecoTask.add(hgcalLocalRecoTask)
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toReplaceWith(calolocalrecoTask, _phase2_calolocalrecoTask)
electronSeedsTask = cms.Task(trackerDrivenElectronSeeds,
ecalDrivenElectronSeeds, electronMergedSeeds)
electronSeeds = cms.Sequence(electronSeedsTask)
_electronSeedsTaskFromMultiCl = electronSeedsTask.copy()
_electronSeedsTaskFromMultiCl.add(
cms.Task(ecalDrivenElectronSeedsFromMultiCl,
electronMergedSeedsFromMultiCl))
_electronSeedsFromMultiCl = cms.Sequence(_electronSeedsTaskFromMultiCl)
from Configuration.ProcessModifiers.pp_on_AA_cff import pp_on_AA
pp_on_AA.toReplaceWith(
electronSeedsTask,
electronSeedsTask.copyAndExclude([trackerDrivenElectronSeeds]))
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toReplaceWith(electronSeedsTask, _electronSeedsTaskFromMultiCl)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
_fastSim_electronSeedsTask = electronSeedsTask.copy()
_fastSim_electronSeedsTask.replace(
trackerDrivenElectronSeeds,
cms.Task(trackerDrivenElectronSeedsTmp, trackerDrivenElectronSeeds))
fastSim.toReplaceWith(electronSeedsTask, _fastSim_electronSeedsTask)
# replace the ECAL driven electron track candidates with the FastSim emulated ones
import FastSimulation.Tracking.electronCkfTrackCandidates_cff
fastElectronCkfTrackCandidates = FastSimulation.Tracking.electronCkfTrackCandidates_cff.electronCkfTrackCandidates.clone(
)
from TrackingTools.GsfTracking.CkfElectronCandidateMaker_cff import *
from TrackingTools.GsfTracking.GsfElectronGsfFit_cff import *
electronGsfTrackingTask = cms.Task(electronSeedsTask,
allMetFilterPaths=['HBHENoiseFilter','HBHENoiseIsoFilter','CSCTightHaloFilter','CSCTightHaloTrkMuUnvetoFilter','CSCTightHalo2015Filter','globalTightHalo2016Filter','globalSuperTightHalo2016Filter','HcalStripHaloFilter','hcalLaserEventFilter','EcalDeadCellTriggerPrimitiveFilter','EcalDeadCellBoundaryEnergyFilter','goodVertices','eeBadScFilter',
'ecalLaserCorrFilter','trkPOGFilters','chargedHadronTrackResolutionFilter','muonBadTrackFilter',
'BadChargedCandidateFilter','BadPFMuonFilter','BadChargedCandidateSummer16Filter','BadPFMuonSummer16Filter',
'trkPOG_manystripclus53X','trkPOG_toomanystripclus53X','trkPOG_logErrorTooManyClusters','METFilters']
def miniAOD_customizeMETFiltersFastSim(process):
"""Replace some MET filters that don't work in FastSim with trivial bools"""
for X in 'CSCTightHaloFilter', 'CSCTightHaloTrkMuUnvetoFilter','CSCTightHalo2015Filter','globalTightHalo2016Filter','globalSuperTightHalo2016Filter','HcalStripHaloFilter':
process.globalReplace(X, cms.EDFilter("HLTBool", result=cms.bool(True)))
for X in 'manystripclus53X', 'toomanystripclus53X', 'logErrorTooManyClusters':
process.globalReplace(X, cms.EDFilter("HLTBool", result=cms.bool(False)))
return process
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toReplaceWith( Flag_HBHENoiseFilter, cms.Path() )
phase2_hgcal.toReplaceWith( Flag_HBHENoiseIsoFilter, cms.Path() )
phase2_hgcal.toReplaceWith( Flag_eeBadScFilter, cms.Path() )
metFilterPathsTask = cms.Task(
HBHENoiseFilterResultProducer,
HBHENoiseFilter,
HBHENoiseIsoFilter,
hcalLaserEventFilter,
EcalDeadCellTriggerPrimitiveFilter,
eeBadScFilter,
ecalLaserCorrFilter,
EcalDeadCellBoundaryEnergyFilter,
primaryVertexFilter,
CSCTightHaloFilter,
CSCTightHaloTrkMuUnvetoFilter,
particleFlowCluster = cms.Sequence(
pfClusteringPS*
pfClusteringECAL*
pfClusteringHBHEHF*
pfClusteringHO
)
#HGCal
from RecoParticleFlow.PFClusterProducer.particleFlowRecHitHGC_cfi import particleFlowRecHitHGC
pfClusteringHGCal = cms.Sequence(particleFlowRecHitHGC)
_phase2_hgcal_particleFlowCluster = particleFlowCluster.copy()
_phase2_hgcal_particleFlowCluster += pfClusteringHGCal
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toReplaceWith( particleFlowCluster, _phase2_hgcal_particleFlowCluster )
#timing
from RecoParticleFlow.PFClusterProducer.particleFlowClusterTimeAssigner_cfi import particleFlowTimeAssignerECAL
from RecoParticleFlow.PFSimProducer.ecalBarrelClusterFastTimer_cfi import ecalBarrelClusterFastTimer
_phase2_timing_particleFlowClusterECALSequence = cms.Sequence(ecalBarrelClusterFastTimer*
particleFlowTimeAssignerECAL*
particleFlowClusterECALSequence.copy())
from Configuration.Eras.Modifier_phase2_timing_cff import phase2_timing
phase2_timing.toReplaceWith(particleFlowClusterECALSequence,
_phase2_timing_particleFlowClusterECALSequence)
phase2_timing.toModify(particleFlowClusterECAL,
inputECAL = cms.InputTag('particleFlowTimeAssignerECAL'))
# GEM settings _gem_RawToDigiTask = RawToDigiTask.copy() _gem_RawToDigiTask.add(muonGEMDigis) from Configuration.Eras.Modifier_run2_GEM_2017_cff import run2_GEM_2017 run2_GEM_2017.toReplaceWith(RawToDigiTask, _gem_RawToDigiTask) from Configuration.Eras.Modifier_run3_GEM_cff import run3_GEM run3_GEM.toReplaceWith(RawToDigiTask, _gem_RawToDigiTask) from EventFilter.HGCalRawToDigi.HGCalRawToDigi_cfi import * _hgcal_RawToDigiTask = RawToDigiTask.copy() _hgcal_RawToDigiTask.add(hgcalDigis) from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal phase2_hgcal.toReplaceWith(RawToDigiTask,_hgcal_RawToDigiTask) # RPC New Readout Validation from Configuration.Eras.Modifier_stage2L1Trigger_2017_cff import stage2L1Trigger_2017 _rpc_NewReadoutVal_RawToDigiTask = RawToDigiTask.copy() _rpc_NewReadoutVal_RawToDigiTask_noTk = RawToDigiTask_noTk.copy() _rpc_NewReadoutVal_RawToDigiTask.add(muonRPCNewDigis) _rpc_NewReadoutVal_RawToDigiTask_noTk.add(muonRPCNewDigis) stage2L1Trigger_2017.toReplaceWith(RawToDigiTask, _rpc_NewReadoutVal_RawToDigiTask) stage2L1Trigger_2017.toReplaceWith(RawToDigiTask_noTk, _rpc_NewReadoutVal_RawToDigiTask) from Configuration.Eras.Modifier_fastSim_cff import fastSim fastSim.toReplaceWith(RawToDigiTask, RawToDigiTask.copyAndExclude([muonRPCNewDigis])) fastSim.toReplaceWith(RawToDigiTask_noTk, RawToDigiTask_noTk.copyAndExclude([muonRPCNewDigis])) _hfnose_RawToDigiTask = RawToDigiTask.copy()
import FWCore.ParameterSet.Config as cms
from RecoParticleFlow.PFTracking.trackerDrivenElectronSeeds_cff import *
from RecoEgamma.EgammaElectronProducers.ecalDrivenElectronSeeds_cfi import *
from RecoParticleFlow.PFTracking.mergedElectronSeeds_cfi import *
electronSeeds = cms.Sequence(trackerDrivenElectronSeeds *
ecalDrivenElectronSeeds * electronMergedSeeds)
_electronSeedsFromMultiCl = electronSeeds.copy()
_electronSeedsFromMultiCl += cms.Sequence(ecalDrivenElectronSeedsFromMultiCl *
electronMergedSeedsFromMultiCl)
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toReplaceWith(electronSeeds, _electronSeedsFromMultiCl)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
_fastSim_electronSeeds = electronSeeds.copy()
_fastSim_electronSeeds.replace(
trackerDrivenElectronSeeds,
trackerDrivenElectronSeedsTmp + trackerDrivenElectronSeeds)
fastSim.toReplaceWith(electronSeeds, _fastSim_electronSeeds)
# replace the ECAL driven electron track candidates with the FastSim emulated ones
import FastSimulation.Tracking.electronCkfTrackCandidates_cff
fastElectronCkfTrackCandidates = FastSimulation.Tracking.electronCkfTrackCandidates_cff.electronCkfTrackCandidates.clone(
)
from TrackingTools.GsfTracking.CkfElectronCandidateMaker_cff import *
from TrackingTools.GsfTracking.GsfElectronGsfFit_cff import *
dqmElectronClientSelectionEtIso = dqmElectronOfflineClient.clone()
dqmElectronClientSelectionEtIso.InputFolderName = cms.string(
"Egamma/Electrons/Ele4_Et10TkIso1")
dqmElectronClientSelectionEtIso.OutputFolderName = cms.string(
"Egamma/Electrons/Ele4_Et10TkIso1")
#dqmElectronClientSelectionEtIsoElID = dqmElectronOfflineClient.clone() ;
#dqmElectronClientSelectionEtIsoElID.InputFolderName = cms.string("Egamma/Electrons/Ele4_Et10TkIso1ElID") ;
#dqmElectronClientSelectionEtIsoElID.OutputFolderName = cms.string("Egamma/Electrons/Ele4_Et10TkIso1ElID") ;
dqmElectronClientTagAndProbe = dqmElectronOfflineClient.clone()
dqmElectronClientTagAndProbe.InputFolderName = cms.string(
"Egamma/Electrons/Ele5_TagAndProbe")
dqmElectronClientTagAndProbe.OutputFolderName = cms.string(
"Egamma/Electrons/Ele5_TagAndProbe")
dqmElectronClientTagAndProbe.EffHistoTitle = cms.string("")
electronOfflineClientSequence = cms.Sequence(
dqmElectronClientAllElectrons * dqmElectronClientSelectionEt *
dqmElectronClientSelectionEtIso
# * dqmElectronClientSelectionEtIsoElID
* dqmElectronClientTagAndProbe)
_electronOfflineClientSequenceFromMultiCl = electronOfflineClientSequence.copy(
)
_electronOfflineClientSequenceFromMultiCl += dqmElectronClientAllElectronsFromMultiCl
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toReplaceWith(electronOfflineClientSequence,
_electronOfflineClientSequenceFromMultiCl)
gsfElectronTask = cms.Task(ecalDrivenGsfElectronCores, ecalDrivenGsfElectrons,
gsfElectronCores, gsfElectrons)
gsfElectronSequence = cms.Sequence(gsfElectronTask)
gsfEcalDrivenElectronTask = cms.Task(ecalDrivenGsfElectronCores,
ecalDrivenGsfElectrons)
gsfEcalDrivenElectronSequence = cms.Sequence(gsfEcalDrivenElectronTask)
_gsfEcalDrivenElectronTaskFromMultiCl = gsfEcalDrivenElectronTask.copy()
_gsfEcalDrivenElectronTaskFromMultiCl.add(
cms.Task(ecalDrivenGsfElectronCoresFromMultiCl,
ecalDrivenGsfElectronsFromMultiCl))
_gsfEcalDrivenElectronSequenceFromMultiCl = cms.Sequence(
_gsfEcalDrivenElectronTaskFromMultiCl)
#gsfElectronMergingSequence = cms.Sequence(gsfElectronCores*gsfElectrons)
from RecoEgamma.EgammaElectronProducers.edBasedElectronIso_cff import *
from RecoEgamma.EgammaElectronProducers.pfBasedElectronIso_cff import *
electronIsoTask = cms.Task(edBasedElectronIsoTask, pfBasedElectronIsoTask)
electronIsoSequence = cms.Sequence(electronIsoTask)
gsfElectronMergingTask = cms.Task(electronIsoTask, gsfElectronCores,
gsfElectrons)
gsfElectronMergingSequence = cms.Sequence(gsfElectronMergingTask)
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toReplaceWith(gsfEcalDrivenElectronTask,
_gsfEcalDrivenElectronTaskFromMultiCl)
import FWCore.ParameterSet.Config as cms
#==============================================================================
# Sequence to make final electrons.
# In the past, this was including the seeding, but this one is directly
# imported in the reco sequences since the integration with pflow.
#==============================================================================
from RecoEgamma.EgammaElectronProducers.gsfElectronModules_cff import *
gsfElectronSequence = cms.Sequence(ecalDrivenGsfElectronCores*ecalDrivenGsfElectrons*gsfElectronCores*gsfElectrons)
gsfEcalDrivenElectronSequence = cms.Sequence(ecalDrivenGsfElectronCores*ecalDrivenGsfElectrons)
_gsfEcalDrivenElectronSequenceFromMultiCl = gsfEcalDrivenElectronSequence.copy()
_gsfEcalDrivenElectronSequenceFromMultiCl += cms.Sequence(ecalDrivenGsfElectronCoresFromMultiCl*ecalDrivenGsfElectronsFromMultiCl)
#gsfElectronMergingSequence = cms.Sequence(gsfElectronCores*gsfElectrons)
from RecoEgamma.EgammaElectronProducers.edBasedElectronIso_cff import *
from RecoEgamma.EgammaElectronProducers.pfBasedElectronIso_cff import *
electronIsoSequence = cms.Sequence(
edBasedElectronIsoSequence+
pfBasedElectronIsoSequence
)
gsfElectronMergingSequence = cms.Sequence(electronIsoSequence*gsfElectronCores*gsfElectrons)
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toReplaceWith(
gsfEcalDrivenElectronSequence, _gsfEcalDrivenElectronSequenceFromMultiCl
)
#from RecoEcal.EgammaClusterProducers.correctedDynamicHybridSuperClusters_cfi import *
# PFECAL super clusters, either hybrid-clustering clone (Box) or mustache.
particleFlowSuperClusteringSequence = cms.Sequence(
particleFlowSuperClusterECAL)
particleFlowSuperClusterHGCal = particleFlowSuperClusterECAL.clone()
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toModify(
particleFlowSuperClusterHGCal,
PFClusters=cms.InputTag('particleFlowClusterHGCal'),
useRegression=cms.bool(False), #no HGCal regression yet
use_preshower=cms.bool(False),
PFBasicClusterCollectionEndcap=cms.string(""),
PFSuperClusterCollectionEndcap=cms.string(""),
PFSuperClusterCollectionEndcapWithPreshower=cms.string(""),
thresh_PFClusterEndcap=cms.double(1.5e-1) # 150 MeV threshold
)
particleFlowSuperClusterHGCalFromMultiCl = particleFlowSuperClusterHGCal.clone(
)
phase2_hgcal.toModify(
particleFlowSuperClusterHGCalFromMultiCl,
PFClusters=cms.InputTag('particleFlowClusterHGCalFromMultiCl'))
_phase2_hgcal_particleFlowSuperClusteringSequence = particleFlowSuperClusteringSequence.copy(
)
_phase2_hgcal_particleFlowSuperClusteringSequence += particleFlowSuperClusterHGCal
_phase2_hgcal_particleFlowSuperClusteringSequence += particleFlowSuperClusterHGCalFromMultiCl
phase2_hgcal.toReplaceWith(particleFlowSuperClusteringSequence,
_phase2_hgcal_particleFlowSuperClusteringSequence)
dqmElectronClientSelectionEt.OutputFolderName = cms.string("Egamma/Electrons/Ele3_Et10") ;
dqmElectronClientSelectionEtIso = dqmElectronOfflineClient.clone() ;
dqmElectronClientSelectionEtIso.InputFolderName = cms.string("Egamma/Electrons/Ele4_Et10TkIso1") ;
dqmElectronClientSelectionEtIso.OutputFolderName = cms.string("Egamma/Electrons/Ele4_Et10TkIso1") ;
#dqmElectronClientSelectionEtIsoElID = dqmElectronOfflineClient.clone() ;
#dqmElectronClientSelectionEtIsoElID.InputFolderName = cms.string("Egamma/Electrons/Ele4_Et10TkIso1ElID") ;
#dqmElectronClientSelectionEtIsoElID.OutputFolderName = cms.string("Egamma/Electrons/Ele4_Et10TkIso1ElID") ;
dqmElectronClientTagAndProbe = dqmElectronOfflineClient.clone() ;
dqmElectronClientTagAndProbe.InputFolderName = cms.string("Egamma/Electrons/Ele5_TagAndProbe") ;
dqmElectronClientTagAndProbe.OutputFolderName = cms.string("Egamma/Electrons/Ele5_TagAndProbe") ;
dqmElectronClientTagAndProbe.EffHistoTitle = cms.string("")
electronOfflineClientSequence = cms.Sequence(
dqmElectronClientAllElectrons
* dqmElectronClientSelectionEt
* dqmElectronClientSelectionEtIso
# * dqmElectronClientSelectionEtIsoElID
* dqmElectronClientTagAndProbe
)
_electronOfflineClientSequenceFromMultiCl = electronOfflineClientSequence.copy()
_electronOfflineClientSequenceFromMultiCl += dqmElectronClientAllElectronsFromMultiCl
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toReplaceWith(
electronOfflineClientSequence, _electronOfflineClientSequenceFromMultiCl
)
import FWCore.ParameterSet.Config as cms # # sequence to make photons from clusters in ECAL # # photon producer from RecoEgamma.EgammaPhotonProducers.photonCore_cfi import * from RecoEgamma.EgammaPhotonProducers.photons_cfi import * photonSequence = cms.Sequence( photonCore + photons ) _photonSequenceFromMultiCl = photonSequence.copy() _photonSequenceFromMultiCl += ( photonCoreFromMultiCl + photonsFromMultiCl) _photonSequenceWithIsland = photonSequence.copy() _photonSequenceWithIsland += ( photonCoreIsland + photonsIsland) from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal phase2_hgcal.toReplaceWith( photonSequence, _photonSequenceFromMultiCl ) from Configuration.Eras.Modifier_pp_on_XeXe_2017_cff import pp_on_XeXe_2017 pp_on_XeXe_2017.toReplaceWith( photonSequence, _photonSequenceWithIsland )
import FWCore.ParameterSet.Config as cms #============================================================================== # Sequence to make final electrons. # In the past, this was including the seeding, but this one is directly # imported in the reco sequences since the integration with pflow. #============================================================================== from RecoEgamma.EgammaElectronProducers.ecalDrivenGsfElectronCores_cfi import ecalDrivenGsfElectronCores from RecoEgamma.EgammaElectronProducers.ecalDrivenGsfElectronCoresHGC_cff import ecalDrivenGsfElectronCoresHGC from RecoEgamma.EgammaElectronProducers.gsfElectrons_cfi import * gsfEcalDrivenElectronTask = cms.Task(ecalDrivenGsfElectronCores,ecalDrivenGsfElectrons) gsfEcalDrivenElectronSequence = cms.Sequence(gsfEcalDrivenElectronTask) _gsfEcalDrivenElectronTaskHGC = gsfEcalDrivenElectronTask.copy() _gsfEcalDrivenElectronTaskHGC.add(cms.Task(ecalDrivenGsfElectronCoresHGC,ecalDrivenGsfElectronsHGC)) from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal phase2_hgcal.toReplaceWith( gsfEcalDrivenElectronTask, _gsfEcalDrivenElectronTaskHGC )
import FWCore.ParameterSet.Config as cms from RecoParticleFlow.PFTracking.trackerDrivenElectronSeeds_cff import * from RecoEgamma.EgammaElectronProducers.ecalDrivenElectronSeeds_cfi import * from RecoParticleFlow.PFTracking.mergedElectronSeeds_cfi import * electronSeeds = cms.Sequence(trackerDrivenElectronSeeds*ecalDrivenElectronSeeds*electronMergedSeeds) _electronSeedsFromMultiCl = electronSeeds.copy() _electronSeedsFromMultiCl += cms.Sequence(ecalDrivenElectronSeedsFromMultiCl*electronMergedSeedsFromMultiCl) from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal phase2_hgcal.toReplaceWith( electronSeeds, _electronSeedsFromMultiCl ) from TrackingTools.GsfTracking.CkfElectronCandidateMaker_cff import * from TrackingTools.GsfTracking.GsfElectronGsfFit_cff import * electronGsfTracking = cms.Sequence(electronSeeds*electronCkfTrackCandidates*electronGsfTracks) _electronGsfTracking = electronGsfTracking.copy() _electronGsfTracking += cms.Sequence(electronCkfTrackCandidatesFromMultiCl*electronGsfTracksFromMultiCl) from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal phase2_hgcal.toReplaceWith( electronGsfTracking, _electronGsfTracking )
from RecoParticleFlow.PFTracking.trackerDrivenElectronSeeds_cff import * from RecoEgamma.EgammaElectronProducers.ecalDrivenElectronSeeds_cfi import * from RecoParticleFlow.PFTracking.mergedElectronSeeds_cfi import * electronSeedsTask = cms.Task(trackerDrivenElectronSeeds,ecalDrivenElectronSeeds,electronMergedSeeds) electronSeeds = cms.Sequence(electronSeedsTask) _electronSeedsTaskFromMultiCl = electronSeedsTask.copy() _electronSeedsTaskFromMultiCl.add(cms.Task(ecalDrivenElectronSeedsFromMultiCl,electronMergedSeedsFromMultiCl)) _electronSeedsFromMultiCl = cms.Sequence(_electronSeedsTaskFromMultiCl) from Configuration.Eras.Modifier_pp_on_AA_2018_cff import pp_on_AA_2018 pp_on_AA_2018.toReplaceWith(electronSeedsTask, electronSeedsTask.copyAndExclude([trackerDrivenElectronSeeds])) from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal phase2_hgcal.toReplaceWith( electronSeedsTask, _electronSeedsTaskFromMultiCl ) from Configuration.Eras.Modifier_fastSim_cff import fastSim _fastSim_electronSeedsTask = electronSeedsTask.copy() _fastSim_electronSeedsTask.replace(trackerDrivenElectronSeeds, cms.Task(trackerDrivenElectronSeedsTmp,trackerDrivenElectronSeeds)) fastSim.toReplaceWith(electronSeedsTask, _fastSim_electronSeedsTask) # replace the ECAL driven electron track candidates with the FastSim emulated ones import FastSimulation.Tracking.electronCkfTrackCandidates_cff fastElectronCkfTrackCandidates = FastSimulation.Tracking.electronCkfTrackCandidates_cff.electronCkfTrackCandidates.clone() from TrackingTools.GsfTracking.CkfElectronCandidateMaker_cff import * from TrackingTools.GsfTracking.GsfElectronGsfFit_cff import * electronGsfTrackingTask = cms.Task(electronSeedsTask,electronCkfTrackCandidates,electronGsfTracks) electronGsfTracking = cms.Sequence(electronGsfTrackingTask) _electronGsfTrackingTask = electronGsfTrackingTask.copy()
_fastSim_electronSeedsTask = electronSeedsTask.copy() _fastSim_electronSeedsTask.replace(trackerDrivenElectronSeeds, cms.Task(trackerDrivenElectronSeedsTmp,trackerDrivenElectronSeeds)) fastSim.toReplaceWith(electronSeedsTask, _fastSim_electronSeedsTask) # replace the ECAL driven electron track candidates with the FastSim emulated ones import FastSimulation.Tracking.electronCkfTrackCandidates_cff fastElectronCkfTrackCandidates = FastSimulation.Tracking.electronCkfTrackCandidates_cff.electronCkfTrackCandidates.clone() from TrackingTools.GsfTracking.CkfElectronCandidateMaker_cff import * from TrackingTools.GsfTracking.GsfElectronGsfFit_cff import * electronGsfTrackingTask = cms.Task(electronSeedsTask,electronCkfTrackCandidates,electronGsfTracks) electronGsfTracking = cms.Sequence(electronGsfTrackingTask) _electronGsfTrackingTask = electronGsfTrackingTask.copy() _fastSim_electronGsfTrackingTask = electronGsfTrackingTask.copy() _fastSim_electronGsfTrackingTask.replace(electronCkfTrackCandidates,fastElectronCkfTrackCandidates) fastSim.toReplaceWith(electronGsfTrackingTask,_fastSim_electronGsfTrackingTask) from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal phase2_hgcal.toReplaceWith( electronGsfTrackingTask, _electronGsfTrackingTask ) from SimTracker.TrackAssociation.trackTimeValueMapProducer_cfi import trackTimeValueMapProducer gsfTrackTimeValueMapProducer = trackTimeValueMapProducer.clone(trackSrc = 'electronGsfTracks') electronGsfTrackingWithTimingTask = cms.Task(electronGsfTrackingTask.copy(),gsfTrackTimeValueMapProducer) electronGsfTrackingWithTiming = cms.Sequence(electronGsfTrackingWithTimingTask) from Configuration.Eras.Modifier_phase2_timing_cff import phase2_timing phase2_timing.toReplaceWith(electronGsfTrackingTask, electronGsfTrackingWithTimingTask)
mergedSuperClustersFromMultiCl = mergedSuperClusters.clone()
mergedSuperClustersFromMultiCl.src = cms.VInputTag(
cms.InputTag("particleFlowSuperClusterECAL",
"particleFlowSuperClusterECALBarrel"),
cms.InputTag("particleFlowSuperClusterHGCalFromMultiCl", ""))
dqmElectronAnalysisAllElectronsFromMultiCl = dqmElectronAnalysisAllElectrons.clone(
)
dqmElectronAnalysisAllElectronsFromMultiCl.OutputFolderName = 'Egamma/Electrons/Ele2FromMultiCl_All'
dqmElectronAnalysisAllElectronsFromMultiCl.MaxAbsEtaMatchingObject = 3.0
dqmElectronAnalysisAllElectronsFromMultiCl.EtaMax = 3.0
dqmElectronAnalysisAllElectronsFromMultiCl.EtaMin = -3.0
dqmElectronAnalysisAllElectronsFromMultiCl.MaxAbsEta = 3.0
dqmElectronAnalysisAllElectronsFromMultiCl.ElectronCollection = 'ecalDrivenGsfElectronsFromMultiCl'
dqmElectronAnalysisAllElectronsFromMultiCl.MatchingObjectCollection = 'mergedSuperClustersFromMultiCl'
_electronAnalyzerSequenceFromMultiCl = electronAnalyzerSequence.copy()
_electronAnalyzerSequenceFromMultiCl += cms.Sequence(
mergedSuperClustersFromMultiCl +
dqmElectronAnalysisAllElectronsFromMultiCl)
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toModify(mergedSuperClusters,
src=cms.VInputTag(
cms.InputTag("particleFlowSuperClusterECAL",
"particleFlowSuperClusterECALBarrel"),
cms.InputTag("particleFlowSuperClusterHGCal", "")))
phase2_hgcal.toReplaceWith(electronAnalyzerSequence,
_electronAnalyzerSequenceFromMultiCl)
_phase2_hgcal_particleFlowTmp = cms.EDProducer("PFCandidateListMerger",
src=cms.VInputTag(
"particleFlowTmpBarrel",
"simPFProducer"))
_phase2_hgcal_simPFTask = cms.Task(pfTrack, hgcalTrackCollection,
tpClusterProducer,
quickTrackAssociatorByHits, simPFProducer)
_phase2_hgcal_simPFSequence = cms.Sequence(_phase2_hgcal_simPFTask)
_phase2_hgcal_particleFlowRecoTask = cms.Task(_phase2_hgcal_simPFTask,
particleFlowRecoTask.copy())
_phase2_hgcal_particleFlowRecoTask.replace(
particleFlowTmpTask, cms.Task(particleFlowTmpBarrel, particleFlowTmp))
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toReplaceWith(particleFlowTmp, _phase2_hgcal_particleFlowTmp)
phase2_hgcal.toReplaceWith(particleFlowRecoTask,
_phase2_hgcal_particleFlowRecoTask)
from Configuration.Eras.Modifier_pp_on_XeXe_2017_cff import pp_on_XeXe_2017
from Configuration.ProcessModifiers.pp_on_AA_cff import pp_on_AA
for e in [pp_on_XeXe_2017, pp_on_AA]:
e.toModify(particleFlowDisplacedVertexCandidate,
tracksSelectorParameters=dict(pt_min=999999.0,
nChi2_max=0.0,
pt_min_prim=999999.0,
dxy=999999.0))
e.toModify(particleFlowBlock, useNuclear=cms.bool(False))
"Egamma/Electrons/Ele3_Et10")
dqmElectronClientSelectionEtIso = dqmElectronOfflineClient.clone()
dqmElectronClientSelectionEtIso.InputFolderName = cms.string(
"Egamma/Electrons/Ele4_Et10TkIso1")
dqmElectronClientSelectionEtIso.OutputFolderName = cms.string(
"Egamma/Electrons/Ele4_Et10TkIso1")
#dqmElectronClientSelectionEtIsoElID = dqmElectronOfflineClient.clone() ;
#dqmElectronClientSelectionEtIsoElID.InputFolderName = cms.string("Egamma/Electrons/Ele4_Et10TkIso1ElID") ;
#dqmElectronClientSelectionEtIsoElID.OutputFolderName = cms.string("Egamma/Electrons/Ele4_Et10TkIso1ElID") ;
dqmElectronClientTagAndProbe = dqmElectronOfflineClient.clone()
dqmElectronClientTagAndProbe.InputFolderName = cms.string(
"Egamma/Electrons/Ele5_TagAndProbe")
dqmElectronClientTagAndProbe.OutputFolderName = cms.string(
"Egamma/Electrons/Ele5_TagAndProbe")
dqmElectronClientTagAndProbe.EffHistoTitle = cms.string("")
electronOfflineClientSequence = cms.Sequence(
dqmElectronClientAllElectrons * dqmElectronClientSelectionEt *
dqmElectronClientSelectionEtIso
# * dqmElectronClientSelectionEtIsoElID
* dqmElectronClientTagAndProbe)
_electronOfflineClientSequenceHGC = electronOfflineClientSequence.copy()
_electronOfflineClientSequenceHGC += dqmElectronClientAllElectronsHGC
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toReplaceWith(electronOfflineClientSequence,
_electronOfflineClientSequenceHGC)
# Preshower clusters
'keep recoPreshowerClusters_multi5x5SuperClustersWithPreshower_*_*',
'keep recoPreshowerClusterShapes_multi5x5PreshowerClusterShape_*_*',
# Particle Flow superclusters (only SuperCluster and CaloCluster outputs, not association map from PFClusters)
'keep recoSuperClusters_particleFlowSuperClusterECAL_*_*',
'keep recoCaloClusters_particleFlowSuperClusterECAL_*_*'
)
)
_phase2_hgcal_scCommands = ['keep *_particleFlowSuperClusterHGCal_*_*']
_phase2_hgcal_scCommandsAOD = ['keep recoSuperClusters_particleFlowSuperClusterHGCal__*',
'keep recoCaloClusters_particleFlowSuperClusterHGCal__*']
_phase2_hgcal_RecoEcalFEVT = RecoEcalFEVT.clone()
_phase2_hgcal_RecoEcalFEVT.outputCommands += _phase2_hgcal_scCommands
_phase2_hgcal_RecoEcalRECO = RecoEcalRECO.clone()
_phase2_hgcal_RecoEcalRECO.outputCommands += _phase2_hgcal_scCommands
_phase2_hgcal_RecoEcalAOD = RecoEcalAOD.clone()
_phase2_hgcal_RecoEcalAOD.outputCommands += _phase2_hgcal_scCommandsAOD
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toReplaceWith( RecoEcalFEVT, _phase2_hgcal_RecoEcalFEVT )
phase2_hgcal.toReplaceWith( RecoEcalRECO, _phase2_hgcal_RecoEcalRECO )
phase2_hgcal.toReplaceWith( RecoEcalAOD , _phase2_hgcal_RecoEcalAOD )
from Configuration.Eras.Modifier_pA_2016_cff import pA_2016
from Configuration.Eras.Modifier_peripheralPbPb_cff import peripheralPbPb
#HI-specific products needed in pp scenario special configurations
for e in [pA_2016, peripheralPbPb]:
for ec in [RecoEcalRECO.outputCommands, RecoEcalFEVT.outputCommands]:
e.toModify( ec, func=lambda outputCommands: outputCommands.extend(['keep recoCaloClusters_islandBasicClusters_*_*'])
)
import FWCore.ParameterSet.Config as cms
#
# sequence to make photons from clusters in ECAL
#
# photon producer
from RecoEgamma.EgammaPhotonProducers.photonCore_cfi import *
from RecoEgamma.EgammaPhotonProducers.photons_cfi import *
photonTask = cms.Task(photonCore, photons)
photonSequence = cms.Sequence(photonTask)
_photonTaskFromMultiCl = photonTask.copy()
_photonTaskFromMultiCl.add(photonCoreFromMultiCl, photonsFromMultiCl)
_photonTaskWithIsland = photonTask.copy()
_photonTaskWithIsland.add(islandPhotonCore, islandPhotons)
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toReplaceWith(photonTask, _photonTaskFromMultiCl)
from Configuration.Eras.Modifier_pA_2016_cff import pA_2016
from Configuration.Eras.Modifier_peripheralPbPb_cff import peripheralPbPb
from Configuration.Eras.Modifier_pp_on_AA_2018_cff import pp_on_AA_2018
from Configuration.Eras.Modifier_pp_on_XeXe_2017_cff import pp_on_XeXe_2017
from Configuration.Eras.Modifier_ppRef_2017_cff import ppRef_2017
for e in [pA_2016, peripheralPbPb, pp_on_AA_2018, pp_on_XeXe_2017, ppRef_2017]:
e.toReplaceWith(photonTask, _photonTaskWithIsland)
simPFProducer )
_phase2_hgcal_particleFlowReco = cms.Sequence( _phase2_hgcal_simPFSequence * particleFlowReco.copy() )
_phase2_hgcal_particleFlowReco.replace( particleFlowTmpSeq, cms.Sequence( particleFlowTmpBarrel * particleFlowTmp ) )
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toModify( quickTrackAssociatorByHits,
pixelSimLinkSrc = cms.InputTag("simSiPixelDigis","Pixel"),
stripSimLinkSrc = cms.InputTag("simSiPixelDigis","Tracker")
)
phase2_hgcal.toModify( tpClusterProducer,
pixelSimLinkSrc = cms.InputTag("simSiPixelDigis", "Pixel"),
phase2OTSimLinkSrc = cms.InputTag("simSiPixelDigis","Tracker")
)
phase2_hgcal.toReplaceWith( particleFlowTmp, _phase2_hgcal_particleFlowTmp )
phase2_hgcal.toReplaceWith( particleFlowReco, _phase2_hgcal_particleFlowReco )
from Configuration.Eras.Modifier_pp_on_XeXe_2017_cff import pp_on_XeXe_2017
pp_on_XeXe_2017.toModify(particleFlowDisplacedVertexCandidate,
tracksSelectorParameters = dict(pt_min = 999999.0,
nChi2_max = 0.0,
pt_min_prim = 999999.0,
dxy = 999999.0)
)
pp_on_XeXe_2017.toModify(particleFlowBlock, useNuclear = cms.bool(False))
pp_on_XeXe_2017.toModify(pfNoPileUpIso, enable = cms.bool(False))
pp_on_XeXe_2017.toModify(pfPileUpIso, enable = cms.bool(False))
# Hcal part # # calo geometry # # changed by tommaso. now the calibrations are read from Configuration/StaqndardSequences/data/*Conditions.cff # #HCAL reconstruction from RecoLocalCalo.Configuration.hcalLocalReco_cff import * from RecoLocalCalo.Configuration.hcalGlobalReco_cff import * # # sequence CaloLocalReco and CaloGlobalReco # calolocalreco = cms.Sequence(ecalLocalRecoSequence+hcalLocalRecoSequence) caloglobalreco = cms.Sequence(hcalGlobalRecoSequence) from RecoLocalCalo.HcalRecProducers.HcalHitSelection_cfi import * reducedHcalRecHitsSequence = cms.Sequence( reducedHcalRecHits ) # # R.Ofierzynski (29.Oct.2009): add NZS sequence # from RecoLocalCalo.Configuration.hcalLocalRecoNZS_cff import * calolocalrecoNZS = cms.Sequence(ecalLocalRecoSequence+hcalLocalRecoSequence+hcalLocalRecoSequenceNZS) from RecoLocalCalo.Configuration.hgcalLocalReco_cff import * _phase2_calolocalreco = calolocalreco.copy() _phase2_calolocalreco += hgcalLocalRecoSequence from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal phase2_hgcal.toReplaceWith( calolocalreco , _phase2_calolocalreco )
# GEM settings _gem_RawToDigi = RawToDigi.copy() _gem_RawToDigi.insert(-1,muonGEMDigis) from Configuration.Eras.Modifier_run2_GEM_2017_cff import run2_GEM_2017 run2_GEM_2017.toReplaceWith(RawToDigi, _gem_RawToDigi) from Configuration.Eras.Modifier_run3_GEM_cff import run3_GEM run3_GEM.toReplaceWith(RawToDigi, _gem_RawToDigi) from EventFilter.HGCalRawToDigi.HGCalRawToDigi_cfi import * _hgcal_RawToDigi = RawToDigi.copy() _hgcal_RawToDigi += hgcalDigis from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal phase2_hgcal.toReplaceWith(RawToDigi,_hgcal_RawToDigi) # RPC New Readout Validation from Configuration.Eras.Modifier_stage2L1Trigger_2017_cff import stage2L1Trigger_2017 _rpc_NewReadoutVal_RawToDigi = RawToDigi.copy() _rpc_NewReadoutVal_RawToDigi_noTk = RawToDigi_noTk.copy() _rpc_NewReadoutVal_RawToDigi += muonRPCNewDigis _rpc_NewReadoutVal_RawToDigi_noTk += muonRPCNewDigis stage2L1Trigger_2017.toReplaceWith(RawToDigi, _rpc_NewReadoutVal_RawToDigi) stage2L1Trigger_2017.toReplaceWith(RawToDigi_noTk, _rpc_NewReadoutVal_RawToDigi) from Configuration.Eras.Modifier_fastSim_cff import fastSim fastSim.toReplaceWith(RawToDigi, RawToDigi.copyAndExclude([muonRPCNewDigis])) fastSim.toReplaceWith(RawToDigi_noTk, RawToDigi_noTk.copyAndExclude([muonRPCNewDigis])) _hfnose_RawToDigi = RawToDigi.copy()
#
# calo geometry
#
# changed by tommaso. now the calibrations are read from Configuration/StaqndardSequences/data/*Conditions.cff
#
#HCAL reconstruction
from RecoLocalCalo.Configuration.hcalLocalReco_cff import *
from RecoLocalCalo.Configuration.hcalGlobalReco_cff import *
#
# sequence CaloLocalReco and CaloGlobalReco
#
calolocalreco = cms.Sequence(ecalLocalRecoSequence + hcalLocalRecoSequence)
caloglobalreco = cms.Sequence(hcalGlobalRecoSequence)
from RecoLocalCalo.HcalRecProducers.HcalHitSelection_cfi import *
reducedHcalRecHitsSequence = cms.Sequence(reducedHcalRecHits)
#
# R.Ofierzynski (29.Oct.2009): add NZS sequence
#
from RecoLocalCalo.Configuration.hcalLocalRecoNZS_cff import *
calolocalrecoNZS = cms.Sequence(ecalLocalRecoSequence + hcalLocalRecoSequence +
hcalLocalRecoSequenceNZS)
from RecoLocalCalo.Configuration.hgcalLocalReco_cff import *
_phase2_calolocalreco = calolocalreco.copy()
_phase2_calolocalreco += hgcalLocalRecoSequence
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toReplaceWith(calolocalreco, _phase2_calolocalreco)
# GEM settings
_gem_RawToDigi = RawToDigi.copy()
_gem_RawToDigi.insert(-1, muonGEMDigis)
from Configuration.Eras.Modifier_run2_GEM_2017_cff import run2_GEM_2017
run2_GEM_2017.toReplaceWith(RawToDigi, _gem_RawToDigi)
from Configuration.Eras.Modifier_run3_GEM_cff import run3_GEM
run3_GEM.toReplaceWith(RawToDigi, _gem_RawToDigi)
from EventFilter.HGCalRawToDigi.HGCalRawToDigi_cfi import *
_hgcal_RawToDigi = RawToDigi.copy()
_hgcal_RawToDigi += hgcalDigis
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toReplaceWith(RawToDigi, _hgcal_RawToDigi)
# RPC New Readout Validation
from Configuration.Eras.Modifier_stage2L1Trigger_2017_cff import stage2L1Trigger_2017
_rpc_NewReadoutVal_RawToDigi = RawToDigi.copy()
_rpc_NewReadoutVal_RawToDigi_noTk = RawToDigi_noTk.copy()
_rpc_NewReadoutVal_RawToDigi += muonRPCNewDigis
_rpc_NewReadoutVal_RawToDigi_noTk += muonRPCNewDigis
stage2L1Trigger_2017.toReplaceWith(RawToDigi, _rpc_NewReadoutVal_RawToDigi)
stage2L1Trigger_2017.toReplaceWith(RawToDigi_noTk,
_rpc_NewReadoutVal_RawToDigi)
from Configuration.Eras.Modifier_fastSim_cff import fastSim
fastSim.toReplaceWith(RawToDigi, RawToDigi.copyAndExclude([muonRPCNewDigis]))
fastSim.toReplaceWith(RawToDigi_noTk,
RawToDigi_noTk.copyAndExclude([muonRPCNewDigis]))
# PFECAL super clusters, either hybrid-clustering clone (Box) or mustache.
particleFlowSuperClusteringTask = cms.Task(particleFlowSuperClusterECAL)
particleFlowSuperClusteringSequence = cms.Sequence(
particleFlowSuperClusteringTask)
particleFlowSuperClusterHGCal = particleFlowSuperClusterECAL.clone()
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toModify(
particleFlowSuperClusterHGCal,
PFClusters='particleFlowClusterHGCal',
useRegression=True,
PFBasicClusterCollectionEndcap="",
PFSuperClusterCollectionEndcap="",
PFSuperClusterCollectionEndcapWithPreshower="",
thresh_PFClusterEndcap=1.5e-1, # 150 MeV threshold
dropUnseedable=True,
)
phase2_hgcal.toModify(particleFlowSuperClusterHGCal.regressionConfig,
regressionKeyEE="superclus_hgcal_mean_offline",
uncertaintyKeyEE="superclus_hgcal_sigma_offline",
isPhaseII=True,
hgcalRecHits="particleFlowRecHitHGC")
_phase2_hgcal_particleFlowSuperClusteringTask = particleFlowSuperClusteringTask.copy(
)
_phase2_hgcal_particleFlowSuperClusteringTask.add(
particleFlowSuperClusterHGCal)
phase2_hgcal.toReplaceWith(particleFlowSuperClusteringTask,
_phase2_hgcal_particleFlowSuperClusteringTask)
metFilters = cms.Sequence(
HBHENoiseFilterResultProducer * HBHENoiseFilter * primaryVertexFilter *
# HBHENoiseIsoFilter*
# HcalStripHaloFilter *
CSCTightHaloFilter *
# hcalLaserEventFilter *
#Various proposals for updated halo filters.
##2015 proposals:
#CSCTightHaloTrkMuUnvetoFilter *
#CSCTightHalo2015Filter *
##2016 proposals
#globalTightHalo2016Filter*
#globalSuperTightHalo2016Filter*
EcalDeadCellTriggerPrimitiveFilter *
# *goodVertices * trackingFailureFilter *
eeBadScFilter *
# ecalLaserCorrFilter *
# trkPOGFilters
chargedHadronTrackResolutionFilter * muonBadTrackFilter)
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toReplaceWith(
metFilters,
metFilters.copyAndExclude([
HBHENoiseFilterResultProducer,
HBHENoiseFilter, # No hcalnoise for hgcal
eeBadScFilter # No EE
]))
SimL1TGlobal
)
SimL1Emulator = cms.Sequence( SimL1EmulatorCore )
#
# Emulators are configured from DB (GlobalTags)
#
from L1Trigger.L1TGlobal.GlobalParameters_cff import *
# 2017 EMTF and TwinMux emulator use payloads from DB, not yet in GT,
# soon to be removed when availble in GTs
from L1Trigger.L1TTwinMux.fakeTwinMuxParams_cff import *
# Customisation for the phase2_hgcal era. Includes the HGCAL L1 trigger
from L1Trigger.L1THGCal.hgcalTriggerPrimitives_cff import *
_phase2_siml1emulator = SimL1Emulator.copy()
_phase2_siml1emulator += hgcalTriggerPrimitives
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toReplaceWith( SimL1Emulator , _phase2_siml1emulator )
# If PreMixing, don't run these modules during first step
from Configuration.ProcessModifiers.premix_stage1_cff import premix_stage1
premix_stage1.toReplaceWith(SimL1Emulator, SimL1Emulator.copyAndExclude([
SimL1TCalorimeter,
SimL1TechnicalTriggers,
SimL1TGlobal
]))
result=cms.bool(False)))
return process
from Configuration.Eras.Modifier_phase2_common_cff import phase2_common
phase2_common.toReplaceWith(Flag_trkPOG_manystripclus53X, cms.Path())
phase2_common.toReplaceWith(Flag_trkPOG_toomanystripclus53X, cms.Path())
phase2_common.toReplaceWith(Flag_trkPOGFilters,
cms.Path(~logErrorTooManyClusters))
from Configuration.Eras.Modifier_run2_common_cff import run2_common
run2_common.toReplaceWith(Flag_ecalBadCalibFilter,
cms.Path(ecalBadCalibFilter))
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toReplaceWith(Flag_HBHENoiseFilter, cms.Path())
phase2_hgcal.toReplaceWith(Flag_HBHENoiseIsoFilter, cms.Path())
phase2_hgcal.toReplaceWith(Flag_eeBadScFilter, cms.Path())
metFilterPathsTask = cms.Task(
HBHENoiseFilterResultProducer, HBHENoiseFilter, HBHENoiseIsoFilter,
hcalLaserEventFilter, EcalDeadCellTriggerPrimitiveFilter, eeBadScFilter,
ecalLaserCorrFilter, EcalDeadCellBoundaryEnergyFilter, ecalBadCalibFilter,
primaryVertexFilter, CSCTightHaloFilter, CSCTightHaloTrkMuUnvetoFilter,
CSCTightHalo2015Filter, globalTightHalo2016Filter,
globalSuperTightHalo2016Filter, HcalStripHaloFilter, goodVertices,
trackingFailureFilter, manystripclus53X, toomanystripclus53X,
logErrorTooManyClusters, chargedHadronTrackResolutionFilter,
muonBadTrackFilter, BadChargedCandidateFilter, BadPFMuonFilter,
BadPFMuonDzFilter, BadChargedCandidateSummer16Filter,
BadPFMuonSummer16Filter)
# Sequence to make final electrons.
# In the past, this was including the seeding, but this one is directly
# imported in the reco sequences since the integration with pflow.
#==============================================================================
from RecoEgamma.EgammaElectronProducers.gsfElectronModules_cff import *
gsfElectronSequence = cms.Sequence(ecalDrivenGsfElectronCores *
ecalDrivenGsfElectrons * gsfElectronCores *
gsfElectrons)
gsfEcalDrivenElectronSequence = cms.Sequence(ecalDrivenGsfElectronCores *
ecalDrivenGsfElectrons)
_gsfEcalDrivenElectronSequenceFromMultiCl = gsfEcalDrivenElectronSequence.copy(
)
_gsfEcalDrivenElectronSequenceFromMultiCl += cms.Sequence(
ecalDrivenGsfElectronCoresFromMultiCl * ecalDrivenGsfElectronsFromMultiCl)
#gsfElectronMergingSequence = cms.Sequence(gsfElectronCores*gsfElectrons)
from RecoEgamma.EgammaElectronProducers.edBasedElectronIso_cff import *
from RecoEgamma.EgammaElectronProducers.pfBasedElectronIso_cff import *
electronIsoSequence = cms.Sequence(edBasedElectronIsoSequence +
pfBasedElectronIsoSequence)
gsfElectronMergingSequence = cms.Sequence(electronIsoSequence *
gsfElectronCores * gsfElectrons)
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toReplaceWith(gsfEcalDrivenElectronSequence,
_gsfEcalDrivenElectronSequenceFromMultiCl)
pfClusteringECALTask,
pfClusteringHBHEHFTask,
pfClusteringHOTask
)
particleFlowCluster = cms.Sequence(particleFlowClusterTask)
#HGCal
from RecoParticleFlow.PFClusterProducer.particleFlowRecHitHGC_cfi import particleFlowRecHitHGC
pfClusteringHGCal = cms.Sequence(particleFlowRecHitHGC)
_phase2_hgcal_particleFlowCluster = particleFlowCluster.copy()
_phase2_hgcal_particleFlowCluster += pfClusteringHGCal
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toReplaceWith( particleFlowCluster, _phase2_hgcal_particleFlowCluster )
#timing
from RecoParticleFlow.PFClusterProducer.particleFlowClusterTimeAssigner_cfi import particleFlowTimeAssignerECAL
from RecoParticleFlow.PFSimProducer.ecalBarrelClusterFastTimer_cfi import ecalBarrelClusterFastTimer
_phase2_timing_particleFlowClusterECALTask = particleFlowClusterECALTask.copy()
_phase2_timing_particleFlowClusterECALTask.add(cms.Task(ecalBarrelClusterFastTimer,
particleFlowTimeAssignerECAL))
from Configuration.Eras.Modifier_phase2_timing_cff import phase2_timing
phase2_timing.toReplaceWith(particleFlowClusterECALTask,
_phase2_timing_particleFlowClusterECALTask)
phase2_timing.toModify(particleFlowClusterECAL,
inputECAL = cms.InputTag('particleFlowTimeAssignerECAL'))
piZeroAnalysis.OutputMEsInRootFile = cms.bool(False)
piZeroAnalysis.Verbosity = cms.untracked.int32(0)
piZeroAnalysis.standAlone = cms.bool(False)
zmumugammaOldAnalysis = DQMOffline.EGamma.zmumugammaAnalyzer_cfi.zmumugammaAnalysis.clone()
zmumugammaOldAnalysis.ComponentName = cms.string('zmumugammaOldAnalysis')
zmumugammaOldAnalysis.analyzerName = cms.string('zmumugammaOldValidation')
zmumugammaOldAnalysis.phoProducer = cms.InputTag('photons')
# HGCal customizations
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
stdPhotonAnalysisHGCal = stdPhotonAnalysis.clone()
stdPhotonAnalysisHGCal.ComponentName = 'stdPhotonAnalyzerHGCalFromMultiCl'
stdPhotonAnalysisHGCal.analyzerName = 'stdPhotonAnalyzerHGCalFromMultiCl'
stdPhotonAnalysisHGCal.phoProducer = 'photonsFromMultiCl'
stdPhotonAnalysisHGCal.isolationStrength = 2
stdPhotonAnalysisHGCal.etaMin = -3.0
stdPhotonAnalysisHGCal.etaMax = 3.0
stdPhotonAnalysisHGCal.maxPhoEta = 3.0
egammaDQMOffline = cms.Sequence(photonAnalysis*stdPhotonAnalysis*zmumugammaOldAnalysis*zmumugammaAnalysis*piZeroAnalysis*electronAnalyzerSequence)
_egammaDQMOfflineHGCal = egammaDQMOffline.copy()
_egammaDQMOfflineHGCal += stdPhotonAnalysisHGCal
phase2_hgcal.toReplaceWith(
egammaDQMOffline, _egammaDQMOfflineHGCal
)
from RecoEgamma.EgammaElectronProducers.gsfElectronModules_cff import *
gsfElectronTask = cms.Task(ecalDrivenGsfElectronCores,ecalDrivenGsfElectrons,gsfElectronCores,gsfElectrons)
gsfElectronSequence = cms.Sequence(gsfElectronTask)
gsfEcalDrivenElectronTask = cms.Task(ecalDrivenGsfElectronCores,ecalDrivenGsfElectrons)
gsfEcalDrivenElectronSequence = cms.Sequence(gsfEcalDrivenElectronTask)
_gsfEcalDrivenElectronTaskFromMultiCl = gsfEcalDrivenElectronTask.copy()
_gsfEcalDrivenElectronTaskFromMultiCl.add(cms.Task(ecalDrivenGsfElectronCoresFromMultiCl,ecalDrivenGsfElectronsFromMultiCl))
_gsfEcalDrivenElectronSequenceFromMultiCl = cms.Sequence(_gsfEcalDrivenElectronTaskFromMultiCl)
#gsfElectronMergingSequence = cms.Sequence(gsfElectronCores*gsfElectrons)
from RecoEgamma.EgammaElectronProducers.edBasedElectronIso_cff import *
from RecoEgamma.EgammaElectronProducers.pfBasedElectronIso_cff import *
electronIsoTask = cms.Task(
edBasedElectronIsoTask,
pfBasedElectronIsoTask
)
electronIsoSequence = cms.Sequence(electronIsoTask)
gsfElectronMergingTask = cms.Task(electronIsoTask,gsfElectronCores,gsfElectrons)
gsfElectronMergingSequence = cms.Sequence(gsfElectronMergingTask)
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toReplaceWith(
gsfEcalDrivenElectronTask, _gsfEcalDrivenElectronTaskFromMultiCl
)
postValidationCosmics = cms.Sequence(
postProcessorMuonTrack
)
postValidationMiniAOD = cms.Sequence(
electronPostValidationSequenceMiniAOD
)
_phase1_postValidation = postValidation.copy()
_phase1_postValidation += siPixelPhase1OfflineDQM_harvestingV
from Configuration.Eras.Modifier_phase1Pixel_cff import phase1Pixel
phase1Pixel.toReplaceWith( postValidation, _phase1_postValidation )
_run3_postValidation = postValidation.copy()
_run3_postValidation += MuonGEMHitsPostProcessors
_run3_postValidation += MuonGEMDigisPostProcessors
_run3_postValidation += MuonGEMRecHitsPostProcessors
_phase2_postValidation = _run3_postValidation.copy()
_phase2_postValidation += hgcalPostProcessor
_phase2_postValidation += MuonME0DigisPostProcessors
_phase2_postValidation += MuonME0SegPostProcessors
from Configuration.Eras.Modifier_run2_GEM_2017_MCTest_cff import run2_GEM_2017_MCTest
run2_GEM_2017_MCTest.toReplaceWith( postValidation, _run3_postValidation )
from Configuration.Eras.Modifier_run3_GEM_cff import run3_GEM
run3_GEM.toReplaceWith( postValidation, _run3_postValidation )
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toReplaceWith( postValidation, _phase2_postValidation )
from L1Trigger.Configuration.SimL1TechnicalTriggers_cff import *
from L1Trigger.L1TCalorimeter.simDigis_cff import *
from L1Trigger.L1TMuon.simDigis_cff import *
from L1Trigger.L1TGlobal.simDigis_cff import *
# define a core which can be extented in customizations:
SimL1EmulatorCore = cms.Sequence(SimL1TCalorimeter + SimL1TMuon +
SimL1TechnicalTriggers + SimL1TGlobal)
SimL1Emulator = cms.Sequence(SimL1EmulatorCore)
#
# Emulators are configured from DB (GlobalTags)
#
from L1Trigger.L1TGlobal.GlobalParameters_cff import *
# 2017 EMTF and TwinMux emulator use payloads from DB, not yet in GT,
# soon to be removed when availble in GTs
from L1Trigger.L1TTwinMux.fakeTwinMuxParams_cff import *
# Customisation for the phase2_hgcal era. Includes the HGCAL L1 trigger
from L1Trigger.L1THGCal.hgcalTriggerPrimitives_cff import *
_phase2_siml1emulator = SimL1Emulator.copy()
_phase2_siml1emulator += hgcalTriggerPrimitives
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
from Configuration.Eras.Modifier_phase2_hgcalV9_cff import phase2_hgcalV9
phase2_hgcal.toReplaceWith(SimL1Emulator, _phase2_siml1emulator)
# Producer for energy corrections
#from RecoEcal.EgammaClusterProducers.correctedDynamicHybridSuperClusters_cfi import *
# PFECAL super clusters, either hybrid-clustering clone (Box) or mustache.
particleFlowSuperClusteringSequence = cms.Sequence(particleFlowSuperClusterECAL)
particleFlowSuperClusterHGCal = particleFlowSuperClusterECAL.clone()
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toModify(
particleFlowSuperClusterHGCal,
PFClusters = cms.InputTag('particleFlowClusterHGCal'),
useRegression = cms.bool(False), #no HGCal regression yet
use_preshower = cms.bool(False),
PFBasicClusterCollectionEndcap = cms.string(""),
PFSuperClusterCollectionEndcap = cms.string(""),
PFSuperClusterCollectionEndcapWithPreshower = cms.string(""),
thresh_PFClusterEndcap = cms.double(1.5e-1), # 150 MeV threshold
dropUnseedable = cms.bool(True),
)
particleFlowSuperClusterHGCalFromMultiCl = particleFlowSuperClusterHGCal.clone()
phase2_hgcal.toModify(
particleFlowSuperClusterHGCalFromMultiCl,
PFClusters = cms.InputTag('particleFlowClusterHGCalFromMultiCl')
)
_phase2_hgcal_particleFlowSuperClusteringSequence = particleFlowSuperClusteringSequence.copy()
_phase2_hgcal_particleFlowSuperClusteringSequence += particleFlowSuperClusterHGCal
_phase2_hgcal_particleFlowSuperClusteringSequence += particleFlowSuperClusterHGCalFromMultiCl
phase2_hgcal.toReplaceWith( particleFlowSuperClusteringSequence, _phase2_hgcal_particleFlowSuperClusteringSequence )
_run3_globalrecoTask = globalrecoTask.copyAndExclude([CastorFullRecoTask])
run3_common.toReplaceWith(globalrecoTask, _run3_globalrecoTask)
_fastSim_globalrecoTask = globalrecoTask.copyAndExclude(
[CastorFullRecoTask, muoncosmicrecoTask])
# insert the few tracking modules to be run after mixing back in the globalreco sequence
_fastSim_globalrecoTask.add(newCombinedSeeds, trackExtrapolator,
caloTowerForTrk, firstStepPrimaryVerticesUnsorted,
ak4CaloJetsForTrk, initialStepTrackRefsForJets,
firstStepPrimaryVertices)
fastSim.toReplaceWith(globalrecoTask, _fastSim_globalrecoTask)
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
_phase2HGALRecoTask = globalrecoTask.copy()
_phase2HGALRecoTask.add(iterTICLTask)
phase2_hgcal.toReplaceWith(globalrecoTask, _phase2HGALRecoTask)
globalreco_plusPLTask = cms.Task(globalrecoTask, ctfTracksPixelLessTask)
globalreco_plusPL = cms.Sequence(globalreco_plusPLTask)
reducedRecHitsTask = cms.Task(reducedEcalRecHitsTask, reducedHcalRecHitsTask)
reducedRecHits = cms.Sequence(reducedRecHitsTask)
highlevelrecoTask = cms.Task(
egammaHighLevelRecoPrePFTask, particleFlowRecoTask,
egammaHighLevelRecoPostPFTask, muoncosmichighlevelrecoTask,
muonshighlevelrecoTask, particleFlowLinksTask, jetHighLevelRecoTask,
metrecoPlusHCALNoiseTask, btaggingTask, recoPFMETTask, PFTauTask,
reducedRecHitsTask, cosmicDCTracksSeqTask, lowPtGsfElectronTask,
conversionOpenTrackTask, gsfTracksOpenConversions)
highlevelreco = cms.Sequence(highlevelrecoTask)
metFilters = cms.Sequence(
HBHENoiseFilterResultProducer *
HBHENoiseFilter *
primaryVertexFilter*
# HBHENoiseIsoFilter*
# HcalStripHaloFilter *
CSCTightHaloFilter *
# hcalLaserEventFilter *
#Various proposals for updated halo filters.
##2015 proposals:
#CSCTightHaloTrkMuUnvetoFilter *
#CSCTightHalo2015Filter *
##2016 proposals
#globalTightHalo2016Filter*
#globalSuperTightHalo2016Filter*
EcalDeadCellTriggerPrimitiveFilter*
# *goodVertices * trackingFailureFilter *
eeBadScFilter*
# ecalLaserCorrFilter *
# trkPOGFilters
chargedHadronTrackResolutionFilter *
muonBadTrackFilter
)
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toReplaceWith(metFilters, metFilters.copyAndExclude([
HBHENoiseFilterResultProducer, HBHENoiseFilter, # No hcalnoise for hgcal
eeBadScFilter # No EE
]))
_ctpps_2016_RawToDigiTask = RawToDigiTask.copy() _ctpps_2016_RawToDigiTask.add(ctppsRawToDigiTask) ctpps_2016.toReplaceWith(RawToDigiTask, _ctpps_2016_RawToDigiTask) _ctpps_2016_RawToDigiTask_noTk = RawToDigiTask_noTk.copy() _ctpps_2016_RawToDigiTask_noTk.add(ctppsRawToDigiTask) ctpps_2016.toReplaceWith(RawToDigiTask_noTk, _ctpps_2016_RawToDigiTask_noTk) # GEM settings _gem_RawToDigiTask = RawToDigiTask.copy() _gem_RawToDigiTask.add(muonGEMDigis) from Configuration.Eras.Modifier_run2_GEM_2017_cff import run2_GEM_2017 run2_GEM_2017.toReplaceWith(RawToDigiTask, _gem_RawToDigiTask) from Configuration.Eras.Modifier_run3_GEM_cff import run3_GEM run3_GEM.toReplaceWith(RawToDigiTask, _gem_RawToDigiTask) from EventFilter.HGCalRawToDigi.HGCalRawToDigi_cfi import * _hgcal_RawToDigiTask = RawToDigiTask.copy() _hgcal_RawToDigiTask.add(hgcalDigis) from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal phase2_hgcal.toReplaceWith(RawToDigiTask, _hgcal_RawToDigiTask) _hfnose_RawToDigiTask = RawToDigiTask.copy() _hfnose_RawToDigiTask.add(hfnoseDigis) from Configuration.Eras.Modifier_phase2_hfnose_cff import phase2_hfnose phase2_hfnose.toReplaceWith(RawToDigiTask, _hfnose_RawToDigiTask)
# * dqmElectronAnalysisSelectionEtIsoElID
* dqmElectronTagProbeAnalysis
)
mergedSuperClustersHGC = mergedSuperClusters.clone(
src = (
("particleFlowSuperClusterECAL:particleFlowSuperClusterECALBarrel"),
("particleFlowSuperClusterHGCal")
)
)
dqmElectronAnalysisAllElectronsHGC = dqmElectronAnalysisAllElectrons.clone(
OutputFolderName = 'Egamma/Electrons/Ele2HGC_All',
MaxAbsEtaMatchingObject = 3.0,
EtaMax = 3.0,
EtaMin = -3.0,
MaxAbsEta = 3.0,
ElectronCollection = 'ecalDrivenGsfElectronsHGC',
MatchingObjectCollection = 'mergedSuperClustersHGC'
)
_electronAnalyzerSequenceHGC = electronAnalyzerSequence.copy()
_electronAnalyzerSequenceHGC += cms.Sequence(mergedSuperClustersHGC+dqmElectronAnalysisAllElectronsHGC)
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toModify( mergedSuperClusters, src = cms.VInputTag( cms.InputTag("particleFlowSuperClusterECAL","particleFlowSuperClusterECALBarrel"), cms.InputTag("particleFlowSuperClusterHGCal","") ) )
phase2_hgcal.toReplaceWith(
electronAnalyzerSequence, _electronAnalyzerSequenceHGC
)
import FWCore.ParameterSet.Config as cms
from DQMOffline.EGamma.photonOfflineClient_cfi import *
import DQMOffline.EGamma.photonOfflineClient_cfi
stdPhotonOfflineClient = DQMOffline.EGamma.photonOfflineClient_cfi.photonOfflineClient.clone(
ComponentName = 'stdPhotonOfflineClient',
analyzerName = 'stdPhotonAnalyzer'
)
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
stdPhotonOfflineClientHGCal = stdPhotonOfflineClient.clone(
analyzerName = 'stdPhotonAnalyzerHGCal'
)
from DQMOffline.EGamma.egammaDQMOffline_cff import stdPhotonAnalysisHGCal
stdPhotonOfflineClientHGCal.etaBin = stdPhotonAnalysisHGCal.etaBin
stdPhotonOfflineClientHGCal.etaMin = stdPhotonAnalysisHGCal.etaMin
stdPhotonOfflineClientHGCal.etaMax = stdPhotonAnalysisHGCal.etaMax
photonOfflineDQMClient = cms.Sequence(photonOfflineClient*stdPhotonOfflineClient)
_photonOfflineDQMClientHGCal = photonOfflineDQMClient.copy()
_photonOfflineDQMClientHGCal += stdPhotonOfflineClientHGCal
phase2_hgcal.toReplaceWith(
photonOfflineDQMClient, _photonOfflineDQMClientHGCal
)
def miniAOD_customizeMETFiltersFastSim(process):
"""Replace some MET filters that don't work in FastSim with trivial bools"""
for X in 'CSCTightHaloFilter', 'CSCTightHaloTrkMuUnvetoFilter','CSCTightHalo2015Filter','globalTightHalo2016Filter','globalSuperTightHalo2016Filter','HcalStripHaloFilter':
process.globalReplace(X, cms.EDFilter("HLTBool", result=cms.bool(True)))
for X in 'manystripclus53X', 'toomanystripclus53X', 'logErrorTooManyClusters':
process.globalReplace(X, cms.EDFilter("HLTBool", result=cms.bool(False)))
return process
from Configuration.Eras.Modifier_phase2_common_cff import phase2_common
phase2_common.toReplaceWith( Flag_trkPOG_manystripclus53X, cms.Path() )
phase2_common.toReplaceWith( Flag_trkPOG_toomanystripclus53X, cms.Path() )
phase2_common.toReplaceWith( Flag_trkPOGFilters, cms.Path(~logErrorTooManyClusters) )
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toReplaceWith( Flag_HBHENoiseFilter, cms.Path() )
phase2_hgcal.toReplaceWith( Flag_HBHENoiseIsoFilter, cms.Path() )
phase2_hgcal.toReplaceWith( Flag_eeBadScFilter, cms.Path() )
metFilterPathsTask = cms.Task(
HBHENoiseFilterResultProducer,
HBHENoiseFilter,
HBHENoiseIsoFilter,
hcalLaserEventFilter,
EcalDeadCellTriggerPrimitiveFilter,
eeBadScFilter,
ecalLaserCorrFilter,
EcalDeadCellBoundaryEnergyFilter,
primaryVertexFilter,
CSCTightHaloFilter,
CSCTightHaloTrkMuUnvetoFilter,
postValidation_gen = cms.Sequence(EventGeneratorPostProcessor) postValidationCosmics = cms.Sequence(postProcessorMuonTrack) postValidationMiniAOD = cms.Sequence(electronPostValidationSequenceMiniAOD) postValidationOuterTracker = cms.Sequence(OuterTracker_harvestingV) _phase1_postValidation = postValidation.copy() _phase1_postValidation += siPixelPhase1OfflineDQM_harvestingV from Configuration.Eras.Modifier_phase1Pixel_cff import phase1Pixel phase1Pixel.toReplaceWith(postValidation, _phase1_postValidation) _run3_postValidation = postValidation.copy() _run3_postValidation += MuonGEMHitsPostProcessors _run3_postValidation += MuonGEMDigisPostProcessors _run3_postValidation += MuonGEMRecHitsPostProcessors _phase2_postValidation = _run3_postValidation.copy() _phase2_postValidation += hgcalPostProcessor _phase2_postValidation += MuonME0DigisPostProcessors _phase2_postValidation += MuonME0SegPostProcessors from Configuration.Eras.Modifier_run2_GEM_2017_cff import run2_GEM_2017 run2_GEM_2017.toReplaceWith(postValidation, _run3_postValidation) from Configuration.Eras.Modifier_run3_GEM_cff import run3_GEM run3_GEM.toReplaceWith(postValidation, _run3_postValidation) from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal phase2_hgcal.toReplaceWith(postValidation, _phase2_postValidation)
* dqmElectronAnalysisSelectionEtIso
# * dqmElectronAnalysisSelectionEtIsoElID
* dqmElectronTagProbeAnalysis
)
mergedSuperClustersFromMultiCl = mergedSuperClusters.clone()
mergedSuperClustersFromMultiCl.src = cms.VInputTag(
cms.InputTag("particleFlowSuperClusterECAL","particleFlowSuperClusterECALBarrel"),
cms.InputTag("particleFlowSuperClusterHGCalFromMultiCl","")
)
dqmElectronAnalysisAllElectronsFromMultiCl = dqmElectronAnalysisAllElectrons.clone()
dqmElectronAnalysisAllElectronsFromMultiCl.OutputFolderName = 'Egamma/Electrons/Ele2FromMultiCl_All'
dqmElectronAnalysisAllElectronsFromMultiCl.MaxAbsEtaMatchingObject = 3.0
dqmElectronAnalysisAllElectronsFromMultiCl.EtaMax = 3.0
dqmElectronAnalysisAllElectronsFromMultiCl.EtaMin = -3.0
dqmElectronAnalysisAllElectronsFromMultiCl.MaxAbsEta = 3.0
dqmElectronAnalysisAllElectronsFromMultiCl.ElectronCollection = 'ecalDrivenGsfElectronsFromMultiCl'
dqmElectronAnalysisAllElectronsFromMultiCl.MatchingObjectCollection = 'mergedSuperClustersFromMultiCl'
_electronAnalyzerSequenceFromMultiCl = electronAnalyzerSequence.copy()
_electronAnalyzerSequenceFromMultiCl += cms.Sequence(mergedSuperClustersFromMultiCl+dqmElectronAnalysisAllElectronsFromMultiCl)
from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal
phase2_hgcal.toModify( mergedSuperClusters, src = cms.VInputTag( cms.InputTag("particleFlowSuperClusterECAL","particleFlowSuperClusterECALBarrel"), cms.InputTag("particleFlowSuperClusterHGCal","") ) )
phase2_hgcal.toReplaceWith(
electronAnalyzerSequence, _electronAnalyzerSequenceFromMultiCl
)
