JetAnalyzers/test/GenJetAna.py

0001 # PYTHON configuration file.
0002 # Description:  Example of analysis pythia produced partons & jets
0003 # Author: R. Harris
0004 # Date:  28 - October - 2008
0005 import FWCore.ParameterSet.Config as cms
0006
0007 process = cms.Process("Ana")
0008 process.load("FWCore.MessageService.MessageLogger_cfi")
0009 #############   Set the number of events #############
0010 process.maxEvents = cms.untracked.PSet(
0011     input = cms.untracked.int32(-1)
0012 )
0013 #############   Define the source file ###############
0014 #process.load("RecoJets.JetAnalyzers.QCD_GenJets_cfi")
0015 process.source = cms.Source("PoolSource",
0016 #    noEventSort = cms.untracked.bool(True),
0017     duplicateCheckMode = cms.untracked.string('noDuplicateCheck'),
0018     fileNames = cms.untracked.vstring(
0019 'file:/uscms_data/d2/rharris/CMSSW_3_3_2/src/PYTHIA6_QCDpt_0_15_900GeV_cff_py_GEN.root',
0020 'file:/uscms_data/d2/rharris/CMSSW_3_3_2/src/PYTHIA6_QCDpt_15_20_900GeV_cff_py_GEN.root',
0021 'file:/uscms_data/d2/rharris/CMSSW_3_3_2/src/PYTHIA6_QCDpt_20_30_900GeV_cff_py_GEN.root',
0022 'file:/uscms_data/d2/rharris/CMSSW_3_3_2/src/PYTHIA6_QCDpt_30_50_900GeV_cff_py_GEN.root',
0023 'file:/uscms_data/d2/rharris/CMSSW_3_3_2/src/PYTHIA6_QCDpt_50_80_900GeV_cff_py_GEN.root',
0024 'file:/uscms_data/d2/rharris/CMSSW_3_3_2/src/PYTHIA6_QCDpt_80_120_900GeV_cff_py_GEN.root',
0025 'file:/uscms_data/d2/rharris/CMSSW_3_3_2/src/PYTHIA6_QCDpt_120_inf_900GeV_cff_py_GEN.root')
0026 )
0027 #############   User analyzer  #####
0028 # Jet Algorithm: which jets to use for making histograms and ana root tree.
0029 # HistoFleName:  name of file which contains histograms and ana root tree.
0030 # debug:         set true for printout each event
0031 # NJets:         Currently only the value 2 is supported for dijet ana.
0032 # eventsGen:     Number of events generated, or will be generated, per pt_hat bin.
0033 # anaLevel:      Analysis level string which steers analysis choices
0034 #               "PtHatOnly": only get PtHat and make PtHat histos. FAST.
0035 #               "Jets":  do analysis of jets, but not partons.
0036 #                        Input file only needs to have GenJets.
0037 #               "all":   do analysis of everything and make histos and root tree
0038 #                        Input file needs to have GenParticles as well as GenJets
0039 #               "generating": prior modes all assumd we were running on file.
0040 #                             This mode assumes we are generating pythia here.
0041 #                            Analysis of everything, make histos and root tree
0042 # xsecGen:       Generator cross section.
0043 #                   Just one value if analysis level is "generating"
0044 #                   Multiple values if analysis level is "all"
0045 # ptHatEdges:    Edges of ptHatBins that correspond to xsecGen values above
0046 #
0047 process.gen = cms.EDAnalyzer("GenJetAnaPythia",
0048     JetAlgorithm    = cms.string('ak7GenJets'),
0049     HistoFileName   = cms.string('QCD900GeVGenJets700Kevts.root'),
0050     debug           = cms.bool(False),
0051     NJets           = cms.int32(2),
0052     eventsGen       = cms.int32(100000),
0053     anaLevel        = cms.string('all'),
0054     xsecGen         = cms.vdouble(3.439E+10, 1.811E+07, 5.469E+06, 7.391E+05, 4.321E+04, 2.200E+03, 1.080E+02),
0055     ptHatEdges      = cms.vdouble(0, 15, 20, 30, 50, 80, 120, 9999)
0056 )
0057 #############   Path       ###########################
0058 process.p = cms.Path(process.gen)
0059 #############   Format MessageLogger #################
0060 process.MessageLogger.cerr.FwkReport.reportEvery = 1000
0061