from Configuration.ProcessModifiers.pp_on_AA_cff import pp_on_AA #products from regular pp which does not fit the normal AOD for e in [pA_2016, peripheralPbPb, pp_on_XeXe_2017, pp_on_AA]: e.toModify(RecoJetsAOD.outputCommands, func=lambda outputCommands: outputCommands.extend( ['keep *_towerMaker_*_*'])) for e in [pp_on_XeXe_2017, pp_on_AA]: e.toModify(RecoJetsAOD.outputCommands, func=lambda outputCommands: outputCommands.extend([ 'keep recoCentrality*_hiCentrality_*_*', 'keep recoClusterCompatibility*_hiClusterCompatibility_*_*' ])) #HI-specific products: needed in AOD, propagate to more inclusive tiers as well pA_2016.toModify( RecoJetsAOD.outputCommands, func=lambda outputCommands: outputCommands.extend([ 'keep recoCentrality*_pACentrality_*_*', 'keep *_hiFJGridEmptyAreaCalculator_*_*', 'keep *_hiFJRhoProducer_*_*' ])) #HI-specific products: needed in AOD, propagate to more inclusive tiers as well peripheralPbPb.toModify(RecoJetsAOD.outputCommands, func=lambda outputCommands: outputCommands.extend( ['keep recoCentrality*_pACentrality_*_*'])) pp_on_AA.toModify(RecoJetsAOD.outputCommands, func=lambda outputCommands: outputCommands.extend([ 'keep *_hiCentrality_*_*', 'keep *_hiFJRhoProducer_*_*', 'keep *_akPu3PFJets_*_*', 'keep *_akPu4PFJets_*_*', 'keep *_kt4PFJetsForRho_*_*', 'keep *_akCs4PFJets_*_*', 'keep *_akPu4CaloJets_*_*', 'drop *_caloTowers_*_*' ])) #RECO content
import FWCore.ParameterSet.Config as cms hiFJGridEmptyAreaCalculator = cms.EDProducer('HiFJGridEmptyAreaCalculator', gridWidth = cms.double(0.05), bandWidth = cms.double(0.2), mapEtaEdges = cms.InputTag('hiFJRhoProducer','mapEtaEdges'), mapToRho = cms.InputTag('hiFJRhoProducer','mapToRho'), mapToRhoM = cms.InputTag('hiFJRhoProducer','mapToRhoM'), pfCandSource = cms.InputTag('particleFlow'), jetSource = cms.InputTag('kt4PFJetsForRho'), doCentrality = cms.bool(True), hiBinCut = cms.int32(100), CentralityBinSrc = cms.InputTag("centralityBin","HFtowers"), keepGridInfo = cms.bool(False), ) from Configuration.Eras.Modifier_pA_2016_cff import pA_2016 pA_2016.toModify(hiFJGridEmptyAreaCalculator, doCentrality = False)
from RecoJets.Configuration.JetIDProducers_cff import * from RecoJets.Configuration.RecoTrackJets_cff import * from RecoJets.Configuration.RecoJetAssociations_cff import * from RecoJets.Configuration.RecoPFJets_cff import * from RecoJets.Configuration.RecoJPTJets_cff import * from JetMETCorrections.Configuration.JetCorrectorsForReco_cff import * jetGlobalReco = cms.Sequence(recoJets*recoJetIds*recoTrackJets) jetHighLevelReco = cms.Sequence(recoPFJets*jetCorrectorsForReco*recoJetAssociations*recoJetAssociationsExplicit*recoJPTJets) from Configuration.Eras.Modifier_pA_2016_cff import pA_2016 #HI-specific algorithms needed in pp scenario special configurations from RecoHI.HiJetAlgos.hiFJRhoProducer import hiFJRhoProducer from RecoHI.HiJetAlgos.hiFJGridEmptyAreaCalculator_cff import hiFJGridEmptyAreaCalculator pA_2016.toModify(hiFJGridEmptyAreaCalculator, doCentrality = False) kt4PFJetsForRho = kt4PFJets.clone(doAreaFastjet = True, jetPtMin = 0.0, GhostArea = 0.005) from RecoHI.HiCentralityAlgos.pACentrality_cfi import pACentrality pA_2016.toModify(pACentrality, producePixelTracks = False) _jetHighLevelReco_pA = jetHighLevelReco.copy() _jetHighLevelReco_pA += kt4PFJetsForRho _jetHighLevelReco_pA += hiFJRhoProducer _jetHighLevelReco_pA += hiFJGridEmptyAreaCalculator _jetHighLevelReco_pA += pACentrality pA_2016.toReplaceWith(jetHighLevelReco, _jetHighLevelReco_pA)
) ) RecoLocalCaloFEVT.outputCommands.extend(ecalLocalRecoFEVT.outputCommands) RecoLocalCaloRECO.outputCommands.extend(ecalLocalRecoRECO.outputCommands) RecoLocalCaloAOD.outputCommands.extend(ecalLocalRecoAOD.outputCommands) def _updateOutput( era, outputPSets, commands): for o in outputPSets: era.toModify( o, outputCommands = o.outputCommands + commands ) # mods for HGCAL from Configuration.Eras.Modifier_phase2_hgcal_cff import phase2_hgcal phase2_hgcal.toModify( RecoLocalCaloFEVT, outputCommands = RecoLocalCaloFEVT.outputCommands + [ 'keep *_HGCalRecHit_*_*', 'keep *_HGCalUncalibRecHit_*_*', 'keep *_hgcalLayerClusters_*_*' ] ) phase2_hgcal.toModify( RecoLocalCaloRECO, outputCommands = RecoLocalCaloRECO.outputCommands + ['keep *_HGCalRecHit_*_*','keep *_hgcalLayerClusters_*_*'] ) # don't modify AOD for HGCal yet, need "reduced" rechits collection first (i.e. requires reconstruction) phase2_hgcal.toModify( RecoLocalCaloAOD, outputCommands = RecoLocalCaloAOD.outputCommands + ['keep *_HGCalRecHit_*_*','keep *_hgcalLayerClusters_*_*'] ) from Configuration.Eras.Modifier_pA_2016_cff import pA_2016 pA_2016.toModify( RecoLocalCaloAOD.outputCommands, func=lambda outputCommands: outputCommands.extend(['keep *_zdcreco_*_*', 'keep ZDCDataFramesSorted_hcalDigis_*_*', 'keep ZDCDataFramesSorted_castorDigis_*_*', 'keep QIE10DataFrameHcalDataFrameContainer_hcalDigis_ZDC_*' ]) )
hltMuonOfflineAnalyzers = cms.Sequence( globalAnalyzerTnP * trackerAnalyzerTnP * tightAnalyzerTnP * looseAnalyzerTnP * globalAnalyzerRef * trackerAnalyzerRef * tightAnalyzerRef * looseAnalyzerRef * globalAnalyzerRef19 * trackerAnalyzerRef19 * tightAnalyzerRef19 * looseAnalyzerRef19 ) from Configuration.Eras.Modifier_pA_2016_cff import pA_2016 for muAna in [globalAnalyzerTnP.targetParams, trackerAnalyzerTnP.targetParams, tightAnalyzerTnP.targetParams, looseAnalyzerTnP.targetParams, globalAnalyzerRef.targetParams, trackerAnalyzerRef.targetParams, tightAnalyzerRef.targetParams, looseAnalyzerRef.targetParams,]: pA_2016.toModify(muAna, ptCut_Jpsi = cms.untracked.double( 5.0)) for muAna in [globalAnalyzerTnP.binParams, trackerAnalyzerTnP.binParams, tightAnalyzerTnP.binParams, looseAnalyzerTnP.binParams, globalAnalyzerRef.binParams, trackerAnalyzerRef.binParams, tightAnalyzerRef.binParams, looseAnalyzerRef.binParams]: pA_2016.toModify(muAna, ptCoarse = cms.untracked.vdouble(0.,1.,2.,3.,4.,5.,7.,9.,12.,15.,20.,30.,40.))
import FWCore.ParameterSet.Config as cms hiFJRhoProducer = cms.EDProducer('HiFJRhoProducer', jetSource=cms.InputTag('kt4PFJetsForRho'), nExcl=cms.int32(2), etaMaxExcl=cms.double(2.), ptMinExcl=cms.double(20.), nExcl2=cms.int32(1), etaMaxExcl2=cms.double(3.), ptMinExcl2=cms.double(20.), etaRanges=cms.vdouble(-5., -3., -2.1, -1.3, 1.3, 2.1, 3., 5.)) from Configuration.Eras.Modifier_pA_2016_cff import pA_2016 pA_2016.toModify(hiFJRhoProducer, etaRanges=cms.vdouble(-5., -3., -2., -1.5, -1., 1., 1.5, 2., 3., 5.))
import FWCore.ParameterSet.Config as cms hiFJGridEmptyAreaCalculator = cms.EDProducer( 'HiFJGridEmptyAreaCalculator', gridWidth=cms.double(0.05), bandWidth=cms.double(0.2), mapEtaEdges=cms.InputTag('hiFJRhoProducer', 'mapEtaEdges'), mapToRho=cms.InputTag('hiFJRhoProducer', 'mapToRho'), mapToRhoM=cms.InputTag('hiFJRhoProducer', 'mapToRhoM'), pfCandSource=cms.InputTag('particleFlow'), jetSource=cms.InputTag('kt4PFJetsForRho'), doCentrality=cms.bool(True), hiBinCut=cms.int32(100), CentralityBinSrc=cms.InputTag("centralityBin", "HFtowers"), keepGridInfo=cms.bool(False), ) from Configuration.Eras.Modifier_pA_2016_cff import pA_2016 pA_2016.toModify(hiFJGridEmptyAreaCalculator, doCentrality=False)
'keep *_ak4GenJetsNoNu_*_*', 'keep *_ak8GenJetsNoNu_*_*', 'keep *_genParticle_*_*' ) ) 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 #products from regular pp which does not fit the normal AOD for e in [pA_2016, peripheralPbPb, pp_on_XeXe_2017]: e.toModify( RecoJetsAOD.outputCommands, func=lambda outputCommands: outputCommands.extend(['keep *_towerMaker_*_*']) ) #HI-specific products: needed in AOD, propagate to more inclusive tiers as well for ec in [RecoJetsAOD.outputCommands, RecoJetsRECO.outputCommands, RecoJetsFEVT.outputCommands]: pA_2016.toModify( ec, func=lambda outputCommands: outputCommands.extend(['keep recoCentrality*_pACentrality_*_*', 'keep *_hiFJGridEmptyAreaCalculator_*_*', 'keep *_hiFJRhoProducer_*_*' ]) ) #HI-specific products: needed in AOD, propagate to more inclusive tiers as well for ec in [RecoJetsAOD.outputCommands, RecoJetsRECO.outputCommands, RecoJetsFEVT.outputCommands]: peripheralPbPb.toModify( ec, func=lambda outputCommands: outputCommands.extend(['keep recoCentrality*_pACentrality_*_*']) )
import FWCore.ParameterSet.Config as cms hiFJRhoProducer = cms.EDProducer('HiFJRhoProducer', jetSource = cms.InputTag('kt4PFJetsForRho'), nExcl = cms.int32(2), etaMaxExcl = cms.double(2.), ptMinExcl = cms.double(20.), nExcl2 = cms.int32(1), etaMaxExcl2 = cms.double(3.), ptMinExcl2 = cms.double(20.), etaRanges = cms.vdouble(-5., -3., -2.1, -1.3, 1.3, 2.1, 3., 5.) ) from Configuration.Eras.Modifier_pA_2016_cff import pA_2016 pA_2016.toModify(hiFJRhoProducer, etaRanges = cms.vdouble(-5., -3., -2., -1.5, -1., 1., 1.5, 2., 3., 5.))
"HLT_Mu17_TrkIsoVVL_v") #tightAnalyzer.probeParams = cms.PSet() hltMuonOfflineAnalyzers = cms.Sequence(globalAnalyzerTnP * trackerAnalyzerTnP * tightAnalyzerTnP * looseAnalyzerTnP * globalAnalyzerRef * trackerAnalyzerRef * tightAnalyzerRef * looseAnalyzerRef) from Configuration.Eras.Modifier_pA_2016_cff import pA_2016 for muAna in [ globalAnalyzerTnP.targetParams, trackerAnalyzerTnP.targetParams, tightAnalyzerTnP.targetParams, looseAnalyzerTnP.targetParams, globalAnalyzerRef.targetParams, trackerAnalyzerRef.targetParams, tightAnalyzerRef.targetParams, looseAnalyzerRef.targetParams, ]: pA_2016.toModify(muAna, ptCut_Jpsi=cms.untracked.double(5.0)) for muAna in [ globalAnalyzerTnP.binParams, trackerAnalyzerTnP.binParams, tightAnalyzerTnP.binParams, looseAnalyzerTnP.binParams, globalAnalyzerRef.binParams, trackerAnalyzerRef.binParams, tightAnalyzerRef.binParams, looseAnalyzerRef.binParams ]: pA_2016.toModify(muAna, ptCoarse=cms.untracked.vdouble(0., 1., 2., 3., 4., 5., 7., 9., 12., 15., 20., 30., 40.))
jetGlobalRecoTask = cms.Task(recoJetsTask, recoJetIdsTask, recoTrackJetsTask) jetGlobalReco = cms.Sequence(jetGlobalRecoTask) jetHighLevelRecoTask = cms.Task(recoPFJetsTask, jetCorrectorsForRecoTask, recoJetAssociationsTask, recoJetAssociationsExplicitTask, recoJPTJetsTask) jetHighLevelReco = cms.Sequence(jetHighLevelRecoTask) from RecoHI.HiJetAlgos.hiFJGridEmptyAreaCalculator_cff import hiFJGridEmptyAreaCalculator from RecoHI.HiJetAlgos.hiFJRhoProducer import hiFJRhoProducer from RecoHI.HiJetAlgos.HiRecoPFJets_cff import kt4PFJetsForRho from Configuration.Eras.Modifier_pA_2016_cff import pA_2016 from Configuration.Eras.Modifier_pp_on_AA_2018_cff import pp_on_AA_2018 from RecoHI.HiCentralityAlgos.pACentrality_cfi import pACentrality pA_2016.toModify(pACentrality, producePixelTracks=False) _jetHighLevelReco_pATask = jetHighLevelRecoTask.copy() _jetHighLevelReco_pATask.add(kt4PFJetsForRho) _jetHighLevelReco_pATask.add(hiFJRhoProducer) _jetHighLevelReco_pATask.add(hiFJGridEmptyAreaCalculator) _jetHighLevelReco_pATask.add(pACentrality) pA_2016.toReplaceWith(jetHighLevelRecoTask, _jetHighLevelReco_pATask) _jetGlobalReco_HI = cms.Sequence(recoJetsHI * recoJetIds) _jetHighLevelReco_HI = cms.Sequence(recoPFJetsHI * jetCorrectorsForReco * recoJetAssociations) pp_on_AA_2018.toReplaceWith(jetGlobalReco, _jetGlobalReco_HI) pp_on_AA_2018.toReplaceWith(jetHighLevelReco, _jetHighLevelReco_HI)
VertexLabel = cms.InputTag("offlineSlimmedPrimaryVertices"), BeamSpotLabel = cms.InputTag("offlineBeamSpot"), etaBin = cms.int32(40), etaBBin = cms.int32(40), etaEBin = cms.int32(40), etaBMin = cms.double(0.), etaBMax = cms.double(1.1), etaECMin = cms.double(0.9), etaECMax = cms.double(2.4), LowMassMin = cms.double(2.0), LowMassMax = cms.double(12.0), HighMassMin = cms.double(70.0), HighMassMax = cms.double(110.0), folder = cms.string("Muons_miniAOD/diMuonHistograms") ) from Configuration.Eras.Modifier_pA_2016_cff import pA_2016 pA_2016.toModify(diMuonHistos, etaBin = 350, etaBBin = 350, etaEBin = 350, LowMassMin = 2.0, LowMassMax = 51.0, HighMassMin = 55.0, HighMassMax = 125.0 )
import FWCore.ParameterSet.Config as cms hiFJRhoProducer = cms.EDProducer('HiFJRhoProducer', jetSource=cms.InputTag('kt4PFJetsForRho'), nExcl=cms.int32(2), etaMaxExcl=cms.double(2.), ptMinExcl=cms.double(20.), nExcl2=cms.int32(1), etaMaxExcl2=cms.double(3.), ptMinExcl2=cms.double(20.), etaRanges=cms.vdouble(-5., -3., -2.1, -1.3, 1.3, 2.1, 3., 5.)) from Configuration.Eras.Modifier_pA_2016_cff import pA_2016 pA_2016.toModify(hiFJRhoProducer, etaRanges=[-5., -3., -2., -1.5, -1., 1., 1.5, 2., 3., 5.])
from RecoJets.Configuration.JetIDProducers_cff import * from RecoJets.Configuration.RecoTrackJets_cff import * from RecoJets.Configuration.RecoJetAssociations_cff import * from RecoJets.Configuration.RecoPFJets_cff import * from RecoJets.Configuration.RecoJPTJets_cff import * from JetMETCorrections.Configuration.JetCorrectorsForReco_cff import * jetGlobalReco = cms.Sequence(recoJets*recoJetIds*recoTrackJets) jetHighLevelReco = cms.Sequence(recoPFJets*jetCorrectorsForReco*recoJetAssociations*recoJetAssociationsExplicit*recoJPTJets) from RecoHI.HiJetAlgos.hiFJGridEmptyAreaCalculator_cff import hiFJGridEmptyAreaCalculator from RecoHI.HiJetAlgos.hiFJRhoProducer import hiFJRhoProducer from RecoHI.HiJetAlgos.HiRecoPFJets_cff import kt4PFJetsForRho from Configuration.Eras.Modifier_pA_2016_cff import pA_2016 from Configuration.Eras.Modifier_pp_on_AA_2018_cff import pp_on_AA_2018 from RecoHI.HiCentralityAlgos.pACentrality_cfi import pACentrality pA_2016.toModify(pACentrality, producePixelTracks = False) _jetHighLevelReco_pA = jetHighLevelReco.copy() _jetHighLevelReco_pA += kt4PFJetsForRho _jetHighLevelReco_pA += hiFJRhoProducer _jetHighLevelReco_pA += hiFJGridEmptyAreaCalculator _jetHighLevelReco_pA += pACentrality pA_2016.toReplaceWith(jetHighLevelReco, _jetHighLevelReco_pA) _jetGlobalReco_HI = cms.Sequence(recoJetsHI*recoJetIds) _jetHighLevelReco_HI = cms.Sequence(recoPFJetsHI*jetCorrectorsForReco*recoJetAssociations) pp_on_AA_2018.toReplaceWith(jetGlobalReco,_jetGlobalReco_HI) pp_on_AA_2018.toReplaceWith(jetHighLevelReco,_jetHighLevelReco_HI)
from RecoJets.Configuration.RecoJetAssociations_cff import * from RecoJets.Configuration.RecoPFJets_cff import * from RecoJets.Configuration.RecoJPTJets_cff import * from JetMETCorrections.Configuration.JetCorrectorsForReco_cff import * jetGlobalReco = cms.Sequence(recoJets * recoJetIds * recoTrackJets) jetHighLevelReco = cms.Sequence(recoPFJets * jetCorrectorsForReco * recoJetAssociations * recoJetAssociationsExplicit * recoJPTJets) from Configuration.Eras.Modifier_pA_2016_cff import pA_2016 #HI-specific algorithms needed in pp scenario special configurations from RecoHI.HiJetAlgos.hiFJRhoProducer import hiFJRhoProducer from RecoHI.HiJetAlgos.hiFJGridEmptyAreaCalculator_cff import hiFJGridEmptyAreaCalculator pA_2016.toModify(hiFJGridEmptyAreaCalculator, doCentrality=False) kt4PFJetsForRho = kt4PFJets.clone(doAreaFastjet=True, jetPtMin=0.0, GhostArea=0.005) from RecoHI.HiCentralityAlgos.pACentrality_cfi import pACentrality pA_2016.toModify(pACentrality, producePixelTracks=False) _jetHighLevelReco_pA = jetHighLevelReco.copy() _jetHighLevelReco_pA += kt4PFJetsForRho _jetHighLevelReco_pA += hiFJRhoProducer _jetHighLevelReco_pA += hiFJGridEmptyAreaCalculator _jetHighLevelReco_pA += pACentrality pA_2016.toReplaceWith(jetHighLevelReco, _jetHighLevelReco_pA)