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.])
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
)
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)