beamtime/sts-lasertest/go4/STS/TSTSHiterizer.cxx (r4864/r3065)

Go to the documentation of this file.

/*
 * TSTSHiterizer.cxx
 *
 *  Created on: Jan 9, 2012
 *      Author: isorokin
 */

#include "TSTSHiterizer.h"

TSTSHiterizer::TSTSHiterizer( Int_t det )
{
        fDetector = det;
        fTimeTolerance = kDefTimeTolerance;
}

TSTSHiterizer::~TSTSHiterizer()
{
}


void TSTSHiterizer::ReconstructHits( Bool_t noTrigger )
{
   if( fClusters[0].size() < 2 ) return;
   if( fClusters[1].size() < 2 ) return;

   Int_t oldestClusterSide = ( fClusters[0].front()->GetTime() < fClusters[1].front()->GetTime() ) ? 0 : 1;
   TSTSCluster* oldestCluster = fClusters[oldestClusterSide].front();

   Int_t oppositeSide = ! oldestClusterSide;

   if( noTrigger ) {
           if( fClusters[oppositeSide].back()->GetTime() - oldestCluster->GetTime() <= 2 * fTimeTolerance ) return;
   }

   std::deque<TSTSCluster*>::iterator matchClustIter = fClusters[oppositeSide].begin();
   while( matchClustIter != fClusters[oppositeSide].end() )
   {
      if( ( * matchClustIter )->GetTime() - oldestCluster->GetTime() <= fTimeTolerance ) {
         TSTSCluster* matchedClusters[2];
         matchedClusters[oldestClusterSide] = new TSTSCluster(*oldestCluster);
         matchedClusters[oppositeSide] = new TSTSCluster( *(*matchClustIter) );
         fHits.push_back( new TSTSHit( matchedClusters[0], matchedClusters[1] ) );
      }
      matchClustIter++;
   }

   delete fClusters[oldestClusterSide].front();
   fClusters[oldestClusterSide].pop_front();
}


TSTSHit * TSTSHiterizer::FetchHit( Bool_t noTrigger )
{
   ReconstructHits( noTrigger );
   if( fHits.size() > 0 ) {
      TSTSHit* returnHit = new TSTSHit( *(fHits.front()) );
      fHits.pop_front();
      return returnHit;
   }
   return 0;
}

void TSTSHiterizer::Clear()
      {
         while(!fHits.empty())
              {
                  TSTSHit* hit=fHits.front();
                  if(hit) delete hit;
                  fHits.pop_front();
              }
         for(int side=0;side<TSTSTopology::kNSides;++side){
            while(!fClusters[side].empty())
               {
                  TSTSCluster* cluster=fClusters[side].front();
                  if(cluster) delete cluster;
                  fClusters[side].pop_front();
               }
         }

      }