ATLAS Offline Software
Loading...
Searching...
No Matches
JetCalibToolsConfig Namespace Reference

Functions

 getJetCalibTool (jetdef, context, data_type, calibseq="", rhoname="", pvname="PrimaryVertices", gscdepth="auto")
 defineJetCalibTool (jetcollection, context, configfile, calibarea, calibseq, data_type, rhoname, pvname, gscdepth)
 getJetCalibToolPrereqs (modspec, jetdef)
 getCalibSpecsFromString (jetdef, modspec)
 getJetCalibToolFromString (jetdef, modspec)

Variables

 jetcaliblog = Logging.logging.getLogger('JetCalibToolsConfig')
list all = ['getJetCalibTool']
dict pflowcontexts
dict topocontexts
dict ufocontexts
dict rscanlc2
dict rscanlc6
dict fatjetcontexts
dict af2configs
dict calibcontexts
list hasInSitu = ["AntiKt4LCTopo", "AntiKt4EMTopo", "AntiKt4EMPFlow", "AntiKt4EMPFlow_noElectrons", "AntiKt4EMPFlow_noMuons", "AntiKt4EMPFlow_noLeptons", "AntiKt4EMPFlow_inclMuons","AntiKt4EMPFlow_tauSeedEleRM","TrigAntiKt4EMTopo"]

Function Documentation

◆ defineJetCalibTool()

defineJetCalibTool ( jetcollection,
context,
configfile,
calibarea,
calibseq,
data_type,
rhoname,
pvname,
gscdepth )

Definition at line 190 of file JetCalibToolsConfig.py.

190def defineJetCalibTool(jetcollection, context, configfile, calibarea, calibseq, data_type, rhoname, pvname, gscdepth):
191 # Abbreviate the calib sequence
192 calibseqshort = ''.join([ step[0] for step in calibseq.split('_') ])
193 toolname = "jetcalib_{0}_{1}_{2}".format(jetcollection,calibseqshort,context)
194 #
195 from AthenaConfiguration.ComponentFactory import CompFactory
196 jct = CompFactory.JetCalibrationTool(toolname,
197 JetCollection = jetcollection,
198 ConfigFile = configfile,
199 CalibArea = calibarea,
200 CalibSequence = calibseq,
201 IsData = (data_type == "data"),
202 RhoKey = rhoname,
203 PrimaryVerticesContainerName = pvname,
204 GSCDepth = gscdepth
205 )
206 return jct
207
208# This method extends the basic config getter to specify the requisite jet
209# moments or other inputs

◆ getCalibSpecsFromString()

getCalibSpecsFromString ( jetdef,
modspec )

Definition at line 241 of file JetCalibToolsConfig.py.

241def getCalibSpecsFromString(jetdef, modspec):
242 calibseq = ""
243 # for the LCTopo EleRM context used for EleRM taus, we need alternative EventShape container
244 rhoname = "EleRM_Kt4LCTopoOriginEventShape" if jetdef.context == "EleRM" else "auto"
245 pvname = "PrimaryVertices"
246 gscdepth = "auto"
247 calibspecs = modspec.split(':')
248 # Probably want data_type to come from elsewhere
249 calibcontext, data_type = calibspecs[:2]
250 if len(calibspecs)>2: calibseq = calibspecs[2]
251 if len(calibspecs)>3: rhoname = calibspecs[3]
252 if len(calibspecs)>4: pvname = calibspecs[4]
253 if len(calibspecs)>5: gscdepth = calibspecs[5]
254
255 return calibcontext, data_type, calibseq, rhoname, pvname, gscdepth
256
257# This method instantiates the JetCalibTool given the input mod specification

◆ getJetCalibTool()

getJetCalibTool ( jetdef,
context,
data_type,
calibseq = "",
rhoname = "",
pvname = "PrimaryVertices",
gscdepth = "auto" )

Definition at line 117 of file JetCalibToolsConfig.py.

