HiTracking/python/hiRegitPixelPairStep_cff.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 ###################################################################################
0004 # pp iterative tracking modified for hiOffline reco (the vertex is the one reconstructed in HI)
0005 ################################### 2nd step: pixel pairs
0006
0007 from RecoHI.HiTracking.HITrackingRegionProducer_cfi import *
0008
0009 ###################################
0010 import RecoTracker.IterativeTracking.PixelPairStep_cff
0011 from RecoTracker.IterativeTracking.PixelPairStep_cff import pixelPairStepTrajectoryBuilder,pixelPairStepTrajectoryFilter,pixelPairStepTrajectoryFilterBase,pixelPairStepTrajectoryFilterInOut,pixelPairStepTrajectoryFilterShape
0012
0013 # NEW CLUSTERS (remove previously used clusters)
0014 hiRegitPixelPairStepClusters = cms.EDProducer("HITrackClusterRemover",
0015                                      clusterLessSolution= cms.bool(True),
0016                                      oldClusterRemovalInfo = cms.InputTag("hiRegitLowPtTripletStepClusters"),
0017                                      trajectories = cms.InputTag("hiRegitLowPtTripletStepTracks"),
0018                                      overrideTrkQuals = cms.InputTag('hiRegitLowPtTripletStepSelector','hiRegitLowPtTripletStep'),
0019                                      TrackQuality = cms.string('highPurity'),
0020                                      pixelClusters = cms.InputTag("siPixelClusters"),
0021                                      stripClusters = cms.InputTag("siStripClusters"),
0022                                      Common = cms.PSet(
0023                             maxChi2 = cms.double(9.0),
0024                             ),
0025                          Strip = cms.PSet(
0026                             maxChi2 = cms.double(9.0),
0027                             #Yen-Jie's mod to preserve merged clusters
0028                             maxSize = cms.uint32(2)
0029                             )
0030                                      )
0031
0032
0033 # SEEDING LAYERS
0034 hiRegitPixelPairStepSeedLayers =  RecoTracker.IterativeTracking.PixelPairStep_cff.pixelPairStepSeedLayers.clone(
0035     BPix = dict(skipClusters = 'hiRegitPixelPairStepClusters'),
0036     FPix = dict(skipClusters = 'hiRegitPixelPairStepClusters')
0037 )
0038
0039
0040 # seeding
0041 hiRegitPixelPairStepSeeds     = RecoTracker.IterativeTracking.PixelPairStep_cff.pixelPairStepSeeds.clone(
0042     RegionFactoryPSet = HiTrackingRegionFactoryFromJetsBlock.clone(
0043     RegionPSet = dict(ptMin = 1.2)
0044     ),
0045     ClusterCheckPSet = dict(doClusterCheck = False), # do not check for max number of clusters pixel or strips
0046     OrderedHitsFactoryPSet = dict(SeedingLayers = 'hiRegitPixelPairStepSeedLayers'),
0047 )
0048
0049 # building: feed the new-named seeds
0050 hiRegitPixelPairStepTrajectoryFilter = RecoTracker.IterativeTracking.PixelPairStep_cff.pixelPairStepTrajectoryFilterBase.clone()
0051
0052 hiRegitPixelPairStepTrajectoryBuilder = RecoTracker.IterativeTracking.PixelPairStep_cff.pixelPairStepTrajectoryBuilder.clone(
0053     trajectoryFilter     = dict(refToPSet_ = 'hiRegitPixelPairStepTrajectoryFilter'),
0054     clustersToSkip       = 'hiRegitPixelPairStepClusters',
0055 )
0056
0057 # trackign candidate
0058 hiRegitPixelPairStepTrackCandidates        =  RecoTracker.IterativeTracking.PixelPairStep_cff._pixelPairStepTrackCandidatesCkf.clone(
0059     src               = 'hiRegitPixelPairStepSeeds',
0060     TrajectoryBuilderPSet = cms.PSet(refToPSet_ = cms.string('hiRegitPixelPairStepTrajectoryBuilder')),
0061     maxNSeeds = 100000
0062 )
0063
0064 # fitting: feed new-names
0065 hiRegitPixelPairStepTracks                 = RecoTracker.IterativeTracking.PixelPairStep_cff.pixelPairStepTracks.clone(
0066     src                 = 'hiRegitPixelPairStepTrackCandidates',
0067     AlgorithmName = 'pixelPairStep',
0068 )
0069
0070
0071 # Track selection
0072 import RecoHI.HiTracking.hiMultiTrackSelector_cfi
0073 hiRegitPixelPairStepSelector = RecoHI.HiTracking.hiMultiTrackSelector_cfi.hiMultiTrackSelector.clone(
0074     src = 'hiRegitPixelPairStepTracks',
0075     trackSelectors = cms.VPSet(
0076        RecoHI.HiTracking.hiMultiTrackSelector_cfi.hiLooseMTS.clone(
0077            name = 'hiRegitPixelPairStepLoose',
0078            d0_par2 = [9999.0, 0.0],
0079            dz_par2 = [9999.0, 0.0],
0080            applyAdaptedPVCuts = False
0081        ), #end of pset
0082        RecoHI.HiTracking.hiMultiTrackSelector_cfi.hiTightMTS.clone(
0083            name = 'hiRegitPixelPairStepTight',
0084            preFilterName = 'hiRegitPixelPairStepLoose',
0085            d0_par2 = [9999.0, 0.0],
0086            dz_par2 = [9999.0, 0.0],
0087            applyAdaptedPVCuts = False
0088        ),
0089        RecoHI.HiTracking.hiMultiTrackSelector_cfi.hiHighpurityMTS.clone(
0090            name = 'hiRegitPixelPairStep',
0091            preFilterName = 'hiRegitPixelPairStepTight',
0092            d0_par2 = [9999.0, 0.0],
0093            dz_par2 = [9999.0, 0.0],
0094            applyAdaptedPVCuts = False
0095        ),
0096     ) #end of vpset
0097 ) #end of clone
0098
0099 hiRegitPixelPairStepTask = cms.Task(hiRegitPixelPairStepClusters,
0100                                     hiRegitPixelPairStepSeedLayers,
0101                                     hiRegitPixelPairStepSeeds,
0102                                     hiRegitPixelPairStepTrackCandidates,
0103                                     hiRegitPixelPairStepTracks,
0104                                     hiRegitPixelPairStepSelector)
0105 hiRegitPixelPairStep = cms.Sequence(hiRegitPixelPairStepTask)