PATTauDiscriminationByMVAIsolationRun2_cff.py

Go to the documentation of this file.

import FWCore.ParameterSet.Config as cms
 
from RecoTauTag.RecoTau.TauDiscriminatorTools import noPrediscriminants
from RecoTauTag.RecoTau.PATTauDiscriminantCutMultiplexer_cfi import *
 
# make sure to load the database containing the mva inputs before using the producers below
# e.g. process.load('RecoTauTag.Configuration.loadRecoTauTagMVAsFromPrepDB_cfi') as in
# RecoTauTag.Configuration.HPSPFTaus_cff
 
patDiscriminationByIsolationMVArun2v1raw = cms.EDProducer("PATTauDiscriminationByMVAIsolationRun2",
 
    # tau collection to discriminate
    PATTauProducer = cms.InputTag('replaceMeByTauCollectionToBeUsed'), # in MiniAOD: slimmedTaus
    Prediscriminants = noPrediscriminants,
    loadMVAfromDB = cms.bool(True),
    inputFileName = cms.FileInPath("RecoTauTag/RecoTau/data/emptyMVAinputFile"), # the filename for MVA if it is not loaded from DB
    mvaName = cms.string("replaceMeByNameOfMVATraining"), # e.g. RecoTauTag_tauIdMVADBoldDMwLTv1
    mvaOpt = cms.string("replaceMeByMVAOption"), # e.g. DBoldDMwLT
    
    # change these only if input isolation sums changed for the MVA training you want to use
    srcChargedIsoPtSum = cms.string('chargedIsoPtSum'),
    srcNeutralIsoPtSum = cms.string('neutralIsoPtSum'),
    srcPUcorrPtSum = cms.string('puCorrPtSum'),
    srcPhotonPtSumOutsideSignalCone = cms.string('photonPtSumOutsideSignalCone'),
    srcFootprintCorrection = cms.string('footprintCorrection'),
 
    verbosity = cms.int32(0)
)
 
patDiscriminationByIsolationMVArun2v1 = patTauDiscriminantCutMultiplexer.clone(
    PATTauProducer = 'replaceMeByTauCollectionToBeUsed', # in MiniAOD: slimmedTaus
    Prediscriminants = noPrediscriminants,
    toMultiplex = 'patDiscriminationByIsolationMVArun2v1raw',
    loadMVAfromDB = True,
    mvaOutput_normalization = "replaceMeByNormalizationToBeUsedIfAny", # e.g. RecoTauTag_tauIdMVADBoldDMwLTv1_mvaOutput_normalization
    mapping = cms.VPSet(
        cms.PSet(
            category = cms.uint32(0),
            cut = cms.string("replaceMeByCut"), # e.g. RecoTauTag_tauIdMVADBoldDMwLTv1_WPEff90
            variable = cms.string("pt"),
        )
    ),
    workingPoints = [
        "Eff80",
        "Eff70",
        "Eff60",
        "Eff50",
        "Eff40"
    ]
)
 
mvaIsolation2TaskRun2 = cms.Task(
   patDiscriminationByIsolationMVArun2v1raw
   , patDiscriminationByIsolationMVArun2v1
)
mvaIsolation2SeqRun2 = cms.Sequence(mvaIsolation2TaskRun2)