SiStripTools/test/MultiplicityMonitor_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002 import FWCore.ParameterSet.VarParsing as VarParsing
0003
0004 process = cms.Process("MultiplicityMonitor")
0005
0006 #prepare options
0007
0008 options = VarParsing.VarParsing("analysis")
0009
0010 options.register ('globalTag',
0011                   "DONOTEXIST",
0012                   VarParsing.VarParsing.multiplicity.singleton, # singleton or list
0013                   VarParsing.VarParsing.varType.string,          # string, int, or float
0014                   "GlobalTag")
0015 options.register ('fromRAW',
0016                   "0",
0017                   VarParsing.VarParsing.multiplicity.singleton, # singleton or list
0018                   VarParsing.VarParsing.varType.int,          # string, int, or float
0019                   "=1 if from RAW")
0020 options.register ('testTag',
0021                   "0",
0022                   VarParsing.VarParsing.multiplicity.singleton, # singleton or list
0023                   VarParsing.VarParsing.varType.int,          # string, int, or float
0024                   "=1 if test tag to be used")
0025
0026 options.parseArguments()
0027
0028 #
0029
0030 process.options = cms.untracked.PSet(
0031     wantSummary = cms.untracked.bool(True),
0032     fileMode = cms.untracked.string("FULLMERGE")
0033     )
0034
0035 process.load("FWCore.MessageService.MessageLogger_cfi")
0036
0037 process.MessageLogger.cout.enable = cms.untracked.bool(True)
0038 process.MessageLogger.cout.threshold = cms.untracked.string("INFO")
0039 process.MessageLogger.cout.default = cms.untracked.PSet(
0040     limit = cms.untracked.int32(10000000)
0041     )
0042 process.MessageLogger.cout.FwkReport = cms.untracked.PSet(
0043     reportEvery = cms.untracked.int32(10000)
0044     )
0045
0046 process.MessageLogger.cerr.enable = cms.untracked.bool(True)
0047 process.MessageLogger.cerr.threshold = cms.untracked.string("WARNING")
0048 process.MessageLogger.cerr.default = cms.untracked.PSet(
0049     limit = cms.untracked.int32(10000000)
0050     )
0051 process.MessageLogger.cerr.FwkReport = cms.untracked.PSet(
0052     reportEvery = cms.untracked.int32(100000)
0053     )
0054
0055 #------------------------------------------------------------------
0056
0057 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(options.maxEvents) )
0058
0059 process.source = cms.Source("PoolSource",
0060                     fileNames = cms.untracked.vstring(options.inputFiles),
0061 #                    skipBadFiles = cms.untracked.bool(True),
0062                     inputCommands = cms.untracked.vstring("keep *", "drop *_MEtoEDMConverter_*_*")
0063                     )
0064
0065 #--------------------------------------
0066 process.seqRECO = cms.Sequence()
0067
0068 if options.fromRAW == 1:
0069     process.load("Configuration.StandardSequences.RawToDigi_Data_cff")
0070     process.load("Configuration.StandardSequences.MagneticField_cff")
0071     process.load("Configuration.StandardSequences.GeometryDB_cff")
0072     process.load("Configuration.StandardSequences.Reconstruction_cff")
0073     process.load("Configuration.StandardSequences.L1Reco_cff")
0074     process.seqRECO = cms.Sequence(process.gtEvmDigis + process.L1Reco
0075                                    + process.siStripDigis + process.siStripZeroSuppression + process.siStripClusters
0076                                    + process.siPixelDigis + process.siPixelClusters )
0077
0078
0079 #
0080
0081 process.load("DPGAnalysis.SiStripTools.sipixelclustermultiplicityprod_cfi")
0082 process.load("DPGAnalysis.SiStripTools.sistripclustermultiplicityprod_cfi")
0083
0084 #process.spclustermultprod.wantedSubDets.append = cms.VPSet(
0085 #    cms.PSet(detSelection = cms.uint32(101),detLabel = cms.string("BPIX"),selection=cms.untracked.vstring("0x1e000000-0x12000000")),
0086 #    cms.PSet(detSelection = cms.uint32(102),detLabel = cms.string("FPIX"),selection=cms.untracked.vstring("0x1e000000-0x14000000"))
0087 #)
0088
0089 #process.ssclustermultprod.wantedSubDets.append = cms.VPSet(
0090 #    cms.PSet(detSelection = cms.uint32(103),detLabel = cms.string("TIB"),selection=cms.untracked.vstring("0x1e000000-0x16000000")),
0091 #    cms.PSet(detSelection = cms.uint32(105),detLabel = cms.string("TOB"),selection=cms.untracked.vstring("0x1e000000-0x1a000000")),
0092 #    cms.PSet(detSelection = cms.uint32(114),detLabel = cms.string("TIDm"),selection=cms.untracked.vstring("0x1e006000-0x18002000")),
0093 #    cms.PSet(detSelection = cms.uint32(124),detLabel = cms.string("TIDp"),selection=cms.untracked.vstring("0x1e006000-0x18004000")),
0094 #    cms.PSet(detSelection = cms.uint32(116),detLabel = cms.string("TECm"),selection=cms.untracked.vstring("0x1e0c0000-0x1c040000")),
0095 #    cms.PSet(detSelection = cms.uint32(126),detLabel = cms.string("TECp"),selection=cms.untracked.vstring("0x1e0c0000-0x1c080000"))
0096 #)
0097
0098
0099
0100 process.seqMultProd = cms.Sequence(process.spclustermultprod+process.ssclustermultprod)
0101
0102 process.goodVertices = cms.EDFilter("VertexSelector",
0103    src = cms.InputTag("offlinePrimaryVertices"),
0104    cut = cms.string("!isFake && ndof > 4 && abs(z) <= 24 && position.Rho <= 2"),
0105    filter = cms.bool(False),   # otherwise it won't filter the events, just produce an empty vertex collection.
0106 )
0107
0108
0109
0110 process.seqProducers = cms.Sequence(process.seqRECO + process.seqMultProd)
0111                                     #+process.goodVertices)
0112
0113 from HLTrigger.HLTfilters.triggerResultsFilter_cfi import *
0114 process.hltSelection = triggerResultsFilter.clone(
0115                                           triggerConditions = cms.vstring("HLT_ZeroBias_*"),
0116                                           hltResults = cms.InputTag( "TriggerResults", "", "HLT" ),
0117                                           l1tResults = cms.InputTag( "" ),
0118                                           throw = cms.bool(False)
0119                                           )
0120
0121 process.manypixelclus = cms.EDFilter('BySiPixelClusterVsSiStripClusterMultiplicityEventFilter',
0122                                      multiplicityConfig = cms.PSet(
0123     firstMultiplicityConfig = cms.PSet(
0124     collectionName = cms.InputTag("siPixelClusters"),
0125     moduleThreshold = cms.untracked.int32(-1),
0126     useQuality = cms.untracked.bool(False),
0127     qualityLabel = cms.untracked.string("")
0128     ),
0129     secondMultiplicityConfig = cms.PSet(
0130     collectionName = cms.InputTag("siStripClusters"),
0131     moduleThreshold = cms.untracked.int32(-1),
0132     useQuality = cms.untracked.bool(False),
0133     qualityLabel = cms.untracked.string("")
0134     )
0135     ),
0136                                      cut = cms.string("( mult1 > 5000+0.1*mult2)")
0137                                      )
0138
0139 process.load("DPGAnalysis.SiStripTools.ssclusmultinvestigator_cfi")
0140
0141 process.ssclusmultinvestigator.scaleFactor=cms.untracked.int32(1)
0142
0143 process.ssclusmultinvestigator.runHisto = cms.untracked.bool(True)
0144 #process.ssclusmultinvestigator.wantedSubDets.append = cms.untracked.VPSet(
0145 #    cms.PSet(detSelection = cms.uint32(103),detLabel = cms.string("TIB"), binMax = cms.int32(1787904/32)),
0146 #    cms.PSet(detSelection = cms.uint32(105),detLabel = cms.string("TOB"), binMax = cms.int32(1787904/32)),
0147 #    cms.PSet(detSelection = cms.uint32(114),detLabel = cms.string("TIDm"), binMax = cms.int32(565248/32)),
0148 #    cms.PSet(detSelection = cms.uint32(124),detLabel = cms.string("TIDp"), binMax = cms.int32(565248/32)),
0149 #    cms.PSet(detSelection = cms.uint32(116),detLabel = cms.string("TECm"), binMax = cms.int32(3866624/64)),
0150 #    cms.PSet(detSelection = cms.uint32(126),detLabel = cms.string("TECp"), binMax = cms.int32(3866624/64))
0151 #    )
0152
0153 process.load("DPGAnalysis.SiStripTools.spclusmultinvestigator_cfi")
0154 process.spclusmultinvestigator.scaleFactor=cms.untracked.int32(5)
0155
0156 #process.spclusmultinvestigator.wantedSubDets.append = cms.untracked.VPSet(
0157 #    cms.PSet(detSelection = cms.uint32(101),detLabel = cms.string("BPIX"), binMax = cms.int32(100000)),
0158 #    cms.PSet(detSelection = cms.uint32(102),detLabel = cms.string("FPIX"), binMax = cms.int32(100000))
0159 #    )
0160
0161 process.load("DPGAnalysis.SiStripTools.multiplicitycorr_cfi")
0162 process.multiplicitycorr.correlationConfigurations = cms.VPSet(
0163    cms.PSet(xMultiplicityMap = cms.InputTag("spclustermultprod"),
0164             xDetSelection = cms.uint32(0), xDetLabel = cms.string("Pixel"), xBins = cms.uint32(2000), xMax=cms.double(50000),
0165             yMultiplicityMap = cms.InputTag("ssclustermultprod"),
0166             yDetSelection = cms.uint32(0), yDetLabel = cms.string("Tracker"), yBins = cms.uint32(2000), yMax=cms.double(500000),
0167             rBins = cms.uint32(200), scaleFactor = cms.untracked.double(.25),
0168             runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False)),
0169    cms.PSet(xMultiplicityMap = cms.InputTag("spclustermultprod"),
0170             xDetSelection = cms.uint32(0), xDetLabel = cms.string("Pixel"), xBins = cms.uint32(2000), xMax=cms.double(50000),
0171             yMultiplicityMap = cms.InputTag("ssclustermultprod"),
0172             yDetSelection = cms.uint32(3), yDetLabel = cms.string("TIB"), yBins = cms.uint32(2000), yMax=cms.double(500000),
0173             rBins = cms.uint32(200), scaleFactor = cms.untracked.double(.25),
0174             runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False)),
0175    cms.PSet(xMultiplicityMap = cms.InputTag("spclustermultprod"),
0176             xDetSelection = cms.uint32(0), xDetLabel = cms.string("Pixel"), xBins = cms.uint32(2000), xMax=cms.double(50000),
0177             yMultiplicityMap = cms.InputTag("ssclustermultprod"),
0178             yDetSelection = cms.uint32(5), yDetLabel = cms.string("TOB"), yBins = cms.uint32(2000), yMax=cms.double(500000),
0179             rBins = cms.uint32(200), scaleFactor = cms.untracked.double(.25),
0180             runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False)),
0181    )
0182
0183 process.multcorrmanypixel = process.multiplicitycorr.clone()
0184
0185 process.seqClusMultInvest = cms.Sequence(process.ssclusmultinvestigator +process.spclusmultinvestigator + process.multiplicitycorr)
0186
0187
0188
0189
0190 process.p0 = cms.Path(
0191     #    process.hltSelection +
0192     process.seqProducers +
0193     process.seqClusMultInvest +
0194     process.manypixelclus +
0195     process.multcorrmanypixel)
0196
0197 #----GlobalTag ------------------------
0198
0199 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0200 from Configuration.AlCa.GlobalTag import GlobalTag
0201 process.GlobalTag = GlobalTag(process.GlobalTag, options.globalTag, '')
0202
0203
0204 process.TFileService = cms.Service('TFileService',
0205                                    fileName = cms.string('MultiplicityMonitor.root')
0206                                    )
0207
0208 print(process.dumpPython())