SiPixelDigitizer/test/testPixelDigitizer.py

0001 ##############################################################################
0002
0003 import FWCore.ParameterSet.Config as cms
0004
0005 process = cms.Process("Test")
0006
0007 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0008 # process.load("Configuration.StandardSequences.Geometry_cff")
0009 process.load("Configuration.Geometry.GeometryIdeal_cff")
0010 process.load("Configuration.StandardSequences.MagneticField_38T_cff")
0011 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0012 process.load("Configuration.StandardSequences.Services_cff")
0013 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0014
0015 # from v7
0016 #process.load("SimGeneral.MixingModule.pixelDigitizer_cfi")
0017 # process.load("SimTracker.Configuration.SimTracker_cff")
0018 process.load("SimG4Core.Configuration.SimG4Core_cff")
0019
0020 # process.load("SimGeneral.MixingModule.mixNoPU_cfi")
0021 from SimGeneral.MixingModule.aliases_cfi import *
0022 from SimGeneral.MixingModule.mixObjects_cfi import *
0023 # from SimGeneral.MixingModule.digitizers_cfi import *
0024 from SimGeneral.MixingModule.pixelDigitizer_cfi import *
0025 from SimGeneral.MixingModule.stripDigitizer_cfi import *
0026 from SimGeneral.MixingModule.trackingTruthProducer_cfi import *
0027
0028 process.simSiPixelDigis = cms.EDProducer("MixingModule",
0029 #    digitizers = cms.PSet(theDigitizers),
0030 #    digitizers = cms.PSet(
0031 #      mergedtruth = cms.PSet(
0032 #            trackingParticles
0033 #      )
0034 #    ),
0035
0036   digitizers = cms.PSet(
0037    pixel = cms.PSet(
0038     pixelDigitizer
0039    ),
0040 #  strip = cms.PSet(
0041 #    stripDigitizer
0042 #  ),
0043   ),
0044
0045 #theDigitizersValid = cms.PSet(
0046 #  pixel = cms.PSet(
0047 #    pixelDigitizer
0048 #  ),
0049 #  strip = cms.PSet(
0050 #    stripDigitizer
0051 #  ),
0052 #  ecal = cms.PSet(
0053 #    ecalDigitizer
0054 #  ),
0055 #  hcal = cms.PSet(
0056 #    hcalDigitizer
0057 #  ),
0058 #  castor = cms.PSet(
0059 #    castorDigitizer
0060 #  ),
0061 #  mergedtruth = cms.PSet(
0062 #    trackingParticles
0063 #  )
0064 #),
0065
0066     LabelPlayback = cms.string(' '),
0067     maxBunch = cms.int32(3),
0068     minBunch = cms.int32(-5), ## in terms of 25 ns
0069
0070     bunchspace = cms.int32(25),
0071     mixProdStep1 = cms.bool(False),
0072     mixProdStep2 = cms.bool(False),
0073
0074     playback = cms.untracked.bool(False),
0075     useCurrentProcessOnly = cms.bool(False),
0076     mixObjects = cms.PSet(
0077         mixTracks = cms.PSet(
0078             mixSimTracks
0079         ),
0080         mixVertices = cms.PSet(
0081             mixSimVertices
0082         ),
0083         mixSH = cms.PSet(
0084 #            mixPixSimHits
0085 # mixPixSimHits = cms.PSet(
0086     input = cms.VInputTag(cms.InputTag("g4SimHits","TrackerHitsPixelBarrelHighTof"),
0087                           cms.InputTag("g4SimHits","TrackerHitsPixelBarrelLowTof"),
0088                           cms.InputTag("g4SimHits","TrackerHitsPixelEndcapHighTof"),
0089                           cms.InputTag("g4SimHits","TrackerHitsPixelEndcapLowTof"),
0090 #                          cms.InputTag("g4SimHits","TrackerHitsTECHighTof"),
0091 #                          cms.InputTag("g4SimHits","TrackerHitsTECLowTof"),
0092 #                          cms.InputTag("g4SimHits","TrackerHitsTIBHighTof"),
0093 #                          cms.InputTag("g4SimHits","TrackerHitsTIBLowTof"),
0094 #                          cms.InputTag("g4SimHits","TrackerHitsTIDHighTof"),
0095 #                          cms.InputTag("g4SimHits","TrackerHitsTIDLowTof"),
0096 #                          cms.InputTag("g4SimHits","TrackerHitsTOBHighTof"),
0097 #                          cms.InputTag("g4SimHits","TrackerHitsTOBLowTof")
0098     ),
0099     type = cms.string('PSimHit'),
0100     subdets = cms.vstring(
0101         'TrackerHitsPixelBarrelHighTof',
0102         'TrackerHitsPixelBarrelLowTof',
0103         'TrackerHitsPixelEndcapHighTof',
0104         'TrackerHitsPixelEndcapLowTof',
0105 #        'TrackerHitsTECHighTof',
0106 #        'TrackerHitsTECLowTof',
0107 #        'TrackerHitsTIBHighTof',
0108 #        'TrackerHitsTIBLowTof',
0109 #        'TrackerHitsTIDHighTof',
0110 #        'TrackerHitsTIDLowTof',
0111 #        'TrackerHitsTOBHighTof',
0112 #        'TrackerHitsTOBLowTof'
0113     ),
0114     crossingFrames = cms.untracked.vstring(),
0115 #        'MuonCSCHits',
0116 #        'MuonDTHits',
0117 #        'MuonRPCHits'),
0118 #)
0119         ),
0120         mixHepMC = cms.PSet(
0121             mixHepMCProducts
0122         )
0123     )
0124 )
0125
0126 process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
0127 #     simMuonCSCDigis = cms.PSet(
0128 #        initialSeed = cms.untracked.uint32(1234567),
0129 #        engineName = cms.untracked.string('TRandom3')
0130 #    ),
0131      simSiPixelDigis = cms.PSet(
0132         initialSeed = cms.untracked.uint32(1234567),
0133         engineName = cms.untracked.string('TRandom3')
0134    )
0135 )
0136
0137
0138 process.maxEvents = cms.untracked.PSet(
0139     input = cms.untracked.int32(-1)
0140 )
0141
0142 process.source = cms.Source("PoolSource", fileNames = cms.untracked.vstring(
0143        'file:/afs/cern.ch/work/d/dkotlins/public//MC/mu/pt100_71_pre5/simhits/simHits2.root'
0144   )
0145 )
0146
0147 # Choose the global tag here:
0148 # for v7.0
0149 #process.GlobalTag.globaltag = 'MC_70_V1::All'
0150 #process.GlobalTag.globaltag = "START70_V1::All"
0151 #process.GlobalTag.globaltag = "START71_V1::All"
0152 #process.GlobalTag.globaltag = 'MC_71_V1::All'
0153 process.GlobalTag.globaltag = 'POSTLS171_V1::All'
0154 #process.GlobalTag.globaltag = "PRE_MC_71_V2::All"
0155
0156 process.o1 = cms.OutputModule("PoolOutputModule",
0157             outputCommands = cms.untracked.vstring('drop *','keep *_*_*_Test'),
0158       fileName = cms.untracked.string('file:/afs/cern.ch/work/d/dkotlins/public/MC/mu/pt100_71_pre7/digis/digis2_postls171.root')
0159 #      fileName = cms.untracked.string('file:dummy.root')
0160 )
0161
0162 process.g4SimHits.Generator.HepMCProductLabel = 'source'
0163 process.simSiPixelDigis.digitizers.pixel.AddPixelInefficiency = cms.bool(False)
0164
0165 # modify digitizer parameters
0166 #process.simSiPixelDigis.ThresholdInElectrons_BPix = 3500.0
0167 #process.mix.digitizers.pixel.ThresholdInElectrons_BPix = 3500.0
0168
0169 #This process is to run the digitizer, pixel gitizer is now clled by the mix module
0170 process.p1 = cms.Path(process.simSiPixelDigis)
0171
0172 process.outpath = cms.EndPath(process.o1)
0173
0174