Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 13:09:12

0001 import FWCore.ParameterSet.Config as cms
0002 import FWCore.ParameterSet.VarParsing as VarParsing
0003 
0004 process = cms.Process("MultiplicityProducerTest")
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 process.MessageLogger.debugModules = cms.untracked.vstring("clustsummmultprod","spclustermultprod","ssclustermultprod")
0056 
0057 #----Remove too verbose PrimaryVertexProducer
0058 
0059 process.MessageLogger.suppressInfo.append("pixelVerticesAdaptive")
0060 process.MessageLogger.suppressInfo.append("pixelVerticesAdaptiveNoBS")
0061 
0062 #----Remove too verbose BeamSpotOnlineProducer
0063 
0064 process.MessageLogger.suppressInfo.append("testBeamSpot")
0065 process.MessageLogger.suppressInfo.append("onlineBeamSpot")
0066 process.MessageLogger.suppressWarning.append("testBeamSpot")
0067 process.MessageLogger.suppressWarning.append("onlineBeamSpot")
0068 
0069 #----Remove too verbose TrackRefitter
0070 
0071 process.MessageLogger.suppressInfo.append("newTracksFromV0")
0072 process.MessageLogger.suppressInfo.append("newTracksFromOtobV0")
0073 
0074 #------------------------------------------------------------------
0075 
0076 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(options.maxEvents) )
0077 
0078 process.source = cms.Source("PoolSource",
0079                     fileNames = cms.untracked.vstring(options.inputFiles),
0080 #                    skipBadFiles = cms.untracked.bool(True),
0081                     inputCommands = cms.untracked.vstring("keep *", "drop *_MEtoEDMConverter_*_*")
0082                     )
0083 
0084 #--------------------------------------
0085 process.seqRECO = cms.Sequence()
0086 
0087 if options.fromRAW == 1:
0088     process.load("Configuration.StandardSequences.RawToDigi_Data_cff")
0089     process.load("Configuration.StandardSequences.MagneticField_cff")
0090     process.load("Configuration.StandardSequences.GeometryDB_cff")
0091     process.load("Configuration.StandardSequences.Reconstruction_cff")
0092     process.load("Configuration.StandardSequences.L1Reco_cff")
0093     process.seqRECO = cms.Sequence(process.gtEvmDigis + process.L1Reco
0094                                    + process.siStripDigis + process.siStripZeroSuppression + process.siStripClusters
0095                                    + process.siPixelDigis + process.siPixelClusters )
0096 
0097 
0098 #
0099 
0100 process.load("DPGAnalysis.SiStripTools.sipixelclustermultiplicityprod_cfi")
0101 process.load("DPGAnalysis.SiStripTools.sistripclustermultiplicityprod_cfi")
0102 
0103 #process.spclustermultprod.withClusterSize=cms.untracked.bool(True)
0104 #process.ssclustermultprod.withClusterSize=cms.untracked.bool(True)
0105 
0106 process.spclustermultprodnew = process.spclustermultprod.clone(wantedSubDets = cms.VPSet(    
0107     cms.PSet(detSelection = cms.uint32(101),detLabel = cms.string("BPIX"),selection=cms.untracked.vstring("0x1e000000-0x12000000")),
0108     cms.PSet(detSelection = cms.uint32(102),detLabel = cms.string("FPIX"),selection=cms.untracked.vstring("0x1e000000-0x14000000")),
0109     cms.PSet(detSelection = cms.uint32(11),detLabel = cms.string("BPIXL1"),selection=cms.untracked.vstring("0x1e0f0000-0x12010000")),
0110     cms.PSet(detSelection = cms.uint32(12),detLabel = cms.string("BPIXL2"),selection=cms.untracked.vstring("0x1e0f0000-0x12020000")),
0111     cms.PSet(detSelection = cms.uint32(13),detLabel = cms.string("BPIXL3"),selection=cms.untracked.vstring("0x1e0f0000-0x12030000")),
0112     cms.PSet(detSelection = cms.uint32(21),detLabel = cms.string("FPIXm"),selection=cms.untracked.vstring("0x1f800000-0x14800000")),
0113     cms.PSet(detSelection = cms.uint32(22),detLabel = cms.string("FPIXp"),selection=cms.untracked.vstring("0x1f800000-0x15000000"))
0114     )
0115 )
0116 process.ssclustermultprodnew = process.ssclustermultprod.clone(wantedSubDets = cms.VPSet(    
0117     cms.PSet(detSelection = cms.uint32(103),detLabel = cms.string("TIB"),selection=cms.untracked.vstring("0x1e000000-0x16000000")),
0118     cms.PSet(detSelection = cms.uint32(114),detLabel = cms.string("TIDm"),selection=cms.untracked.vstring("0x1e006000-0x18002000")),
0119     cms.PSet(detSelection = cms.uint32(124),detLabel = cms.string("TIDp"),selection=cms.untracked.vstring("0x1e006000-0x18004000")),
0120     cms.PSet(detSelection = cms.uint32(116),detLabel = cms.string("TECm"),selection=cms.untracked.vstring("0x1e0c0000-0x1c040000")),
0121     cms.PSet(detSelection = cms.uint32(126),detLabel = cms.string("TECp"),selection=cms.untracked.vstring("0x1e0c0000-0x1c080000")),
0122     cms.PSet(detSelection = cms.uint32(999),detLabel = cms.string("Module"),selection=cms.untracked.vstring("0x1fffffff-0x1600546d"))
0123     )
0124 )
0125 
0126 process.load("DPGAnalysis.SiStripTools.clustersummarymultiplicityprod_cfi")
0127 
0128 process.clustsummmultprod.wantedSubDets = cms.VPSet(
0129     cms.PSet(detSelection = cms.uint32(0),detLabel = cms.string("TK"), subDetEnum = cms.int32(0)        ),
0130     cms.PSet(detSelection = cms.uint32(1),detLabel = cms.string("TIB"), subDetEnum = cms.int32(1)       ),
0131     cms.PSet(detSelection = cms.uint32(2),detLabel = cms.string("TOB"), subDetEnum = cms.int32(2)       ),
0132     cms.PSet(detSelection = cms.uint32(3),detLabel = cms.string("TID"), subDetEnum = cms.int32(3)       ),
0133     cms.PSet(detSelection = cms.uint32(4),detLabel = cms.string("TEC"), subDetEnum = cms.int32(4)       ),
0134     cms.PSet(detSelection = cms.uint32(1005),detLabel = cms.string("Pixel"), subDetEnum = cms.int32(5)  ),
0135     cms.PSet(detSelection = cms.uint32(1006),detLabel = cms.string("FPIX"), subDetEnum = cms.int32(6)   ),
0136     cms.PSet(detSelection = cms.uint32(1007),detLabel = cms.string("BPIX"), subDetEnum = cms.int32(7)   ),
0137     )
0138 
0139 #process.ssclustermultprod.clusterdigiCollection = cms.InputTag("calZeroBiasClusters")
0140 #process.ssclustermultprodnew.clusterdigiCollection = cms.InputTag("calZeroBiasClusters")
0141 
0142 
0143 #process.seqMultProd = cms.Sequence(process.spclustermultprod+process.ssclustermultprod
0144 #                                   +process.spclustermultprodnew+process.ssclustermultprodnew
0145 #                                   +process.clustsummmultprod)
0146 process.seqMultProd = cms.Sequence(process.ssclustermultprod+process.ssclustermultprodnew)
0147 #process.seqMultProd = cms.Sequence(process.clustsummmultprod)
0148 
0149 process.goodVertices = cms.EDFilter("VertexSelector",
0150    src = cms.InputTag("offlinePrimaryVertices"),
0151    cut = cms.string("!isFake && ndof > 4 && abs(z) <= 24 && position.Rho <= 2"),  
0152    filter = cms.bool(False),   # otherwise it won't filter the events, just produce an empty vertex collection.
0153 )
0154 
0155 process.load("RecoLocalTracker.SiStripClusterizer.SiStripClusterToDigiProducer_cfi")
0156 
0157 
0158 process.seqProducers = cms.Sequence(process.seqRECO + process.seqMultProd + process.siStripClustersToDigis
0159                                     #+process.goodVertices
0160                                     )
0161 
0162 from HLTrigger.HLTfilters.triggerResultsFilter_cfi import *
0163 process.hltSelection = triggerResultsFilter.clone(
0164                                           triggerConditions = cms.vstring("HLT_ZeroBias_*"),
0165                                           hltResults = cms.InputTag( "TriggerResults", "", "HLT" ),
0166                                           l1tResults = cms.InputTag( "" ),
0167                                           throw = cms.bool(False)
0168                                           )
0169 
0170 process.manystripclus = cms.EDFilter('BySiPixelClusterVsSiStripClusterMultiplicityEventFilter',
0171                                      multiplicityConfig = cms.PSet(
0172     firstMultiplicityConfig = cms.PSet(
0173     collectionName = cms.InputTag("siPixelClusters"),
0174     moduleThreshold = cms.untracked.int32(-1),
0175     useQuality = cms.untracked.bool(False),
0176     qualityLabel = cms.untracked.string("")
0177     ),
0178     secondMultiplicityConfig = cms.PSet(
0179     collectionName = cms.InputTag("siStripClusters"),
0180     moduleThreshold = cms.untracked.int32(-1),
0181     useQuality = cms.untracked.bool(False),
0182     qualityLabel = cms.untracked.string("")
0183     )
0184     ),
0185                                      cut = cms.string("( mult2 > 20000+7*mult1)")
0186                                      )
0187 
0188 process.manystripclus53X = cms.EDFilter('ByClusterSummaryMultiplicityPairEventFilter',
0189                                         multiplicityConfig = cms.PSet(
0190     firstMultiplicityConfig = cms.PSet(
0191     clusterSummaryCollection = cms.InputTag("clusterSummaryProducer"),
0192     subDetEnum = cms.int32(5),
0193     subDetVariable = cms.string("pHits")
0194     ),
0195     secondMultiplicityConfig = cms.PSet(
0196     clusterSummaryCollection = cms.InputTag("clusterSummaryProducer"),
0197     subDetEnum = cms.int32(0),
0198     subDetVariable = cms.string("cHits")
0199     ),
0200     ),
0201                                         cut = cms.string("( mult2 > 20000+7*mult1)")
0202                                         )
0203 
0204 process.load("DPGAnalysis.SiStripTools.ssclusmultinvestigator_cfi")
0205 #process.ssclusmultinvestigator.vertexCollection = cms.InputTag("goodVertices")
0206 #process.ssclusmultinvestigator.wantInvestHist = cms.bool(True)
0207 #process.ssclusmultinvestigator.wantVtxCorrHist = cms.bool(True)
0208 #process.ssclusmultinvestigator.digiVtxCorrConfig = cms.PSet(
0209 #    wantedSubDets = process.ssclusmultinvestigator.wantedSubDets,
0210 #    hitName = cms.untracked.string("cluster"),
0211 #    numberOfBins = cms.untracked.int32(100),   
0212 #    scaleFactor = cms.untracked.int32(1),
0213 #    maxNvtx = cms.untracked.int32(100)
0214 #    )
0215 #process.ssclusmultinvestigator.wantLumiCorrHist = cms.bool(True)
0216 #process.ssclusmultinvestigator.digiLumiCorrConfig = cms.PSet(
0217 #    lumiProducer = cms.InputTag("lumiProducer"),
0218 #    wantedSubDets = process.ssclusmultinvestigator.wantedSubDets,
0219 #    hitName = cms.untracked.string("cluster"),
0220 #    numberOfBins = cms.untracked.int32(100),   
0221 #    scaleFactor = cms.untracked.int32(1),
0222 #    maxLumi = cms.untracked.double(25)
0223 #    )
0224 
0225 
0226 #process.ssclusmultinvestigator.runHisto = cms.untracked.bool(True)
0227 process.ssclusmultinvestigator.scaleFactor=cms.untracked.int32(1)
0228 
0229 process.ssclusmultinvestigatornew = process.ssclusmultinvestigator.clone()
0230 process.ssclusmultinvestigatornew.runHisto = cms.untracked.bool(True)
0231 process.ssclusmultinvestigatornew.wantedSubDets = cms.untracked.VPSet(    
0232     cms.PSet(detSelection = cms.uint32(103),detLabel = cms.string("TIB"), binMax = cms.int32(1787904/32)),
0233     cms.PSet(detSelection = cms.uint32(114),detLabel = cms.string("TIDm"), binMax = cms.int32(565248/32)),
0234     cms.PSet(detSelection = cms.uint32(124),detLabel = cms.string("TIDp"), binMax = cms.int32(565248/32)),
0235     cms.PSet(detSelection = cms.uint32(116),detLabel = cms.string("TECm"), binMax = cms.int32(3866624/64)),
0236     cms.PSet(detSelection = cms.uint32(126),detLabel = cms.string("TECp"), binMax = cms.int32(3866624/64)),
0237     cms.PSet(detSelection = cms.uint32(999),detLabel = cms.string("Module"), binMax = cms.int32(200))
0238     )
0239 #process.ssclusmultinvestigatornew.digiVtxCorrConfig.wantedSubDets = cms.untracked.VPSet(    
0240 #    cms.PSet(detSelection = cms.uint32(103),detLabel = cms.string("TIB"), binMax = cms.int32(1787904/32)),
0241 #    cms.PSet(detSelection = cms.uint32(114),detLabel = cms.string("TIDm"), binMax = cms.int32(565248/32)),
0242 #    cms.PSet(detSelection = cms.uint32(124),detLabel = cms.string("TIDp"), binMax = cms.int32(565248/32)),
0243 #    cms.PSet(detSelection = cms.uint32(116),detLabel = cms.string("TECm"), binMax = cms.int32(3866624/64)),
0244 #    cms.PSet(detSelection = cms.uint32(126),detLabel = cms.string("TECp"), binMax = cms.int32(3866624/64)),
0245 #    cms.PSet(detSelection = cms.uint32(999),detLabel = cms.string("Module"), binMax = cms.int32(200))
0246 #    )
0247 #process.ssclusmultinvestigatornew.digiLumiCorrConfig.wantedSubDets = cms.untracked.VPSet(    
0248 #    cms.PSet(detSelection = cms.uint32(103),detLabel = cms.string("TIB"), binMax = cms.int32(1787904/32)),
0249 #    cms.PSet(detSelection = cms.uint32(114),detLabel = cms.string("TIDm"), binMax = cms.int32(565248/32)),
0250 #    cms.PSet(detSelection = cms.uint32(124),detLabel = cms.string("TIDp"), binMax = cms.int32(565248/32)),
0251 #    cms.PSet(detSelection = cms.uint32(116),detLabel = cms.string("TECm"), binMax = cms.int32(3866624/64)),
0252 #    cms.PSet(detSelection = cms.uint32(126),detLabel = cms.string("TECp"), binMax = cms.int32(3866624/64)),
0253 #    cms.PSet(detSelection = cms.uint32(999),detLabel = cms.string("Module"), binMax = cms.int32(200))
0254 #    )
0255 process.ssclusmultinvestigatornew.multiplicityMap = cms.InputTag("ssclustermultprodnew")
0256 
0257 
0258 process.load("DPGAnalysis.SiStripTools.spclusmultinvestigator_cfi")
0259 process.spclusmultinvestigator.vertexCollection = cms.InputTag("goodVertices")
0260 process.spclusmultinvestigator.wantInvestHist = cms.bool(True)
0261 process.spclusmultinvestigator.wantVtxCorrHist = cms.bool(True)
0262 process.spclusmultinvestigator.digiVtxCorrConfig = cms.PSet(
0263     wantedSubDets = process.spclusmultinvestigator.wantedSubDets,
0264     hitName = cms.untracked.string("cluster"),
0265     numberOfBins = cms.untracked.int32(100),   
0266     scaleFactor = cms.untracked.int32(5)
0267     )
0268 process.spclusmultinvestigator.wantLumiCorrHist = cms.bool(True)
0269 process.spclusmultinvestigator.digiLumiCorrConfig = cms.PSet(
0270     lumiProducer = cms.InputTag("lumiProducer"),
0271     wantedSubDets = process.spclusmultinvestigator.wantedSubDets,
0272     hitName = cms.untracked.string("cluster"),
0273     numberOfBins = cms.untracked.int32(100),   
0274     scaleFactor = cms.untracked.int32(5)
0275     )
0276 #process.spclusmultinvestigator.runHisto = cms.untracked.bool(True)
0277 process.spclusmultinvestigator.scaleFactor=cms.untracked.int32(5)
0278 
0279 process.spclusmultinvestigatornew = process.spclusmultinvestigator.clone()
0280 process.spclusmultinvestigatornew.wantedSubDets = cms.untracked.VPSet(    
0281     cms.PSet(detSelection = cms.uint32(101),detLabel = cms.string("BPIX"), binMax = cms.int32(100000)),
0282     cms.PSet(detSelection = cms.uint32(11),detLabel = cms.string("BPIXL1"), binMax = cms.int32(100000)),
0283     cms.PSet(detSelection = cms.uint32(12),detLabel = cms.string("BPIXL2"), binMax = cms.int32(100000)),
0284     cms.PSet(detSelection = cms.uint32(13),detLabel = cms.string("BPIXL3"), binMax = cms.int32(100000)),
0285     cms.PSet(detSelection = cms.uint32(102),detLabel = cms.string("FPIX"), binMax = cms.int32(100000)),
0286     cms.PSet(detSelection = cms.uint32(21),detLabel = cms.string("FPIXm"), binMax = cms.int32(100000)),
0287     cms.PSet(detSelection = cms.uint32(22),detLabel = cms.string("FPIXp"), binMax = cms.int32(100000))
0288     )
0289 process.spclusmultinvestigatornew.digiVtxCorrConfig.wantedSubDets = cms.untracked.VPSet(    
0290     cms.PSet(detSelection = cms.uint32(101),detLabel = cms.string("BPIX"), binMax = cms.int32(100000)),
0291     cms.PSet(detSelection = cms.uint32(11),detLabel = cms.string("BPIXL1"), binMax = cms.int32(100000)),
0292     cms.PSet(detSelection = cms.uint32(12),detLabel = cms.string("BPIXL2"), binMax = cms.int32(100000)),
0293     cms.PSet(detSelection = cms.uint32(13),detLabel = cms.string("BPIXL3"), binMax = cms.int32(100000)),
0294     cms.PSet(detSelection = cms.uint32(102),detLabel = cms.string("FPIX"), binMax = cms.int32(100000)),
0295     cms.PSet(detSelection = cms.uint32(21),detLabel = cms.string("FPIXm"), binMax = cms.int32(100000)),
0296     cms.PSet(detSelection = cms.uint32(22),detLabel = cms.string("FPIXp"), binMax = cms.int32(100000))
0297     )
0298 process.spclusmultinvestigatornew.digiLumiCorrConfig.wantedSubDets = cms.untracked.VPSet(    
0299     cms.PSet(detSelection = cms.uint32(101),detLabel = cms.string("BPIX"), binMax = cms.int32(100000)),
0300     cms.PSet(detSelection = cms.uint32(11),detLabel = cms.string("BPIXL1"), binMax = cms.int32(100000)),
0301     cms.PSet(detSelection = cms.uint32(12),detLabel = cms.string("BPIXL2"), binMax = cms.int32(100000)),
0302     cms.PSet(detSelection = cms.uint32(13),detLabel = cms.string("BPIXL3"), binMax = cms.int32(100000)),
0303     cms.PSet(detSelection = cms.uint32(102),detLabel = cms.string("FPIX"), binMax = cms.int32(100000)),
0304     cms.PSet(detSelection = cms.uint32(21),detLabel = cms.string("FPIXm"), binMax = cms.int32(100000)),
0305     cms.PSet(detSelection = cms.uint32(22),detLabel = cms.string("FPIXp"), binMax = cms.int32(100000))
0306     )
0307 process.spclusmultinvestigatornew.multiplicityMap = cms.InputTag("spclustermultprodnew")
0308 
0309 process.load("DPGAnalysis.SiStripTools.multiplicitycorr_cfi")
0310 process.multiplicitycorr.correlationConfigurations = cms.VPSet(
0311    cms.PSet(xMultiplicityMap = cms.InputTag("spclustermultprodnew"),
0312             xDetSelection = cms.uint32(12), xDetLabel = cms.string("BPixL2"), xBins = cms.uint32(1000), xMax=cms.double(10000), 
0313             yMultiplicityMap = cms.InputTag("ssclustermultprodnew"),
0314             yDetSelection = cms.uint32(103), yDetLabel = cms.string("TIB"), yBins = cms.uint32(1000), yMax=cms.double(50000),
0315             rBins = cms.uint32(200), scaleFactor = cms.untracked.double(.25),
0316             runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False)),
0317    cms.PSet(xMultiplicityMap = cms.InputTag("spclustermultprodnew"),
0318             xDetSelection = cms.uint32(12), xDetLabel = cms.string("BPixL2"), xBins = cms.uint32(1000), xMax=cms.double(10000), 
0319             yMultiplicityMap = cms.InputTag("ssclustermultprodnew"),
0320             yDetSelection = cms.uint32(114), yDetLabel = cms.string("TIDm"), yBins = cms.uint32(1000), yMax=cms.double(20000),
0321             rBins = cms.uint32(200), scaleFactor = cms.untracked.double(.5),
0322             runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False)),
0323    cms.PSet(xMultiplicityMap = cms.InputTag("spclustermultprodnew"),
0324             xDetSelection = cms.uint32(12), xDetLabel = cms.string("BPixL2"), xBins = cms.uint32(1000), xMax=cms.double(10000), 
0325             yMultiplicityMap = cms.InputTag("ssclustermultprodnew"),
0326             yDetSelection = cms.uint32(124), yDetLabel = cms.string("TIDp"), yBins = cms.uint32(1000), yMax=cms.double(20000),
0327             rBins = cms.uint32(200), scaleFactor = cms.untracked.double(.5),
0328             runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False)),
0329    cms.PSet(xMultiplicityMap = cms.InputTag("spclustermultprodnew"),
0330             xDetSelection = cms.uint32(12), xDetLabel = cms.string("BPixL2"), xBins = cms.uint32(1000), xMax=cms.double(10000), 
0331             yMultiplicityMap = cms.InputTag("spclustermultprodnew"),
0332             yDetSelection = cms.uint32(11), yDetLabel = cms.string("BPixL1"), yBins = cms.uint32(1000), yMax=cms.double(10000),
0333             rBins = cms.uint32(200), scaleFactor = cms.untracked.double(1.),
0334             runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False))
0335    )
0336 
0337 process.multcorrclustsumm = process.multiplicitycorr.clone()
0338 process.multcorrclustsumm.correlationConfigurations = cms.VPSet(
0339    cms.PSet(yMultiplicityMap = cms.InputTag("spclustermultprod"),
0340             yDetSelection = cms.uint32(0), yDetLabel = cms.string("Pixel"), yBins = cms.uint32(1000), yMax=cms.double(1000), 
0341             xMultiplicityMap = cms.InputTag("clustsummmultprod"),
0342             xDetSelection = cms.uint32(1005), xDetLabel = cms.string("Pixelcs"), xBins = cms.uint32(1000), xMax=cms.double(1000),
0343             rBins = cms.uint32(200), scaleFactor = cms.untracked.double(1.),
0344             runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False)),
0345    cms.PSet(xMultiplicityMap = cms.InputTag("spclustermultprodnew"),
0346             xDetSelection = cms.uint32(101), xDetLabel = cms.string("BPIX"), xBins = cms.uint32(1000), xMax=cms.double(1000), 
0347             yMultiplicityMap = cms.InputTag("clustsummmultprod"),
0348             yDetSelection = cms.uint32(1007), yDetLabel = cms.string("BPIXcs"), yBins = cms.uint32(1000), yMax=cms.double(1000),
0349             rBins = cms.uint32(200), scaleFactor = cms.untracked.double(1.),
0350             runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False)),
0351    cms.PSet(xMultiplicityMap = cms.InputTag("spclustermultprodnew"),
0352             xDetSelection = cms.uint32(102), xDetLabel = cms.string("FPIX"), xBins = cms.uint32(1000), xMax=cms.double(1000), 
0353             yMultiplicityMap = cms.InputTag("clustsummmultprod"),
0354             yDetSelection = cms.uint32(1006), yDetLabel = cms.string("FPIXcs"), yBins = cms.uint32(1000), yMax=cms.double(1000),
0355             rBins = cms.uint32(200), scaleFactor = cms.untracked.double(1.),
0356             runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False)),
0357    cms.PSet(xMultiplicityMap = cms.InputTag("ssclustermultprod"),
0358             xDetSelection = cms.uint32(3), xDetLabel = cms.string("TIB"), xBins = cms.uint32(1000), xMax=cms.double(50000), 
0359             yMultiplicityMap = cms.InputTag("clustsummmultprod"),
0360             yDetSelection = cms.uint32(1), yDetLabel = cms.string("TIBcs"), yBins = cms.uint32(1000), yMax=cms.double(50000),
0361             rBins = cms.uint32(200), scaleFactor = cms.untracked.double(1.),
0362             runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False)),
0363    cms.PSet(xMultiplicityMap = cms.InputTag("ssclustermultprod"),
0364             xDetSelection = cms.uint32(4), xDetLabel = cms.string("TID"), xBins = cms.uint32(1000), xMax=cms.double(50000), 
0365             yMultiplicityMap = cms.InputTag("clustsummmultprod"),
0366             yDetSelection = cms.uint32(3), yDetLabel = cms.string("TIDcs"), yBins = cms.uint32(1000), yMax=cms.double(50000),
0367             rBins = cms.uint32(200), scaleFactor = cms.untracked.double(1.),
0368             runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False)),
0369    cms.PSet(xMultiplicityMap = cms.InputTag("ssclustermultprod"),
0370             xDetSelection = cms.uint32(5), xDetLabel = cms.string("TOB"), xBins = cms.uint32(1000), xMax=cms.double(50000), 
0371             yMultiplicityMap = cms.InputTag("clustsummmultprod"),
0372             yDetSelection = cms.uint32(2), yDetLabel = cms.string("TOBcs"), yBins = cms.uint32(1000), yMax=cms.double(50000),
0373             rBins = cms.uint32(200), scaleFactor = cms.untracked.double(1.),
0374             runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False)),
0375    cms.PSet(xMultiplicityMap = cms.InputTag("ssclustermultprod"),
0376             xDetSelection = cms.uint32(6), xDetLabel = cms.string("TEC"), xBins = cms.uint32(1000), xMax=cms.double(50000), 
0377             yMultiplicityMap = cms.InputTag("clustsummmultprod"),
0378             yDetSelection = cms.uint32(4), yDetLabel = cms.string("TECcs"), yBins = cms.uint32(1000), yMax=cms.double(50000),
0379             rBins = cms.uint32(200), scaleFactor = cms.untracked.double(1.),
0380             runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False))
0381    )
0382 
0383 
0384 process.clustsummmultcorr =  cms.EDAnalyzer('MultiplicityCorrelator',
0385                                             correlationConfigurations = cms.VPSet(
0386     cms.PSet(xMultiplicityMap = cms.InputTag("clustsummmultprod"), xDetSelection = cms.uint32(0), xDetLabel = cms.string("TK"), xBins = cms.uint32(1000), xMax=cms.double(150000),
0387              yMultiplicityMap = cms.InputTag("clustsummmultprod"), yDetSelection = cms.uint32(1005), yDetLabel = cms.string("Pixel"), yBins = cms.uint32(1000), yMax=cms.double(20000),
0388              rBins = cms.uint32(200), scaleFactor = cms.untracked.double(5.),
0389              runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False))
0390     )
0391                                             )
0392 #process.clustsummmultcorr =  cms.EDAnalyzer('MultiplicityCorrelator',
0393 #                                            correlationConfigurations = cms.VPSet(
0394 #    cms.PSet(xMultiplicityMap = cms.InputTag("ssclustermultprod"), xDetSelection = cms.uint32(0), xDetLabel = cms.string("TK"), xBins = cms.uint32(1000), xMax=cms.double(150000),
0395 #             yMultiplicityMap = cms.InputTag("spclustermultprod"), yDetSelection = cms.uint32(0), yDetLabel = cms.string("Pixel"), yBins = cms.uint32(1000), yMax=cms.double(20000),
0396 #             rBins = cms.uint32(200), scaleFactor = cms.untracked.double(5.),
0397 #             runHisto=cms.bool(False),runHistoBXProfile=cms.bool(False),runHistoBX=cms.bool(False),runHisto2D=cms.bool(False))
0398 #    )
0399 #                                            )
0400 
0401 process.clustsummmultcorrmanystripclus = process.clustsummmultcorr.clone()
0402 
0403 process.load("TrackingPFG.Utilities.logerroranalyzer_cfi")
0404 
0405 process.load("DPGAnalysis.SiStripTools.clusterbigeventsdebugger_cfi")
0406 process.clusterbigeventsdebugger.selections.append(cms.PSet(label=cms.string("Module"),selection=cms.untracked.vstring("0x1fffffff-0x1600546d")))
0407 #process.clusterbigeventsdebugger.collection = cms.InputTag("calZeroBiasClusters")
0408 #process.clusterbigeventsdebugger.foldedStrips = cms.untracked.bool(True)
0409 #process.clusterbigeventsdebugger.singleEvents = cms.bool(True)
0410 
0411 process.load("DPGAnalysis.SiStripTools.digibigeventsdebugger_cfi")
0412 process.digibigeventsdebugger.selections = cms.VPSet(
0413     cms.PSet(label=cms.string("TIB"),selection=cms.untracked.vstring("0x1e000000-0x16000000")),
0414     cms.PSet(label=cms.string("TIBD"),selection=cms.untracked.vstring("0x1e000000-0x16000000","0x1e000000-0x18000000")),
0415     cms.PSet(label=cms.string("TECp"),selection=cms.untracked.vstring("0x1e0c0000-0x1c080000")),
0416     cms.PSet(label=cms.string("TECm"),selection=cms.untracked.vstring("0x1e0c0000-0x1c040000")),
0417     cms.PSet(label=cms.string("TOB"),selection=cms.untracked.vstring("0x1e000000-0x1a000000")),
0418     cms.PSet(label=cms.string("TID"),selection=cms.untracked.vstring("0x1e000000-0x18000000"))
0419     )
0420 
0421 process.digibigeventsdebugger.selections.append(cms.PSet(label=cms.string("Module"),selection=cms.untracked.vstring("0x1fffffff-0x1600546d")))
0422 process.digibigeventsdebugger.collection = cms.InputTag("siStripClustersToDigis","ZeroSuppressed")
0423 process.digibigeventsdebugger.foldedStrips = cms.untracked.bool(True)
0424 #process.digibigeventsdebugger.singleEvents = cms.bool(True)
0425 
0426 process.load("DPGAnalysis.SiStripTools.fedbadmodulefilter_cfi")
0427 process.fedbadmodulefilter.badModThr = 0
0428 process.fedbadmodulefilter.wantedHisto = cms.untracked.bool(True)
0429 process.fedbadmodulefilter.moduleSelection = cms.untracked.PSet(selection=cms.untracked.vstring("0x1fffffff-0x1600546d"))
0430 
0431 process.fedbadmoduleTECmD8 = process.fedbadmodulefilter.clone(moduleSelection = cms.untracked.PSet(selection=cms.untracked.vstring("0x1e0fc000-0x1c060000")))
0432 process.fedbadmoduleTECmD4 = process.fedbadmodulefilter.clone(moduleSelection = cms.untracked.PSet(selection=cms.untracked.vstring("0x1e0fc000-0x1c050000")))
0433 process.fedbadmoduleTECmD7 = process.fedbadmodulefilter.clone(moduleSelection = cms.untracked.PSet(selection=cms.untracked.vstring("0x1e0fc000-0x1c05c000")))
0434 
0435 
0436 #process.seqClusMultInvest = cms.Sequence(process.spclusmultinvestigator + process.ssclusmultinvestigator +
0437 #                                         process.spclusmultinvestigatornew + process.ssclusmultinvestigatornew
0438 #                                         + process.multiplicitycorr + process.multcorrclustsumm) 
0439 #process.seqClusMultInvest = cms.Sequence(process.clustsummmultcorr) 
0440 process.seqClusMultInvest = cms.Sequence(process.ssclusmultinvestigatornew +
0441                                          process.clusterbigeventsdebugger + process.digibigeventsdebugger +
0442                                          process.fedbadmodulefilter + process.fedbadmoduleTECmD8 + process.fedbadmoduleTECmD4 + process.fedbadmoduleTECmD7 
0443                                          ) 
0444 
0445 
0446 
0447 process.p0 = cms.Path(
0448     #    process.hltSelection +
0449     process.seqProducers +
0450     process.seqClusMultInvest)
0451 
0452 #process.pmanystripclus = cms.Path(process.seqProducers +
0453 #                                  process.manystripclus53X +
0454 ##                                  process.manystripclus +
0455 #                                  process.clustsummmultcorrmanystripclus +
0456 #                                  process.logerroranalyzer
0457 #                                  )
0458 
0459 #process.outoddevent = cms.OutputModule("PoolOutputModule",
0460 #                                       fileName = cms.untracked.string("manystripclusevents.root"),
0461 #                                       outputCommands = cms.untracked.vstring("keep *"),
0462 #                                       SelectEvents = cms.untracked.PSet(SelectEvents = cms.vstring( "pmanystripclus"))
0463 #                                       )                                       
0464 
0465 #process.e = cms.EndPath(process.outoddevent)
0466 
0467 #----GlobalTag ------------------------
0468 
0469 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0470 from Configuration.AlCa.GlobalTag import GlobalTag
0471 process.GlobalTag = GlobalTag(process.GlobalTag, options.globalTag, '')
0472 
0473 if options.testTag==1:
0474     process.stripConditions = cms.ESSource(
0475         "PoolDBESSource",
0476         process.CondDBSetup,
0477         connect = cms.string('frontier://FrontierProd/CMS_COND_31X_STRIP'),
0478         toGet = cms.VPSet(
0479         cms.PSet(
0480         record = cms.string('SiStripNoisesRcd'), tag = cms.string('SiStripNoise_GR10_v2_offline')
0481         ),
0482         ),
0483         )
0484 
0485     process.es_prefer_strips = cms.ESPrefer("PoolDBESSource","stripConditions")
0486 
0487 process.TFileService = cms.Service('TFileService',
0488                                    fileName = cms.string('MultiplicityProducerTest.root')
0489                                    )
0490 
0491 #print process.dumpPython()