359 def makeAlgs (self, config) :
360
361 log = logging.getLogger('RScanJetAnalysisConfig')
362
363 jetCollectionName=self.jetCollection
364 if(self.jetCollection=="AnalysisJets") :
365 jetCollectionName="AntiKt4EMPFlowJets"
366 if(self.jetCollection=="AnalysisLargeRJets") :
367 jetCollectionName="AntiKt10LCTopoTrimmedPtFrac5SmallR20Jets"
368
369 if not config.isPhyslite() or self.recalibratePhyslite:
370 if self.jetInput not in ["LCTopo", "HI"]:
371 raise ValueError(
372 "Unsupported input type '{0}' for R-scan jets!".format(self.jetInput) )
373
374 calibToolName = 'JetCalibTool_' + jetCollectionName[:-4]
375 calibTool = config.createPublicTool( 'JetCalibrationTool', calibToolName )
376 calibTool.JetCollection = jetCollectionName[:-4]
377 if self.jetInput=="LCTopo":
378 calibTool.ConfigFile = \
379 "JES_MC16Recommendation_Rscan{0}LC_Feb2022_R21.config".format(self.radius)
380 if config.dataType() is DataType.Data:
381 calibTool.CalibSequence = "JetArea_Residual_EtaJES_GSC_Insitu"
382 else:
383 calibTool.CalibSequence = "JetArea_Residual_EtaJES_GSC_Smear"
384 elif self.jetInput=="HI":
385 calibTool.ConfigFile = \
386 "JES_MC16_HI_Jan2021_5TeV.config"
387 if config.dataType() is DataType.Data:
388 calibTool.CalibSequence = "EtaJES_Insitu"
389 else:
390 calibTool.CalibSequence = "EtaJES"
391 calibTool.IsData = (config.dataType() is DataType.Data)
392
393 alg = config.createAlgorithm( 'CP::JetCalibrationAlg', 'JetCalibrationAlg' )
394 alg.HIsetup = self.jetInput == "HI"
395 alg.calibrationTool = f'{calibTool.getType()}/{calibTool.getName()}'
396 alg.jets = config.readName (self.containerName)
397
398 log.warning("Uncertainties for R-Scan jets are not yet released!")
399
400