TracksterLinkingbySkeletons.cc File Reference

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <Math/VectorUtil.h>
#include "DataFormats/Math/interface/deltaR.h"
#include "DataFormats/HGCalReco/interface/Common.h"
#include "DataFormats/GeometrySurface/interface/BoundDisk.h"
#include "DataFormats/HGCalReco/interface/Trackster.h"
#include "TrackingTools/TrajectoryState/interface/TrajectoryStateTransform.h"
#include "RecoParticleFlow/PFProducer/interface/PFMuonAlgo.h"
#include "RecoHGCal/TICL/interface/TracksterLinkingAlgoBase.h"
#include "RecoHGCal/TICL/plugins/TracksterLinkingbySkeletons.h"
#include "TICLGraph.h"

Go to the source code of this file.

Functions
bool	isInCylinder (const std::array< ticl::Vector, 3 > &mySkeleton, const std::array< ticl::Vector, 3 > &otherSkeleton, const float radius_sqr)

Function Documentation

isInCylinder()

bool isInCylinder	(	const std::array< ticl::Vector, 3 > &	mySkeleton,
		const std::array< ticl::Vector, 3 > &	otherSkeleton,
		const float	radius_sqr
	)

Definition at line 172 of file TracksterLinkingbySkeletons.cc.

References HLT_2024v14_cff::distance, dqmdumpme::first, dqmdumpme::last, and amptDefault_cfi::proj.

Referenced by ticl::TracksterLinkingbySkeletons::areCompatible().

                                          {
  const auto &first = mySkeleton[0];
  const auto &last = mySkeleton[2];
  const auto &pointToCheck = otherSkeleton[0];

  const auto &cylAxis = last - first;
  const auto &vecToPoint = pointToCheck - first;

  auto axisNorm = cylAxis.Dot(cylAxis);
  auto projLength = vecToPoint.Dot(cylAxis) / axisNorm;
  bool isWithinLength = projLength >= 0 && projLength <= 1;

  if (!isWithinLength)
    return false;

  const auto &proj = cylAxis * projLength;

  const auto &pointOnAxis = first + proj;

  const auto &distance = pointToCheck - pointOnAxis;
  auto distance2 = distance.Dot(distance);

  bool isWithinRadius = distance2 <= radius_sqr;

  return isWithinRadius;
}