Reference Manual

Go to the documentation of this file.
00001 
00002 # computing of isolation from deposits,
00003 # as provided by Matthew LeBourgeois
00004 
00005 from RecoEgamma.EgammaIsolationAlgos.egammaIsolationSequence_cff import *
00006 from RecoEgamma.EgammaIsolationAlgos.egammaIsolationSequencePAT_cff import *
00007 
00008 # create both IsoDeposits from the full collection and the reduced collection
00009 
00010 eleIsoDepositEcalFromHitsFull = eleIsoDepositEcalFromHits.clone()
00011 eleIsoDepositEcalFromHitsReduced = eleIsoDepositEcalFromHits.clone()
00012 
00013 eleIsoDepositEcalFromHitsFull.ExtractorPSet.barrelEcalHits = cms.InputTag("ecalRecHit","EcalRecHitsEB")
00014 eleIsoDepositEcalFromHitsFull.ExtractorPSet.endcapEcalHits = cms.InputTag("ecalRecHit","EcalRecHitsEE")
00015 eleIsoDepositEcalFromHitsReduced.ExtractorPSet.barrelEcalHits = cms.InputTag("reducedEcalRecHitsEB")
00016 eleIsoDepositEcalFromHitsReduced.ExtractorPSet.endcapEcalHits = cms.InputTag("reducedEcalRecHitsEE")
00017 
00018 # clone the value map producers for each DR
00019 
00020 eleIsoFromDepsTk03                           = eleIsoFromDepsTk.clone()
00021 eleIsoFromDepsTk04                           = eleIsoFromDepsTk.clone()
00022 eleIsoFromDepsEcalFromHitsByCrystalFull03    = eleIsoFromDepsEcalFromHitsByCrystal.clone()
00023 eleIsoFromDepsEcalFromHitsByCrystalFull04    = eleIsoFromDepsEcalFromHitsByCrystal.clone()
00024 eleIsoFromDepsEcalFromHitsByCrystalReduced03 = eleIsoFromDepsEcalFromHitsByCrystal.clone()
00025 eleIsoFromDepsEcalFromHitsByCrystalReduced04 = eleIsoFromDepsEcalFromHitsByCrystal.clone()
00026 eleIsoFromDepsHcalFromTowers03               = eleIsoFromDepsHcalFromTowers.clone()
00027 eleIsoFromDepsHcalFromTowers04               = eleIsoFromDepsHcalFromTowers.clone()
00028 
00029 # set the correct delta R
00030 eleIsoFromDepsTk03.deposits[0].deltaR                         = 0.3
00031 eleIsoFromDepsTk04.deposits[0].deltaR                         = 0.4
00032 eleIsoFromDepsEcalFromHitsByCrystalFull03.deposits[0].deltaR  = 0.3
00033 eleIsoFromDepsEcalFromHitsByCrystalFull04.deposits[0].deltaR  = 0.4
00034 eleIsoFromDepsEcalFromHitsByCrystalReduced03.deposits[0].deltaR  = 0.3
00035 eleIsoFromDepsEcalFromHitsByCrystalReduced04.deposits[0].deltaR  = 0.4
00036 eleIsoFromDepsHcalFromTowers03.deposits[0].deltaR             = 0.3
00037 eleIsoFromDepsHcalFromTowers04.deposits[0].deltaR             = 0.4
00038 
00039 # change the source on the ECAL hits to make sure to get the reduced or full collection
00040 
00041 eleIsoFromDepsEcalFromHitsByCrystalFull03.deposits[0].src = "eleIsoDepositEcalFromHitsFull"
00042 eleIsoFromDepsEcalFromHitsByCrystalFull04.deposits[0].src = "eleIsoDepositEcalFromHitsFull"
00043 eleIsoFromDepsEcalFromHitsByCrystalReduced03.deposits[0].src = "eleIsoDepositEcalFromHitsReduced"
00044 eleIsoFromDepsEcalFromHitsByCrystalReduced04.deposits[0].src = "eleIsoDepositEcalFromHitsReduced"
00045 
00046 # the sequence
00047 
00048 electronIsoFromDeps = cms.Sequence(
00049   eleIsoDepositTk*eleIsoDepositEcalFromHitsFull*
00050   eleIsoDepositEcalFromHitsReduced*
00051   eleIsoDepositHcalFromTowers*      
00052   eleIsoFromDepsTk03*
00053   eleIsoFromDepsTk04*
00054   eleIsoFromDepsEcalFromHitsByCrystalFull03*
00055   eleIsoFromDepsEcalFromHitsByCrystalFull04*
00056   eleIsoFromDepsEcalFromHitsByCrystalReduced03*
00057   eleIsoFromDepsEcalFromHitsByCrystalReduced04*
00058   eleIsoFromDepsHcalFromTowers03*
00059   eleIsoFromDepsHcalFromTowers04
00060 )
/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_14/src/Validation/RecoEgamma/python/electronIsoFromDeps_cff.py
