TrackerElectrons/python/ThStepGsf_cff.py

# $Id: ThStepGsf_cff.py,v 1.4 2009/03/21 12:46:57 loizides Exp $

import FWCore.ParameterSet.Config as cms

import RecoLocalTracker.SiPixelRecHits.SiPixelRecHits_cfi
#TRACKER HITS
thPixelRecHitsGsf = RecoLocalTracker.SiPixelRecHits.SiPixelRecHits_cfi.siPixelRecHits.clone()
import RecoLocalTracker.SiStripRecHitConverter.SiStripRecHitConverter_cfi
thStripRecHitsGsf = RecoLocalTracker.SiStripRecHitConverter.SiStripRecHitConverter_cfi.siStripMatchedRecHits.clone()
thPixelRecHitsGsf.src = 'thClustersGsf'
thStripRecHitsGsf.ClusterProducer = 'thClustersGsf'

import RecoTracker.TkSeedGenerator.GlobalMixedSeeds_cff
#SEEDS
thPLSeedsGsf = RecoTracker.TkSeedGenerator.GlobalMixedSeeds_cff.globalMixedSeeds.clone()
import RecoTracker.MeasurementDet.MeasurementTrackerESProducer_cfi
thPLSeedsGsf.OrderedHitsFactoryPSet.SeedingLayers = 'ThLayerPairsGsf'
thPLSeedsGsf.RegionFactoryPSet.RegionPSet.ptMin = 0.6
thPLSeedsGsf.RegionFactoryPSet.RegionPSet.originHalfLength = 7.0
thPLSeedsGsf.RegionFactoryPSet.RegionPSet.originRadius = 0.7

#TRAJECTORY MEASUREMENT
thMeasurementTrackerGsf = RecoTracker.MeasurementDet.MeasurementTrackerESProducer_cfi.MeasurementTracker.clone()
import TrackingTools.TrajectoryFiltering.TrajectoryFilterESProducer_cfi
thMeasurementTrackerGsf.ComponentName = 'thMeasurementTrackerGsf'
thMeasurementTrackerGsf.pixelClusterProducer = 'thClustersGsf'
thMeasurementTrackerGsf.stripClusterProducer = 'thClustersGsf'

#TRAJECTORY FILTER
thCkfTrajectoryFilterGsf = TrackingTools.TrajectoryFiltering.TrajectoryFilterESProducer_cfi.trajectoryFilterESProducer.clone()
import RecoTracker.CkfPattern.GroupedCkfTrajectoryBuilderESProducer_cfi
thCkfTrajectoryFilterGsf.ComponentName = 'thCkfTrajectoryFilterGsf'
thCkfTrajectoryFilterGsf.filterPset.maxLostHits = 0
thCkfTrajectoryFilterGsf.filterPset.minimumNumberOfHits = 3
thCkfTrajectoryFilterGsf.filterPset.minPt = 0.3

# #TRAJECTORY BUILDER
# thCkfTrajectoryBuilderGsf = RecoTracker.CkfPattern.GroupedCkfTrajectoryBuilderESProducer_cfi.GroupedCkfTrajectoryBuilder.clone()
# import RecoTracker.CkfPattern.CkfTrackCandidates_cfi
# thCkfTrajectoryBuilderGsf.ComponentName = 'thCkfTrajectoryBuilderGsf'
# thCkfTrajectoryBuilderGsf.MeasurementTrackerName = 'thMeasurementTrackerGsf'
# thCkfTrajectoryBuilderGsf.trajectoryFilterName = 'thCkfTrajectoryFilterGsf'
# 
# 
# #TRACK CANDIDATES
# thTrackCandidatesGsf = RecoTracker.CkfPattern.CkfTrackCandidates_cfi.ckfTrackCandidates.clone()
# import RecoTracker.TrackProducer.CTFFinalFitWithMaterial_cfi
# thTrackCandidatesGsf.SeedProducer = 'thPLSeedsGsf'
# thTrackCandidatesGsf.TrajectoryBuilder = 'thCkfTrajectoryBuilderGsf'
# thTrackCandidatesGsf.doSeedingRegionRebuilding = True
# thTrackCandidatesGsf.useHitsSplitting = True

import RecoEgamma.EgammaElectronProducers.gsfElectronCkfTrackCandidateMaker_cff
#TRAJECTORY BUILDER
thCkfTrajectoryBuilderGsf = RecoEgamma.EgammaElectronProducers.gsfElectronCkfTrackCandidateMaker_cff.TrajectoryBuilderForPixelMatchGsfElectrons.clone()
thCkfTrajectoryBuilderGsf.ComponentName = 'thCkfTrajectoryBuilderGsf'
thCkfTrajectoryBuilderGsf.MeasurementTrackerName = 'thMeasurementTrackerGsf'
thCkfTrajectoryBuilderGsf.trajectoryFilterName = 'thCkfTrajectoryFilterGsf'


