#include <DAQQuickClustering_module.h>

Public Member Functions
	ClusterHitsInTime (double TimeWindow)

void	DoIt (std::vector< recoHit > cHitVector)

double	GetTimeWindow () const

int	getNClusters () const

std::vector< cluster >	getClusterVector () const

void	SetTimeWindow (double inTimeWindow=1)

Private Attributes
int	fNClusters = 0

double	fTimeWindow = 0

std::vector< recoHit >	fHitVector

std::vector< cluster >	fVecClusters

Detailed Description

Definition at line 170 of file DAQQuickClustering_module.h.

Constructor & Destructor Documentation

ClusterHitsInTime::ClusterHitsInTime ( double TimeWindow )

inline

Definition at line 173 of file DAQQuickClustering_module.h.

                                       : fNClusters(0),
                                            fTimeWindow(TimeWindow),
                                            fHitVector(),
                                            fVecClusters(){};

Member Function Documentation

void ClusterHitsInTime::DoIt ( std::vector< recoHit > cHitVector )

Definition at line 3 of file DAQQuickClustering_module.cc.

 {
   fHitVector.clear();
   fHitVector = cHitVector;
 
   //ORDER IN THE HIT VECTOR BY TIME.
   std::sort(fHitVector.begin(), fHitVector.end(),
             [](const recoHit& lhs, const recoHit& rhs){return lhs<rhs;});
 
   for(unsigned int i = 0; i < fHitVector.size()-1; i++)
   {
     //std::cout << "fHitVector.at(" << i << ").getHitTime() " << fHitVector.at(i).getHitTime()
     //         << std::endl;
     std::vector<recoHit> vec_TempHits;
     if(std::abs(fHitVector.at(i).getHitTime()-fHitVector.at(i+1).getHitTime())<=fTimeWindow)
     {
       int  timeCount = 1;
 
       vec_TempHits.push_back(fHitVector.at(i));
       vec_TempHits.push_back(fHitVector.at(i+1));
 
       while((i+timeCount+1)<fHitVector.size() &&
             std::abs(fHitVector.at(i+timeCount).getHitTime()-fHitVector.at(i+timeCount+1).getHitTime()) <= fTimeWindow)
       {
         vec_TempHits.push_back(fHitVector.at(i + timeCount + 1));
         timeCount++;
       }
 
       i = i + timeCount;
       cluster temp(fHitVector.at(0).getEvent(), vec_TempHits);
       fVecClusters.push_back(temp);
     }
   }
 
   fNClusters = fVecClusters.size();
 }