| 1 |
|
import FWCore.ParameterSet.Config as cms |
| 2 |
|
|
| 3 |
< |
### this needs some genjet cleaning or improved matching criteria |
| 4 |
< |
### to do proper matching in the heavy-ion environment |
| 3 |
> |
### genjet cleaning for improved matching in HI environment |
| 4 |
> |
|
| 5 |
> |
iterativeCone5HiCleanedGenJets = cms.EDProducer('HiGenJetCleaner', |
| 6 |
> |
src = cms.untracked.string('iterativeCone5HiGenJets'), |
| 7 |
> |
deltaR = cms.untracked.double(0.25), |
| 8 |
> |
ptCut = cms.untracked.double(20), |
| 9 |
> |
createNewCollection = cms.untracked.bool(True), |
| 10 |
> |
fillDummyEntries = cms.untracked.bool(True) |
| 11 |
> |
) |
| 12 |
> |
|
| 13 |
> |
iterativeCone7HiCleanedGenJets = cms.EDProducer('HiGenJetCleaner', |
| 14 |
> |
src = cms.untracked.string('iterativeCone7HiGenJets'), |
| 15 |
> |
deltaR = cms.untracked.double(0.25), |
| 16 |
> |
ptCut = cms.untracked.double(20), |
| 17 |
> |
createNewCollection = cms.untracked.bool(True), |
| 18 |
> |
fillDummyEntries = cms.untracked.bool(True) |
| 19 |
> |
) |
| 20 |
> |
|
| 21 |
> |
ak5HiCleanedGenJets = cms.EDProducer('HiGenJetCleaner', |
| 22 |
> |
src = cms.untracked.string('ak5HiGenJets'), |
| 23 |
> |
deltaR = cms.untracked.double(0.25), |
| 24 |
> |
ptCut = cms.untracked.double(20), |
| 25 |
> |
createNewCollection = cms.untracked.bool(True), |
| 26 |
> |
fillDummyEntries = cms.untracked.bool(True) |
| 27 |
> |
) |
| 28 |
> |
|
| 29 |
> |
ak7HiCleanedGenJets = cms.EDProducer('HiGenJetCleaner', |
| 30 |
> |
src = cms.untracked.string('ak7HiGenJets'), |
| 31 |
> |
deltaR = cms.untracked.double(0.25), |
| 32 |
> |
ptCut = cms.untracked.double(20), |
| 33 |
> |
createNewCollection = cms.untracked.bool(True), |
| 34 |
> |
fillDummyEntries = cms.untracked.bool(True) |
| 35 |
> |
) |
| 36 |
> |
|
| 37 |
> |
### jet analyzer for two radii (0.5, 0.7) and three algorithms: |
| 38 |
> |
### iterative cone with PU, anti-kt with PU, anti-kt with fastjet PU |
| 39 |
|
|
| 40 |
|
JetAnalyzerICPU5Calo = cms.EDAnalyzer("CaloJetTester", |
| 41 |
|
src = cms.InputTag("iterativeConePu5CaloJets"), |
| 42 |
< |
srcGen = cms.InputTag("iterativeCone5HiGenJets"), |
| 42 |
> |
srcGen = cms.InputTag("iterativeCone5HiCleanedGenJets"), |
| 43 |
|
genEnergyFractionThreshold = cms.double(0.05), |
| 44 |
|
genPtThreshold = cms.double(1.0), |
| 45 |
|
RThreshold = cms.double(0.3), |
| 48 |
|
|
| 49 |
|
JetAnalyzerICPU7Calo = cms.EDAnalyzer("CaloJetTester", |
| 50 |
|
src = cms.InputTag("iterativeConePu7CaloJets"), |
| 51 |
< |
srcGen = cms.InputTag("iterativeCone7HiGenJets"), |
| 51 |
> |
srcGen = cms.InputTag("iterativeCone7HiCleanedGenJets"), |
| 52 |
|
genEnergyFractionThreshold = cms.double(0.05), |
| 53 |
|
genPtThreshold = cms.double(1.0), |
| 54 |
|
RThreshold = cms.double(0.3), |
| 57 |
|
|
| 58 |
|
JetAnalyzerAkPU5Calo = cms.EDAnalyzer("CaloJetTester", |
| 59 |
|
src = cms.InputTag("akPu5CaloJets"), |
| 60 |
< |
srcGen = cms.InputTag("ak5HiGenJets"), |
| 60 |
> |
srcGen = cms.InputTag("ak5HiCleanedGenJets"), |
| 61 |
|
genEnergyFractionThreshold = cms.double(0.05), |
| 62 |
|
genPtThreshold = cms.double(1.0), |
| 63 |
|
RThreshold = cms.double(0.3), |
| 66 |
|
|
| 67 |
|
JetAnalyzerAkPU7Calo = cms.EDAnalyzer("CaloJetTester", |
| 68 |
|
src = cms.InputTag("akPu7CaloJets"), |
| 69 |
< |
srcGen = cms.InputTag("ak7HiGenJets"), |
| 69 |
> |
srcGen = cms.InputTag("ak7HiCleanedGenJets"), |
| 70 |
|
genEnergyFractionThreshold = cms.double(0.05), |
| 71 |
|
genPtThreshold = cms.double(1.0), |
| 72 |
|
RThreshold = cms.double(0.3), |
| 75 |
|
|
| 76 |
|
JetAnalyzerAkFastPU5Calo = cms.EDAnalyzer("CaloJetTester", |
| 77 |
|
src = cms.InputTag("akFastPu5CaloJets"), |
| 78 |
< |
srcGen = cms.InputTag("ak5HiGenJets"), |
| 78 |
> |
srcGen = cms.InputTag("ak5HiCleanedGenJets"), |
| 79 |
|
genEnergyFractionThreshold = cms.double(0.05), |
| 80 |
|
genPtThreshold = cms.double(1.0), |
| 81 |
|
RThreshold = cms.double(0.3), |
| 84 |
|
|
| 85 |
|
JetAnalyzerAkFastPU7Calo = cms.EDAnalyzer("CaloJetTester", |
| 86 |
|
src = cms.InputTag("akFastPu7CaloJets"), |
| 87 |
< |
srcGen = cms.InputTag("ak7HiGenJets"), |
| 87 |
> |
srcGen = cms.InputTag("ak7HiCleanedGenJets"), |
| 88 |
|
genEnergyFractionThreshold = cms.double(0.05), |
| 89 |
|
genPtThreshold = cms.double(1.0), |
| 90 |
|
RThreshold = cms.double(0.3), |
| 91 |
|
reverseEnergyFractionThreshold = cms.double(0.5) |
| 92 |
|
) |
| 93 |
|
|
| 94 |
< |
hiJetValidation = cms.Sequence(JetAnalyzerICPU5Calo |
| 95 |
< |
#*JetAnalyzerICPU7Calo |
| 96 |
< |
#*JetAnalyzerAkPU5Calo |
| 97 |
< |
#*JetAnalyzerAkPU7Calo |
| 98 |
< |
#*JetAnalyzerAkFastPU5Calo |
| 99 |
< |
#*JetAnalyzerAkFastPU7Calo |
| 66 |
< |
) |
| 94 |
> |
hiJetValidation = cms.Sequence( |
| 95 |
> |
(iterativeCone5HiCleanedGenJets * JetAnalyzerICPU5Calo) |
| 96 |
> |
#+ (iterativeCone7HiCleanedGenJets * JetAnalyzerICPU7Calo) |
| 97 |
> |
#+ (ak5HiCleanedGenJets * JetAnalyzerAkPU5Calo * JetAnalyzerAkFastPU5Calo |
| 98 |
> |
#+ (ak7HiCleanedGenJets*JetAnalyzerAkPU7Calo *JetAnalyzerAkFastPU7Calo) |
| 99 |
> |
) |
