Skip to content

Commit

Permalink
Merge pull request #34571 from jeongeun/recoMET12X
Browse files Browse the repository at this point in the history 
Migrate module configuration in RecoMET to use default cfipython
  • Loading branch information
@cmsbuild
cmsbuild committed Aug 11, 2021
2 parents 55b543e + 0620244 commit ac8d8f2
Showing 1 changed file with 62 additions and 64 deletions.
126 changes: 62 additions & 64 deletions RecoMET/METPUSubtraction/python/mvaPFMET_leptons_cfi.py
View file
Expand Up @@ -4,100 +4,98 @@
isomuons = cms.EDFilter(
"MuonSelector",
src = cms.InputTag('muons'),
cut = cms.string( "(isTrackerMuon) && std::abs(eta) < 2.5 && pt > 9.5"+#17. "+
"&& isPFMuon"+
"&& globalTrack.isNonnull"+
"&& innerTrack.hitPattern.numberOfValidPixelHits > 0"+
"&& innerTrack.normalizedChi2 < 10"+
"&& numberOfMatches > 0"+
"&& innerTrack.hitPattern.numberOfValidTrackerHits>5"+
"&& globalTrack.hitPattern.numberOfValidHits>0"+
"&& (pfIsolationR03.sumChargedHadronPt+pfIsolationR03.sumNeutralHadronEt+pfIsolationR03.sumPhotonEt)/pt < 0.3"+
"&& std::abs(innerTrack().dxy)<2.0"
),
cut = cms.string("(isTrackerMuon) && std::abs(eta) < 2.5 && pt > 9.5"+#17. "+
"&& isPFMuon"+
"&& globalTrack.isNonnull"+
"&& innerTrack.hitPattern.numberOfValidPixelHits > 0"+
"&& innerTrack.normalizedChi2 < 10"+
"&& numberOfMatches > 0"+
"&& innerTrack.hitPattern.numberOfValidTrackerHits>5"+
"&& globalTrack.hitPattern.numberOfValidHits>0"+
"&& (pfIsolationR03.sumChargedHadronPt+pfIsolationR03.sumNeutralHadronEt+pfIsolationR03.sumPhotonEt)/pt < 0.3"+
"&& std::abs(innerTrack().dxy)<2.0"
),
filter = cms.bool(False)
)
)


isoelectrons = cms.EDFilter(
"GsfElectronSelector",
src = cms.InputTag('gsfElectrons'),
cut = cms.string(
"std::abs(eta) < 2.5 && pt > 9.5" +
"&& gsfTrack.trackerExpectedHitsInner.numberOfHits == 0" +
# "&& (pfIsolationVariables.chargedHadronIso+pfIsolationVariables.neutralHadronIso)/et < 0.3" +
"&& (isolationVariables03.tkSumPt)/et < 0.2" +
"&& ((std::abs(eta) < 1.4442 " +
"&& std::abs(deltaEtaSuperClusterTrackAtVtx) < 0.007"+
"&& std::abs(deltaPhiSuperClusterTrackAtVtx) < 0.8" +
"&& sigmaIetaIeta < 0.01" +
"&& hcalOverEcal < 0.15" +
"&& std::abs(1./superCluster.energy - 1./p) < 0.05)"+
"|| (std::abs(eta) > 1.566 "+
"&& std::abs(deltaEtaSuperClusterTrackAtVtx) < 0.009"+
"&& std::abs(deltaPhiSuperClusterTrackAtVtx) < 0.10" +
"&& sigmaIetaIeta < 0.03" +
"&& hcalOverEcal < 0.10" +
"&& std::abs(1./superCluster.energy - 1./p) < 0.05))"
),
filter = cms.bool(False)
)
src = cms.InputTag('gsfElectrons'),
cut = cms.string("std::abs(eta) < 2.5 && pt > 9.5" +
"&& gsfTrack.trackerExpectedHitsInner.numberOfHits == 0" +
# "&& (pfIsolationVariables.chargedHadronIso+pfIsolationVariables.neutralHadronIso)/et < 0.3" +
"&& (isolationVariables03.tkSumPt)/et < 0.2" +
"&& ((std::abs(eta) < 1.4442 " +
"&& std::abs(deltaEtaSuperClusterTrackAtVtx) < 0.007"+
"&& std::abs(deltaPhiSuperClusterTrackAtVtx) < 0.8" +
"&& sigmaIetaIeta < 0.01" +
"&& hcalOverEcal < 0.15" +
"&& std::abs(1./superCluster.energy - 1./p) < 0.05)"+
"|| (std::abs(eta) > 1.566 "+
"&& std::abs(deltaEtaSuperClusterTrackAtVtx) < 0.009"+
"&& std::abs(deltaPhiSuperClusterTrackAtVtx) < 0.10" +
"&& sigmaIetaIeta < 0.03" +
"&& hcalOverEcal < 0.10" +
"&& std::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
)
doRhoFastjet = True,
voronoiRfact = 0.9
)

