Skip to content

Commit c99cf57

Browse files
hqucmshqucms
committed
Update HGCal testbeam relval (wf 77.0).
1 parent 8226edd commit c99cf57

File tree

1 file changed

+88
-23
lines changed

1 file changed

+88
-23
lines changed

‎RecoLocalCalo/Configuration/python/hgcalTestBeamLocalReco_cff.py‎

Lines changed: 88 additions & 23 deletions
Original file line numberDiff line numberDiff line change
@@ -1,54 +1,116 @@
11
import FWCore.ParameterSet.Config as cms
2+
import os
3+
24

35
def runRecoForSep2024TB(process):
46
process.load('Configuration.StandardSequences.Accelerators_cff')
57

68
process.load(f"Configuration.Geometry.GeometryExtendedRun4D104Reco_cff")
79
process.load(f"Configuration.Geometry.GeometryExtendedRun4D104_cff")
810
from Geometry.HGCalMapping.hgcalmapping_cff import customise_hgcalmapper
9-
process = customise_hgcalmapper(process, modules='Geometry/HGCalMapping/data/ModuleMaps/modulelocator_Sep2024TBv2.txt')
11+
process = customise_hgcalmapper(
12+
process, modules='Geometry/HGCalMapping/data/ModuleMaps/modulelocator_P5v1.txt')
1013

11-
# Exclude rawMetaDataCollector from input TB file
12-
process.source.inputCommands = cms.untracked.vstring('keep *', 'drop *_rawMetaDataCollector_*_*')
14+
# Setup DTH source
15+
inputdir = '/store/group/dpg_hgcal/comm_hgcal/relval'
16+
files = [
17+
'run20000000_ls0000_EoR_source2000.jsn',
18+
'run20000000_ls0000_EoR_source2001.jsn',
19+
'run20000000_ls0001_EoLS_source2000.jsn',
20+
'run20000000_ls0001_EoLS_source2001.jsn',
21+
'run20000000_ls0001_index000000_source2000.raw',
22+
'run20000000_ls0001_index000000_source2001.raw',
23+
'config_feds_v1.json',
24+
'config_econds_v1.json',
25+
]
26+
localdir = 'local_daq/ramdisk/run20000000'
27+
os.makedirs(localdir, exist_ok=True)
28+
for f in files:
29+
os.system(f'xrdcp --silent -f root://cms-xrd-global.cern.ch//{inputdir}/{f} {localdir}')
1330

14-
# Keep HGCal output
15-
process.RecoLocalCaloRECO.outputCommands = ['keep *_hgcal*_*_*']
31+
process.EvFDaqDirector = cms.Service("EvFDaqDirector",
32+
baseDir=cms.untracked.string('local_daq/fu'),
33+
buBaseDir=cms.untracked.string('local_daq/ramdisk'),
34+
buBaseDirsAll=cms.untracked.vstring('local_daq/ramdisk'),
35+
buBaseDirsNumStreams=cms.untracked.vint32(2),
36+
buBaseDirsStreamIDs=cms.untracked.vint32(2000, 2001),
37+
directorIsBU=cms.untracked.bool(False),
38+
fileBrokerHost=cms.untracked.string('htcp40.cern.ch'),
39+
fileBrokerHostFromCfg=cms.untracked.bool(False),
40+
runNumber=cms.untracked.uint32(20000000),
41+
sourceIdentifier=cms.untracked.string('source'),
42+
useFileBroker=cms.untracked.bool(True)
43+
)
44+
45+
process.source = cms.Source("DAQSource",
46+
dataMode=cms.untracked.string('DTH'),
47+
eventChunkBlock=cms.untracked.uint32(100),
48+
eventChunkSize=cms.untracked.uint32(100),
49+
fileDiscoveryMode=cms.untracked.bool(True),
50+
keepRawFiles=cms.untracked.bool(True),
51+
maxBufferedFiles=cms.untracked.uint32(2),
52+
maxChunkSize=cms.untracked.uint32(1000),
53+
numBuffers=cms.untracked.uint32(3),
54+
testing=cms.untracked.bool(True),
55+
useL1EventID=cms.untracked.bool(False),
56+
verifyChecksum=cms.untracked.bool(False)
57+
)
58+
59+
process.FastMonitoringService = cms.Service("FastMonitoringService",
60+
sleepTime=cms.untracked.int32(1)
61+
)
62+
63+
process.hgcalConfigESProducer = cms.ESSource(
64+
"HGCalConfigurationESProducer", bePassthroughMode=cms.int32(-1),
65+
cbHeaderMarker=cms.int32(-1),
66+
charMode=cms.int32(-1),
67+
econdHeaderMarker=cms.int32(-1),
68+
fedjson=cms.string(os.path.abspath('local_daq/ramdisk/run20000000/config_feds_v1.json')),
69+
indexSource=cms.ESInputTag("hgCalMappingESProducer", ""),
70+
modjson=cms.string(os.path.abspath('local_daq/ramdisk/run20000000/config_econds_v1.json')),
71+
slinkHeaderMarker=cms.int32(-1))
72+
73+
# Setup HGCal unpacker
74+
process.hgcalDigis = cms.EDProducer("HGCalRawToDigi",
75+
src=cms.InputTag("rawDataCollector")
76+
)
1677

