IntermediateTrajectoryCleaner Class Reference

#include <IntermediateTrajectoryCleaner.h>

List of all members.

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 cond::ecalcond::all, and TempTrajectory::isValid().

                                                                                                          {

  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::not1(std::mem_fun_ref(&TempTrajectory::isValid))),
 //                                     boost::bind(&TempTrajectory::isValid,_1)), 
                        theTrajectories.end());
}