test/python/simHitValidation_cfg.py

0001 ###############################################################################
0002 # Way to use this:
0003 #   cmsRun simHitValidation_cfg.py geometry=D92
0004 #
0005 #   Options for geometry D88, D92, D93
0006 #
0007 ###############################################################################
0008 import FWCore.ParameterSet.Config as cms
0009 import os, sys, imp, re
0010 import FWCore.ParameterSet.VarParsing as VarParsing
0011
0012 ####################################################################
0013 ### SETUP OPTIONS
0014 options = VarParsing.VarParsing('standard')
0015 options.register('geometry',
0016                  "D88",
0017                   VarParsing.VarParsing.multiplicity.singleton,
0018                   VarParsing.VarParsing.varType.string,
0019                   "geometry of operations: D88, D92, D93")
0020
0021 ### get and parse the command line arguments
0022 options.parseArguments()
0023
0024 print(options)
0025
0026 ####################################################################
0027 # Use the options
0028 from Configuration.Eras.Era_Phase2C17I13M9_cff import Phase2C17I13M9
0029 process = cms.Process('HGCalSimHitLocal',Phase2C17I13M9)
0030
0031 geomFile = "Configuration.Geometry.GeometryExtended2026" + options.geometry + "Reco_cff"
0032 fileName = "file:SimHitVal" + options.geometry + ".root"
0033 outFile = "file:step1" + options.geometry + ".root"
0034
0035 print("Geometry file:          ", geomFile)
0036 print("SIM Output file:        ", outFile)
0037 print("Validation Output file: ", fileName)
0038
0039 # import of standard configurations
0040 process.load(geomFile)
0041 process.load('Configuration.StandardSequences.Services_cff')
0042 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0043 process.load('FWCore.MessageService.MessageLogger_cfi')
0044 process.load('Configuration.EventContent.EventContent_cff')
0045 process.load('FWCore.MessageService.MessageLogger_cfi')
0046 process.load('Configuration.EventContent.EventContent_cff')
0047 process.load('SimGeneral.MixingModule.mixNoPU_cfi')
0048 process.load('Configuration.StandardSequences.MagneticField_cff')
0049 process.load('Configuration.StandardSequences.Generator_cff')
0050 process.load('IOMC.EventVertexGenerators.VtxSmearedGauss_cfi')
0051 process.load('GeneratorInterface.Core.genFilterSummary_cff')
0052 process.load('Configuration.StandardSequences.SimIdeal_cff')
0053 process.load('Configuration.StandardSequences.Digi_cff')
0054 process.load('Configuration.StandardSequences.SimL1Emulator_cff')
0055 process.load('Configuration.StandardSequences.DigiToRaw_cff')
0056 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0057 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0058 from Configuration.AlCa.GlobalTag import GlobalTag
0059 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:phase2_realistic_T21', '')
0060
0061 process.maxEvents = cms.untracked.PSet(
0062         input = cms.untracked.int32(10)
0063 )
0064
0065 # Input source
0066 process.source = cms.Source("EmptySource")
0067
0068 process.options = cms.untracked.PSet()
0069
0070 # Production Info
0071 process.configurationMetadata = cms.untracked.PSet(
0072     version    = cms.untracked.string('$Revision: 1.20 $'),
0073     annotation = cms.untracked.string('SingleMuonPt10_cfi nevts:10'),
0074     name       = cms.untracked.string('Applications')
0075     )
0076
0077 # Output definition
0078
0079 process.FEVTDEBUGHLToutput = cms.OutputModule("PoolOutputModule",
0080                           splitLevel                   = cms.untracked.int32(0),
0081                           eventAutoFlushCompressedSize = cms.untracked.int32(5242880),
0082                           outputCommands               = process.FEVTDEBUGHLTEventContent.outputCommands,
0083                           fileName                     = cms.untracked.string(outFile),
0084                           dataset                      = cms.untracked.PSet(
0085     filterName = cms.untracked.string(''),
0086     dataTier   = cms.untracked.string('GEN-SIM-DIGI-RAW')
0087     ),
0088                           SelectEvents                 = cms.untracked.PSet(SelectEvents = cms.vstring('generation_step'))
0089                           )
0090
0091 # Additional output definition
0092
0093 # Other statements
0094 process.genstepfilter.triggerConditions=cms.vstring("generation_step")
0095
0096 process.generator = cms.EDProducer("FlatRandomPtGunProducer",
0097    PGunParameters = cms.PSet(
0098     MaxPt  = cms.double(10.01),
0099     MinPt  = cms.double(9.99),
0100     PartID = cms.vint32(13),
0101     MaxEta = cms.double(2.50),
0102     MaxPhi = cms.double(3.14159265359),
0103     MinEta = cms.double(1.75),
0104     MinPhi = cms.double(-3.14159265359)
0105     ),
0106    Verbosity       = cms.untracked.int32(0),
0107    psethack        = cms.string('single electron pt 10'),
0108    AddAntiParticle = cms.bool(True),
0109    firstRun        = cms.untracked.uint32(1)
0110 )
0111
0112 process.mix.digitizers = cms.PSet(process.theDigitizersValid)
0113
0114 # Path and EndPath definitions
0115 process.generation_step = cms.Path(process.pgen)
0116 process.simulation_step = cms.Path(process.psim)
0117 process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
0118 process.digitisation_step = cms.Path(process.pdigi_valid)
0119 process.L1simulation_step = cms.Path(process.SimL1Emulator)
0120 process.digi2raw_step = cms.Path(process.DigiToRaw)
0121 process.endjob_step = cms.EndPath(process.endOfProcess)
0122 process.FEVTDEBUGHLToutput_step = cms.EndPath(process.FEVTDEBUGHLToutput)
0123
0124 # Output definition
0125 process.ValidationOutput = cms.OutputModule("PoolOutputModule",
0126     outputCommands = cms.untracked.vstring('drop *', 'keep *_MEtoEDMConverter_*_*'),
0127                                             fileName = cms.untracked.string(fileName),
0128 )
0129
0130 process.load("DQMServices.Core.DQM_cfg")
0131 process.load("DQMServices.Components.MEtoEDMConverter_cfi")
0132
0133 process.load("Validation.HGCalValidation.digiValidation_cff")
0134
0135 ## path and endpath deffinition
0136 process.p1 = cms.Path(process.hgcalDigiValidationEE+
0137               process.hgcalDigiValidationHEF+
0138               process.hgcalDigiValidationHEB)
0139 process.p2 = cms.Path(process.MEtoEDMConverter)
0140 process.output_step = cms.EndPath(process.ValidationOutput)
0141
0142 # Schedule definition
0143 process.schedule = cms.Schedule(process.generation_step,process.genfiltersummary_step,
0144                 process.simulation_step,process.digitisation_step,
0145                 process.L1simulation_step,process.digi2raw_step,
0146                 process.endjob_step,process.p1,process.p2,process.output_step)
0147
0148 # filter all path with the production filter sequence
0149 for path in process.paths:
0150         getattr(process,path)._seq = process.generator * getattr(process,path)._seq
0151