hpsPFTauDiscriminationByLooseIsolationMVA2 = hpsPFTauDiscriminationByDecayModeFinding.clone(
Prediscriminants=cms.PSet(BooleanOperator=cms.string("and"),
mva=cms.PSet(Producer=cms.InputTag(
'hpsPFTauDiscriminationByIsolationMVA2raw'),
cut=cms.double(0.85))))
hpsPFTauDiscriminationByMediumIsolationMVA2 = copy.deepcopy(
hpsPFTauDiscriminationByLooseIsolationMVA2)
hpsPFTauDiscriminationByMediumIsolationMVA2.Prediscriminants.mva.cut = cms.double(
0.90)
hpsPFTauDiscriminationByTightIsolationMVA2 = copy.deepcopy(
hpsPFTauDiscriminationByLooseIsolationMVA2)
hpsPFTauDiscriminationByTightIsolationMVA2.Prediscriminants.mva.cut = cms.double(
0.94)
from RecoJets.Configuration.RecoPFJets_cff import kt6PFJets as _dummy
kt6PFJetsForRhoComputationVoronoi = _dummy.clone(doRhoFastjet=True,
voronoiRfact=0.9)
hpsPFTauDiscriminationByMVAIsolationSeq = cms.Sequence(
kt6PFJetsForRhoComputationVoronoi *
hpsPFTauDiscriminationByIsolationMVAraw *
hpsPFTauDiscriminationByLooseIsolationMVA *
hpsPFTauDiscriminationByMediumIsolationMVA *
hpsPFTauDiscriminationByTightIsolationMVA *
hpsPFTauDiscriminationByIsolationMVA2raw *
hpsPFTauDiscriminationByLooseIsolationMVA2 *
hpsPFTauDiscriminationByMediumIsolationMVA2 *
hpsPFTauDiscriminationByTightIsolationMVA2)
#copying discriminator against electrons and muons
hpsPFTauDiscriminationByLooseElectronRejection = pfRecoTauDiscriminationAgainstElectron.clone(
PFTauProducer=cms.InputTag('hpsPFTauProducer'),
"&& sigmaIetaIeta < 0.01" +
"&& hcalOverEcal < 0.15" +
"&& abs(1./superCluster.energy - 1./p) < 0.05)"+
"|| (abs(eta) > 1.566 "+
"&& abs(deltaEtaSuperClusterTrackAtVtx) < 0.009"+
"&& abs(deltaPhiSuperClusterTrackAtVtx) < 0.10" +
"&& sigmaIetaIeta < 0.03" +
"&& hcalOverEcal < 0.10" +
"&& abs(1./superCluster.energy - 1./p) < 0.05))"
),
filter = cms.bool(False)
)
from RecoJets.Configuration.RecoPFJets_cff import kt6PFJets as dummy
kt6PFJetsForRhoComputationVoronoiMet = dummy.clone(
doRhoFastjet = True,
voronoiRfact = 0.9
)
from RecoTauTag.RecoTau.PFRecoTauDiscriminationByHPSSelection_cfi import hpsSelectionDiscriminator
hpsPFTauDiscriminationByDecayModeFindingMVAMET = hpsSelectionDiscriminator.clone(
PFTauProducer = cms.InputTag('hpsPFTauProducer')
)
from RecoTauTag.RecoTau.TauDiscriminatorTools import requireLeadTrack
# Define decay mode prediscriminant
requireDecayMode = cms.PSet(
BooleanOperator = cms.string("and"),
decayMode = cms.PSet(
Producer = cms.InputTag('hpsPFTauDiscriminationByDecayModeFindingMVAMET'),
cut = cms.double(0.5)
)
# minNumber = cms.uint32(0),
# maxNumber = cms.uint32(0),
# src = cms.InputTag("trackerDrivenOnlyElectrons")
# )
#trackerDrivenRemoverSeq = cms.Sequence( trackerDrivenOnlyElectrons * trackerDrivenRemover )
#trackerDrivenRemoverSeq = cms.Sequence( trackerDrivenOnlyElectrons)
from Calibration.EcalAlCaRecoProducers.alCaIsolatedElectrons_cfi import *
from Calibration.EcalAlCaRecoProducers.AlCaElectronTracksReducer_cfi import *
from Calibration.EcalAlCaRecoProducers.eleIsoSequence_cff import *
from Calibration.EcalAlCaRecoProducers.WZElectronSkims_cff import *
from RecoJets.Configuration.RecoPFJets_cff import kt6PFJets
kt6PFJetsForRhoCorrection = kt6PFJets.clone(doRhoFastjet = True)
kt6PFJetsForRhoCorrection.Rho_EtaMax = cms.double(2.5)
#list of SCs to be used for the recHit reduction
#GsfMatchedPhotonCands = cms.EDProducer("ElectronMatchedCandidateProducer",
# src = cms.InputTag("goodPhotons"),
# ReferenceElectronCollection = cms.untracked.InputTag("goodElectrons"),
# deltaR = cms.untracked.double(0.3)
#)
alcarecoEcalRecHitReducerSeq = cms.Sequence(alCaIsolatedElectrons)
alcarecoElectronTracksReducerSeq = cms.Sequence(alcaElectronTracksReducer)
# sequence that reduces the RECO format (only ECAL part) into ALCARECO
ALCARECOEcalCalElectronECALSeq = cms.Sequence( alCaIsolatedElectrons)
applyPostfix( process, 'pfElectronsFromVertex' , postfix ).vertices = cms.InputTag( pfVertices )
applyPostfix( process, 'pfElectronsFromVertex' , postfix ).d0Cut = pfD0Cut
applyPostfix( process, 'pfElectronsFromVertex' , postfix ).dzCut = pfDzCut
applyPostfix( process, 'pfSelectedElectrons' , postfix ).cut = pfElectronSelectionCut
#applyPostfix( process, 'isoValElectronWithCharged', postfix ).deposits[0].deltaR = pfElectronIsoConeR
#applyPostfix( process, 'isoValElectronWithNeutral', postfix ).deposits[0].deltaR = pfElectronIsoConeR
#applyPostfix( process, 'isoValElectronWithPhotons', postfix ).deposits[0].deltaR = pfElectronIsoConeR
#applyPostfix( process, 'pfIsolatedElectrons' , postfix ).combinedIsolationCut = pfElectronCombIsoCut
applyPostfix( process, 'pfIsolatedElectrons' , postfix ).isolationCut = pfElectronCombIsoCut
applyPostfix( process, 'patElectrons', postfix ).pvSrc = cms.InputTag( pfVertices )
applyPostfix( process, 'patMuons', postfix ).pvSrc = cms.InputTag( pfVertices )
from RecoJets.Configuration.RecoPFJets_cff import kt6PFJets
kt6PFJetsChs = kt6PFJets.clone(
rParam = cms.double( 0.6 )
, src = cms.InputTag( 'pfNoElectron' )
, doAreaFastjet = cms.bool( True )
, doRhoFastjet = cms.bool( True )
, voronoiRfact = cms.double( -0.9 )
)
applyPostfix( process, 'patJetCorrFactors', postfix ).primaryVertices = cms.InputTag( pfVertices )
kt6PFJetsPFChs = kt6PFJetsChs.clone( src = cms.InputTag( 'pfNoElectron' + postfix ) )
setattr( process, 'kt6PFJetsChs' + postfix, kt6PFJetsPFChs )
getattr( process, 'patPF2PATSequence' + postfix).replace( getattr( process, 'patJetCorrFactors' + postfix )
, getattr( process, 'kt6PFJetsChs' + postfix ) * getattr( process, 'patJetCorrFactors' + postfix )
)
applyPostfix( process, 'patJetCorrFactors', postfix ).rho = cms.InputTag( 'kt6PFJetsChs' + postfix, 'rho' )
# In order to have a coherent semileptonic channel also, add
# some "loose" leptons to do QCD estimates.
process.pfIsolatedMuonsLoosePFlow = process.pfIsolatedMuonsPFlow.clone(
combinedIsolationCut = cms.double(999.0)
process.pfPileUpPFlow.Vertices = cms.InputTag('goodOfflinePrimaryVertices')
process.pfPileUpPFlow.verbose = True
process.pfJetsPFlow.doAreaFastjet = True
process.pfJetsPFlow.doRhoFastjet = False
from PhysicsTools.SelectorUtils.pvSelector_cfi import pvSelector
process.goodOfflinePrimaryVertices = cms.EDFilter(
"PrimaryVertexObjectFilter",
filterParams=pvSelector.clone(minNdof=cms.double(4.0),
maxZ=cms.double(24.0)),
src=cms.InputTag('offlinePrimaryVertices'))
# Compute the mean pt per unit area (rho) from the PFCHS inputs
from RecoJets.Configuration.RecoPFJets_cff import kt6PFJets, ak5PFJets
process.kt6PFJetsPFlow = kt6PFJets.clone(src=cms.InputTag('pfNoElectron' +
postfix),
doAreaFastjet=cms.bool(True),
doRhoFastjet=cms.bool(True))
process.patJetCorrFactorsPFlow.rho = cms.InputTag("kt6PFJetsPFlow", "rho")
process.ak5PFJetsPileUp = ak5PFJets.clone(src=cms.InputTag('pfPileUp' +
postfix))
"""
process.pileUpN1PrimaryVertices = cms.EDFilter("RecoTauPileUpVertexSelector",
filter = cms.bool(False),
src = cms.InputTag('goodOfflinePrimaryVertices'),
minTrackSumPt = cms.double(0.)
#minTrackSumPt = cms.double(5.)
)
process.pileUpN2PrimaryVertices = process.pileUpN1PrimaryVertices.clone(
src = cms.InputTag('pileUpN1PrimaryVertices')
def rhoFor2012Aeff(process,postfix): from RecoJets.Configuration.RecoPFJets_cff import kt6PFJets process.kt6PFJetsForIsolation = kt6PFJets.clone(doRhoFastjet = True,Rho_EtaMax = 2.5) process.p += process.kt6PFJetsForIsolation
process.pfJetsPFlow.doAreaFastjet = True
process.pfJetsPFlow.doRhoFastjet = False
from PhysicsTools.SelectorUtils.pvSelector_cfi import pvSelector
process.goodOfflinePrimaryVertices = cms.EDFilter(
"PrimaryVertexObjectFilter",
filterParams=pvSelector.clone(minNdof=cms.double(4.0), maxZ=cms.double(24.0)),
src=cms.InputTag("offlinePrimaryVertices"),
)
# Compute the mean pt per unit area (rho) from the PFCHS inputs
from RecoJets.Configuration.RecoPFJets_cff import kt6PFJets, ak5PFJets
process.kt6PFJetsPFlow = kt6PFJets.clone(
src=cms.InputTag("pfNoElectron" + postfix), doAreaFastjet=cms.bool(True), doRhoFastjet=cms.bool(True)
)
process.patJetCorrFactorsPFlow.rho = cms.InputTag("kt6PFJetsPFlow", "rho")
process.ak5PFJetsPileUp = ak5PFJets.clone(src=cms.InputTag("pfPileUp" + postfix))
"""
process.pileUpN1PrimaryVertices = cms.EDFilter("RecoTauPileUpVertexSelector",
filter = cms.bool(False),
src = cms.InputTag('goodOfflinePrimaryVertices'),
minTrackSumPt = cms.double(0.)
#minTrackSumPt = cms.double(5.)
)
process.pileUpN2PrimaryVertices = process.pileUpN1PrimaryVertices.clone(
src = cms.InputTag('pileUpN1PrimaryVertices')
import FWCore.ParameterSet.Config as cms from RecoJets.Configuration.RecoPFJets_cff import kt6PFJets kt6PFJets = kt6PFJets.clone( rParam = 0.6, doRhoFastjet = True )
#trackerDrivenRemover = cms.EDFilter("PATCandViewCountFilter",
# minNumber = cms.uint32(0),
# maxNumber = cms.uint32(0),
# src = cms.InputTag("trackerDrivenOnlyElectrons")
# )
#trackerDrivenRemoverSeq = cms.Sequence( trackerDrivenOnlyElectrons * trackerDrivenRemover )
#trackerDrivenRemoverSeq = cms.Sequence( trackerDrivenOnlyElectrons)
from Calibration.EcalAlCaRecoProducers.alCaIsolatedElectrons_cfi import *
from Calibration.EcalAlCaRecoProducers.AlCaElectronTracksReducer_cfi import *
from Calibration.EcalAlCaRecoProducers.eleIsoSequence_cff import *
from Calibration.EcalAlCaRecoProducers.WZElectronSkims_cff import *
from RecoJets.Configuration.RecoPFJets_cff import kt6PFJets
kt6PFJetsForRhoCorrection = kt6PFJets.clone(doRhoFastjet=True)
kt6PFJetsForRhoCorrection.Rho_EtaMax = cms.double(2.5)
#list of SCs to be used for the recHit reduction
#GsfMatchedPhotonCands = cms.EDProducer("ElectronMatchedCandidateProducer",
# src = cms.InputTag("goodPhotons"),
# ReferenceElectronCollection = cms.untracked.InputTag("goodElectrons"),
# deltaR = cms.untracked.double(0.3)
#)
alcarecoEcalRecHitReducerSeq = cms.Sequence(alCaIsolatedElectrons)
alcarecoElectronTracksReducerSeq = cms.Sequence(alcaElectronTracksReducer)
# sequence that reduces the RECO format (only ECAL part) into ALCARECO
ALCARECOEcalCalElectronECALSeq = cms.Sequence(alCaIsolatedElectrons)
import FWCore.ParameterSet.Config as cms
#needed for MLP tagger
from PhysicsTools.SelectorUtils.pvSelector_cfi import pvSelector
goodOfflinePrimaryVerticesQG = cms.EDFilter(
"PrimaryVertexObjectFilter",
filterParams=pvSelector.clone(minNdof=cms.double(4.0),
maxZ=cms.double(24.0)),
src=cms.InputTag('offlinePrimaryVertices'))
from RecoJets.Configuration.RecoPFJets_cff import kt6PFJets
kt6PFJetsQG = kt6PFJets.clone()
#kt6PFJetsQG = kt6PFJetsPFlow.clone()
#kt6PFJetsQG = kt6PFJets.clone( src = cms.InputTag( 'particleFlow' )
# , doRhoFastjet = True
# )
#needed for Likelihood tagger
from RecoJets.JetProducers.kt4PFJets_cfi import *
#from RecoJets.Configuration.RecoPFJets_cff import kt6PFJets
kt6PFJetsIsoQG = kt4PFJets.clone(rParam=0.6, doRhoFastjet=True)
kt6PFJetsIsoQG.Rho_EtaMax = cms.double(2.5)
QGTagger = cms.EDProducer(
'QGTagger',
# srcRho = cms.InputTag('kt6PFJetsQG','rho'),
srcRho=cms.InputTag('kt6PFJetsIsoQG', 'rho'),
srcRhoIso=cms.InputTag('kt6PFJetsIsoQG', 'rho'),
)
QuarkGluonTagger = cms.Sequence(kt6PFJetsQG + kt6PFJetsIsoQG + QGTagger)
import FWCore.ParameterSet.Config as cms
#needed for MLP tagger
from PhysicsTools.SelectorUtils.pvSelector_cfi import pvSelector
goodOfflinePrimaryVerticesQG = cms.EDFilter("PrimaryVertexObjectFilter",
filterParams = pvSelector.clone( minNdof = cms.double(4.0), maxZ = cms.double(24.0) ),
src = cms.InputTag('offlinePrimaryVertices')
)
from RecoJets.Configuration.RecoPFJets_cff import kt6PFJets
kt6PFJetsQG = kt6PFJets.clone()
#needed for Likelihood tagger
from RecoJets.JetProducers.kt4PFJets_cfi import *
kt6PFJetsIsoQG = kt4PFJets.clone( rParam = 0.6, doRhoFastjet = True )
kt6PFJetsIsoQG.Rho_EtaMax = cms.double(2.5)
QGTagger = cms.EDProducer('QGTagger',
srcRho = cms.InputTag('kt6PFJetsQG','rho'),
srcRhoIso = cms.InputTag('kt6PFJetsIsoQG','rho'),
)
QuarkGluonTagger = cms.Sequence(goodOfflinePrimaryVerticesQG + kt6PFJetsQG + kt6PFJetsIsoQG + QGTagger)