CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
RunI_MixedTripletStep_cff.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
2 
3 ###############################################################
4 # Large impact parameter Tracking using mixed-triplet seeding #
5 ###############################################################
6 from RecoLocalTracker.SubCollectionProducers.trackClusterRemover_cfi import *
7 mixedTripletStepClusters = trackClusterRemover.clone(
8  maxChi2 = cms.double(9.0),
9  trajectories = cms.InputTag("detachedTripletStepTracks"),
10  pixelClusters = cms.InputTag("siPixelClusters"),
11  stripClusters = cms.InputTag("siStripClusters"),
12  oldClusterRemovalInfo = cms.InputTag("detachedTripletStepClusters"),
13  overrideTrkQuals = cms.InputTag('detachedTripletStep'),
14  TrackQuality = cms.string('highPurity'),
15  minNumberOfLayersWithMeasBeforeFiltering = cms.int32(0),
16 )
17 
18 # SEEDING LAYERS
19 mixedTripletStepSeedLayersA = cms.EDProducer("SeedingLayersEDProducer",
20  layerList = cms.vstring('BPix1+BPix2+BPix3',
21  'BPix1+BPix2+FPix1_pos', 'BPix1+BPix2+FPix1_neg',
22  'BPix1+FPix1_pos+FPix2_pos', 'BPix1+FPix1_neg+FPix2_neg',
23  'BPix2+FPix1_pos+FPix2_pos', 'BPix2+FPix1_neg+FPix2_neg',
24  'FPix1_pos+FPix2_pos+TEC1_pos', 'FPix1_neg+FPix2_neg+TEC1_neg',
25  'FPix2_pos+TEC2_pos+TEC3_pos', 'FPix2_neg+TEC2_neg+TEC3_neg'),
26  BPix = cms.PSet(
27  TTRHBuilder = cms.string('WithTrackAngle'),
28  HitProducer = cms.string('siPixelRecHits'),
29  skipClusters = cms.InputTag('mixedTripletStepClusters')
30  ),
31  FPix = cms.PSet(
32  TTRHBuilder = cms.string('WithTrackAngle'),
33  HitProducer = cms.string('siPixelRecHits'),
34  skipClusters = cms.InputTag('mixedTripletStepClusters')
35  ),
36  TEC = cms.PSet(
37  matchedRecHits = cms.InputTag("siStripMatchedRecHits","matchedRecHit"),
38  useRingSlector = cms.bool(True),
39  TTRHBuilder = cms.string('WithTrackAngle'), minGoodCharge = cms.double(2069),
40  minRing = cms.int32(1),
41  maxRing = cms.int32(1),
42  skipClusters = cms.InputTag('mixedTripletStepClusters')
43  )
44 )
45 
46 # SEEDS
47 from RecoPixelVertexing.PixelTriplets.PixelTripletLargeTipGenerator_cfi import *
48 PixelTripletLargeTipGenerator.extraHitRZtolerance = 0.0
49 PixelTripletLargeTipGenerator.extraHitRPhitolerance = 0.0
50 import RecoTracker.TkSeedGenerator.GlobalSeedsFromTriplets_cff
51 mixedTripletStepSeedsA = RecoTracker.TkSeedGenerator.GlobalSeedsFromTriplets_cff.globalSeedsFromTriplets.clone()
52 mixedTripletStepSeedsA.OrderedHitsFactoryPSet.SeedingLayers = 'mixedTripletStepSeedLayersA'
53 mixedTripletStepSeedsA.OrderedHitsFactoryPSet.GeneratorPSet = cms.PSet(PixelTripletLargeTipGenerator)
54 mixedTripletStepSeedsA.SeedCreatorPSet.ComponentName = 'SeedFromConsecutiveHitsTripletOnlyCreator'
55 mixedTripletStepSeedsA.RegionFactoryPSet.RegionPSet.ptMin = 0.4
56 mixedTripletStepSeedsA.RegionFactoryPSet.RegionPSet.originHalfLength = 15.0
57 mixedTripletStepSeedsA.RegionFactoryPSet.RegionPSet.originRadius = 1.5
58 
59 mixedTripletStepSeedsA.SeedComparitorPSet = cms.PSet(
60  ComponentName = cms.string('PixelClusterShapeSeedComparitor'),
61  FilterAtHelixStage = cms.bool(False),
62  FilterPixelHits = cms.bool(True),
63  FilterStripHits = cms.bool(True),
64  ClusterShapeHitFilterName = cms.string('ClusterShapeHitFilter'),
65  ClusterShapeCacheSrc = cms.InputTag('siPixelClusterShapeCache')
66  )
67 
68 # SEEDING LAYERS
69 mixedTripletStepSeedLayersB = cms.EDProducer("SeedingLayersEDProducer",
70  layerList = cms.vstring('BPix2+BPix3+TIB1', 'BPix2+BPix3+TIB2'),
71  BPix = cms.PSet(
72  TTRHBuilder = cms.string('WithTrackAngle'),
73  HitProducer = cms.string('siPixelRecHits'),
74  skipClusters = cms.InputTag('mixedTripletStepClusters')
75  ),
76  TIB = cms.PSet(
77  matchedRecHits = cms.InputTag("siStripMatchedRecHits","matchedRecHit"),
78  TTRHBuilder = cms.string('WithTrackAngle'), minGoodCharge = cms.double(2069),
79  skipClusters = cms.InputTag('mixedTripletStepClusters')
80  )
81 )
82 
83 # SEEDS
84 from RecoPixelVertexing.PixelTriplets.PixelTripletLargeTipGenerator_cfi import *
85 PixelTripletLargeTipGenerator.extraHitRZtolerance = 0.0
86 PixelTripletLargeTipGenerator.extraHitRPhitolerance = 0.0
87 import RecoTracker.TkSeedGenerator.GlobalSeedsFromTriplets_cff
88 mixedTripletStepSeedsB = RecoTracker.TkSeedGenerator.GlobalSeedsFromTriplets_cff.globalSeedsFromTriplets.clone()
89 mixedTripletStepSeedsB.OrderedHitsFactoryPSet.SeedingLayers = 'mixedTripletStepSeedLayersB'
90 mixedTripletStepSeedsB.OrderedHitsFactoryPSet.GeneratorPSet = cms.PSet(PixelTripletLargeTipGenerator)
91 mixedTripletStepSeedsB.SeedCreatorPSet.ComponentName = 'SeedFromConsecutiveHitsTripletOnlyCreator'
92 mixedTripletStepSeedsB.RegionFactoryPSet.RegionPSet.ptMin = 0.6
93 mixedTripletStepSeedsB.RegionFactoryPSet.RegionPSet.originHalfLength = 10.0
94 mixedTripletStepSeedsB.RegionFactoryPSet.RegionPSet.originRadius = 1.5
95 
96 mixedTripletStepSeedsB.SeedComparitorPSet = cms.PSet(
97  ComponentName = cms.string('PixelClusterShapeSeedComparitor'),
98  FilterAtHelixStage = cms.bool(False),
99  FilterPixelHits = cms.bool(True),
100  FilterStripHits = cms.bool(True),
101  ClusterShapeHitFilterName = cms.string('ClusterShapeHitFilter')
102  )
103 
104 import RecoTracker.TkSeedGenerator.GlobalCombinedSeeds_cfi
105 mixedTripletStepSeeds = RecoTracker.TkSeedGenerator.GlobalCombinedSeeds_cfi.globalCombinedSeeds.clone()
106 mixedTripletStepSeeds.seedCollections = cms.VInputTag(
107  cms.InputTag('mixedTripletStepSeedsA'),
108  cms.InputTag('mixedTripletStepSeedsB'),
109  )
110 
111 # QUALITY CUTS DURING TRACK BUILDING
112 import TrackingTools.TrajectoryFiltering.TrajectoryFilter_cff
113 mixedTripletStepTrajectoryFilter = TrackingTools.TrajectoryFiltering.TrajectoryFilter_cff.CkfBaseTrajectoryFilter_block.clone(
114  maxLostHits = 0,
115  minimumNumberOfHits = 3,
116  minPt = 0.1
117  )
118 
119 # Propagator taking into account momentum uncertainty in multiple scattering calculation.
120 import TrackingTools.MaterialEffects.MaterialPropagator_cfi
121 mixedTripletStepPropagator = TrackingTools.MaterialEffects.MaterialPropagator_cfi.MaterialPropagator.clone(
122  ComponentName = 'mixedTripletStepPropagator',
123  ptMin = 0.1
124  )
125 import TrackingTools.MaterialEffects.OppositeMaterialPropagator_cfi
126 mixedTripletStepPropagatorOpposite = TrackingTools.MaterialEffects.OppositeMaterialPropagator_cfi.OppositeMaterialPropagator.clone(
127  ComponentName = 'mixedTripletStepPropagatorOpposite',
128  ptMin = 0.1
129  )
130 
131 import TrackingTools.KalmanUpdators.Chi2MeasurementEstimatorESProducer_cfi
132 mixedTripletStepChi2Est = TrackingTools.KalmanUpdators.Chi2MeasurementEstimatorESProducer_cfi.Chi2MeasurementEstimator.clone(
133  ComponentName = cms.string('mixedTripletStepChi2Est'),
134  nSigma = cms.double(3.0),
135  MaxChi2 = cms.double(16.0)
136 )
137 
138 # TRACK BUILDING
139 import RecoTracker.CkfPattern.GroupedCkfTrajectoryBuilder_cfi
140 mixedTripletStepTrajectoryBuilder = RecoTracker.CkfPattern.GroupedCkfTrajectoryBuilder_cfi.GroupedCkfTrajectoryBuilder.clone(
141  MeasurementTrackerName = '',
142  trajectoryFilter = cms.PSet(refToPSet_ = cms.string('mixedTripletStepTrajectoryFilter')),
143  propagatorAlong = cms.string('mixedTripletStepPropagator'),
144  propagatorOpposite = cms.string('mixedTripletStepPropagatorOpposite'),
145  maxCand = 2,
146  estimator = cms.string('mixedTripletStepChi2Est'),
147  maxDPhiForLooperReconstruction = cms.double(2.0),
148  maxPtForLooperReconstruction = cms.double(0.7)
149  )
150 
151 # MAKING OF TRACK CANDIDATES
152 import RecoTracker.CkfPattern.CkfTrackCandidates_cfi
153 mixedTripletStepTrackCandidates = RecoTracker.CkfPattern.CkfTrackCandidates_cfi.ckfTrackCandidates.clone(
154  src = cms.InputTag('mixedTripletStepSeeds'),
155  clustersToSkip = cms.InputTag('mixedTripletStepClusters'),
156  ### these two parameters are relevant only for the CachingSeedCleanerBySharedInput
157  numHitsForSeedCleaner = cms.int32(50),
158  #onlyPixelHitsForSeedCleaner = cms.bool(True),
159 
160  TrajectoryBuilderPSet = cms.PSet(refToPSet_ = cms.string('mixedTripletStepTrajectoryBuilder')),
161  doSeedingRegionRebuilding = True,
162  useHitsSplitting = True
163 )
164 
165 
166 from TrackingTools.TrajectoryCleaning.TrajectoryCleanerBySharedHits_cfi import trajectoryCleanerBySharedHits
167 mixedTripletStepTrajectoryCleanerBySharedHits = trajectoryCleanerBySharedHits.clone(
168  ComponentName = cms.string('mixedTripletStepTrajectoryCleanerBySharedHits'),
169  fractionShared = cms.double(0.11),
170  allowSharedFirstHit = cms.bool(True)
171  )
172 mixedTripletStepTrackCandidates.TrajectoryCleaner = 'mixedTripletStepTrajectoryCleanerBySharedHits'
173 
174 
175 # TRACK FITTING
176 import RecoTracker.TrackProducer.TrackProducer_cfi
177 mixedTripletStepTracks = RecoTracker.TrackProducer.TrackProducer_cfi.TrackProducer.clone(
178  AlgorithmName = cms.string('mixedTripletStep'),
179  src = 'mixedTripletStepTrackCandidates',
180  Fitter = cms.string('FlexibleKFFittingSmoother')
181 )
182 
183 # TRACK SELECTION AND QUALITY FLAG SETTING.
184 import RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi
185 mixedTripletStepSelector = RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.multiTrackSelector.clone(
186  src='mixedTripletStepTracks',
187  useAnyMVA = cms.bool(True),
188  GBRForestLabel = cms.string('MVASelectorIter4'),
189  trackSelectors= cms.VPSet(
190  RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.looseMTS.clone(
191  name = 'mixedTripletStepVtxLoose',
192  chi2n_par = 1.0,
193  res_par = ( 0.003, 0.001 ),
194  minNumberLayers = 3,
195  maxNumberLostLayers = 1,
196  minNumber3DLayers = 2,
197  d0_par1 = ( 1.2, 3.0 ),
198  dz_par1 = ( 1.2, 3.0 ),
199  d0_par2 = ( 1.3, 3.0 ),
200  dz_par2 = ( 1.3, 3.0 )
201  ),
202  RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.looseMTS.clone(
203  name = 'mixedTripletStepTrkLoose',
204  chi2n_par = 0.6,
205  res_par = ( 0.003, 0.001 ),
206  minNumberLayers = 4,
207  maxNumberLostLayers = 1,
208  minNumber3DLayers = 3,
209  d0_par1 = ( 1.1, 4.0 ),
210  dz_par1 = ( 1.1, 4.0 ),
211  d0_par2 = ( 1.1, 4.0 ),
212  dz_par2 = ( 1.1, 4.0 )
213  ),
214  RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.tightMTS.clone(
215  name = 'mixedTripletStepVtxTight',
216  preFilterName = 'mixedTripletStepVtxLoose',
217  chi2n_par = 0.6,
218  res_par = ( 0.003, 0.001 ),
219  minNumberLayers = 3,
220  maxNumberLostLayers = 1,
221  minNumber3DLayers = 3,
222  d0_par1 = ( 1.1, 3.0 ),
223  dz_par1 = ( 1.1, 3.0 ),
224  d0_par2 = ( 1.2, 3.0 ),
225  dz_par2 = ( 1.2, 3.0 )
226  ),
227  RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.tightMTS.clone(
228  name = 'mixedTripletStepTrkTight',
229  preFilterName = 'mixedTripletStepTrkLoose',
230  chi2n_par = 0.4,
231  res_par = ( 0.003, 0.001 ),
232  minNumberLayers = 5,
233  maxNumberLostLayers = 1,
234  minNumber3DLayers = 4,
235  d0_par1 = ( 1.0, 4.0 ),
236  dz_par1 = ( 1.0, 4.0 ),
237  d0_par2 = ( 1.0, 4.0 ),
238  dz_par2 = ( 1.0, 4.0 )
239  ),
240  RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.highpurityMTS.clone(
241  name = 'mixedTripletStepVtx',
242  preFilterName = 'mixedTripletStepVtxTight',
243  chi2n_par = 0.4,
244  res_par = ( 0.003, 0.001 ),
245  minNumberLayers = 3,
246  maxNumberLostLayers = 1,
247  minNumber3DLayers = 3,
248  max_minMissHitOutOrIn = 1,
249  max_lostHitFraction = 1.0,
250  d0_par1 = ( 1.0, 3.0 ),
251  dz_par1 = ( 1.0, 3.0 ),
252  d0_par2 = ( 1.1, 3.0 ),
253  dz_par2 = ( 1.1, 3.0 )
254  ),
255  RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.highpurityMTS.clone(
256  name = 'mixedTripletStepTrk',
257  preFilterName = 'mixedTripletStepTrkTight',
258  chi2n_par = 0.25,
259  res_par = ( 0.003, 0.001 ),
260  minNumberLayers = 5,
261  maxNumberLostLayers = 0,
262  minNumber3DLayers = 4,
263  max_minMissHitOutOrIn = 1,
264  max_lostHitFraction = 1.0,
265  d0_par1 = ( 0.8, 4.0 ),
266  dz_par1 = ( 0.8, 4.0 ),
267  d0_par2 = ( 0.8, 4.0 ),
268  dz_par2 = ( 0.8, 4.0 )
269  )
270  ) #end of vpset
271  ) #end of clone
272 
273 import RecoTracker.FinalTrackSelectors.trackListMerger_cfi
274 mixedTripletStep = RecoTracker.FinalTrackSelectors.trackListMerger_cfi.trackListMerger.clone(
275  TrackProducers = cms.VInputTag(cms.InputTag('mixedTripletStepTracks'),
276  cms.InputTag('mixedTripletStepTracks')),
277  hasSelector=cms.vint32(1,1),
278  shareFrac=cms.double(0.11),
279  indivShareFrac=cms.vdouble(0.11,0.11),
280  selectedTrackQuals = cms.VInputTag(cms.InputTag("mixedTripletStepSelector","mixedTripletStepVtx"),
281  cms.InputTag("mixedTripletStepSelector","mixedTripletStepTrk")),
282  setsToMerge = cms.VPSet( cms.PSet( tLists=cms.vint32(0,1), pQual=cms.bool(True) )),
283  writeOnlyTrkQuals=cms.bool(True)
284 )
285 
286 
287 MixedTripletStep = cms.Sequence(mixedTripletStepClusters*
288  mixedTripletStepSeedLayersA*
289  mixedTripletStepSeedsA*
290  mixedTripletStepSeedLayersB*
291  mixedTripletStepSeedsB*
292  mixedTripletStepSeeds*
293  mixedTripletStepTrackCandidates*
294  mixedTripletStepTracks*
295  mixedTripletStepSelector*
296  mixedTripletStep)
297 
OppositeMaterialPropagator_cfi
PixelTripletLargeTipGenerator_cfi
CkfTrackCandidates_cfi
MaterialPropagator_cfi
trackClusterRemover_cfi
Chi2MeasurementEstimatorESProducer_cfi
TrajectoryFilter_cff
TrajectoryCleanerBySharedHits_cfi
TrackProducer_cfi
GroupedCkfTrajectoryBuilder_cfi
GlobalCombinedSeeds_cfi
GlobalSeedsFromTriplets_cff
Generated for CMSSW Reference Manual by doxygen 1.8.5