Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-04-06 12:06:44

0001 import FWCore.ParameterSet.Config as cms
0002 import FWCore.ParameterSet.VarParsing as VarParsing
0003 
0004 
0005 from PhysicsTools.PatAlgos.tools.helpers import cloneProcessingSnippet
0006 
0007 process = cms.Process("OOTmultiplicity")
0008 
0009 #prepare options
0010 
0011 options = VarParsing.VarParsing("analysis")
0012 
0013 options.register ('globalTag',
0014                   "DONOTEXIST",
0015                   VarParsing.VarParsing.multiplicity.singleton, # singleton or list
0016                   VarParsing.VarParsing.varType.string,          # string, int, or float
0017                   "GlobalTag")
0018 
0019 options.parseArguments()
0020 
0021 process.options = cms.untracked.PSet(
0022     wantSummary = cms.untracked.bool(True),
0023     fileMode = cms.untracked.string("FULLMERGE")
0024     )
0025 
0026 process.load("FWCore.MessageService.MessageLogger_cfi")
0027 
0028 #----------------------------------------------------------------
0029 
0030 
0031 process.MessageLogger.cout.enable = cms.untracked.bool(True)
0032 process.MessageLogger.cout.threshold = cms.untracked.string("INFO")
0033 process.MessageLogger.cout.default = cms.untracked.PSet(
0034     limit = cms.untracked.int32(10000000)
0035     )
0036 process.MessageLogger.cout.CondDBESSource = cms.untracked.PSet(
0037     limit = cms.untracked.int32(0)
0038     )
0039 process.MessageLogger.cout.FwkReport = cms.untracked.PSet(
0040     reportEvery = cms.untracked.int32(10000)
0041     )
0042 
0043 process.MessageLogger.cerr.enable = cms.untracked.bool(True)
0044 process.MessageLogger.cerr.threshold = cms.untracked.string("WARNING")
0045 process.MessageLogger.cerr.default = cms.untracked.PSet(
0046     limit = cms.untracked.int32(10000000)
0047     )
0048 process.MessageLogger.cerr.FwkReport = cms.untracked.PSet(
0049     reportEvery = cms.untracked.int32(100000)
0050     )
0051 
0052 process.MessageLogger.suppressWarning.append("consecutiveHEs")
0053 
0054 
0055 #------------------------------------------------------------------
0056 
0057 #process.MessageLogger.cout.threshold = cms.untracked.string("DEBUG")
0058 #process.MessageLogger.debugModules = cms.untracked.vstring("eventtimedistribution")
0059 
0060 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(options.maxEvents) )
0061 
0062 process.source = cms.Source("PoolSource",
0063                     fileNames = cms.untracked.vstring(options.inputFiles),
0064 #                    skipBadFiles = cms.untracked.bool(True),
0065                     inputCommands = cms.untracked.vstring("keep *", "drop *_MEtoEDMConverter_*_*")
0066                     )
0067 
0068 
0069 
0070 process.load("DPGAnalysis.SiStripTools.eventwithhistoryproducerfroml1abc_cfi")
0071 
0072 import DPGAnalysis.SiStripTools.apvcyclephaseproducerfroml1ts2011_cfi 
0073 process.APVPhases = DPGAnalysis.SiStripTools.apvcyclephaseproducerfroml1ts2011_cfi.APVPhases 
0074 
0075 process.APVPhases.wantHistos = cms.untracked.bool(True)
0076 
0077 # Filter for no noisy bins
0078 
0079 process.load("DPGAnalysis.SiStripTools.filters.Potential_TIBTEC_HugeEvents_AlCaReco_cfi")
0080 process.load("DPGAnalysis.SiStripTools.filters.Potential_TIBTEC_FrameHeaderEvents_firstpeak_AlCaReco_cfi")
0081 process.load("DPGAnalysis.SiStripTools.filters.Potential_TIBTEC_FrameHeaderEvents_additionalpeak_AlCaReco_cfi")
0082 
0083 process.latencyPlusOne = cms.EDFilter('EventWithHistoryEDFilter',
0084                                       commonConfiguration = cms.untracked.PSet(historyProduct = cms.untracked.InputTag("consecutiveHEs"),
0085                                                                                APVPhaseLabel  = cms.untracked.string("APVPhases"),
0086                                                                                ),
0087                                       filterConfigurations = cms.untracked.VPSet(cms.PSet(apvModes = cms.untracked.vint32(47),
0088                                                                                           dbxRangeLtcyAware = cms.untracked.vint32(1,1)
0089                                                                                           )
0090                                                                                  )
0091                                       )
0092 
0093 process.PotentialNewHugeEvents = cms.EDFilter('EventWithHistoryEDFilter',
0094                                               commonConfiguration = cms.untracked.PSet(historyProduct = cms.untracked.InputTag("consecutiveHEs"),
0095                                                                                        APVPhaseLabel  = cms.untracked.string("APVPhases"),
0096                                                                                        ),
0097                                               filterConfigurations = cms.untracked.VPSet(cms.PSet(apvModes = cms.untracked.vint32(37),
0098                                                                                                   partitionName = cms.untracked.string("Any"),
0099                                                                                                   absBXInCycleRangeLtcyAware = cms.untracked.vint32(12,12)
0100                                                                                                   )
0101                                                                                          )
0102                                               )
0103 
0104 process.PotentialAdditionalHugeEvents = cms.EDFilter('EventWithHistoryEDFilter',
0105                                                      commonConfiguration = cms.untracked.PSet(historyProduct = cms.untracked.InputTag("consecutiveHEs"),
0106                                                                                               APVPhaseLabel  = cms.untracked.string("APVPhases"),
0107                                                                                               ),
0108                                                      filterConfigurations = cms.untracked.VPSet(cms.PSet(apvModes = cms.untracked.vint32(47),
0109                                                                                                   partitionName = cms.untracked.string("Any"),
0110                                                                                                   absBXInCycleRangeLtcyAware = cms.untracked.vint32(9,10)
0111                                                                                                   )
0112                                                                                          )
0113                                               )
0114 
0115 
0116 process.seqNoNoisyBins = cms.Sequence(~process.PotentialTIBTECHugeEvents
0117                                       + ~process.PotentialTIBTECFrameHeaderEventsFPeak
0118                                       + ~process.PotentialTIBTECFrameHeaderEventsAdditionalPeak
0119                                       + ~process.latencyPlusOne
0120                                       + ~process.PotentialNewHugeEvents
0121                                       + ~process.PotentialAdditionalHugeEvents
0122                                       )
0123 
0124 
0125 process.load("DPGAnalysis.SiStripTools.apvcyclephasemonitor_cfi")
0126 
0127 process.load("DPGAnalysis.SiStripTools.eventtimedistribution_cfi")
0128 process.eventtimedistribution.wantEWHDepthHisto = cms.untracked.bool(True)
0129 process.eventtimedistrRandom = process.eventtimedistribution.clone()
0130 process.eventtimedistrnonoisybins = process.eventtimedistribution.clone()
0131 process.eventtimedistrnonoisybinsRandom = process.eventtimedistribution.clone()
0132 
0133 
0134 # cluster analysis
0135 
0136 process.load("DPGAnalysis.SiStripTools.sistripclustermultiplicityprod_cfi")
0137 process.ssclustermultprod.wantedSubDets.extend(cms.VPSet(
0138     cms.PSet(detSelection = cms.uint32(101),detLabel = cms.string("TIBL1"),selection=cms.untracked.vstring("0x1e01c000-0x16004000")),
0139     cms.PSet(detSelection = cms.uint32(102),detLabel = cms.string("TIBL2"),selection=cms.untracked.vstring("0x1e01c000-0x16008000")),
0140     cms.PSet(detSelection = cms.uint32(103),detLabel = cms.string("TIBL3"),selection=cms.untracked.vstring("0x1e01c000-0x1600c000")),
0141     cms.PSet(detSelection = cms.uint32(104),detLabel = cms.string("TIBL4"),selection=cms.untracked.vstring("0x1e01c000-0x16010000")),
0142     cms.PSet(detSelection = cms.uint32(201),detLabel = cms.string("TOBL1"),selection=cms.untracked.vstring("0x1e01c000-0x1a004000")),
0143     cms.PSet(detSelection = cms.uint32(202),detLabel = cms.string("TOBL2"),selection=cms.untracked.vstring("0x1e01c000-0x1a008000")),
0144     cms.PSet(detSelection = cms.uint32(203),detLabel = cms.string("TOBL3"),selection=cms.untracked.vstring("0x1e01c000-0x1a00c000")),
0145     cms.PSet(detSelection = cms.uint32(204),detLabel = cms.string("TOBL4"),selection=cms.untracked.vstring("0x1e01c000-0x1a010000")),
0146     cms.PSet(detSelection = cms.uint32(205),detLabel = cms.string("TOBL5"),selection=cms.untracked.vstring("0x1e01c000-0x1a014000")),
0147     cms.PSet(detSelection = cms.uint32(206),detLabel = cms.string("TOBL6"),selection=cms.untracked.vstring("0x1e01c000-0x1a018000"))
0148 ))
0149 
0150 process.load("DPGAnalysis.SiStripTools.sipixelclustermultiplicityprod_cfi")
0151 
0152 process.load("DPGAnalysis.SiStripTools.ssclusmultinvestigator_cfi")
0153 process.ssclusmultinvestigator.scaleFactor = cms.untracked.int32(2)
0154 process.ssclusmultinvestigator.fillHisto = cms.untracked.bool(True)
0155 process.ssclusmultinvestigator.wantedSubDets.extend(cms.untracked.VPSet(    
0156     cms.PSet(detSelection = cms.uint32(101),detLabel = cms.string("TIBL1"),  binMax = cms.int32(10000)),
0157     cms.PSet(detSelection = cms.uint32(102),detLabel = cms.string("TIBL2"),  binMax = cms.int32(10000)),
0158     cms.PSet(detSelection = cms.uint32(103),detLabel = cms.string("TIBL3"),  binMax = cms.int32(10000)),
0159     cms.PSet(detSelection = cms.uint32(104),detLabel = cms.string("TIBL4"),  binMax = cms.int32(10000)),
0160     cms.PSet(detSelection = cms.uint32(201),detLabel = cms.string("TOBL1"),  binMax = cms.int32(10000)),
0161     cms.PSet(detSelection = cms.uint32(202),detLabel = cms.string("TOBL2"),  binMax = cms.int32(10000)),
0162     cms.PSet(detSelection = cms.uint32(203),detLabel = cms.string("TOBL3"),  binMax = cms.int32(10000)),
0163     cms.PSet(detSelection = cms.uint32(204),detLabel = cms.string("TOBL4"),  binMax = cms.int32(10000)),
0164     cms.PSet(detSelection = cms.uint32(205),detLabel = cms.string("TOBL5"),  binMax = cms.int32(10000)),
0165     cms.PSet(detSelection = cms.uint32(206),detLabel = cms.string("TOBL6"),  binMax = cms.int32(10000))
0166 ))
0167 process.ssclusmultinvestRandom = process.ssclusmultinvestigator.clone()
0168 process.ssclusmultinvestnonoisybins = process.ssclusmultinvestigator.clone()
0169 process.ssclusmultinvestnonoisybinsRandom = process.ssclusmultinvestigator.clone()
0170 
0171 process.load("DPGAnalysis.SiStripTools.spclusmultinvestigator_cfi")
0172 process.spclusmultinvestigator.scaleFactor = cms.untracked.int32(2)
0173 process.spclusmultinvestigator.fillHisto = cms.untracked.bool(True)
0174 process.spclusmultinvestRandom = process.spclusmultinvestigator.clone()
0175 process.spclusmultinvestnonoisybins = process.spclusmultinvestigator.clone()
0176 process.spclusmultinvestnonoisybinsRandom = process.spclusmultinvestigator.clone()
0177 
0178 process.load("DPGAnalysis.SiStripTools.ssclusmulttimecorrelations_cfi")
0179 process.ssclusmulttimecorrelations.historyProduct = cms.InputTag("consecutiveHEs")
0180 process.ssclusmulttimecorrelations.scaleFactors = cms.untracked.vint32()
0181 process.ssclusmulttimecorrelations.runHisto = cms.untracked.bool(True)
0182 process.ssclusmulttimecorrelations.wantedSubDets = cms.untracked.VPSet(    
0183     cms.PSet(detSelection = cms.uint32(0),detLabel = cms.string("TK"),  binMax = cms.int32(9523712/64)),
0184     cms.PSet(detSelection = cms.uint32(3),detLabel = cms.string("TIB"), binMax = cms.int32(1787904/64)),
0185     cms.PSet(detSelection = cms.uint32(4),detLabel = cms.string("TID"), binMax = cms.int32( 565248/64)),
0186     cms.PSet(detSelection = cms.uint32(5),detLabel = cms.string("TOB"), binMax = cms.int32(3303936/64)),
0187     cms.PSet(detSelection = cms.uint32(6),detLabel = cms.string("TEC"), binMax = cms.int32(3866624/64)),
0188     cms.PSet(detSelection = cms.uint32(101),detLabel = cms.string("TIBL1"),  binMax = cms.int32(10000)),
0189     cms.PSet(detSelection = cms.uint32(102),detLabel = cms.string("TIBL2"),  binMax = cms.int32(10000)),
0190     cms.PSet(detSelection = cms.uint32(103),detLabel = cms.string("TIBL3"),  binMax = cms.int32(10000)),
0191     cms.PSet(detSelection = cms.uint32(104),detLabel = cms.string("TIBL4"),  binMax = cms.int32(10000)),
0192     cms.PSet(detSelection = cms.uint32(201),detLabel = cms.string("TOBL1"),  binMax = cms.int32(10000)),
0193     cms.PSet(detSelection = cms.uint32(202),detLabel = cms.string("TOBL2"),  binMax = cms.int32(10000)),
0194     cms.PSet(detSelection = cms.uint32(203),detLabel = cms.string("TOBL3"),  binMax = cms.int32(10000)),
0195     cms.PSet(detSelection = cms.uint32(204),detLabel = cms.string("TOBL4"),  binMax = cms.int32(10000)),
0196     cms.PSet(detSelection = cms.uint32(205),detLabel = cms.string("TOBL5"),  binMax = cms.int32(10000)),
0197     cms.PSet(detSelection = cms.uint32(206),detLabel = cms.string("TOBL6"),  binMax = cms.int32(10000))
0198 )
0199 process.ssclusmulttimecorrRandom = process.ssclusmulttimecorrelations.clone()
0200 process.ssclusmulttimecorrnonoisybins = process.ssclusmulttimecorrelations.clone()
0201 process.ssclusmulttimecorrnonoisybinsRandom = process.ssclusmulttimecorrelations.clone()
0202 
0203 process.load("DPGAnalysis.SiStripTools.spclusmulttimecorrelations_cfi")
0204 process.spclusmulttimecorrelations.historyProduct = cms.InputTag("consecutiveHEs")
0205 process.spclusmulttimecorrelations.scaleFactors = cms.untracked.vint32()
0206 process.spclusmulttimecorrelations.runHisto = cms.untracked.bool(True)
0207 process.spclusmulttimecorrRandom = process.spclusmulttimecorrelations.clone()
0208 process.spclusmulttimecorrnonoisybins = process.spclusmulttimecorrelations.clone()
0209 process.spclusmulttimecorrnonoisybinsRandom = process.spclusmulttimecorrelations.clone()
0210 
0211 #
0212 from HLTrigger.HLTfilters.triggerResultsFilter_cfi import *
0213 process.hltZeroBias = triggerResultsFilter.clone(
0214     triggerConditions = cms.vstring("HLT_ZeroBias_*"),
0215     hltResults = cms.InputTag( "TriggerResults", "", "HLT" ),
0216     l1tResults = cms.InputTag( "" ),
0217     throw = cms.bool(False)
0218                                           )
0219 process.hltRandom = triggerResultsFilter.clone(
0220     triggerConditions = cms.vstring("HLT_Random_*"),
0221     hltResults = cms.InputTag( "TriggerResults", "", "HLT" ),
0222     l1tResults = cms.InputTag( "" ),
0223     throw = cms.bool(False)
0224     )
0225 
0226 
0227 
0228 process.TFileService = cms.Service('TFileService',
0229                                    fileName = cms.string('OOTmultiplicity.root')
0230                                    )
0231 
0232 #----GlobalTag ------------------------
0233 
0234 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0235 from Configuration.AlCa.GlobalTag import GlobalTag
0236 process.GlobalTag = GlobalTag(process.GlobalTag, options.globalTag, '')
0237 
0238 process.pzerobias = cms.Path(process.hltZeroBias
0239                              + process.consecutiveHEs
0240                              + process.APVPhases
0241                              + process.apvcyclephasemonitor
0242                              + process.eventtimedistribution
0243                              + process.ssclustermultprod + process.ssclusmultinvestigator
0244 #                             + process.ssclusmulttimecorrelations 
0245                              + process.spclustermultprod + process.spclusmultinvestigator
0246 #                             + process.spclusmulttimecorrelations
0247                              )
0248 
0249 process.prandom = cms.Path(process.hltRandom
0250                            + process.consecutiveHEs
0251                            + process.APVPhases
0252                            + process.apvcyclephasemonitor
0253                            + process.eventtimedistrRandom
0254                            + process.ssclustermultprod + process.ssclusmultinvestRandom
0255 #                           + process.ssclusmulttimecorrRandom 
0256                            + process.spclustermultprod + process.spclusmultinvestRandom
0257 #                           + process.spclusmulttimecorrRandom
0258                            )
0259 
0260 
0261 process.pnonoisybins = cms.Path(process.hltZeroBias
0262                                 + process.consecutiveHEs
0263                                 + process.APVPhases
0264                                 + process.seqNoNoisyBins
0265                                 + process.ssclustermultprod 
0266                                 + process.spclustermultprod 
0267                                 + process.eventtimedistrnonoisybins
0268                                 + process.ssclusmultinvestnonoisybins
0269 #                                + process.ssclusmulttimecorrnonoisybins
0270                                 + process.spclusmultinvestnonoisybins
0271 #                                + process.spclusmulttimecorrnonoisybins
0272                                 )
0273 
0274 process.pnonoisybinsRandom = cms.Path(process.hltRandom
0275                                       + process.consecutiveHEs
0276                                       + process.APVPhases
0277                                       + process.seqNoNoisyBins
0278                                       + process.ssclustermultprod 
0279                                       + process.spclustermultprod 
0280                                       + process.eventtimedistrnonoisybinsRandom
0281                                       + process.ssclusmultinvestnonoisybinsRandom
0282 #                                      + process.ssclusmulttimecorrnonoisybinsRandom
0283                                       + process.spclusmultinvestnonoisybinsRandom
0284 #                                      + process.spclusmulttimecorrnonoisybinsRandom
0285                                 )
0286