1 |
bendavid |
1.5 |
## $Id: generalGsfTracking_cfi.py,v 1.4 2009/03/21 12:46:57 loizides Exp $
|
2 |
loizides |
1.3 |
|
3 |
bendavid |
1.5 |
#import FWCore.ParameterSet.Config as cms
|
4 |
bendavid |
1.1 |
|
5 |
bendavid |
1.5 |
#import RecoEgamma.EgammaElectronProducers.gsfElectronCkfTrackCandidateMaker_cff
|
6 |
bendavid |
1.2 |
|
7 |
bendavid |
1.5 |
#TrajectoryBuilderForGeneralGsfTracks = RecoEgamma.EgammaElectronProducers.gsfElectronCkfTrackCandidateMaker_cff.TrajectoryBuilderForPixelMatchGsfElectrons.clone()
|
8 |
|
|
#TrajectoryBuilderForGeneralGsfTracks.ComponentName = 'TrajectoryBuilderForGeneralGsfTracks'
|
9 |
|
|
#TrajectoryBuilderForGeneralGsfTracks.trajectoryFilterName = 'newTrajectoryFilter'
|
10 |
|
|
## TrajectoryBuilderForGeneralGsfTracks.trajectoryFilterName = 'TrajectoryFilterForGeneralGsfTracks'
|
11 |
|
|
|
12 |
|
|
## TrajectoryFilterForGeneralGsfTracks = RecoEgamma.EgammaElectronProducers.gsfElectronCkfTrackCandidateMaker_cff.TrajectoryFilterForPixelMatchGsfElectrons.clone()
|
13 |
|
|
## TrajectoryFilterForGeneralGsfTracks.ComponentName = 'TrajectoryFilterForGeneralGsfTracks'
|
14 |
|
|
## TrajectoryFilterForGeneralGsfTracks.filterPset.minPt = 0.0
|
15 |
|
|
## TrajectoryFilterForGeneralGsfTracks.filterPset.minPt.nSigmaMinPt = 0.0
|
16 |
|
|
|
17 |
|
|
#newTrackCandidateMakerGsf = RecoEgamma.EgammaElectronProducers.gsfElectronCkfTrackCandidateMaker_cff.egammaCkfTrackCandidates.clone()
|
18 |
|
|
#newTrackCandidateMakerGsf.TrajectoryBuilder = 'TrajectoryBuilderForGeneralGsfTracks'
|
19 |
|
|
##newTrackCandidateMakerGsf.SeedProducer = 'newCombinedSeeds'
|
20 |
|
|
|
21 |
|
|
##normal iterative tracking procedure reimplemented using gsf fits for all tracks
|
22 |
|
|
#import TrackingTools.GsfTracking.GsfElectronFit_cfi
|
23 |
|
|
#preFilterFirstStepTracksGsf = TrackingTools.GsfTracking.GsfElectronFit_cfi.GsfGlobalElectronTest.clone()
|
24 |
|
|
#preFilterFirstStepTracksGsf.src = 'newTrackCandidateMakerGsf'
|
25 |
|
|
#preFilterFirstStepTracksGsf.TrajectoryInEvent = True
|
26 |
|
|
|
27 |
|
|
#from MitEdm.TrackerElectrons.iterativeTkGsf_cff import *
|
28 |
|
|
#from MitEdm.TrackerElectrons.MergeTrackCollectionsGsf_cff import *
|
29 |
|
|
|
30 |
|
|
#import MitEdm.Producers.gsfTrackAssociator_cff
|
31 |
|
|
#generalGsfTrackAssociator = MitEdm.Producers.gsfTrackAssociator_cff.gsfTrackAssociator.clone()
|
32 |
|
|
#generalGsfTrackAssociator.fromTracks = 'generalTracksGsf'
|
33 |
|
|
|
34 |
|
|
#import SimTracker.TrackAssociation.trackingParticleRecoTrackAsssociation_cfi
|
35 |
|
|
#gsfMcMatch = SimTracker.TrackAssociation.trackingParticleRecoTrackAsssociation_cfi.trackingParticleRecoTrackAsssociation.clone()
|
36 |
|
|
#gsfMcMatch.label_tr = 'generalTracksGsf'
|
37 |
|
|
|
38 |
|
|
##configure tracker-seeded electron producer with very loose cuts
|
39 |
|
|
#import RecoEgamma.EgammaElectronProducers.gsfElectrons_cfi
|
40 |
|
|
#globalGeneralGsfElectrons = RecoEgamma.EgammaElectronProducers.gsfElectrons_cfi.gsfElectrons.clone()
|
41 |
|
|
#globalGeneralGsfElectrons.tracks = cms.InputTag("generalTracksGsf")
|
42 |
|
|
#globalGeneralGsfElectrons.minEOverPEndcaps = 0.0
|
43 |
|
|
#globalGeneralGsfElectrons.minEOverPBarrel = 0.0
|
44 |
|
|
#globalGeneralGsfElectrons.maxEOverPEndcaps = 99999.9
|
45 |
|
|
#globalGeneralGsfElectrons.maxEOverPBarrel = 99999.9
|
46 |
|
|
##globalGeneralGsfElectrons.maxHOverE = 99999.9
|
47 |
|
|
##globalGeneralGsfElectrons.maxDeltaPhi = 0.3
|
48 |
|
|
##globalGeneralGsfElectrons.maxDeltaEta = 0.06
|
49 |
|
|
##globalGeneralGsfElectrons.applyAmbResolution = True
|
50 |
|
|
|
51 |
|
|
##run isolation and electron id algorithms on new electron collection
|
52 |
|
|
#import RecoEgamma.EgammaIsolationAlgos.eleIsolationSequence_cff
|
53 |
|
|
#eleIsoDepositTkGsf = RecoEgamma.EgammaIsolationAlgos.eleIsolationSequence_cff.eleIsoDepositTk.clone()
|
54 |
|
|
#eleIsoDepositTkGsf.src = cms.InputTag("globalGeneralGsfElectrons")
|
55 |
|
|
#eleIsoDepositEcalFromHitsGsf = RecoEgamma.EgammaIsolationAlgos.eleIsolationSequence_cff.eleIsoDepositEcalFromHits.clone()
|
56 |
|
|
#eleIsoDepositEcalFromHitsGsf.src = cms.InputTag("globalGeneralGsfElectrons")
|
57 |
|
|
##eleIsoDepositHcalFromHitsGsf = RecoEgamma.EgammaIsolationAlgos.eleIsolationSequence_cff.eleIsoDepositHcalFromHits.clone()
|
58 |
|
|
##eleIsoDepositHcalFromHitsGsf.src = cms.InputTag("globalGeneralGsfElectrons")
|
59 |
|
|
#eleIsoFromDepsTkGsf = RecoEgamma.EgammaIsolationAlgos.eleIsolationSequence_cff.eleIsoFromDepsTk.clone()
|
60 |
|
|
#eleIsoFromDepsTkGsf.deposits[0].src = cms.InputTag("eleIsoDepositEcalFromHitsGsf")
|
61 |
|
|
#eleIsoFromDepsEcalFromHitsGsf = RecoEgamma.EgammaIsolationAlgos.eleIsolationSequence_cff.eleIsoFromDepsEcalFromHits.clone()
|
62 |
|
|
#eleIsoFromDepsEcalFromHitsGsf.deposits[0].src = cms.InputTag("eleIsoDepositEcalFromHitsGsf")
|
63 |
|
|
#eleIsoFromDepsHcalFromHitsGsf = RecoEgamma.EgammaIsolationAlgos.eleIsolationSequence_cff.eleIsoFromDepsHcalFromHits.clone()
|
64 |
|
|
#eleIsoFromDepsHcalFromHitsGsf.deposits[0].src = cms.InputTag("eleIsoDepositHcalFromHitsGsf")
|
65 |
|
|
|
66 |
|
|
#generalGsfElectronsIso = cms.Sequence(eleIsoDepositTkGsf*eleIsoDepositEcalFromHitsGsf*eleIsoDepositHcalFromHitsGsf
|
67 |
|
|
#*eleIsoFromDepsTkGsf*eleIsoFromDepsEcalFromHitsGsf*eleIsoFromDepsHcalFromHitsGsf)
|
68 |
|
|
|
69 |
|
|
#import RecoEgamma.ElectronIdentification.electronIdSequence_cff
|
70 |
|
|
#eidRobustLooseGsf = RecoEgamma.ElectronIdentification.electronIdSequence_cff.eidRobustLoose.clone()
|
71 |
|
|
#eidRobustLooseGsf.src = cms.InputTag("globalGeneralGsfElectrons")
|
72 |
|
|
#eidRobustTightGsf = RecoEgamma.ElectronIdentification.electronIdSequence_cff.eidRobustTight.clone()
|
73 |
|
|
#eidRobustTightGsf.src = cms.InputTag("globalGeneralGsfElectrons")
|
74 |
|
|
#eidLooseGsf = RecoEgamma.ElectronIdentification.electronIdSequence_cff.eidLoose.clone()
|
75 |
|
|
#eidLooseGsf.src = cms.InputTag("globalGeneralGsfElectrons")
|
76 |
|
|
#eidTightGsf = RecoEgamma.ElectronIdentification.electronIdSequence_cff.eidTight.clone()
|
77 |
|
|
#eidTightGsf.src = cms.InputTag("globalGeneralGsfElectrons")
|
78 |
|
|
|
79 |
|
|
#import RecoEgamma.ElectronIdentification.electronIdNeuralNet_cfi
|
80 |
|
|
#eidNeuralNetGsf = RecoEgamma.ElectronIdentification.electronIdNeuralNet_cfi.eidNeuralNet.clone()
|
81 |
|
|
#eidNeuralNetGsf.src = cms.InputTag("globalGeneralGsfElectrons")
|
82 |
|
|
|
83 |
|
|
#import RecoEgamma.ElectronIdentification.electronIdLikelihoodExt_cfi
|
84 |
|
|
#eidLikelihoodExtGsf = RecoEgamma.ElectronIdentification.electronIdLikelihoodExt_cfi.eidLikelihoodExt.clone()
|
85 |
|
|
#eidLikelihoodExtGsf.src = cms.InputTag("globalGeneralGsfElectrons")
|
86 |
|
|
|
87 |
|
|
#generalGsfElectronId = cms.Sequence(eidRobustLooseGsf*eidRobustTightGsf*eidLooseGsf*eidTightGsf*eidNeuralNetGsf*eidLikelihoodExtGsf)
|
88 |
|
|
|
89 |
|
|
#generalGsfTracking = cms.Sequence(newTrackCandidateMakerGsf*preFilterFirstStepTracksGsf*iterTrackingGsf*trackCollectionMergingGsf*
|
90 |
|
|
#gsfMcMatch*generalGsfTrackAssociator*
|
91 |
|
|
#globalGeneralGsfElectrons*generalGsfElectronsIso*generalGsfElectronId)
|