Back to home page

Project CMSSW displayed by LXR

 
 

    


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

0001 # PYTHON configuration file.
0002 # Description:  Example of dijet mass plot
0003 #               with corrected and uncorrected jets
0004 #               for various jet algorithms
0005 # Author: Kalanand Mishra
0006 # Date:  22 - November - 2009 
0007 
0008 import FWCore.ParameterSet.Config as cms
0009 
0010 process = cms.Process("Ana")
0011 process.load("FWCore.MessageService.MessageLogger_cfi")
0012 #############   Set the number of events #############
0013 process.maxEvents = cms.untracked.PSet(
0014     input = cms.untracked.int32(1000)
0015 )
0016 #############   Define the source file ###############
0017 process.source = cms.Source("PoolSource",
0018     fileNames = cms.untracked.vstring(
0019     '/store/mc/Summer09/QCDFlat_Pt15to3000/GEN-SIM-RECO/MC_31X_V9_7TeV-v1/0000/FABD2A94-C0D3-DE11-B6FD-00237DA13C2E.root')
0020 )
0021 process.source.inputCommands = cms.untracked.vstring("keep *","drop *_MEtoEDMConverter_*_*")
0022 
0023 #############   Include the jet corrections ##########
0024 process.load("JetMETCorrections.Configuration.L2L3Corrections_Summer09_7TeV_ReReco332_cff")
0025 # set the record's IOV. Must be defined once. Choose ANY correction service. #
0026 process.prefer("L2L3JetCorrectorAK5Calo")
0027 
0028 #############   User analyzer (calo jets) ##
0029 process.DijetMassCaloJets = cms.EDAnalyzer("DijetMassCaloJets",
0030     # Uncorrected CaloJets
0031     AKJets           = cms.string('ak5CaloJets'),
0032     ICJets           = cms.string('iterativeCone5CaloJets'),
0033     SCJets           = cms.string('sisCone5CaloJets'),
0034     KTJets           = cms.string('kt4CaloJets'),
0035     # Corrected CaloJets                                          
0036     AKCorrectedJets  = cms.string('L2L3CorJetAK5Calo'), 
0037     ICCorrectedJets  = cms.string('L2L3CorJetIC5Calo'), 
0038     SCCorrectedJets  = cms.string('L2L3CorJetSC5Calo'), 
0039     KTCorrectedJets  = cms.string('L2L3CorJetKT4Calo'),
0040     # Name of the output ROOT file containing the histograms 
0041     HistoFileName = cms.untracked.string('DijetMassCaloJets.root')
0042 )
0043 
0044 #############   User analyzer (PF jets) ##
0045 process.DijetMassPFJets = cms.EDAnalyzer("DijetMassPFJets",
0046     # Uncorrected PFJets
0047     AKJets           = cms.string('ak5PFJets'),
0048     ICJets           = cms.string('iterativeCone5PFJets'),
0049     SCJets           = cms.string('sisCone5PFJets'),
0050     KTJets           = cms.string('kt4PFJets'),
0051     # Corrected PFJets                                          
0052     AKCorrectedJets  = cms.string('L2L3CorJetAK5PF'), 
0053     ICCorrectedJets  = cms.string('L2L3CorJetIC5PF'), 
0054     SCCorrectedJets  = cms.string('L2L3CorJetSC5PF'), 
0055     KTCorrectedJets  = cms.string('L2L3CorJetKT4PF'),
0056     # Name of the output ROOT file containing the histograms 
0057     HistoFileName = cms.untracked.string('DijetMassPFJets.root')
0058 )
0059 
0060 
0061 #############   User analyzer (gen jets) ##
0062 # ak5GenJets are NOT there: First load the needed modules
0063 process.load("RecoJets.Configuration.GenJetParticles_cff")
0064 process.load("RecoJets.JetProducers.ak5GenJets_cfi")
0065 process.DijetMassGenJets = cms.EDAnalyzer("DijetMassGenJets",
0066     # Uncorrected GenJets
0067     AKJets           = cms.string('ak5GenJets'),
0068     ICJets           = cms.string('iterativeCone5GenJets'),
0069     SCJets           = cms.string('sisCone5GenJets'),
0070     KTJets           = cms.string('kt4GenJets'),
0071     # Corrected GenJets  == Uncorrected GenJets                                        
0072     AKCorrectedJets  = cms.string('ak5GenJets'), 
0073     ICCorrectedJets  = cms.string('iterativeCone5GenJets'), 
0074     SCCorrectedJets  = cms.string('sisCone5GenJets'), 
0075     KTCorrectedJets  = cms.string('kt4GenJets'),
0076     # Name of the output ROOT file containing the histograms 
0077     HistoFileName = cms.untracked.string('DijetMassGenJets.root')
0078 )
0079 
0080 
0081 #############   Path       ###########################
0082 process.p = cms.Path(process.L2L3CorJetAK5Calo * process.L2L3CorJetIC5Calo * process.L2L3CorJetSC5Calo * process.L2L3CorJetKT4Calo * process.DijetMassCaloJets)
0083 process.p2 = cms.Path(process.L2L3CorJetAK5PF * process.L2L3CorJetIC5PF * process.L2L3CorJetSC5PF * process.L2L3CorJetKT4PF * process.DijetMassPFJets)
0084 process.p3 = cms.Path(process.genParticlesForJets * process.ak5GenJets * process.DijetMassGenJets)
0085 
0086 #############   Format MessageLogger #################
0087 process.MessageLogger.cerr.FwkReport.reportEvery = 10
0088