test/python/testHitRelabeller_cfg.py

0001 ###############################################################################
0002 # Way to use this:
0003 #   cmsRun testHitRelabeller_cfg.py geometry=2021
0004 #
0005 #   Options for geometry 2016, 2017, 2018, 2021, 2023
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                  "2021",
0017                   VarParsing.VarParsing.multiplicity.singleton,
0018                   VarParsing.VarParsing.varType.string,
0019                   "geometry of operations: 2016, 2017, 2018, 2021, 2023")
0020
0021 ### get and parse the command line arguments
0022 options.parseArguments()
0023
0024 print(options)
0025
0026 ####################################################################
0027 # Use the options
0028
0029 geomName = "Configuration.Geometry.GeometryExtended" + options.geometry + "Reco_cff"
0030
0031 if (options.geometry == "2016"):
0032     from Configuration.Eras.Era_Run2_2016_cff import Run2_2016
0033     process = cms.Process('HitRelabeller',Run2_2016)
0034     baseName = "auto:run2_data"
0035 elif (options.geometry == "2017"):
0036     from Configuration.Eras.Era_Run2_2017_cff import Run2_2017
0037     process = cms.Process('HitRelabeller',Run2_2017)
0038     baseName = "auto:run2_data"
0039 elif (options.geometry == "2018"):
0040     from Configuration.Eras.Era_Run2_2018_cff import Run2_2018
0041     process = cms.Process('HitRelabeller',Run2_2018)
0042     baseName = "auto:run2_data"
0043 else:
0044     from Configuration.Eras.Era_Run3_DDD_cff import Run3_DDD
0045     process = cms.Process('HitRelabeller',Run3_DDD)
0046     baseName = "auto:run3_data"
0047
0048 print("Base file Name: ", baseName)
0049 print("Geom file Name: ", geomName)
0050
0051
0052 # import of standard configurations
0053 process.load('Configuration.StandardSequences.Services_cff')
0054 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0055 process.load('FWCore.MessageService.MessageLogger_cfi')
0056 process.load('Configuration.EventContent.EventContent_cff')
0057 process.load('SimGeneral.MixingModule.mixNoPU_cfi')
0058 process.load('Configuration.StandardSequences.MagneticField_cff')
0059 process.load('Configuration.StandardSequences.Generator_cff')
0060 process.load('IOMC.EventVertexGenerators.VtxSmearedRealistic25ns13p6TeVEarly2022Collision_cfi')
0061 process.load('GeneratorInterface.Core.genFilterSummary_cff')
0062 process.load('Configuration.StandardSequences.SimIdeal_cff')
0063 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0064 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0065 process.load(geomName)
0066
0067 if hasattr(process,'MessageLogger'):
0068     process.MessageLogger.HCalGeom=dict()
0069
0070 process.maxEvents = cms.untracked.PSet(
0071     input = cms.untracked.int32(1),
0072     output = cms.optional.untracked.allowed(cms.int32,cms.PSet)
0073 )
0074
0075 # Input source
0076 process.source = cms.Source("EmptySource")
0077
0078 process.options = cms.untracked.PSet(
0079     IgnoreCompletely = cms.untracked.vstring(),
0080     Rethrow = cms.untracked.vstring(),
0081     TryToContinue = cms.untracked.vstring(),
0082     allowUnscheduled = cms.obsolete.untracked.bool,
0083     canDeleteEarly = cms.untracked.vstring(),
0084     emptyRunLumiMode = cms.obsolete.untracked.string,
0085     eventSetup = cms.untracked.PSet(
0086         forceNumberOfConcurrentIOVs = cms.untracked.PSet(
0087
0088         ),
0089         numberOfConcurrentIOVs = cms.untracked.uint32(1)
0090     ),
0091     fileMode = cms.untracked.string('FULLMERGE'),
0092     forceEventSetupCacheClearOnNewRun = cms.untracked.bool(False),
0093     makeTriggerResults = cms.obsolete.untracked.bool,
0094     numberOfConcurrentLuminosityBlocks = cms.untracked.uint32(1),
0095     numberOfConcurrentRuns = cms.untracked.uint32(1),
0096     numberOfStreams = cms.untracked.uint32(0),
0097     numberOfThreads = cms.untracked.uint32(1),
0098     printDependencies = cms.untracked.bool(False),
0099     sizeOfStackForThreadsInKB = cms.optional.untracked.uint32,
0100     throwIfIllegalParameter = cms.untracked.bool(True),
0101     wantSummary = cms.untracked.bool(False)
0102 )
0103
0104 # Production Info
0105 process.configurationMetadata = cms.untracked.PSet(
0106     annotation = cms.untracked.string('ZMM_14TeV_TuneCP5_cfi nevts:10'),
0107     name = cms.untracked.string('Applications'),
0108     version = cms.untracked.string('$Revision: 1.19 $')
0109 )
0110
0111 # Other statements
0112 process.genstepfilter.triggerConditions=cms.vstring("generation_step")
0113 from Configuration.AlCa.GlobalTag import GlobalTag
0114 process.GlobalTag = GlobalTag(process.GlobalTag, baseName, '')
0115
0116 process.generator = cms.EDFilter("Pythia8GeneratorFilter",
0117                                  pythiaHepMCVerbosity = cms.untracked.bool(False),
0118                                  maxEventsToPrint = cms.untracked.int32(0),
0119                                  pythiaPylistVerbosity = cms.untracked.int32(0),
0120                                  filterEfficiency = cms.untracked.double(1.0),
0121                                  comEnergy = cms.double(14000.0),
0122                                  PythiaParameters = cms.PSet(
0123                                      pythia8CommonSettings = cms.vstring(
0124                                          'Tune:preferLHAPDF = 2',
0125                                          'Main:timesAllowErrors = 10000',
0126                                          'Check:epTolErr = 0.01',
0127                                          'Beams:setProductionScalesFromLHEF = off',
0128                                          'SLHA:minMassSM = 1000.',
0129                                          'ParticleDecays:limitTau0 = on',
0130                                          'ParticleDecays:tau0Max = 10',
0131                                         'ParticleDecays:allowPhotonRadiation = on',
0132                                      ),
0133                                      pythia8CP5Settings = cms.vstring(
0134                                          'Tune:pp 14',
0135                                          'Tune:ee 7',
0136                                          'MultipartonInteractions:ecmPow=0.03344',
0137                                          'MultipartonInteractions:bProfile=2',
0138                                          'MultipartonInteractions:pT0Ref=1.41',
0139                                          'MultipartonInteractions:coreRadius=0.7634',
0140                                          'MultipartonInteractions:coreFraction=0.63',
0141                                          'ColourReconnection:range=5.176',
0142                                          'SigmaTotal:zeroAXB=off',
0143                                          'SpaceShower:alphaSorder=2',
0144                                          'SpaceShower:alphaSvalue=0.118',
0145                                          'SigmaProcess:alphaSvalue=0.118',
0146                                          'SigmaProcess:alphaSorder=2',
0147                                          'MultipartonInteractions:alphaSvalue=0.118',
0148                                          'MultipartonInteractions:alphaSorder=2',
0149                                          'TimeShower:alphaSorder=2',
0150                                          'TimeShower:alphaSvalue=0.118',
0151                                          'SigmaTotal:mode = 0',
0152                                          'SigmaTotal:sigmaEl = 21.89',
0153                                          'SigmaTotal:sigmaTot = 100.309',
0154                                          'PDF:pSet=LHAPDF6:NNPDF31_nnlo_as_0118',
0155                                      ),
0156                                      processParameters = cms.vstring(
0157                                          'WeakSingleBoson:ffbar2gmZ = on',
0158                                          '23:onMode = off',
0159                                          '23:onIfAny = 13',
0160                                          'PhaseSpace:mHatMin = 75.',
0161                                      ),
0162                                      parameterSets = cms.vstring('pythia8CommonSettings',
0163                                                                  'pythia8CP5Settings',
0164                                                                  'processParameters',
0165                                                              )
0166                                  )
0167 )
0168
0169 process.mumugenfilter = cms.EDFilter("MCParticlePairFilter",
0170     MaxEta = cms.untracked.vdouble(4.0, 4.0),
0171     MinEta = cms.untracked.vdouble(-4.0, -4.0),
0172     MinPt = cms.untracked.vdouble(2.5, 2.5),
0173     ParticleCharge = cms.untracked.int32(-1),
0174     ParticleID1 = cms.untracked.vint32(13),
0175     ParticleID2 = cms.untracked.vint32(13),
0176     Status = cms.untracked.vint32(1, 1)
0177 )
0178
0179
0180 process.ProductionFilterSequence = cms.Sequence(process.generator+process.mumugenfilter)
0181 process.load("Geometry.HcalCommonData.hcalHitRelabellerTester_cfi")
0182
0183 # Path and EndPath definitions
0184 process.generation_step = cms.Path(process.pgen)
0185 process.simulation_step = cms.Path(process.psim)
0186 process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
0187 process.endjob_step = cms.EndPath(process.endOfProcess)
0188 process.analysis_step = cms.EndPath(process.hcalHitRelabellerTester)
0189
0190 # Schedule definition
0191 process.schedule = cms.Schedule(process.generation_step,
0192                                 process.genfiltersummary_step,
0193                                 process.simulation_step,
0194                                 process.endjob_step,
0195                                 process.analysis_step)
0196 from PhysicsTools.PatAlgos.tools.helpers import associatePatAlgosToolsTask
0197 associatePatAlgosToolsTask(process)
0198 # filter all path with the production filter sequence
0199 for path in process.paths:
0200     getattr(process,path).insert(0, process.ProductionFilterSequence)
0201
0202
0203 # Customisation from command line
0204
0205 # Add early deletion of temporary data products to reduce peak memory need
0206 from Configuration.StandardSequences.earlyDeleteSettings_cff import customiseEarlyDelete
0207 process = customiseEarlyDelete(process)
0208 # End adding early deletion