CP::PartonHistoryUtils Namespace Reference

Functions
const xAOD::TruthParticle *	findAfterFSR (const xAOD::TruthParticle *particle)
	Return particle after FSR (before the decay vertex).
bool	isChildOf (const xAOD::TruthParticle *parent, const xAOD::TruthParticle *potentialChild)
bool	isBrokenTop (const xAOD::TruthParticle *particle)
	Looking for tops without children -> must be broken.
bool	isAfterFSR (const xAOD::TruthParticle *particle)
	Determine whether particle is afterFSR.
bool	hasParticleIdenticalParent (const xAOD::TruthParticle *particle)
	Return true when particle is a top before FSR.
bool	hasParentPdgId (const xAOD::TruthParticle *particle, int PdgId)
bool	hasParentPdgId (const xAOD::TruthParticle *particle)
bool	hasParentAbsPdgId (const xAOD::TruthParticle *particle, int absPdgId)
bool	hasIdenticalChild (const xAOD::TruthParticle *particle)
bool	isQuarkFromPDF (const xAOD::TruthParticle *particle)
std::string	getretrievalstring (const std::string &parent, const std::string &suffix, const std::string &postfix)

Function Documentation

findAfterFSR()

const xAOD::TruthParticle * CP::PartonHistoryUtils::findAfterFSR

(

const xAOD::TruthParticle *

particle

)

Return particle after FSR (before the decay vertex).

Definition at line 93 of file PartonHistoryUtils.cxx.

                                                                  {
  bool isAfter(false);
  const int pdgId = p->pdgId();
  int i = 0;
  while (!isAfter) {
    i = 0;
    for (size_t j = 0; j < p->nChildren(); j++) {
      const xAOD::TruthParticle* tmp_children = p->child(j);
      if (tmp_children && tmp_children->pdgId() == pdgId) {
        p = p->child(j);
        i++;
        break;
      }  // if
    }  // for
    if (i == 0)
      isAfter = true;
  }  // while
  return p;
}

getretrievalstring()

std::string CP::PartonHistoryUtils::getretrievalstring	(	const std::string &	parent,
		const std::string &	suffix,
		const std::string &	postfix )

Definition at line 122 of file PartonHistoryUtils.cxx.

                                                         {
  std::string baseprefix = "MC_";
  return baseprefix + (parent.empty() ? "" : parent + "_") + suffix + postfix;
}

hasIdenticalChild()

bool CP::PartonHistoryUtils::hasIdenticalChild

(

const xAOD::TruthParticle *

particle

)

Definition at line 31 of file PartonHistoryUtils.cxx.

                                                   {
  // Checks if the given particle has at least one child with an identical PDG
  // ID.
  bool check = false;
  for (size_t i = 0; i < p->nChildren(); i++) {
    if (p->child(i) && p->child(i)->pdgId() == p->pdgId())
      check = true;
  }
  return check;
}

hasParentAbsPdgId()

bool CP::PartonHistoryUtils::hasParentAbsPdgId	(	const xAOD::TruthParticle *	particle,
		int	absPdgId )

Definition at line 42 of file PartonHistoryUtils.cxx.

                                                                 {
  // Checks if the parent of the given particle has a specific absolute PDG ID.
  return p->parent(0) && p->parent(0)->absPdgId() == absPdgId;
}

hasParentPdgId() [1/2]

bool CP::PartonHistoryUtils::hasParentPdgId

(

const xAOD::TruthParticle *

particle

)

Definition at line 26 of file PartonHistoryUtils.cxx.

                                                {
  // Checks if the parent of the given particle has a specific PDG ID.
  return hasParentPdgId(p, p->pdgId());
}

hasParentPdgId() [2/2]

bool CP::PartonHistoryUtils::hasParentPdgId	(	const xAOD::TruthParticle *	particle,
		int	PdgId )

Definition at line 17 of file PartonHistoryUtils.cxx.

                                                           {
  // Checks if the particle or its parent are nullptr.
  if (!p || !p->parent(0)) {
    return false;
  }
  // Checks if the parent of the given particle has a specific PDG ID.
  return p->parent(0)->pdgId() == pdgId;
}

hasParticleIdenticalParent()

bool CP::PartonHistoryUtils::hasParticleIdenticalParent

(

const xAOD::TruthParticle *

particle

)

Return true when particle is a top before FSR.

Definition at line 47 of file PartonHistoryUtils.cxx.

                                                            {
  // Checks if particle and any of ity parents are identical.
  bool skipit(false);
  for (size_t i = 0; i < p->nParents(); i++) {
    const xAOD::TruthParticle* parent = p->parent(i);
    if (parent && parent->pdgId() == p->pdgId()) {
      skipit = true;
      break;
    }  // if
  }  // for
  return skipit;
}

isAfterFSR()

bool CP::PartonHistoryUtils::isAfterFSR

(

const xAOD::TruthParticle *

particle

)

Determine whether particle is afterFSR.

Definition at line 84 of file PartonHistoryUtils.cxx.

                                            {
  bool isAfter = true;
  for (size_t j = 0; j < p->nChildren(); j++) {
    if (p->child(j) && p->child(j)->pdgId() == p->pdgId())
      isAfter = false;
  }
  return isAfter;
}

isBrokenTop()

bool CP::PartonHistoryUtils::isBrokenTop

(

const xAOD::TruthParticle *

particle

)

Looking for tops without children -> must be broken.

Definition at line 117 of file PartonHistoryUtils.cxx.

                                             {
  // check if particle is a top without children.
  return (p->pdgId() == 6 && p->nChildren() == 0);
}

isChildOf()

bool CP::PartonHistoryUtils::isChildOf	(	const xAOD::TruthParticle *	parent,
		const xAOD::TruthParticle *	potentialChild )

Definition at line 60 of file PartonHistoryUtils.cxx.

                                               {
  // Checks if child is child of parent.
  if (!(parent && child))
    return false;
  if (parent->uid() == child->uid()) {
    return true;
  }
  // Loop through all children of the parent
  for (size_t i = 0; i < parent->nChildren(); ++i) {
    // Recursively check if child is a child of the current child
    if (isChildOf(parent->child(i), child)) {
      return true;
    }
  }
  // If child is not found in the children or their descendants, return
  // false
  return false;
}

isQuarkFromPDF()

bool CP::PartonHistoryUtils::isQuarkFromPDF

(

const xAOD::TruthParticle *

particle

)

Definition at line 129 of file PartonHistoryUtils.cxx.

                                                       {
  // In principle we could use some status codes here, e.g. 31/41/42/53/61 for
  // Pythia 8. But that is not guaranteed to be compatible across generators,
  // so instead we just check if it's a massless quark.
  bool isQuark = 1 <= particle->absPdgId() && particle->absPdgId() <= 5;
  bool isMassless = particle->m() == 0;
  return isQuark && isMassless;
}