|
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
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
|
import FWCore.ParameterSet.Config as cms
import FWCore.ParameterSet.VarParsing as VarParsing
options = VarParsing.VarParsing()
options.register("isUnitTest",
False,
VarParsing.VarParsing.multiplicity.singleton,
VarParsing.VarParsing.varType.bool,
"are we running the unit test")
options.parseArguments()
process = cms.Process("HitEff")
process.load("Configuration.StandardSequences.MagneticField_cff")
process.load("Configuration.StandardSequences.GeometryRecoDB_cff")
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
process.load('FWCore.MessageService.MessageLogger_cfi')
from Configuration.AlCa.GlobalTag import GlobalTag
process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:run2_data', '')
process.source = cms.Source("PoolSource", fileNames=cms.untracked.vstring(
# 10 random files, will need the rest later
"/store/express/Run2018D/StreamExpress/ALCARECO/SiStripCalMinBias-Express-v1/000/325/172/00000/E3D6AECF-3F12-6540-97DC-4A6994CFEBF3.root",
"/store/express/Run2018D/StreamExpress/ALCARECO/SiStripCalMinBias-Express-v1/000/325/172/00000/242566C3-0540-8C43-8D6E-BB42C1FE0BB5.root",
"/store/express/Run2018D/StreamExpress/ALCARECO/SiStripCalMinBias-Express-v1/000/325/172/00000/BC8C0839-F645-B948-9040-15FCB5D50472.root",
"/store/express/Run2018D/StreamExpress/ALCARECO/SiStripCalMinBias-Express-v1/000/325/172/00000/3A806401-2CBC-4345-A5CB-593AABD4BE4E.root",
"/store/express/Run2018D/StreamExpress/ALCARECO/SiStripCalMinBias-Express-v1/000/325/172/00000/852C3C1E-2BD4-A843-A65B-51110A503FBD.root",
"/store/express/Run2018D/StreamExpress/ALCARECO/SiStripCalMinBias-Express-v1/000/325/172/00000/B795F9A0-4681-A34A-B879-E33A0DEC8720.root",
"/store/express/Run2018D/StreamExpress/ALCARECO/SiStripCalMinBias-Express-v1/000/325/172/00000/3A0884F2-A395-C541-8EFB-740C45A57CCE.root",
"/store/express/Run2018D/StreamExpress/ALCARECO/SiStripCalMinBias-Express-v1/000/325/172/00000/D274E7C1-5A9D-A544-B9B3-6A30166FC16C.root",
"/store/express/Run2018D/StreamExpress/ALCARECO/SiStripCalMinBias-Express-v1/000/325/172/00000/C4D243DC-2A09-CF42-A050-7678EF4A90D7.root",
"/store/express/Run2018D/StreamExpress/ALCARECO/SiStripCalMinBias-Express-v1/000/325/172/00000/7946A89D-8AC5-6B4F-BAD2-AE3B971865C5.root",
))
if(options.isUnitTest):
process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(20))
else:
process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(100000))
process.load("RecoVertex.BeamSpotProducer.BeamSpot_cfi")
process.load("RecoTracker.TrackProducer.TrackRefitters_cff")
process.refitTracks = process.TrackRefitter.clone(src=cms.InputTag("ALCARECOSiStripCalMinBias"))
tracks = cms.InputTag("refitTracks")
process.hiteff = cms.EDProducer("SiStripHitEfficiencyWorker",
lumiScalers =cms.InputTag("scalersRawToDigi"),
addLumi = cms.untracked.bool(True),
commonMode = cms.InputTag("siStripDigis", "CommonMode"),
addCommonMode = cms.untracked.bool(False),
combinatorialTracks = tracks,
trajectories = tracks,
siStripClusters = cms.InputTag("siStripClusters"),
siStripDigis = cms.InputTag("siStripDigis"),
trackerEvent = cms.InputTag("MeasurementTrackerEvent"),
# part 2
Layer = cms.int32(0), # =0 means do all layers
Debug = cms.untracked.bool(True),
# do not cut on the total number of tracks
cutOnTracks = cms.bool(True),
trackMultiplicity = cms.uint32(100),
# use or not first and last measurement of a trajectory (biases), default is false
useFirstMeas = cms.bool(False),
useLastMeas = cms.bool(False),
useAllHitsFromTracksWithMissingHits = cms.bool(False),
## non-default settings
ClusterMatchingMethod = cms.int32(4), # default 0 case0,1,2,3,4
ClusterTrajDist = cms.double(15), # default 64
)
process.load("DQM.SiStripCommon.TkHistoMap_cff")
## OLD HITEFF
from CalibTracker.SiStripHitEfficiency.SiStripHitEff_cff import anEff
process.anEff = anEff
process.anEff.Debug = True
process.anEff.combinatorialTracks = tracks
process.anEff.trajectories = tracks
process.TFileService = cms.Service("TFileService",
fileName = cms.string('HitEffTree.root')
)
process.load("CalibTracker.SiStripCommon.ShallowEventDataProducer_cfi")
process.eventInfo = cms.EDAnalyzer(
"ShallowTree",
outputCommands = cms.untracked.vstring(
'drop *',
'keep *_shallowEventRun_*_*',
)
)
## END OLD HITEFF
## TODO double-check in main CalibTree config if hiteff also takes refitted tracks
process.allPath = cms.Path(process.MeasurementTrackerEvent*process.offlineBeamSpot*process.refitTracks
*process.anEff*process.shallowEventRun*process.eventInfo
*process.hiteff)
# save the DQM plots in the DQMIO format
process.dqmOutput = cms.OutputModule("DQMRootOutputModule",
fileName = cms.untracked.string("DQM.root")
)
# also save in legacy format, for easier comparison
process.dqmSaver = cms.EDAnalyzer("DQMFileSaver",
convention = cms.untracked.string('Offline'),
fileFormat = cms.untracked.string('ROOT'),
producer = cms.untracked.string('DQM'),
workflow = cms.untracked.string('/Harvesting/SiStripHitEfficiency/All'),
dirName = cms.untracked.string('.'),
saveByRun = cms.untracked.int32(1),
saveAtJobEnd = cms.untracked.bool(False),
)
process.HitEffOutput = cms.EndPath(process.dqmOutput*process.dqmSaver)
if(options.isUnitTest):
process.MessageLogger.cerr.enable = False
process.MessageLogger.TkHistoMap = dict()
process.MessageLogger.SiStripHitEfficiency = dict()
process.MessageLogger.SiStripHitEfficiencyWorker = dict()
process.MessageLogger.cout = cms.untracked.PSet(
enable = cms.untracked.bool(True),
threshold = cms.untracked.string("INFO"),
enableStatistics = cms.untracked.bool(True),
default = cms.untracked.PSet(limit = cms.untracked.int32(0)),
FwkReport = cms.untracked.PSet(limit = cms.untracked.int32(-1),
reportEvery = cms.untracked.int32(1)),
TkHistoMap = cms.untracked.PSet( limit = cms.untracked.int32(-1)),
SiStripHitEfficiency = cms.untracked.PSet( limit = cms.untracked.int32(-1)),
SiStripHitEfficiencyWorker = cms.untracked.PSet( limit = cms.untracked.int32(-1))
)
else:
process.MessageLogger = cms.Service(
"MessageLogger",
destinations = cms.untracked.vstring("log_tkhistomap"),
debugModules = cms.untracked.vstring("hiteff", "anEff"),
log_tkhistomap = cms.untracked.PSet(threshold = cms.untracked.string("DEBUG"),
default = cms.untracked.PSet(limit = cms.untracked.int32(-1))),
categories=cms.untracked.vstring("TkHistoMap",
"SiStripHitEfficiency:HitEff",
"SiStripHitEfficiency",
"SiStripHitEfficiencyWorker")
)
# Run the rest of the CT-based sequence with
# cmsRun test/testSiStripHitEffFromCalibTree_cfg.py inputFiles=HitEffTree.root runNumber=325172
|