TrackAndECALLinker Class Reference

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 edm::ParameterSet &conf)
	BlockElementLinkerBase (const BlockElementLinkerBase &)=delete
const std::string &	name () const
BlockElementLinkerBase &	operator= (const BlockElementLinkerBase &)=delete

Private Attributes
const bool	_debug
const bool	_useKDTree

Detailed Description

Definition at line 7 of file TrackAndECALLinker.cc.

Constructor & Destructor Documentation

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

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

overridevirtual

Reimplemented from BlockElementLinkerBase.

Definition at line 29 of file TrackAndECALLinker.cc.

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

                                                   {  
  bool result = false;
  switch( elem1->type() ){ 
  case reco::PFBlockElement::TRACK:
    result = (elem1->isMultilinksValide() && elem1->getMultilinks().size() > 0);
    break;
  case reco::PFBlockElement::ECAL:
    result = (elem2->isMultilinksValide() && elem2->getMultilinks().size() > 0);
  default:
    break;
  } 
  return (_useKDTree ? result : true);  
}

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

overridevirtual

Implements BlockElementLinkerBase.

Definition at line 45 of file TrackAndECALLinker.cc.

References _debug, _useKDTree, reco::PFBlockElementCluster::clusterRef(), LinkByRecHit::computeDist(), constexpr, gather_cfg::cout, reco::PFTrajectoryPoint::ECALShowerMax, reco::PFBlockElement::getMultilinks(), reco::PFBlockElement::isMultilinksValide(), reco::PFTrajectoryPoint::isValid(), NULL, reco::PFTrajectoryPoint::positionREP(), LinkByRecHit::testTrackAndClusterByRecHit(), reco::PFBlockElementTrack::trackRefPF(), and reco::PFBlockElement::type().

                                              {  
  constexpr reco::PFTrajectoryPoint::LayerType ECALShowerMax =
    reco::PFTrajectoryPoint::ECALShowerMax;
  const reco::PFBlockElementCluster *ecalelem(NULL);
  const reco::PFBlockElementTrack   *tkelem(NULL);
  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;
}

Member Data Documentation

const bool TrackAndECALLinker::_debug

private

Definition at line 21 of file TrackAndECALLinker.cc.

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

const bool TrackAndECALLinker::_useKDTree

private

Definition at line 21 of file TrackAndECALLinker.cc.

Referenced by linkPrefilter(), and testLink().