EcalDigis/test/SRValidationTest_cfg.py

0001 # The following comments couldn't be translated into the new config version:
0002
0003 #,
0004 import FWCore.ParameterSet.Config as cms
0005
0006 process = cms.Process("EcalSelectiveReadoutValid")
0007 # initialize  MessageLogger
0008 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0009
0010 # initialize magnetic field
0011 process.load("Configuration.StandardSequences.MagneticField_cff")
0012
0013 # geometry (Only Ecal)
0014 #process.load("Geometry.EcalCommonData.EcalOnly_cfi")
0015 # geometry (Only Ecal)
0016 process.load("Geometry.EcalCommonData.EcalOnly_cfi")
0017 process.load("Geometry.CaloEventSetup.CaloGeometry_cff")
0018 process.load("Geometry.CaloEventSetup.EcalTrigTowerConstituents_cfi")
0019 process.load("Geometry.EcalMapping.EcalMapping_cfi")
0020 process.load("Geometry.EcalMapping.EcalMappingRecord_cfi")
0021
0022
0023 # DQM services:
0024 process.load("DQMServices.Core.DQM_cfg")
0025
0026 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0027
0028 # event vertex smearing - applies only once (internal check)
0029 # Note : all internal generators will always do (0,0,0) vertex
0030 #
0031 process.load("IOMC.EventVertexGenerators.VtxSmearedGauss_cfi")
0032
0033 # run simulation, with EcalHits Validation specific watcher
0034 process.load("SimG4Core.Application.g4SimHits_cfi")
0035
0036 #  replace g4SimHits.Watchers = {
0037 #       { string type = "EcalSimHitsValidProducer"
0038 #         untracked string instanceLabel="EcalValidInfo"
0039 #         untracked bool verbose = false
0040 #       }
0041 #  }
0042 # Mixing Module
0043 process.load("SimGeneral.MixingModule.mixNoPU_cfi")
0044
0045 process.load("CalibCalorimetry.Configuration.Ecal_FakeConditions_cff")
0046
0047 # ECAL digitization sequence
0048 process.load("SimCalorimetry.Configuration.ecalDigiSequence_cff")
0049
0050 # ECAL digis validation sequence
0051 #include "Validation/EcalDigis/data/ecalDigisValidationSequence.cff"
0052 # Defines Ecal seletive readout validation module, ecalSelectiveReadoutValidation:
0053 process.load("Validation.EcalDigis.ecalSelectiveReadoutValidation_cfi")
0054 process.ecalSelectiveReadoutValidation.outputFile = 'srvalid_hists.root'
0055
0056 #ECAL reco sequence:
0057 process.load("RecoLocalCalo.Configuration.ecalLocalRecoSequence_cff")
0058 process.ecalWeightUncalibRecHit.EBdigiCollection = cms.InputTag("simEcalDigis", "ebDigis")
0059 process.ecalWeightUncalibRecHit.EEdigiCollection = cms.InputTag("simEcalDigis", "eeDigis")
0060
0061 process.maxEvents = cms.untracked.PSet(
0062     input = cms.untracked.int32(1)
0063 )
0064 process.o1 = cms.OutputModule("PoolOutputModule",
0065     outputCommands = cms.untracked.vstring('keep *'),
0066     fileName = cms.untracked.string('QCD_pt30_50_all_SRValidation.root')
0067 )
0068
0069 process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
0070     moduleSeeds = cms.PSet(
0071         g4SimHits = cms.untracked.uint32(9876),
0072         simEcalUnsuppressedDigis = cms.untracked.uint32(9876),
0073         VtxSmeared = cms.untracked.uint32(123456789)
0074     ),
0075     sourceSeed = cms.untracked.uint32(135799753)
0076 )
0077
0078 #Pythia configuration to generate multijet event with pt_hat between 30
0079 #and 50 GeV/c
0080 process.load("Configuration.Generator.QCD_Pt_30_50_cfi")
0081
0082 process.tpparams12 = cms.ESSource("EmptyESSource",
0083     recordName = cms.string('EcalTPGPhysicsConstRcd'),
0084     iovIsRunNotTime = cms.bool(True),
0085     firstValid = cms.vuint32(1)
0086 )
0087
0088 process.DQM.collectorHost = ''
0089
0090 process.g4SimHits.Generator.HepMCProductLabel = 'source'
0091
0092 process.simEcalDigis.writeSrFlags = True
0093
0094 # detector response simulation path:
0095 process.detSim = cms.Sequence(process.VtxSmeared*process.g4SimHits)
0096
0097 # processing path:
0098 process.p1 = cms.Path(process.detSim*process.mix*process.simEcalUnsuppressedDigis*process.simEcalTriggerPrimitiveDigis*process.simEcalDigis*process.ecalWeightUncalibRecHit*process.ecalRecHit*process.ecalSelectiveReadoutValidation)
0099
0100 process.outpath = cms.EndPath(process.o1)
0101
0102
0103 process.simEcalDigis.srpEndcapLowInterestChannelZS = -1.e9 #-0.06
0104 process.simEcalDigis.srpBarrelLowInterestChannelZS = -1.e9 #-0.035
0105