CMS 3D CMS Logo

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