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