test/python/runHcalCellCountRun4_cfg.py

0001 ###############################################################################
0002 # Way to use this:
0003 #   cmsRun runHcalCellCountRun4_cfg.py geometry=D110
0004 #
0005 #   Options for geometry D95, D96, D98, D99, D100, D101, D102, D103, D104,
0006 #                        D105, D106, D107, D108, D109, D110, D111, D112, D113,
0007 #                        D114, D115, D116
0008 #
0009 ###############################################################################
0010 import FWCore.ParameterSet.Config as cms
0011 import os, sys, importlib, re
0012 import FWCore.ParameterSet.VarParsing as VarParsing
0013
0014 ####################################################################
0015 ### SETUP OPTIONS
0016 options = VarParsing.VarParsing('standard')
0017 options.register('geometry',
0018                  "D110",
0019                   VarParsing.VarParsing.multiplicity.singleton,
0020                   VarParsing.VarParsing.varType.string,
0021                   "geometry of operations: D95, D96, D98, D99, D100, D101, D102, D103, D104, D105, D106, D107, D108, D109, D110, D111, D112, D113, D114, D115, D116")
0022 ### get and parse the command line arguments
0023 options.parseArguments()
0024
0025 print(options)
0026
0027 ####################################################################
0028 # Use the options
0029
0030 geomFile = "Configuration.Geometry.GeometryExtendedRun4" + options.geometry + "Reco_cff"
0031 geomName = "Run4" + options.geometry
0032
0033 print("Geometry Name:  ", geomName)
0034 print("Geom file Name: ", geomFile)
0035
0036 import Configuration.Geometry.defaultPhase2ConditionsEra_cff as _settings
0037 GLOBAL_TAG, ERA = _settings.get_era_and_conditions(geomName)
0038
0039 print("Global Tag Name: ", GLOBAL_TAG)
0040 print("Era Name:        ", ERA)
0041
0042
0043 process = cms.Process('HcalCellCount',ERA)
0044
0045 process.load('Configuration.StandardSequences.Services_cff')
0046 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0047 process.load('FWCore.MessageService.MessageLogger_cfi')
0048 process.load('Configuration.EventContent.EventContent_cff')
0049 process.load(geomFile)
0050 process.load('Configuration.StandardSequences.MagneticField_cff')
0051 process.load('Configuration.StandardSequences.Generator_cff')
0052 process.load('IOMC.EventVertexGenerators.VtxSmearedRealisticHLLHC_cfi')
0053 process.load('GeneratorInterface.Core.genFilterSummary_cff')
0054 process.load('Configuration.StandardSequences.SimIdeal_cff')
0055 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0056 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0057 process.load('Geometry.HcalTowerAlgo.hcalCellCount_cfi')
0058
0059 process.MessageLogger.G4cout=dict()
0060
0061 process.maxEvents = cms.untracked.PSet(
0062     input = cms.untracked.int32(1)
0063 )
0064
0065 if hasattr(process,'MessageLogger'):
0066     process.MessageLogger.HCalGeom=dict()
0067
0068
0069 process.source = cms.Source("EmptySource")
0070
0071 process.generator = cms.EDProducer("FlatRandomPtGunProducer",
0072     PGunParameters = cms.PSet(
0073         PartID = cms.vint32(13),
0074         MinEta = cms.double(-2.5),
0075         MaxEta = cms.double(2.5),
0076         MinPhi = cms.double(-3.14159265359),
0077         MaxPhi = cms.double(3.14159265359),
0078         MinPt  = cms.double(9.99),
0079         MaxPt  = cms.double(10.01)
0080     ),
0081     AddAntiParticle = cms.bool(False),
0082     Verbosity       = cms.untracked.int32(0),
0083     firstRun        = cms.untracked.uint32(1)
0084 )
0085
0086 process.options = cms.untracked.PSet(
0087     IgnoreCompletely = cms.untracked.vstring(),
0088     Rethrow = cms.untracked.vstring(),
0089     TryToContinue = cms.untracked.vstring(),
0090     accelerators = cms.untracked.vstring('*'),
0091     allowUnscheduled = cms.obsolete.untracked.bool,
0092     canDeleteEarly = cms.untracked.vstring(),
0093     deleteNonConsumedUnscheduledModules = cms.untracked.bool(True),
0094     dumpOptions = cms.untracked.bool(False),
0095     emptyRunLumiMode = cms.obsolete.untracked.string,
0096     eventSetup = cms.untracked.PSet(
0097         forceNumberOfConcurrentIOVs = cms.untracked.PSet(
0098             allowAnyLabel_=cms.required.untracked.uint32
0099         ),
0100         numberOfConcurrentIOVs = cms.untracked.uint32(0)
0101     ),
0102     fileMode = cms.untracked.string('FULLMERGE'),
0103     forceEventSetupCacheClearOnNewRun = cms.untracked.bool(False),
0104     holdsReferencesToDeleteEarly = cms.untracked.VPSet(),
0105     makeTriggerResults = cms.obsolete.untracked.bool,
0106     modulesToCallForTryToContinue = cms.untracked.vstring(),
0107     modulesToIgnoreForDeleteEarly = cms.untracked.vstring(),
0108     numberOfConcurrentLuminosityBlocks = cms.untracked.uint32(0),
0109     numberOfConcurrentRuns = cms.untracked.uint32(1),
0110     numberOfStreams = cms.untracked.uint32(0),
0111     numberOfThreads = cms.untracked.uint32(1),
0112     printDependencies = cms.untracked.bool(False),
0113     sizeOfStackForThreadsInKB = cms.optional.untracked.uint32,
0114     throwIfIllegalParameter = cms.untracked.bool(True),
0115     wantSummary = cms.untracked.bool(False)
0116 )
0117
0118 process.ProductionFilterSequence = cms.Sequence(process.generator)
0119 from Configuration.AlCa.GlobalTag import GlobalTag
0120 process.GlobalTag = GlobalTag(process.GlobalTag, GLOBAL_TAG, '')
0121
0122 # Path and EndPath definitions
0123 process.generation_step = cms.Path(process.pgen)
0124 process.simulation_step = cms.Path(process.psim)
0125 process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
0126 process.endjob_step = cms.EndPath(process.endOfProcess)
0127 process.analysis_step = cms.EndPath(process.hcalCellCount)
0128
0129 # Schedule definition
0130 process.schedule = cms.Schedule(process.generation_step,
0131                                 process.genfiltersummary_step,
0132                                 process.simulation_step,
0133                                 process.analysis_step,
0134                                 process.endjob_step)
0135
0136 from PhysicsTools.PatAlgos.tools.helpers import associatePatAlgosToolsTask
0137 associatePatAlgosToolsTask(process)
0138 # filter all path with the production filter sequence
0139 for path in process.paths:
0140         getattr(process,path).insert(0, process.ProductionFilterSequence)
0141
0142 process.g4SimHits.UseMagneticField        = False
0143 process.g4SimHits.Physics.DefaultCutValue = 10.
0144