117def getJetCalibTool(jetdef, context, data_type, calibseq = "", rhoname = "", pvname = "PrimaryVertices", gscdepth = "auto"):
118
119 jetcollection = jetdef.basename
120
121 # In principle we could autoconfigure
122 if data_type not in ['data','mc','afii']:
123 jetcaliblog.error("JetCalibConfig accepts data_type values: 'data', 'mc', 'afii'")
124 raise ValueError("Unsupported data_type provided: '{0}".format(data_type))
125
126 if jetcollection.endswith("Jets"):
127 jetcaliblog.error("Jet collection '{0}'should be specified w/o 'Jets' in the name.".format(jetcollection))
128 raise ValueError("Bad jet collection formatting in getJetCalibTool.")
129
130 jetcaliblog.debug("Preparing calibration for {0}, in context {1} on sample type {2}".format(jetcollection,context,data_type))
131
132 jetcontexts = calibcontexts[jetcollection]
133
134 if jetcollection == "AntiKt4EMTopo" and context == "T0":
135 if jetdef._cflags.GeoModel.Run >= LHCPeriod.Run4:
136 context = "HLLHC"
137
138 try:
139 configfile, calibarea, calibseq_def = jetcontexts[context]
140 _calibseq = calibseq if calibseq else calibseq_def
141
142 # Check that the calib sequence requests something sensible for the in situ calibration
143 # Leave other checks for the tool code.
144 # Might need to specialise if we decide MC trigger jets should also have in situ.
145 if _calibseq.endswith("Insitu"):
146 if data_type == 'data':
147 if jetcollection not in hasInSitu:
148 raise ValueError("In situ calibration does not exist for {0}, context {1}".format(jetcollection,context))
149 else:
150 raise ValueError("In situ calibration requested for MC on {0}, context {1}".format(jetcollection,context))
151 _configfile = configfile
152 if data_type == "afii" and jetcollection in af2configs.keys():
153 _configfile = af2configs[jetcollection]
154 _data_type = data_type
155 # Most likely an oversight, but R20/21 JetRecCalibrationFinder did not set the data flag
156 # This affects the residual correction by scaling mu.
157 # We should revert this later on, set up now for validation purposes
158 if context == "T0":
159 _data_type = "data"
160 _pvname = ""
161 if "Residual" in _calibseq or "GSC" in _calibseq and gscdepth!="EM3" or "LargeRDNN" in _calibseq:
162 _pvname = pvname
163 # HACK: For Trigger HI UPC who want to do pflow but avoid track-GSC
164 if context == "TrigHIUPC":
165 _pvname = pvname
166 # HACK: For testing while we don't have finalised calibrations for trigger PF jets
167 _jetcollection = jetcollection
168 if "PFlow" in jetcollection and (context=="TrigSoftDrop" or context=="TrigHLTSoftDrop"):
169 _jetcollection = jetcollection.replace("EMPFlow","UFO")
170 # Use standard calibrations for certain jet collections
171 # with minor changes to their constituents that should
172 # not invalidate the JES
173 ignore_substr = [
174 '_noElectrons', '_noMuons', '_noLeptons', '_inclMuons',
175 '_tauSeedEleRM',
176 'ByVertex', 'CustomVtxGNN', 'CustomVtx',
177 ]
178 for text in ignore_substr:
179 if text in jetcollection :
180 _jetcollection = _jetcollection.replace(text,"")
181
182 return defineJetCalibTool(_jetcollection, context, _configfile, calibarea, _calibseq, _data_type, rhoname, _pvname, gscdepth)
183 except KeyError as e:
184 jetcaliblog.error("Context '{0}' not found for jet collection '{1}'".format(context,jetcollection))
185 jetcaliblog.error("Options are '{0}".format(','.join(jetcontexts.keys())))
186 raise e
187 return None
188
189# This method actually sets up the tool

◆ getJetCalibToolFromString()

getJetCalibToolFromString ( jetdef,
modspec )

Definition at line 258 of file JetCalibToolsConfig.py.

258def getJetCalibToolFromString(jetdef, modspec):
259 calibcontext, data_type, calibseq, rhoname, pvname, gscdepth = getCalibSpecsFromString(jetdef, modspec)
260 return getJetCalibTool(jetdef,calibcontext,data_type,calibseq,rhoname,pvname,gscdepth)

◆ getJetCalibToolPrereqs()

getJetCalibToolPrereqs ( modspec,
jetdef )

Definition at line 210 of file JetCalibToolsConfig.py.

