IterativeTracking/python/HighPtTripletStep_cff.py

0001 import FWCore.ParameterSet.Config as cms
0002 import RecoTracker.IterativeTracking.iterativeTkConfig as _cfg
0003 from Configuration.Eras.Modifier_fastSim_cff import fastSim
0004
0005 # for dnn classifier
0006 from Configuration.ProcessModifiers.trackdnn_cff import trackdnn
0007 from RecoTracker.IterativeTracking.dnnQualityCuts import qualityCutDictionary
0008
0009 # for no-loopers
0010 from Configuration.ProcessModifiers.trackingNoLoopers_cff import trackingNoLoopers
0011
0012 ### high-pT triplets ###
0013
0014 # NEW CLUSTERS (remove previously used clusters)
0015 highPtTripletStepClusters = _cfg.clusterRemoverForIter('HighPtTripletStep')
0016 for _eraName, _postfix, _era in _cfg.nonDefaultEras():
0017     _era.toReplaceWith(highPtTripletStepClusters, _cfg.clusterRemoverForIter('HighPtTripletStep', _eraName, _postfix))
0018
0019
0020 # SEEDING LAYERS
0021 import RecoTracker.TkSeedingLayers.PixelLayerTriplets_cfi as _PixelLayerTriplets_cfi
0022 highPtTripletStepSeedLayers = _PixelLayerTriplets_cfi.PixelLayerTriplets.clone(
0023     layerList = [
0024         'BPix1+BPix2+BPix3',
0025         'BPix2+BPix3+BPix4',
0026         'BPix1+BPix3+BPix4',
0027         'BPix1+BPix2+BPix4',
0028         'BPix2+BPix3+FPix1_pos', 'BPix2+BPix3+FPix1_neg',
0029         'BPix1+BPix2+FPix1_pos', 'BPix1+BPix2+FPix1_neg',
0030         'BPix1+BPix3+FPix1_pos', 'BPix1+BPix3+FPix1_neg',
0031         'BPix2+FPix1_pos+FPix2_pos', 'BPix2+FPix1_neg+FPix2_neg',
0032         'BPix1+FPix1_pos+FPix2_pos', 'BPix1+FPix1_neg+FPix2_neg',
0033         'BPix1+BPix2+FPix2_pos', 'BPix1+BPix2+FPix2_neg',
0034         'FPix1_pos+FPix2_pos+FPix3_pos', 'FPix1_neg+FPix2_neg+FPix3_neg',
0035         'BPix1+FPix2_pos+FPix3_pos', 'BPix1+FPix2_neg+FPix3_neg',
0036         'BPix1+FPix1_pos+FPix3_pos', 'BPix1+FPix1_neg+FPix3_neg'
0037     ],
0038     BPix = dict(skipClusters = cms.InputTag('highPtTripletStepClusters')),
0039     FPix = dict(skipClusters = cms.InputTag('highPtTripletStepClusters'))
0040 )
0041
0042 from Configuration.Eras.Modifier_trackingPhase2PU140_cff import trackingPhase2PU140
0043 trackingPhase2PU140.toModify(highPtTripletStepSeedLayers,
0044 # combination with gap removed as only source of fakes in current geometry (kept for doc)
0045     layerList = ['BPix1+BPix2+BPix3', 'BPix2+BPix3+BPix4',
0046                  'BPix1+BPix3+BPix4', 'BPix1+BPix2+BPix4',
0047                  'BPix2+BPix3+FPix1_pos', 'BPix2+BPix3+FPix1_neg',
0048                  'BPix1+BPix2+FPix1_pos', 'BPix1+BPix2+FPix1_neg',
0049                  'BPix2+FPix1_pos+FPix2_pos', 'BPix2+FPix1_neg+FPix2_neg',
0050                  'BPix1+FPix1_pos+FPix2_pos', 'BPix1+FPix1_neg+FPix2_neg',
0051 #                 'BPix1+BPix2+FPix2_pos', 'BPix1+BPix2+FPix2_neg',
0052                  'FPix1_pos+FPix2_pos+FPix3_pos', 'FPix1_neg+FPix2_neg+FPix3_neg',
0053                  'BPix1+FPix2_pos+FPix3_pos', 'BPix1+FPix2_neg+FPix3_neg',
0054 #                 'BPix1+FPix1_pos+FPix3_pos', 'BPix1+FPix1_neg+FPix3_neg',
0055                  'FPix2_pos+FPix3_pos+FPix4_pos', 'FPix2_neg+FPix3_neg+FPix4_neg',
0056                  'FPix3_pos+FPix4_pos+FPix5_pos', 'FPix3_neg+FPix4_neg+FPix5_neg',
0057                  'FPix4_pos+FPix5_pos+FPix6_pos', 'FPix4_neg+FPix5_neg+FPix6_neg',
0058                  'FPix5_pos+FPix6_pos+FPix7_pos', 'FPix5_neg+FPix6_neg+FPix7_neg',
0059                  'FPix6_pos+FPix7_pos+FPix8_pos', 'FPix6_neg+FPix7_neg+FPix8_neg',
0060 #  removed as redunant and covering effectively only eta>4   (here for documentation, to be optimized after TDR)
0061 #                 'FPix6_pos+FPix7_pos+FPix9_pos', 'FPix6_neg+FPix7_neg+FPix9_neg']
0062      ]
0063 )
0064 # TrackingRegion
0065 from RecoTracker.TkTrackingRegions.globalTrackingRegionFromBeamSpot_cfi import globalTrackingRegionFromBeamSpot as _globalTrackingRegionFromBeamSpot
0066 highPtTripletStepTrackingRegions = _globalTrackingRegionFromBeamSpot.clone(RegionPSet = dict(
0067     ptMin        = 0.55,
0068     originRadius = 0.02,
0069     nSigmaZ      = 4.0
0070 ))
0071 trackingPhase2PU140.toModify(highPtTripletStepTrackingRegions, RegionPSet = dict(ptMin = 0.7, originRadius = 0.02))
0072
0073 from Configuration.Eras.Modifier_pp_on_XeXe_2017_cff import pp_on_XeXe_2017
0074 from Configuration.ProcessModifiers.pp_on_AA_cff import pp_on_AA
0075 from RecoTracker.TkTrackingRegions.globalTrackingRegionWithVertices_cff import globalTrackingRegionWithVertices as _globalTrackingRegionWithVertices
0076 (pp_on_XeXe_2017 | pp_on_AA).toReplaceWith(highPtTripletStepTrackingRegions,
0077                 _globalTrackingRegionWithVertices.clone(RegionPSet=dict(
0078                     fixedError   = 0.2,
0079                     ptMin        = 0.7,
0080                     originRadius = 0.02
0081                 )
0082                                                                       )
0083 )
0084 from Configuration.Eras.Modifier_highBetaStar_cff import highBetaStar
0085 highBetaStar.toModify(highPtTripletStepTrackingRegions,RegionPSet = dict(
0086      ptMin        = 0.05,
0087      originRadius = 0.2
0088 ))
0089
0090
0091 # seeding
0092 from RecoTracker.TkHitPairs.hitPairEDProducer_cfi import hitPairEDProducer as _hitPairEDProducer
0093 highPtTripletStepHitDoublets = _hitPairEDProducer.clone(
0094     seedingLayers   = 'highPtTripletStepSeedLayers',
0095     trackingRegions = 'highPtTripletStepTrackingRegions',
0096     layerPairs      = [0,1], # layer pairs (0,1), (1,2)
0097     maxElement      = 50000000,
0098     produceIntermediateHitDoublets = True,
0099 )
0100 from RecoTracker.PixelSeeding.caHitTripletEDProducer_cfi import caHitTripletEDProducer as _caHitTripletEDProducer
0101 from RecoTracker.PixelSeeding.pixelTripletHLTEDProducer_cfi import pixelTripletHLTEDProducer as _pixelTripletHLTEDProducer
0102 from RecoTracker.PixelLowPtUtilities.ClusterShapeHitFilterESProducer_cfi import *
0103 import RecoTracker.PixelLowPtUtilities.LowPtClusterShapeSeedComparitor_cfi
0104 highPtTripletStepHitTriplets = _caHitTripletEDProducer.clone(
0105     doublets = 'highPtTripletStepHitDoublets',
0106     extraHitRPhitolerance = _pixelTripletHLTEDProducer.extraHitRPhitolerance,
0107     SeedComparitorPSet = RecoTracker.PixelLowPtUtilities.LowPtClusterShapeSeedComparitor_cfi.LowPtClusterShapeSeedComparitor.clone(),
0108     maxChi2 = dict(
0109         pt1    = 0.8, pt2    = 8,
0110         value1 = 100, value2 = 6,
0111     ),
0112     useBendingCorrection = True,
0113     CAThetaCut           = 0.004,
0114     CAPhiCut             = 0.07,
0115     CAHardPtCut          = 0.3,
0116 )
0117
0118 trackingPhase2PU140.toModify(highPtTripletStepHitTriplets,CAThetaCut = 0.003,CAPhiCut = 0.06,CAHardPtCut = 0.5)
0119 highBetaStar.toModify(highPtTripletStepHitTriplets,CAThetaCut = 0.008,CAPhiCut = 0.14,CAHardPtCut = 0)
0120
0121 from RecoTracker.TkSeedGenerator.seedCreatorFromRegionConsecutiveHitsEDProducer_cff import seedCreatorFromRegionConsecutiveHitsEDProducer as _seedCreatorFromRegionConsecutiveHitsEDProducer
0122 highPtTripletStepSeeds = _seedCreatorFromRegionConsecutiveHitsEDProducer.clone(
0123     seedingHitSets = 'highPtTripletStepHitTriplets',
0124 )
0125
0126 #For FastSim phase1 tracking
0127 import FastSimulation.Tracking.TrajectorySeedProducer_cfi
0128 from FastSimulation.Tracking.SeedingMigration import _hitSetProducerToFactoryPSet
0129 _fastSim_highPtTripletStepSeeds = FastSimulation.Tracking.TrajectorySeedProducer_cfi.trajectorySeedProducer.clone(
0130     trackingRegions = 'highPtTripletStepTrackingRegions',
0131     hitMasks        = cms.InputTag('highPtTripletStepMasks'),
0132     seedFinderSelector = dict( CAHitTripletGeneratorFactory = _hitSetProducerToFactoryPSet(highPtTripletStepHitTriplets),
0133                                layerList = highPtTripletStepSeedLayers.layerList.value(),
0134                                #new parameters required for phase1 seeding
0135                                BPix = dict(TTRHBuilder = 'WithoutRefit', HitProducer = 'TrackingRecHitProducer',),
0136                                FPix = dict(TTRHBuilder = 'WithoutRefit', HitProducer = 'TrackingRecHitProducer',),
0137                                layerPairs = highPtTripletStepHitDoublets.layerPairs.value()
0138                                ))
0139
0140 _fastSim_highPtTripletStepSeeds.seedFinderSelector.CAHitTripletGeneratorFactory.SeedComparitorPSet.ComponentName = 'none'
0141 fastSim.toReplaceWith(highPtTripletStepSeeds,_fastSim_highPtTripletStepSeeds)
0142
0143 # QUALITY CUTS DURING TRACK BUILDING
0144 import TrackingTools.TrajectoryFiltering.TrajectoryFilter_cff as _TrajectoryFilter_cff
0145 _highPtTripletStepTrajectoryFilterBase = _TrajectoryFilter_cff.CkfBaseTrajectoryFilter_block.clone(
0146     minimumNumberOfHits = 3,
0147     minPt               = 0.2,
0148 )
0149 highPtTripletStepTrajectoryFilterBase = _highPtTripletStepTrajectoryFilterBase.clone(
0150     maxCCCLostHits     = 0,
0151     minGoodStripCharge = dict(refToPSet_ = 'SiStripClusterChargeCutLoose')
0152 )
0153 trackingPhase2PU140.toReplaceWith(highPtTripletStepTrajectoryFilterBase, _highPtTripletStepTrajectoryFilterBase)
0154
0155 (pp_on_XeXe_2017 | pp_on_AA).toModify(highPtTripletStepTrajectoryFilterBase, minPt=0.7)
0156 highBetaStar.toModify(highPtTripletStepTrajectoryFilterBase, minPt=0.05)
0157
0158 highPtTripletStepTrajectoryFilter = _TrajectoryFilter_cff.CompositeTrajectoryFilter_block.clone(
0159     filters = [cms.PSet(refToPSet_ = cms.string('highPtTripletStepTrajectoryFilterBase'))]
0160 )
0161
0162 trackingPhase2PU140.toModify(highPtTripletStepTrajectoryFilter,
0163     filters = highPtTripletStepTrajectoryFilter.filters + [cms.PSet(refToPSet_ = cms.string('ClusterShapeTrajectoryFilter'))]
0164 )
0165
0166
0167 highPtTripletStepTrajectoryFilterInOut = highPtTripletStepTrajectoryFilterBase.clone(
0168     minPt               = 0.4,
0169     minimumNumberOfHits = 4,
0170     seedExtension       = 1,
0171     strictSeedExtension = False, # allow inactive
0172     pixelSeedExtension  = False,
0173 )
0174 highBetaStar.toModify(highPtTripletStepTrajectoryFilterInOut, minPt=0.05)
0175
0176 import RecoTracker.MeasurementDet.Chi2ChargeMeasurementEstimator_cfi
0177 highPtTripletStepChi2Est = RecoTracker.MeasurementDet.Chi2ChargeMeasurementEstimator_cfi.Chi2ChargeMeasurementEstimator.clone(
0178     ComponentName = 'highPtTripletStepChi2Est',
0179     nSigma        = 3.0,
0180     MaxChi2       = 30.0,
0181     clusterChargeCut = dict(refToPSet_ = 'SiStripClusterChargeCutLoose'),
0182     pTChargeCutThreshold = 15.
0183 )
0184 trackingPhase2PU140.toModify(highPtTripletStepChi2Est,
0185     clusterChargeCut = dict(refToPSet_ = 'SiStripClusterChargeCutNone'),
0186     MaxChi2 = 20.0
0187 )
0188
0189
0190 # TRACK BUILDING
0191 import RecoTracker.CkfPattern.GroupedCkfTrajectoryBuilder_cfi as _GroupedCkfTrajectoryBuilder_cfi
0192 highPtTripletStepTrajectoryBuilder = _GroupedCkfTrajectoryBuilder_cfi.GroupedCkfTrajectoryBuilderIterativeDefault.clone(
0193     trajectoryFilter     = dict(refToPSet_ = 'highPtTripletStepTrajectoryFilter'),
0194     alwaysUseInvalidHits = True,
0195     maxCand              = 3,
0196     estimator            = 'highPtTripletStepChi2Est',
0197     maxDPhiForLooperReconstruction = 2.0,
0198     # 0.63 GeV is the maximum pT for a charged particle to loop within the 1.1m radius
0199     # of the outermost Tracker barrel layer (with B=3.8T)
0200     maxPtForLooperReconstruction = 0.7,
0201 )
0202 trackingNoLoopers.toModify(highPtTripletStepTrajectoryBuilder,
0203                            maxPtForLooperReconstruction = 0.0)
0204 trackingPhase2PU140.toModify(highPtTripletStepTrajectoryBuilder,
0205     inOutTrajectoryFilter = dict(refToPSet_ = 'highPtTripletStepTrajectoryFilterInOut'),
0206     useSameTrajFilter = False,
0207     maxCand = 3,
0208 )
0209
0210 # MAKING OF TRACK CANDIDATES
0211 import RecoTracker.CkfPattern.CkfTrackCandidates_cfi as _CkfTrackCandidates_cfi
0212 highPtTripletStepTrackCandidates = _CkfTrackCandidates_cfi.ckfTrackCandidatesIterativeDefault.clone(
0213     src = 'highPtTripletStepSeeds',
0214     clustersToSkip = 'highPtTripletStepClusters',
0215     ### these two parameters are relevant only for the CachingSeedCleanerBySharedInput
0216     numHitsForSeedCleaner = 50,
0217     onlyPixelHitsForSeedCleaner = True,
0218     TrajectoryBuilderPSet = dict(refToPSet_ = 'highPtTripletStepTrajectoryBuilder'),
0219     doSeedingRegionRebuilding = True,
0220     useHitsSplitting = True,
0221 )
0222
0223 from Configuration.ProcessModifiers.trackingMkFitHighPtTripletStep_cff import trackingMkFitHighPtTripletStep
0224 import RecoTracker.MkFit.mkFitSeedConverter_cfi as mkFitSeedConverter_cfi
0225 import RecoTracker.MkFit.mkFitIterationConfigESProducer_cfi as mkFitIterationConfigESProducer_cfi
0226 import RecoTracker.MkFit.mkFitProducer_cfi as mkFitProducer_cfi
0227 import RecoTracker.MkFit.mkFitOutputConverter_cfi as mkFitOutputConverter_cfi
0228 highPtTripletStepTrackCandidatesMkFitSeeds = mkFitSeedConverter_cfi.mkFitSeedConverter.clone(
0229     seeds = 'highPtTripletStepSeeds',
0230 )
0231 highPtTripletStepTrackCandidatesMkFitConfig = mkFitIterationConfigESProducer_cfi.mkFitIterationConfigESProducer.clone(
0232     ComponentName = 'highPtTripletStepTrackCandidatesMkFitConfig',
0233     config = 'RecoTracker/MkFit/data/mkfit-phase1-highPtTripletStep.json',
0234 )
0235 highPtTripletStepTrackCandidatesMkFit = mkFitProducer_cfi.mkFitProducer.clone(
0236     seeds = 'highPtTripletStepTrackCandidatesMkFitSeeds',
0237     config = ('', 'highPtTripletStepTrackCandidatesMkFitConfig'),
0238     clustersToSkip = 'highPtTripletStepClusters',
0239 )
0240 trackingMkFitHighPtTripletStep.toReplaceWith(highPtTripletStepTrackCandidates, mkFitOutputConverter_cfi.mkFitOutputConverter.clone(
0241     seeds = 'highPtTripletStepSeeds',
0242     mkFitSeeds = 'highPtTripletStepTrackCandidatesMkFitSeeds',
0243     tracks = 'highPtTripletStepTrackCandidatesMkFit',
0244     candMVASel = True,
0245     candWP = -0.3,
0246 ))
0247 (pp_on_XeXe_2017 | pp_on_AA).toModify(highPtTripletStepTrackCandidatesMkFitConfig, minPt=0.7)
0248
0249 # For Phase2PU140
0250 from TrackingTools.TrajectoryCleaning.TrajectoryCleanerBySharedHits_cfi import trajectoryCleanerBySharedHits as _trajectoryCleanerBySharedHits
0251 highPtTripletStepTrajectoryCleanerBySharedHits = _trajectoryCleanerBySharedHits.clone(
0252     ComponentName       = 'highPtTripletStepTrajectoryCleanerBySharedHits',
0253     fractionShared      = 0.16,
0254     allowSharedFirstHit = True
0255 )
0256 trackingPhase2PU140.toModify(highPtTripletStepTrackCandidates,
0257     TrajectoryCleaner    = 'highPtTripletStepTrajectoryCleanerBySharedHits',
0258     clustersToSkip       = '',
0259     phase2clustersToSkip = 'highPtTripletStepClusters'
0260 )
0261
0262 from Configuration.ProcessModifiers.trackingLST_cff import trackingLST
0263 from RecoTracker.LST.lstOutputConverter_cfi import lstOutputConverter as _lstOutputConverter
0264 (trackingPhase2PU140 & trackingLST).toReplaceWith(highPtTripletStepTrackCandidates, _lstOutputConverter.clone())
0265
0266 #For FastSim phase1 tracking
0267 import FastSimulation.Tracking.TrackCandidateProducer_cfi
0268 _fastSim_highPtTripletStepTrackCandidates = FastSimulation.Tracking.TrackCandidateProducer_cfi.trackCandidateProducer.clone(
0269     src                      = 'highPtTripletStepSeeds',
0270     MinNumberOfCrossedLayers = 3,
0271     hitMasks                 = cms.InputTag('highPtTripletStepMasks')
0272 )
0273 fastSim.toReplaceWith(highPtTripletStepTrackCandidates,_fastSim_highPtTripletStepTrackCandidates)
0274
0275 # TRACK FITTING
0276 import RecoTracker.TrackProducer.TrackProducerIterativeDefault_cfi
0277 highPtTripletStepTracks = RecoTracker.TrackProducer.TrackProducerIterativeDefault_cfi.TrackProducerIterativeDefault.clone(
0278     src           = 'highPtTripletStepTrackCandidates',
0279     AlgorithmName = 'highPtTripletStep',
0280     Fitter        = 'FlexibleKFFittingSmoother',
0281 )
0282 fastSim.toModify(highPtTripletStepTracks,TTRHBuilder = 'WithoutRefit')
0283
0284 from Configuration.Eras.Modifier_phase2_timing_layer_cff import phase2_timing_layer
0285 phase2_timing_layer.toModify(highPtTripletStepTracks, TrajectoryInEvent = True)
0286
0287 highPtTripletStepLSTpTracks = highPtTripletStepTracks.clone(
0288      src           = 'highPtTripletStepTrackCandidates:pTCsLST'
0289 )
0290 highPtTripletStepLSTT5Tracks = highPtTripletStepTracks.clone(
0291      src           = 'highPtTripletStepTrackCandidates:t5TCsLST'
0292 )
0293 _highPtTripletStepTracks_LST = RecoTracker.FinalTrackSelectors.trackListMerger_cfi.trackListMerger.clone(
0294     TrackProducers     = ['highPtTripletStepLSTpTracks',
0295                           'highPtTripletStepLSTT5Tracks'],
0296     hasSelector        = [1,0],
0297     indivShareFrac     = [0.1,0.1],
0298     selectedTrackQuals = ['highPtTripletStepSelector:highPtTripletStep',
0299                           'highPtTripletStepSelectorLSTT5:highPtTripletStepLSTT5'],
0300     copyExtras         = True,
0301     copyMVA            = False,
0302     setsToMerge        = [cms.PSet( tLists=cms.vint32(0,1), pQual=cms.bool(True) )]
0303 )
0304 (trackingPhase2PU140 & trackingLST).toReplaceWith(highPtTripletStepTracks, _highPtTripletStepTracks_LST)
0305
0306 # Final selection
0307 from RecoTracker.FinalTrackSelectors.TrackMVAClassifierPrompt_cfi import *
0308 highPtTripletStep = TrackMVAClassifierPrompt.clone(
0309      mva         = dict(GBRForestLabel = 'MVASelectorHighPtTripletStep_Phase1'),
0310      src         = 'highPtTripletStepTracks',
0311      qualityCuts = [0.2,0.3,0.4]
0312 )
0313 pp_on_AA.toModify(highPtTripletStep,
0314         mva = dict(GBRForestLabel = 'HIMVASelectorHighPtTripletStep_Phase1'),
0315         qualityCuts = [-0.9, -0.3, 0.85],
0316 )
0317
0318 from RecoTracker.FinalTrackSelectors.trackTfClassifier_cfi import *
0319 from RecoTracker.FinalTrackSelectors.trackSelectionTf_cfi import *
0320 from RecoTracker.FinalTrackSelectors.trackSelectionTf_CKF_cfi import *
0321 trackdnn.toReplaceWith(highPtTripletStep, trackTfClassifier.clone(
0322     src = 'highPtTripletStepTracks',
0323     qualityCuts = qualityCutDictionary.HighPtTripletStep.value()
0324 ))
0325 highBetaStar.toModify(highPtTripletStep,qualityCuts = [-0.2,0.3,0.4])
0326
0327 fastSim.toModify(highPtTripletStep,vertices = 'firstStepPrimaryVerticesBeforeMixing')
0328
0329 (pp_on_AA & trackdnn).toModify(highPtTripletStep, qualityCuts =  [-0.5, 0.5, 0.94] )
0330
0331 # For Phase2PU140
0332 import RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi
0333 highPtTripletStepSelector = RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.multiTrackSelector.clone(
0334     src = 'highPtTripletStepTracks',
0335     trackSelectors = [
0336         RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.looseMTS.clone(
0337             name = 'highPtTripletStepLoose',
0338             chi2n_par = 2.0,
0339             res_par = ( 0.003, 0.002 ),
0340             minNumberLayers = 3,
0341             maxNumberLostLayers = 3,
0342             minNumber3DLayers = 3,
0343             d0_par1 = ( 0.7, 4.0 ),
0344             dz_par1 = ( 0.8, 4.0 ),
0345             d0_par2 = ( 0.6, 4.0 ),
0346             dz_par2 = ( 0.6, 4.0 )
0347             ), #end of pset
0348         RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.tightMTS.clone(
0349             name = 'highPtTripletStepTight',
0350             preFilterName = 'highPtTripletStepLoose',
0351             chi2n_par = 1.0,
0352             res_par = ( 0.003, 0.002 ),
0353             minNumberLayers = 3,
0354             maxNumberLostLayers = 2,
0355             minNumber3DLayers = 3,
0356             d0_par1 = ( 0.6, 4.0 ),
0357             dz_par1 = ( 0.7, 4.0 ),
0358             d0_par2 = ( 0.5, 4.0 ),
0359             dz_par2 = ( 0.6, 4.0 )
0360             ),
0361         RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.highpurityMTS.clone(
0362             name = 'highPtTripletStep',
0363             preFilterName = 'highPtTripletStepTight',
0364             min_eta = -4.1,
0365             max_eta = 4.1,
0366             chi2n_par = 0.8,
0367             res_par = ( 0.003, 0.001 ),
0368             min_nhits = 4,
0369             minNumberLayers = 4,
0370             maxNumberLostLayers = 2,
0371             minNumber3DLayers = 4,
0372             d0_par1 = ( 0.6, 4.0 ),
0373             dz_par1 = ( 0.7, 4.0 ),
0374             d0_par2 = ( 0.45, 4.0 ),
0375             dz_par2 = ( 0.55, 4.0 )
0376             ),
0377     ] #end of vpset
0378 ) #end of clone
0379
0380 from Configuration.ProcessModifiers.vectorHits_cff import vectorHits
0381 vectorHits.toModify(highPtTripletStepSelector.trackSelectors[2], minNumberLayers = 3, minNumber3DLayers = 3, d0_par1 = ( 0.5, 4.0 ), dz_par1 = ( 0.6, 4.0 ))
0382
0383 (trackingPhase2PU140 & trackingLST).toModify(highPtTripletStepSelector, src = 'highPtTripletStepLSTpTracks')
0384 # Passthrough selector to satisfy the TrackListMerger requirement for selector values
0385 highPtTripletStepSelectorLSTT5 = RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.multiTrackSelector.clone(
0386     src = 'highPtTripletStepLSTT5Tracks',
0387     trackSelectors = [
0388         RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.looseMTS.clone(
0389             name = 'highPtTripletStepLSTT5Loose',
0390             minHitsToBypassChecks = 0
0391             ), #end of pset
0392         RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.tightMTS.clone(
0393             name = 'highPtTripletStepLSTT5Tight',
0394             preFilterName = 'highPtTripletStepLSTT5Loose',
0395             minHitsToBypassChecks = 0
0396             ),
0397         RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.highpurityMTS.clone(
0398             name = 'highPtTripletStepLSTT5',
0399             preFilterName = 'highPtTripletStepLSTT5Tight',
0400             minHitsToBypassChecks = 0
0401             ),
0402     ] #end of vpset
0403 ) #end of clone
0404
0405 # Final sequence
0406 HighPtTripletStepTask = cms.Task(highPtTripletStepClusters,
0407                                  highPtTripletStepSeedLayers,
0408                                  highPtTripletStepTrackingRegions,
0409                                  highPtTripletStepHitDoublets,
0410                                  highPtTripletStepHitTriplets,
0411                                  highPtTripletStepSeeds,
0412                                  highPtTripletStepTrackCandidates,
0413                                  highPtTripletStepTracks,
0414                                  highPtTripletStep)
0415 HighPtTripletStep = cms.Sequence(HighPtTripletStepTask)
0416
0417 _HighPtTripletStepTask_trackingMkFit = HighPtTripletStepTask.copy()
0418 _HighPtTripletStepTask_trackingMkFit.add(highPtTripletStepTrackCandidatesMkFitSeeds, highPtTripletStepTrackCandidatesMkFit, highPtTripletStepTrackCandidatesMkFitConfig)
0419 trackingMkFitHighPtTripletStep.toReplaceWith(HighPtTripletStepTask, _HighPtTripletStepTask_trackingMkFit)
0420
0421 _HighPtTripletStepTask_Phase2PU140 = HighPtTripletStepTask.copy()
0422 _HighPtTripletStepTask_Phase2PU140.replace(highPtTripletStep, highPtTripletStepSelector)
0423 _HighPtTripletStep_Phase2PU140 = cms.Sequence(_HighPtTripletStepTask_Phase2PU140)
0424 trackingPhase2PU140.toReplaceWith(HighPtTripletStepTask, _HighPtTripletStepTask_Phase2PU140)
0425
0426 _HighPtTripletStepTask_LST = HighPtTripletStepTask.copy()
0427 from RecoLocalTracker.Phase2TrackerRecHits.Phase2TrackerRecHits_cfi import siPhase2RecHits
0428 from RecoTracker.LST.lstSeedTracks_cff import lstInitialStepSeedTracks,lstHighPtTripletStepSeedTracks
0429 from RecoTracker.LST.lstPixelSeedInputProducer_cfi import lstPixelSeedInputProducer
0430 from RecoTracker.LST.lstPhase2OTHitsInputProducer_cfi import lstPhase2OTHitsInputProducer
0431 from RecoTracker.LST.lstProducerTask_cff import *
0432
0433 _HighPtTripletStepTask_LST.add(siPhase2RecHits, lstInitialStepSeedTracks, lstHighPtTripletStepSeedTracks, lstPixelSeedInputProducer, lstPhase2OTHitsInputProducer,
0434                                lstProducerTask, highPtTripletStepLSTpTracks, highPtTripletStepLSTT5Tracks, highPtTripletStepSelectorLSTT5)
0435 (trackingPhase2PU140 & trackingLST).toReplaceWith(HighPtTripletStepTask, _HighPtTripletStepTask_LST)
0436
0437 # fast tracking mask producer
0438 from FastSimulation.Tracking.FastTrackerRecHitMaskProducer_cfi import maskProducerFromClusterRemover
0439 highPtTripletStepMasks = maskProducerFromClusterRemover(highPtTripletStepClusters)
0440 fastSim.toReplaceWith(HighPtTripletStepTask,
0441                       cms.Task(highPtTripletStepMasks
0442                                ,highPtTripletStepTrackingRegions
0443                                ,highPtTripletStepSeeds
0444                                ,highPtTripletStepTrackCandidates
0445                                ,highPtTripletStepTracks
0446                                ,highPtTripletStep
0447                                ) )