python/test/Validate_DigiToRawToClusters_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 Mode = str("ZS")    # Options: "ZS", "VR", "PR", "FK"
0004 Write = bool(False) # Write output to disk
0005
0006 process = cms.Process("DigiToRawToClusters")
0007
0008 # ---- Data source ----
0009 process.source = cms.Source(
0010    "PoolSource",
0011    fileNames = cms.untracked.vstring(
0012    '/store/relval/CMSSW_3_1_1/RelValQCD_FlatPt_15_3000/GEN-SIM-DIGI-RAW-HLTDEBUG/MC_31X_V2-v1/0002/ECAD7ED7-966B-DE11-B4FE-000423D99CEE.root'
0013     )
0014 )
0015
0016 # ---- Services ----
0017 process.load("DQM.SiStripCommon.MessageLogger_cfi")
0018 process.MessageLogger.debugModules = cms.untracked.vstring()
0019 process.Timing = cms.Service("Timing")
0020 process.Tracer = cms.Service(
0021     "Tracer",
0022     sourceSeed = cms.untracked.string("$$")
0023     )
0024
0025 # ---- Conditions ----
0026 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0027 process.GlobalTag.globaltag = "MC_31X_V2::All"
0028
0029 # ---- Region cabling ----
0030 process.load('Configuration.StandardSequences.Geometry_cff')
0031 process.load('Configuration.StandardSequences.MagneticField_cff')
0032 process.load('CalibTracker.SiStripESProducers.SiStripRegionConnectivity_cfi')
0033
0034 # ---- Reference clusters ----
0035 from RecoLocalTracker.SiStripClusterizer.SiStripClusterizer_cfi import *
0036 process.referenceSiStripClusters = siStripClusters.clone()
0037 process.referenceSiStripClusters.DigiProducersList = cms.VInputTag(cms.InputTag('simSiStripDigis:ZeroSuppressed'))
0038
0039 # ---- DigiToRaw ----
0040 process.load("EventFilter.SiStripRawToDigi.SiStripDigiToRaw_cfi")
0041 process.SiStripDigiToRaw.InputDigis = cms.InputTag('simSiStripDigis', 'ZeroSuppressed')
0042
0043 # ---- RawToClusters ----
0044 process.load('EventFilter.SiStripRawToDigi.SiStripRawToClusters_cfi')
0045 process.SiStripRawToClustersFacility.ProductLabel = cms.InputTag("SiStripDigiToRaw")
0046 process.load('EventFilter.SiStripRawToDigi.SiStripRawToClustersRoI_cfi')
0047 process.SiStripRoI.SiStripLazyGetter = cms.InputTag("SiStripRawToClustersFacility")
0048 process.load('EventFilter.SiStripRawToDigi.test.SiStripClustersDSVBuilder_cfi')
0049 process.siStripClustersDSV.SiStripLazyGetter = cms.InputTag("SiStripRawToClustersFacility")
0050 process.siStripClustersDSV.SiStripRefGetter = cms.InputTag("SiStripRoI")
0051 process.siStripClustersDSV.DetSetVectorNew = True
0052 process.SiStripRawToClusters = cms.Sequence( process.SiStripRawToClustersFacility * process.SiStripRoI * process.siStripClustersDSV )
0053
0054 # ---- Validation ----
0055 process.load('EventFilter.SiStripRawToDigi.test.SiStripClusterValidator_cfi')
0056 process.ValidateSiStripClusters.Collection1 = cms.untracked.InputTag("referenceSiStripClusters")
0057 process.ValidateSiStripClusters.Collection2 = cms.untracked.InputTag("siStripClustersDSV")
0058 process.ValidateSiStripClusters.DetSetVectorNew = True
0059
0060 # ----- FedReadoutMode -----
0061 if Mode == str("ZS") :
0062     process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(1000) )
0063     process.SiStripDigiToRaw.FedReadoutMode = cms.string('ZERO_SUPPRESSED')
0064     process.SiStripDigiToRaw.PacketCode = cms.string('ZERO_SUPPRESSED')
0065 elif Mode == str("VR") :
0066     process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(10) )
0067     process.SiStripDigiToRaw.FedReadoutMode = cms.string('VIRGIN_RAW')
0068     process.SiStripDigiToRaw.PacketCode = cms.string('VIRGIN_RAW')
0069 elif Mode == str("PR") :
0070     process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(10) )
0071     process.SiStripDigiToRaw.FedReadoutMode = cms.string('PROCESSED_RAW')
0072     process.SiStripDigiToRaw.PacketCode = cms.string('PROCESSED_RAW')
0073 else :
0074     print("UNKNOWN FED READOUT MODE!")
0075     import sys
0076     sys.exit()
0077
0078 # ---- Sequence ----
0079 process.p = cms.Path(
0080     process.referenceSiStripClusters *
0081     process.SiStripDigiToRaw *
0082     process.SiStripRawToClusters *
0083     process.ValidateSiStripClusters
0084     )
0085
0086 # ----- WriteToDisk -----
0087 process.output = cms.OutputModule(
0088     "PoolOutputModule",
0089     fileName = cms.untracked.string('output.root'),
0090     outputCommands = cms.untracked.vstring(
0091     'drop *',
0092     'keep SiStrip*_simSiStripDigis_*_*', # (to drop SimLinks)
0093     'keep *_*_*_DigiToRawToClusters'
0094     )
0095     )
0096 process.output.fileName = "DigiToRawToClusters"+Mode+".root"
0097 if Write == bool(True) :
0098     process.e = cms.EndPath( process.output )
0099 else :
0100     print("Event content not written to disk!")
0101