CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
HiRegitMuonPixelPairStep_cff.py
Go to the documentation of this file.
2 
3 ###################################################################################
4 # pp iterative tracking modified for hiOffline reco (the vertex is the one reconstructed in HI)
5 ################################### 2nd step: pixel pairs
6 
8 HiTrackingRegionFactoryFromSTAMuonsBlock.MuonTrackingRegionBuilder.vertexCollection = cms.InputTag("hiSelectedVertex")
9 HiTrackingRegionFactoryFromSTAMuonsBlock.MuonSrc= cms.InputTag("standAloneMuons","UpdatedAtVtx")
10 
11 HiTrackingRegionFactoryFromSTAMuonsBlock.MuonTrackingRegionBuilder.UseVertex = True
12 
13 HiTrackingRegionFactoryFromSTAMuonsBlock.MuonTrackingRegionBuilder.Phi_fixed = True
14 HiTrackingRegionFactoryFromSTAMuonsBlock.MuonTrackingRegionBuilder.Eta_fixed = True
15 HiTrackingRegionFactoryFromSTAMuonsBlock.MuonTrackingRegionBuilder.DeltaPhi = 0.3
16 HiTrackingRegionFactoryFromSTAMuonsBlock.MuonTrackingRegionBuilder.DeltaEta = 0.2
17 
18 ###################################
20 
21 # NEW CLUSTERS (remove previously used clusters)
22 hiRegitMuPixelPairStepClusters = RecoTracker.IterativeTracking.PixelPairStep_cff.pixelPairStepClusters.clone(
23  trajectories = cms.InputTag("hiRegitMuInitialStepTracks"),
24  overrideTrkQuals = cms.InputTag('hiRegitMuInitialStepSelector','hiRegitMuInitialStep'),
25  oldClusterRemovalInfo = cms.InputTag(""),
26  TrackQuality = cms.string('tight')
27 )
28 
29 
30 # SEEDING LAYERS
31 hiRegitMuPixelPairStepSeedLayers = RecoTracker.IterativeTracking.PixelPairStep_cff.pixelPairStepSeedLayers.clone()
32 hiRegitMuPixelPairStepSeedLayers.BPix.skipClusters = cms.InputTag('hiRegitMuPixelPairStepClusters')
33 hiRegitMuPixelPairStepSeedLayers.FPix.skipClusters = cms.InputTag('hiRegitMuPixelPairStepClusters')
34 
35 
36 
37 # seeding
38 hiRegitMuPixelPairStepSeeds = RecoTracker.IterativeTracking.PixelPairStep_cff.pixelPairStepSeeds.clone()
39 hiRegitMuPixelPairStepSeeds.RegionFactoryPSet = HiTrackingRegionFactoryFromSTAMuonsBlock.clone()
40 hiRegitMuPixelPairStepSeeds.ClusterCheckPSet.doClusterCheck = False # do not check for max number of clusters pixel or strips
41 hiRegitMuPixelPairStepSeeds.RegionFactoryPSet.MuonTrackingRegionBuilder.Pt_min = 1.0
42 hiRegitMuPixelPairStepSeeds.RegionFactoryPSet.MuonTrackingRegionBuilder.DeltaR = 0.01 # default = 0.2
43 hiRegitMuPixelPairStepSeeds.RegionFactoryPSet.MuonTrackingRegionBuilder.DeltaZ = 0.09 # this give you the length
44 hiRegitMuPixelPairStepSeeds.RegionFactoryPSet.MuonTrackingRegionBuilder.Rescale_Dz = 0. # max(DeltaZ_Region,Rescale_Dz*vtx->zError())
45 hiRegitMuPixelPairStepSeeds.OrderedHitsFactoryPSet.SeedingLayers = 'hiRegitMuPixelPairStepSeedLayers'
46 
47 
48 # building: feed the new-named seeds
49 hiRegitMuPixelPairStepTrajectoryFilterBase = RecoTracker.IterativeTracking.PixelPairStep_cff.pixelPairStepTrajectoryFilterBase.clone()
50 hiRegitMuPixelPairStepTrajectoryFilterBase.minPt = 0.8
51 hiRegitMuPixelPairStepTrajectoryFilterBase.minimumNumberOfHits = 6
52 hiRegitMuPixelPairStepTrajectoryFilterBase.minHitsMinPt = 4
53 
54 hiRegitMuPixelPairStepTrajectoryFilter = RecoTracker.IterativeTracking.PixelPairStep_cff.pixelPairStepTrajectoryFilter.clone()
55 hiRegitMuPixelPairStepTrajectoryFilter.filters = cms.VPSet(
56  cms.PSet( refToPSet_ = cms.string('hiRegitMuPixelPairStepTrajectoryFilterBase')),
57  cms.PSet( refToPSet_ = cms.string('pixelPairStepTrajectoryFilterShape')))
58 
59 
60 hiRegitMuPixelPairStepTrajectoryBuilder = RecoTracker.IterativeTracking.PixelPairStep_cff.pixelPairStepTrajectoryBuilder.clone(
61  trajectoryFilter = cms.PSet(
62  refToPSet_ = cms.string('hiRegitMuPixelPairStepTrajectoryFilter')
63  ),
64  minNrOfHitsForRebuild = 6 #change from default 4
65 )
66 
67 # trackign candidate
68 hiRegitMuPixelPairStepTrackCandidates = RecoTracker.IterativeTracking.PixelPairStep_cff.pixelPairStepTrackCandidates.clone(
69  src = cms.InputTag('hiRegitMuPixelPairStepSeeds'),
70  TrajectoryBuilder = 'hiRegitMuPixelPairStepTrajectoryBuilder',
71  clustersToSkip = cms.InputTag("hiRegitMuPixelPairStepClusters"),
72  maxNSeeds = cms.uint32(1000000)
73  )
74 
75 # fitting: feed new-names
76 hiRegitMuPixelPairStepTracks = RecoTracker.IterativeTracking.PixelPairStep_cff.pixelPairStepTracks.clone(
77  AlgorithmName = cms.string('hiRegitMuPixelPairStep'),
78  src = 'hiRegitMuPixelPairStepTrackCandidates',
79  clustersToSkip = cms.InputTag('hiRegitMuPixelPairStepClusters'),
80 )
81 
82 
83 import RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi
85 hiRegitMuPixelPairStepSelector = RecoHI.HiTracking.hiMultiTrackSelector_cfi.hiMultiTrackSelector.clone(
86  src ='hiRegitMuPixelPairStepTracks',
87  vertices = cms.InputTag("hiSelectedVertex"),
88  useAnyMVA = cms.bool(True),
89  GBRForestLabel = cms.string('HIMVASelectorIter6'),
90  GBRForestVars = cms.vstring(['chi2perdofperlayer', 'dxyperdxyerror', 'dzperdzerror', 'nhits', 'nlayers', 'eta']),
91  trackSelectors= cms.VPSet(
92  RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.looseMTS.clone(
93  name = 'hiRegitMuPixelPairStepLoose',
94  min_nhits = cms.uint32(8)
95  ), #end of pset
96  RecoHI.HiTracking.hiMultiTrackSelector_cfi.hiTightMTS.clone(
97  name = 'hiRegitMuPixelPairStepTight',
98  preFilterName = 'hiRegitMuPixelPairStepLoose',
99  min_nhits = cms.uint32(8),
100  useMVA = cms.bool(True),
101  minMVA = cms.double(-0.58)
102  ),
103  RecoHI.HiTracking.hiMultiTrackSelector_cfi.hiHighpurityMTS.clone(
104  name = 'hiRegitMuPixelPairStep',
105  preFilterName = 'hiRegitMuPixelPairStepTight',
106  min_nhits = cms.uint32(8),
107  useMVA = cms.bool(True),
108  minMVA = cms.double(0.77)
109  ),
110  ) #end of vpset
111 )
112 
113 hiRegitMuonPixelPairStep = cms.Sequence(hiRegitMuPixelPairStepClusters*
114  hiRegitMuPixelPairStepSeedLayers*
115  hiRegitMuPixelPairStepSeeds*
116  hiRegitMuPixelPairStepTrackCandidates*
117  hiRegitMuPixelPairStepTracks*
118  hiRegitMuPixelPairStepSelector)
pp iterative tracking modified for hiOffline reco (the vertex is the one reconstructed in HI) 3rd ste...