pfChargedAll = cms.InputTag("elPFIsoValueChargedAll03PFIdPFIso")
)
process.patElectrons.isolationValuesNoPFId = cms.PSet(
pfPUChargedHadrons = cms.InputTag("elPFIsoValuePU03NoPFIdPFIso"),
pfChargedHadrons = cms.InputTag("elPFIsoValueCharged03NoPFIdPFIso"),
pfNeutralHadrons = cms.InputTag("elPFIsoValueNeutral03NoPFIdPFIso"),
pfPhotons = cms.InputTag("elPFIsoValueGamma03NoPFIdPFIso"),
pfChargedAll = cms.InputTag("elPFIsoValueChargedAll03NoPFIdPFIso")
)
process.eleIsoSequence = cms.Sequence(process.pfElectronIsolationSequencePFIso)
process.makePatElectrons = cms.Sequence(process.eleIsoSequence+
process.patElectrons)
###photons
process.phoIsoSequence = setupPFPhotonIso(process, 'photons')
process.photonPFIsolationDepositsSequencePFIso = cms.Sequence(process.phPFIsoDepositChargedPFIso+
process.phPFIsoDepositChargedAllPFIso+
process.phPFIsoDepositGammaPFIso+
process.phPFIsoDepositNeutralPFIso+
process.phPFIsoDepositPUPFIso)
process.pfPhotonIsolationSequencePFIso = cms.Sequence(process.photonPFIsolationDepositsSequencePFIso+
process.phPFIsoValueCharged03PFIdPFIso+
process.phPFIsoValueChargedAll03PFIdPFIso+
process.phPFIsoValueGamma03PFIdPFIso+
process.phPFIsoValueNeutral03PFIdPFIso+
process.phPFIsoValuePU03PFIdPFIso)
process.load('PhysicsTools.PatAlgos.producersLayer1.photonProducer_cfi')
process.patPhotons.isoDeposits = cms.PSet(
user = cms.VInputTag(
#process.selectedPatJetsAK5Calo.cut = cms.string("pt > 15")
# Flavor history stuff
process.load("PhysicsTools.HepMCCandAlgos.flavorHistoryPaths_cfi")
process.flavorHistoryFilter.pathToSelect = cms.int32(-1)
process.cFlavorHistoryProducer.matchedSrc = cms.InputTag("ak5GenJetsNoNu")
process.bFlavorHistoryProducer.matchedSrc = cms.InputTag("ak5GenJetsNoNu")
process.prePathCounter = cms.EDProducer("EventCountProducer")
process.postPathCounter = cms.EDProducer("EventCountProducer")
process.load('TopBrussels.TopTreeProducer.eventCleaning.eventCleaning_cff')
### photon sequence ###
from CommonTools.ParticleFlow.Tools.pfIsolation import setupPFElectronIso, setupPFMuonIso, setupPFPhotonIso
process.phoIsoSequence = setupPFPhotonIso(process, 'selectedPatPhotons')
if runOnMC is False:
process.makePatPhotons.remove(process.photonMatch)
process.photonSequence = cms.Sequence (
process.makePatPhotons+
process.selectedPatPhotons+
process.phoIsoSequence
)
# let it run
process.patseq = cms.Sequence(
process.prePathCounter*
# process.kt6PFJetsForIsolation*
process.goodOfflinePrimaryVertices*
addPfMET(process, 'PF')
process.patMETsPF.metSource = cms.InputTag("pfMet")
process.patPFMETsTypeIcorrected = process.patMETs.clone(
metSource=cms.InputTag('pfType1CorrectedMet'),
addMuonCorrections=cms.bool(False),
genMETSource=cms.InputTag('genMetTrue'),
addGenMET=cms.bool(False))
process.metAnalysisSequence = cms.Sequence(process.patPFMETsTypeIcorrected)
process.cleanPatPhotons.checkOverlaps.electrons.requireNoOverlaps = cms.bool(
False)
from CommonTools.ParticleFlow.Tools.pfIsolation import setupPFElectronIso, setupPFPhotonIso
process.eleIsoSequence = setupPFElectronIso(process, 'gsfElectrons')
process.phoIsoSequence = setupPFPhotonIso(process, 'photons')
process.load("ElectroWeakAnalysis.MultiBosons.VgAnalyzerKit_cfi")
process.VgAnalyzerKit.doGenParticles = cms.bool(False)
process.VgAnalyzerKit.doStoreJets = cms.bool(True)
process.VgAnalyzerKit.doJetHLTMatch = cms.bool(True)
process.VgAnalyzerKit.doStoreSCs = cms.bool(False)
process.VgAnalyzerKit.triggerResults = cms.InputTag("TriggerResults::HLT")
process.VgAnalyzerKit.rhoLabel = cms.InputTag("kt6PFJets25", "rho")
process.VgAnalyzerKit.sigmaLabel = cms.InputTag("kt6PFJets25", "sigma")
process.VgAnalyzerKit.rho2011Label = cms.InputTag("kt6PFJetsForIsolation",
"rho")
process.VgAnalyzerKit.rho2012Label = cms.InputTag("kt6PFJets", "rho")
process.VgAnalyzerKit.TypeIpfMETLabel = cms.InputTag("patPFMETsTypeIcorrected")
process.VgAnalyzerKit.doSkim = cms.bool(False)
#process.selectedPatJetsAK5Calo.cut = cms.string("pt > 15")
# Flavor history stuff
process.load("PhysicsTools.HepMCCandAlgos.flavorHistoryPaths_cfi")
process.flavorHistoryFilter.pathToSelect = cms.int32(-1)
process.cFlavorHistoryProducer.matchedSrc = cms.InputTag("ak5GenJetsNoNu")
process.bFlavorHistoryProducer.matchedSrc = cms.InputTag("ak5GenJetsNoNu")
process.prePathCounter = cms.EDProducer("EventCountProducer")
process.postPathCounter = cms.EDProducer("EventCountProducer")
process.load('TopBrussels.TopTreeProducer.eventCleaning.eventCleaning_cff')
### photon sequence ###
from CommonTools.ParticleFlow.Tools.pfIsolation import setupPFElectronIso, setupPFMuonIso, setupPFPhotonIso
process.phoIsoSequence = setupPFPhotonIso(process, 'selectedPatPhotons')
process.photonSequence = cms.Sequence (
process.makePatPhotons+
process.selectedPatPhotons+
process.phoIsoSequence
)
# let it run
process.patseq = cms.Sequence(
process.prePathCounter*
# process.kt6PFJetsForIsolation*
process.goodOfflinePrimaryVertices*
# process.ak5GenJetsSeq*
process.primaryVertexFilter * #removes events with no good pv (but if cuts to determine good pv change...)
process.eventCleaning*
