-
Notifications
You must be signed in to change notification settings - Fork 4.2k
/
customiseHybrid.py
44 lines (37 loc) · 2.13 KB
/
customiseHybrid.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
import FWCore.ParameterSet.Config as cms
## Adjust the local reco sequence for running on hybrid zero-suppressed data
def runOnHybridZS(process):
process.load("RecoLocalTracker.SiStripZeroSuppression.SiStripZeroSuppression_cfi")
process.load("RecoLocalTracker.SiStripClusterizer.SiStripClusterizer_cfi")
process.siStripZeroSuppression.Algorithms.APVInspectMode = "Hybrid"
zsInputs = process.siStripZeroSuppression.RawDigiProducersList
clusInputs = process.siStripClusters.DigiProducersList
unpackedZS = cms.InputTag("siStripDigis", "ZeroSuppressed")
zsInputs.append(unpackedZS)
clusInputs.remove(unpackedZS)
clusInputs.append(cms.InputTag("siStripZeroSuppression","ZeroSuppressed"))
# for on-demand clusterizer
from FWCore.ParameterSet.MassReplace import massReplaceParameter
massReplaceParameter(process, "HybridZeroSuppressed", cms.bool(False), cms.bool(True))
return process
## Change the (normal, ZS) repacker to use zero-suppressed hybrid data
def repackZSHybrid(process):
process.SiStripDigiToZSRaw.InputDigis = cms.InputTag("siStripZeroSuppression", "ZeroSuppressed")
process.DigiToRawRepack.insert(0, process.siStripZeroSuppression)
return process
## Add the ZS algorithm (in hybrid emulation mode) before repacking, to produce emulated hybrid samples with
## cmsDriver --step RAW2DIGI,REPACK:DigiToHybridRawRepack --customiseRecoLocalTracker/SiStripZeroSuppression/customiseHybrid.addHybridEmulationBeforeRepack ...
def addHybridEmulationBeforeRepack(process):
process.load("RecoLocalTracker.SiStripZeroSuppression.SiStripZeroSuppression_cfi")
zs = process.siStripZeroSuppression
zs.produceRawDigis = False
zs.produceHybridFormat = True
zs.Algorithms.APVInspectMode = "HybridEmulation"
zs.Algorithms.APVRestoreMode = ""
zs.Algorithms.CommonModeNoiseSubtractionMode = 'Median'
zs.Algorithms.MeanCM = 0
zs.Algorithms.DeltaCMThreshold = 20
zs.Algorithms.Use10bitsTruncation = True
zs.RawDigiProducersList = cms.VInputTag(cms.InputTag("siStripDigis", "VirginRaw"))
process.DigiToHybridRawRepack.insert(0, zs) ## insert before repacking
return process