1778
# ESProducer to load calibration parameters from JSON file, like pedestals
18-
process.hgcalCalibParamESProducer = cms.ESProducer('hgcalrechit::HGCalCalibrationESProducer@alpaka',
19-
filename=cms.FileInPath('RecoLocalCalo/HGCalRecProducers/data/testbeam/level0_calib_Relay1727033054.json'),
79+
process.hgcalCalibParamESProducer = cms.ESProducer(
80+
'hgcalrechit::HGCalCalibrationESProducer@alpaka',
81+
filename=cms.FileInPath('RecoLocalCalo/HGCalRecProducers/data/testbeam/level0_calib_params_test.json'),
2082
filenameEnergyLoss=cms.FileInPath('RecoLocalCalo/HGCalRecProducers/data/testbeam/hgcal_energyloss_v16.json'),
21-
indexSource=cms.ESInputTag('hgCalMappingESProducer',''),
22-
mapSource=cms.ESInputTag('hgCalMappingModuleESProducer','')
23-
)
83+
indexSource=cms.ESInputTag('hgCalMappingESProducer', ''),
84+
mapSource=cms.ESInputTag('hgCalMappingModuleESProducer', ''))
2485

25-
process.hgcalSoARecHits = cms.EDProducer('alpaka_serial_sync::HGCalRecHitsProducer',
26-
digis=cms.InputTag('hgcalDigis', ''),
27-
calibSource=cms.ESInputTag('hgcalCalibParamESProducer', ''),
28-
n_hits_scale=cms.int32(1),
29-
n_blocks=cms.int32(1024),
30-
n_threads=cms.int32(4096)
31-
)
86+
process.hgcalSoARecHits = cms.EDProducer('HGCalRecHitsProducer@alpaka',
87+
digis=cms.InputTag('hgcalDigis', ''),
88+
calibSource=cms.ESInputTag('hgcalCalibParamESProducer', ''),
89+
n_hits_scale=cms.int32(1),
90+
n_blocks=cms.int32(1024),
91+
n_threads=cms.int32(1024)
92+
)
3293

3394
from RecoLocalCalo.HGCalRecProducers.hgCalSoARecHitsLayerClustersProducer_cfi import hgCalSoARecHitsLayerClustersProducer
3495
process.hgcalSoARecHitsLayerClusters = hgCalSoARecHitsLayerClustersProducer.clone(
35-
hgcalRecHitsSoA = "hgcalSoARecHits"
96+
hgcalRecHitsSoA="hgcalSoARecHits"
3697
)
3798

3899
from RecoLocalCalo.HGCalRecProducers.hgCalSoALayerClustersProducer_cfi import hgCalSoALayerClustersProducer
39100
process.hgcalSoALayerClusters = hgCalSoALayerClustersProducer.clone(
40-
hgcalRecHitsLayerClustersSoA = "hgcalSoARecHitsLayerClusters",
41-
hgcalRecHitsSoA = "hgcalSoARecHits"
101+
hgcalRecHitsLayerClustersSoA="hgcalSoARecHitsLayerClusters",
102+
hgcalRecHitsSoA="hgcalSoARecHits"
42103
)
43104

44105
from RecoLocalCalo.HGCalRecProducers.hgCalLayerClustersFromSoAProducer_cfi import hgCalLayerClustersFromSoAProducer
45106
process.hgcalMergeLayerClusters = hgCalLayerClustersFromSoAProducer.clone(
46-
hgcalRecHitsLayerClustersSoA = "hgcalSoARecHitsLayerClusters",
47-
hgcalRecHitsSoA = "hgcalSoARecHits",
48-
src = "hgcalSoALayerClusters"
107+
hgcalRecHitsLayerClustersSoA="hgcalSoARecHitsLayerClusters",
108+
hgcalRecHitsSoA="hgcalSoARecHits",
109+
src="hgcalSoALayerClusters"
49110
)
50111

51112
process.reco_task = cms.Task(
113+
process.hgcalDigis,
52114
process.hgcalSoARecHits,
53115
process.hgcalSoARecHitsLayerClusters,
54116
process.hgcalSoALayerClusters,
@@ -57,4 +119,7 @@ def runRecoForSep2024TB(process):
57119
process.hgcalTestBeamLocalRecoSequence = cms.Path(process.reco_task)
58120
process.schedule.insert(0, process.hgcalTestBeamLocalRecoSequence)
59121

122+
# Keep HGCal output
123+
process.FEVTDEBUGoutput.outputCommands = ['keep *_hgcal*_*_*']
124+
60125
return process

0 commit comments

Comments
 (0)