CMS 3D CMS Logo

/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_14/src/RecoJets/Configuration/python/RecoJets_cff.py

Go to the documentation of this file.
00001 import FWCore.ParameterSet.Config as cms
00002 
00003 from RecoJets.JetProducers.sc5CaloJets_cfi import sisCone5CaloJets
00004 from RecoJets.JetProducers.ic5CaloJets_cfi import iterativeCone5CaloJets
00005 from RecoJets.JetProducers.ak5CaloJets_cfi import ak5CaloJets
00006 from RecoJets.JetProducers.gk5CaloJets_cfi import gk5CaloJets
00007 from RecoJets.JetProducers.kt4CaloJets_cfi import kt4CaloJets
00008 from RecoJets.JetProducers.ca4CaloJets_cfi import ca4CaloJets
00009 
00010 
00011 sisCone7CaloJets = sisCone5CaloJets.clone( rParam = 0.7 )
00012 ak7CaloJets = ak5CaloJets.clone( rParam = 0.7 )
00013 gk7CaloJets = gk5CaloJets.clone( rParam = 0.7 )
00014 kt6CaloJets = kt4CaloJets.clone( rParam = 0.6 )
00015 ca6CaloJets = ca4CaloJets.clone( rParam = 0.6 )
00016 
00017 
00018 # Restrict SISCone algorithm to 1000 towers input
00019 sisCone5CaloJets.restrictInputs = cms.bool(True)
00020 sisCone5CaloJets.maxInputs = cms.uint32(1000)
00021 
00022 sisCone7CaloJets.restrictInputs = cms.bool(True)
00023 sisCone7CaloJets.maxInputs = cms.uint32(1000)
00024 
00025 doPileup = cms.bool(True)
00026 
00027 sisCone5CaloJetsPUCorr      =sisCone5CaloJets.clone      (doPUOffsetCorr = doPileup)
00028 sisCone7CaloJetsPUCorr      =sisCone7CaloJets.clone      (doPUOffsetCorr = doPileup)
00029 kt4CaloJetsPUCorr           =kt4CaloJets.clone           (doPUOffsetCorr = doPileup)
00030 kt6CaloJetsPUCorr           =kt6CaloJets.clone           (doPUOffsetCorr = doPileup)
00031 iterativeCone5CaloJetsPUCorr=iterativeCone5CaloJets.clone(doPUOffsetCorr = doPileup)
00032 ak5CaloJetsPUCorr           =ak5CaloJets.clone           (doPUOffsetCorr = doPileup)
00033 ak7CaloJetsPUCorr           =ak7CaloJets.clone           (doPUOffsetCorr = doPileup)
00034 gk5CaloJetsPUCorr           =gk5CaloJets.clone           (doPUOffsetCorr = doPileup)
00035 gk7CaloJetsPUCorr           =gk7CaloJets.clone           (doPUOffsetCorr = doPileup)
00036 ca4CaloJetsPUCorr           =ca4CaloJets.clone           (doPUOffsetCorr = doPileup)
00037 ca6CaloJetsPUCorr           =ca6CaloJets.clone           (doPUOffsetCorr = doPileup)
00038 
00039 #compute areas for Fastjet PU subtraction  
00040 kt6CaloJets.doRhoFastjet = True
00041 kt6CaloJets.doAreaFastjet = True
00042 #use active areas and not Voronoi tessellation for the moment
00043 #kt6CaloJets.voronoiRfact = 0.9
00044 ak5CaloJets.doAreaFastjet = True
00045 ak7CaloJets.doAreaFastjet = True
00046 
00047 
00048 kt6CaloJetsCentral = kt6CaloJets.clone(
00049     Ghost_EtaMax = cms.double(3.1),
00050     Rho_EtaMax = cms.double(2.5)
00051     )
00052 
00053 kt6CaloJetsCentralPUCorr           =kt6CaloJetsCentral.clone           (doPUOffsetCorr = doPileup)
00054 
00055 
00056 RecoJets   =cms.Sequence(kt4CaloJets+kt6CaloJets+kt6CaloJetsCentral+
00057                          iterativeCone5CaloJets+
00058                          ak5CaloJets+ak7CaloJets)
00059 
00060 recoAllJets=cms.Sequence(sisCone5CaloJets+sisCone7CaloJets+
00061                          kt4CaloJets+kt6CaloJets+kt6CaloJetsCentral+
00062                          iterativeCone5CaloJets+
00063                          ak5CaloJets+ak7CaloJets+
00064                          gk5CaloJets+gk7CaloJets+
00065                          ca4CaloJets+ca6CaloJets)
00066 
00067 
00068 recoAllJetsPUOffsetCorr=cms.Sequence(sisCone5CaloJetsPUCorr+sisCone7CaloJetsPUCorr+
00069                                      kt4CaloJetsPUCorr+kt6CaloJetsPUCorr+kt6CaloJetsCentralPUCorr+
00070                                      iterativeCone5CaloJetsPUCorr+
00071                                      ak5CaloJetsPUCorr+ak7CaloJetsPUCorr+
00072                                      gk5CaloJetsPUCorr+gk7CaloJetsPUCorr+
00073                                      ca4CaloJetsPUCorr+ca6CaloJetsPUCorr)