| 1 |
## $Id: generalGsfTracking_cfi.py,v 1.4 2009/03/21 12:46:57 loizides Exp $
|
| 2 |
|
| 3 |
#import FWCore.ParameterSet.Config as cms
|
| 4 |
|
| 5 |
#import RecoEgamma.EgammaElectronProducers.gsfElectronCkfTrackCandidateMaker_cff
|
| 6 |
|
| 7 |
#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)
|
