ttFullHadEvtBuilder_cff.py

Go to the documentation of this file.

import FWCore.ParameterSet.Config as cms

#
# produce ttFullHadEvent structure with all necessary ingredients,
# needs ttGenEvent as input
#

## std sequence to produce the ttFullHadEventHypotheses
from TopQuarkAnalysis.TopEventProducers.sequences.ttFullHadEvtHypotheses_cff import *

## configure ttFullHadEventBuilder
from TopQuarkAnalysis.TopEventProducers.producers.TtFullHadEvtBuilder_cfi import *

### make ttFullHadEvent
makeTtFullHadEventTask = cms.Task(
    makeTtFullHadHypothesesTask,
    ttFullHadEvent
)
#makeTtFullHadEvent = cms.Sequence(makeTtFullHadEventTask)

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

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

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

    ### include hypotheses in the standard sequence
    #sequence = getattr(process, "makeTtFullHadHypotheses")
    #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 removeTtFullHadHypGenMatch(process):
    #process.makeTtFullHadHypotheses.remove(process.makeHypothesis_genMatch)
    process.ttFullHadEvent.hypotheses.remove("ttFullHadHypGenMatch")


## set a specific attribute for all hypotheses to a given value
## -> this works for "jets", "maxNJets", "jetCorrectionLevel"
def setForAllTtFullHadHypotheses(process, attribute, value):
    modules = ["ttFullHadJetPartonMatch",
               "ttFullHadHypGenMatch",
               "kinFitTtFullHadEventHypothesis",
               "ttFullHadHypKinFit"]
    for obj in range(len(modules)):
        object = getattr(process, modules[obj])
        if hasattr(object, attribute):
            setattr(object, attribute, value)