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