from RecoTauTag.RecoTau.PFRecoTauDiscriminationByHPSSelection_cfi import hpsSelectionDiscriminator
hpsPFTauDiscriminationByDecayModeFinding = hpsSelectionDiscriminator.clone(
PFTauProducer = 'hpsPFTauProducer'
)
PFTauProducer = 'hpsPFTauProducer'
)

from RecoTauTag.RecoTau.TauDiscriminatorTools import requireLeadTrack
# Define decay mode prediscriminant
requireDecayMode = cms.PSet(
BooleanOperator = cms.string("and"),
decayMode = cms.PSet(
Producer = cms.InputTag('hpsPFTauDiscriminationByDecayModeFinding'),
cut = cms.double(0.5)
)
)
BooleanOperator = cms.string("and"),
decayMode = cms.PSet(
Producer = cms.InputTag('hpsPFTauDiscriminationByDecayModeFinding'),
cut = cms.double(0.5)
)
)

from RecoTauTag.Configuration.HPSPFTaus_cff import hpsPFTauDiscriminationByLooseCombinedIsolationDBSumPtCorr3Hits
import RecoTauTag.RecoTau.pfRecoTauDiscriminationAgainstMuon2_cfi as _mod

hpsPFTauDiscriminationAgainstMuon2 = cms.EDProducer("PFRecoTauDiscriminationAgainstMuon2",
PFTauProducer = cms.InputTag('hpsPFTauProducer'),
Prediscriminants = requireDecayMode.clone(),
discriminatorOption = cms.string('loose'), # available options are: 'loose', 'medium', 'tight'
HoPMin = cms.double(0.2)
)

hpsPFTauDiscriminationAgainstMuon2 = _mod.pfRecoTauDiscriminationAgainstMuon2.clone(
PFTauProducer = 'hpsPFTauProducer',
Prediscriminants = requireDecayMode.clone(),
discriminatorOption = 'loose', # available options are: 'loose', 'medium', 'tight'
)

hpsPFTauDiscriminationByMVAIsolation = cms.EDProducer(
"PFRecoTauDiscriminationByMVAIsolation",
PFTauProducer = cms.InputTag('hpsPFTauProducer'),
rhoProducer = cms.InputTag('kt6PFJetsForRhoComputationVoronoiMet','rho'),
Prediscriminants = requireDecayMode.clone(),
gbrfFilePath = cms.FileInPath('RecoTauTag/RecoTau/data/gbrfTauIso_v2.root'),
returnMVA = cms.bool(False),
mvaMin = cms.double(0.8),
)
PFTauProducer = cms.InputTag('hpsPFTauProducer'),
rhoProducer = cms.InputTag('kt6PFJetsForRhoComputationVoronoiMet','rho'),
Prediscriminants = requireDecayMode.clone(),
gbrfFilePath = cms.FileInPath('RecoTauTag/RecoTau/data/gbrfTauIso_v2.root'),
returnMVA = cms.bool(False),
mvaMin = cms.double(0.8),
)

isotaus = cms.EDFilter(
"PFTauSelector",
src = cms.InputTag('hpsPFTauProducer'),
BooleanOperator = cms.string("and"),
discriminators = cms.VPSet(
cms.PSet( discriminator=cms.InputTag("hpsPFTauDiscriminationByDecayModeFinding"), selectionCut=cms.double(0.5)),
#cms.PSet( discriminator=cms.InputTag("hpsPFTauDiscriminationByMVAIsolation"), selectionCut=cms.double(0.5)),
cms.PSet( discriminator=cms.InputTag("hpsPFTauDiscriminationByLooseCombinedIsolationDBSumPtCorr3Hits"), selectionCut=cms.double(0.5)),
cms.PSet( discriminator=cms.InputTag("hpsPFTauDiscriminationByLooseElectronRejection"), selectionCut=cms.double(0.5)),
cms.PSet( discriminator=cms.InputTag("hpsPFTauDiscriminationAgainstMuon2"), selectionCut=cms.double(0.5))
cms.PSet( discriminator=cms.InputTag("hpsPFTauDiscriminationByDecayModeFinding"), selectionCut=cms.double(0.5)),
#cms.PSet( discriminator=cms.InputTag("hpsPFTauDiscriminationByMVAIsolation"), selectionCut=cms.double(0.5)),
cms.PSet( discriminator=cms.InputTag("hpsPFTauDiscriminationByLooseCombinedIsolationDBSumPtCorr3Hits"), selectionCut=cms.double(0.5)),
cms.PSet( discriminator=cms.InputTag("hpsPFTauDiscriminationByLooseElectronRejection"), selectionCut=cms.double(0.5)),
cms.PSet( discriminator=cms.InputTag("hpsPFTauDiscriminationAgainstMuon2"), selectionCut=cms.double(0.5))
),
cut = cms.string("std::abs(eta) < 2.3 && pt > 19.0 "),
filter = cms.bool(False)
)
)

isomuonTask = cms.Task(isomuons)
isomuonseq = cms.Sequence(isomuonsTask)
Expand Down

0 comments on commit ac8d8f2

Please sign in to comment.