TrackAndECALLinker Class Reference

Inheritance diagram for TrackAndECALLinker:

Public Member Functions
bool	linkPrefilter (const reco::PFBlockElement *, const reco::PFBlockElement *) const override
double	testLink (const reco::PFBlockElement *, const reco::PFBlockElement *) const override
	TrackAndECALLinker (const edm::ParameterSet &conf)
Public Member Functions inherited from BlockElementLinkerBase
	BlockElementLinkerBase (const edm::ParameterSet &conf)
	BlockElementLinkerBase (const BlockElementLinkerBase &)=delete
const std::string &	name () const
BlockElementLinkerBase &	operator= (const BlockElementLinkerBase &)=delete
virtual	~BlockElementLinkerBase ()=default

Private Attributes
const bool	debug_
const bool	useKDTree_

Detailed Description

Definition at line 7 of file TrackAndECALLinker.cc.

Constructor & Destructor Documentation

TrackAndECALLinker::TrackAndECALLinker ( const edm::ParameterSet &  conf )

inline

Definition at line 9 of file TrackAndECALLinker.cc.

      : BlockElementLinkerBase(conf),
        useKDTree_(conf.getParameter<bool>("useKDTree")),
        debug_(conf.getUntrackedParameter<bool>("debug", false)) {}

Member Function Documentation

bool TrackAndECALLinker::linkPrefilter ( const reco::PFBlockElement *  elem1, const reco::PFBlockElement *  elem2  ) const

overridevirtual

Reimplemented from BlockElementLinkerBase.

Definition at line 24 of file TrackAndECALLinker.cc.

References reco::PFBlockElement::ECAL, relativeConstraints::empty, reco::PFBlockElement::getMultilinks(), reco::PFBlockElement::isMultilinksValide(), mps_fire::result, reco::PFBlockElement::TRACK, reco::PFBlockElement::type(), and useKDTree_.

                                                                                                             {
  bool result = false;
  // Track-ECAL KDTree multilinks are stored to ecal's elem
  switch (elem1->type()) {
    case reco::PFBlockElement::TRACK:
      result = (elem2->isMultilinksValide(elem1->type()) && !elem2->getMultilinks(elem1->type()).empty() &&
                elem1->isMultilinksValide(elem2->type()));
      break;
    case reco::PFBlockElement::ECAL:
      result = (elem1->isMultilinksValide(elem2->type()) && !elem1->getMultilinks(elem2->type()).empty() &&
                elem2->isMultilinksValide(elem1->type()));
    default:
      break;
  }
  return (useKDTree_ ? result : true);
}

double TrackAndECALLinker::testLink ( const reco::PFBlockElement *  elem1, const reco::PFBlockElement *  elem2  ) const

overridevirtual

Implements BlockElementLinkerBase.

Definition at line 41 of file TrackAndECALLinker.cc.

References reco::PFBlockElementCluster::clusterRef(), LinkByRecHit::computeDist(), gather_cfg::cout, debug_, reco::PFTrajectoryPoint::ECALShowerMax, reco::PFBlockElement::getMultilinks(), reco::PFBlockElement::isMultilinksValide(), reco::PFTrajectoryPoint::isValid(), reco::PFTrajectoryPoint::positionREP(), LinkByRecHit::testTrackAndClusterByRecHit(), reco::PFBlockElementTrack::trackRefPF(), reco::PFBlockElement::type(), and useKDTree_.

                                                                                                          {
  constexpr reco::PFTrajectoryPoint::LayerType ECALShowerMax = reco::PFTrajectoryPoint::ECALShowerMax;
  const reco::PFBlockElementCluster* ecalelem(nullptr);
  const reco::PFBlockElementTrack* tkelem(nullptr);
  double dist(-1.0);
  if (elem1->type() < elem2->type()) {
    tkelem = static_cast<const reco::PFBlockElementTrack*>(elem1);
    ecalelem = static_cast<const reco::PFBlockElementCluster*>(elem2);
  } else {
    tkelem = static_cast<const reco::PFBlockElementTrack*>(elem2);
    ecalelem = static_cast<const reco::PFBlockElementCluster*>(elem1);
  }
  const reco::PFRecTrackRef& trackref = tkelem->trackRefPF();
  const reco::PFClusterRef& clusterref = ecalelem->clusterRef();
  const reco::PFCluster::REPPoint& ecalreppos = clusterref->positionREP();
  const reco::PFTrajectoryPoint& tkAtECAL = trackref->extrapolatedPoint(ECALShowerMax);

  // Check if the linking has been done using the KDTree algo
  // Glowinski & Gouzevitch
  if (useKDTree_ && ecalelem->isMultilinksValide(tkelem->type())) {  //KDTree Algo
    const reco::PFMultilinksType& multilinks = ecalelem->getMultilinks(tkelem->type());
    const double tracketa = tkAtECAL.positionREP().Eta();
    const double trackphi = tkAtECAL.positionREP().Phi();
    // Check if the link Track/Ecal exist
    reco::PFMultilinksType::const_iterator mlit = multilinks.begin();
    for (; mlit != multilinks.end(); ++mlit)
      if (mlit->trackRef == trackref)
        break;

    // If the link exist, we fill dist and linktest.
    if (mlit != multilinks.end()) {
      dist = LinkByRecHit::computeDist(ecalreppos.Eta(), ecalreppos.Phi(), tracketa, trackphi);
    }

  } else {  // Old algorithm
    if (tkAtECAL.isValid())
      dist = LinkByRecHit::testTrackAndClusterByRecHit(*trackref, *clusterref, false, debug_);
  }

  if (debug_) {
    if (dist > 0.) {
      std::cout << " Here a link has been established"
                << " between a track an Ecal with dist  " << dist << std::endl;
    } else
      std::cout << " No link found " << std::endl;
  }
  return dist;
}

Member Data Documentation

const bool TrackAndECALLinker::debug_

private

Definition at line 19 of file TrackAndECALLinker.cc.

Referenced by testLink().

const bool TrackAndECALLinker::useKDTree_

private

Definition at line 19 of file TrackAndECALLinker.cc.

Referenced by linkPrefilter(), and testLink().