IntermediateTrajectoryCleaner Class Reference

#include <IntermediateTrajectoryCleaner.h>

Static Public Member Functions
static void	clean (TempTrajectoryContainer &tracks)

Private Types
typedef TransientTrackingRecHit::ConstRecHitPointer	ConstRecHitPointer
typedef BaseCkfTrajectoryBuilder::TempTrajectoryContainer	TempTrajectoryContainer

Detailed Description

Definition at line 8 of file IntermediateTrajectoryCleaner.h.

Member Typedef Documentation

typedef TransientTrackingRecHit::ConstRecHitPointer IntermediateTrajectoryCleaner::ConstRecHitPointer

private

Definition at line 10 of file IntermediateTrajectoryCleaner.h.

typedef BaseCkfTrajectoryBuilder::TempTrajectoryContainer IntermediateTrajectoryCleaner::TempTrajectoryContainer

private

Definition at line 9 of file IntermediateTrajectoryCleaner.h.

Member Function Documentation

void IntermediateTrajectoryCleaner::clean ( IntermediateTrajectoryCleaner::TempTrajectoryContainer &  theTrajectories )

static

Definition at line 14 of file IntermediateTrajectoryCleaner.cc.

References TrackingRecHit::all, and TempTrajectory::isValid().

Referenced by GroupedCkfTrajectoryBuilder::groupedLimitedCandidates(), and CkfTrajectoryBuilder::limitedCandidates().

                                                                                                          {

  if (theTrajectories.empty()) return;
  if (theTrajectories[0].measurements().size()<4) return;

  for (TempTrajectoryContainer::iterator firstTraj=theTrajectories.begin(), firstEnd=theTrajectories.end() - 1;
     firstTraj != firstEnd; ++firstTraj) {

    if ( (!firstTraj->isValid()) ||
         (!firstTraj->lastMeasurement().recHit()->isValid()) ) continue;
    
    TempTrajectory::DataContainer::const_iterator itFirst = firstTraj->measurements().rbegin();
    ConstRecHitPointer first_hit1 = itFirst->recHit(); --itFirst; 
    ConstRecHitPointer first_hit2 = itFirst->recHit(); --itFirst; 
    ConstRecHitPointer first_hit3 = itFirst->recHit();

    bool fh2Valid = first_hit2->isValid();

    for (TempTrajectoryContainer::iterator secondTraj = (firstTraj+1), secondEnd = theTrajectories.end();
       secondTraj != secondEnd; ++secondTraj) {

      if ( (!secondTraj->isValid()) ||
           (!secondTraj->lastMeasurement().recHit()->isValid()) ) continue;

        TempTrajectory::DataContainer::const_iterator itSecond = secondTraj->measurements().rbegin();
        ConstRecHitPointer second_hit1 = itSecond->recHit(); --itSecond; 
        ConstRecHitPointer second_hit2 = itSecond->recHit(); --itSecond; 
        ConstRecHitPointer second_hit3 = itSecond->recHit();

    if ( &(*first_hit3) == &(*second_hit3) ) {
      if (fh2Valid ^ second_hit2->isValid()) { // ^ = XOR!
        if ( first_hit1->hit()->sharesInput( second_hit1->hit(), TrackingRecHit::all ) ){
          
          if (!fh2Valid) {
        firstTraj->invalidate();
        break;
          }
          // else // if (!second_hit2->isValid())  // should be always true, as we did XOR !
          secondTraj->invalidate();
        }
      }
    }
    }
  }
  theTrajectories.erase(std::remove_if( theTrajectories.begin(),theTrajectories.end(),
                                        std::not_fn(&TempTrajectory::isValid)),
                        theTrajectories.end());
}