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 #   connect = cms.string('frontier://cmsfrontier.cern.ch:8000/FrontierPrep/CMS_COND_30X_HCAL'),
0064     connect = cms.string('sqlite_file:testExample.db'),
0065      toGet = cms.VPSet(
0066          cms.PSet(
0067              record = cms.string('CastorPedestalsRcd'),
0068              tag = cms.string('castor_pedestals_v1.0_test')
0069              ),
0070          cms.PSet(
0071              record = cms.string('CastorPedestalWidthsRcd'),
0072              tag = cms.string('castor_widths_v1.0_test')
0073              ),
0074          cms.PSet(
0075              record = cms.string('CastorGainsRcd'),
0076              tag = cms.string('castor_gains_v1.0_test')
0077              ),
0078          cms.PSet(
0079              record = cms.string('CastorQIEDataRcd'),
0080              tag = cms.string('castor_qie_v1.0_test')
0081              ),
0082          cms.PSet(
0083              record = cms.string('CastorElectronicsMapRcd'),
0084              tag = cms.string('castor_emap_v1.0_test')
0085               ),
0086          cms.PSet(
0087              record = cms.string('CastorChannelQualityRcd'),
0088              tag = cms.string('castor_channelstatus_v1.0_test')
0089
0090              )
0091
0092      )
0093 )
0094 process.es_hardcode = cms.ESSource("CastorHardcodeCalibrations",
0095      toGet = cms.untracked.vstring('GainWidths')
0096  )
0097
0098 process.common_maximum_timex = cms.PSet( # need to be localy redefined
0099    MaxTrackTime  = cms.double(500.0),  # need to be localy redefined
0100    MaxTimeNames  = cms.vstring(), # need to be localy redefined
0101    MaxTrackTimes = cms.vdouble()  # need to be localy redefined
0102 )
0103
0104
0105
0106 process.CaloSD = cms.PSet(
0107     DetailedTiming = cms.bool(False),
0108     EminTrack = cms.double(1.0),
0109     Verbosity = cms.int32(0),
0110     UseMap = cms.bool(True),
0111     CheckHits = cms.int32(25)
0112 )
0113
0114
0115
0116
0117
0118 process.o1 = cms.OutputModule("PoolOutputModule",
0119     fileName = cms.untracked.string('gen-sim-digi-reco_muon.root')
0120 )
0121
0122 process.p1 = cms.Path(process.generator*process.pgen*process.VtxSmeared*process.g4SimHits*process.mix*process.simCastorDigis*process.castorreco)
0123 #*process.CastorFullReco)
0124 process.outpath=cms.EndPath(process.o1)
0125
0126 process.g4SimHits.UseMagneticField = False
0127 process.g4SimHits.Physics.DefaultCutValue = 10.
0128
0129 process.g4SimHits.Generator.ApplyEtaCuts = False
0130 process.g4SimHits.UseMagneticField = False
0131
0132 process.g4SimHits.CaloTrkProcessing.TestBeam = True
0133
0134 process.g4SimHits.StackingAction = cms.PSet(
0135    process.common_heavy_suppression,
0136    process.common_maximum_timex,        # need to be localy redefined
0137    TrackNeutrino = cms.bool(False),
0138    KillHeavy     = cms.bool(False),
0139    SaveFirstLevelSecondary = cms.untracked.bool(True),
0140    SavePrimaryDecayProductsAndConversionsInTracker = cms.untracked.bool(True),
0141    SavePrimaryDecayProductsAndConversionsInCalo    = cms.untracked.bool(True),
0142    SavePrimaryDecayProductsAndConversionsInMuon    = cms.untracked.bool(True)
0143 )
0144
0145 process.g4SimHits.SteppingAction = cms.PSet(
0146    process.common_maximum_timex, # need to be localy redefined
0147    KillBeamPipe            = cms.bool(True),
0148    CriticalEnergyForVacuum = cms.double(2.0),
0149    CriticalDensity         = cms.double(1e-15),
0150    EkinNames               = cms.vstring(),
0151    EkinThresholds          = cms.vdouble(),
0152    EkinParticles           = cms.vstring(),
0153    Verbosity               = cms.untracked.int32(0)
0154 )
0155
0156 process.g4SimHits.Watchers = cms.VPSet(cms.PSet(
0157     type = cms.string('CastorTestAnalysis'),
0158     CastorTestAnalysis = cms.PSet(
0159         EventNtupleFileName = cms.string('eventNtuple_muon.root'),
0160         Verbosity = cms.int32(0),
0161         StepNtupleFileName = cms.string('stepNtuple_electron.root'),
0162         StepNtupleFlag = cms.int32(0),
0163         EventNtupleFlag = cms.int32(1)
0164     )
0165 ))
0166
0167