CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
reco_skim_cfg_mod.py
Go to the documentation of this file.
2 
3 process = cms.Process('RERECO')
4 
5 # this is to avoid the postpathendrun probem with same process name (only with http reader)
6 process.options = cms.untracked.PSet(
7  IgnoreCompletely = cms.untracked.vstring('Configuration')
8 # SkipEvent = cms.untracked.vstring('Configuration')
9 )
10 
11 
12 # for ispy
13 process.add_(
14  cms.Service("ISpyService",
15  outputFileName = cms.untracked.string('Ispy.ig'),
16  outputMaxEvents = cms.untracked.int32 (1000),
17  online = cms.untracked.bool(True),
18  debug = cms.untracked.bool(True)
19  )
20 )
21 
22 
23 # import of standard configurations
24 process.load('Configuration/StandardSequences/Services_cff')
25 process.load('FWCore/MessageService/MessageLogger_cfi')
26 process.load('Configuration/StandardSequences/GeometryIdeal_cff')
27 process.load('Configuration/StandardSequences/MagneticField_AutoFromDBCurrent_cff')
28 process.load('Configuration/StandardSequences/RawToDigi_Data_cff')
29 process.load('Configuration/StandardSequences/Reconstruction_cff')
30 process.load('DQMOffline/Configuration/DQMOffline_cff')
31 process.load('Configuration/StandardSequences/EndOfProcess_cff')
32 process.load('Configuration/StandardSequences/FrontierConditions_GlobalTag_cff')
33 process.load('Configuration/EventContent/EventContent_cff')
34 process.load('ISpy/Analyzers/ISpy_Producer_cff')
35 
36 
37 
38 ######### FILTERING Section #############################
39 
40 # this is for filtering on HLT path
41 process.hltHighLevel = cms.EDFilter("HLTHighLevel",
42  TriggerResultsTag = cms.InputTag("TriggerResults","","HLT"),
43 # HLTPaths = cms.vstring('HLT_Activity_L1A'), # provide list of HLT paths (or patterns) you want
44  HLTPaths = cms.vstring('HLT_MinBiasBSC'), # provide list of HLT paths (or patterns) you want
45  eventSetupPathsKey = cms.string(''), # not empty => use read paths from AlCaRecoTriggerBitsRcd via this key
46  andOr = cms.bool(True), # how to deal with multiple triggers: True (OR) accept if ANY is true, False (AND) accept if ALL are true
47  throw = cms.bool(True), # throw exception on unknown path names
48  saveTags = cms.bool(False)
49  )
50 
51 # this is for filtering based on reco variables
52 process.skimming = cms.EDFilter("BeamSplash",
53  energycuttot = cms.untracked.double(1000.0),
54  energycutecal = cms.untracked.double(700.0),
55  energycuthcal = cms.untracked.double(700.0),
56  ebrechitcollection = cms.InputTag("ecalRecHit","EcalRecHitsEB"),
57  eerechitcollection = cms.InputTag("ecalRecHit","EcalRecHitsEE"),
58  hbherechitcollection = cms.InputTag("hbhereco"),
59  applyfilter = cms.untracked.bool(False)
60 )
61 
62 # this is for filtering on trigger type
63 
64 process.load("HLTrigger.special.HLTTriggerTypeFilter_cfi")
65 # 0=random, 1=physics, 2=calibration, 3=technical
66 process.hltTriggerTypeFilter.SelectedTriggerType = 1
67 
68 # this is for filtering on L1 technical trigger bit
69 process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
70 process.load('HLTrigger/HLTfilters/hltLevel1GTSeed_cfi')
71 process.hltLevel1GTSeed.L1TechTriggerSeeding = cms.bool(True)
72 process.hltLevel1GTSeed.L1SeedsLogicalExpression = cms.string('32 OR 33 OR 40 OR 41')
73 
74 #this is for filtering/tagging PhysDecl bit
75 process.physdecl = cms.EDFilter("PhysDecl",
76  applyfilter = cms.untracked.bool(False),
77  debugOn = cms.untracked.bool(True)
78  )
79 
80 
81 process.configurationMetadata = cms.untracked.PSet(
82  version = cms.untracked.string('$Revision: 1.5 $'),
83  annotation = cms.untracked.string('promptReco nevts:1'),
84  name = cms.untracked.string('PyReleaseValidation')
85 )
86 process.maxEvents = cms.untracked.PSet(
87  input = cms.untracked.int32(NUMEVENTS)
88 )
89 process.options = cms.untracked.PSet(
90  Rethrow = cms.untracked.vstring('ProductNotFound'),
91  wantSummary = cms.untracked.bool(True)
92 )
93 
94 process.source = cms.Source("EventStreamHttpReader",
95 
96 # streaming##################################################
97 # in p5
98 # sourceURL = cms.string('http://srv-c2d05-14.cms:22100/urn:xdaq-application:lid=30'),
99 # consumerName = cms.untracked.string('DQM Source'),
100 
101 # tunnel to proxy
102 # THIS SHOULD BE THE CORRECT FOR OFFLINE ACCESSING THE REVERSE PROXY
103 # sourceURL = cms.string('http://cmsdaq0.cern.ch/event-server/urn:xdaq-application:lid=30'),
104 
105 # special tunnel configuration, need to setup an external tunnel
106 # sourceURL = cms.string('http://localhost:22100/urn:xdaq-application:lid=30'),
107  sourceURL = SOURCE,
108  consumerName = cms.untracked.string('Event Display'),
109 
110 # direct storage manager
111 # sourceURL = cms.string('http://localhost:22100/urn:xdaq-application:service=storagemanager'),
112 # consumerName = cms.untracked.string('Event Display'),
113 
114 # playback###################################################
115 # in pt5
116 # sourceURL = cms.string('http://srv-c2d05-05:50082/urn:xdaq-application:lid=29'),
117 
118 # tunnel
119 # sourceURL = cms.string('http://localhost:50082/urn:xdaq-application:lid=29'),
120 # consumerName = cms.untracked.string('Playback Source'),
121 #################################################################
122 
123  consumerPriority = cms.untracked.string('normal'),
124  max_event_size = cms.int32(7000000),
125  SelectHLTOutput = SELECTHLT,
126 # SelectHLTOutput = cms.untracked.string('hltOutputDQM'),
127 # SelectHLTOutput = cms.untracked.string('hltOutputExpress'),
128  max_queue_depth = cms.int32(5),
129  maxEventRequestRate = cms.untracked.double(2.0),
130  SelectEvents = cms.untracked.PSet(
131 # SelectEvents = cms.vstring('*DQM')
132  SelectEvents = cms.vstring('*')
133 # SelectEvents = cms.vstring('PhysicsPath')
134  ),
135  headerRetryInterval = cms.untracked.int32(3)
136 )
137 
138 
139 
140 #process.source = cms.Source("PoolSource",
141 # debugVerbosity = cms.untracked.uint32(0),
142 # debugFlag = cms.untracked.bool(False),
143 # fileNames = cms.untracked.vstring(
144 ##'/store/data/Commissioning08/BeamHalo/RECO/StuffAlmostToP5_v1/000/061/642/10A0FE34-A67D-DD11-AD05-000423D94E1C.root'
145 ##
146 ##'/store/express/CRAFT09/ExpressMuon/FEVT/v1/000/110/835/FED0EFCD-AB87-DE11-9B72-000423D99658.root'
147 ##'/store/express/CRAFT09/ExpressMuon/FEVT/v1/000/110/835/FC629BD2-CF87-DE11-9077-001D09F25438.root',
148 ##'/store/express/CRAFT09/ExpressMuon/FEVT/v1/000/110/835/FC38EE75-BD87-DE11-822A-001D09F253C0.root',
149 ##'/store/express/CRAFT09/ExpressMuon/FEVT/v1/000/110/835/FC1CB101-A487-DE11-9F10-000423D99660.root'
150 #))
151 
152 
153 process.FEVT = cms.OutputModule("PoolOutputModule",
154  maxSize = cms.untracked.int32(1000),
155  fileName = cms.untracked.string('EVDISPSM_SUFFIX.root'),
156  outputCommands = cms.untracked.vstring('keep *','drop *_MEtoEDMConverter_*_*'),
157  dataset = cms.untracked.PSet(
158  dataTier = cms.untracked.string('RAW-RECO'),
159 # filterName = cms.untracked.string(''))
160  filterName = cms.untracked.string('EVDISP')),
161  SelectEvents = cms.untracked.PSet(
162  SelectEvents = cms.vstring('fullpath')
163  )
164 )
165 
166 # Other statements
167 #process.GlobalTag.connect = 'sqlite_file:/afs/cern.ch/user/m/malgeri/public/gtfirstcoll.db'
168 process.GlobalTag.globaltag = 'GR09_P_V6::All'
169 
170 process.fifthCkfTrajectoryFilter.filterPset.minimumNumberOfHits = 2
171 process.fifthCkfTrajectoryFilter.filterPset.maxLostHits = 4
172 process.fifthCkfTrajectoryFilter.filterPset.maxConsecLostHits = 2
173 process.fifthCkfInOutTrajectoryFilter.filterPset.minimumNumberOfHits = 2
174 process.fifthCkfInOutTrajectoryFilter.filterPset.maxLostHits = 4
175 process.fifthCkfInOutTrajectoryFilter.filterPset.maxConsecLostHits = 2
176 process.fifthCkfTrajectoryBuilder.minNrOfHitsForRebuild = 2
177 process.fifthRKTrajectorySmoother.minHits = 2
178 process.fifthRKTrajectoryFitter.minHits = 2
179 process.fifthFittingSmootherWithOutlierRejection.MinNumberOfHits = 2
180 process.tobtecStepLoose.minNumberLayers = 2
181 process.tobtecStepLoose.maxNumberLostLayers = 2
182 process.tobtecStepLoose.dz_par1 = cms.vdouble(10.5, 4.0)
183 process.tobtecStepLoose.dz_par2 = cms.vdouble(10.5, 4.0)
184 process.tobtecStepLoose.d0_par1 = cms.vdouble(10.5, 4.0)
185 process.tobtecStepLoose.d0_par2 = cms.vdouble(10.5, 4.0)
186 process.tobtecStepLoose.chi2n_par = cms.double(100.0)
187 process.fifthSeeds.RegionFactoryPSet.RegionPSet.originHalfLength = 100
188 process.fifthSeeds.RegionFactoryPSet.RegionPSet.originRadius = 10
189 process.Chi2MeasurementEstimator.MaxChi2 = 100
190 
191 
192 # to filter on MinBias...
193 #process.fullpath = cms.Path(process.hltTriggerTypeFilter+process.hltHighLevel+process.RawToDigi+process.reconstruction)
194 
195 # to filter on trigger type only
196 #process.fullpath = cms.Path(process.hltTriggerTypeFilter+process.hltHighLevel+process.RawToDigi+process.reconstruction)
197 
198 #process.fullpath = cms.Path(process.RawToDigi+process.reconstruction+process.skimming+process.iSpy_sequence)
199 #process.fullpath = cms.Path(process.hltTriggerTypeFilter+process.RawToDigi+process.reconstruction+process.skimming+process.iSpy_sequence)
200 # added physdecl in tagging mode to catch physdeclared bit in log files
201 # process.fullpath = cms.Path(process.hltTriggerTypeFilter+process.RawToDigi+process.physdecl+process.reconstruction+process.skimming+process.iSpy_sequence)
202 process.fullpath = cms.Path(process.RawToDigi+process.physdecl+process.reconstruction+process.skimming+process.iSpy_sequence)
203 
204 process.out_step = cms.EndPath(process.FEVT)
205 
206 # Schedule definition
207 
208 process.schedule = cms.Schedule(process.fullpath,process.out_step)
209 
210 
211 #process.e = cms.EndPath(process.out)
212