/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_14/src/TopQuarkAnalysis/TopEventProducers/python/sequences/ttSemiLepEvtBuilder_cff.py

Go to the documentation of this file.

import FWCore.ParameterSet.Config as cms


################################################################################
# produce ttSemiLepEvent structure with all necessary ingredients
################################################################################

## std sequence to produce the ttSemiLepEventHypotheses
from TopQuarkAnalysis.TopEventProducers.sequences.ttSemiLepEvtHypotheses_cff import *

## configure ttSemiLepEventBuilder
from TopQuarkAnalysis.TopEventProducers.producers.TtSemiLepEvtBuilder_cfi import *

## make ttSemiLepEvent
makeTtSemiLepEvent = cms.Sequence(makeTtSemiLepHypotheses *
                                  ttSemiLepEvent
                                  )


################################################################################
## helper functions
## (examples of usage can be found in the ttSemiLepEvtBuilder_cfg.py)
################################################################################

## add hypotheses to the process
def addTtSemiLepHypotheses(process,
                           names
                           ):

    ## edit list of input hypotheses for the TtSemiLepEventBuilder
    labels =  getattr(process.ttSemiLepEvent, "hypotheses")
    for obj in range(len(names)):
        ## create correct label from HypoClassKey string (stripping the leading "k")
        ## e.g. kKinFit -> ttSemiLepHypKinFit
        label = "ttSemiLepHyp" + names[obj][1:]
        ## add it to the list
        labels.append(label)
    process.ttSemiLepEvent.hypotheses = labels

    ## include hypotheses in the standard sequence
    sequence = getattr(process, "makeTtSemiLepHypotheses")
    for obj in range(len(names)):
        ## create correct label from HypoClassKey string (stripping the leading "k")
        ## e.g. kKinFit -> makeHypothesis_kinFit
        if names[obj][1:4] == "MVA":
            label = "makeHypothesis_" + names[obj][1:4].lower() + names[obj][4:]
        else:
            label = "makeHypothesis_" + names[obj][1:2].lower() + names[obj][2:]
        ## add it to the sequence
        sequence += getattr(process, label)


## remove genMatch hypothesis from the process
def removeTtSemiLepHypGenMatch(process):
    process.makeTtSemiLepHypotheses.remove(process.makeHypothesis_genMatch)
    process.ttSemiLepEvent.hypotheses.remove("ttSemiLepHypGenMatch")


## set a specific attribute for all hypotheses to a given value
## -> this works for "jets", "leps", "mets", "maxNJets"
def setForAllTtSemiLepHypotheses(process, attribute, value):
    modules = ["findTtSemiLepJetCombGeom",
               "findTtSemiLepJetCombMaxSumPtWMass",
               "findTtSemiLepJetCombMVA",
               "findTtSemiLepJetCombWMassDeltaTopMass",
               "findTtSemiLepJetCombWMassMaxSumPt",
               "hitFitTtSemiLepEventHypothesis",
               "kinFitTtSemiLepEventHypothesis",
               "ttSemiLepJetPartonMatch",
               "ttSemiLepHypGenMatch",
               "ttSemiLepHypGeom",
               "ttSemiLepHypHitFit",
               "ttSemiLepHypKinFit",
               "ttSemiLepHypMaxSumPtWMass",
               "ttSemiLepHypMVADisc",               
               "ttSemiLepHypWMassDeltaTopMass",
               "ttSemiLepHypWMassMaxSumPt"
               ]
    for obj in range(len(modules)):
        object = getattr(process, modules[obj])
        if hasattr(object, attribute):
            setattr(object, attribute, value)

## use electrons instead of muons for the hypotheses
def useElectronsForAllTtSemiLepHypotheses(process, elecLabel = "selectedPatElectrons"):
    ## use correct KinFitter module
    import TopQuarkAnalysis.TopKinFitter.TtSemiLepKinFitProducer_Electrons_cfi
    process.kinFitTtSemiLepEventHypothesis = TopQuarkAnalysis.TopKinFitter.TtSemiLepKinFitProducer_Electrons_cfi.kinFitTtSemiLepEvent.clone()
    import TopQuarkAnalysis.TopHitFit.TtSemiLepHitFitProducer_Electrons_cfi
    process.hitFitTtSemiLepEventHypothesis = TopQuarkAnalysis.TopHitFit.TtSemiLepHitFitProducer_Electrons_cfi.hitFitTtSemiLepEvent.clone()
    ## replace lepton InputTags in all modules
    setForAllTtSemiLepHypotheses(process, "leps", elecLabel)