#TRACK CANDIDATES
thTrackCandidatesGsf = RecoEgamma.EgammaElectronProducers.gsfElectronCkfTrackCandidateMaker_cff.egammaCkfTrackCandidates.clone()
#thTrackCandidatesGsf.SeedProducer = 'thPLSeedsGsf'
thTrackCandidatesGsf.TrajectoryBuilder = 'thCkfTrajectoryBuilderGsf'
thTrackCandidatesGsf.doSeedingRegionRebuilding = True
thTrackCandidatesGsf.useHitsSplitting = True

import TrackingTools.GsfTracking.GsfElectronFit_cfi
#TRACKS
thWithMaterialTracksGsf = TrackingTools.GsfTracking.GsfElectronFit_cfi.GsfGlobalElectronTest.clone()
thWithMaterialTracksGsf.src = 'thTrackCandidatesGsf'
thWithMaterialTracksGsf.clusterRemovalInfo = cms.InputTag('thClustersGsf')
thWithMaterialTracksGsf.TrajectoryInEvent = True


#HIT REMOVAL
thClustersGsf = cms.EDFilter("TrackClusterRemover",
    oldClusterRemovalInfo = cms.InputTag("secClustersGsf"),
    trajectories = cms.InputTag("secWithMaterialTracksGsf"),
    pixelClusters = cms.InputTag("secClustersGsf"),
    Common = cms.PSet(
        maxChi2 = cms.double(30.0)
    ),
    stripClusters = cms.InputTag("secClustersGsf")
)

#SEEDING LAYERS
thlayerpairsGsf = cms.ESProducer("MixedLayerPairsESProducer",
    ComponentName = cms.string('ThLayerPairsGsf'),
    layerList = cms.vstring('BPix1+BPix2', 
        'BPix2+BPix3',
        'BPix1+FPix1_pos',
        'BPix1+FPix1_neg',
        'FPix1_pos+FPix2_pos',
        'FPix1_neg+FPix2_neg',
        'FPix2_pos+TEC2_pos',
        'FPix2_neg+TEC2_neg'),
    TEC = cms.PSet(
        matchedRecHits = cms.InputTag("thStripRecHitsGsf","matchedRecHit"),
        useRingSlector = cms.untracked.bool(True),
        TTRHBuilder = cms.string('WithTrackAngle'),
        minRing = cms.int32(1),
        maxRing = cms.int32(2)
    ),
    BPix = cms.PSet(
        useErrorsFromParam = cms.untracked.bool(True),
        hitErrorRPhi = cms.double(0.0027),
        TTRHBuilder = cms.string('TTRHBuilderWithoutAngle4MixedPairs'),
        HitProducer = cms.string('thPixelRecHitsGsf'),
        hitErrorRZ = cms.double(0.006)
    ),
    FPix = cms.PSet(
        useErrorsFromParam = cms.untracked.bool(True),
        hitErrorRPhi = cms.double(0.0051),
        TTRHBuilder = cms.string('TTRHBuilderWithoutAngle4MixedPairs'),
        HitProducer = cms.string('thPixelRecHitsGsf'),
        hitErrorRZ = cms.double(0.0036)
    )
)

##from RecoTracker.IterativeTracking.ThVxFilter_cff import *
#import RecoTracker.FinalTrackSelectors.selectHighPurity_cfi
#thStepVtx = RecoTracker.FinalTrackSelectors.selectHighPurity_cfi.selectHighPurity.clone()
#thStepVtx.src = 'thWithMaterialTracks'
#thStepVtx.copyTrajectories = True
#thStepVtx.chi2n_par = 0.9
#thStepVtx.res_par = ( 0.003, 0.001 )
#thStepVtx.d0_par1 = ( 0.9, 3.0 )
#thStepVtx.dz_par1 = ( 0.9, 3.0 )
#thStepVtx.d0_par2 = ( 1.0, 3.0 )
#thStepVtx.dz_par2 = ( 1.0, 3.0 )

#thStepTrk = RecoTracker.FinalTrackSelectors.selectHighPurity_cfi.selectHighPurity.clone()
#thStepTrk.src = 'thWithMaterialTracks'
#thStepTrk.copyTrajectories = True
#thStepTrk.chi2n_par = 0.5
#thStepTrk.res_par = ( 0.003, 0.001 )
#thStepTrk.minNumberLayers = 5
#thStepTrk.d0_par1 = ( 1.0, 4.0 )
#thStepTrk.dz_par1 = ( 1.0, 4.0 )
#thStepTrk.d0_par2 = ( 1.0, 4.0 )
#thStepTrk.dz_par2 = ( 1.0, 4.0 )
#import RecoTracker.FinalTrackSelectors.ctfrsTrackListMerger_cfi

