Inheritance diagram for TrackAndHCALLinker:

Public Member Functions
double	testLink (const reco::PFBlockElement , const reco::PFBlockElement ) const override

	TrackAndHCALLinker (const edm::ParameterSet &conf)

Public Member Functions inherited from BlockElementLinkerBase
	BlockElementLinkerBase (const edm::ParameterSet &conf)

	BlockElementLinkerBase (const BlockElementLinkerBase &)=delete

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

const std::string &	name () const

BlockElementLinkerBase &	operator= (const BlockElementLinkerBase &)=delete

Private Attributes
bool	_debug

bool	_useKDTree

Detailed Description

Definition at line 9 of file TrackAndHCALLinker.cc.

Constructor & Destructor Documentation

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

inline

Definition at line 11 of file TrackAndHCALLinker.cc.

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

Member Function Documentation

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

overridevirtual

Implements BlockElementLinkerBase.

Definition at line 29 of file TrackAndHCALLinker.cc.

References reco::PFBlockElementCluster::clusterRef(), LinkByRecHit::computeDist(), constexpr, reco::deltaPhi(), reco::PFBlockElement::getMultilinks(), reco::PFTrajectoryPoint::HCALEntrance, reco::PFTrajectoryPoint::HCALExit, reco::PFBlockElement::isMultilinksValide(), reco::PFTrajectoryPoint::isValid(), NULL, reco::PFTrajectoryPoint::position(), reco::PFTrajectoryPoint::positionREP(), LinkByRecHit::testTrackAndClusterByRecHit(), reco::PFBlockElementTrack::trackRefPF(), and reco::PFBlockElement::type().

                                            {  
   constexpr reco::PFTrajectoryPoint::LayerType HCALEntrance =
     reco::PFTrajectoryPoint::HCALEntrance;
   constexpr reco::PFTrajectoryPoint::LayerType HCALExit =
     reco::PFTrajectoryPoint::HCALExit;
   const reco::PFBlockElementCluster *hcalelem(NULL);
   const reco::PFBlockElementTrack   *tkelem(NULL);
   double dist(-1.0);
   if( elem1->type() < elem2->type() ) {
     tkelem = static_cast<const reco::PFBlockElementTrack*>(elem1);
     hcalelem = static_cast<const reco::PFBlockElementCluster*>(elem2);
   } else {
     tkelem = static_cast<const reco::PFBlockElementTrack*>(elem2);
     hcalelem = static_cast<const reco::PFBlockElementCluster*>(elem1);
   }
   const reco::PFRecTrackRef& trackref = tkelem->trackRefPF();
   const reco::PFClusterRef& clusterref = hcalelem->clusterRef();
   const reco::PFCluster::REPPoint& hcalreppos = clusterref->positionREP(); 
   const reco::PFTrajectoryPoint& tkAtHCALEnt =
     trackref->extrapolatedPoint( HCALEntrance );
   const reco::PFTrajectoryPoint& tkAtHCALEx =
     trackref->extrapolatedPoint( HCALExit );
   const double dHEta = ( tkAtHCALEx.positionREP().Eta() - 
              tkAtHCALEnt.positionREP().Eta()  );
   const double dHPhi = reco::deltaPhi( tkAtHCALEx.positionREP().Phi(), 
                        tkAtHCALEnt.positionREP().Phi() );
   if ( _useKDTree && hcalelem->isMultilinksValide() ) { //KDTree Algo
     const reco::PFMultilinksType& multilinks = hcalelem->getMultilinks();
     const double tracketa = tkAtHCALEnt.positionREP().Eta();
     const double trackphi = tkAtHCALEnt.positionREP().Phi();
     
     // Check if the link Track/Hcal exist
     reco::PFMultilinksType::const_iterator mlit = multilinks.begin();
     for (; mlit != multilinks.end(); ++mlit)
       if ((mlit->first == trackphi) && (mlit->second == tracketa))
     break;
     
     // If the link exist, we fill dist and linktest.     
 
 
 
     if (mlit != multilinks.end()){     
 
 
       //special case ! A looper  can exit the barrel inwards and hit the endcap
       //In this case calculate the distance based on the first crossing since
       //the looper will probably never make it to the endcap
       if (tkAtHCALEx.position().R()<tkAtHCALEnt.position().R()) {
     dist = LinkByRecHit::computeDist(hcalreppos.Eta(), 
                      hcalreppos.Phi(), 
                      tracketa, 
                      trackphi);
     
     edm::LogWarning("TrackHCALLinker ") <<"Special case of linking with track hitting HCAL and looping back in the tracker ";
       }
       else {
     dist = LinkByRecHit::computeDist(hcalreppos.Eta(), 
                        hcalreppos.Phi(), 
                        tracketa + 0.1 * dHEta, 
                        trackphi + 0.1 * dHPhi);
       }
 
     }
 
   }    else {// Old algorithm
     if ( tkAtHCALEnt.isValid() )
       dist = LinkByRecHit::testTrackAndClusterByRecHit( *trackref, 
                             *clusterref,
                             false, _debug );
   }
   return dist;
 }

Member Data Documentation

bool TrackAndHCALLinker::_debug

private

Definition at line 21 of file TrackAndHCALLinker.cc.

Referenced by cmsPerfSuite.PerfSuite::optionParse(), and cmsPerfSuite.PerfSuite::runCmsReport().

bool TrackAndHCALLinker::_useKDTree