EcalAlCaRecoProducers/python/eleIsoSequence_cff.py

0001 #from RecoParticleFlow.PFProducer.electronPFIsolationValues_cff import *
0002 #from CommonTools.ParticleFlow.PFBRECO_cff import *
0003 from CommonTools.ParticleFlow.Isolation.tools_cfi import *
0004 import FWCore.ParameterSet.Config as cms
0005 #Now prepare the iso deposits
0006 elPFIsoDepositChargedGsf=isoDepositReplace('gedGsfElectrons','pfAllChargedHadrons')
0007 elPFIsoDepositChargedAllGsf=isoDepositReplace('gedGsfElectrons','pfAllChargedParticles')
0008 elPFIsoDepositNeutralGsf=isoDepositReplace('gedGsfElectrons','pfAllNeutralHadrons')
0009 elPFIsoDepositPUGsf=isoDepositReplace('gedGsfElectrons','pfPileUpAllChargedParticles')
0010 #elPFIsoDepositGamma=isoDepositReplace('pfSelectedElectrons','pfAllPhotons')
0011 elPFIsoDepositGammaGsf= cms.EDProducer("CandIsoDepositProducer",
0012                                        src = cms.InputTag("gedGsfElectrons"),
0013                                        MultipleDepositsFlag = cms.bool(False),
0014                                        trackType = cms.string('candidate'),
0015                                        ExtractorPSet = cms.PSet(
0016     Diff_z = cms.double(99999.99),
0017     ComponentName = cms.string('PFCandWithSuperClusterExtractor'),
0018     DR_Max = cms.double(1.0),
0019     Diff_r = cms.double(99999.99),
0020     inputCandView = cms.InputTag("pfAllPhotons"),
0021     DR_Veto = cms.double(0),
0022     SCMatch_Veto = cms.bool(False),
0023     MissHitSCMatch_Veto = cms.bool(True),
0024     DepositLabel = cms.untracked.string('')
0025     )
0026                                        )
0027
0028 # elPFIsoDepositChargedGsf= elPFIsoDepositCharged.clone()
0029 # elPFIsoDepositChargedGsf.src = 'gedGsfElectrons'
0030 # elPFIsoDepositChargedAllGsf = elPFIsoDepositChargedAll.clone()
0031 # elPFIsoDepositChargedAllGsf.src = 'gedGsfElectrons'
0032 # elPFIsoDepositNeutralGsf = elPFIsoDepositNeutral.clone()
0033 # elPFIsoDepositNeutralGsf.src = 'gedGsfElectrons'
0034 # elPFIsoDepositGammaGsf = elPFIsoDepositGamma.clone()
0035 # elPFIsoDepositGammaGsf.src = 'gedGsfElectrons'
0036 # elPFIsoDepositPUGsf = elPFIsoDepositPU.clone()
0037 # elPFIsoDepositPUGsf.src = 'gedGsfElectrons'
0038
0039 elPFIsoValueCharged03PFIdGsf = cms.EDProducer("PFCandIsolatorFromDeposits",
0040                                            deposits = cms.VPSet(
0041     cms.PSet(
0042     src = cms.InputTag("elPFIsoDepositChargedGsf"),
0043     deltaR = cms.double(0.3),
0044     weight = cms.string('1'),
0045     vetos = cms.vstring('EcalEndcaps:ConeVeto(0.015)'),
0046     skipDefaultVeto = cms.bool(True),
0047     mode = cms.string('sum'),
0048     PivotCoordinatesForEBEE = cms.bool(True)
0049     )
0050     )
0051                                            )
0052 elPFIsoValueChargedAll03PFIdGsf = cms.EDProducer("PFCandIsolatorFromDeposits",
0053                                                  deposits = cms.VPSet(
0054           cms.PSet(
0055           src = cms.InputTag("elPFIsoDepositChargedAllGsf"),
0056           deltaR = cms.double(0.3),
0057           weight = cms.string('1'),
0058           vetos = cms.vstring('EcalEndcaps:ConeVeto(0.015)'),
0059           skipDefaultVeto = cms.bool(True),
0060           mode = cms.string('sum'),
0061           PivotCoordinatesForEBEE = cms.bool(True)
0062    )
0063  )
0064                                               )
0065
0066 elPFIsoValueGamma03PFIdGsf = cms.EDProducer("PFCandIsolatorFromDeposits",
0067                                          deposits = cms.VPSet(
0068            cms.PSet(
0069            src = cms.InputTag("elPFIsoDepositGammaGsf"),
0070            deltaR = cms.double(0.3),
0071            weight = cms.string('1'),
0072            vetos = cms.vstring('EcalEndcaps:ConeVeto(0.08)'),
0073            skipDefaultVeto = cms.bool(True),
0074            mode = cms.string('sum'),
0075            PivotCoordinatesForEBEE = cms.bool(True)
0076      )
0077   )
0078 )
0079
0080 elPFIsoValueNeutral03PFIdGsf = cms.EDProducer("PFCandIsolatorFromDeposits",
0081                                            deposits = cms.VPSet(
0082              cms.PSet(
0083              src = cms.InputTag("elPFIsoDepositNeutralGsf"),
0084              deltaR = cms.double(0.3),
0085              weight = cms.string('1'),
0086              vetos = cms.vstring(),
0087              skipDefaultVeto = cms.bool(True),
0088              mode = cms.string('sum'),
0089              PivotCoordinatesForEBEE = cms.bool(True)
0090              )
0091          )
0092                                            )
0093
0094 elPFIsoValuePU03PFIdGsf = cms.EDProducer("PFCandIsolatorFromDeposits",
0095                                       deposits = cms.VPSet(
0096            cms.PSet(
0097            src = cms.InputTag("elPFIsoDepositPUGsf"),
0098            deltaR = cms.double(0.3),
0099            weight = cms.string('1'),
0100            vetos = cms.vstring('EcalEndcaps:ConeVeto(0.015)'),
0101            skipDefaultVeto = cms.bool(True),
0102            mode = cms.string('sum'),
0103            PivotCoordinatesForEBEE = cms.bool(True)
0104            )
0105   )
0106                                          )
0107
0108
0109 # elPFIsoValueCharged03PFIdGsf = elPFIsoValueCharged03PFId.clone()
0110 # #elPFIsoValueCharged03PFIdGsf.deposits[0].src = cms.InputTag(elPFIsoDepositChargedGsf)
0111
0112 # elPFIsoValueChargedAll03PFIdGsf = elPFIsoValueChargedAll03PFId.clone()
0113 # #elPFIsoValueChargedAll03PFIdGsf.deposits[0].src = cms.InputTag(elPFIsoDepositChargedAllGsf)
0114
0115 # elPFIsoValueGamma03PFIdGsf = elPFIsoValueGamma03PFId.clone()
0116 # #elPFIsoValueGamma03PFIdGsf.deposits[0].src = cms.InputTag(elPFIsoDepositGammaGsf)
0117
0118
0119 # elPFIsoValueNeutral03PFIdGsf = elPFIsoValueNeutral03PFId.clone()
0120 # #elPFIsoValueNeutral03PFIdGsf.deposits[0].src = cms.InputTag(elPFIsoDepositNeutralGsf)
0121
0122 # elPFIsoValuePU03PFIdGsf = elPFIsoValuePU03PFId.clone()
0123 # #elPFIsoValuePU03PFIdGsf.deposits[0].src = cms.InputTag(elPFIsoDepositPUGsf)
0124
0125 eleIsoSequence = cms.Sequence((elPFIsoDepositChargedGsf + elPFIsoDepositChargedAllGsf + elPFIsoDepositNeutralGsf + elPFIsoDepositGammaGsf + elPFIsoDepositPUGsf))
0126 eleIsoSequence *= cms.Sequence(elPFIsoValueCharged03PFIdGsf+elPFIsoValueChargedAll03PFIdGsf+elPFIsoValueGamma03PFIdGsf+elPFIsoValueNeutral03PFIdGsf+elPFIsoValuePU03PFIdGsf )
0127 pfisoALCARECO = cms.Sequence(eleIsoSequence) #pfParticleSelectionSequence + eleIsoSequence)