Python QcdD3PDMakerFlags.doRecJet Examples

Example #1

File: myJets_rel172_1.py Project: rushioda/PIXELVALID_athena

from IOVDbSvc.CondDB import conddb

#---
# Jet Width
WidthT = JetWidthTool("JetWidth")
#---

from JetMomentTools.JetMomentToolsConf import LowEtConstituentsFracTool

ConstMomTool = LowEtConstituentsFracTool("LowEtConstituentsFrac")

#set inputFileType
if rec.readAOD():
    jetFlags.inputFileType = 'AOD'

if QcdD3PDMakerFlags.doRecJet():
    if QcdD3PDMakerFlags.doTopoJet():  # off in default
        #TopoJets (GCW+JES)
        make_StandardJetGetter('AntiKt',
                               0.4,
                               'Topo',
                               doEMCut=False,
                               minPt=10 * GeV,
                               doCalib=True,
                               calibName='EM:GCW_ORIGIN_ETAJES',
                               GCWScheme='Snowmass',
                               addJetQualityMoments=True,
                               writeMoments=True,
                               outputCollectionName='AntiKt4TopoNewJets')
        make_StandardJetGetter('AntiKt',
                               0.6,

Example #2

from BTagging.BTaggingFlags import BTaggingFlags
from JetRec.JetRecFlags import jetFlags
from RecExConfig.RecFlags import rec

BTaggingFlags.Jets = []
#BTaggingFlags.Jets += ['AntiKt4TopoEMJets','AntiKt6TopoEMJets',
#                       'AntiKt4LCTopoJets','AntiKt6LCTopoJets'] # for the default jets
if rec.doTruth():
    BTaggingFlags.Jets += ['AntiKt4Truth']
if rec.doTruth() and QcdD3PDMakerFlags.doRecTruthJet():
    BTaggingFlags.Jets += [
        'AntiKt4TruthNew', 'AntiKt5TruthNew', 'AntiKt6TruthNew',
        'AntiKt4TruthWithNoInt', 'AntiKt6TruthWithNoInt',
        'AntiKt4TruthWithMuNoInt', 'AntiKt6TruthWithMuNoInt'
    ]
if QcdD3PDMakerFlags.doRecJet() and QcdD3PDMakerFlags.doTopoEMJet():
    BTaggingFlags.Jets += ['AntiKt4TopoEMNew', 'AntiKt6TopoEMNew']
if QcdD3PDMakerFlags.doRecJet() and QcdD3PDMakerFlags.doLCTopoJet():
    BTaggingFlags.Jets += [
        'AntiKt4LCTopoNew', 'AntiKt5LCTopoNew', 'AntiKt6LCTopoNew'
    ]

BTaggingFlags.JetsWithInfoPlus = BTaggingFlags.Jets[:]
#BTaggingFlags.CalibrationTag="BTagCalibALL-07-02" # up to now it is not linked to the global flag

BTaggingFlags.CalibrationChannelAliases += [
    'AntiKt4Truth->AntiKt4TopoEM',
    'AntiKt4TruthNew->AntiKt4TopoEM',
    'AntiKt5TruthNew->AntiKt4TopoEM',
    'AntiKt6TruthNew->AntiKt6TopoEM',
    'AntiKt4TruthWithNoInt->AntiKt4TopoEM',

Example #3

File: QcdD3PDMakerPreSetup.calib.py Project: rushioda/PIXELVALID_athena

#==================================
# PrimaryDPDMaker specific options
#==================================

#no skimming
primDPD.ApplySkimming = False

#not make DESD for default
primDPD.MakePrimaryDPD = False

#==================================
# RecExCommon specific options
#==================================

if QcdD3PDMakerFlags.doRecJet() or QcdD3PDMakerFlags.doRecMET(
) or QcdD3PDMakerFlags.doRecTruthJet():
    rec.oldFlagTopSteering.set_Value_and_Lock(True)

#==================================
#Setting AOD flags. Not much needs to be done if you read an ESD 'cause everything's already built.
#See https://twiki.cern.ch/twiki/bin/view/Atlas/UserAnalysisTest#The_AOD_Production_Flags
#==================================
#other settings
jetFlags.noStandardConfig = True

if rec.readAOD() or not (QcdD3PDMakerFlags.doRecJet()
                         or QcdD3PDMakerFlags.doRecMET()):
    jobproperties.CaloRecFlags.Enabled.set_Value_and_Lock(False)
    jobproperties.CaloRecFlags.doCaloCluster.set_Value_and_Lock(False)
    jobproperties.CaloRecFlags.doCaloTopoCluster.set_Value_and_Lock(False)

Example #4

#==================================
# PrimaryDPDMaker specific options
#==================================

#no skimming
primDPD.ApplySkimming = False

#not make DESD for default
primDPD.MakePrimaryDPD = False

#==================================
# RecExCommon specific options
#==================================

if QcdD3PDMakerFlags.doRecJet() or QcdD3PDMakerFlags.doRecMET(
) or QcdD3PDMakerFlags.doRecTruthJet():
    rec.oldFlagTopSteering.set_Value_and_Lock(True)

#==================================
#Setting AOD flags. Not much needs to be done if you read an ESD 'cause everything's already built.
#See https://twiki.cern.ch/twiki/bin/view/Atlas/UserAnalysisTest#The_AOD_Production_Flags
#==================================
#other settings
jetFlags.noStandardConfig = True

#if rec.readAOD() or not (QcdD3PDMakerFlags.doRecJet() or QcdD3PDMakerFlags.doRecMET()):
#  jobproperties.CaloRecFlags.Enabled.set_Value_and_Lock(False)
#  jobproperties.CaloRecFlags.doCaloCluster.set_Value_and_Lock(False)
#  jobproperties.CaloRecFlags.doCaloTopoCluster.set_Value_and_Lock(False)
#else:

Example #5

File: metAlgs_rel166.py Project: rushioda/PIXELVALID_athena

from MissingET.METRefGetter_new import *
from QcdD3PDMaker.QcdD3PDMakerFlags import QcdD3PDMakerFlags

if QcdD3PDMakerFlags.doRecJet(
) and not QcdD3PDMakerFlags.doRecJetOnlyAdditional():
    my_alg = make_METRefAlg(_suffix='_GCW_pt20')
    my_alg.jet_JetInputCollectionKey = 'AntiKt6TopoNewJets'
    my_alg.jet_JetPtCut = 20.0 * GeV
    my_alg.jet_ApplyJetScale = "Yes"
    my_alg.jet_UseJetMomentForScale = True
    my_alg.jet_JetMomentForScale = "GCWJES"
    my_alg.jet_RunSoftJetsTool = True
    my_alg.jet_SoftJetsPtCut = 7.0 * GeV
    my_alg.jet_SoftJetsMaxPtCut = 20.0 * GeV
    my_alg.jet_ApplySoftJetsScale = "No"
    my_alg.jet_calibType = 'H1'
    my_alg.ele_calibType = 'RefCalib'
    my_alg.gamma_calibType = 'EmScale'
    my_alg.cellout_calibType = 'H1'
    my_alg.tau_calibType = 'H1'
    my_alg.cryo_ApplyCorrection = "On"
    my_alg.muon_algorithm = "Staco"
    my_alg.muon_isolationAlg = "dRJet"
    my_alg()

    my_alg3 = make_METRefAlg(_suffix='_LCW_pt20')
    my_alg3.jet_JetInputCollectionKey = 'AntiKt6LCTopoNewJets'
    my_alg3.jet_JetPtCut = 20.0 * GeV
    my_alg3.jet_ApplyJetScale = "No"
    my_alg3.jet_UseJetMomentForScale = True
    my_alg3.jet_JetMomentForScale = "LCJES"