Exemplo n.º 1
0
                               writeMoments=True,
                               outputCollectionName='AntiKt4TopoNewJets')
        make_StandardJetGetter('AntiKt',
                               0.6,
                               'Topo',
                               doEMCut=False,
                               minPt=10 * GeV,
                               doCalib=True,
                               calibName='EM:GCW_ORIGIN_ETAJES',
                               GCWScheme='Snowmass',
                               addJetQualityMoments=True,
                               writeMoments=True,
                               outputCollectionName='AntiKt6TopoNewJets')

    if QcdD3PDMakerFlags.doTopoEMJet(
    ) and not QcdD3PDMakerFlags.doRecJetOnlyAdditional():
        #TopoJets (EM+JES)
        make_StandardJetGetter('AntiKt',
                               0.4,
                               'Topo',
                               doEMCut=False,
                               minPt=10 * GeV,
                               doCalib=True,
                               calibName='EM:OFFSET_ORIGIN_ETAJES',
                               addJetQualityMoments=True,
                               writeMoments=True,
                               outputCollectionName='AntiKt4TopoEMNewJets')
        make_StandardJetGetter('AntiKt',
                               0.6,
                               'Topo',
                               doEMCut=False,
Exemplo n.º 2
0
from JetCalibTools.JetCalibToolsConf import *
from CaloRec.CaloTowerCmbGetter import CaloTowerCmbGetter
from AthenaCommon.SystemOfUnits import GeV
from RecExConfig.RecFlags import rec
from PrimaryDPDMaker.PrimaryDPDFlags import primDPD
from QcdD3PDMaker.QcdD3PDMakerFlags import QcdD3PDMakerFlags

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

if QcdD3PDMakerFlags.doRecJet():
  #B-tag setting
  include("BTagging/BTagging_LoadTools.py")

  if QcdD3PDMakerFlags.doTopoJet() and not QcdD3PDMakerFlags.doRecJetOnlyAdditional():
    #TopoJets
    make_StandardJetGetter('AntiKt',0.4,'H1Topo', doEMCut=False, minPt=7*GeV, doCalib=True, calibName='H1NumInv', addJetQualityMoments=True, outputCollectionName='AntiKt4H1TopoNewJets')
    make_StandardJetGetter('AntiKt',0.6,'H1Topo', doEMCut=False, minPt=7*GeV, doCalib=True, calibName='H1NumInv', addJetQualityMoments=True, outputCollectionName='AntiKt6H1TopoNewJets')

  if QcdD3PDMakerFlags.doTopoJet():
    #TopoJets (EM)
    make_StandardJetGetter('AntiKt',0.6,'H1Topo', minPt=4*GeV, calibTools=[], inputSuff="EMscale", outputCollectionName='AntiKt6H1TopoEMJets')

  if not rec.readAOD():
    if QcdD3PDMakerFlags.doTowerJet() and not QcdD3PDMakerFlags.doRecJetOnlyAdditional():
      #TopoTowerJets
      make_StandardJetGetter('AntiKt',0.4,'H1TopoTower', doEMCut=False, minPt=7*GeV, doCalib=True, calibName='H1NumInv', addJetQualityMoments=True, outputCollectionName='AntiKt4H1TowerNewJets', GCWScheme='Snowmass')
      make_StandardJetGetter('AntiKt',0.6,'H1TopoTower', doEMCut=False, minPt=7*GeV, doCalib=True, calibName='H1NumInv', addJetQualityMoments=True, outputCollectionName='AntiKt6H1TowerNewJets', GCWScheme='Snowmass')

    if QcdD3PDMakerFlags.doGhostJet():