CastorSim/data/CastorDigiValidation.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("CastorTest")
0004
0005 process.load("SimGeneral.HepPDTESSource.pdt_cfi")
0006 process.load("IOMC.EventVertexGenerators.VtxSmearedGauss_cfi")
0007 process.load("Configuration.StandardSequences.GeometryExtended_cff")
0008 process.load("Configuration.EventContent.EventContent_cff")
0009
0010 process.load("SimG4Core.Application.g4SimHits_cfi")
0011 process.load("Configuration.StandardSequences.Generator_cff")
0012 process.load("SimGeneral.MixingModule.mixNoPU_cfi")
0013 process.load("Geometry.CaloEventSetup.CaloGeometry_cff")
0014
0015 process.load("SimCalorimetry.CastorSim.castordigi_cfi")
0016
0017 process.load("RecoLocalCalo.CastorReco.CastorSimpleReconstructor_cfi")
0018
0019
0020 process.load("CondCore.DBCommon.CondDBSetup_cfi")
0021 process.load("CondCore.DBCommon.CondDBCommon_cfi")
0022
0023 process.load('RecoLocalCalo.Castor.Castor_cff')
0024 process.castor_db_producer = cms.ESProducer("CastorDbProducer")
0025
0026 process.MessageLogger = cms.Service("MessageLogger",
0027     destinations = cms.untracked.vstring('cout')
0028 )
0029
0030 process.load("Configuration.StandardSequences.SimulationRandomNumberGeneratorSeeds_cff")
0031
0032 process.RandomNumberGeneratorService.theSource.initialSeed = 15875319
0033 process.RandomNumberGeneratorService.generator.initialSeed = 9792301
0034 process.RandomNumberGeneratorService.VtxSmeared.initialSeed = 2263014
0035 process.RandomNumberGeneratorService.g4SimHits.initialSeed = 124
0036 process.RandomNumberGeneratorService.mix.initialSeed = 22159
0037 process.RandomNumberGeneratorService.simCastorDigis.initialSeed = 210
0038
0039
0040 process.maxEvents = cms.untracked.PSet(
0041     input = cms.untracked.int32(10)
0042 )
0043
0044 process.source = cms.Source("EmptySource")
0045
0046 process.generator = cms.EDProducer("FlatRandomEGunProducer",
0047     PGunParameters = cms.PSet(
0048         PartID = cms.vint32(13),
0049         MinEta = cms.double(5.2),
0050         MaxEta = cms.double(6.6),
0051         MinPhi = cms.double(-3.14),
0052         MaxPhi = cms.double(3.14),
0053         MinE = cms.double(50.00),
0054         MaxE = cms.double(50.00)
0055     ),
0056     AddAntiParticle = cms.bool(True),
0057     Verbosity = cms.untracked.int32(0)
0058 )
0059
0060
0061 process.es_pool = cms.ESSource( "PoolDBESSource",
0062      process.CondDBSetup,
0063      timetype = cms.string('runnumber'),
0064 #   connect = cms.string('frontier://cmsfrontier.cern.ch:8000/FrontierPrep/CMS_COND_30X_HCAL'),
0065     connect = cms.string('sqlite_file:testExample.db'),
0066      authenticationMethod = cms.untracked.uint32(0),
0067      toGet = cms.VPSet(
0068          cms.PSet(
0069              record = cms.string('CastorPedestalsRcd'),
0070              tag = cms.string('castor_pedestals_v1.0_test')
0071              ),
0072          cms.PSet(
0073              record = cms.string('CastorPedestalWidthsRcd'),
0074              tag = cms.string('castor_widths_v1.0_test')
0075              ),
0076          cms.PSet(
0077              record = cms.string('CastorGainsRcd'),
0078              tag = cms.string('castor_gains_v1.0_test')
0079              ),
0080          cms.PSet(
0081              record = cms.string('CastorQIEDataRcd'),
0082              tag = cms.string('castor_qie_v1.0_test')
0083              ),
0084          cms.PSet(
0085              record = cms.string('CastorElectronicsMapRcd'),
0086              tag = cms.string('castor_emap_v1.0_test')
0087               ),
0088          cms.PSet(
0089              record = cms.string('CastorChannelQualityRcd'),
0090              tag = cms.string('castor_channelstatus_v1.0_test')
0091
0092              )
0093
0094      )
0095 )
0096 process.es_hardcode = cms.ESSource("CastorHardcodeCalibrations",
0097      toGet = cms.untracked.vstring('GainWidths')
0098  )
0099
0100 process.common_maximum_timex = cms.PSet( # need to be localy redefined
0101    MaxTrackTime  = cms.double(500.0),  # need to be localy redefined
0102    MaxTimeNames  = cms.vstring(), # need to be localy redefined
0103    MaxTrackTimes = cms.vdouble()  # need to be localy redefined
0104 )
0105
0106
0107
0108 process.CaloSD = cms.PSet(
0109     DetailedTiming = cms.bool(False),
0110     EminTrack = cms.double(1.0),
0111     Verbosity = cms.int32(0),
0112     UseMap = cms.bool(True),
0113     CheckHits = cms.int32(25)
0114 )
0115
0116
0117
0118
0119
0120 process.o1 = cms.OutputModule("PoolOutputModule",
0121     fileName = cms.untracked.string('gen-sim-digi-reco_muon.root')
0122 )
0123
0124 process.p1 = cms.Path(process.generator*process.pgen*process.VtxSmeared*process.g4SimHits*process.mix*process.simCastorDigis*process.castorreco)
0125 #*process.CastorFullReco)
0126 process.outpath=cms.EndPath(process.o1)
0127
0128 process.g4SimHits.UseMagneticField = False
0129 process.g4SimHits.Physics.DefaultCutValue = 10.
0130
0131 process.g4SimHits.Generator.ApplyEtaCuts = False
0132 process.g4SimHits.UseMagneticField = False
0133
0134 process.g4SimHits.CaloTrkProcessing.TestBeam = True
0135
0136 process.g4SimHits.StackingAction = cms.PSet(
0137    process.common_heavy_suppression,
0138    process.common_maximum_timex,        # need to be localy redefined
0139    TrackNeutrino = cms.bool(False),
0140    KillHeavy     = cms.bool(False),
0141    SaveFirstLevelSecondary = cms.untracked.bool(True),
0142    SavePrimaryDecayProductsAndConversionsInTracker = cms.untracked.bool(True),
0143    SavePrimaryDecayProductsAndConversionsInCalo    = cms.untracked.bool(True),
0144    SavePrimaryDecayProductsAndConversionsInMuon    = cms.untracked.bool(True)
0145 )
0146
0147 process.g4SimHits.SteppingAction = cms.PSet(
0148    process.common_maximum_timex, # need to be localy redefined
0149    KillBeamPipe            = cms.bool(True),
0150    CriticalEnergyForVacuum = cms.double(2.0),
0151    CriticalDensity         = cms.double(1e-15),
0152    EkinNames               = cms.vstring(),
0153    EkinThresholds          = cms.vdouble(),
0154    EkinParticles           = cms.vstring(),
0155    Verbosity               = cms.untracked.int32(0)
0156 )
0157
0158 process.g4SimHits.Watchers = cms.VPSet(cms.PSet(
0159     type = cms.string('CastorTestAnalysis'),
0160     CastorTestAnalysis = cms.PSet(
0161         EventNtupleFileName = cms.string('eventNtuple_muon.root'),
0162         Verbosity = cms.int32(0),
0163         StepNtupleFileName = cms.string('stepNtuple_electron.root'),
0164         StepNtupleFlag = cms.int32(0),
0165         EventNtupleFlag = cms.int32(1)
0166     )
0167 ))
0168
0169