210def getJetCalibToolPrereqs(modspec,jetdef):
211 calibcontext, data_type, calibseq, rhoname, pvname, gscdepth = getCalibSpecsFromString(jetdef, modspec)
212 if calibseq=="":
213 cfg, calibarea, calibseq = calibcontexts[jetdef.basename][calibcontext]
214
215 # For now, only dependent on calibseq -- can ignore Insitu, which is
216 # added when getting the concrete tool
217 prereqs = []
218 prereqs.append("mod:ConstitFourMom")
219 if "JetArea" in calibseq: # Will not insert a prefix here
220 if calibcontext.startswith("Trig"): prereqs.append("input:HLT_EventDensity")
221 elif pvname == "PrimaryVertices_initial": prereqs.append("input:EventDensityCustomVtxGNN")
222 elif pvname != "PrimaryVertices": prereqs.append("input:EventDensityCustomVtx")
223 else: prereqs.append(inputsFromContext("EventDensity")(jetdef))
224 if "GSC" in calibseq:
225 prereqs += ["mod:CaloEnergies"]
226 if calibcontext != "TrigRun2": # No track/MS GSC for trigger w/o FTK
227 prereqs += ["mod:TrackMoments",
228 "ghost:MuonSegment"]
229 if "GNNC" in calibseq:
230 prereqs += ["mod:CaloQuality","mod:TrackMoments"]
231 if "CombinedMass" in calibcontext:
232 prereqs += ["mod:TrackSumMoments"]
233 if "LargeRDNN" in calibseq:
234 prereqs += ["mod:CaloEnergiesLargeR","mod:ConstitFrac","mod:groomMRatio","mod:Width",
235 "mod:nsubjettiness","mod:nsubjettinessR","mod:ktsplitter","mod:ecorr",
236 "mod:ecorrR","mod:qw"]
237 jetcaliblog.debug("Prereqs for calibseq '{0}': {1}".format(calibseq,str(prereqs)))
238 return prereqs
239
240# This method translates the mod specification string into calibration specifications

Variable Documentation

◆ af2configs

dict JetCalibToolsConfig.af2configs
Initial value:
= {
"AntiKt4EMPFlow": "JES_MC16Recommendation_AFII_PFlow_Apr2019_Rel21.config",
"AntiKt4GPFlow": "JES_MC16Recommendation_AFII_PFlow_Apr2019_Rel21.config",
"AntiKt4EMTopo": "JES_MC16Recommendation_AFII_EMTopo_Apr2019_Rel21.config",
"AntiKt4LCTopo": "JES_MC16Recommendation_AFII_EMTopo_Apr2019_Rel21.config",
}

Definition at line 70 of file JetCalibToolsConfig.py.

◆ all

list JetCalibToolsConfig.all = ['getJetCalibTool']

Definition at line 19 of file JetCalibToolsConfig.py.

◆ calibcontexts

dict JetCalibToolsConfig.calibcontexts

Definition at line 77 of file JetCalibToolsConfig.py.

◆ fatjetcontexts

dict JetCalibToolsConfig.fatjetcontexts
Initial value:
= {
"CombinedMass": ("JES_MC16recommendation_FatJet_Trimmed_JMS_comb_17Oct2018.config","00-04-82","EtaJES_JMS"),
"CaloMass": ("JES_MC16recommendation_FatJet_Trimmed_JMS_calo_12Oct2018.config","00-04-82","EtaJES_JMS"),
"TAMass": ("JES_MC16recommendation_FatJet_Trimmed_JMS_TA_12Oct2018.config","00-04-82","EtaJES_JMS"),
"SoftDrop": ("JES_MC20PreRecommendation_R10_UFO_CSSK_SoftDrop_JMS_R21Insitu_26Nov2024.config","00-04-83","EtaJES_JMS"),
"TrigUngroomed": ("JES_Full2012dataset_Rscan_June2014.config","00-04-77","JetArea_EtaJES"),
"TrigTrimmed": ("JES_MC15recommendation_FatJet_June2015_PtFrac4.config","00-04-82","EtaJES_JMS"),
"TrigSoftDrop": ("JES_MC16recommendation_R10_UFO_CSSK_SoftDrop_JMS_01April2020.config","00-04-82","EtaJES_JMS"),
"TrigHLTSoftDrop":("JES_JMS_MC21_R10_PFlow_CSSK_SoftDrop_Dec2023_Trigger.config", "00-04-83","EtaJES_JMS"),
"LargeRDNN": ("JES_JMS_MC20dnnc_R10_UFO_CSSK_SoftDrop_Jun2023.config","00-04-83","LargeRDNN"),
}

Definition at line 57 of file JetCalibToolsConfig.py.

◆ hasInSitu

