Skims/python/reco_skim_cfg_mod.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process('RERECO')
0004
0005 # this is to avoid the postpathendrun probem with same process name (only with http reader)
0006 process.options = cms.untracked.PSet(
0007     IgnoreCompletely = cms.untracked.vstring('Configuration')
0008 #    TryToContinue = cms.untracked.vstring('Configuration')
0009 )
0010
0011
0012 # for ispy
0013 process.add_(
0014     cms.Service("ISpyService",
0015     outputFileName = cms.untracked.string('Ispy.ig'),
0016     outputMaxEvents = cms.untracked.int32 (1000),
0017     online = cms.untracked.bool(True),
0018     debug = cms.untracked.bool(True)
0019     )
0020 )
0021
0022
0023 # import of standard configurations
0024 process.load('Configuration/StandardSequences/Services_cff')
0025 process.load('FWCore/MessageService/MessageLogger_cfi')
0026 process.load('Configuration/StandardSequences/GeometryIdeal_cff')
0027 process.load('Configuration/StandardSequences/MagneticField_AutoFromDBCurrent_cff')
0028 process.load('Configuration/StandardSequences/RawToDigi_Data_cff')
0029 process.load('Configuration/StandardSequences/Reconstruction_cff')
0030 process.load('DQMOffline/Configuration/DQMOffline_cff')
0031 process.load('Configuration/StandardSequences/EndOfProcess_cff')
0032 process.load('Configuration/StandardSequences/FrontierConditions_GlobalTag_cff')
0033 process.load('Configuration/EventContent/EventContent_cff')
0034 process.load('ISpy/Analyzers/ISpy_Producer_cff')
0035
0036
0037
0038 ######### FILTERING Section #############################
0039
0040 # this is for filtering on HLT path
0041 process.hltHighLevel = cms.EDFilter("HLTHighLevel",
0042      TriggerResultsTag = cms.InputTag("TriggerResults","","HLT"),
0043 #     HLTPaths = cms.vstring('HLT_Activity_L1A'),           # provide list of HLT paths (or patterns) you want
0044      HLTPaths = cms.vstring('HLT_MinBiasBSC'),           # provide list of HLT paths (or patterns) you want
0045      eventSetupPathsKey = cms.string(''), # not empty => use read paths from AlCaRecoTriggerBitsRcd via this key
0046      andOr = cms.bool(True),             # how to deal with multiple triggers: True (OR) accept if ANY is true, False (AND) accept if ALL are true
0047      throw = cms.bool(True),    # throw exception on unknown path names
0048      saveTags = cms.bool(False)
0049  )
0050
0051 # this is for filtering based on reco variables
0052 process.skimming = cms.EDFilter("BeamSplash",
0053     energycuttot = cms.untracked.double(1000.0),
0054     energycutecal = cms.untracked.double(700.0),
0055     energycuthcal = cms.untracked.double(700.0),
0056     ebrechitcollection =   cms.InputTag("ecalRecHit","EcalRecHitsEB"),
0057     eerechitcollection =   cms.InputTag("ecalRecHit","EcalRecHitsEE"),
0058     hbherechitcollection =   cms.InputTag("hbhereco"),
0059     applyfilter = cms.untracked.bool(False)
0060 )
0061
0062 # this is for filtering on trigger type
0063
0064 process.load("HLTrigger.special.HLTTriggerTypeFilter_cfi")
0065 # 0=random, 1=physics, 2=calibration, 3=technical
0066 process.hltTriggerTypeFilter.SelectedTriggerType = 1
0067
0068 # this is for filtering on L1 technical trigger bit
0069 process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
0070 process.load('HLTrigger/HLTfilters/hltLevel1GTSeed_cfi')
0071 process.hltLevel1GTSeed.L1TechTriggerSeeding = cms.bool(True)
0072 process.hltLevel1GTSeed.L1SeedsLogicalExpression = cms.string('32 OR 33 OR 40 OR 41')
0073
0074 #this is for filtering/tagging PhysDecl bit
0075 process.physdecl = cms.EDFilter("PhysDecl",
0076      applyfilter = cms.untracked.bool(False),
0077      debugOn = cms.untracked.bool(True)
0078     )
0079
0080
0081 process.configurationMetadata = cms.untracked.PSet(
0082     version = cms.untracked.string('$Revision: 1.4 $'),
0083     annotation = cms.untracked.string('promptReco nevts:1'),
0084     name = cms.untracked.string('PyReleaseValidation')
0085 )
0086 process.maxEvents = cms.untracked.PSet(
0087     input = cms.untracked.int32(NUMEVENTS)
0088 )
0089 process.options = cms.untracked.PSet(
0090     Rethrow = cms.untracked.vstring('ProductNotFound'),
0091     wantSummary = cms.untracked.bool(True)
0092 )
0093
0094 process.source = cms.Source("EventStreamHttpReader",
0095
0096 # streaming##################################################
0097 # in p5
0098 #   sourceURL = cms.string('http://srv-c2d05-14.cms:22100/urn:xdaq-application:lid=30'),
0099 #   consumerName = cms.untracked.string('DQM Source'),
0100
0101 # tunnel to proxy
0102 # THIS SHOULD BE THE CORRECT FOR OFFLINE ACCESSING THE REVERSE PROXY
0103 #   sourceURL = cms.string('http://cmsdaq0.cern.ch/event-server/urn:xdaq-application:lid=30'),
0104
0105 # special tunnel configuration, need to setup an external tunnel
0106 #   sourceURL = cms.string('http://localhost:22100/urn:xdaq-application:lid=30'),
0107    sourceURL = SOURCE,
0108    consumerName = cms.untracked.string('Event Display'),
0109
0110 # direct storage manager
0111 #   sourceURL = cms.string('http://localhost:22100/urn:xdaq-application:service=storagemanager'),
0112 #   consumerName = cms.untracked.string('Event Display'),
0113
0114 # playback###################################################
0115 # in pt5
0116 #    sourceURL = cms.string('http://srv-c2d05-05:50082/urn:xdaq-application:lid=29'),
0117
0118 # tunnel
0119 #   sourceURL = cms.string('http://localhost:50082/urn:xdaq-application:lid=29'),
0120 #   consumerName = cms.untracked.string('Playback Source'),
0121 #################################################################
0122
0123    consumerPriority = cms.untracked.string('normal'),
0124    max_event_size = cms.int32(7000000),
0125    SelectHLTOutput = SELECTHLT,
0126 #   SelectHLTOutput = cms.untracked.string('hltOutputDQM'),
0127 #   SelectHLTOutput = cms.untracked.string('hltOutputExpress'),
0128    max_queue_depth = cms.int32(5),
0129    maxEventRequestRate = cms.untracked.double(2.0),
0130    SelectEvents = cms.untracked.PSet(
0131 #       SelectEvents = cms.vstring('*DQM')
0132        SelectEvents = cms.vstring('*')
0133 #       SelectEvents = cms.vstring('PhysicsPath')
0134    ),
0135    headerRetryInterval = cms.untracked.int32(3)
0136 )
0137
0138
0139
0140 #process.source = cms.Source("PoolSource",
0141 #    debugVerbosity = cms.untracked.uint32(0),
0142 #    debugFlag = cms.untracked.bool(False),
0143 #    fileNames = cms.untracked.vstring(
0144 ##'/store/data/Commissioning08/BeamHalo/RECO/StuffAlmostToP5_v1/000/061/642/10A0FE34-A67D-DD11-AD05-000423D94E1C.root'
0145 ##
0146 ##'/store/express/CRAFT09/ExpressMuon/FEVT/v1/000/110/835/FED0EFCD-AB87-DE11-9B72-000423D99658.root'
0147 ##'/store/express/CRAFT09/ExpressMuon/FEVT/v1/000/110/835/FC629BD2-CF87-DE11-9077-001D09F25438.root',
0148 ##'/store/express/CRAFT09/ExpressMuon/FEVT/v1/000/110/835/FC38EE75-BD87-DE11-822A-001D09F253C0.root',
0149 ##'/store/express/CRAFT09/ExpressMuon/FEVT/v1/000/110/835/FC1CB101-A487-DE11-9F10-000423D99660.root'
0150 #))
0151
0152
0153 process.FEVT = cms.OutputModule("PoolOutputModule",
0154     maxSize = cms.untracked.int32(1000),
0155     fileName = cms.untracked.string('EVDISPSM_SUFFIX.root'),
0156     outputCommands = cms.untracked.vstring('keep *','drop *_MEtoEDMConverter_*_*'),
0157     dataset = cms.untracked.PSet(
0158               dataTier = cms.untracked.string('RAW-RECO'),
0159 #             filterName = cms.untracked.string(''))
0160               filterName = cms.untracked.string('EVDISP')),
0161                                 SelectEvents = cms.untracked.PSet(
0162         SelectEvents = cms.vstring('fullpath')
0163     )
0164 )
0165
0166 # Other statements
0167 #process.GlobalTag.connect = 'sqlite_file:/afs/cern.ch/user/m/malgeri/public/gtfirstcoll.db'
0168 process.GlobalTag.globaltag = 'GR09_P_V6::All'
0169
0170 process.fifthCkfTrajectoryFilter.filterPset.minimumNumberOfHits = 2
0171 process.fifthCkfTrajectoryFilter.filterPset.maxLostHits = 4
0172 process.fifthCkfTrajectoryFilter.filterPset.maxConsecLostHits = 2
0173 process.fifthCkfInOutTrajectoryFilter.filterPset.minimumNumberOfHits = 2
0174 process.fifthCkfInOutTrajectoryFilter.filterPset.maxLostHits = 4
0175 process.fifthCkfInOutTrajectoryFilter.filterPset.maxConsecLostHits = 2
0176 process.fifthCkfTrajectoryBuilder.minNrOfHitsForRebuild = 2
0177 process.fifthRKTrajectorySmoother.minHits = 2
0178 process.fifthRKTrajectoryFitter.minHits = 2
0179 process.fifthFittingSmootherWithOutlierRejection.MinNumberOfHits = 2
0180 process.tobtecStepLoose.minNumberLayers = 2
0181 process.tobtecStepLoose.maxNumberLostLayers = 2
0182 process.tobtecStepLoose.dz_par1 = cms.vdouble(10.5, 4.0)
0183 process.tobtecStepLoose.dz_par2 = cms.vdouble(10.5, 4.0)
0184 process.tobtecStepLoose.d0_par1 = cms.vdouble(10.5, 4.0)
0185 process.tobtecStepLoose.d0_par2 = cms.vdouble(10.5, 4.0)
0186 process.tobtecStepLoose.chi2n_par = cms.double(100.0)
0187 process.fifthSeeds.RegionFactoryPSet.RegionPSet.originHalfLength = 100
0188 process.fifthSeeds.RegionFactoryPSet.RegionPSet.originRadius     = 10
0189 process.Chi2MeasurementEstimator.MaxChi2 = 100
0190
0191
0192 # to filter on MinBias...
0193 #process.fullpath = cms.Path(process.hltTriggerTypeFilter+process.hltHighLevel+process.RawToDigi+process.reconstruction)
0194
0195 # to filter on trigger type only
0196 #process.fullpath = cms.Path(process.hltTriggerTypeFilter+process.hltHighLevel+process.RawToDigi+process.reconstruction)
0197
0198 #process.fullpath = cms.Path(process.RawToDigi+process.reconstruction+process.skimming+process.iSpy_sequence)
0199 #process.fullpath = cms.Path(process.hltTriggerTypeFilter+process.RawToDigi+process.reconstruction+process.skimming+process.iSpy_sequence)
0200 # added physdecl in tagging mode to catch physdeclared bit in log files
0201 # process.fullpath = cms.Path(process.hltTriggerTypeFilter+process.RawToDigi+process.physdecl+process.reconstruction+process.skimming+process.iSpy_sequence)
0202 process.fullpath = cms.Path(process.RawToDigi+process.physdecl+process.reconstruction+process.skimming+process.iSpy_sequence)
0203
0204 process.out_step = cms.EndPath(process.FEVT)
0205
0206 # Schedule definition
0207
0208 process.schedule = cms.Schedule(process.fullpath,process.out_step)
0209
0210
0211 #process.e = cms.EndPath(process.out)
0212