HiJetAlgos/python/HiRecoPFJets_cff.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 ## Default Parameter Sets
0004 from RecoJets.JetProducers.AnomalousCellParameters_cfi import *
0005 from RecoHI.HiJetAlgos.HiPFJetParameters_cff import *
0006
0007 #pseudo towers for noise suppression background subtraction
0008 from RecoHI.HiJetAlgos.PFTowers_cfi import PFTowers
0009
0010 #dummy sequence to speed-up reconstruction in pp_on_AA era
0011 pfEmptyCollection = cms.EDFilter('GenericPFCandidateSelector',
0012                                  src = cms.InputTag('particleFlow'),
0013                                  cut = cms.string("pt<0")
0014                                 )
0015
0016 ak5PFJets = cms.EDProducer(
0017     "FastjetJetProducer",
0018     HiPFJetParameters,
0019     AnomalousCellParameters,
0020     MultipleAlgoIteratorBlock,
0021     jetAlgorithm = cms.string("AntiKt"),
0022     rParam       = cms.double(0.5)
0023 )
0024 ak5PFJets.src = 'particleFlow'
0025
0026 akPu5PFJets = ak5PFJets.clone(
0027     jetType        = 'BasicJet',
0028     doPVCorrection = False,
0029     doPUOffsetCorr = True,
0030     subtractorName = "MultipleAlgoIterator",
0031     src            = 'PFTowers',
0032     doAreaFastjet  = False,
0033     puPtMin        = cms.double(25)
0034 )
0035
0036 akPu1PFJets = akPu5PFJets.clone(rParam = 0.1, puPtMin = 10)
0037 akPu2PFJets = akPu5PFJets.clone(rParam = 0.2, puPtMin = 10)
0038 akPu3PFJets = akPu5PFJets.clone(rParam = 0.3, puPtMin = 15)
0039 akPu4PFJets = akPu5PFJets.clone(rParam = 0.4, puPtMin = 20)
0040 akPu6PFJets = akPu5PFJets.clone(rParam = 0.6, puPtMin = 30)
0041 akPu7PFJets = akPu5PFJets.clone(rParam = 0.7, puPtMin = 35)
0042
0043 hiPFCandCleanerforJets = cms.EDFilter('GenericPFCandidateSelector',
0044                                 src = cms.InputTag('particleFlow'),
0045                                 cut = cms.string("pt>5 && abs(eta)< 2")
0046                                 )
0047
0048 ak4PFJetsForFlow = akPu5PFJets.clone(
0049    Ghost_EtaMax = 5.0,
0050    Rho_EtaMax   = 4.4,
0051    doRhoFastjet = False,
0052    jetPtMin     = 15.0,
0053    nSigmaPU     = 1.0,
0054    rParam       = 0.4,
0055    radiusPU     = 0.5,
0056    src          = "hiPFCandCleanerforJets",
0057 )
0058
0059 kt4PFJetsForRho = cms.EDProducer(
0060     "FastjetJetProducer",
0061     HiPFJetParameters,
0062     AnomalousCellParameters,
0063     jetAlgorithm = cms.string("Kt"),
0064     rParam       = cms.double(0.4)
0065 )
0066 kt4PFJetsForRho.src           = 'particleFlow'
0067 kt4PFJetsForRho.doAreaFastjet = True
0068 kt4PFJetsForRho.jetPtMin      = 0.0
0069 kt4PFJetsForRho.GhostArea     = 0.005
0070
0071 from RecoHI.HiJetAlgos.hiFJRhoProducer import hiFJRhoProducer
0072
0073 import RecoHI.HiJetAlgos.hiFJRhoFlowModulationProducer_cfi as _mod
0074 hiFJRhoFlowModulation = _mod.hiFJRhoFlowModulationProducer.clone()
0075
0076 import RecoHI.HiJetAlgos.hiPuRhoProducer_cfi as _mod
0077 hiPuRho = _mod.hiPuRhoProducer.clone()
0078
0079 akCs4PFJets = cms.EDProducer(
0080     "CSJetProducer",
0081     HiPFJetParameters,
0082     AnomalousCellParameters,
0083     jetAlgorithm  = cms.string("AntiKt"),
0084     rParam        = cms.double(0.4),
0085     etaMap = cms.InputTag('hiPuRho', 'mapEtaEdges'),
0086     rho = cms.InputTag('hiPuRho', 'mapToRho'),
0087     rhom = cms.InputTag('hiPuRho', 'mapToRhoM'),
0088     csRParam  = cms.double(-1.),
0089     csAlpha   = cms.double(2.),
0090     writeJetsWithConst = cms.bool(True),
0091     useModulatedRho = cms.bool(False),
0092     rhoFlowFitParams = cms.InputTag('hiFJRhoFlowModulation', 'rhoFlowFitParams'),
0093     jetCollInstanceName = cms.string("pfParticlesCs"),
0094 )
0095 akCs4PFJets.src               = 'particleFlow'
0096 akCs4PFJets.doAreaFastjet     = True
0097 akCs4PFJets.jetPtMin          = 0.0
0098 akCs4PFJets.useExplicitGhosts = cms.bool(True)
0099 akCs4PFJets.GhostArea         = 0.005
0100
0101 akCs3PFJets = akCs4PFJets.clone(rParam = 0.3)
0102
0103 hiRecoPFJetsTask = cms.Task(
0104                            PFTowers,
0105                            akPu3PFJets,
0106                            akPu4PFJets,
0107                            akPu5PFJets,
0108                            hiPFCandCleanerforJets,
0109                            kt4PFJetsForRho,
0110                            ak4PFJetsForFlow,
0111                            hiFJRhoProducer,
0112                            hiPuRho,
0113                            hiFJRhoFlowModulation,
0114                            akCs3PFJets,
0115                            akCs4PFJets
0116     )
0117 hiRecoPFJets = cms.Sequence(hiRecoPFJetsTask)
0118
0119 from Configuration.ProcessModifiers.run2_miniAOD_pp_on_AA_103X_cff import run2_miniAOD_pp_on_AA_103X
0120 run2_miniAOD_pp_on_AA_103X.toModify(akCs4PFJets,src = 'cleanedParticleFlow')
0121 run2_miniAOD_pp_on_AA_103X.toModify(PFTowers,src = 'cleanedParticleFlow')
0122 run2_miniAOD_pp_on_AA_103X.toModify(hiPuRho,src = "PackedPFTowers")