SiStripZeroSuppression/test/SiStripZS_EmulateHybrid_RealData.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process('TEST1')
0004
0005 # import of standard configurations
0006 process.load('Configuration.StandardSequences.Services_cff')
0007 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0008 process.load('FWCore.MessageService.MessageLogger_cfi')
0009 process.load('Configuration.StandardSequences.GeometryDB_cff')
0010 process.load('Configuration.StandardSequences.MagneticField_38T_cff')
0011 process.load('Configuration.StandardSequences.RawToDigi_Data_cff')
0012 process.load('Configuration.StandardSequences.L1Reco_cff')
0013 process.load('Configuration.StandardSequences.Reconstruction_cff')
0014 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0015 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0016 process.load('Configuration.EventContent.EventContentHeavyIons_cff')
0017
0018 process.configurationMetadata = cms.untracked.PSet(
0019     annotation = cms.untracked.string('hiRecoDM nevts:2'),
0020     name = cms.untracked.string('PyReleaseValidation')
0021 )
0022 process.maxEvents = cms.untracked.PSet(
0023     input = cms.untracked.int32(1)
0024 )
0025
0026 # Input source
0027 process.source = cms.Source("PoolSource",
0028     fileNames = cms.untracked.vstring(
0029     #'file:/data/abaty/VirginRaw_CentralitySkims/VirginRAW_2010_HICorePhysics_SKIM_Cent_0_5_1.root',
0030     #'file:/data/abaty/HLT_Emulated_2010Data/outputHIPhysicsVirginRaw.root'
0031     #'file:/data/abaty/VirginRaw_CentralitySkims/VirginRAW_2010_HICorePhysics_SKIM_Cent_0_5_10.root',
0032     #'/store/hidata/HIRun2015/HITrackerVirginRaw/RAW/v1/000/263/400/00000/40322926-4AA3-E511-95F7-02163E0146A8.root'
0033      'file:/data/TrackerStudies/2015_VR_forBaselineFollowerStudies_Aug2017.root'
0034     # 'root://cms-xrd-global.cern.ch////store/hidata/HIRun2015/HITrackerVirginRaw/RAW/v1/000/262/495/00000/14634975-5493-E511-AC10-02163E0138B3.root'
0035    )
0036 )
0037
0038 #process.source = cms.Source("NewEventStreamFileReader",
0039 #fileNames = cms.untracked.vstring('root://eoscms.cern.ch//store/t0streamer/Data/HIPhysicsVirginRaw/000/262/296/run262296_ls0223_streamHIPhysicsVirginRaw_StorageManager.dat'
0040 #)
0041 #)
0042
0043 # Output definition
0044 process.RECOoutput = cms.OutputModule("PoolOutputModule",
0045     splitLevel = cms.untracked.int32(0),
0046     fileName = cms.untracked.string('RECO.root'),
0047     dataset = cms.untracked.PSet(
0048         filterName = cms.untracked.string(''),
0049         dataTier = cms.untracked.string('GEN-SIM-RECO')
0050    ),
0051
0052     outputCommands = cms.untracked.vstring('keep FEDRawDataCollection_myRawDataCollector_*_*')
0053
0054 )
0055
0056
0057
0058 process.load('Configuration/StandardSequences/FrontierConditions_GlobalTag_cff')
0059 process.GlobalTag.globaltag = '75X_dataRun2_HLTHI_v4'
0060
0061 ## Offline Silicon Tracker Zero Suppression
0062 from RecoLocalTracker.SiStripZeroSuppression.DefaultAlgorithms_cff import *
0063 process.siStripZeroSuppression.produceRawDigis = False
0064 process.siStripZeroSuppression.produceCalculatedBaseline = False
0065 process.siStripZeroSuppression.produceBaselinePoints = False
0066 process.siStripZeroSuppression.storeCM = True
0067 process.siStripZeroSuppression.produceHybridFormat = True
0068 process.siStripZeroSuppression.Algorithms.APVInspectMode = "HybridEmulation"
0069 process.siStripZeroSuppression.Algorithms.CommonModeNoiseSubtractionMode = 'Median'
0070 process.siStripZeroSuppression.Algorithms.MeanCM = 512
0071 process.siStripZeroSuppression.Algorithms.DeltaCMThreshold = 20
0072 process.siStripZeroSuppression.Algorithms.Use10bitsTruncation = True
0073
0074 process.TFileService = cms.Service("TFileService",
0075         fileName=cms.string("Baselines.root"))
0076
0077 process.hybridAna = cms.EDAnalyzer("SiStripHybridFormatAnalyzer",
0078
0079     srcDigis =  cms.InputTag('siStripZeroSuppression','VirginRaw'),
0080     srcAPVCM =  cms.InputTag('siStripZeroSuppression','APVCMVirginRaw'),
0081     nModuletoDisplay = cms.uint32(10000),
0082     plotAPVCM   = cms.bool(True)
0083 )
0084
0085 process.siStripDigiToRawData = cms.EDProducer(
0086     "SiStripDigiToRawModule",
0087     InputModuleLabel = cms.string('siStripZeroSuppression'),
0088     InputDigiLabel = cms.string('VirginRaw'),
0089     FedReadoutMode = cms.string('ZERO_SUPPRESSED'),
0090     UseFedKey = cms.bool(False),
0091     UseWrongDigiType = cms.bool(False),
0092     CopyBufferHeader = cms.bool(False),
0093     RawDataTag = cms.InputTag('rawDataCollector')
0094 )
0095
0096 process.myRawDataCollector = cms.EDProducer("RawDataCollectorByLabel",
0097     verbose = cms.untracked.int32(0),     # 0 = quiet, 1 = collection list, 2 = FED list
0098     RawCollectionList = cms.VInputTag( #cms.InputTag('hcalRawData'),
0099                                        cms.InputTag('siStripDigiToRawData'),
0100                    #                     cms.InputTag('rawDataRepacker')
0101                                        cms.InputTag('rawDataCollector')
0102                    #                    cms.InputTag('siPixelRawData'),
0103     ),
0104 )
0105
0106 # Path and EndPath definitions
0107 process.raw2digi_step = cms.Path(process.siStripDigis)
0108 #process.reconstruction_step = cms.Path(process.striptrackerlocalreco+process.moddedZS+process.moddedClust+process.baselineAna+process.moddedbaselineAna+process.clusterMatching)
0109 process.reconstruction_step = cms.Path(process.striptrackerlocalreco*process.hybridAna*process.siStripDigiToRawData*process.myRawDataCollector)
0110 #process.reconstruction_step = cms.Path(process.striptrackerlocalreco);
0111 #process.endjob_step = cms.EndPath(process.endOfProcess)
0112 process.RECOoutput_step = cms.EndPath(process.RECOoutput)
0113
0114
0115 # Schedule definition
0116 process.schedule = cms.Schedule(process.raw2digi_step,process.reconstruction_step, process.RECOoutput_step)