CMS 3D CMS Logo

/data/refman/pasoursint/CMSSW_5_3_3/src/RecoJets/Configuration/python/RecoPFJets_cff.py

Go to the documentation of this file.
00001 import FWCore.ParameterSet.Config as cms
00002 
00003 from RecoJets.JetProducers.sc5PFJets_cfi import sisCone5PFJets
00004 from RecoJets.JetProducers.ic5PFJets_cfi import iterativeCone5PFJets
00005 from RecoJets.JetProducers.ak5PFJets_cfi import ak5PFJets
00006 from RecoJets.JetProducers.ak5PFJetsTrimmed_cfi import ak5PFJetsTrimmed
00007 from RecoJets.JetProducers.ak5PFJetsFiltered_cfi import ak5PFJetsFiltered, ak5PFJetsMassDropFiltered
00008 from RecoJets.JetProducers.ak5PFJetsPruned_cfi import ak5PFJetsPruned
00009 from RecoJets.JetProducers.gk5PFJets_cfi import gk5PFJets
00010 from RecoJets.JetProducers.kt4PFJets_cfi import kt4PFJets
00011 from RecoJets.JetProducers.ca4PFJets_cfi import ca4PFJets
00012 from RecoJets.JetProducers.fixedGridRhoProducer_cfi import fixedGridRhoAll
00013 from RecoJets.JetProducers.fixedGridRhoProducerFastjet_cfi import fixedGridRhoFastjetAll
00014 
00015 sisCone7PFJets = sisCone5PFJets.clone( rParam = 0.7 )
00016 ak7PFJets = ak5PFJets.clone( rParam = 0.7 )
00017 gk7PFJets = gk5PFJets.clone( rParam = 0.7 )
00018 kt6PFJets = kt4PFJets.clone( rParam = 0.6 )
00019 ca6PFJets = ca4PFJets.clone( rParam = 0.6 )
00020 
00021 #compute areas for Fastjet PU subtraction  
00022 kt6PFJets.doRhoFastjet = True
00023 kt6PFJets.doAreaFastjet = True
00024 #use active areas and not Voronoi tessellation for the moment
00025 kt6PFJets.voronoiRfact = 0.9
00026 ak5PFJets.doAreaFastjet = True
00027 ak5PFJetsTrimmed.doAreaFastjet = True
00028 ak7PFJets.doAreaFastjet = True
00029 
00030 
00031 
00032 
00033 kt6PFJetsCentralChargedPileUp = kt6PFJets.clone(
00034     src = cms.InputTag("pfPileUpAllChargedParticles"),
00035     Ghost_EtaMax = cms.double(3.1),
00036     Rho_EtaMax = cms.double(2.5)
00037     )
00038 
00039 kt6PFJetsCentralNeutral = kt6PFJets.clone(
00040     src = cms.InputTag("pfAllNeutralHadronsAndPhotons"),
00041     Ghost_EtaMax = cms.double(3.1),
00042     Rho_EtaMax = cms.double(2.5),
00043     inputEtMin = cms.double(0.5)
00044     )
00045 
00046 
00047 
00048 kt6PFJetsCentralNeutralTight = kt6PFJetsCentralNeutral.clone(
00049     inputEtMin = cms.double(1.0)
00050     )
00051 
00052 
00053 recoPFJets   =cms.Sequence(kt4PFJets+kt6PFJets+
00054                            kt6PFJetsCentralChargedPileUp+
00055                            kt6PFJetsCentralNeutral+
00056                            kt6PFJetsCentralNeutralTight+
00057                            fixedGridRhoAll+
00058                            fixedGridRhoFastjetAll+
00059                            iterativeCone5PFJets+
00060                            ak5PFJets+ak7PFJets)
00061 
00062 recoAllPFJets=cms.Sequence(sisCone5PFJets+sisCone7PFJets+
00063                            kt4PFJets+kt6PFJets+
00064                            kt6PFJetsCentralChargedPileUp+
00065                            kt6PFJetsCentralNeutral+
00066                            kt6PFJetsCentralNeutralTight+
00067                            fixedGridRhoAll+
00068                            fixedGridRhoFastjetAll+
00069                            iterativeCone5PFJets+
00070                            ak5PFJets+ak7PFJets+
00071                            ak5PFJetsTrimmed+ak5PFJetsPruned+
00072                            ak5PFJetsFiltered+ak5PFJetsMassDropFiltered+
00073                            gk5PFJets+gk7PFJets+
00074                            ca4PFJets+ca6PFJets)