electronIsoFromDeps_cff.py

Go to the documentation of this file.

# computing of isolation from deposits,
# as provided by Matthew LeBourgeois

from RecoEgamma.EgammaIsolationAlgos.egammaIsolationSequence_cff import *
from RecoEgamma.EgammaIsolationAlgos.egammaIsolationSequencePAT_cff import *

# create both IsoDeposits from the full collection and the reduced collection

eleIsoDepositEcalFromHitsFull = eleIsoDepositEcalFromHits.clone()
eleIsoDepositEcalFromHitsReduced = eleIsoDepositEcalFromHits.clone()

eleIsoDepositEcalFromHitsFull.ExtractorPSet.barrelEcalHits = cms.InputTag("ecalRecHit","EcalRecHitsEB")
eleIsoDepositEcalFromHitsFull.ExtractorPSet.endcapEcalHits = cms.InputTag("ecalRecHit","EcalRecHitsEE")
eleIsoDepositEcalFromHitsReduced.ExtractorPSet.barrelEcalHits = cms.InputTag("reducedEcalRecHitsEB")
eleIsoDepositEcalFromHitsReduced.ExtractorPSet.endcapEcalHits = cms.InputTag("reducedEcalRecHitsEE")

# clone the value map producers for each DR

eleIsoFromDepsTk03                           = eleIsoFromDepsTk.clone()
eleIsoFromDepsTk04                           = eleIsoFromDepsTk.clone()
eleIsoFromDepsEcalFromHitsByCrystalFull03    = eleIsoFromDepsEcalFromHitsByCrystal.clone()
eleIsoFromDepsEcalFromHitsByCrystalFull04    = eleIsoFromDepsEcalFromHitsByCrystal.clone()
eleIsoFromDepsEcalFromHitsByCrystalReduced03 = eleIsoFromDepsEcalFromHitsByCrystal.clone()
eleIsoFromDepsEcalFromHitsByCrystalReduced04 = eleIsoFromDepsEcalFromHitsByCrystal.clone()
eleIsoFromDepsHcalFromTowers03               = eleIsoFromDepsHcalFromTowers.clone()
eleIsoFromDepsHcalFromTowers04               = eleIsoFromDepsHcalFromTowers.clone()

# set the correct delta R
eleIsoFromDepsTk03.deposits[0].deltaR                         = 0.3
eleIsoFromDepsTk04.deposits[0].deltaR                         = 0.4
eleIsoFromDepsEcalFromHitsByCrystalFull03.deposits[0].deltaR  = 0.3
eleIsoFromDepsEcalFromHitsByCrystalFull04.deposits[0].deltaR  = 0.4
eleIsoFromDepsEcalFromHitsByCrystalReduced03.deposits[0].deltaR  = 0.3
eleIsoFromDepsEcalFromHitsByCrystalReduced04.deposits[0].deltaR  = 0.4
eleIsoFromDepsHcalFromTowers03.deposits[0].deltaR             = 0.3
eleIsoFromDepsHcalFromTowers04.deposits[0].deltaR             = 0.4

# change the source on the ECAL hits to make sure to get the reduced or full collection

eleIsoFromDepsEcalFromHitsByCrystalFull03.deposits[0].src = "eleIsoDepositEcalFromHitsFull"
eleIsoFromDepsEcalFromHitsByCrystalFull04.deposits[0].src = "eleIsoDepositEcalFromHitsFull"
eleIsoFromDepsEcalFromHitsByCrystalReduced03.deposits[0].src = "eleIsoDepositEcalFromHitsReduced"
eleIsoFromDepsEcalFromHitsByCrystalReduced04.deposits[0].src = "eleIsoDepositEcalFromHitsReduced"

# the sequence

electronIsoFromDeps = cms.Sequence(
  eleIsoDepositTk*eleIsoDepositEcalFromHitsFull*
  eleIsoDepositEcalFromHitsReduced*
  eleIsoDepositHcalFromTowers*      
  eleIsoFromDepsTk03*
  eleIsoFromDepsTk04*
  eleIsoFromDepsEcalFromHitsByCrystalFull03*
  eleIsoFromDepsEcalFromHitsByCrystalFull04*
  eleIsoFromDepsEcalFromHitsByCrystalReduced03*
  eleIsoFromDepsEcalFromHitsByCrystalReduced04*
  eleIsoFromDepsHcalFromTowers03*
  eleIsoFromDepsHcalFromTowers04
)