Making JPT jets from SisCone jets

JPT corrected jets can be made starting from IC5 jets or SisCone jets. To make the JPT jets from SisCone I edited two new configuation files in NtupleMakers /python.

ZSPJetCorrections219SisCone _cff.py

import FWCore.ParameterSet.Config as cms

# Jet corrections.
#
# Define the correction services for each algorithm.
#
ZSPJetCorrectorSisCone = cms.ESSource("ZSPJetCorrectionService",
    tagName = cms.string('ZSP_CMSSW219_Iterative_Cone_05'),
    label = cms.string('ZSPJetCorrectorSisCone')
#    tagName = cms.string('ZSP_CMSSW152_Iterative_Cone_05'),
)

#   
#   Define the producers of corrected jet collections for each algorithm.
#
ZSPJetCorJetSisCone = cms.EDProducer("CaloJetCorrectionProducer",
    src = cms.InputTag("sisCone5CaloJets"),
    correctors = cms.vstring('ZSPJetCorrectorSisCone'),
    alias = cms.untracked.string('ZSPJetCorSisCone')
)

#
#  Define a sequence to make all corrected jet collections at once.
#
ZSPJetCorrectionsSisCone = cms.Sequence(ZSPJetCorJetSisCone)

JetPlusTrackCorrectionsSisCone _cff.py

import FWCore.ParameterSet.Config as cms

# File: JetCorrections.cff
# Author: O. Kodolova
# Date: 1/24/07
#
# Jet corrections with tracks for the icone5 jets with ZSP corrections.
# ... adapted by DLE for SisCone jets ...
# 

from JetMETCorrections.Configuration.JetCorrectionsRecord_cfi import *
from RecoJets.Configuration.RecoJetAssociations_cff import *

JetPlusTrackZSPCorrectorSisCone = cms.ESSource("JetPlusTrackCorrectionService",
    JetTrackCollectionAtCalo = cms.InputTag("ZSPsisConeJetTracksAssociatorAtCaloFace"),
    respalgo = cms.int32(5),
    JetTrackCollectionAtVertex = cms.InputTag("ZSPsisConeJetTracksAssociatorAtVertex"),
    muonSrc = cms.InputTag("muons"),
    AddOutOfConeTracks = cms.bool(True),
    NonEfficiencyFile = cms.string('CMSSW_167_TrackNonEff'),
    NonEfficiencyFileResp = cms.string('CMSSW_167_TrackLeakage'),
    ResponseFile = cms.string('CMSSW_167_response'),
    label = cms.string('JetPlusTrackZSPCorrectorSisCone'),
    TrackQuality = cms.string('highPurity'),
    UseQuality = cms.bool(True)
)

JetPlusTrackZSPCorJetSisCone = cms.EDProducer("CaloJetCorrectionProducer",
    src = cms.InputTag("ZSPJetCorJetSisCone"),
    correctors = cms.vstring('JetPlusTrackZSPCorrectorSisCone'),
    alias = cms.untracked.string('JetPlusTrackZSPCorJetSisCone')
)

from RecoJets.JetAssociationProducers.iterativeCone5JTA_cff import*

ZSPsisConeJetTracksAssociatorAtVertex = iterativeCone5JetTracksAssociatorAtVertex.clone()
ZSPsisConeJetTracksAssociatorAtVertex.jets = cms.InputTag("ZSPJetCorJetSisCone")

ZSPsisConeJetTracksAssociatorAtCaloFace = iterativeCone5JetTracksAssociatorAtCaloFace.clone()
ZSPsisConeJetTracksAssociatorAtCaloFace.jets = cms.InputTag("ZSPJetCorJetSisCone")

ZSPsisConeJetExtender = iterativeCone5JetExtender.clone()
ZSPsisConeJetExtender.jets = cms.InputTag("ZSPJetCorJetSisCone")
ZSPsisConeJetExtender.jet2TracksAtCALO = cms.InputTag("ZSPsisConeJetTracksAssociatorAtCaloFace")
ZSPsisConeJetExtender.jet2TracksAtVX = cms.InputTag("ZSPsisConeJetTracksAssociatorAtVertex")


ZSPrecoJetAssociationsSisCone = cms.Sequence(ZSPsisConeJetTracksAssociatorAtVertex*ZSPsisConeJetTracksAssociatorAtCaloFace*ZSPsisConeJetExtender)

JetPlusTrackCorrectionsSisCone = cms.Sequence(ZSPrecoJetAssociationsSisCone*JetPlusTrackZSPCorJetSisCone)

Then in the cfg file

#-------------------------------------------------
# load JPT producer
#-------------------------------------------------
process.load("CMS2.NtupleMaker.JetPlusTrackCorrectionsSisCone_cff")
process.load("CMS2.NtupleMaker.ZSPJetCorrections219SisCone_cff")
#process.load("JetMETCorrections.Configuration.JetPlusTrackCorrections_cff")
#process.load("JetMETCorrections.Configuration.ZSPJetCorrections219_cff")

and change the sequence also

#process.JPTCorrection = cms.Sequence(process.ZSPJetCorrections*process.JetPlusTrackCorrections*process.jptMaker)
process.JPTCorrection = cms.Sequence(process.ZSPJetCorrectionsSisCone*process.JetPlusTrackCorrectionsSisCone*process.jptMaker)

Also need to change the python/jptMaker_cfi.py as appropriate for the jptInputTag and the calojetInputTag according to the new files above.

Validation

/store/mc/Summer08/ZeeJet_Pt80to120/GEN-SIM-RECO/IDEAL_V9_v1/0003/E071D34D-C2D1-DD11-915A-003048C185DC.root

Distribution of the ET of JPT corrected jets made from IC5 jets (black) compared with distribution of JPT corrected jets made from SisCone jets (red)

The ratio of the above with statistical errors

-- DaveEvans - 17 Apr 2009