L1TZDC/test/l1ZDCProducerTest.py

0001 #Copied rom GM zdcDigiAnalyze.py 2023.08.04
0002 #Found here https://github.com/ginnocen/UPCopenHFanalysis/blob/zdc_calibrationcode/zdc_calibration/newZDCAnalyzer/test/zdcDigiAnalyze.py
0003 #CMcGinn it modifying to test the l1zdc producer, see comments below
0004 #Bugs, contact christopher.mc.ginn@cern.ch or cfmcginn on github
0005
0006 import FWCore.ParameterSet.Config as cms
0007 import FWCore.Utilities.FileUtils as FileUtils
0008
0009 from Configuration.StandardSequences.Eras import eras
0010
0011 process = cms.Process("Demo",eras.Run2_2018_pp_on_AA)
0012
0013 process.load('Configuration.StandardSequences.Services_cff')
0014 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0015 process.load('FWCore.MessageService.MessageLogger_cfi')
0016 process.load('Configuration.EventContent.EventContent_cff')
0017 process.load('SimGeneral.MixingModule.mixNoPU_cfi')
0018 process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
0019 process.load('Configuration.Geometry.GeometrySimDB_cff')
0020 process.load('Configuration.StandardSequences.MagneticField_38T_cff')
0021 process.load('Configuration.StandardSequences.Generator_cff')
0022 process.load('IOMC.EventVertexGenerators.VtxSmearedNominalCollision2015_cfi')
0023 process.load('GeneratorInterface.Core.genFilterSummary_cff')
0024 process.load('Configuration.StandardSequences.SimIdeal_cff')
0025 process.load('Configuration.StandardSequences.Digi_cff')
0026 process.load('Configuration.StandardSequences.SimL1Emulator_cff')
0027 process.load('Configuration.StandardSequences.DigiToRaw_cff')
0028 process.load('HLTrigger.Configuration.HLT_GRun_cff')
0029 process.load('Configuration.StandardSequences.RawToDigi_cff')
0030 process.load('Configuration.StandardSequences.Reconstruction_cff')
0031 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0032
0033 #process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_condDBv2_cff")
0034 #process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(-1))
0035 #process.MessageLogger.cerr.FwkReport.reportEvery = 10000
0036
0037 process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(1000))
0038 process.MessageLogger.cerr.FwkReport.reportEvery = 10
0039 #process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(10))
0040 #process.MessageLogger.cerr.FwkReport.reportEvery = 1
0041
0042 process.TFileService = cms.Service("TFileService",fileName=cms.string("zdcdigitree_DoDUMMYPyLUT_327524.root"))
0043
0044 #mylist = FileUtils.loadListFromFile('files_temp.txt')
0045 mylist = FileUtils.loadListFromFile('files_327524.txt')
0046 readFiles = cms.untracked.vstring(*mylist)
0047
0048 process.source = cms.Source("PoolSource",
0049                             fileNames = cms.untracked.vstring(
0050                             'file:/afs/cern.ch/user/m/mcsanad/public/CMSSW_10_3_1/src/zdc/newZDCAnalyzer/test/ED0B7A21-B558-924C-A57E-B1651E8BFFA3.root'
0051 #                           *mylist
0052                            )
0053 )
0054
0055 #Try some real basic replacement - producer and analyzer
0056 process.etSumZdcProducer = cms.EDProducer('L1TZDCProducer',
0057                                           zdcDigis = cms.InputTag("hcalDigis", "ZDC", "reRECO"),
0058                                           bxFirst = cms.int32(-2),
0059                                           bxLast = cms.int32(2),
0060                                           sampleToCenterBX = cms.int32(4)
0061
0062 )
0063
0064 process.etSumZdcAnalyzer = cms.EDAnalyzer('L1TZDCAnalyzer',
0065                                           etSumTag = cms.InputTag("etSumZdcProducer")
0066 )
0067
0068
0069 process.skimOutput = cms.OutputModule("PoolOutputModule",
0070                                       fileName = cms.untracked.string("comp.root"),
0071                                       outputCommands = cms.untracked.vstring(
0072                                           'keep *'
0073                                       ),
0074 )
0075
0076
0077 #process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0078 #process.GlobalTag.globaltag = '92X_upgrade2017_realistic_v10'
0079 #process.GlobalTag.globaltag = '103X_dataRun2_Express_v2'
0080 #process.GlobalTag.globaltag = '103X_dataRun2_Prompt_v3'
0081
0082 #CM Note - below global tag file doesnt appear to exist, swapping in the file from runEmulator-CaloStage2.py for now
0083 #process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_condDBv2_cff')
0084 #the below is the swap from runEmulator-CaloStage2.py
0085 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0086 process.GlobalTag.globaltag = cms.string('103X_dataRun2_v6')
0087
0088 #For lookup tables
0089 # Automatic addition of the customisation function from L1Trigger.Configuration.customiseSettings
0090 from L1Trigger.Configuration.customiseSettings import L1TSettingsToCaloParams_2018_v1_4_1
0091 #call to customisation function L1TSettingsToCaloParams_2018_v1_4_1 imported from L1Trigger.Configuration.customiseSettings
0092
0093 #process.caloStage2Params
0094 #L1TSettingsToCaloParams_2018_v1_4_1.caloStage2Params.zdcLUTFile        = cms.FileInPath("L1Trigger/L1TCalorimeter/data/zdcLUT_HI_DUMMY.txt"),
0095
0096 process = L1TSettingsToCaloParams_2018_v1_4_1(process)
0097 #process.caloStage2Params.zdcLUTFile        = cms.FileInPath("L1Trigger/L1TCalorimeter/data/zdcLUT_HI_DUMMY.txt"),
0098
0099
0100 process.produce_step = cms.Path(process.etSumZdcProducer)
0101 process.analyze_step = cms.Path(process.etSumZdcAnalyzer)
0102 process.output_step = cms.EndPath(process.skimOutput)
0103
0104
0105 process.schedule = cms.Schedule(
0106                     process.produce_step,
0107                     process.analyze_step,
0108                     process.output_step
0109 )
0110
0111
0112 #uncomment the output for debugging/saving the etsums