JetValidationHeavyIons_cff.py

Go to the documentation of this file.

import FWCore.ParameterSet.Config as cms

### genjet cleaning for improved matching in HI environment

from RecoHI.HiJetAlgos.HiGenCleaner_cff import *

iterativeCone5HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = cms.InputTag('iterativeCone5HiGenJets'))
iterativeCone7HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = cms.InputTag('iterativeCone7HiGenJets'))
ak5HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = cms.InputTag('ak5HiGenJets'))
ak7HiCleanedGenJets = heavyIonCleanedGenJets.clone( src = cms.InputTag('ak7HiGenJets'))

### jet analyzer for two radii (0.5, 0.7) and three algorithms:
### iterative cone with PU, anti-kt with PU, anti-kt with fastjet PU


JetAnalyzerICPU5Calo = cms.EDAnalyzer("JetTester",
    JetType = cms.untracked.string('calo'),
    OutputFile = cms.untracked.string(''),
    src = cms.InputTag("iterativeConePu5CaloJets"),
    srcGen = cms.InputTag("iterativeCone5HiCleanedGenJets"),
    JetCorrections = cms.string(""),
    recoJetPtThreshold = cms.double(40),
    genEnergyFractionThreshold     = cms.double(0.05),
    matchGenPtThreshold                 = cms.double(20.0),
    RThreshold                     = cms.double(0.3)
)

JetAnalyzerICPU7Calo = cms.EDAnalyzer("JetTester",
    JetType = cms.untracked.string('calo'),
    OutputFile = cms.untracked.string(''),
    src = cms.InputTag("iterativeConePu7CaloJets"),
    srcGen = cms.InputTag("iterativeCone7HiCleanedGenJets"), 
    JetCorrections = cms.string(""),
    recoJetPtThreshold = cms.double(40),
    genEnergyFractionThreshold     = cms.double(0.05),
    matchGenPtThreshold                 = cms.double(20.0),
    RThreshold                     = cms.double(0.3) 
)

JetAnalyzerAkPU5Calo = cms.EDAnalyzer("JetTester",
    JetType = cms.untracked.string('calo'),
    OutputFile = cms.untracked.string(''),
    src = cms.InputTag("akPu5CaloJets"),
    srcGen = cms.InputTag("ak5HiCleanedGenJets"),  
    JetCorrections = cms.string(""),
    recoJetPtThreshold = cms.double(40),
    genEnergyFractionThreshold     = cms.double(0.05),
    matchGenPtThreshold                 = cms.double(20.0),
    RThreshold                     = cms.double(0.3) 
)

JetAnalyzerAkPU7Calo = cms.EDAnalyzer("JetTester",
    JetType = cms.untracked.string('calo'),
    OutputFile = cms.untracked.string(''),
    src = cms.InputTag("akPu7CaloJets"),
    srcGen = cms.InputTag("ak7HiCleanedGenJets"),       
    JetCorrections = cms.string(""),
    recoJetPtThreshold = cms.double(40),
    genEnergyFractionThreshold     = cms.double(0.05),
    matchGenPtThreshold                 = cms.double(20.0),
    RThreshold                     = cms.double(0.3) 
)

JetAnalyzerAkFastPU5Calo = cms.EDAnalyzer("JetTester",
    JetType = cms.untracked.string('calo'),
    OutputFile = cms.untracked.string(''),
    src = cms.InputTag("akFastPu5CaloJets"),
    srcGen = cms.InputTag("ak5HiCleanedGenJets"),  
    JetCorrections = cms.string(""),
    recoJetPtThreshold = cms.double(40),
    genEnergyFractionThreshold     = cms.double(0.05),
    matchGenPtThreshold                 = cms.double(20.0),
    RThreshold                     = cms.double(0.3) 
)

JetAnalyzerAkFastPU7Calo = cms.EDAnalyzer("JetTester",
    JetType = cms.untracked.string('calo'),
    OutputFile = cms.untracked.string(''),
    src = cms.InputTag("akFastPu7CaloJets"),
    srcGen = cms.InputTag("ak7HiCleanedGenJets"),       
    JetCorrections = cms.string(""),
    recoJetPtThreshold = cms.double(40),
    genEnergyFractionThreshold     = cms.double(0.05),
    matchGenPtThreshold                 = cms.double(20.0),
    RThreshold                     = cms.double(0.3) 
)

hiJetValidation = cms.Sequence(
    (iterativeCone5HiCleanedGenJets * JetAnalyzerICPU5Calo) 
    #+ (iterativeCone7HiCleanedGenJets * JetAnalyzerICPU7Calo)
    #+ (ak5HiCleanedGenJets * JetAnalyzerAkPU5Calo * JetAnalyzerAkFastPU5Calo
    #+ (ak7HiCleanedGenJets*JetAnalyzerAkPU7Calo *JetAnalyzerAkFastPU7Calo)
    )