Skip to content

Refactor truth information

Adds truth information about S and H coming from X, and quarks coming from S and H. This information is added as decorations on EventInfo, and retrieved by the tree maker and recorded as new branches:

  • truth_b_fromH_pt
  • truth_b_fromH_eta
  • truth_b_fromH_phi
  • truth_b_fromH_m
  • truth_b_fromS_pt
  • truth_b_fromS_eta
  • truth_b_fromS_phi
  • truth_b_fromS_m
  • truth_H_pt
  • truth_H_eta
  • truth_H_phi
  • truth_H_m
  • truth_S_pt
  • truth_S_eta
  • truth_S_phi
  • truth_S_m

Also returns a View container for the selected truth particles to event store for other algs (DiHiggs, SH) to use:

# Select truth S, H, and children b quarks

if dataType != "data":
    log.info("Adding truth particle info seq")
    cfg.merge(
        TruthParticleInformationAlgCfg(
            flags,
            inputContainerName=containers["inputs"]["truthParticles"],
            outputContainerName=containers["outputs"]["truthParticles"],
        )
    )

# Other algs that uses this information

if do_dihiggs_analysis:
    cfg.merge(
        DiHiggsAnalysisAlgCfg(
            flags,
            SmallJetKey=containers["outputs"]["reco4Jet"].replace("%SYS%", "NOSYS"),
            LargeJetKey=containers["outputs"]["reco10Jet"].replace(
                "%SYS%", "NOSYS"
            ),
            TruthParticleInformationKey=containers["outputs"]["truthParticles"],
            btag_wps=btag_wps,
            vr_btag_wps=vr_btag_wps,
        )
    )

// in DiHiggsAnalysisAlg.h

SG::ReadHandleKey<xAOD::TruthParticleContainer> m_truthParticleInfoInKey{
    this, "TruthParticleInformationKey", "",

related to #21 (closed)

Edited by Frederic Renner

Merge request reports

Loading