list JetCalibToolsConfig.hasInSitu = ["AntiKt4LCTopo", "AntiKt4EMTopo", "AntiKt4EMPFlow", "AntiKt4EMPFlow_noElectrons", "AntiKt4EMPFlow_noMuons", "AntiKt4EMPFlow_noLeptons", "AntiKt4EMPFlow_inclMuons","AntiKt4EMPFlow_tauSeedEleRM","TrigAntiKt4EMTopo"]

Definition at line 111 of file JetCalibToolsConfig.py.

◆ jetcaliblog

JetCalibToolsConfig.jetcaliblog = Logging.logging.getLogger('JetCalibToolsConfig')

Definition at line 17 of file JetCalibToolsConfig.py.

◆ pflowcontexts

dict JetCalibToolsConfig.pflowcontexts
Initial value:
= {
"T0":("JES_MC15cRecommendation_PFlow_Aug2016_rel21.config","00-04-77","JetArea_Residual_EtaJES"),
# Omit smearing, to avoid any efficiency loss
"AnalysisLatest":("PreRec_R22_PFlow_ResPU_EtaJES_GSC_February23_230215.config","00-04-82","JetArea_Residual_EtaJES_GSC_Insitu"),
"TrigLS2":("JES_MC16Recommendation_Consolidated_PFlow_Apr2019_Rel21_Trigger.config","00-04-82","JetArea_Residual_EtaJES_GSC_Insitu"),
"Trigger":("JES_MC16Recommendation_Consolidated_PFlow_30May2022_Rel22_Trigger.config","00-04-82","JetArea_Residual_EtaJES_GSC_Insitu"),
"TrigR22Prerec":("PreRec_R22_PFlow_ResPU_EtaJES_GSC_February23_230215.config", "00-04-82", "JetArea_Residual_EtaJES_GSC_Insitu"),
"TrigHIUPC" : ("JES_MC16Recommendation_LowMu1718_MCJES_GSC_nTrkOn_PFlow_Sep2023_Rel21.config", "00-04-82", "EtaJES_GSC"),
}

Definition at line 24 of file JetCalibToolsConfig.py.

◆ rscanlc2

dict JetCalibToolsConfig.rscanlc2
Initial value:
= {
"RScanLatest":("JES_MC16Recommendation_Rscan2LC_22Feb2018_rel21.config","00-04-81","JetArea_Residual_EtaJES_GSC_Insitu")
}

Definition at line 49 of file JetCalibToolsConfig.py.

◆ rscanlc6

dict JetCalibToolsConfig.rscanlc6
Initial value:
= {
"RScanLatest":("JES_MC16Recommendation_Rscan6LC_22Feb2018_rel21.config","00-04-81","JetArea_Residual_EtaJES_GSC_Insitu")
}

Definition at line 53 of file JetCalibToolsConfig.py.

◆ topocontexts

dict JetCalibToolsConfig.topocontexts
Initial value:
= {
"T0":("JES_MC15cRecommendation_May2016_rel21.config","00-04-77","JetArea_Residual_EtaJES"),
# Omit smearing, to avoid any efficiency loss
"AnalysisLatest":("JES_MC16Recommendation_Consolidated_EMTopo_Apr2019_Rel21.config","00-04-82","JetArea_Residual_EtaJES_GSC_Insitu"),
"TrigRun2":("JES_MC15cRecommendation_May2016_Trigger.config","00-04-77","JetArea_EtaJES_GSC_Insitu"),
"TrigRun2GSC":("JES_data2016_data2015_Recommendation_Dec2016_rel21.config","00-04-77","JetArea_EtaJES_GSC_Insitu"),
"TrigLS2":("JES_MC16Recommendation_Consolidated_EMTopo_Apr2019_Rel21_Trigger.config","00-04-82","JetArea_Residual_EtaJES_GSC_Insitu"),
"Trigger":("JES_MC16Recommendation_Consolidated_EMTopo_Apr2019_Rel21_Trigger.config","00-04-82","JetArea_Residual_EtaJES_GSC_Insitu"),
"HLLHC":("HLLHC/JES_MC16_HLLHC_Aug2021_rel21.config","00-04-82","JetArea_Residual_EtaJES"),
}

Definition at line 34 of file JetCalibToolsConfig.py.

◆ ufocontexts

dict JetCalibToolsConfig.ufocontexts
Initial value:
= {
"T0":("Consolidated_R22_CSSKUFO_ResPU_EtaJES_GNNC_20231103.config", "00-04-83", "Residual_EtaJES_GNNC"),
}

Definition at line 45 of file JetCalibToolsConfig.py.