IterativeTobTecStep_cff.py

Go to the documentation of this file.

import FWCore.ParameterSet.Config as cms

# step 5

# seeding
#from FastSimulation.Tracking.IterativeFifthSeedProducer_cff import *
import FastSimulation.Tracking.TrajectorySeedProducer_cfi
iterativeTobTecSeeds = FastSimulation.Tracking.TrajectorySeedProducer_cfi.trajectorySeedProducer.clone()
iterativeTobTecSeeds.firstHitSubDetectorNumber = [2]
iterativeTobTecSeeds.firstHitSubDetectors = [5, 6]
iterativeTobTecSeeds.secondHitSubDetectorNumber = [2]
iterativeTobTecSeeds.secondHitSubDetectors = [5, 6]
iterativeTobTecSeeds.thirdHitSubDetectorNumber = [0]
iterativeTobTecSeeds.thirdHitSubDetectors = []
iterativeTobTecSeeds.seedingAlgo = ['TobTecLayerPairs']
iterativeTobTecSeeds.minRecHits = [4]
iterativeTobTecSeeds.pTMin = [0.3]
#cut on fastsim simtracks. I think it should be removed for the 5th step
iterativeTobTecSeeds.maxD0 = [99.]
iterativeTobTecSeeds.maxZ0 = [99.]
#-----
iterativeTobTecSeeds.numberOfHits = [2]
#values for the seed compatibility constraint
iterativeTobTecSeeds.originRadius = [6.0] # was 5.0
iterativeTobTecSeeds.originHalfLength = [30.0] # was 10.0
iterativeTobTecSeeds.originpTMin = [0.6] # was 0.5
iterativeTobTecSeeds.zVertexConstraint = [-1.0]
iterativeTobTecSeeds.primaryVertices = ['none']

# candidate producer
#from FastSimulation.Tracking.IterativeFifthCandidateProducer_cff import *
import FastSimulation.Tracking.TrackCandidateProducer_cfi
iterativeTobTecTrackCandidates = FastSimulation.Tracking.TrackCandidateProducer_cfi.trackCandidateProducer.clone()
iterativeTobTecTrackCandidates.SeedProducer = cms.InputTag("iterativeTobTecSeeds","TobTecLayerPairs")
iterativeTobTecTrackCandidates.TrackProducers = ['pixelPairStepTracks','detachedTripletStepTracks','mixedTripletStepTracks','pixelLessStepTracks'] # add 0 and 0.5?
iterativeTobTecTrackCandidates.KeepFittedTracks = False
iterativeTobTecTrackCandidates.MinNumberOfCrossedLayers = 3


# track producer
#from FastSimulation.Tracking.IterativeFifthTrackProducer_cff import *
import RecoTracker.TrackProducer.CTFFinalFitWithMaterial_cfi
iterativeTobTecTracks = RecoTracker.TrackProducer.CTFFinalFitWithMaterial_cfi.ctfWithMaterialTracks.clone()
iterativeTobTecTracks.src = 'iterativeTobTecTrackCandidates'
iterativeTobTecTracks.TTRHBuilder = 'WithoutRefit'
iterativeTobTecTracks.Fitter = 'KFFittingSmootherFifth'
iterativeTobTecTracks.Propagator = 'PropagatorWithMaterial'


# track merger
#from FastSimulation.Tracking.IterativeFifthTrackMerger_cfi import *
tobTecStepTracks = cms.EDProducer("FastTrackMerger",
                                  TrackProducers = cms.VInputTag(cms.InputTag("iterativeTobTecTrackCandidates"),
                                                                 cms.InputTag("iterativeTobTecTracks")),
                                  RemoveTrackProducers =  cms.untracked.VInputTag(cms.InputTag("initialStepTracks"),
                                                                                  cms.InputTag("lowPtTripletStepTracks"),  
                                                                                  cms.InputTag("pixelPairStepTracks"),  
                                                                                  cms.InputTag("detachedTripletStepTracks"),    
                                                                                  cms.InputTag("mixedTripletStepTracks"),     
                                                                                  cms.InputTag("pixelLessStepTracks")),   
                                  trackAlgo = cms.untracked.uint32(10), # iter6
                                  MinNumberOfTrajHits = cms.untracked.uint32(6), # was 4
                                  MaxLostTrajHits = cms.untracked.uint32(0)
                                  )


# track selection
import RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi
tobTecStepSelector = RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.multiTrackSelector.clone(
        src='tobTecStepTracks',
            trackSelectors= cms.VPSet(
            RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.looseMTS.clone(
                name = 'tobTecStepLoose',
                            chi2n_par = 0.4,
                            res_par = ( 0.003, 0.001 ),
                            minNumberLayers = 5,
                            maxNumberLostLayers = 1,
                            minNumber3DLayers = 2,
                            d0_par1 = ( 2.0, 4.0 ),
                            dz_par1 = ( 1.8, 4.0 ),
                            d0_par2 = ( 2.0, 4.0 ),
                            dz_par2 = ( 1.8, 4.0 )
                            ),
                    RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.tightMTS.clone(
                name = 'tobTecStepTight',
                            preFilterName = 'tobTecStepLoose',
                            chi2n_par = 0.3,
                            res_par = ( 0.003, 0.001 ),
                            minNumberLayers = 5,
                            maxNumberLostLayers = 0,
                            minNumber3DLayers = 2,
                            d0_par1 = ( 1.5, 4.0 ),
                            dz_par1 = ( 1.4, 4.0 ),
                            d0_par2 = ( 1.5, 4.0 ),
                            dz_par2 = ( 1.4, 4.0 )
                            ),
                    RecoTracker.FinalTrackSelectors.multiTrackSelector_cfi.highpurityMTS.clone(
                name = 'tobTecStep',
                            preFilterName = 'tobTecStepTight',
                            chi2n_par = 0.2,
                            res_par = ( 0.003, 0.001 ),
                            minNumberLayers = 5,
                            maxNumberLostLayers = 0,
                            minNumber3DLayers = 2,
                            d0_par1 = ( 1.4, 4.0 ),
                            dz_par1 = ( 1.3, 4.0 ),
                            d0_par2 = ( 1.4, 4.0 ),
                            dz_par2 = ( 1.3, 4.0 )
                            ),
                    ) #end of vpset
            ) #end of clone

# sequence
iterativeTobTecStep = cms.Sequence(iterativeTobTecSeeds
                                      +iterativeTobTecTrackCandidates
                                      +iterativeTobTecTracks
                                      +tobTecStepTracks
                                      +tobTecStepSelector)