1 |
# $Id: PixelLessStepGsf_cff.py,v 1.3 2009/03/20 17:18:39 loizides Exp $
|
2 |
|
3 |
import FWCore.ParameterSet.Config as cms
|
4 |
|
5 |
#HIT REMOVAL
|
6 |
fourthClustersGsf = cms.EDFilter("TrackClusterRemover",
|
7 |
oldClusterRemovalInfo = cms.InputTag("thClustersGsf"),
|
8 |
trajectories = cms.InputTag("thWithMaterialTracksGsf"),
|
9 |
pixelClusters = cms.InputTag("thClustersGsf"),
|
10 |
Common = cms.PSet(
|
11 |
maxChi2 = cms.double(30.0)
|
12 |
),
|
13 |
stripClusters = cms.InputTag("thClustersGsf")
|
14 |
)
|
15 |
|
16 |
#TRACKER HITS
|
17 |
import RecoLocalTracker.SiPixelRecHits.SiPixelRecHits_cfi
|
18 |
import RecoLocalTracker.SiStripRecHitConverter.SiStripRecHitConverter_cfi
|
19 |
fourthPixelRecHitsGsf = RecoLocalTracker.SiPixelRecHits.SiPixelRecHits_cfi.siPixelRecHits.clone()
|
20 |
fourthStripRecHitsGsf = RecoLocalTracker.SiStripRecHitConverter.SiStripRecHitConverter_cfi.siStripMatchedRecHits.clone()
|
21 |
fourthPixelRecHitsGsf.src = 'fourthClustersGsf'
|
22 |
fourthStripRecHitsGsf.ClusterProducer = 'fourthClustersGsf'
|
23 |
|
24 |
#SEEDS
|
25 |
import RecoTracker.TkSeedGenerator.GlobalMixedSeeds_cfi
|
26 |
fourthPLSeedsGsf = RecoTracker.TkSeedGenerator.GlobalMixedSeeds_cfi.globalMixedSeeds.clone()
|
27 |
import RecoTracker.MeasurementDet.MeasurementTrackerESProducer_cfi
|
28 |
fourthPLSeedsGsf.OrderedHitsFactoryPSet.SeedingLayers = 'FourthLayerPairsGsf'
|
29 |
fourthPLSeedsGsf.RegionFactoryPSet.RegionPSet.ptMin = 0.6
|
30 |
fourthPLSeedsGsf.RegionFactoryPSet.RegionPSet.originHalfLength = 10.0
|
31 |
fourthPLSeedsGsf.RegionFactoryPSet.RegionPSet.originRadius = 2.0
|
32 |
|
33 |
#TRAJECTORY MEASUREMENT
|
34 |
fourthMeasurementTrackerGsf = RecoTracker.MeasurementDet.MeasurementTrackerESProducer_cfi.MeasurementTracker.clone()
|
35 |
import TrackingTools.TrajectoryFiltering.TrajectoryFilterESProducer_cfi
|
36 |
fourthMeasurementTrackerGsf.ComponentName = 'fourthMeasurementTrackerGsf'
|
37 |
fourthMeasurementTrackerGsf.pixelClusterProducer = 'fourthClustersGsf'
|
38 |
fourthMeasurementTrackerGsf.stripClusterProducer = 'fourthClustersGsf'
|
39 |
|
40 |
#TRAJECTORY FILTER
|
41 |
fourthCkfTrajectoryFilterGsf = TrackingTools.TrajectoryFiltering.TrajectoryFilterESProducer_cfi.trajectoryFilterESProducer.clone()
|
42 |
import RecoTracker.CkfPattern.GroupedCkfTrajectoryBuilderESProducer_cfi
|
43 |
fourthCkfTrajectoryFilterGsf.ComponentName = 'fourthCkfTrajectoryFilterGsf'
|
44 |
fourthCkfTrajectoryFilterGsf.filterPset.maxLostHits = 0
|
45 |
fourthCkfTrajectoryFilterGsf.filterPset.minimumNumberOfHits = 5
|
46 |
fourthCkfTrajectoryFilterGsf.filterPset.minPt = 0.3
|
47 |
|
48 |
# #TRAJECTORY BUILDER
|
49 |
# fourthCkfTrajectoryBuilderGsf = RecoTracker.CkfPattern.GroupedCkfTrajectoryBuilderESProducer_cfi.GroupedCkfTrajectoryBuilder.clone()
|
50 |
# import RecoTracker.CkfPattern.CkfTrackCandidates_cfi
|
51 |
# fourthCkfTrajectoryBuilderGsf.ComponentName = 'fourthCkfTrajectoryBuilderGsf'
|
52 |
# fourthCkfTrajectoryBuilderGsf.MeasurementTrackerName = 'fourthMeasurementTrackerGsf'
|
53 |
# fourthCkfTrajectoryBuilderGsf.trajectoryFilterName = 'fourthCkfTrajectoryFilterGsf'
|
54 |
#
|
55 |
#
|
56 |
# #TRACK CANDIDATES
|
57 |
# fourthTrackCandidatesGsf = RecoTracker.CkfPattern.CkfTrackCandidates_cfi.ckfTrackCandidates.clone()
|
58 |
# import RecoTracker.TrackProducer.CTFFinalFitWithMaterial_cfi
|
59 |
# fourthTrackCandidatesGsf.SeedProducer = 'fourthPLSeedsGsf'
|
60 |
# fourthTrackCandidatesGsf.TrajectoryBuilder = 'fourthCkfTrajectoryBuilderGsf'
|
61 |
|
62 |
import RecoEgamma.EgammaElectronProducers.gsfElectronCkfTrackCandidateMaker_cff
|
63 |
#TRAJECTORY BUILDER
|
64 |
fourthCkfTrajectoryBuilderGsf = RecoEgamma.EgammaElectronProducers.gsfElectronCkfTrackCandidateMaker_cff.TrajectoryBuilderForPixelMatchGsfElectrons.clone()
|
65 |
fourthCkfTrajectoryBuilderGsf.ComponentName = 'fourthCkfTrajectoryBuilderGsf'
|
66 |
fourthCkfTrajectoryBuilderGsf.MeasurementTrackerName = 'fourthMeasurementTrackerGsf'
|
67 |
fourthCkfTrajectoryBuilderGsf.trajectoryFilterName = 'fourthCkfTrajectoryFilterGsf'
|
68 |
|
69 |
|
70 |
#TRACK CANDIDATES
|
71 |
fourthTrackCandidatesGsf = RecoEgamma.EgammaElectronProducers.gsfElectronCkfTrackCandidateMaker_cff.egammaCkfTrackCandidates.clone()
|
72 |
fourthTrackCandidatesGsf.SeedProducer = 'fourthPLSeedsGsf'
|
73 |
fourthTrackCandidatesGsf.TrajectoryBuilder = 'fourthCkfTrajectoryBuilderGsf'
|
74 |
|
75 |
|
76 |
import TrackingTools.GsfTracking.GsfElectronFit_cfi
|
77 |
#TRACKS
|
78 |
fourthWithMaterialTracksGsf = TrackingTools.GsfTracking.GsfElectronFit_cfi.GsfGlobalElectronTest.clone()
|
79 |
fourthWithMaterialTracksGsf.src = 'fourthTrackCandidatesGsf'
|
80 |
fourthWithMaterialTracksGsf.clusterRemovalInfo = cms.InputTag('fourthClustersGsf')
|
81 |
fourthWithMaterialTracksGsf.AlgorithmName = cms.string('iter4')
|
82 |
fourthWithMaterialTracksGsf.TrajectoryInEvent = True
|
83 |
|
84 |
#SEEDING LAYERS
|
85 |
fourthlayerpairsGsf = cms.ESProducer("PixelLessLayerPairsESProducer",
|
86 |
ComponentName = cms.string('FourthLayerPairsGsf'),
|
87 |
layerList = cms.vstring('TIB1+TIB2',
|
88 |
'TIB1+TID1_pos','TIB1+TID1_neg',
|
89 |
'TID1_pos+TID2_pos','TID2_pos+TID3_pos','TID3_pos+TEC1_pos',
|
90 |
'TEC1_pos+TEC2_pos','TEC2_pos+TEC3_pos','TEC3_pos+TEC4_pos','TEC3_pos+TEC5_pos','TEC4_pos+TEC5_pos',
|
91 |
'TID1_neg+TID2_neg','TID2_neg+TID3_neg','TID3_neg+TEC1_neg',
|
92 |
'TEC1_neg+TEC2_neg','TEC2_neg+TEC3_neg','TEC3_neg+TEC4_neg','TEC3_neg+TEC5_neg','TEC4_neg+TEC5_neg'),
|
93 |
TIB = cms.PSet(
|
94 |
TTRHBuilder = cms.string('WithTrackAngle'),
|
95 |
matchedRecHits = cms.InputTag("fourthStripRecHitsGsf","matchedRecHit")
|
96 |
),
|
97 |
TID = cms.PSet(
|
98 |
matchedRecHits = cms.InputTag("fourthStripRecHitsGsf","matchedRecHit"),
|
99 |
useRingSlector = cms.untracked.bool(True),
|
100 |
TTRHBuilder = cms.string('WithTrackAngle'),
|
101 |
minRing = cms.int32(1),
|
102 |
maxRing = cms.int32(2)
|
103 |
),
|
104 |
TEC = cms.PSet(
|
105 |
matchedRecHits = cms.InputTag("fourthStripRecHitsGsf","matchedRecHit"),
|
106 |
useRingSlector = cms.untracked.bool(True),
|
107 |
TTRHBuilder = cms.string('WithTrackAngle'),
|
108 |
minRing = cms.int32(1),
|
109 |
maxRing = cms.int32(2)
|
110 |
)
|
111 |
)
|
112 |
|
113 |
#import RecoTracker.FinalTrackSelectors.selectHighPurity_cfi
|
114 |
#pixellessStep = RecoTracker.FinalTrackSelectors.selectHighPurity_cfi.selectHighPurity.clone()
|
115 |
#pixellessStep.src = 'fourthWithMaterialTracks'
|
116 |
#pixellessStep.copyTrajectories = True
|
117 |
#pixellessStep.chi2n_par = 0.3
|
118 |
#pixellessStep.res_par = ( 0.003, 0.001 )
|
119 |
#pixellessStep.minNumberLayers = 5
|
120 |
#pixellessStep.d0_par1 = ( 1.0, 4.0 )
|
121 |
#pixellessStep.dz_par1 = ( 1.0, 4.0 )
|
122 |
#pixellessStep.d0_par2 = ( 1.0, 4.0 )
|
123 |
#pixellessStep.dz_par2 = ( 1.0, 4.0 )
|
124 |
|
125 |
fourthStepGsf = cms.Sequence(fourthClustersGsf*
|
126 |
fourthPixelRecHitsGsf*fourthStripRecHitsGsf*
|
127 |
fourthPLSeedsGsf*
|
128 |
fourthTrackCandidatesGsf*
|
129 |
fourthWithMaterialTracksGsf)
|