HiJetAlgos/python/HiRecoJets_cff.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 from RecoJets.Configuration.CaloTowersRec_cff import *
0004
0005 ## Default Parameter Sets
0006 from RecoJets.JetProducers.AnomalousCellParameters_cfi import *
0007 from RecoHI.HiJetAlgos.HiCaloJetParameters_cff import *
0008
0009 ## Calo Towers
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 ## Noise reducing PU subtraction algos
0025
0026 ## Iterative Cone
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 ## kT
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 ## anti-kT
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 ## Default Sequence
0096 hiRecoJetsTask = cms.Task(
0097     caloTowersRecTask,caloTowers,
0098     iterativeConePu5CaloJets,
0099     akPu3CaloJets,akPu4CaloJets,akPu5CaloJets
0100     )
0101 hiRecoJets = cms.Sequence(hiRecoJetsTask)
0102
0103 ## Extended Sequence
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)