CSCSegment/test/run_on_simdigi_sk.py

0001 ## Process sim digi events with CSC rechit & segment builders - Tim Cox - 11.02.2015
0002 ## This version runs in 7_4_0_preX on a 7_3_0 simulated data DIGI relval sample.
0003 ##     -- USING OLD ALGO 'SK' --
0004 ## Run on  100  events of a 25ns PU TTbar sample
0005
0006 import FWCore.ParameterSet.Config as cms
0007
0008 process = cms.Process("TEST")
0009
0010 process.load("Configuration.StandardSequences.GeometryRecoDB_cff")
0011 process.load("Configuration.StandardSequences.MagneticField_cff")
0012 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0013 ##process.load("Configuration.StandardSequences.RawToDigi_Data_cff")
0014 process.load("Configuration.StandardSequences.Reconstruction_cff")
0015 process.load("Configuration.StandardSequences.EndOfProcess_cff")
0016
0017 process.load("CalibMuon.CSCCalibration.CSCChannelMapper_cfi")
0018 process.load("CalibMuon.CSCCalibration.CSCIndexer_cfi")
0019 process.CSCIndexerESProducer.AlgoName = cms.string("CSCIndexerPostls1")
0020 process.CSCChannelMapperESProducer.AlgoName = cms.string("CSCChannelMapperPostls1")
0021
0022 # --- MATCH GT TO RELEASE AND DATA SAMPLE
0023
0024 process.GlobalTag.globaltag = "MCRUN2_73_V5::All"
0025
0026 # --- NUMBER OF EVENTS
0027
0028 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(100) )
0029
0030 process.options   = cms.untracked.PSet( SkipEvent = cms.untracked.vstring("ProductNotFound") )
0031 process.options   = cms.untracked.PSet( wantSummary = cms.untracked.bool(True) )
0032
0033 ## ttbar+pu is 200 events per file so need 5 for 1000 events
0034 process.source    = cms.Source("PoolSource",
0035     fileNames = cms.untracked.vstring(
0036             "/store/relval/CMSSW_7_3_0/RelValTTbar_13/GEN-SIM-DIGI-RAW-HLTDEBUG/PU25ns_MCRUN2_73_V7_71XGENSIM-v1/00000/044157C8-A181-E411-AC04-002354EF3BD2.root"
0037  ,
0038 "/store/relval/CMSSW_7_3_0/RelValTTbar_13/GEN-SIM-DIGI-RAW-HLTDEBUG/PU25ns_MCRUN2_73_V7_71XGENSIM-v1/00000/0A963931-A181-E411-B4C5-0026189438DC.root",
0039 "/store/relval/CMSSW_7_3_0/RelValTTbar_13/GEN-SIM-DIGI-RAW-HLTDEBUG/PU25ns_MCRUN2_73_V7_71XGENSIM-v1/00000/145BE1DC-A181-E411-816D-0025905A609E.root",
0040 "/store/relval/CMSSW_7_3_0/RelValTTbar_13/GEN-SIM-DIGI-RAW-HLTDEBUG/PU25ns_MCRUN2_73_V7_71XGENSIM-v1/00000/14DDEDD0-A181-E411-9476-0026189438F8.root",
0041 "/store/relval/CMSSW_7_3_0/RelValTTbar_13/GEN-SIM-DIGI-RAW-HLTDEBUG/PU25ns_MCRUN2_73_V7_71XGENSIM-v1/00000/18F85F36-A181-E411-8AF4-0025905B85D6.root"
0042     )
0043 )
0044
0045 # ME1/1A is  u n g a n g e d  postls1
0046
0047 process.CSCGeometryESModule.useGangedStripsInME1a = False
0048 ##process.CSCGeometryESModule.debugV = True
0049 ##process.idealForDigiCSCGeometry.useGangedStripsInME1a = False
0050
0051 # Turn off some flags for CSCRecHitD that are turned ON in default config
0052
0053 process.csc2DRecHits.readBadChannels = cms.bool(False)
0054 process.csc2DRecHits.CSCUseGasGainCorrections = cms.bool(False)
0055 # Already defaults OFF...
0056 ## process.csc2DRecHits.CSCUseTimingCorrections = cms.bool(False)
0057
0058 # Switch input for CSCRecHitD to  s i m u l a t e d  digis
0059
0060 process.csc2DRecHits.wireDigiTag  = cms.InputTag("simMuonCSCDigis","MuonCSCWireDigi")
0061 process.csc2DRecHits.stripDigiTag = cms.InputTag("simMuonCSCDigis","MuonCSCStripDigi")
0062
0063 # -- ACCESSING "DEEP" PARAMETERS OF THE ALGO IS TRICKY
0064 # THE FOLLOWING FOUND BY EXPLORING CONFIG WITH python -i
0065
0066 # switch to CSCSegAlgoSK
0067 process.cscSegments.algo_type = cms.int32(1)
0068
0069 # --- Activate LogVerbatim IN CSCSegment
0070 process.MessageLogger.cerr.enable = False
0071 process.MessageLogger.cout = cms.untracked.PSet(
0072     enable    = cms.untracked.bool(True),
0073     threshold = cms.untracked.string("INFO"),
0074     default   = cms.untracked.PSet( limit = cms.untracked.int32(0)  ),
0075     FwkReport = cms.untracked.PSet( limit = cms.untracked.int32(-1) ),
0076     CSCSegment = cms.untracked.PSet( limit = cms.untracked.int32(-1) )
0077 )
0078
0079 # Path and EndPath def
0080 process.reco = cms.Path(process.csc2DRecHits * process.cscSegments)
0081 process.endjob = cms.EndPath(process.endOfProcess)
0082
0083 # Schedule definition
0084 process.schedule = cms.Schedule(process.reco, process.endjob)