CMS 3D CMS Logo
MuonSeededStep_cff.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
2 
3 
4 from RecoMuon.MuonIdentification.earlyMuons_cfi import earlyMuons
5 
6 
7 import RecoTracker.SpecialSeedGenerators.outInSeedsFromStandaloneMuons_cfi
8 import RecoTracker.SpecialSeedGenerators.inOutSeedsFromTrackerMuons_cfi
9 muonSeededSeedsOutIn = RecoTracker.SpecialSeedGenerators.outInSeedsFromStandaloneMuons_cfi.outInSeedsFromStandaloneMuons.clone(
10  src = "earlyMuons",
11 )
12 muonSeededSeedsInOut = RecoTracker.SpecialSeedGenerators.inOutSeedsFromTrackerMuons_cfi.inOutSeedsFromTrackerMuons.clone(
13  src = "earlyMuons",
14 )
15 
16 from RecoTracker.SpecialSeedGenerators.outInSeedsFromStandaloneMuons_cfi import hitCollectorForOutInMuonSeeds
17 from Configuration.Eras.Modifier_tracker_apv_vfp30_2016_cff import tracker_apv_vfp30_2016 as _tracker_apv_vfp30_2016
18 _tracker_apv_vfp30_2016.toModify(hitCollectorForOutInMuonSeeds, MinPtForHitRecoveryInGluedDet=1e9)
19 
20 
22 import TrackingTools.TrajectoryCleaning.TrajectoryCleanerBySharedHits_cfi
23 muonSeededTrajectoryCleanerBySharedHits = TrackingTools.TrajectoryCleaning.TrajectoryCleanerBySharedHits_cfi.trajectoryCleanerBySharedHits.clone(
24  ComponentName = cms.string('muonSeededTrajectoryCleanerBySharedHits'),
25  fractionShared = cms.double(0.1),
26  ValidHitBonus = cms.double(1000.0),
27  MissingHitPenalty = cms.double(1.0),
28  ComponentType = cms.string('TrajectoryCleanerBySharedHits'),
29  allowSharedFirstHit = cms.bool(True)
30 )
31 
32 
33 
34 from TrackingTools.KalmanUpdators.Chi2MeasurementEstimator_cfi import Chi2MeasurementEstimator as _Chi2MeasurementEstimator
35 _muonSeededMeasurementEstimatorForInOutBase = _Chi2MeasurementEstimator.clone(
36  ComponentName = cms.string('muonSeededMeasurementEstimatorForInOut'),
37  MaxChi2 = cms.double(80.0),
38  nSigma = cms.double(4.),
39 )
40 muonSeededMeasurementEstimatorForInOut = _muonSeededMeasurementEstimatorForInOutBase.clone(
41  MaxSagitta = cms.double(-1.)
42 )
43 from Configuration.Eras.Modifier_trackingPhase2PU140_cff import trackingPhase2PU140
44 trackingPhase2PU140.toModify(muonSeededMeasurementEstimatorForInOut, MaxChi2 = 400.0, MaxSagitta = 2)
45 
46 _muonSeededMeasurementEstimatorForOutInBase = _Chi2MeasurementEstimator.clone(
47  ComponentName = cms.string('muonSeededMeasurementEstimatorForOutIn'),
48  MaxChi2 = cms.double(30.0),
49  nSigma = cms.double(3.),
50 )
51 from Configuration.Eras.Modifier_tracker_apv_vfp30_2016_cff import tracker_apv_vfp30_2016 as _tracker_apv_vfp30_2016
52 _tracker_apv_vfp30_2016.toModify(_muonSeededMeasurementEstimatorForOutInBase, MinPtForHitRecoveryInGluedDet=1e9)
53 muonSeededMeasurementEstimatorForOutIn = _muonSeededMeasurementEstimatorForOutInBase.clone(
54  MaxSagitta = cms.double(-1.)
55 )
56 
57 
58 import TrackingTools.TrajectoryFiltering.TrajectoryFilter_cff
59 muonSeededTrajectoryFilterForInOut = TrackingTools.TrajectoryFiltering.TrajectoryFilter_cff.CkfBaseTrajectoryFilter_block.clone()
60 muonSeededTrajectoryFilterForInOut.constantValueForLostHitsFractionFilter = 10
61 muonSeededTrajectoryFilterForInOut.minimumNumberOfHits = 3
62 
63 muonSeededTrajectoryFilterForOutIn = muonSeededTrajectoryFilterForInOut.clone()
64 muonSeededTrajectoryFilterForOutIn.constantValueForLostHitsFractionFilter = 10
65 muonSeededTrajectoryFilterForOutIn.minimumNumberOfHits = 5
66 
67 
68 import RecoTracker.CkfPattern.GroupedCkfTrajectoryBuilder_cfi
69 muonSeededTrajectoryBuilderForInOut = RecoTracker.CkfPattern.GroupedCkfTrajectoryBuilder_cfi.GroupedCkfTrajectoryBuilder.clone(
70  foundHitBonus = cms.double(1000.0),
71  lostHitPenalty = cms.double(1.0),
72  maxCand = cms.int32(5),
73  estimator = cms.string('muonSeededMeasurementEstimatorForInOut'),
74  trajectoryFilter = cms.PSet(refToPSet_ = cms.string('muonSeededTrajectoryFilterForInOut')),
75  inOutTrajectoryFilter = cms.PSet(refToPSet_ = cms.string('muonSeededTrajectoryFilterForInOut')), # not sure if it is used
76  minNrOfHitsForRebuild = cms.int32(2),
77  requireSeedHitsInRebuild = cms.bool(True),
78  keepOriginalIfRebuildFails = cms.bool(True),
79 )
80 muonSeededTrajectoryBuilderForOutIn = RecoTracker.CkfPattern.GroupedCkfTrajectoryBuilder_cfi.GroupedCkfTrajectoryBuilder.clone(
81  foundHitBonus = cms.double(1000.0),
82  lostHitPenalty = cms.double(1.0),
83  maxCand = cms.int32(3),
84  estimator = cms.string('muonSeededMeasurementEstimatorForOutIn'),
85  trajectoryFilter = cms.PSet(refToPSet_ = cms.string('muonSeededTrajectoryFilterForOutIn')),
86  inOutTrajectoryFilter = cms.PSet(refToPSet_ = cms.string('muonSeededTrajectoryFilterForOutIn')), # not sure if it is used
87  minNrOfHitsForRebuild = cms.int32(5),
88  requireSeedHitsInRebuild = cms.bool(True),
89  keepOriginalIfRebuildFails = cms.bool(False),
90 )
91 
92 
93 import TrackingTools.TrackFitters.RungeKuttaFitters_cff
94 muonSeededFittingSmootherWithOutliersRejectionAndRK = TrackingTools.TrackFitters.RungeKuttaFitters_cff.KFFittingSmootherWithOutliersRejectionAndRK.clone(
95  ComponentName = cms.string("muonSeededFittingSmootherWithOutliersRejectionAndRK"),
96  BreakTrajWith2ConsecutiveMissing = cms.bool(False),
97  EstimateCut = cms.double(50.),
98 )
99 
100 
101 import RecoTracker.CkfPattern.CkfTrackCandidates_cfi
102 muonSeededTrackCandidatesInOut = RecoTracker.CkfPattern.CkfTrackCandidates_cfi.ckfTrackCandidates.clone(
103  src = cms.InputTag("muonSeededSeedsInOut"),
104  TrajectoryBuilderPSet = cms.PSet(refToPSet_ = cms.string("muonSeededTrajectoryBuilderForInOut")),
105  TrajectoryCleaner = cms.string('muonSeededTrajectoryCleanerBySharedHits'),
106  RedundantSeedCleaner = cms.string("none"),
107 )
108 muonSeededTrackCandidatesOutIn = RecoTracker.CkfPattern.CkfTrackCandidates_cfi.ckfTrackCandidates.clone(
109  src = cms.InputTag("muonSeededSeedsOutIn"),
110  TrajectoryBuilderPSet = cms.PSet(refToPSet_ = cms.string("muonSeededTrajectoryBuilderForOutIn")),
111  TrajectoryCleaner = cms.string('muonSeededTrajectoryCleanerBySharedHits'),
112  numHitsForSeedCleaner = cms.int32(50),
113  onlyPixelHitsForSeedCleaner = cms.bool(False),
114 )
115 
116 
117 import RecoTracker.TrackProducer.TrackProducer_cfi
118 muonSeededTracksOutIn = RecoTracker.TrackProducer.TrackProducer_cfi.TrackProducer.clone(
119  src = cms.InputTag("muonSeededTrackCandidatesOutIn"),
120  AlgorithmName = cms.string('muonSeededStepOutIn'),
121  Fitter = cms.string("muonSeededFittingSmootherWithOutliersRejectionAndRK"),
122 )
123 muonSeededTracksInOut = RecoTracker.TrackProducer.TrackProducer_cfi.TrackProducer.clone(
124  src = cms.InputTag("muonSeededTrackCandidatesInOut"),
125  AlgorithmName = cms.string('muonSeededStepInOut'),
126  Fitter = cms.string("muonSeededFittingSmootherWithOutliersRejectionAndRK"),
127 )
128 
129 
130 # Final Classifier
131 from RecoTracker.FinalTrackSelectors.TrackCutClassifier_cff import *
132 muonSeededTracksInOutClassifier = TrackCutClassifier.clone()
133 muonSeededTracksInOutClassifier.src='muonSeededTracksInOut'
134 muonSeededTracksInOutClassifier.mva.minPixelHits = [0,0,0]
135 muonSeededTracksInOutClassifier.mva.maxChi2 = [9999.,9999.,9999.]
136 muonSeededTracksInOutClassifier.mva.maxChi2n = [10.0,1.0,0.4]
137 muonSeededTracksInOutClassifier.mva.minLayers = [3,5,5]
138 muonSeededTracksInOutClassifier.mva.min3DLayers = [1,2,2]
139 muonSeededTracksInOutClassifier.mva.maxLostLayers = [4,3,2]
140 
141 
142 muonSeededTracksOutInClassifier = TrackCutClassifier.clone()
143 muonSeededTracksOutInClassifier.src='muonSeededTracksOutIn'
144 muonSeededTracksOutInClassifier.mva.minPixelHits = [0,0,0]
145 muonSeededTracksOutInClassifier.mva.maxChi2 = [9999.,9999.,9999.]
146 muonSeededTracksOutInClassifier.mva.maxChi2n = [10.0,1.0,0.4]
147 muonSeededTracksOutInClassifier.mva.minLayers = [3,5,5]
148 muonSeededTracksOutInClassifier.mva.min3DLayers = [1,2,2]
149 muonSeededTracksOutInClassifier.mva.maxLostLayers = [4,3,2]
150 
151 
152 # For Phase2PU140
153 import RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi
154 muonSeededTracksInOutSelector = RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.multiTrackSelector.clone(
155  src='muonSeededTracksInOut',
156  trackSelectors= cms.VPSet(
157  RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.looseMTS.clone(
158  name = 'muonSeededTracksInOutLoose',
159  applyAdaptedPVCuts = cms.bool(False),
160  chi2n_par = 10.0,
161  minNumberLayers = 3,
162  min_nhits = 5,
163  maxNumberLostLayers = 4,
164  minNumber3DLayers = 0,
165  minHitsToBypassChecks = 7
166  ),
167  RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.tightMTS.clone(
168  name = 'muonSeededTracksInOutTight',
169  preFilterName = 'muonSeededTracksInOutLoose',
170  applyAdaptedPVCuts = cms.bool(False),
171  chi2n_par = 1.0,
172  minNumberLayers = 5,
173  min_nhits = 6,
174  maxNumberLostLayers = 3,
175  minNumber3DLayers = 2,
176  minHitsToBypassChecks = 10
177  ),
178  RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.highpurityMTS.clone(
179  name = 'muonSeededTracksInOutHighPurity',
180  preFilterName = 'muonSeededTracksInOutTight',
181  applyAdaptedPVCuts = cms.bool(False),
182  chi2n_par = 0.4,
183  minNumberLayers = 5,
184  min_nhits = 7,
185  maxNumberLostLayers = 2,
186  minNumber3DLayers = 2,
187  minHitsToBypassChecks = 20
188  ),
189  ) #end of vpset
190  ) #end of clone
191 muonSeededTracksOutInSelector = RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.multiTrackSelector.clone(
192  src='muonSeededTracksOutIn',
193  trackSelectors= cms.VPSet(
194  RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.looseMTS.clone(
195  name = 'muonSeededTracksOutInLoose',
196  applyAdaptedPVCuts = cms.bool(False),
197  chi2n_par = 10.0,
198  minNumberLayers = 3,
199  min_nhits = 5,
200  maxNumberLostLayers = 4,
201  minNumber3DLayers = 0,
202  minHitsToBypassChecks = 7
203  ),
204  RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.tightMTS.clone(
205  name = 'muonSeededTracksOutInTight',
206  preFilterName = 'muonSeededTracksOutInLoose',
207  applyAdaptedPVCuts = cms.bool(False),
208  chi2n_par = 1.0,
209  minNumberLayers = 5,
210  min_nhits = 6,
211  maxNumberLostLayers = 3,
212  minNumber3DLayers = 2,
213  minHitsToBypassChecks = 10
214  ),
215  RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.highpurityMTS.clone(
216  name = 'muonSeededTracksOutInHighPurity',
217  preFilterName = 'muonSeededTracksOutInTight',
218  applyAdaptedPVCuts = cms.bool(False),
219  chi2n_par = 0.4,
220  minNumberLayers = 5,
221  min_nhits = 7,
222  maxNumberLostLayers = 2,
223  minNumber3DLayers = 2,
224  minHitsToBypassChecks = 20
225  ),
226  ) #end of vpset
227  ) #end of clone
228 
229 
230 
231 muonSeededStepCoreInOutTask = cms.Task(
232  muonSeededSeedsInOut , muonSeededTrackCandidatesInOut , muonSeededTracksInOut
233 )
234 muonSeededStepCoreInOut = cms.Sequence(muonSeededStepCoreInOutTask)
235 
236 muonSeededStepCoreOutInTask = cms.Task(
237  muonSeededSeedsOutIn , muonSeededTrackCandidatesOutIn , muonSeededTracksOutIn
238 )
239 muonSeededStepCoreOutIn = cms.Sequence(muonSeededStepCoreOutInTask)
240 
241 muonSeededStepCoreTask = cms.Task(
242  muonSeededStepCoreInOutTask ,
243  muonSeededStepCoreOutInTask
244 )
245 muonSeededStepCore = cms.Sequence(muonSeededStepCoreTask)
246 #Phase2 : just muon Seed InOut is used in this moment
247 #trackingPhase2PU140.toReplaceWith(muonSeededStepCore, muonSeededStepCoreInOut)
248 muonSeededStepExtraInOutTask = cms.Task(
249  muonSeededTracksInOutClassifier
250 )
251 muonSeededStepExtraInOut = cms.Sequence(muonSeededStepExtraInOutTask)
252 
253 trackingPhase2PU140.toReplaceWith(muonSeededStepExtraInOutTask, cms.Task(
254  muonSeededTracksInOutSelector
255 ))
256 
257 muonSeededStepExtraTask = cms.Task(
258  muonSeededStepExtraInOutTask ,
259  muonSeededTracksOutInClassifier
260 )
261 
262 muonSeededStepExtra = cms.Sequence(muonSeededStepExtraTask)
263 trackingPhase2PU140.toReplaceWith(muonSeededStepExtraTask, cms.Task(
264  muonSeededStepExtraInOutTask ,
265  muonSeededTracksOutInSelector
266 ))
267 
268 muonSeededStepTask = cms.Task(
269  earlyMuons,
270  muonSeededStepCoreTask,
271  muonSeededStepExtraTask
272 )
273 muonSeededStep = cms.Sequence(muonSeededStepTask)
274 
275 
276 
277 muonSeededSeedsInOutAsTracks = cms.EDProducer("FakeTrackProducerFromSeed", src = cms.InputTag("muonSeededSeedsInOut"))
278 muonSeededSeedsOutInAsTracks = cms.EDProducer("FakeTrackProducerFromSeed", src = cms.InputTag("muonSeededSeedsOutIn"))
279 muonSeededTrackCandidatesInOutAsTracks = cms.EDProducer("FakeTrackProducerFromCandidate", src = cms.InputTag("muonSeededTrackCandidatesInOut"))
280 muonSeededTrackCandidatesOutInAsTracks = cms.EDProducer("FakeTrackProducerFromCandidate", src = cms.InputTag("muonSeededTrackCandidatesOutIn"))
281 muonSeededStepDebugInOutTask = cms.Task(
282  muonSeededSeedsInOutAsTracks , muonSeededTrackCandidatesInOutAsTracks
283 )
284 muonSeededStepDebugInOut = cms.Sequence(muonSeededStepDebugInOutTask)
285 muonSeededStepDebugTask = cms.Task(
286  muonSeededSeedsOutInAsTracks , muonSeededTrackCandidatesOutInAsTracks ,
287  muonSeededStepDebugInOutTask
288 )
289 muonSeededStepDebug = cms.Sequence(muonSeededStepDebugTask)
TrajectoryFilter_cff
GroupedCkfTrajectoryBuilder_cfi
TrackCutClassifier_cff
TrackProducer_cfi
RungeKuttaFitters_cff
TrajectoryCleanerBySharedHits_cfi
CkfTrackCandidates_cfi
Chi2MeasurementEstimator_cfi
Generated for CMSSW Reference Manual by doxygen 1.8.16