CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
Phase1PU70_LowPtTripletStep_cff.py
Go to the documentation of this file.
2 
3 # NEW CLUSTERS (remove previously used clusters)
4 lowPtTripletStepClusters = cms.EDProducer("TrackClusterRemover",
5  clusterLessSolution= cms.bool(True),
6  trajectories = cms.InputTag("lowPtQuadStepTracks"),
7  overrideTrkQuals = cms.InputTag('lowPtQuadStepSelector','lowPtQuadStep'),
8  TrackQuality = cms.string('highPurity'),
9  minNumberOfLayersWithMeasBeforeFiltering = cms.int32(0),
10  pixelClusters = cms.InputTag("siPixelClusters"),
11  stripClusters = cms.InputTag("siStripClusters"),
12  Common = cms.PSet(
13  maxChi2 = cms.double(9.0)
14  )
15 )
16 
17 # SEEDING LAYERS
19 lowPtTripletStepSeedLayers = RecoTracker.TkSeedingLayers.PixelLayerTriplets_cfi.PixelLayerTriplets.clone(
20  layerList = cms.vstring('BPix1+BPix2+BPix3', 'BPix2+BPix3+BPix4',
21  'BPix1+BPix3+BPix4', 'BPix1+BPix2+BPix4',
22  'BPix2+BPix3+FPix1_pos', 'BPix2+BPix3+FPix1_neg',
23  'BPix1+BPix2+FPix1_pos', 'BPix1+BPix2+FPix1_neg',
24  'BPix1+BPix3+FPix1_pos', 'BPix1+BPix3+FPix1_neg',
25  'BPix2+FPix1_pos+FPix2_pos', 'BPix2+FPix1_neg+FPix2_neg',
26  'BPix1+FPix1_pos+FPix2_pos', 'BPix1+FPix1_neg+FPix2_neg',
27  'BPix1+BPix2+FPix2_pos', 'BPix1+BPix2+FPix2_neg',
28  'FPix1_pos+FPix2_pos+FPix3_pos', 'FPix1_neg+FPix2_neg+FPix3_neg',
29  'BPix1+FPix2_pos+FPix3_pos', 'BPix1+FPix2_neg+FPix3_neg',
30  'BPix1+FPix1_pos+FPix3_pos', 'BPix1+FPix1_neg+FPix3_neg')
31  )
32 lowPtTripletStepSeedLayers.BPix.skipClusters = cms.InputTag('lowPtTripletStepClusters')
33 lowPtTripletStepSeedLayers.FPix.skipClusters = cms.InputTag('lowPtTripletStepClusters')
34 
35 # SEEDS
37 from RecoTracker.TkTrackingRegions.GlobalTrackingRegionFromBeamSpot_cfi import RegionPsetFomBeamSpotBlock
38 lowPtTripletStepSeeds = RecoTracker.TkSeedGenerator.GlobalSeedsFromTriplets_cff.globalSeedsFromTriplets.clone(
39  RegionFactoryPSet = RegionPsetFomBeamSpotBlock.clone(
40  ComponentName = cms.string('GlobalRegionProducerFromBeamSpot'),
41  RegionPSet = RegionPsetFomBeamSpotBlock.RegionPSet.clone(
42  ptMin = 0.35,
43  originRadius = 0.015,
44  nSigmaZ = 4.0
45  )
46  )
47  )
48 lowPtTripletStepSeeds.OrderedHitsFactoryPSet.SeedingLayers = 'lowPtTripletStepSeedLayers'
49 
51 lowPtTripletStepSeeds.OrderedHitsFactoryPSet.GeneratorPSet.SeedComparitorPSet.ComponentName = 'LowPtClusterShapeSeedComparitor'
52 lowPtTripletStepSeeds.ClusterCheckPSet.doClusterCheck = cms.bool(False)
53 lowPtTripletStepSeeds.OrderedHitsFactoryPSet.GeneratorPSet.maxElement = cms.uint32(0)
54 
55 # QUALITY CUTS DURING TRACK BUILDING
56 import TrackingTools.TrajectoryFiltering.TrajectoryFilterESProducer_cfi
57 lowPtTripletStepStandardTrajectoryFilter = TrackingTools.TrajectoryFiltering.TrajectoryFilterESProducer_cfi.trajectoryFilterESProducer.clone(
58  ComponentName = 'lowPtTripletStepStandardTrajectoryFilter',
59  filterPset = TrackingTools.TrajectoryFiltering.TrajectoryFilterESProducer_cfi.trajectoryFilterESProducer.filterPset.clone(
60  minimumNumberOfHits = 3,
61  minPt = 0.075
62  )
63  )
64 
66 # Composite filter
68 lowPtTripletStepTrajectoryFilter = TrackingTools.TrajectoryFiltering.CompositeTrajectoryFilterESProducer_cfi.compositeTrajectoryFilterESProducer.clone(
69  ComponentName = cms.string('lowPtTripletStepTrajectoryFilter'),
70  filterNames = cms.vstring('lowPtTripletStepStandardTrajectoryFilter',
71  'clusterShapeTrajectoryFilter')
72  )
73 
75 lowPtTripletStepChi2Est = TrackingTools.KalmanUpdators.Chi2MeasurementEstimatorESProducer_cfi.Chi2MeasurementEstimator.clone(
76  ComponentName = cms.string('lowPtTripletStepChi2Est'),
77  nSigma = cms.double(3.0),
78  MaxChi2 = cms.double(9.0)
79 )
80 
81 # TRACK BUILDING
83 lowPtTripletStepTrajectoryBuilder = RecoTracker.CkfPattern.GroupedCkfTrajectoryBuilderESProducer_cfi.GroupedCkfTrajectoryBuilder.clone(
84  ComponentName = 'lowPtTripletStepTrajectoryBuilder',
85  MeasurementTrackerName = '',
86  trajectoryFilterName = 'lowPtTripletStepTrajectoryFilter',
87  clustersToSkip = cms.InputTag('lowPtTripletStepClusters'),
88  maxCand = 4,
89  estimator = cms.string('lowPtTripletStepChi2Est'),
90  maxDPhiForLooperReconstruction = cms.double(2.0),
91  # 0.63 GeV is the maximum pT for a charged particle to loop within the 1.1m radius
92  # of the outermost Tracker barrel layer (with B=3.8T)
93  maxPtForLooperReconstruction = cms.double(0.7)
94  )
95 
96 # MAKING OF TRACK CANDIDATES
98 lowPtTripletStepTrackCandidates = RecoTracker.CkfPattern.CkfTrackCandidates_cfi.ckfTrackCandidates.clone(
99  src = cms.InputTag('lowPtTripletStepSeeds'),
100  ### these two parameters are relevant only for the CachingSeedCleanerBySharedInput
101  numHitsForSeedCleaner = cms.int32(50),
102  onlyPixelHitsForSeedCleaner = cms.bool(True),
103  TrajectoryBuilder = 'lowPtTripletStepTrajectoryBuilder',
104  doSeedingRegionRebuilding = True,
105  useHitsSplitting = True
106 )
107 
108 # TRACK FITTING
110 lowPtTripletStepTracks = RecoTracker.TrackProducer.TrackProducer_cfi.TrackProducer.clone(
111  src = 'lowPtTripletStepTrackCandidates',
112  AlgorithmName = cms.string('iter3'),
113  Fitter = cms.string('FlexibleKFFittingSmoother'),
114  TTRHBuilder=cms.string('WithTrackAngle')
115  )
116 
117 from TrackingTools.TrajectoryCleaning.TrajectoryCleanerBySharedHits_cfi import trajectoryCleanerBySharedHits
118 lowPtTripletStepTrajectoryCleanerBySharedHits = trajectoryCleanerBySharedHits.clone(
119  ComponentName = cms.string('lowPtTripletStepTrajectoryCleanerBySharedHits'),
120  fractionShared = cms.double(0.09),
121  allowSharedFirstHit = cms.bool(True)
122  )
123 lowPtTripletStepTrackCandidates.TrajectoryCleaner = 'lowPtTripletStepTrajectoryCleanerBySharedHits'
124 
125 # Final selection
126 import RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi
127 lowPtTripletStepSelector = RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.multiTrackSelector.clone(
128  src='lowPtTripletStepTracks',
129  trackSelectors= cms.VPSet(
130  RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.looseMTS.clone(
131  name = 'lowPtTripletStepLoose',
132  chi2n_par = 2.0,
133  res_par = ( 0.003, 0.002 ),
134  minNumberLayers = 3,
135  maxNumberLostLayers = 2,
136  minNumber3DLayers = 3,
137  d0_par1 = ( 0.8, 4.0 ),
138  dz_par1 = ( 0.7, 4.0 ),
139  d0_par2 = ( 0.5, 4.0 ),
140  dz_par2 = ( 0.5, 4.0 )
141  ), #end of pset
142  RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.tightMTS.clone(
143  name = 'lowPtTripletStepTight',
144  preFilterName = 'lowPtTripletStepLoose',
145  chi2n_par = 0.9,
146  res_par = ( 0.003, 0.002 ),
147  minNumberLayers = 3,
148  maxNumberLostLayers = 2,
149  minNumber3DLayers = 3,
150  d0_par1 = ( 0.7, 4.0 ),
151  dz_par1 = ( 0.6, 4.0 ),
152  d0_par2 = ( 0.4, 4.0 ),
153  dz_par2 = ( 0.4, 4.0 )
154  ),
155  RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.highpurityMTS.clone(
156  name = 'lowPtTripletStep',
157  preFilterName = 'lowPtTripletStepTight',
158  chi2n_par = 0.5,
159  res_par = ( 0.003, 0.001 ),
160  minNumberLayers = 3,
161  maxNumberLostLayers = 2,
162  minNumber3DLayers = 3,
163  d0_par1 = ( 0.6, 4.0 ),
164  dz_par1 = ( 0.45, 4.0 ),
165  d0_par2 = ( 0.3, 4.0 ),
166  dz_par2 = ( 0.4, 4.0 )
167  ),
168  ) #end of vpset
169  ) #end of clone
170 
171 # Final sequence
172 LowPtTripletStep = cms.Sequence(lowPtTripletStepClusters*
173  lowPtTripletStepSeedLayers*
174  lowPtTripletStepSeeds*
175  lowPtTripletStepTrackCandidates*
176  lowPtTripletStepTracks*
177  lowPtTripletStepSelector)