Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-11-26 02:34:12

0001 ##-----
0002 # Set job-specific inputs based on shell
0003 # the following enviromental variables
0004 # export COSMIC_MODE=True # or False
0005 # export JOB_NAME="CosmicStream112220"
0006 # export NEVENTS=1000
0007 # export ALL_HISTS=True
0008 # export TRIGGER_SET=HLT
0009 # export READ_LIST_FROM_FILE=False # or True
0010 # export INPUTFILES='' # root file(s)
0011 # export INPUTFILES_LIST='inputfile_BeamHaloExpress_120015.txt'
0012 # for details see https://twiki.cern.ch/twiki/bin/view/CMS/JetMETDQMPromptAnalysis
0013 ##-----
0014 import os
0015 import FWCore.ParameterSet.Config as cms
0016 #
0017 # --- [cosmic sequence (default=True)?]
0018 iscosmics = (os.environ.get('COSMIC_MODE','False'))
0019 print('iscosmics (default=True) = '+str(iscosmics))
0020 #
0021 # --- [name of job & output file (default=test)?]
0022 jobname = (os.environ.get('JOB_NAME','test'))
0023 print('jobname (default=test) = '+str(jobname))
0024 #
0025 # --- [number of events (default=1000)]
0026 nevents = int(os.environ.get('NEVENTS','1000'))
0027 print('nevents (default=1000)  = '+str(nevents))
0028 #
0029 # --- [turn on all histograms (default=True)?]
0030 allhist = (os.environ.get('ALL_HISTS','True'))
0031 print('allhist (default=True) = '+str(allhist))
0032 #
0033 #--- [read list of input files from a text file? or not (default=False)]
0034 read_from_file = (os.environ.get('READ_LIST_FROM_FILE','True'))
0035 print('read list of input files from a text file (default=False) = '+str(read_from_file))
0036 #
0037 #--- [trigger set (default=HLT)]
0038 trigger_set = (os.environ.get('TRIGGER_SET','HLT'))
0039 print('trigger set name (default=HLT) = '+str(trigger_set))
0040 #
0041 #--- [define list of input files]
0042 inputfiles = []
0043 if read_from_file=="True":
0044   #--- [name of the text file (default=inputfile_list_default.txt)]
0045   filename = (os.environ.get('INPUTFILES_LIST','inputfile_list_default.txt'))
0046   file=open(filename)
0047   print(file.read())
0048   f = open(filename)
0049   try:
0050     for line in f:
0051         inputfiles.append(line)
0052   finally:
0053     f.close()
0054 else:
0055   inputfiles = os.environ.get('INPUTFILES',
0056   '/store/hidata/HIRun2010/HIAllPhysics/RECO/PromptReco-v2/000/150/314/D28FA8E9-4EEA-DF11-9B2A-000423D987E0.root').split(",")
0057   #'/store/data/Commissioning09/MinimumBias/RECO/v4/000/102/347/F85D1BC6-A06A-DE11-BDF8-0019B9F581C9.root').split(",")
0058   #'/store/data/CRAFT09/Calo/RECO/v1/000/112/220/F0B768A4-5E93-DE11-B222-000423D94524.root').split(",")
0059 
0060 print('List of input files')
0061 print(inputfiles)
0062 #-----
0063 
0064 #
0065 #-----
0066 process = cms.Process("test")
0067 process.load("CondCore.DBCommon.CondDBSetup_cfi")
0068 
0069 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0070 
0071 #
0072 # DQM
0073 #
0074 process.load("DQMServices.Core.DQM_cfg")
0075 
0076 process.load("DQMServices.Components.MEtoEDMConverter_cfi")
0077 
0078 ###for HI running
0079 #process.load('Configuration.StandardSequences.MagneticField_38T_cff')
0080 process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
0081 process.load('Configuration.StandardSequences.Reconstruction_cff')
0082 process.load('RecoJets.JetAssociationProducers.trackExtrapolator_cfi')
0083 process.load('RecoJets.JetAssociationProducers.ak7JTA_cff') 
0084 process.load('RecoJets.JetAssociationProducers.ak5JTA_cff')
0085 process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
0086 process.load('HLTrigger/HLTfilters/hltLevel1GTSeed_cfi')
0087 process.load('Configuration/StandardSequences/FrontierConditions_GlobalTag_cff')
0088 
0089 process.GlobalTag.globaltag ='GR_P_V27A::All'
0090 
0091 # the task - JetMET objects
0092 if iscosmics =="True":
0093   process.load("DQMOffline.JetMET.jetMETDQMOfflineSourceCosmic_cff")
0094 else:
0095   process.load("DQMOffline.JetMET.jetMETDQMOfflineSourceHI_cff")
0096 
0097 process.jetMETAnalyzer.OutputMEsInRootFile = cms.bool(True)
0098 process.jetMETAnalyzer.OutputFileName = cms.string("jetMETMonitoring_%s.root" % jobname)
0099 process.jetMETAnalyzer.TriggerResultsLabel = cms.InputTag("TriggerResults","",trigger_set)
0100 process.jetMETAnalyzer.processname = cms.string(trigger_set)
0101 
0102 # check # of bins
0103 process.load("DQMServices.Components.DQMStoreStats_cfi")
0104 
0105 # for igprof
0106 #process.IgProfService = cms.Service("IgProfService",
0107 #  reportFirstEvent            = cms.untracked.int32(0),
0108 #  reportEventInterval         = cms.untracked.int32(25),
0109 #  reportToFileAtPostEvent     = cms.untracked.string("| gzip -c > igdqm.%I.gz")
0110 #)
0111 
0112 #
0113 # /Wmunu/Summer09-MC_31X_V3-v1/GEN-SIM-RECO
0114 process.source = cms.Source("PoolSource",
0115     fileNames = cms.untracked.vstring(*inputfiles))
0116 
0117 #
0118 process.source.inputCommands = cms.untracked.vstring('keep *', 'drop *_MEtoEDMConverter_*_*')
0119 
0120 #
0121 process.maxEvents = cms.untracked.PSet(
0122     input = cms.untracked.int32( nevents )
0123 )
0124 process.Timing = cms.Service("Timing")
0125 
0126 ## # Comment this out or reconfigure to see error messages 
0127 process.MessageLogger = cms.Service("MessageLogger",
0128     debugModules = cms.untracked.vstring('jetMETAnalyzer'),
0129     cout = cms.untracked.PSet(
0130         default = cms.untracked.PSet(
0131             limit = cms.untracked.int32(0)
0132         ),
0133         jetMETAnalyzer = cms.untracked.PSet(
0134             limit = cms.untracked.int32(1)
0135         ),
0136         noLineBreaks = cms.untracked.bool(True),
0137         DEBUG = cms.untracked.PSet(
0138             limit = cms.untracked.int32(0)
0139         ),
0140         threshold = cms.untracked.string('DEBUG')
0141     ),
0142     categories = cms.untracked.vstring('jetMETAnalyzer'),
0143     destinations = cms.untracked.vstring('cout')
0144 )
0145 
0146 
0147 process.options = cms.untracked.PSet(
0148     wantSummary = cms.untracked.bool(True)
0149 )
0150 
0151 process.FEVT = cms.OutputModule("PoolOutputModule",
0152     outputCommands = cms.untracked.vstring('keep *_MEtoEDMConverter_*_*'),
0153     #outputCommands = cms.untracked.vstring('keep *'),
0154     fileName = cms.untracked.string("reco_DQM_%s.root" % jobname)
0155 )
0156 
0157 process.options = cms.untracked.PSet(
0158     wantSummary = cms.untracked.bool(True) ## default is false
0159 
0160 )
0161 
0162 if iscosmics=="True":
0163   process.p = cms.Path(process.jetMETDQMOfflineSourceCosmic
0164                      * process.dqmStoreStats
0165                      * process.MEtoEDMConverter
0166                      )
0167 else:
0168   process.p = cms.Path(process.jetMETDQMOfflineSource
0169                      * process.dqmStoreStats
0170                      * process.MEtoEDMConverter
0171                      )
0172 
0173 process.outpath = cms.EndPath(process.FEVT)
0174 process.DQM.collectorHost = ''
0175 
0176