SiStripTools/test/ByMultiplicityFilterTest_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002 import FWCore.ParameterSet.VarParsing as VarParsing
0003
0004 process = cms.Process("ByMultiplicityFilterTest")
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
0016 options.parseArguments()
0017
0018 #
0019
0020 process.options = cms.untracked.PSet(
0021     wantSummary = cms.untracked.bool(True),
0022     fileMode = cms.untracked.string("FULLMERGE")
0023     )
0024
0025 process.load("FWCore.MessageService.MessageLogger_cfi")
0026
0027 process.MessageLogger.cout.enable = cms.untracked.bool(True)
0028 process.MessageLogger.cout.threshold = cms.untracked.string("INFO")
0029 process.MessageLogger.cout.default = cms.untracked.PSet(
0030     limit = cms.untracked.int32(10000000)
0031     )
0032 process.MessageLogger.cout.FwkReport = cms.untracked.PSet(
0033     reportEvery = cms.untracked.int32(10000)
0034     )
0035
0036 process.MessageLogger.cerr.enable = cms.untracked.bool(True)
0037 process.MessageLogger.cerr.threshold = cms.untracked.string("WARNING")
0038 process.MessageLogger.cerr.default = cms.untracked.PSet(
0039     limit = cms.untracked.int32(10000000)
0040     )
0041 process.MessageLogger.cerr.FwkReport = cms.untracked.PSet(
0042     reportEvery = cms.untracked.int32(100000)
0043     )
0044
0045 #----Remove too verbose PrimaryVertexProducer
0046
0047 process.MessageLogger.suppressInfo.append("pixelVerticesAdaptive")
0048 process.MessageLogger.suppressInfo.append("pixelVerticesAdaptiveNoBS")
0049
0050 #----Remove too verbose BeamSpotOnlineProducer
0051
0052 process.MessageLogger.suppressInfo.append("testBeamSpot")
0053 process.MessageLogger.suppressInfo.append("onlineBeamSpot")
0054 process.MessageLogger.suppressWarning.append("testBeamSpot")
0055 process.MessageLogger.suppressWarning.append("onlineBeamSpot")
0056
0057 #----Remove too verbose TrackRefitter
0058
0059 process.MessageLogger.suppressInfo.append("newTracksFromV0")
0060 process.MessageLogger.suppressInfo.append("newTracksFromOtobV0")
0061
0062 #------------------------------------------------------------------
0063
0064 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )
0065
0066 process.source = cms.Source("PoolSource",
0067                     fileNames = cms.untracked.vstring(options.inputFiles),
0068 #                    skipBadFiles = cms.untracked.bool(True),
0069                     inputCommands = cms.untracked.vstring("keep *", "drop *_MEtoEDMConverter_*_*")
0070                     )
0071
0072 #--------------------------------------
0073
0074 process.load("DPGAnalysis.SiStripTools.sipixelclustermultiplicityprod_cfi")
0075 process.load("DPGAnalysis.SiStripTools.sistripclustermultiplicityprod_cfi")
0076 process.load("DPGAnalysis.SiStripTools.clustersummarymultiplicityprod_cfi")
0077 process.seqMultProd = cms.Sequence(process.spclustermultprod+process.ssclustermultprod + process.clustsummmultprod)
0078 #process.seqMultProd = cms.Sequence(process.clustsummmultprod)
0079
0080 process.load("DPGAnalysis.SiStripTools.multiplicitycorr_cfi")
0081 process.multiplicitycorr.correlationConfigurations = cms.VPSet(
0082    cms.PSet(xMultiplicityMap = cms.InputTag("ssclustermultprod"),
0083             xDetSelection = cms.uint32(0), xDetLabel = cms.string("TK"), xBins = cms.uint32(3000), xMax=cms.double(100000),
0084             yMultiplicityMap = cms.InputTag("spclustermultprod"),
0085             yDetSelection = cms.uint32(0), yDetLabel = cms.string("Pixel"), yBins = cms.uint32(1000), yMax=cms.double(30000),
0086             rBins = cms.uint32(200), scaleFactor = cms.untracked.double(5.),
0087             runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False))
0088    )
0089
0090 process.clustsummmultiplicitycorr = process.multiplicitycorr.clone()
0091 process.clustsummmultiplicitycorr.correlationConfigurations = cms.VPSet(
0092    cms.PSet(xMultiplicityMap = cms.InputTag("clustsummmultprod"),
0093             xDetSelection = cms.uint32(0), xDetLabel = cms.string("TK"), xBins = cms.uint32(3000), xMax=cms.double(100000),
0094             yMultiplicityMap = cms.InputTag("clustsummmultprod"),
0095             yDetSelection = cms.uint32(1005), yDetLabel = cms.string("Pixel"), yBins = cms.uint32(1000), yMax=cms.double(30000),
0096             rBins = cms.uint32(200), scaleFactor = cms.untracked.double(5.),
0097             runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False))
0098    )
0099
0100 process.multiplicitycorrtest1 = process.multiplicitycorr.clone()
0101 process.multiplicitycorrtest2 = process.multiplicitycorr.clone()
0102 process.multiplicitycorrtest1not = process.multiplicitycorr.clone()
0103 process.multiplicitycorrtest2not = process.multiplicitycorr.clone()
0104 process.clustsummmultiplicitycorrtest1 = process.multiplicitycorr.clone()
0105 process.clustsummmultiplicitycorrtest2 = process.multiplicitycorr.clone()
0106 process.clustsummmultiplicitycorrtest1not = process.multiplicitycorr.clone()
0107 process.clustsummmultiplicitycorrtest2not = process.multiplicitycorr.clone()
0108 process.multiplicitycorrstripconsistencytest1 = process.multiplicitycorr.clone()
0109 process.multiplicitycorrstripconsistencytest2 = process.multiplicitycorr.clone()
0110 process.multiplicitycorrpixelconsistencytest1 = process.multiplicitycorr.clone()
0111 process.multiplicitycorrpixelconsistencytest2 = process.multiplicitycorr.clone()
0112 process.multiplicitycorrstripconsistencytestnew1 = process.multiplicitycorr.clone()
0113 process.multiplicitycorrstripconsistencytestnew2 = process.multiplicitycorr.clone()
0114 process.multiplicitycorrpixelconsistencytestnew1 = process.multiplicitycorr.clone()
0115 process.multiplicitycorrpixelconsistencytestnew2 = process.multiplicitycorr.clone()
0116 process.multiplicitycorrpixelstripconsistencytest1 = process.multiplicitycorr.clone()
0117 process.multiplicitycorrpixelstripconsistencytestnot1 = process.multiplicitycorr.clone()
0118 process.multiplicitycorrpixelstripconsistencytest2 = process.multiplicitycorr.clone()
0119 process.multiplicitycorrpixelstripconsistencytestnot2 = process.multiplicitycorr.clone()
0120
0121 process.seqClusMultInvest = cms.Sequence(process.multiplicitycorr + process.clustsummmultiplicitycorr)
0122 #--------------------------------------------------------------------
0123
0124 process.load("DPGAnalysis.SiStripTools.largesipixelclusterevents_cfi")
0125 process.largeSiPixelClusterEvents.absoluteThreshold = 3000
0126 process.largeSiPixelClusterEvents.moduleThreshold = -1
0127
0128 process.load("DPGAnalysis.SiStripTools.largesistripclusterevents_cfi")
0129 process.largeSiStripClusterEvents.absoluteThreshold = 20000
0130 process.largeSiStripClusterEvents.moduleThreshold = -1
0131
0132 process.load("DPGAnalysis.SiStripTools.bysipixelclustmulteventfilter_cfi")
0133 process.bysipixelclustmulteventfilter.multiplicityConfig.moduleThreshold = -1
0134 process.bysipixelclustmulteventfilter.cut = cms.string("mult > 3000")
0135
0136 process.load("DPGAnalysis.SiStripTools.bysistripclustmulteventfilter_cfi")
0137 process.bysistripclustmulteventfilter.multiplicityConfig.moduleThreshold = -1
0138 process.bysistripclustmulteventfilter.cut = cms.string("mult > 20000")
0139
0140 process.load("DPGAnalysis.SiStripTools.byclustsummsipixelmulteventfilter_cfi")
0141 process.byclustsummsipixelmulteventfilter.cut = cms.string("mult > 3000")
0142
0143 process.load("DPGAnalysis.SiStripTools.byclustsummsistripmulteventfilter_cfi")
0144 process.byclustsummsistripmulteventfilter.cut = cms.string("mult > 20000")
0145
0146 process.stripfiltertest1 = cms.Sequence(process.largeSiStripClusterEvents + ~process.bysistripclustmulteventfilter)
0147 process.stripfiltertest2 = cms.Sequence(~process.largeSiStripClusterEvents + process.bysistripclustmulteventfilter)
0148
0149 process.pixelfiltertest1 = cms.Sequence(process.largeSiPixelClusterEvents + ~process.bysipixelclustmulteventfilter)
0150 process.pixelfiltertest2 = cms.Sequence(~process.largeSiPixelClusterEvents + process.bysipixelclustmulteventfilter)
0151
0152 process.stripfiltertestnew1 = cms.Sequence(process.byclustsummsistripmulteventfilter + ~process.bysistripclustmulteventfilter)
0153 process.stripfiltertestnew2 = cms.Sequence(~process.byclustsummsistripmulteventfilter + process.bysistripclustmulteventfilter)
0154
0155 process.pixelfiltertestnew1 = cms.Sequence(process.byclustsummsipixelmulteventfilter + ~process.bysipixelclustmulteventfilter)
0156 process.pixelfiltertestnew2 = cms.Sequence(~process.byclustsummsipixelmulteventfilter + process.bysipixelclustmulteventfilter)
0157
0158 process.load("DPGAnalysis.SiStripTools.bysipixelvssistripclustmulteventfilter_cfi")
0159 process.pixelvsstripfilter1 = process.bysipixelvssistripclustmulteventfilter.clone(cut=cms.string("(mult2 > 10000) && ( mult2 > 2000+7*mult1)"))
0160 process.pixelvsstripfilter2 = process.bysipixelvssistripclustmulteventfilter.clone(cut=cms.string("(mult1 > 1000) && (mult2 <30000) && ( mult2 < -2000+7*mult1)"))
0161
0162 process.load("DPGAnalysis.SiStripTools.byclustsummsipixelvssistripmulteventfilter_cfi")
0163 process.clustsummpixelvsstripfilter1 = process.byclustsummsipixelvssistripmulteventfilter.clone(cut=cms.string("(mult2 > 10000) && ( mult2 > 2000+7*mult1)"))
0164 process.clustsummpixelvsstripfilter2 = process.byclustsummsipixelvssistripmulteventfilter.clone(cut=cms.string("(mult1 > 1000) && (mult2 <30000) && ( mult2 < -2000+7*mult1)"))
0165
0166 process.pixelvsstripfiltertest1 = cms.Sequence(process.clustsummpixelvsstripfilter1 + ~process.pixelvsstripfilter1)
0167 process.pixelvsstripfiltertestnot1 = cms.Sequence(~process.clustsummpixelvsstripfilter1 + process.pixelvsstripfilter1)
0168 process.pixelvsstripfiltertest2 = cms.Sequence(process.clustsummpixelvsstripfilter2 + ~process.pixelvsstripfilter2)
0169 process.pixelvsstripfiltertestnot2 = cms.Sequence(~process.clustsummpixelvsstripfilter2 + process.pixelvsstripfilter2)
0170
0171
0172
0173 #-------------------------------------------------------------------------------------------
0174
0175 process.seqProducers = cms.Sequence(process.seqMultProd)
0176
0177 process.pstripfiltertest1 = cms.Path(process.seqProducers + process.stripfiltertest1 + process.multiplicitycorrstripconsistencytest1)
0178 process.pstripfiltertest2 = cms.Path(process.seqProducers + process.stripfiltertest2 + process.multiplicitycorrstripconsistencytest2)
0179 process.ppixelfiltertest1 = cms.Path(process.seqProducers + process.pixelfiltertest1 + process.multiplicitycorrpixelconsistencytest1)
0180 process.ppixelfiltertest2 = cms.Path(process.seqProducers + process.pixelfiltertest2 + process.multiplicitycorrpixelconsistencytest2)
0181
0182 process.pstripfiltertestnew1 = cms.Path(process.seqProducers + process.stripfiltertestnew1 + process.multiplicitycorrstripconsistencytestnew1)
0183 process.pstripfiltertestnew2 = cms.Path(process.seqProducers + process.stripfiltertestnew2 + process.multiplicitycorrstripconsistencytestnew2)
0184 process.ppixelfiltertestnew1 = cms.Path(process.seqProducers + process.pixelfiltertestnew1 + process.multiplicitycorrpixelconsistencytestnew1)
0185 process.ppixelfiltertestnew2 = cms.Path(process.seqProducers + process.pixelfiltertestnew2 + process.multiplicitycorrpixelconsistencytestnew2)
0186
0187 process.ppixelstripfiltertest1 = cms.Path(process.seqProducers + process.pixelvsstripfiltertest1 + process.multiplicitycorrpixelstripconsistencytest1)
0188 process.ppixelstripfiltertestnot1 = cms.Path(process.seqProducers + process.pixelvsstripfiltertestnot1 + process.multiplicitycorrpixelstripconsistencytestnot1)
0189 process.ppixelstripfiltertest2 = cms.Path(process.seqProducers + process.pixelvsstripfiltertest2 + process.multiplicitycorrpixelstripconsistencytest2)
0190 process.ppixelstripfiltertestnot2 = cms.Path(process.seqProducers + process.pixelvsstripfiltertestnot2 + process.multiplicitycorrpixelstripconsistencytestnot2)
0191
0192 process.p0 = cms.Path(
0193    process.seqProducers +
0194    process.seqClusMultInvest
0195    )
0196 process.pfiltertest1 = cms.Path(process.seqProducers + process.pixelvsstripfilter1 + process.multiplicitycorrtest1)
0197 process.pfiltertest2 = cms.Path(process.seqProducers + process.pixelvsstripfilter2 + process.multiplicitycorrtest2)
0198 process.pfiltertest1not = cms.Path(process.seqProducers + ~process.pixelvsstripfilter1 + process.multiplicitycorrtest1not)
0199 process.pfiltertest2not = cms.Path(process.seqProducers + ~process.pixelvsstripfilter2 + process.multiplicitycorrtest2not)
0200
0201 process.pclustsummfiltertest1 = cms.Path(process.seqProducers + process.clustsummpixelvsstripfilter1 + process.clustsummmultiplicitycorrtest1)
0202 process.pclustsummfiltertest2 = cms.Path(process.seqProducers + process.clustsummpixelvsstripfilter2 + process.clustsummmultiplicitycorrtest2)
0203 process.pclustsummfiltertest1not = cms.Path(process.seqProducers + ~process.clustsummpixelvsstripfilter1 + process.clustsummmultiplicitycorrtest1not)
0204 process.pclustsummfiltertest2not = cms.Path(process.seqProducers + ~process.clustsummpixelvsstripfilter2 + process.clustsummmultiplicitycorrtest2not)
0205
0206 #----GlobalTag ------------------------
0207
0208 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0209 from Configuration.AlCa.GlobalTag import GlobalTag
0210 process.GlobalTag = GlobalTag(process.GlobalTag, options.globalTag, '')
0211
0212
0213 process.TFileService = cms.Service('TFileService',
0214                                    fileName = cms.string('ByMultiplicityFilterTest.root')
0215                                    )
0216
0217 #print process.dumpPython()