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 BlockElementLinkerBase &)=delete

	BlockElementLinkerBase (const edm::ParameterSet &conf)

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::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

◆ linkPrefilter()

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

overridevirtual

Reimplemented from BlockElementLinkerBase.

Definition at line 24 of file TrackAndECALLinker.cc.

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

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

◆ testLink()

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

overridevirtual

Implements BlockElementLinkerBase.

Definition at line 39 of file TrackAndECALLinker.cc.

                                                                                                           {
   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);
   const reco::PFCluster::REPPoint& tkreppos = tkAtECAL.positionREP();
  
   // Check if the linking has been done using the KDTree algo
   // Glowinski & Gouzevitch
   if (useKDTree_ && tkelem->isMultilinksValide()) {  //KDTree Algo
     const reco::PFMultilinksType& multilinks = tkelem->getMultilinks();
     const double ecalphi = ecalreppos.Phi();
     const double ecaleta = ecalreppos.Eta();
  
     // Check if the link Track/Ecal exist
     reco::PFMultilinksType::const_iterator mlit = multilinks.begin();
     for (; mlit != multilinks.end(); ++mlit)
       if ((mlit->first == ecalphi) && (mlit->second == ecaleta))
         break;
  
     // If the link exist, we fill dist and linktest.
     if (mlit != multilinks.end()) {
       dist = LinkByRecHit::computeDist(ecaleta, ecalphi, tkreppos.Eta(), tkreppos.Phi());
     }
  
   } 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;
 }

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_.

Member Data Documentation

◆ debug_

const bool TrackAndECALLinker::debug_

private

Definition at line 19 of file TrackAndECALLinker.cc.

Referenced by testLink().

◆ useKDTree_

const bool TrackAndECALLinker::useKDTree_