RecoTauCleanerPlugins.py

Go to the documentation of this file.

import FWCore.ParameterSet.Config as cms

'''

Plugins for ranking PFTau candidates

'''

tolerance_default = cms.double(0)

matchingConeCut = cms.PSet(
    name = cms.string("MatchingCone"),
    plugin = cms.string("RecoTauStringCleanerPlugin"),
    # Prefer taus that are within DR<0.1 of the jet axis
    selection = cms.string("deltaR(eta, phi, jetRef().eta, jetRef().phi) < 0.1"),
    selectionPassFunction = cms.string("0"),
    selectionFailValue = cms.double(1e3),
    tolerance = tolerance_default,
)

# Prefer taus with charge == 1 (no three prongs with charge = 3)
unitCharge = cms.PSet(
    name = cms.string("UnitCharge"),
    plugin = cms.string("RecoTauStringCleanerPlugin"),
    # Only effects three prongs
    selection = cms.string("signalChargedHadrCands().size() = 3"),
    # As 1 is lower than 3, this will always prefer those with unit charge
    selectionPassFunction = cms.string("abs(charge())-1"),
    # If it is a one prong, consider it just as good as a
    # three prong with unit charge
    selectionFailValue = cms.double(0),
    tolerance = tolerance_default,
)

# similar to unitCharge but handles also cases where tau is made up of
# a combination of tracks and pf charged hadrons
charge = cms.PSet(
    name = cms.string("Charge"),
    plugin = cms.string("RecoTauChargeCleanerPlugin"),
    # cleaner is applied to decay modes with the number of prongs given here
    nprongs = cms.vuint32(1,3),
    # taus with charge != 1 are rejected
    passForCharge = cms.int32(1),
    selectionFailValue = cms.double(0),
    tolerance = tolerance_default,
)

# Prefer taus with pt greater 15
ptGt15 = cms.PSet(
    name = cms.string("PtGt15"),
    plugin = cms.string("RecoTauStringCleanerPlugin"),
    selection = cms.string("pt > 15."),
    selectionPassFunction = cms.string("0"),
    selectionFailValue = cms.double(1e3),
    tolerance = tolerance_default,
)

leadPionFinding = cms.PSet(
    name = cms.string("LeadPion"),
    plugin = cms.string("RecoTauDiscriminantCleanerPlugin"),
    src = cms.InputTag("DISCRIMINATOR_SRC"),
    tolerance = tolerance_default,
)

pt = cms.PSet(
    name = cms.string("Pt"),
    plugin = cms.string("RecoTauStringCleanerPlugin"),
    # Require that cones were built by ensuring the a leadCand exits
    selection = cms.string("leadCand().isNonnull()"),
    selectionPassFunction = cms.string("-pt()"), # CV: negative sign means that we prefer candidates of high pT
    selectionFailValue = cms.double(1e3),
    tolerance = cms.double(1.e-2) # CV: consider candidates with almost equal pT to be of the same rank (to avoid sensitivity to rounding errors)
)

chargedHadronMultiplicity = cms.PSet(
    name = cms.string("ChargedHadronMultiplicity"),
    plugin = cms.string("RecoTauChargedHadronMultiplicityCleanerPlugin"),
    tolerance = tolerance_default,
)

stripMultiplicity = cms.PSet(
    name = cms.string("StripMultiplicity"),
    plugin = cms.string("RecoTauStringCleanerPlugin"),
    # Require that cones were built by ensuring the a leadCand exits
    selection = cms.string("leadCand().isNonnull()"),
    selectionPassFunction = cms.string("-signalPiZeroCandidates().size()"),
    selectionFailValue = cms.double(1e3),
    tolerance = tolerance_default,
)

combinedIsolation = cms.PSet(
    name = cms.string("CombinedIsolation"),
    plugin = cms.string("RecoTauStringCleanerPlugin"),
    # Require that cones were built by ensuring the a leadCand exits
    selection = cms.string("leadCand().isNonnull()"),
    selectionPassFunction = cms.string("isolationPFChargedHadrCandsPtSum() + isolationPFGammaCandsEtSum()"),
    selectionFailValue = cms.double(1e3),
    tolerance = tolerance_default,
)

chargeIsolation = cms.PSet(
    name = cms.string("ChargeIsolation"),
    plugin = cms.string("RecoTauStringCleanerPlugin"),
    # Require that cones were built by ensuring the a leadCand exits
    selection = cms.string("leadCand().isNonnull()"),
    # Prefer lower isolation activity
    selectionPassFunction = cms.string("isolationPFChargedHadrCandsPtSum()"),
    selectionFailValue = cms.double(1e3),
    tolerance = tolerance_default,
)

ecalIsolation = cms.PSet(
    name = cms.string("GammaIsolation"),
    plugin = cms.string("RecoTauStringCleanerPlugin"),
    # Require that cones were built by ensuring the a leadCand exits
    selection = cms.string("leadCand().isNonnull()"),
    # Prefer lower isolation activity
    selectionPassFunction = cms.string("isolationPFGammaCandsEtSum()"),
    selectionFailValue = cms.double(1e3),
    tolerance = tolerance_default,
)

# CV: Reject 2-prong candidates in which one of the tracks has low pT,
#     in order to reduce rate of 1-prong taus migrating to 2-prong decay modecanddates
killSoftTwoProngTaus = cms.PSet(
    name = cms.string("killSoftTwoProngTaus"),
    plugin = cms.string("RecoTauSoftTwoProngTausCleanerPlugin"),
    minTrackPt = cms.double(5.),
    tolerance = tolerance_default,
)