JetAnalyzers/test/dumpFFTJetGrids.py

0001 #
0002 # This example configuration file dumps FFTJet energy discretization grids
0003 # in the form of root histograms.
0004 #
0005 # I. Volobouev, June 29 2011
0006 #
0007 import FWCore.ParameterSet.Config as cms
0008
0009 process = cms.Process("FFTJetTest")
0010
0011 # Various standard stuff
0012 process.load('Configuration/StandardSequences/GeometryExtended_cff')
0013 process.load('Configuration/StandardSequences/Reconstruction_cff')
0014 process.load('Configuration.StandardSequences.Services_cff')
0015 process.load('Configuration/StandardSequences/FrontierConditions_GlobalTag_cff')
0016 process.GlobalTag.globaltag = 'FT_R_42_V13A::All'
0017
0018 process.load("FWCore.MessageService.MessageLogger_cfi")
0019 process.MessageLogger.cerr.FwkReport.reportEvery = 100
0020
0021 process.TFileService = cms.Service(
0022     "TFileService",
0023     fileName = cms.string("fftjet_grids.root")
0024 )
0025
0026 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )
0027
0028 process.source = cms.Source("PoolSource",
0029     fileNames = cms.untracked.vstring(
0030         '/store/data/Run2011A/MinimumBias/RECO/May10ReReco-v2/0002/FA15E8AA-F782-E011-9D3C-001BFCDBD1BE.root'
0031     )
0032 )
0033
0034 # process.source = cms.Source("PoolSource",
0035 #     fileNames = cms.untracked.vstring(
0036 #         '/store/data/Run2011A/MinimumBias/RECO/May10ReReco-v2/0002/FE8F09A1-8082-E011-990F-001A92971B7E.root',
0037 #         '/store/data/Run2011A/MinimumBias/RECO/May10ReReco-v2/0002/FA15E8AA-F782-E011-9D3C-001BFCDBD1BE.root',
0038 #         '/store/data/Run2011A/MinimumBias/RECO/May10ReReco-v2/0002/F88FDE18-6782-E011-B5D8-003048678ED4.root',
0039 #         '/store/data/Run2011A/MinimumBias/RECO/May10ReReco-v2/0002/F839BD95-F183-E011-BE30-002618FDA262.root',
0040 #         '/store/data/Run2011A/MinimumBias/RECO/May10ReReco-v2/0002/D89DE790-FC82-E011-8CA6-0026189438CB.root',
0041 #         '/store/data/Run2011A/MinimumBias/RECO/May10ReReco-v2/0002/C4BE9561-6782-E011-A791-0030486792B4.root',
0042 #         '/store/data/Run2011A/MinimumBias/RECO/May10ReReco-v2/0002/C0CE5857-8482-E011-9BD9-002354EF3BDE.root',
0043 #         '/store/data/Run2011A/MinimumBias/RECO/May10ReReco-v2/0002/BC5BF6CB-6682-E011-B5E2-00248C0BE01E.root',
0044 #         '/store/data/Run2011A/MinimumBias/RECO/May10ReReco-v2/0002/AA6C29E6-6F82-E011-8F12-001A92971B62.root',
0045 #         '/store/data/Run2011A/MinimumBias/RECO/May10ReReco-v2/0002/A476C008-8A82-E011-A676-002354EF3BDA.root'
0046 #     )
0047 # )
0048
0049 # Configure the trigger
0050 process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
0051 from HLTrigger.HLTfilters.hltLevel1GTSeed_cfi import hltLevel1GTSeed
0052
0053 # L1 technical triggers:
0054 # 0     -- beam crossing
0055 # 36-39 -- beam halo
0056 process.VetoHardInt = hltLevel1GTSeed.clone(
0057     L1TechTriggerSeeding = cms.bool(True),
0058     L1SeedsLogicalExpression = cms.string('0 AND NOT (36 OR 37 OR 38 OR 39)')
0059 )
0060
0061 # Remove the so-called "scraping" events
0062 process.noScraping = cms.EDFilter("FilterOutScraping",
0063     applyfilter = cms.untracked.bool(True),
0064     debugOn = cms.untracked.bool(False), ## Or 'True' to get some per-event info
0065     numtrack = cms.untracked.uint32(10),
0066     thresh = cms.untracked.double(0.25)
0067 )
0068
0069 # min bias filter
0070 process.HLTZeroBias =cms.EDFilter("HLTHighLevel",
0071     TriggerResultsTag = cms.InputTag("TriggerResults","","HLT"),
0072     HLTPaths = cms.vstring('HLT_ZeroBias','HLT_ZeroBias_v1','HLT_ZeroBias_v2','HLT_ZeroBias_v3'),
0073     eventSetupPathsKey = cms.string(''),
0074     andOr = cms.bool(True), #----- True = OR, False = AND between the HLTPaths
0075     throw = cms.bool(False)
0076 )
0077
0078 # Configure the FFTJet pattern recognition module
0079 from RecoJets.FFTJetProducers.fftjetcommon_cfi import *
0080 from RecoJets.FFTJetProducers.fftjetpatrecoproducer_cfi import *
0081
0082 fftjet_patreco_producer.src = cms.InputTag("particleFlow")
0083 fftjet_patreco_producer.jetType = cms.string("PFJet")
0084 fftjet_patreco_producer.storeDiscretizationGrid = cms.bool(True)
0085 fftjet_patreco_producer.makeClusteringTree = cms.bool(False)
0086 fftjet_patreco_producer.GridConfiguration = fftjet_grid_256_128
0087
0088 # Configure the FFTJet pile-up analyzer module to dump the grids
0089 from RecoJets.JetAnalyzers.fftjetpileupanalyzer_cfi import *
0090
0091 fftjet_pileup_analyzer.collectPileup = cms.bool(False)
0092 fftjet_pileup_analyzer.collectSummaries = cms.bool(False)
0093 fftjet_pileup_analyzer.collectGrids = cms.bool(True)
0094 fftjet_pileup_analyzer.collectVertexInfo = cms.bool(True)
0095
0096 process.fftjetpatreco = fftjet_patreco_producer
0097 process.pileupanalyzer = fftjet_pileup_analyzer
0098
0099 process.p = cms.Path(
0100     process.VetoHardInt *
0101     process.HLTZeroBias *
0102     process.noScraping *
0103     process.fftjetpatreco *
0104     process.pileupanalyzer
0105 )