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,
from MissingET.METRefGetter_plup import *
from QcdD3PDMaker.QcdD3PDMakerFlags import QcdD3PDMakerFlags
if QcdD3PDMakerFlags.doRecJet():
my_alg1 = make_METRefAlg(_suffix='_default')
my_alg1.jet_JetInputCollectionKey ='AntiKt4LCTopoNewJets'
my_alg1.jet_JetPtCut = 20.0*GeV
my_alg1.jet_ApplyJetScale = "Yes"
my_alg1.jet_UseJetMomentForScale = True
my_alg1.jet_JetMomentForScale = "LCJES"
my_alg1.jet_RunSoftJetsTool = False
my_alg1.jet_SoftJetsPtCut = 10.0*GeV
my_alg1.jet_SoftJetsMaxPtCut = 20.0*GeV
my_alg1.jet_ApplySoftJetsScale = "No"
my_alg1.jet_calibType ='LocHad'
my_alg1.ele_calibType ='RefCalib'
my_alg1.gamma_calibType ='EmScale'
my_alg1.cellout_calibType ='Eflow'
my_alg1.tau_calibType ='ExclRefCalib'
my_alg1.cryo_ApplyCorrection = "Off"
my_alg1.muon_algorithm = "Staco"
my_alg1.muon_isolationAlg = "dRJet"
my_alg1()
my_alg2 = make_METRefAlg(_suffix='_LCW_pt20')
my_alg2.jet_JetInputCollectionKey ='AntiKt4LCTopoNewJets'
my_alg2.jet_JetPtCut = 20.0*GeV
my_alg2.jet_ApplyJetScale = "No"
my_alg2.jet_RunSoftJetsTool = False
my_alg2.jet_calibType ='LocHad'
from JetMomentTools.JetMomentToolsConf import JetMomentsFromCalib
from RecExConfig.RecFlags import rec
from QcdD3PDMaker.QcdD3PDMakerFlags import QcdD3PDMakerFlags
mt = JetMomentsFromCalib()
from JetCalibTools.MakeCalibSequences import alternateCalibSequence
from JetRec.JetAlgConfiguration import checkAndUpdateOptions
#Request the GSC calib sequence, starting from the jet *final scale* (i.e. ETAJES for TopoEMJets)
d = checkAndUpdateOptions(input='Topo')
StructSeq = alternateCalibSequence("CAL:STRUCT1", d)
mt.Calibrator = StructSeq # link the calibrator to the moment tool.
# redefine StructSeq, and set options to save only the last
# 2 moments of the GSC sequence.
StructSeq = StructSeq.CalibToolSequence[-1]
StructSeq.MomentTag = "GSCFactor"
StructSeq.SimpleFactorCorrection = True
for t in StructSeq.CalibToolSequence[:
-1]: #switch off all moments writing but the last.
t.WriteMoments = False
# schedule the alg computing jet moments :
make_JetMomentGetter("AntiKt4TopoEMJets", [mt])
make_JetMomentGetter("AntiKt6TopoEMJets", [mt])
if QcdD3PDMakerFlags.doRecJet():
if QcdD3PDMakerFlags.doTopoEMJet():
make_JetMomentGetter("AntiKt4TopoEMNewJets", [mt])
make_JetMomentGetter("AntiKt6TopoEMNewJets", [mt])
if QcdD3PDMakerFlags.doLCTopoJet():
make_JetMomentGetter("AntiKt4LCTopoNewJets", [mt])
make_JetMomentGetter("AntiKt5LCTopoNewJets", [mt])
make_JetMomentGetter("AntiKt6LCTopoNewJets", [mt])
#==================================
# 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:
from QcdD3PDMaker.QcdD3PDMakerFlags import QcdD3PDMakerFlags
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',
from RecExConfig.RecFlags import rec
from JetRec.JetMomentGetter import make_JetMomentGetter
from JetMomentTools.SetupJetMomentTools import getJetVertexAssociationTool
jvatool = getJetVertexAssociationTool('AntiKt', 0.4, 'Topo')
jvatool.StoreQualityVariables = True
# jvatool.JVFMomentName = "JVF2"
make_JetMomentGetter('AntiKt4TopoEMJets', [jvatool])
make_JetMomentGetter('AntiKt6TopoEMJets', [jvatool])
make_JetMomentGetter('AntiKt4TopoJets', [jvatool])
make_JetMomentGetter('AntiKt6TopoJets', [jvatool])
make_JetMomentGetter('AntiKt4LCTopoJets', [jvatool])
make_JetMomentGetter('AntiKt6LCTopoJets', [jvatool])
make_JetMomentGetter('AntiKt4TowerJets', [jvatool])
make_JetMomentGetter('AntiKt6TowerJets', [jvatool])
from QcdD3PDMaker.QcdD3PDMakerFlags import QcdD3PDMakerFlags
if QcdD3PDMakerFlags.doRecJet():
make_JetMomentGetter('AntiKt4TopoEMNewJets', [jvatool])
make_JetMomentGetter('AntiKt6TopoEMNewJets', [jvatool])
make_JetMomentGetter('AntiKt4TopoNewJets', [jvatool])
make_JetMomentGetter('AntiKt6TopoNewJets', [jvatool])
make_JetMomentGetter('AntiKt4LCTopoNewJets', [jvatool])
make_JetMomentGetter('AntiKt6LCTopoNewJets', [jvatool])
if not rec.readAOD():
make_JetMomentGetter('AntiKt4TowerEMNewJets', [jvatool])
make_JetMomentGetter('AntiKt6TowerEMNewJets', [jvatool])
make_JetMomentGetter('AntiKt4TowerNewJets', [jvatool])
make_JetMomentGetter('AntiKt6TowerNewJets', [jvatool])
# Imports
from JetRec.JetGetters import *
from JetRec.JetRecConf import *
from JetRec.JetRecFlags import jetFlags
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
from AthenaCommon.AppMgr import ServiceMgr as svcMgr
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():
#B-tag setting
from BTagging.BTaggingFlags import BTaggingFlags
BTaggingFlags.Jets += ['AntiKt4TopoNew','AntiKt6TopoNew','AntiKt4TopoEMNew','AntiKt6TopoEMNew',
'AntiKt4TowerNew','AntiKt6TowerNew','AntiKt4TowerEMNew','AntiKt6TowerEMNew',
'AntiKt4LCTopoNew','AntiKt6LCTopoNew']
BTaggingFlags.CalibrationChannelAliases += [
'AntiKt4TopoNew->AntiKt4Topo,AntiKt4H1Topo',
'AntiKt6TopoNew->AntiKt6Topo,AntiKt6H1Topo,AntiKt6H1Tower',
'AntiKt4TopoEMNew->AntiKt4TopoEM,AntiKt4H1Topo',
'AntiKt6TopoEMNew->AntiKt6TopoEM,AntiKt6H1Topo,AntiKt6H1Tower',
'AntiKt4TowerNew->AntiKt4Tower,AntiKt4H1Tower',
'AntiKt6TowerNew->AntiKt6Tower,AntiKt6H1Tower',
'AntiKt4TowerEMNew->AntiKt4Tower,AntiKt4H1Tower',
'AntiKt6TowerEMNew->AntiKt6Tower,AntiKt6H1Tower',
#==================================
# 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)
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"