From a890f1301548918950a5810016a19077804b2101 Mon Sep 17 00:00:00 2001 From: Francesco Fiori Date: Tue, 20 Oct 2015 17:46:42 +0200 Subject: [PATCH 1/2] ready to use cfg for shifters --- .../SiStripTools/test/APVShot_Selector_cfg.py | 88 +++++++++++++ .../ClusterMultiplicityRegionSelector_cfg.py | 121 ++++++++++++++++++ 2 files changed, 209 insertions(+) create mode 100644 DPGAnalysis/SiStripTools/test/APVShot_Selector_cfg.py create mode 100644 DPGAnalysis/SiStripTools/test/ClusterMultiplicityRegionSelector_cfg.py diff --git a/DPGAnalysis/SiStripTools/test/APVShot_Selector_cfg.py b/DPGAnalysis/SiStripTools/test/APVShot_Selector_cfg.py new file mode 100644 index 0000000000000..b772ac6af93ad --- /dev/null +++ b/DPGAnalysis/SiStripTools/test/APVShot_Selector_cfg.py @@ -0,0 +1,88 @@ +import FWCore.ParameterSet.Config as cms +import FWCore.ParameterSet.VarParsing as VarParsing + +process = cms.Process("APVShotAnalyzer") + +#prepare options + +process.DQMStore=cms.Service("DQMStore") +process.TkDetMap=cms.Service("TkDetMap") +process.SiStripDetInfoFileReader = cms.Service("SiStripDetInfoFileReader") + +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.parseArguments() + +# + +process.options = cms.untracked.PSet( + wantSummary = cms.untracked.bool(False), + fileMode = cms.untracked.string("FULLMERGE") + ) + +process.load("FWCore.MessageService.MessageLogger_cfi") + + +process.MessageLogger.debugModules=cms.untracked.vstring("apvshotfilter") +#------------------------------------------------------------------ + +process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(10000) ) + +process.source = cms.Source("PoolSource", + fileNames = cms.untracked.vstring('/store/data/Run2015D/ZeroBias/RAW/v1/000/257/805/00000/AA3FF22F-E766-E511-B8F2-02163E01234B.root','/store/data/Run2015D/ZeroBias/RAW/v1/000/257/805/00000/FAEE2C85-C666-E511-8F92-02163E014708.root'), +# skipBadFiles = cms.untracked.bool(True), + inputCommands = cms.untracked.vstring("keep *", "drop *_MEtoEDMConverter_*_*") + ) + +#-------------------------------------- +process.load("Configuration.StandardSequences.RawToDigi_Data_cff") +process.load("Configuration.StandardSequences.MagneticField_AutoFromDBCurrent_cff") +process.load("Configuration.StandardSequences.GeometryRecoDB_cff") +process.load("Configuration.StandardSequences.Reconstruction_cff") + +process.froml1abcHEs = cms.EDProducer("EventWithHistoryProducerFromL1ABC", + l1ABCCollection=cms.InputTag("scalersRawToDigi") + ) +process.load("DPGAnalysis.SiStripTools.apvcyclephaseproducerfroml1tsDB_cfi") +process.load("DPGAnalysis.SiStripTools.eventtimedistribution_cfi") + +process.seqEventHistoryReco = process.seqEventHistoryReco = cms.Sequence(process.froml1abcHEs + process.APVPhases) +process.seqEventHistory = cms.Sequence(process.eventtimedistribution) + +process.eventtimedistribution.historyProduct = cms.InputTag("froml1abcHEs") + + +process.load("DPGAnalysis.SiStripTools.apvshotsfilter_cfi") +process.apvshotsfilter.useCabling = cms.untracked.bool(True) + +process.load("DPGAnalysis.SiStripTools.apvshotsanalyzer_cfi") +process.apvshotsanalyzer.historyProduct = cms.InputTag("froml1abcHEs") +process.apvshotsanalyzer.useCabling = cms.untracked.bool(True) + +#----GlobalTag ------------------------ + +process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_condDBv2_cff") +from Configuration.AlCa.GlobalTag_condDBv2 import GlobalTag +process.GlobalTag = GlobalTag(process.GlobalTag,'74X_dataRun2_Prompt_v2' , '') + + + +process.p0 = cms.Path( + process.siStripDigis + process.siStripZeroSuppression + + process.scalersRawToDigi + + process.seqEventHistoryReco + + process.seqEventHistory + + process.apvshotsfilter + + process.apvshotsanalyzer + ) + +process.TFileService = cms.Service('TFileService', + fileName = cms.string('APVShotAnalyzer.root') + ) + diff --git a/DPGAnalysis/SiStripTools/test/ClusterMultiplicityRegionSelector_cfg.py b/DPGAnalysis/SiStripTools/test/ClusterMultiplicityRegionSelector_cfg.py new file mode 100644 index 0000000000000..48aa29e02f12a --- /dev/null +++ b/DPGAnalysis/SiStripTools/test/ClusterMultiplicityRegionSelector_cfg.py @@ -0,0 +1,121 @@ +import FWCore.ParameterSet.Config as cms + +process = cms.Process("BeamBackground") + +process.load('Configuration.StandardSequences.Services_cff') +process.load('FWCore.MessageService.MessageLogger_cfi') + +process.maxEvents = cms.untracked.PSet( + input = cms.untracked.int32(500) +) + + +process.source = cms.Source("PoolSource", + fileNames = cms.untracked.vstring('/store/data/Run2015C/ZeroBias/RECO/PromptReco-v1/000/254/906/00000/226DA26D-D84B-E511-9CFD-02163E011C96.root'), + secondaryFileNames = cms.untracked.vstring() +) + + +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.load("DPGAnalysis.SiStripTools.sipixelclustermultiplicityprod_cfi") +process.load("DPGAnalysis.SiStripTools.sistripclustermultiplicityprod_cfi") + +ssclustermultprod = cms.EDProducer("SiStripClusterMultiplicityProducer", + clusterdigiCollection = cms.InputTag("siStripClusters"), + wantedSubDets = cms.VPSet( + cms.PSet(detSelection = cms.uint32(0),detLabel = cms.string("TK")), + cms.PSet(detSelection = cms.uint32(3),detLabel = cms.string("TIB")), + cms.PSet(detSelection = cms.uint32(4),detLabel = cms.string("TID")), + cms.PSet(detSelection = cms.uint32(5),detLabel = cms.string("TOB")), + cms.PSet(detSelection = cms.uint32(6),detLabel = cms.string("TEC")) + ) + ) + + +spclustermultprod = cms.EDProducer("SiPixelClusterMultiplicityProducer", + clusterdigiCollection = cms.InputTag("siPixelClusters"), + wantedSubDets = cms.VPSet( + cms.PSet(detSelection = cms.uint32(0),detLabel = cms.string("Pixel")), + cms.PSet(detSelection = cms.uint32(1),detLabel = cms.string("BPIX")), + cms.PSet(detSelection = cms.uint32(2),detLabel = cms.string("FPIX")) + ) + ) + + +process.load("DPGAnalysis.SiStripTools.multiplicitycorr_cfi") + +process.multiplicitycorr.correlationConfigurations = cms.VPSet( + cms.PSet(xMultiplicityMap = cms.InputTag("ssclustermultprod"), + xDetSelection = cms.uint32(0), xDetLabel = cms.string("TK"), xBins = cms.uint32(3000), xMax=cms.double(100000), + yMultiplicityMap = cms.InputTag("spclustermultprod"), + yDetSelection = cms.uint32(0), yDetLabel = cms.string("Pixel"), yBins = cms.uint32(1000), yMax=cms.double(30000), + rBins = cms.uint32(200), scaleFactor = cms.untracked.double(10.4),#10.4 for 25ns // 7.7 for 50 ns + runHisto=cms.bool(True),runHistoBXProfile=cms.bool(True),runHistoBX=cms.bool(True),runHisto2D=cms.bool(True)) +) + + +process.load("DPGAnalysis.SiStripTools.eventtimedistribution_cfi") + + +import DPGAnalysis.SiStripTools.apvcyclephaseproducerfroml1tsDB_cfi +process.APVPhases = DPGAnalysis.SiStripTools.apvcyclephaseproducerfroml1tsDB_cfi.APVPhases + +process.load("DPGAnalysis.SiStripTools.eventwithhistoryproducerfroml1abc_cfi") + +process.load("DPGAnalysis.SiStripTools.bysipixelvssistripclustmulteventfilter_cfi") + +# mult2= Number of Strip Clusters; mult1= Number of Pixel Clusters + +#filter for BeamBackground events +process.offdiagonal = process.bysipixelvssistripclustmulteventfilter.clone(cut=cms.string("( mult2 < 7.4*(mult1-300))")) #7.4 for 25ns // 5.5 for 50 ns +process.NoZeroSClusters = process.bysipixelvssistripclustmulteventfilter.clone(cut=cms.string("( mult2 > 500)")) + +#Filter for high strip noise +#process.offdiagonal = process.bysipixelvssistripclustmulteventfilter.clone(cut=cms.string("( mult2 > 7.4*(mult1+300))")) #7.4 for 25ns // 5.5 for 50 ns +#process.NoZeroSClusters = process.bysipixelvssistripclustmulteventfilter.clone(cut=cms.string("( mult2 > 26000)")) + +#Filter for strip noise +#process.offdiagonal = process.bysipixelvssistripclustmulteventfilter.clone(cut=cms.string("( mult2 > 7.4*(mult1+300))")) #7.4 for 25ns // 5.5 for 50 ns +#process.NoZeroSClusters = process.bysipixelvssistripclustmulteventfilter.clone(cut=cms.string("( mult2 < 26000)")) + +#Filter for main diagonal +#process.offdiagonal = process.bysipixelvssistripclustmulteventfilter.clone(cut=cms.string("( mult2 > 7.4*(mult1-300))")) #7.4 for 25ns // 5.5 for 50 ns +#process.offdiagonal1 = process.bysipixelvssistripclustmulteventfilter.clone(cut=cms.string("( mult2 < 7.4*(mult1+300))")) #7.4 for 25ns // 5.5 for 50 ns + + + + +process.multiplicitycorrAfter=process.multiplicitycorr.clone() + +process.multiplicitycorrAfter.correlationConfigurations = cms.VPSet( + cms.PSet(xMultiplicityMap = cms.InputTag("ssclustermultprod"), + xDetSelection = cms.uint32(0), xDetLabel = cms.string("TK"), xBins = cms.uint32(3000), xMax=cms.double(100000), + yMultiplicityMap = cms.InputTag("spclustermultprod"), + yDetSelection = cms.uint32(0), yDetLabel = cms.string("Pixel"), yBins = cms.uint32(1000), yMax=cms.double(30000), + rBins = cms.uint32(200), scaleFactor = cms.untracked.double(10.4),#10.4 for 25 ns + runHisto=cms.bool(True),runHistoBXProfile=cms.bool(True),runHistoBX=cms.bool(True),runHisto2D=cms.bool(True)) +) + +process.eventtimedistributionAfter= process.eventtimedistribution.clone() + +process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_condDBv2_cff") + +from Configuration.AlCa.GlobalTag_condDBv2 import GlobalTag + +process.GlobalTag = GlobalTag(process.GlobalTag, '74X_dataRun2_Prompt_v2', '') + +process.TFileService = cms.Service('TFileService', + fileName = cms.string('BeamBackground.root') + ) + +process.MainSeq= cms.Sequence(process.hltSelection+process.consecutiveHEs+process.APVPhases+process.ssclustermultprod+process.spclustermultprod+process.eventtimedistribution+process.multiplicitycorr+process.NoZeroSClusters+process.offdiagonal+process.eventtimedistributionAfter+process.multiplicitycorrAfter) + +process.p0 = cms.Path(process.MainSeq) From 5dec854899b1c53fac4857a1fa9ecb60c0d8b002 Mon Sep 17 00:00:00 2001 From: Francesco Fiori Date: Wed, 21 Oct 2015 13:45:04 +0200 Subject: [PATCH 2/2] add var parsing --- .../SiStripTools/test/APVShot_Selector_cfg.py | 4 ++-- .../ClusterMultiplicityRegionSelector_cfg.py | 16 ++++++++++++++-- 2 files changed, 16 insertions(+), 4 deletions(-) diff --git a/DPGAnalysis/SiStripTools/test/APVShot_Selector_cfg.py b/DPGAnalysis/SiStripTools/test/APVShot_Selector_cfg.py index b772ac6af93ad..e7ef76c32e646 100644 --- a/DPGAnalysis/SiStripTools/test/APVShot_Selector_cfg.py +++ b/DPGAnalysis/SiStripTools/test/APVShot_Selector_cfg.py @@ -35,7 +35,7 @@ process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(10000) ) process.source = cms.Source("PoolSource", - fileNames = cms.untracked.vstring('/store/data/Run2015D/ZeroBias/RAW/v1/000/257/805/00000/AA3FF22F-E766-E511-B8F2-02163E01234B.root','/store/data/Run2015D/ZeroBias/RAW/v1/000/257/805/00000/FAEE2C85-C666-E511-8F92-02163E014708.root'), + fileNames = cms.untracked.vstring(options.inputFiles), # skipBadFiles = cms.untracked.bool(True), inputCommands = cms.untracked.vstring("keep *", "drop *_MEtoEDMConverter_*_*") ) @@ -69,7 +69,7 @@ process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_condDBv2_cff") from Configuration.AlCa.GlobalTag_condDBv2 import GlobalTag -process.GlobalTag = GlobalTag(process.GlobalTag,'74X_dataRun2_Prompt_v2' , '') +process.GlobalTag = GlobalTag(process.GlobalTag,options.globalTag, '') diff --git a/DPGAnalysis/SiStripTools/test/ClusterMultiplicityRegionSelector_cfg.py b/DPGAnalysis/SiStripTools/test/ClusterMultiplicityRegionSelector_cfg.py index 48aa29e02f12a..11b3760572683 100644 --- a/DPGAnalysis/SiStripTools/test/ClusterMultiplicityRegionSelector_cfg.py +++ b/DPGAnalysis/SiStripTools/test/ClusterMultiplicityRegionSelector_cfg.py @@ -1,17 +1,29 @@ import FWCore.ParameterSet.Config as cms +import FWCore.ParameterSet.VarParsing as VarParsing process = cms.Process("BeamBackground") +options = VarParsing.VarParsing("analysis") + process.load('Configuration.StandardSequences.Services_cff') process.load('FWCore.MessageService.MessageLogger_cfi') +options.register ('globalTag', + "DONOTEXIST", + VarParsing.VarParsing.multiplicity.singleton, # singleton or list + VarParsing.VarParsing.varType.string, # string, int, or float + "GlobalTag") + +options.parseArguments() + + process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(500) ) process.source = cms.Source("PoolSource", - fileNames = cms.untracked.vstring('/store/data/Run2015C/ZeroBias/RECO/PromptReco-v1/000/254/906/00000/226DA26D-D84B-E511-9CFD-02163E011C96.root'), + fileNames = cms.untracked.vstring(options.inputFiles), secondaryFileNames = cms.untracked.vstring() ) @@ -110,7 +122,7 @@ from Configuration.AlCa.GlobalTag_condDBv2 import GlobalTag -process.GlobalTag = GlobalTag(process.GlobalTag, '74X_dataRun2_Prompt_v2', '') +process.GlobalTag = GlobalTag(process.GlobalTag,options.globalTag, '') process.TFileService = cms.Service('TFileService', fileName = cms.string('BeamBackground.root')