TrackAndHCALLinker Class Reference

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
virtual	~BlockElementLinkerBase ()=default

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.

References testLink().

                                                  :
    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 _debug, _useKDTree, 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().

Referenced by TrackAndHCALLinker().

                                           {  
  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(), cmsPerfSuite.PerfSuite::runCmsReport(), and testLink().

bool TrackAndHCALLinker::_useKDTree

private

Definition at line 21 of file TrackAndHCALLinker.cc.

Referenced by testLink().