Inheritance diagram for GSFAndHGCalLinker:

Public Member Functions
	GSFAndHGCalLinker (const edm::ParameterSet &conf)

double	testLink (const reco::PFBlockElement , const reco::PFBlockElement ) const override

Public Member Functions inherited from BlockElementLinkerBase
	BlockElementLinkerBase (const BlockElementLinkerBase &)=delete

	BlockElementLinkerBase (const edm::ParameterSet &conf)

virtual bool	linkPrefilter (const reco::PFBlockElement , const reco::PFBlockElement ) const

const std::string &	name () const

BlockElementLinkerBase &	operator= (const BlockElementLinkerBase &)=delete

virtual	~BlockElementLinkerBase ()=default

Private Attributes
bool	debug_

bool	useKDTree_

Detailed Description

Definition at line 7 of file GSFAndHGCalLinker.cc.

Constructor & Destructor Documentation

◆ GSFAndHGCalLinker()

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

inline

Definition at line 9 of file GSFAndHGCalLinker.cc.

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

Member Function Documentation

◆ testLink()

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

overridevirtual

Implements BlockElementLinkerBase.

Definition at line 22 of file GSFAndHGCalLinker.cc.

                                                                                                          {
   constexpr reco::PFTrajectoryPoint::LayerType ECALShowerMax = reco::PFTrajectoryPoint::ECALShowerMax;
   const reco::PFBlockElementCluster* hgcalelem(nullptr);
   const reco::PFBlockElementGsfTrack* gsfelem(nullptr);
   double dist(-1.0);
   if (elem1->type() > elem2->type()) {
     hgcalelem = static_cast<const reco::PFBlockElementCluster*>(elem1);
     gsfelem = static_cast<const reco::PFBlockElementGsfTrack*>(elem2);
   } else {
     hgcalelem = static_cast<const reco::PFBlockElementCluster*>(elem2);
     gsfelem = static_cast<const reco::PFBlockElementGsfTrack*>(elem1);
   }
   const reco::PFRecTrack& track = gsfelem->GsftrackPF();
   const reco::PFClusterRef& clusterref = hgcalelem->clusterRef();
   const reco::PFTrajectoryPoint& tkAtECAL = track.extrapolatedPoint(ECALShowerMax);
   if (tkAtECAL.isValid()) {
     dist = LinkByRecHit::computeDist(tkAtECAL.positionREP().eta(),
                                      tkAtECAL.positionREP().phi(),
                                      clusterref->positionREP().Eta(),
                                      clusterref->positionREP().Phi());
     if (dist > 0.3)
       dist = -1.0;
   }
   if (debug_) {
     if (dist > 0.) {
       std::cout << " Here a link has been established"
                 << " between a GSF track an HGCal with dist  " << dist << std::endl;
     } else {
       if (debug_)
         std::cout << " No link found " << std::endl;
     }
   }
  
   return dist;
 }

References reco::PFBlockElementCluster::clusterRef(), LinkByRecHit::computeDist(), gather_cfg::cout, debug_, reco::PFTrajectoryPoint::ECALShowerMax, reco::PFBlockElementGsfTrack::GsftrackPF(), reco::PFTrajectoryPoint::isValid(), reco::PFTrajectoryPoint::positionREP(), HLT_2018_cff::track, and reco::PFBlockElement::type().

Member Data Documentation

◆ debug_

bool GSFAndHGCalLinker::debug_

private

Definition at line 17 of file GSFAndHGCalLinker.cc.

Referenced by testLink().

◆ useKDTree_

bool GSFAndHGCalLinker::useKDTree_