CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
RunI_PixelLessStep_cff.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
2 
3 ##########################################################################
4 # Large impact parameter tracking using TIB/TID/TEC stereo layer seeding #
5 ##########################################################################
6 
8 pixelLessStepClusters = trackClusterRemover.clone(
9  maxChi2 = cms.double(9.0),
10  trajectories = cms.InputTag("mixedTripletStepTracks"),
11  pixelClusters = cms.InputTag("siPixelClusters"),
12  stripClusters = cms.InputTag("siStripClusters"),
13  oldClusterRemovalInfo = cms.InputTag("mixedTripletStepClusters"),
14  overrideTrkQuals = cms.InputTag('mixedTripletStep'),
15  TrackQuality = cms.string('highPurity'),
16  minNumberOfLayersWithMeasBeforeFiltering = cms.int32(0),
17 )
18 
19 # SEEDING LAYERS
20 pixelLessStepSeedLayers = cms.EDProducer("SeedingLayersEDProducer",
21  layerList = cms.vstring('TIB1+TIB2',
22  'TID1_pos+TID2_pos','TID2_pos+TID3_pos',
23  'TEC1_pos+TEC2_pos','TEC2_pos+TEC3_pos','TEC3_pos+TEC4_pos','TEC3_pos+TEC5_pos','TEC4_pos+TEC5_pos',
24  'TID1_neg+TID2_neg','TID2_neg+TID3_neg',
25  'TEC1_neg+TEC2_neg','TEC2_neg+TEC3_neg','TEC3_neg+TEC4_neg','TEC3_neg+TEC5_neg','TEC4_neg+TEC5_neg'),
26  TIB = cms.PSet(
27  TTRHBuilder = cms.string('WithTrackAngle'), clusterChargeCut = cms.PSet(refToPSet_ = cms.string('SiStripClusterChargeCutNone')),
28  matchedRecHits = cms.InputTag("siStripMatchedRecHits","matchedRecHit"),
29  skipClusters = cms.InputTag('pixelLessStepClusters')
30  ),
31  TID = cms.PSet(
32  matchedRecHits = cms.InputTag("siStripMatchedRecHits","matchedRecHit"),
33  skipClusters = cms.InputTag('pixelLessStepClusters'),
34  useRingSlector = cms.bool(True),
35  TTRHBuilder = cms.string('WithTrackAngle'), clusterChargeCut = cms.PSet(refToPSet_ = cms.string('SiStripClusterChargeCutNone')),
36  minRing = cms.int32(1),
37  maxRing = cms.int32(2)
38  ),
39  TEC = cms.PSet(
40  matchedRecHits = cms.InputTag("siStripMatchedRecHits","matchedRecHit"),
41  skipClusters = cms.InputTag('pixelLessStepClusters'),
42  useRingSlector = cms.bool(True),
43  TTRHBuilder = cms.string('WithTrackAngle'), clusterChargeCut = cms.PSet(refToPSet_ = cms.string('SiStripClusterChargeCutNone')),
44  minRing = cms.int32(1),
45  maxRing = cms.int32(2)
46  )
47 )
48 
49 # SEEDS
51 pixelLessStepSeeds = RecoTracker.TkSeedGenerator.GlobalMixedSeeds_cff.globalMixedSeeds.clone()
52 pixelLessStepSeeds.OrderedHitsFactoryPSet.SeedingLayers = 'pixelLessStepSeedLayers'
53 pixelLessStepSeeds.RegionFactoryPSet.RegionPSet.ptMin = 0.7
54 pixelLessStepSeeds.RegionFactoryPSet.RegionPSet.originHalfLength = 10.0
55 pixelLessStepSeeds.RegionFactoryPSet.RegionPSet.originRadius = 2.0
56 
57 pixelLessStepSeeds.SeedComparitorPSet = cms.PSet(
58  ComponentName = cms.string('PixelClusterShapeSeedComparitor'),
59  FilterAtHelixStage = cms.bool(True),
60  FilterPixelHits = cms.bool(False),
61  FilterStripHits = cms.bool(True),
62  ClusterShapeHitFilterName = cms.string('ClusterShapeHitFilter'),
63  ClusterShapeCacheSrc = cms.InputTag("siPixelClusterShapeCache") # not really needed here since FilterPixelHits=False
64  )
65 
66 # QUALITY CUTS DURING TRACK BUILDING
68 pixelLessStepTrajectoryFilter = TrackingTools.TrajectoryFiltering.TrajectoryFilter_cff.CkfBaseTrajectoryFilter_block.clone(
69  maxLostHits = 0,
70  minimumNumberOfHits = 4,
71  minPt = 0.1
72  )
73 
75 pixelLessStepChi2Est = TrackingTools.KalmanUpdators.Chi2MeasurementEstimator_cfi.Chi2MeasurementEstimator.clone(
76  ComponentName = cms.string('pixelLessStepChi2Est'),
77  nSigma = cms.double(3.0),
78  MaxChi2 = cms.double(16.0),
79  clusterChargeCut = cms.PSet(refToPSet_ = cms.string('SiStripClusterChargeCutNone'))
80 )
81 
82 # TRACK BUILDING
84 pixelLessStepTrajectoryBuilder = RecoTracker.CkfPattern.GroupedCkfTrajectoryBuilder_cfi.GroupedCkfTrajectoryBuilder.clone(
85  MeasurementTrackerName = '',
86  trajectoryFilter = cms.PSet(refToPSet_ = cms.string('pixelLessStepTrajectoryFilter')),
87  minNrOfHitsForRebuild = 4,
88  maxCand = 2,
89  alwaysUseInvalidHits = False,
90  estimator = cms.string('pixelLessStepChi2Est'),
91  maxDPhiForLooperReconstruction = cms.double(2.0),
92  maxPtForLooperReconstruction = cms.double(0.7)
93  )
94 
95 # MAKING OF TRACK CANDIDATES
97 pixelLessStepTrackCandidates = RecoTracker.CkfPattern.CkfTrackCandidates_cfi.ckfTrackCandidates.clone(
98  src = cms.InputTag('pixelLessStepSeeds'),
99  clustersToSkip = cms.InputTag('pixelLessStepClusters'),
100  ### these two parameters are relevant only for the CachingSeedCleanerBySharedInput
101  numHitsForSeedCleaner = cms.int32(50),
102  #onlyPixelHitsForSeedCleaner = cms.bool(True),
103  TrajectoryBuilderPSet = cms.PSet(refToPSet_ = cms.string('pixelLessStepTrajectoryBuilder'))
104 )
105 
106 from TrackingTools.TrajectoryCleaning.TrajectoryCleanerBySharedHits_cfi import trajectoryCleanerBySharedHits
107 pixelLessStepTrajectoryCleanerBySharedHits = trajectoryCleanerBySharedHits.clone(
108  ComponentName = cms.string('pixelLessStepTrajectoryCleanerBySharedHits'),
109  fractionShared = cms.double(0.19),
110  allowSharedFirstHit = cms.bool(True)
111  )
112 pixelLessStepTrackCandidates.TrajectoryCleaner = 'pixelLessStepTrajectoryCleanerBySharedHits'
113 
114 
115 # TRACK FITTING
117 pixelLessStepTracks = RecoTracker.TrackProducer.TrackProducer_cfi.TrackProducer.clone(
118  src = 'pixelLessStepTrackCandidates',
119  AlgorithmName = cms.string('pixelLessStep'),
120  Fitter = cms.string('FlexibleKFFittingSmoother')
121  )
122 
123 import RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi
124 pixelLessStepSelector = RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.multiTrackSelector.clone(
125  src='pixelLessStepTracks',
126  useAnyMVA = cms.bool(False),
127  GBRForestLabel = cms.string('MVASelectorIter5'),
128  trackSelectors= cms.VPSet(
129  RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.looseMTS.clone(
130  name = 'pixelLessStepLoose',
131  chi2n_par = 0.5,
132  res_par = ( 0.003, 0.001 ),
133  minNumberLayers = 4,
134  maxNumberLostLayers = 1,
135  minNumber3DLayers = 3,
136  d0_par1 = ( 1.3, 4.0 ),
137  dz_par1 = ( 1.3, 4.0 ),
138  d0_par2 = ( 1.3, 4.0 ),
139  dz_par2 = ( 1.3, 4.0 )
140  ),
141  RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.tightMTS.clone(
142  name = 'pixelLessStepTight',
143  preFilterName = 'pixelLessStepLoose',
144  chi2n_par = 0.35,
145  res_par = ( 0.003, 0.001 ),
146  minNumberLayers = 4,
147  maxNumberLostLayers = 0,
148  minNumber3DLayers = 3,
149  d0_par1 = ( 1.1, 4.0 ),
150  dz_par1 = ( 1.1, 4.0 ),
151  d0_par2 = ( 1.1, 4.0 ),
152  dz_par2 = ( 1.1, 4.0 )
153  ),
154  RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.highpurityMTS.clone(
155  name = 'QualityMasks',
156  preFilterName = 'pixelLessStepTight',
157  chi2n_par = 0.2,
158  res_par = ( 0.003, 0.001 ),
159  minNumberLayers = 4,
160  maxNumberLostLayers = 0,
161  minNumber3DLayers = 3,
162  d0_par1 = ( 0.9, 4.0 ),
163  dz_par1 = ( 0.9, 4.0 ),
164  d0_par2 = ( 0.9, 4.0 ),
165  dz_par2 = ( 0.9, 4.0 )
166  ),
167  ),
168  vertices = cms.InputTag("pixelVertices")#end of vpset
169  ) #end of clone
170 
171 
172 PixelLessStep = cms.Sequence(pixelLessStepClusters*
173  pixelLessStepSeedLayers*
174  pixelLessStepSeeds*
175  pixelLessStepTrackCandidates*
176  pixelLessStepTracks*
177  pixelLessStepSelector)
178