test/python/testHGCalBHValid_cfg.py

0001 ###############################################################################
0002 # Way to use this:
0003 #   cmsRun runHGCalBHValid_cfg.py geometry=D88
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                  "D92",
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('HGCGeomAnalysis',Phase2C17I13M9)
0030
0031 geomFile = "Configuration.Geometry.GeometryExtended2026" + options.geometry + "Reco_cff"
0032 fileName = "hgcBHValid" + options.geometry + ".root"
0033
0034 print("Geometry file:  ", geomFile)
0035 print("Output file:    ", fileName)
0036
0037 # import of standard configurations
0038 process.load(geomFile)
0039 process.load('Configuration.StandardSequences.Services_cff')
0040 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0041 process.load('FWCore.MessageService.MessageLogger_cfi')
0042 process.load('Configuration.EventContent.EventContent_cff')
0043 process.load('SimGeneral.MixingModule.mixNoPU_cfi')
0044 process.load('Configuration.StandardSequences.MagneticField_cff')
0045 process.load('Configuration.StandardSequences.Generator_cff')
0046 process.load('IOMC.EventVertexGenerators.VtxSmearedRealistic50ns13TeVCollision_cfi')
0047 process.load('GeneratorInterface.Core.genFilterSummary_cff')
0048 process.load('Configuration.StandardSequences.SimIdeal_cff')
0049 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0050 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0051 process.load('Configuration.StandardSequences.Digi_cff')
0052 process.load('Configuration.StandardSequences.SimL1Emulator_cff')
0053 process.load('Configuration.StandardSequences.L1TrackTrigger_cff')
0054 process.load('Configuration.StandardSequences.DigiToRaw_cff')
0055 process.load('HLTrigger.Configuration.HLT_Fake2_cff')
0056 process.load('Configuration.StandardSequences.RawToDigi_cff')
0057 process.load('Configuration.StandardSequences.L1Reco_cff')
0058 process.load('Configuration.StandardSequences.Reconstruction_cff')
0059 process.load('Configuration.StandardSequences.RecoSim_cff')
0060 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0061
0062 process.maxEvents = cms.untracked.PSet(
0063     input = cms.untracked.int32(2000)
0064 )
0065
0066 process.MessageLogger.cerr.FwkReport.reportEvery = 5
0067 # Input source
0068 process.source = cms.Source("EmptySource")
0069
0070 process.options = cms.untracked.PSet(
0071     wantSummary = cms.untracked.bool(True)
0072 )
0073
0074 # Production Info
0075 process.configurationMetadata = cms.untracked.PSet(
0076     version = cms.untracked.string(''),
0077     annotation = cms.untracked.string(''),
0078     name = cms.untracked.string('Applications')
0079 )
0080
0081 if hasattr(process,'MessageLogger'):
0082     process.MessageLogger.ValidHGCal=dict()
0083     process.MessageLogger.HcalSim=dict()
0084
0085 # Additional output definition
0086 process.load('Validation.HGCalValidation.hgcalBHValidation_cfi')
0087
0088 process.TFileService = cms.Service("TFileService",
0089                                    fileName = cms.string(fileName),
0090                                    closeFileFast = cms.untracked.bool(True)
0091                                    )
0092
0093
0094 # Other statements
0095 process.genstepfilter.triggerConditions=cms.vstring("generation_step")
0096 from Configuration.AlCa.GlobalTag import GlobalTag
0097 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:phase2_realistic_T21', '')
0098
0099 process.generator = cms.EDProducer("FlatRandomPtGunProducer",
0100     PGunParameters = cms.PSet(
0101         MaxPt = cms.double(35.0),
0102         MinPt = cms.double(35.0),
0103         PartID = cms.vint32(13), #--->muon
0104         MinPhi = cms.double(-3.14159265359),
0105         MaxPhi = cms.double(3.14159265359),
0106         MinEta = cms.double(1.2),
0107         MaxEta = cms.double(3.0)
0108     ),
0109     Verbosity = cms.untracked.int32(0),
0110     psethack = cms.string('single muon pt 35'),
0111     AddAntiParticle = cms.bool(True),
0112     firstRun = cms.untracked.uint32(1)
0113 )
0114
0115
0116 #Modified to produce hgceedigis
0117 process.mix.digitizers = cms.PSet(process.theDigitizersValid)
0118 process.ProductionFilterSequence = cms.Sequence(process.generator)
0119
0120 #Following Removes Mag Field
0121 process.g4SimHits.UseMagneticField = False
0122 process.g4SimHits.Physics.bField = cms.double(0.0)
0123
0124
0125 # Path and EndPath definitions
0126 process.generation_step = cms.Path(process.pgen)
0127 process.simulation_step = cms.Path(process.psim)
0128 process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
0129 process.digitisation_step = cms.Path(process.pdigi_valid)
0130 process.L1simulation_step = cms.Path(process.SimL1Emulator)
0131 process.L1TrackTrigger_step = cms.Path(process.L1TrackTrigger)
0132 process.digi2raw_step = cms.Path(process.DigiToRaw)
0133 process.raw2digi_step = cms.Path(process.RawToDigi)
0134 process.L1Reco_step = cms.Path(process.L1Reco)
0135 process.reconstruction_step = cms.Path(process.localreco)
0136 process.recosim_step = cms.Path(process.recosim)
0137 process.analysis_step = cms.Path(process.hgcalBHAnalysis)
0138
0139 # Schedule definition
0140 process.schedule = cms.Schedule(process.generation_step,
0141                 process.simulation_step,
0142                 process.digitisation_step,
0143                                 process.L1simulation_step,
0144                                 process.L1TrackTrigger_step,
0145                                 process.digi2raw_step,
0146                                 process.raw2digi_step,
0147                                 process.L1Reco_step,
0148                                 process.reconstruction_step,
0149                                 process.recosim_step,
0150                                 process.analysis_step,
0151                 )
0152
0153 # filter all path with the production filter sequence
0154 for path in process.paths:
0155         if getattr(process,path)._seq is not None: getattr(process,path)._seq = process.ProductionFilterSequence * getattr(process,path)._seq