RecoEgamma/python/electronIsoFromDeps_cff.py

0001
0002 # computing of isolation from deposits,

0003 # as provided by Matthew LeBourgeois

0004
0005 from RecoEgamma.EgammaIsolationAlgos.egammaIsolationSequence_cff import *
0006 from RecoEgamma.EgammaIsolationAlgos.egammaIsolationSequencePAT_cff import *
0007
0008 # create both IsoDeposits from the full collection and the reduced collection

0009
0010 eleIsoDepositEcalFromHitsFull = eleIsoDepositEcalFromHits.clone()
0011 eleIsoDepositEcalFromHitsReduced = eleIsoDepositEcalFromHits.clone()
0012
0013 eleIsoDepositEcalFromHitsFull.ExtractorPSet.barrelEcalHits = "ecalRecHit : EcalRecHitsEB"
0014 eleIsoDepositEcalFromHitsFull.ExtractorPSet.endcapEcalHits = "ecalRecHit : EcalRecHitsEE"
0015 eleIsoDepositEcalFromHitsReduced.ExtractorPSet.barrelEcalHits = "reducedEcalRecHitsEB"
0016 eleIsoDepositEcalFromHitsReduced.ExtractorPSet.endcapEcalHits = "reducedEcalRecHitsEE"
0017
0018
0019 # clone the value map producers for each DR

0020
0021 eleIsoFromDepsTk03                           = eleIsoFromDepsTk.clone()
0022 eleIsoFromDepsTk04                           = eleIsoFromDepsTk.clone()
0023 eleIsoFromDepsEcalFromHitsByCrystalFull03    = eleIsoFromDepsEcalFromHitsByCrystal.clone()
0024 eleIsoFromDepsEcalFromHitsByCrystalFull04    = eleIsoFromDepsEcalFromHitsByCrystal.clone()
0025 eleIsoFromDepsEcalFromHitsByCrystalReduced03 = eleIsoFromDepsEcalFromHitsByCrystal.clone()
0026 eleIsoFromDepsEcalFromHitsByCrystalReduced04 = eleIsoFromDepsEcalFromHitsByCrystal.clone()
0027 eleIsoFromDepsHcalFromTowers03               = eleIsoFromDepsHcalFromTowers.clone()
0028 eleIsoFromDepsHcalFromTowers04               = eleIsoFromDepsHcalFromTowers.clone()
0029
0030 # set the correct delta R

0031 eleIsoFromDepsTk03.deposits[0].deltaR                         = 0.3
0032 eleIsoFromDepsTk04.deposits[0].deltaR                         = 0.4
0033 eleIsoFromDepsEcalFromHitsByCrystalFull03.deposits[0].deltaR  = 0.3
0034 eleIsoFromDepsEcalFromHitsByCrystalFull04.deposits[0].deltaR  = 0.4
0035 eleIsoFromDepsEcalFromHitsByCrystalReduced03.deposits[0].deltaR  = 0.3
0036 eleIsoFromDepsEcalFromHitsByCrystalReduced04.deposits[0].deltaR  = 0.4
0037 eleIsoFromDepsHcalFromTowers03.deposits[0].deltaR             = 0.3
0038 eleIsoFromDepsHcalFromTowers04.deposits[0].deltaR             = 0.4
0039
0040 # change the source on the ECAL hits to make sure to get the reduced or full collection

0041
0042 eleIsoFromDepsEcalFromHitsByCrystalFull03.deposits[0].src = "eleIsoDepositEcalFromHitsFull"
0043 eleIsoFromDepsEcalFromHitsByCrystalFull04.deposits[0].src = "eleIsoDepositEcalFromHitsFull"
0044 eleIsoFromDepsEcalFromHitsByCrystalReduced03.deposits[0].src = "eleIsoDepositEcalFromHitsReduced"
0045 eleIsoFromDepsEcalFromHitsByCrystalReduced04.deposits[0].src = "eleIsoDepositEcalFromHitsReduced"
0046
0047 # the sequence

0048
0049 electronIsoFromDeps = cms.Sequence(
0050   eleIsoDepositTk*eleIsoDepositEcalFromHitsFull*
0051   eleIsoDepositEcalFromHitsReduced*
0052   eleIsoDepositHcalFromTowers*
0053   eleIsoFromDepsTk03*
0054   eleIsoFromDepsTk04*
0055   eleIsoFromDepsEcalFromHitsByCrystalFull03*
0056   eleIsoFromDepsEcalFromHitsByCrystalFull04*
0057   eleIsoFromDepsEcalFromHitsByCrystalReduced03*
0058   eleIsoFromDepsEcalFromHitsByCrystalReduced04*
0059   eleIsoFromDepsHcalFromTowers03*
0060   eleIsoFromDepsHcalFromTowers04
0061 )