FSQHLTOfflineSource_cfi Namespace Reference

Functions
def	getDoublePFJet
def	getFSQAll
def	getFSQHI
def	getHighMultVPSet
def	getPTAveVPSet
def	getSinglePFJet
def	getZeroBias_SinglePixelTrackVPSet

Variables
string	fsqdirname = "HLT/FSQ/"
tuple	fsqHLTOfflineSource
string	processName = "HLT"

Function Documentation

def FSQHLTOfflineSource_cfi.getDoublePFJet	(		thresholds,
			flavour = None,
			etaMin = -1,
			srcType = "genJets"
	)

Definition at line 466 of file FSQHLTOfflineSource_cfi.py.

References SiStripPI.max, min(), and str.

Referenced by getFSQAll().

def getDoublePFJet(thresholds, flavour=None, etaMin=-1, srcType="genJets" ):
    if srcType == "genJets":
        inputCol = cms.InputTag("ak4GenJets")
        handlerType = "FromRecoCandidate"
        label = srcType
    elif srcType == "ak4PFJetsCHS":
        inputCol = cms.InputTag("ak4PFJetsCHS")
        handlerType = "RecoPFJetWithJEC"
        label = srcType
    elif srcType == "hlt":
        inputCol = cms.InputTag("S")
        handlerType = "FromHLT"
        label = srcType
    else:
        raise Exception("Whooops!")

    combinedObjectSortCriteria = "at(0).pt + at(1).pt"
    combinedObjectSelection = "1 == 1"
    if flavour != None and "FB" in flavour :
        combinedObjectSortCriteria = "("+combinedObjectSortCriteria+")*(  ? at(0).eta*at(1).eta < 0 ? 1 : 0 )"
        combinedObjectSelection = "at(0).eta*at(1).eta < 0"
        
    if etaMin == None:
        etaMin = -1

    ret=cms.VPSet()
    for t in thresholds:
        partialPathName = "HLT_DiPFJet"+ str(t)+"_"
        if flavour != None:
            partialPathName += flavour+"_"
        partialPathName += "v"

        marginLow = max(t-t//2, 15)
        ptBinLow  = max(t-marginLow,0)
        marginHigh =  min(max(t//3, 15), 50)
        ptBinHigh = t+marginHigh
        ptBins = min(100, ptBinHigh-ptBinLow)
        fromJets =  cms.PSet(
            triggerSelection = cms.string(partialPathName+"*"),
            handlerType = cms.string(handlerType),
            PFJetCorLabel        = cms.InputTag("ak4PFL1FastL2L3Corrector"),
            inputCol = inputCol,
            partialPathName = cms.string(partialPathName),
            partialFilterName  = cms.string("hltDoublePFJet"),
            dqmhistolabel  = cms.string(label),
            mainDQMDirname = cms.untracked.string(fsqdirname),
            singleObjectsPreselection = cms.string("abs(eta) < 5.5 && abs(eta) > " + str(etaMin) ),
            #singleObjectsPreselection = cms.string("pt > 15 && abs(eta) < 5.5 && abs(eta) > " + str(etaMin) ),
            singleObjectDrawables =  cms.VPSet(),
            combinedObjectSelection =  cms.string(combinedObjectSelection),
            combinedObjectSortCriteria = cms.string(combinedObjectSortCriteria),
            combinedObjectDimension = cms.int32(2),
            combinedObjectDrawables =  cms.VPSet(
                    #cms.PSet (name = cms.string("ptLead"), expression = cms.string("max(at(0).pt, at(1).pt())"), 
                    #          bins = cms.int32(ptBins), min = cms.double(ptBinLow), max = cms.double(ptBinHigh)),
                    #cms.PSet (name = cms.string("ptSublead"), expression = cms.string("min(at(0).pt, at(1).pt())"), 
                    #          bins = cms.int32(ptBins), min = cms.double(ptBinLow), max = cms.double(ptBinHigh)),
                    cms.PSet (name = cms.string("ptMostFwd"), 
                              expression = cms.string("? at(0).eta > at(1).eta ? at(0).pt : at(1).pt"), 
                              bins = cms.int32(ptBins), min = cms.double(ptBinLow), max = cms.double(ptBinHigh)),
                    cms.PSet (name = cms.string("ptMostBkw"), 
                              expression = cms.string("? at(0).eta > at(1).eta ? at(1).pt : at(0).pt"), 
                              bins = cms.int32(ptBins), min = cms.double(ptBinLow), max = cms.double(ptBinHigh)),

                    cms.PSet (name = cms.string("etaMostFwd"), 
                              expression = cms.string("? at(0).eta > at(1).eta ? at(0).eta : at(1).eta"), 
                              bins = cms.int32(52), min = cms.double(-5.2), max = cms.double(5.2)),
                    cms.PSet (name = cms.string("etaMostBkw"), 
                              expression = cms.string("? at(0).eta > at(1).eta ? at(1).eta : at(0).eta"), 
                              bins = cms.int32(52), min = cms.double(-5.2), max = cms.double(5.2)),
                    #cms.PSet (name = cms.string("pt_nominator"), expression = cms.string("min(at(0).pt, at(1).pt)"),
                    #         bins = cms.int32(ptBins), min = cms.double(ptBinLow), max = cms.double(ptBinHigh)  )
            )
        )
        #fromJets.triggerSelection = cms.string("HLT_ZeroBias_v*")
        ret.append(fromJets)
        if srcType != "hlt":
            fromJets.combinedObjectDrawables.append( cms.PSet(name = cms.string("ptm_nominator"), 
                                                              expression = cms.string("min(at(0).pt, at(1).pt)"),
                                                              bins = cms.int32(ptBins), 
                                                              min = cms.double(ptBinLow), 
                                                              max = cms.double(ptBinHigh)  ))

            fromJetsDenom  = fromJets.clone()
            fromJetsDenom.triggerSelection = cms.string("HLT_ZeroBias_v*")
            fromJetsDenom.singleObjectDrawables =  cms.VPSet()
            fromJetsDenom.combinedObjectDrawables =  cms.VPSet(
                cms.PSet (name = cms.string("ptm_denominator"), expression = cms.string("min(at(0).pt, at(1).pt)"),
                          bins = cms.int32(ptBins), min = cms.double(ptBinLow), max = cms.double(ptBinHigh)  )
            )
            ret.append(fromJetsDenom)
            

    return ret