#thStep = RecoTracker.FinalTrackSelectors.ctfrsTrackListMerger_cfi.ctfrsTrackListMerger.clone()
#thStep.TrackProducer1 = 'thStepVtx'
#thStep.TrackProducer2 = 'thStepTrk'

thirdStepGsf = cms.Sequence(thClustersGsf*
                            thPixelRecHitsGsf*thStripRecHitsGsf*
                            thPLSeedsGsf*
                            thTrackCandidatesGsf*
                            thWithMaterialTracksGsf
                           )
Revision:	1.5
Committed:	Fri Jul 17 19:53:50 2009 UTC (15 years, 9 months ago) by bendavid
Content type:	text/x-python
Branch:	MAIN
CVS Tags:	Mit_032, Mit_031, Mit_025c_branch2, Mit_025c_branch1, Mit_030, Mit_029c, Mit_029b, Mit_030_pre1, Mit_029a, Mit_029, Mit_029_pre1, Mit_028a, Mit_025c_branch0, Mit_028, Mit_027a, Mit_027, Mit_026, Mit_025e, Mit_025d, Mit_025c, Mit_025b, Mit_025a, Mit_025, Mit_025pre2, Mit_024b, Mit_025pre1, Mit_024a, Mit_024, Mit_023, Mit_022a, Mit_022, Mit_020d, TMit_020d, Mit_020c, Mit_021, Mit_021pre2, Mit_021pre1, Mit_020b, Mit_020a, Mit_020, Mit_020pre1, Mit_018, Mit_017, Mit_017pre3, Mit_017pre2, Mit_017pre1, V07-05-00, Mit_016, Mit_015b, Mit_015a, Mit_015, Mit_014e, Mit_014d, Mit_014c, Mit_014b, ConvRejection-10-06-09, Mit_014a, Mit_014, Mit_014pre3, Mit_014pre2, Mit_014pre1, Mit_013d, Mit_013c, Mit_013b, Mit_013a, Mit_013, Mit_013pre1, Mit_012i, Mit_012h, Mit_012g, Mit_012f, Mit_012e, Mit_012d, Mit_012c, Mit_012b, Mit_012a, Mit_012, Mit_011a, Mit_011, Mit_010a, Mit_010, HEAD
Branch point for:	Mit_025c_branch
Changes since 1.4:	+4 -4 lines
Log Message:	Silenced python errors, but all the python here needs to be redone from the 31 tracking stuff
#	User	Rev	Content
1	bendavid	1.5	# $Id: ThStepGsf_cff.py,v 1.4 2009/03/21 12:46:57 loizides Exp $
2	loizides	1.3
3	bendavid	1.1	import FWCore.ParameterSet.Config as cms
4
5			import RecoLocalTracker.SiPixelRecHits.SiPixelRecHits_cfi
6			#TRACKER HITS
7			thPixelRecHitsGsf = RecoLocalTracker.SiPixelRecHits.SiPixelRecHits_cfi.siPixelRecHits.clone()
8			import RecoLocalTracker.SiStripRecHitConverter.SiStripRecHitConverter_cfi
9			thStripRecHitsGsf = RecoLocalTracker.SiStripRecHitConverter.SiStripRecHitConverter_cfi.siStripMatchedRecHits.clone()
10			thPixelRecHitsGsf.src = 'thClustersGsf'
11			thStripRecHitsGsf.ClusterProducer = 'thClustersGsf'
12
13	bendavid	1.5	import RecoTracker.TkSeedGenerator.GlobalMixedSeeds_cff
14	bendavid	1.1	#SEEDS
15	bendavid	1.5	thPLSeedsGsf = RecoTracker.TkSeedGenerator.GlobalMixedSeeds_cff.globalMixedSeeds.clone()
16	bendavid	1.1	import RecoTracker.MeasurementDet.MeasurementTrackerESProducer_cfi
17			thPLSeedsGsf.OrderedHitsFactoryPSet.SeedingLayers = 'ThLayerPairsGsf'
18			thPLSeedsGsf.RegionFactoryPSet.RegionPSet.ptMin = 0.6
19			thPLSeedsGsf.RegionFactoryPSet.RegionPSet.originHalfLength = 7.0
20			thPLSeedsGsf.RegionFactoryPSet.RegionPSet.originRadius = 0.7
21
22			#TRAJECTORY MEASUREMENT
23			thMeasurementTrackerGsf = RecoTracker.MeasurementDet.MeasurementTrackerESProducer_cfi.MeasurementTracker.clone()
24			import TrackingTools.TrajectoryFiltering.TrajectoryFilterESProducer_cfi
25			thMeasurementTrackerGsf.ComponentName = 'thMeasurementTrackerGsf'
26			thMeasurementTrackerGsf.pixelClusterProducer = 'thClustersGsf'
27			thMeasurementTrackerGsf.stripClusterProducer = 'thClustersGsf'
28
29			#TRAJECTORY FILTER
30			thCkfTrajectoryFilterGsf = TrackingTools.TrajectoryFiltering.TrajectoryFilterESProducer_cfi.trajectoryFilterESProducer.clone()
31			import RecoTracker.CkfPattern.GroupedCkfTrajectoryBuilderESProducer_cfi
32			thCkfTrajectoryFilterGsf.ComponentName = 'thCkfTrajectoryFilterGsf'
33			thCkfTrajectoryFilterGsf.filterPset.maxLostHits = 0
34			thCkfTrajectoryFilterGsf.filterPset.minimumNumberOfHits = 3
35			thCkfTrajectoryFilterGsf.filterPset.minPt = 0.3
36
37	bendavid	1.2	# #TRAJECTORY BUILDER
38			# thCkfTrajectoryBuilderGsf = RecoTracker.CkfPattern.GroupedCkfTrajectoryBuilderESProducer_cfi.GroupedCkfTrajectoryBuilder.clone()
39			# import RecoTracker.CkfPattern.CkfTrackCandidates_cfi
40			# thCkfTrajectoryBuilderGsf.ComponentName = 'thCkfTrajectoryBuilderGsf'
41			# thCkfTrajectoryBuilderGsf.MeasurementTrackerName = 'thMeasurementTrackerGsf'
42			# thCkfTrajectoryBuilderGsf.trajectoryFilterName = 'thCkfTrajectoryFilterGsf'
43			#
44			#
45			# #TRACK CANDIDATES
46			# thTrackCandidatesGsf = RecoTracker.CkfPattern.CkfTrackCandidates_cfi.ckfTrackCandidates.clone()
47			# import RecoTracker.TrackProducer.CTFFinalFitWithMaterial_cfi
48			# thTrackCandidatesGsf.SeedProducer = 'thPLSeedsGsf'
49			# thTrackCandidatesGsf.TrajectoryBuilder = 'thCkfTrajectoryBuilderGsf'
50			# thTrackCandidatesGsf.doSeedingRegionRebuilding = True
51			# thTrackCandidatesGsf.useHitsSplitting = True
52
53			import RecoEgamma.EgammaElectronProducers.gsfElectronCkfTrackCandidateMaker_cff
54	bendavid	1.1	#TRAJECTORY BUILDER
55	bendavid	1.2	thCkfTrajectoryBuilderGsf = RecoEgamma.EgammaElectronProducers.gsfElectronCkfTrackCandidateMaker_cff.TrajectoryBuilderForPixelMatchGsfElectrons.clone()
56	bendavid	1.1	thCkfTrajectoryBuilderGsf.ComponentName = 'thCkfTrajectoryBuilderGsf'
57			thCkfTrajectoryBuilderGsf.MeasurementTrackerName = 'thMeasurementTrackerGsf'
58			thCkfTrajectoryBuilderGsf.trajectoryFilterName = 'thCkfTrajectoryFilterGsf'
59
60
61			#TRACK CANDIDATES
62	bendavid	1.2	thTrackCandidatesGsf = RecoEgamma.EgammaElectronProducers.gsfElectronCkfTrackCandidateMaker_cff.egammaCkfTrackCandidates.clone()
63	bendavid	1.5	#thTrackCandidatesGsf.SeedProducer = 'thPLSeedsGsf'
64	bendavid	1.1	thTrackCandidatesGsf.TrajectoryBuilder = 'thCkfTrajectoryBuilderGsf'
65			thTrackCandidatesGsf.doSeedingRegionRebuilding = True
66			thTrackCandidatesGsf.useHitsSplitting = True
67
68			import TrackingTools.GsfTracking.GsfElectronFit_cfi
69			#TRACKS
70			thWithMaterialTracksGsf = TrackingTools.GsfTracking.GsfElectronFit_cfi.GsfGlobalElectronTest.clone()
71			thWithMaterialTracksGsf.src = 'thTrackCandidatesGsf'
72			thWithMaterialTracksGsf.clusterRemovalInfo = cms.InputTag('thClustersGsf')
73			thWithMaterialTracksGsf.TrajectoryInEvent = True
74
75
76			#HIT REMOVAL
77			thClustersGsf = cms.EDFilter("TrackClusterRemover",
78			oldClusterRemovalInfo = cms.InputTag("secClustersGsf"),
79			trajectories = cms.InputTag("secWithMaterialTracksGsf"),
80			pixelClusters = cms.InputTag("secClustersGsf"),
81			Common = cms.PSet(
82			maxChi2 = cms.double(30.0)
83			),
84			stripClusters = cms.InputTag("secClustersGsf")
85			)
86
87			#SEEDING LAYERS
88			thlayerpairsGsf = cms.ESProducer("MixedLayerPairsESProducer",
89			ComponentName = cms.string('ThLayerPairsGsf'),
90			layerList = cms.vstring('BPix1+BPix2',
91			'BPix2+BPix3',
92			'BPix1+FPix1_pos',
93			'BPix1+FPix1_neg',
94			'FPix1_pos+FPix2_pos',
95			'FPix1_neg+FPix2_neg',
96			'FPix2_pos+TEC2_pos',
97			'FPix2_neg+TEC2_neg'),
98			TEC = cms.PSet(
99			matchedRecHits = cms.InputTag("thStripRecHitsGsf","matchedRecHit"),
100			useRingSlector = cms.untracked.bool(True),
101			TTRHBuilder = cms.string('WithTrackAngle'),
102			minRing = cms.int32(1),
103			maxRing = cms.int32(2)
104			),
105			BPix = cms.PSet(
106			useErrorsFromParam = cms.untracked.bool(True),
107			hitErrorRPhi = cms.double(0.0027),
108			TTRHBuilder = cms.string('TTRHBuilderWithoutAngle4MixedPairs'),
109			HitProducer = cms.string('thPixelRecHitsGsf'),
110			hitErrorRZ = cms.double(0.006)
111			),
112			FPix = cms.PSet(
113			useErrorsFromParam = cms.untracked.bool(True),
114			hitErrorRPhi = cms.double(0.0051),
115			TTRHBuilder = cms.string('TTRHBuilderWithoutAngle4MixedPairs'),
116			HitProducer = cms.string('thPixelRecHitsGsf'),
117			hitErrorRZ = cms.double(0.0036)
118			)
119			)
120
121			##from RecoTracker.IterativeTracking.ThVxFilter_cff import *
122			#import RecoTracker.FinalTrackSelectors.selectHighPurity_cfi
123			#thStepVtx = RecoTracker.FinalTrackSelectors.selectHighPurity_cfi.selectHighPurity.clone()
124			#thStepVtx.src = 'thWithMaterialTracks'
125			#thStepVtx.copyTrajectories = True
126			#thStepVtx.chi2n_par = 0.9
127			#thStepVtx.res_par = ( 0.003, 0.001 )
128			#thStepVtx.d0_par1 = ( 0.9, 3.0 )
129			#thStepVtx.dz_par1 = ( 0.9, 3.0 )
130			#thStepVtx.d0_par2 = ( 1.0, 3.0 )
131			#thStepVtx.dz_par2 = ( 1.0, 3.0 )
132
133			#thStepTrk = RecoTracker.FinalTrackSelectors.selectHighPurity_cfi.selectHighPurity.clone()
134			#thStepTrk.src = 'thWithMaterialTracks'
135			#thStepTrk.copyTrajectories = True
136			#thStepTrk.chi2n_par = 0.5
137			#thStepTrk.res_par = ( 0.003, 0.001 )
138			#thStepTrk.minNumberLayers = 5
139			#thStepTrk.d0_par1 = ( 1.0, 4.0 )
140			#thStepTrk.dz_par1 = ( 1.0, 4.0 )
141			#thStepTrk.d0_par2 = ( 1.0, 4.0 )
142			#thStepTrk.dz_par2 = ( 1.0, 4.0 )
143			#import RecoTracker.FinalTrackSelectors.ctfrsTrackListMerger_cfi
144
145			#thStep = RecoTracker.FinalTrackSelectors.ctfrsTrackListMerger_cfi.ctfrsTrackListMerger.clone()
146			#thStep.TrackProducer1 = 'thStepVtx'
147			#thStep.TrackProducer2 = 'thStepTrk'
148
149			thirdStepGsf = cms.Sequence(thClustersGsf*
150	loizides	1.4	thPixelRecHitsGsfthStripRecHitsGsf
151			thPLSeedsGsf*
152			thTrackCandidatesGsf*
153			thWithMaterialTracksGsf
154			)