test/python/runHGC11_cfg.py

0001 ###############################################################################
0002 # Way to use this:
0003 #   cmsRun runHGC1_cfg.py geometry=V18
0004 #   Options for geometry D107, D116, D104, D120
0005 #                        (V16) (V17) (V18) (V19)
0006 ###############################################################################
0007 import FWCore.ParameterSet.Config as cms
0008 import os, sys, importlib, re
0009 import FWCore.ParameterSet.VarParsing as VarParsing
0010
0011 ####################################################################
0012 ### SETUP OPTIONS
0013 options = VarParsing.VarParsing('standard')
0014 options.register('geometry',
0015                  "D116",
0016                   VarParsing.VarParsing.multiplicity.singleton,
0017                   VarParsing.VarParsing.varType.string,
0018                   "geometry of operations: D107, D116, D104, D120")
0019
0020 ### get and parse the command line arguments
0021
0022 options.parseArguments()
0023 print(options)
0024
0025 ####################################################################
0026
0027 geomName = "Run4" + options.geometry
0028 geomFile = "Configuration.Geometry.GeometryExtended" + geomName + "Reco_cff"
0029 import Configuration.Geometry.defaultPhase2ConditionsEra_cff as _settings
0030 GLOBAL_TAG, ERA = _settings.get_era_and_conditions(geomName)
0031 print("Geometry Name:   ", geomName)
0032 print("Geom file Name:  ", geomFile)
0033 print("Global Tag Name: ", GLOBAL_TAG)
0034 print("Era Name:        ", ERA)
0035
0036 process = cms.Process("PROD",ERA)
0037
0038 process.load('Configuration.StandardSequences.Services_cff')
0039 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0040 process.load('FWCore.MessageService.MessageLogger_cfi')
0041 process.load('Configuration.EventContent.EventContent_cff')
0042 process.load('SimGeneral.MixingModule.mixNoPU_cfi')
0043 process.load(geomFile)
0044 process.load("Configuration.StandardSequences.MagneticField_cff")
0045 process.load('Configuration.StandardSequences.Generator_cff')
0046 process.load('IOMC.EventVertexGenerators.VtxSmearedRealisticHLLHC_cfi')
0047 process.load('GeneratorInterface.Core.genFilterSummary_cff')
0048 process.load('Configuration.StandardSequences.SimIdeal_cff')
0049 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0050 from Configuration.AlCa.GlobalTag import GlobalTag
0051 process.GlobalTag = GlobalTag(process.GlobalTag, GLOBAL_TAG, '')
0052
0053 if hasattr(process,'MessageLogger'):
0054     process.MessageLogger.HGCalGeom = dict()
0055     process.MessageLogger.HGCSim = dict()
0056 #   process.MessageLogger.CaloSim = dict()
0057
0058 process.load("IOMC.RandomEngine.IOMC_cff")
0059 process.RandomNumberGeneratorService.generator.initialSeed = 456789
0060 process.RandomNumberGeneratorService.g4SimHits.initialSeed = 9876
0061 process.RandomNumberGeneratorService.VtxSmeared.initialSeed = 123456789
0062
0063 process.maxEvents = cms.untracked.PSet(
0064     input = cms.untracked.int32(100)
0065 )
0066
0067 process.source = cms.Source("EmptySource",
0068                             firstRun        = cms.untracked.uint32(1),
0069                             firstEvent      = cms.untracked.uint32(1)
0070 )
0071
0072 process.generator = cms.EDFilter("Pythia8PtGun",
0073                          PGunParameters = cms.PSet(
0074                              MaxPt = cms.double(100.01),
0075                              MinPt = cms.double(99.99),
0076                              ParticleID = cms.vint32(-11),
0077                              AddAntiParticle = cms.bool(True),
0078                              MaxEta = cms.double(3.1),
0079                              MaxPhi = cms.double(3.14159265359),
0080                              MinEta = cms.double(2.5),
0081                              MinPhi = cms.double(-3.14159265359) ## in radians
0082                          ),
0083                          Verbosity = cms.untracked.int32(0), ## set to 1 (or greater)  for printouts
0084                          psethack = cms.string('single mu pt 100'),
0085                          firstRun = cms.untracked.uint32(1),
0086                          PythiaParameters = cms.PSet(parameterSets = cms.vstring())
0087 )
0088
0089 # Other statements
0090 process.genstepfilter.triggerConditions=cms.vstring("generation_step")
0091 process.ProductionFilterSequence = cms.Sequence(process.generator)
0092
0093 process.generation_step = cms.Path(process.pgen)
0094 process.simulation_step = cms.Path(process.psim)
0095 process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
0096
0097 # Schedule definition
0098 process.schedule = cms.Schedule(process.generation_step,
0099                                 process.genfiltersummary_step,
0100                                 process.simulation_step,
0101                                 )
0102
0103 # filter all path with the production filter sequence
0104 for path in process.paths:
0105         getattr(process,path).insert(0, process.ProductionFilterSequence)