beamtime/sts-lasertest/go4/STS/TSTSClusterizer.cxx (r4864/r3054)

Go to the documentation of this file.

#include "TSTSClusterizer.h"
#include "TMath.h"

TSTSClusterizer::TSTSClusterizer( Int_t det, Int_t side ) :
        fDetector( det ),
        fSide( side )
{
        fTimeTolerance = kDefTimeTolerance;
        fSpaceTolerance = kDefSpaceTolerance;
        fTimeWindow = kDefTimeWindow;
}

TSTSClusterizer::~TSTSClusterizer()
{

}


Bool_t TSTSClusterizer::Match( TSTSDigi const digi1, const TSTSDigi digi2 )
{
        return ( TMath::Abs( digi1.GetTime() - digi2.GetTime() ) <= fTimeTolerance )
                        && ( TMath::Abs( digi1.GetStrip() - digi2.GetStrip() ) <= fSpaceTolerance );
}



Bool_t TSTSClusterizer::Match( TSTSCluster const * const cluster, const TSTSDigi digi )
{
        for( Int_t iDigi = 0; iDigi < cluster->GetNDigis(); iDigi++ ) {
                if( Match( digi, cluster->GetDigi( iDigi ) ) ) {
                        return true;
                }
        }
        return false;
}


/*
 * Expects that: cluster exists, itDigi valid
 */
void TSTSClusterizer::MoveDigiToCluster( TSTSCluster* cluster, const std::deque<TSTSDigi>::iterator& itDigi )
{
   cluster->AddDigi(*itDigi);
   fDigis.erase(itDigi);
}


/*
 * Expects that: cluster has at least one digi, itDigi exists
 */
Bool_t TSTSClusterizer::FindDigiFitsCluster( Bool_t noTrigger, TSTSCluster const * const cluster, std::deque<TSTSDigi>::iterator& itDigi )
{
        for( std::deque<TSTSDigi>::iterator it = fDigis.begin(); it != fDigis.end(); it++ ) {
                if( Match( cluster, *it ) ) {
                        if( noTrigger ) {
                                if( (fDigis.back().GetTime() - it->GetTime()) < fTimeTolerance ) {
                                        printf("Error in TClusterizer::FindDigiFitsCluster: Too small time window. Result might be wrong!\n");
                                }
                        }
                        itDigi = it;
                        return true;
                }
        }
        return false;
}



TSTSCluster* TSTSClusterizer::FetchCluster( Bool_t noTrigger )
{
        if( fDigis.size() == 0 ) {
                return 0;
        }
        if( noTrigger ) {
                if( fDigis.back().GetTime() - fDigis.front().GetTime() < fTimeWindow ) {
                        return 0;
                }
        }
        TSTSCluster* cluster = new TSTSCluster();
        MoveDigiToCluster(cluster, fDigis.begin());
        std::deque<TSTSDigi>::iterator itDigi;
        while ( FindDigiFitsCluster( noTrigger, cluster, itDigi ) ) {
                MoveDigiToCluster( cluster, itDigi );
        }
        return cluster;
}