413 def makeAlgs (self, config) :
414
415 if config.dataType() is DataType.Data:
416
417 return
418 log = logging.getLogger('makeGeneratorAnalysisSequence')
419
420 if self.runNumber is None:
421 self.runNumber = config.runNumber()
422
423 if self.saveCutBookkeepers and not self.runNumber:
424 raise ValueError ("invalid run number: " + str(self.runNumber))
425
426
427 if self.saveCutBookkeepers:
428 alg = config.createAlgorithm('CP::AsgCutBookkeeperAlg', 'CutBookkeeperAlg')
429 alg.RootStreamName = self.streamName
430 alg.runNumber = self.runNumber
431 if self.cutBookkeepersSystematics is None:
432 alg.enableSystematics = not config.noSystematics()
433 else:
434 alg.enableSystematics = self.cutBookkeepersSystematics
435 if self.histPattern:
436 alg.histPattern = self.histPattern
437 config.addPrivateTool( 'truthWeightTool', 'PMGTools::PMGTruthWeightTool' )
438
439
440 alg = config.createAlgorithm( 'CP::PMGTruthWeightAlg', 'PMGTruthWeightAlg' )
441 config.addPrivateTool( 'truthWeightTool', 'PMGTools::PMGTruthWeightTool' )
442 alg.decoration = 'generatorWeight_%SYS%'
443 config.addOutputVar ('EventInfo', 'generatorWeight_%SYS%', 'weight_mc')
444
445 if self.detailedPDFinfo:
446 alg = config.createAlgorithm( 'CP::PDFinfoAlg', 'PDFinfoAlg', reentrant=True )
447 for var in ["PDFID1","PDFID2","PDGID1","PDGID2","Q","X1","X2","XF1","XF2"]:
448 config.addOutputVar ('EventInfo', var, 'PDFinfo_' + var, noSys=True)
449
450 if self.doPDFReweighting:
451 alg = config.createAlgorithm( 'CP::PDFReweightAlg', 'PDFReweightAlg', reentrant=True )
452
453 for pdf_set in self.outPDFName:
454 config.addOutputVar('EventInfo', f'PDFReweightSF_{pdf_set.replace("/", "_")}',
455 f'PDFReweightSF_{pdf_set.replace("/", "_")}', noSys=True)
456
457
458 if self.doHFProdFracReweighting:
459 generatorInfo = config.flags.Input.GeneratorsInfo
460 log.info(f"Loaded generator info: {generatorInfo}")
461
462 DSID = "000000"
463
464 if not generatorInfo:
465 log.warning("No generator info found.")
466 DSID = "000000"
467 elif isinstance(generatorInfo, dict):
468 if "Pythia8" in generatorInfo:
469 DSID = "410470"
470 elif "Sherpa" in generatorInfo and "2.2.8" in generatorInfo["Sherpa"]:
471 DSID = "421152"
472 elif "Sherpa" in generatorInfo and "2.2.10" in generatorInfo["Sherpa"]:
473 DSID = "700122"
474 elif "Sherpa" in generatorInfo and "2.2.11" in generatorInfo["Sherpa"]:
475 log.warning("HF production fraction reweighting is not configured for Sherpa 2.2.11. Using weights for Sherpa 2.2.10 instead.")
476 DSID = "700122"
477 elif "Sherpa" in generatorInfo and "2.2.12" in generatorInfo["Sherpa"]:
478 log.warning("HF production fraction reweighting is not configured for Sherpa 2.2.12. Using weights for Sherpa 2.2.10 instead.")
479 DSID = "700122"
480 elif "Sherpa" in generatorInfo and "2.2.14" in generatorInfo["Sherpa"]:
481 log.warning("HF production fraction reweighting is not configured for Sherpa 2.2.14. New weights need to be calculated.")
482 DSID = "000000"
483 elif "Sherpa" in generatorInfo and "2.2.1" in generatorInfo["Sherpa"]:
484 DSID = "410250"
485 elif "Herwig7" in generatorInfo and "7.1.3" in generatorInfo["Herwig7"]:
486 DSID = "411233"
487 elif "Herwig7" in generatorInfo and "7.2.1" in generatorInfo["Herwig7"]:
488 DSID = "600666"
489 elif "Herwig7" in generatorInfo and "7." in generatorInfo["Herwig7"]:
490 DSID = "410558"
491 elif "amc@NLO" in generatorInfo:
492 DSID = "410464"
493 else:
494 log.warning(f"HF production fraction reweighting is not configured for this generator: {generatorInfo}")
495 log.warning("New weights need to be calculated.")
496 DSID = "000000"
497 else:
498 log.warning("Failed to determine generator from metadata")
499 DSID = "000000"
500
501 log.info(f"Using HF production fraction weights calculated using DSID {DSID}")
502 if DSID == "000000":
503 log.warning("HF production fraction reweighting will return dummy weights of 1.0")
504
505 alg = config.createAlgorithm( 'CP::SysTruthWeightAlg', 'SysTruthWeightAlg' + self.streamName )
506 config.addPrivateTool( 'sysTruthWeightTool', 'PMGTools::PMGHFProductionFractionTool' )
507 alg.decoration = 'prodFracWeight_%SYS%'
508 alg.TruthParticleContainer = self.truthParticleContainer
509 alg.sysTruthWeightTool.ShowerGenerator = DSID
510 config.addOutputVar ('EventInfo', 'prodFracWeight_%SYS%', 'weight_HF_prod_frac')
511
1.14.0