Inheritance diagram for PreshowerAndECALLinker:

Public Member Functions
bool	linkPrefilter (const reco::PFBlockElement , const reco::PFBlockElement ) const override

	PreshowerAndECALLinker (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)

const std::string &	name () const

BlockElementLinkerBase &	operator= (const BlockElementLinkerBase &)=delete

virtual	~BlockElementLinkerBase ()=default

Private Attributes
bool	debug_

bool	useKDTree_

Detailed Description

Definition at line 6 of file PreshowerAndECALLinker.cc.

Constructor & Destructor Documentation

◆ PreshowerAndECALLinker()

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

inline

Definition at line 8 of file PreshowerAndECALLinker.cc.

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

Member Function Documentation

◆ linkPrefilter()

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

overridevirtual

Reimplemented from BlockElementLinkerBase.

Definition at line 23 of file PreshowerAndECALLinker.cc.

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

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

◆ testLink()

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

overridevirtual

Implements BlockElementLinkerBase.

Definition at line 42 of file PreshowerAndECALLinker.cc.

                                                                                                               {
   const reco::PFBlockElementCluster *pselem(nullptr), *ecalelem(nullptr);
   double dist(-1.0);
   if (elem1->type() < elem2->type()) {
     pselem = static_cast<const reco::PFBlockElementCluster*>(elem1);
     ecalelem = static_cast<const reco::PFBlockElementCluster*>(elem2);
   } else {
     pselem = static_cast<const reco::PFBlockElementCluster*>(elem2);
     ecalelem = static_cast<const reco::PFBlockElementCluster*>(elem1);
   }
   const reco::PFClusterRef& psref = pselem->clusterRef();
   const reco::PFClusterRef& ecalref = ecalelem->clusterRef();
   if (psref.isNull() || ecalref.isNull()) {
     throw cms::Exception("BadClusterRefs") << "PFBlockElementCluster's refs are null!";
   }
   // Check if the linking has been done using the KDTree algo
   // Glowinski & Gouzevitch
   if (useKDTree_ && pselem->isMultilinksValide(ecalelem->type())) {  // KDTree algo
     const reco::PFMultilinksType& multilinks = pselem->getMultilinks(ecalelem->type());
     const reco::PFCluster::REPPoint& ecalreppos = ecalref->positionREP();
     const math::XYZPoint& ecalxyzpos = ecalref->position();
     const math::XYZPoint& psxyzpos = psref->position();
     const double ecalPhi = ecalreppos.Phi();
     const double ecalEta = ecalreppos.Eta();
  
     // Check if the link PS/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(
           ecalxyzpos.X() / 1000., ecalxyzpos.Y() / 1000., psxyzpos.X() / 1000., psxyzpos.Y() / 1000., false);
     }
   } else {  //Old algorithm
     dist = LinkByRecHit::testECALAndPSByRecHit(*ecalref, *psref, debug_);
   }
   return dist;
 }

References reco::PFBlockElementCluster::clusterRef(), LinkByRecHit::computeDist(), debug_, egammaTools::ecalEta(), egammaTools::ecalPhi(), Exception, edm::Ref< C, T, F >::isNull(), LinkByRecHit::testECALAndPSByRecHit(), reco::PFBlockElement::type(), and useKDTree_.

Member Data Documentation

◆ debug_

bool PreshowerAndECALLinker::debug_

private

Definition at line 18 of file PreshowerAndECALLinker.cc.

Referenced by testLink().

◆ useKDTree_

bool PreshowerAndECALLinker::useKDTree_