|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
|
import FWCore.ParameterSet.Config as cms
import FWCore.ParameterSet.VarParsing as VarParsing
process = cms.Process("MultiplicityProducerTest")
#prepare options
options = VarParsing.VarParsing("analysis")
options.register ('globalTag',
"DONOTEXIST",
VarParsing.VarParsing.multiplicity.singleton, # singleton or list
VarParsing.VarParsing.varType.string, # string, int, or float
"GlobalTag")
options.register ('fromRAW',
"0",
VarParsing.VarParsing.multiplicity.singleton, # singleton or list
VarParsing.VarParsing.varType.int, # string, int, or float
"=1 if from RAW")
options.register ('testTag',
"0",
VarParsing.VarParsing.multiplicity.singleton, # singleton or list
VarParsing.VarParsing.varType.int, # string, int, or float
"=1 if test tag to be used")
options.parseArguments()
#
process.options = cms.untracked.PSet(
wantSummary = cms.untracked.bool(True),
fileMode = cms.untracked.string("FULLMERGE")
)
process.load("FWCore.MessageService.MessageLogger_cfi")
process.MessageLogger.cout.enable = cms.untracked.bool(True)
process.MessageLogger.cout.threshold = cms.untracked.string("INFO")
process.MessageLogger.cout.default = cms.untracked.PSet(
limit = cms.untracked.int32(10000000)
)
process.MessageLogger.cout.FwkReport = cms.untracked.PSet(
reportEvery = cms.untracked.int32(10000)
)
process.MessageLogger.cerr.enable = cms.untracked.bool(True)
process.MessageLogger.cerr.threshold = cms.untracked.string("WARNING")
process.MessageLogger.cerr.default = cms.untracked.PSet(
limit = cms.untracked.int32(10000000)
)
process.MessageLogger.cerr.FwkReport = cms.untracked.PSet(
reportEvery = cms.untracked.int32(100000)
)
process.MessageLogger.debugModules = cms.untracked.vstring("clustsummmultprod","spclustermultprod","ssclustermultprod")
#----Remove too verbose PrimaryVertexProducer
process.MessageLogger.suppressInfo.append("pixelVerticesAdaptive")
process.MessageLogger.suppressInfo.append("pixelVerticesAdaptiveNoBS")
#----Remove too verbose BeamSpotOnlineProducer
process.MessageLogger.suppressInfo.append("testBeamSpot")
process.MessageLogger.suppressInfo.append("onlineBeamSpot")
process.MessageLogger.suppressWarning.append("testBeamSpot")
process.MessageLogger.suppressWarning.append("onlineBeamSpot")
#----Remove too verbose TrackRefitter
process.MessageLogger.suppressInfo.append("newTracksFromV0")
process.MessageLogger.suppressInfo.append("newTracksFromOtobV0")
#------------------------------------------------------------------
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(options.maxEvents) )
process.source = cms.Source("PoolSource",
fileNames = cms.untracked.vstring(options.inputFiles),
# skipBadFiles = cms.untracked.bool(True),
inputCommands = cms.untracked.vstring("keep *", "drop *_MEtoEDMConverter_*_*")
)
#--------------------------------------
process.seqRECO = cms.Sequence()
if options.fromRAW == 1:
process.load("Configuration.StandardSequences.RawToDigi_Data_cff")
process.load("Configuration.StandardSequences.MagneticField_cff")
process.load("Configuration.StandardSequences.GeometryDB_cff")
process.load("Configuration.StandardSequences.Reconstruction_cff")
process.load("Configuration.StandardSequences.L1Reco_cff")
process.seqRECO = cms.Sequence(process.gtEvmDigis + process.L1Reco
+ process.siStripDigis + process.siStripZeroSuppression + process.siStripClusters
+ process.siPixelDigis + process.siPixelClusters )
#
process.load("DPGAnalysis.SiStripTools.sipixelclustermultiplicityprod_cfi")
process.load("DPGAnalysis.SiStripTools.sistripclustermultiplicityprod_cfi")
#process.spclustermultprod.withClusterSize=cms.untracked.bool(True)
#process.ssclustermultprod.withClusterSize=cms.untracked.bool(True)
process.spclustermultprodnew = process.spclustermultprod.clone(wantedSubDets = cms.VPSet(
cms.PSet(detSelection = cms.uint32(101),detLabel = cms.string("BPIX"),selection=cms.untracked.vstring("0x1e000000-0x12000000")),
cms.PSet(detSelection = cms.uint32(102),detLabel = cms.string("FPIX"),selection=cms.untracked.vstring("0x1e000000-0x14000000")),
cms.PSet(detSelection = cms.uint32(11),detLabel = cms.string("BPIXL1"),selection=cms.untracked.vstring("0x1e0f0000-0x12010000")),
cms.PSet(detSelection = cms.uint32(12),detLabel = cms.string("BPIXL2"),selection=cms.untracked.vstring("0x1e0f0000-0x12020000")),
cms.PSet(detSelection = cms.uint32(13),detLabel = cms.string("BPIXL3"),selection=cms.untracked.vstring("0x1e0f0000-0x12030000")),
cms.PSet(detSelection = cms.uint32(21),detLabel = cms.string("FPIXm"),selection=cms.untracked.vstring("0x1f800000-0x14800000")),
cms.PSet(detSelection = cms.uint32(22),detLabel = cms.string("FPIXp"),selection=cms.untracked.vstring("0x1f800000-0x15000000"))
)
)
process.ssclustermultprodnew = process.ssclustermultprod.clone(wantedSubDets = cms.VPSet(
cms.PSet(detSelection = cms.uint32(103),detLabel = cms.string("TIB"),selection=cms.untracked.vstring("0x1e000000-0x16000000")),
cms.PSet(detSelection = cms.uint32(114),detLabel = cms.string("TIDm"),selection=cms.untracked.vstring("0x1e006000-0x18002000")),
cms.PSet(detSelection = cms.uint32(124),detLabel = cms.string("TIDp"),selection=cms.untracked.vstring("0x1e006000-0x18004000")),
cms.PSet(detSelection = cms.uint32(116),detLabel = cms.string("TECm"),selection=cms.untracked.vstring("0x1e0c0000-0x1c040000")),
cms.PSet(detSelection = cms.uint32(126),detLabel = cms.string("TECp"),selection=cms.untracked.vstring("0x1e0c0000-0x1c080000")),
cms.PSet(detSelection = cms.uint32(999),detLabel = cms.string("Module"),selection=cms.untracked.vstring("0x1fffffff-0x1600546d"))
)
)
process.load("DPGAnalysis.SiStripTools.clustersummarymultiplicityprod_cfi")
process.clustsummmultprod.wantedSubDets = cms.VPSet(
cms.PSet(detSelection = cms.uint32(0),detLabel = cms.string("TK"), subDetEnum = cms.int32(0) ),
cms.PSet(detSelection = cms.uint32(1),detLabel = cms.string("TIB"), subDetEnum = cms.int32(1) ),
cms.PSet(detSelection = cms.uint32(2),detLabel = cms.string("TOB"), subDetEnum = cms.int32(2) ),
cms.PSet(detSelection = cms.uint32(3),detLabel = cms.string("TID"), subDetEnum = cms.int32(3) ),
cms.PSet(detSelection = cms.uint32(4),detLabel = cms.string("TEC"), subDetEnum = cms.int32(4) ),
cms.PSet(detSelection = cms.uint32(1005),detLabel = cms.string("Pixel"), subDetEnum = cms.int32(5) ),
cms.PSet(detSelection = cms.uint32(1006),detLabel = cms.string("FPIX"), subDetEnum = cms.int32(6) ),
cms.PSet(detSelection = cms.uint32(1007),detLabel = cms.string("BPIX"), subDetEnum = cms.int32(7) ),
)
#process.ssclustermultprod.clusterdigiCollection = cms.InputTag("calZeroBiasClusters")
#process.ssclustermultprodnew.clusterdigiCollection = cms.InputTag("calZeroBiasClusters")
#process.seqMultProd = cms.Sequence(process.spclustermultprod+process.ssclustermultprod
# +process.spclustermultprodnew+process.ssclustermultprodnew
# +process.clustsummmultprod)
process.seqMultProd = cms.Sequence(process.ssclustermultprod+process.ssclustermultprodnew)
#process.seqMultProd = cms.Sequence(process.clustsummmultprod)
process.goodVertices = cms.EDFilter("VertexSelector",
src = cms.InputTag("offlinePrimaryVertices"),
cut = cms.string("!isFake && ndof > 4 && abs(z) <= 24 && position.Rho <= 2"),
filter = cms.bool(False), # otherwise it won't filter the events, just produce an empty vertex collection.
)
process.load("RecoLocalTracker.SiStripClusterizer.SiStripClusterToDigiProducer_cfi")
process.seqProducers = cms.Sequence(process.seqRECO + process.seqMultProd + process.siStripClustersToDigis
#+process.goodVertices
)
from HLTrigger.HLTfilters.triggerResultsFilter_cfi import *
process.hltSelection = triggerResultsFilter.clone(
triggerConditions = cms.vstring("HLT_ZeroBias_*"),
hltResults = cms.InputTag( "TriggerResults", "", "HLT" ),
l1tResults = cms.InputTag( "" ),
throw = cms.bool(False)
)
process.manystripclus = cms.EDFilter('BySiPixelClusterVsSiStripClusterMultiplicityEventFilter',
multiplicityConfig = cms.PSet(
firstMultiplicityConfig = cms.PSet(
collectionName = cms.InputTag("siPixelClusters"),
moduleThreshold = cms.untracked.int32(-1),
useQuality = cms.untracked.bool(False),
qualityLabel = cms.untracked.string("")
),
secondMultiplicityConfig = cms.PSet(
collectionName = cms.InputTag("siStripClusters"),
moduleThreshold = cms.untracked.int32(-1),
useQuality = cms.untracked.bool(False),
qualityLabel = cms.untracked.string("")
)
),
cut = cms.string("( mult2 > 20000+7*mult1)")
)
process.manystripclus53X = cms.EDFilter('ByClusterSummaryMultiplicityPairEventFilter',
multiplicityConfig = cms.PSet(
firstMultiplicityConfig = cms.PSet(
clusterSummaryCollection = cms.InputTag("clusterSummaryProducer"),
subDetEnum = cms.int32(5),
subDetVariable = cms.string("pHits")
),
secondMultiplicityConfig = cms.PSet(
clusterSummaryCollection = cms.InputTag("clusterSummaryProducer"),
subDetEnum = cms.int32(0),
subDetVariable = cms.string("cHits")
),
),
cut = cms.string("( mult2 > 20000+7*mult1)")
)
process.load("DPGAnalysis.SiStripTools.ssclusmultinvestigator_cfi")
#process.ssclusmultinvestigator.vertexCollection = cms.InputTag("goodVertices")
#process.ssclusmultinvestigator.wantInvestHist = cms.bool(True)
#process.ssclusmultinvestigator.wantVtxCorrHist = cms.bool(True)
#process.ssclusmultinvestigator.digiVtxCorrConfig = cms.PSet(
# wantedSubDets = process.ssclusmultinvestigator.wantedSubDets,
# hitName = cms.untracked.string("cluster"),
# numberOfBins = cms.untracked.int32(100),
# scaleFactor = cms.untracked.int32(1),
# maxNvtx = cms.untracked.int32(100)
# )
#process.ssclusmultinvestigator.wantLumiCorrHist = cms.bool(True)
#process.ssclusmultinvestigator.digiLumiCorrConfig = cms.PSet(
# lumiProducer = cms.InputTag("lumiProducer"),
# wantedSubDets = process.ssclusmultinvestigator.wantedSubDets,
# hitName = cms.untracked.string("cluster"),
# numberOfBins = cms.untracked.int32(100),
# scaleFactor = cms.untracked.int32(1),
# maxLumi = cms.untracked.double(25)
# )
#process.ssclusmultinvestigator.runHisto = cms.untracked.bool(True)
process.ssclusmultinvestigator.scaleFactor=cms.untracked.int32(1)
process.ssclusmultinvestigatornew = process.ssclusmultinvestigator.clone()
process.ssclusmultinvestigatornew.runHisto = cms.untracked.bool(True)
process.ssclusmultinvestigatornew.wantedSubDets = cms.untracked.VPSet(
cms.PSet(detSelection = cms.uint32(103),detLabel = cms.string("TIB"), binMax = cms.int32(1787904/32)),
cms.PSet(detSelection = cms.uint32(114),detLabel = cms.string("TIDm"), binMax = cms.int32(565248/32)),
cms.PSet(detSelection = cms.uint32(124),detLabel = cms.string("TIDp"), binMax = cms.int32(565248/32)),
cms.PSet(detSelection = cms.uint32(116),detLabel = cms.string("TECm"), binMax = cms.int32(3866624/64)),
cms.PSet(detSelection = cms.uint32(126),detLabel = cms.string("TECp"), binMax = cms.int32(3866624/64)),
cms.PSet(detSelection = cms.uint32(999),detLabel = cms.string("Module"), binMax = cms.int32(200))
)
#process.ssclusmultinvestigatornew.digiVtxCorrConfig.wantedSubDets = cms.untracked.VPSet(
# cms.PSet(detSelection = cms.uint32(103),detLabel = cms.string("TIB"), binMax = cms.int32(1787904/32)),
# cms.PSet(detSelection = cms.uint32(114),detLabel = cms.string("TIDm"), binMax = cms.int32(565248/32)),
# cms.PSet(detSelection = cms.uint32(124),detLabel = cms.string("TIDp"), binMax = cms.int32(565248/32)),
# cms.PSet(detSelection = cms.uint32(116),detLabel = cms.string("TECm"), binMax = cms.int32(3866624/64)),
# cms.PSet(detSelection = cms.uint32(126),detLabel = cms.string("TECp"), binMax = cms.int32(3866624/64)),
# cms.PSet(detSelection = cms.uint32(999),detLabel = cms.string("Module"), binMax = cms.int32(200))
# )
#process.ssclusmultinvestigatornew.digiLumiCorrConfig.wantedSubDets = cms.untracked.VPSet(
# cms.PSet(detSelection = cms.uint32(103),detLabel = cms.string("TIB"), binMax = cms.int32(1787904/32)),
# cms.PSet(detSelection = cms.uint32(114),detLabel = cms.string("TIDm"), binMax = cms.int32(565248/32)),
# cms.PSet(detSelection = cms.uint32(124),detLabel = cms.string("TIDp"), binMax = cms.int32(565248/32)),
# cms.PSet(detSelection = cms.uint32(116),detLabel = cms.string("TECm"), binMax = cms.int32(3866624/64)),
# cms.PSet(detSelection = cms.uint32(126),detLabel = cms.string("TECp"), binMax = cms.int32(3866624/64)),
# cms.PSet(detSelection = cms.uint32(999),detLabel = cms.string("Module"), binMax = cms.int32(200))
# )
process.ssclusmultinvestigatornew.multiplicityMap = cms.InputTag("ssclustermultprodnew")
process.load("DPGAnalysis.SiStripTools.spclusmultinvestigator_cfi")
process.spclusmultinvestigator.vertexCollection = cms.InputTag("goodVertices")
process.spclusmultinvestigator.wantInvestHist = cms.bool(True)
process.spclusmultinvestigator.wantVtxCorrHist = cms.bool(True)
process.spclusmultinvestigator.digiVtxCorrConfig = cms.PSet(
wantedSubDets = process.spclusmultinvestigator.wantedSubDets,
hitName = cms.untracked.string("cluster"),
numberOfBins = cms.untracked.int32(100),
scaleFactor = cms.untracked.int32(5)
)
process.spclusmultinvestigator.wantLumiCorrHist = cms.bool(True)
process.spclusmultinvestigator.digiLumiCorrConfig = cms.PSet(
lumiProducer = cms.InputTag("lumiProducer"),
wantedSubDets = process.spclusmultinvestigator.wantedSubDets,
hitName = cms.untracked.string("cluster"),
numberOfBins = cms.untracked.int32(100),
scaleFactor = cms.untracked.int32(5)
)
#process.spclusmultinvestigator.runHisto = cms.untracked.bool(True)
process.spclusmultinvestigator.scaleFactor=cms.untracked.int32(5)
process.spclusmultinvestigatornew = process.spclusmultinvestigator.clone()
process.spclusmultinvestigatornew.wantedSubDets = cms.untracked.VPSet(
cms.PSet(detSelection = cms.uint32(101),detLabel = cms.string("BPIX"), binMax = cms.int32(100000)),
cms.PSet(detSelection = cms.uint32(11),detLabel = cms.string("BPIXL1"), binMax = cms.int32(100000)),
cms.PSet(detSelection = cms.uint32(12),detLabel = cms.string("BPIXL2"), binMax = cms.int32(100000)),
cms.PSet(detSelection = cms.uint32(13),detLabel = cms.string("BPIXL3"), binMax = cms.int32(100000)),
cms.PSet(detSelection = cms.uint32(102),detLabel = cms.string("FPIX"), binMax = cms.int32(100000)),
cms.PSet(detSelection = cms.uint32(21),detLabel = cms.string("FPIXm"), binMax = cms.int32(100000)),
cms.PSet(detSelection = cms.uint32(22),detLabel = cms.string("FPIXp"), binMax = cms.int32(100000))
)
process.spclusmultinvestigatornew.digiVtxCorrConfig.wantedSubDets = cms.untracked.VPSet(
cms.PSet(detSelection = cms.uint32(101),detLabel = cms.string("BPIX"), binMax = cms.int32(100000)),
cms.PSet(detSelection = cms.uint32(11),detLabel = cms.string("BPIXL1"), binMax = cms.int32(100000)),
cms.PSet(detSelection = cms.uint32(12),detLabel = cms.string("BPIXL2"), binMax = cms.int32(100000)),
cms.PSet(detSelection = cms.uint32(13),detLabel = cms.string("BPIXL3"), binMax = cms.int32(100000)),
cms.PSet(detSelection = cms.uint32(102),detLabel = cms.string("FPIX"), binMax = cms.int32(100000)),
cms.PSet(detSelection = cms.uint32(21),detLabel = cms.string("FPIXm"), binMax = cms.int32(100000)),
cms.PSet(detSelection = cms.uint32(22),detLabel = cms.string("FPIXp"), binMax = cms.int32(100000))
)
process.spclusmultinvestigatornew.digiLumiCorrConfig.wantedSubDets = cms.untracked.VPSet(
cms.PSet(detSelection = cms.uint32(101),detLabel = cms.string("BPIX"), binMax = cms.int32(100000)),
cms.PSet(detSelection = cms.uint32(11),detLabel = cms.string("BPIXL1"), binMax = cms.int32(100000)),
cms.PSet(detSelection = cms.uint32(12),detLabel = cms.string("BPIXL2"), binMax = cms.int32(100000)),
cms.PSet(detSelection = cms.uint32(13),detLabel = cms.string("BPIXL3"), binMax = cms.int32(100000)),
cms.PSet(detSelection = cms.uint32(102),detLabel = cms.string("FPIX"), binMax = cms.int32(100000)),
cms.PSet(detSelection = cms.uint32(21),detLabel = cms.string("FPIXm"), binMax = cms.int32(100000)),
cms.PSet(detSelection = cms.uint32(22),detLabel = cms.string("FPIXp"), binMax = cms.int32(100000))
)
process.spclusmultinvestigatornew.multiplicityMap = cms.InputTag("spclustermultprodnew")
process.load("DPGAnalysis.SiStripTools.multiplicitycorr_cfi")
process.multiplicitycorr.correlationConfigurations = cms.VPSet(
cms.PSet(xMultiplicityMap = cms.InputTag("spclustermultprodnew"),
xDetSelection = cms.uint32(12), xDetLabel = cms.string("BPixL2"), xBins = cms.uint32(1000), xMax=cms.double(10000),
yMultiplicityMap = cms.InputTag("ssclustermultprodnew"),
yDetSelection = cms.uint32(103), yDetLabel = cms.string("TIB"), yBins = cms.uint32(1000), yMax=cms.double(50000),
rBins = cms.uint32(200), scaleFactor = cms.untracked.double(.25),
runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False)),
cms.PSet(xMultiplicityMap = cms.InputTag("spclustermultprodnew"),
xDetSelection = cms.uint32(12), xDetLabel = cms.string("BPixL2"), xBins = cms.uint32(1000), xMax=cms.double(10000),
yMultiplicityMap = cms.InputTag("ssclustermultprodnew"),
yDetSelection = cms.uint32(114), yDetLabel = cms.string("TIDm"), yBins = cms.uint32(1000), yMax=cms.double(20000),
rBins = cms.uint32(200), scaleFactor = cms.untracked.double(.5),
runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False)),
cms.PSet(xMultiplicityMap = cms.InputTag("spclustermultprodnew"),
xDetSelection = cms.uint32(12), xDetLabel = cms.string("BPixL2"), xBins = cms.uint32(1000), xMax=cms.double(10000),
yMultiplicityMap = cms.InputTag("ssclustermultprodnew"),
yDetSelection = cms.uint32(124), yDetLabel = cms.string("TIDp"), yBins = cms.uint32(1000), yMax=cms.double(20000),
rBins = cms.uint32(200), scaleFactor = cms.untracked.double(.5),
runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False)),
cms.PSet(xMultiplicityMap = cms.InputTag("spclustermultprodnew"),
xDetSelection = cms.uint32(12), xDetLabel = cms.string("BPixL2"), xBins = cms.uint32(1000), xMax=cms.double(10000),
yMultiplicityMap = cms.InputTag("spclustermultprodnew"),
yDetSelection = cms.uint32(11), yDetLabel = cms.string("BPixL1"), yBins = cms.uint32(1000), yMax=cms.double(10000),
rBins = cms.uint32(200), scaleFactor = cms.untracked.double(1.),
runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False))
)
process.multcorrclustsumm = process.multiplicitycorr.clone()
process.multcorrclustsumm.correlationConfigurations = cms.VPSet(
cms.PSet(yMultiplicityMap = cms.InputTag("spclustermultprod"),
yDetSelection = cms.uint32(0), yDetLabel = cms.string("Pixel"), yBins = cms.uint32(1000), yMax=cms.double(1000),
xMultiplicityMap = cms.InputTag("clustsummmultprod"),
xDetSelection = cms.uint32(1005), xDetLabel = cms.string("Pixelcs"), xBins = cms.uint32(1000), xMax=cms.double(1000),
rBins = cms.uint32(200), scaleFactor = cms.untracked.double(1.),
runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False)),
cms.PSet(xMultiplicityMap = cms.InputTag("spclustermultprodnew"),
xDetSelection = cms.uint32(101), xDetLabel = cms.string("BPIX"), xBins = cms.uint32(1000), xMax=cms.double(1000),
yMultiplicityMap = cms.InputTag("clustsummmultprod"),
yDetSelection = cms.uint32(1007), yDetLabel = cms.string("BPIXcs"), yBins = cms.uint32(1000), yMax=cms.double(1000),
rBins = cms.uint32(200), scaleFactor = cms.untracked.double(1.),
runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False)),
cms.PSet(xMultiplicityMap = cms.InputTag("spclustermultprodnew"),
xDetSelection = cms.uint32(102), xDetLabel = cms.string("FPIX"), xBins = cms.uint32(1000), xMax=cms.double(1000),
yMultiplicityMap = cms.InputTag("clustsummmultprod"),
yDetSelection = cms.uint32(1006), yDetLabel = cms.string("FPIXcs"), yBins = cms.uint32(1000), yMax=cms.double(1000),
rBins = cms.uint32(200), scaleFactor = cms.untracked.double(1.),
runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False)),
cms.PSet(xMultiplicityMap = cms.InputTag("ssclustermultprod"),
xDetSelection = cms.uint32(3), xDetLabel = cms.string("TIB"), xBins = cms.uint32(1000), xMax=cms.double(50000),
yMultiplicityMap = cms.InputTag("clustsummmultprod"),
yDetSelection = cms.uint32(1), yDetLabel = cms.string("TIBcs"), yBins = cms.uint32(1000), yMax=cms.double(50000),
rBins = cms.uint32(200), scaleFactor = cms.untracked.double(1.),
runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False)),
cms.PSet(xMultiplicityMap = cms.InputTag("ssclustermultprod"),
xDetSelection = cms.uint32(4), xDetLabel = cms.string("TID"), xBins = cms.uint32(1000), xMax=cms.double(50000),
yMultiplicityMap = cms.InputTag("clustsummmultprod"),
yDetSelection = cms.uint32(3), yDetLabel = cms.string("TIDcs"), yBins = cms.uint32(1000), yMax=cms.double(50000),
rBins = cms.uint32(200), scaleFactor = cms.untracked.double(1.),
runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False)),
cms.PSet(xMultiplicityMap = cms.InputTag("ssclustermultprod"),
xDetSelection = cms.uint32(5), xDetLabel = cms.string("TOB"), xBins = cms.uint32(1000), xMax=cms.double(50000),
yMultiplicityMap = cms.InputTag("clustsummmultprod"),
yDetSelection = cms.uint32(2), yDetLabel = cms.string("TOBcs"), yBins = cms.uint32(1000), yMax=cms.double(50000),
rBins = cms.uint32(200), scaleFactor = cms.untracked.double(1.),
runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False)),
cms.PSet(xMultiplicityMap = cms.InputTag("ssclustermultprod"),
xDetSelection = cms.uint32(6), xDetLabel = cms.string("TEC"), xBins = cms.uint32(1000), xMax=cms.double(50000),
yMultiplicityMap = cms.InputTag("clustsummmultprod"),
yDetSelection = cms.uint32(4), yDetLabel = cms.string("TECcs"), yBins = cms.uint32(1000), yMax=cms.double(50000),
rBins = cms.uint32(200), scaleFactor = cms.untracked.double(1.),
runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False))
)
process.clustsummmultcorr = cms.EDAnalyzer('MultiplicityCorrelator',
correlationConfigurations = cms.VPSet(
cms.PSet(xMultiplicityMap = cms.InputTag("clustsummmultprod"), xDetSelection = cms.uint32(0), xDetLabel = cms.string("TK"), xBins = cms.uint32(1000), xMax=cms.double(150000),
yMultiplicityMap = cms.InputTag("clustsummmultprod"), yDetSelection = cms.uint32(1005), yDetLabel = cms.string("Pixel"), yBins = cms.uint32(1000), yMax=cms.double(20000),
rBins = cms.uint32(200), scaleFactor = cms.untracked.double(5.),
runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False))
)
)
#process.clustsummmultcorr = cms.EDAnalyzer('MultiplicityCorrelator',
# correlationConfigurations = cms.VPSet(
# cms.PSet(xMultiplicityMap = cms.InputTag("ssclustermultprod"), xDetSelection = cms.uint32(0), xDetLabel = cms.string("TK"), xBins = cms.uint32(1000), xMax=cms.double(150000),
# yMultiplicityMap = cms.InputTag("spclustermultprod"), yDetSelection = cms.uint32(0), yDetLabel = cms.string("Pixel"), yBins = cms.uint32(1000), yMax=cms.double(20000),
# rBins = cms.uint32(200), scaleFactor = cms.untracked.double(5.),
# runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False))
# )
# )
process.clustsummmultcorrmanystripclus = process.clustsummmultcorr.clone()
process.load("TrackingPFG.Utilities.logerroranalyzer_cfi")
process.load("DPGAnalysis.SiStripTools.clusterbigeventsdebugger_cfi")
process.clusterbigeventsdebugger.selections.append(cms.PSet(label=cms.string("Module"),selection=cms.untracked.vstring("0x1fffffff-0x1600546d")))
#process.clusterbigeventsdebugger.collection = cms.InputTag("calZeroBiasClusters")
#process.clusterbigeventsdebugger.foldedStrips = cms.untracked.bool(True)
#process.clusterbigeventsdebugger.singleEvents = cms.bool(True)
process.load("DPGAnalysis.SiStripTools.digibigeventsdebugger_cfi")
process.digibigeventsdebugger.selections = cms.VPSet(
cms.PSet(label=cms.string("TIB"),selection=cms.untracked.vstring("0x1e000000-0x16000000")),
cms.PSet(label=cms.string("TIBD"),selection=cms.untracked.vstring("0x1e000000-0x16000000","0x1e000000-0x18000000")),
cms.PSet(label=cms.string("TECp"),selection=cms.untracked.vstring("0x1e0c0000-0x1c080000")),
cms.PSet(label=cms.string("TECm"),selection=cms.untracked.vstring("0x1e0c0000-0x1c040000")),
cms.PSet(label=cms.string("TOB"),selection=cms.untracked.vstring("0x1e000000-0x1a000000")),
cms.PSet(label=cms.string("TID"),selection=cms.untracked.vstring("0x1e000000-0x18000000"))
)
process.digibigeventsdebugger.selections.append(cms.PSet(label=cms.string("Module"),selection=cms.untracked.vstring("0x1fffffff-0x1600546d")))
process.digibigeventsdebugger.collection = cms.InputTag("siStripClustersToDigis","ZeroSuppressed")
process.digibigeventsdebugger.foldedStrips = cms.untracked.bool(True)
#process.digibigeventsdebugger.singleEvents = cms.bool(True)
process.load("DPGAnalysis.SiStripTools.fedbadmodulefilter_cfi")
process.fedbadmodulefilter.badModThr = 0
process.fedbadmodulefilter.wantedHisto = cms.untracked.bool(True)
process.fedbadmodulefilter.moduleSelection = cms.untracked.PSet(selection=cms.untracked.vstring("0x1fffffff-0x1600546d"))
process.fedbadmoduleTECmD8 = process.fedbadmodulefilter.clone(moduleSelection = cms.untracked.PSet(selection=cms.untracked.vstring("0x1e0fc000-0x1c060000")))
process.fedbadmoduleTECmD4 = process.fedbadmodulefilter.clone(moduleSelection = cms.untracked.PSet(selection=cms.untracked.vstring("0x1e0fc000-0x1c050000")))
process.fedbadmoduleTECmD7 = process.fedbadmodulefilter.clone(moduleSelection = cms.untracked.PSet(selection=cms.untracked.vstring("0x1e0fc000-0x1c05c000")))
#process.seqClusMultInvest = cms.Sequence(process.spclusmultinvestigator + process.ssclusmultinvestigator +
# process.spclusmultinvestigatornew + process.ssclusmultinvestigatornew
# + process.multiplicitycorr + process.multcorrclustsumm)
#process.seqClusMultInvest = cms.Sequence(process.clustsummmultcorr)
process.seqClusMultInvest = cms.Sequence(process.ssclusmultinvestigatornew +
process.clusterbigeventsdebugger + process.digibigeventsdebugger +
process.fedbadmodulefilter + process.fedbadmoduleTECmD8 + process.fedbadmoduleTECmD4 + process.fedbadmoduleTECmD7
)
process.p0 = cms.Path(
# process.hltSelection +
process.seqProducers +
process.seqClusMultInvest)
#process.pmanystripclus = cms.Path(process.seqProducers +
# process.manystripclus53X +
## process.manystripclus +
# process.clustsummmultcorrmanystripclus +
# process.logerroranalyzer
# )
#process.outoddevent = cms.OutputModule("PoolOutputModule",
# fileName = cms.untracked.string("manystripclusevents.root"),
# outputCommands = cms.untracked.vstring("keep *"),
# SelectEvents = cms.untracked.PSet(SelectEvents = cms.vstring( "pmanystripclus"))
# )
#process.e = cms.EndPath(process.outoddevent)
#----GlobalTag ------------------------
process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
from Configuration.AlCa.GlobalTag import GlobalTag
process.GlobalTag = GlobalTag(process.GlobalTag, options.globalTag, '')
if options.testTag==1:
process.stripConditions = cms.ESSource(
"PoolDBESSource",
process.CondDBSetup,
connect = cms.string('frontier://FrontierProd/CMS_COND_31X_STRIP'),
toGet = cms.VPSet(
cms.PSet(
record = cms.string('SiStripNoisesRcd'), tag = cms.string('SiStripNoise_GR10_v2_offline')
),
),
)
process.es_prefer_strips = cms.ESPrefer("PoolDBESSource","stripConditions")
process.TFileService = cms.Service('TFileService',
fileName = cms.string('MultiplicityProducerTest.root')
)
#print process.dumpPython()
|