File indexing completed on 2024-04-06 12:25:19
0001 import FWCore.ParameterSet.Config as cms
0002
0003 from RecoJets.Configuration.CaloTowersRec_cff import *
0004
0005
0006 from RecoJets.JetProducers.AnomalousCellParameters_cfi import *
0007 from RecoHI.HiJetAlgos.HiCaloJetParameters_cff import *
0008
0009
0010 CaloTowerConstituentsMapBuilder = cms.ESProducer("CaloTowerConstituentsMapBuilder",
0011 MapFile = cms.untracked.string('Geometry/CaloTopology/data/CaloTowerEEGeometric.map.gz')
0012 )
0013
0014 caloTowers = cms.EDProducer("CaloTowerCandidateCreator",
0015 src = cms.InputTag("towerMaker"),
0016 e = cms.double(0.0),
0017 verbose = cms.untracked.int32(0),
0018 pt = cms.double(0.0),
0019 minimumE = cms.double(0.0),
0020 minimumEt = cms.double(0.0),
0021 et = cms.double(0.0)
0022 )
0023
0024
0025
0026
0027 iterativeConePu5CaloJets = cms.EDProducer(
0028 "FastjetJetProducer",
0029 HiCaloJetParameters,
0030 AnomalousCellParameters,
0031 MultipleAlgoIteratorBlock,
0032 jetAlgorithm = cms.string("IterativeCone"),
0033 rParam = cms.double(0.5)
0034 )
0035 iterativeConePu5CaloJets.radiusPU = 0.5
0036
0037
0038 ktPu4CaloJets = cms.EDProducer(
0039 "FastjetJetProducer",
0040 HiCaloJetParameters,
0041 AnomalousCellParameters,
0042 MultipleAlgoIteratorBlock,
0043 jetAlgorithm = cms.string("Kt"),
0044 rParam = cms.double(0.4)
0045 )
0046 ktPu4CaloJets.radiusPU = 0.5
0047
0048 ktPu6CaloJets = cms.EDProducer(
0049 "FastjetJetProducer",
0050 HiCaloJetParameters,
0051 AnomalousCellParameters,
0052 MultipleAlgoIteratorBlock,
0053 jetAlgorithm = cms.string("Kt"),
0054 rParam = cms.double(0.6)
0055 )
0056 ktPu6CaloJets.radiusPU = 0.7
0057
0058
0059 akPu5CaloJets = cms.EDProducer(
0060 "FastjetJetProducer",
0061 HiCaloJetParameters,
0062 AnomalousCellParameters,
0063 MultipleAlgoIteratorBlock,
0064 jetAlgorithm = cms.string("AntiKt"),
0065 rParam = cms.double(0.5)
0066 )
0067 akPu5CaloJets.radiusPU = 0.5
0068 akPu5CaloJets.puPtMin = 10
0069
0070 akPu1CaloJets = akPu5CaloJets.clone(rParam = 0.1, puPtMin = 4)
0071 akPu2CaloJets = akPu5CaloJets.clone(rParam = 0.2, puPtMin = 4)
0072 akPu3CaloJets = akPu5CaloJets.clone(rParam = 0.3, puPtMin = 6)
0073 akPu4CaloJets = akPu5CaloJets.clone(rParam = 0.4, puPtMin = 8)
0074 akPu6CaloJets = akPu5CaloJets.clone(rParam = 0.6, puPtMin = 12)
0075 akPu7CaloJets = akPu5CaloJets.clone(rParam = 0.7, puPtMin = 14)
0076
0077 ak5CaloJets = cms.EDProducer(
0078 "FastjetJetProducer",
0079 HiCaloJetParameters,
0080 AnomalousCellParameters,
0081 MultipleAlgoIteratorBlock,
0082 jetAlgorithm = cms.string("AntiKt"),
0083 rParam = cms.double(0.5)
0084 )
0085 ak5CaloJets.doPUOffsetCorr = False
0086
0087 ak1CaloJets = ak5CaloJets.clone(rParam = 0.1)
0088 ak2CaloJets = ak5CaloJets.clone(rParam = 0.2)
0089 ak3CaloJets = ak5CaloJets.clone(rParam = 0.3)
0090 ak4CaloJets = ak5CaloJets.clone(rParam = 0.4)
0091 ak6CaloJets = ak5CaloJets.clone(rParam = 0.6)
0092 ak7CaloJets = ak5CaloJets.clone(rParam = 0.7)
0093
0094
0095
0096 hiRecoJetsTask = cms.Task(
0097 caloTowersRecTask,caloTowers,
0098 iterativeConePu5CaloJets,
0099 akPu3CaloJets,akPu4CaloJets,akPu5CaloJets
0100 )
0101 hiRecoJets = cms.Sequence(hiRecoJetsTask)
0102
0103
0104 hiRecoAllJetsTask = cms.Task(
0105 caloTowersRecTask,caloTowers,iterativeConePu5CaloJets
0106 ,ak1CaloJets,ak2CaloJets,ak3CaloJets,ak4CaloJets,ak5CaloJets,ak6CaloJets,ak7CaloJets
0107 ,akPu1CaloJets,akPu2CaloJets,akPu3CaloJets,akPu4CaloJets,akPu5CaloJets,akPu6CaloJets,akPu7CaloJets,
0108 ktPu4CaloJets,ktPu6CaloJets
0109 )
0110 hiRecoAllJets = cms.Sequence(hiRecoAllJetsTask)