Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 13:34:59

0001 #
0002 # This configuration file collects the data needed to calibrate
0003 # FFTJet pule-up estimator for PFJets
0004 #
0005 # Replace strings INPUTFILE and OUTPUTFILE in the configuration below
0006 #
0007 # I. Volobouev, April 27 2011
0008 #
0009 import FWCore.ParameterSet.Config as cms
0010 
0011 process = cms.Process("FFTJetTest")
0012 
0013 process.load("FWCore.MessageService.MessageLogger_cfi")
0014 process.MessageLogger.cerr.FwkReport.reportEvery = 1000
0015 
0016 process.TFileService = cms.Service(
0017     "TFileService",
0018     fileName = cms.string("OUTPUTFILE")
0019 )
0020 
0021 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )
0022 
0023 process.source = cms.Source("PoolSource",
0024     fileNames = cms.untracked.vstring(
0025         'file:INPUTFILE'
0026     )
0027 )
0028 
0029 from RecoJets.FFTJetProducers.fftjetpileupprocessor_pfprod_cfi import *
0030 from RecoJets.FFTJetProducers.fftjetpileupestimator_pf_cfi import *
0031 from RecoJets.JetAnalyzers.fftjetpileupanalyzer_cfi import *
0032 
0033 fftjet_pileup_analyzer.summaryLabel = cms.InputTag(
0034     "pileupestimator", "FFTJetPileupEstimatePF")
0035 fftjet_pileup_analyzer.fastJetRhoLabel = cms.InputTag("mykt6PFJets", "rho")
0036 fftjet_pileup_analyzer.fastJetSigmaLabel = cms.InputTag("mykt6PFJets", "sigma")
0037 fftjet_pileup_analyzer.collectFastJetRho = cms.bool(True)
0038 
0039 process.pileupprocessor = fftjet_pileup_processor_pf
0040 process.pileupestimator = fftjet_pileup_estimator_pf
0041 process.pileupanalyzer = fftjet_pileup_analyzer
0042 
0043 # Configure FastJet rho reconstruction
0044 from RecoJets.JetProducers.kt4PFJets_cfi import kt4PFJets
0045 mykt6PFJets = kt4PFJets.clone( rParam = 0.6 )
0046 mykt6PFJets.doRhoFastjet = True
0047 mykt6PFJets.doAreaFastjet = True
0048 mykt6PFJets.voronoiRfact = 0.9
0049 process.mykt6PFJets = mykt6PFJets
0050 
0051 process.p = cms.Path(process.mykt6PFJets*process.pileupprocessor*process.pileupestimator*process.pileupanalyzer)