MdcSegData.cxx

Go to the documentation of this file.

//--------------------------------------------------------------------------
// File and Version Information:
//      $Id: MdcSegData.cxx,v 1.5 2011/09/26 01:06:37 zhangy Exp $
//
// Author(s): Steve Schaffner
//------------------------------------------------------------------------
 
#include "MdcTrkRecon/MdcSegData.h"
#include "MdcData/MdcHit.h"
#include "MdcData/MdcHitMap.h"
#include "MdcGeom/MdcDetector.h"
#include "MdcTrkRecon/MdcSegUsage.h"
#include "MdcTrkRecon/mdcWrapWire.h"
 
//------------------------------------------------------------------------
MdcSegData::MdcSegData( bool skip )
//------------------------------------------------------------------------
{
  eventNumber = 0;
  _segUsage   = 0;
  _hitMap     = 0;
  _bunchTime  = 0;
  _skipUsed   = skip;
 
  _theHits = 0;
}
 
//------------------------------------------------------------------------
MdcSegData::~MdcSegData( void ) {
  //------------------------------------------------------------------------
  delete[] _segUsage;
}
 
//------------------------------------------------------------------------
int MdcSegData::nhits() const {
  //------------------------------------------------------------------------
  return _theHits->size();
}
 
//------------------------------------------------------------------------
void MdcSegData::loadevent( MdcHitCol* hitcol, MdcHitMap* hitmap, double tb ) {
  //------------------------------------------------------------------------
  _bunchTime = tb;
  _theHits   = hitcol;
  _hitMap    = hitmap;
  assert( _hitMap != 0 );
 
  _segUsageDict.clear();
  delete[] _segUsage;
  _segUsage = 0;
  eventNumber++;
 
  _segUsage = new MdcSegUsage[nhits()];
 
  MdcHit*             aHit  = 0;
  MdcHitCol::iterator iter  = _theHits->begin();
  int                 index = 0;
  for ( ; iter != _theHits->end(); iter++ )
  {
    aHit = *iter;
    _segUsageDict.put( aHit, &( _segUsage[index] ) );
    if ( skippingUsed() && aHit->usedHit() )
    {
      // Mark for later skipping
      _segUsage[index].killHit();
    }
    index++;
  }
}
 
//------------------------------------------------------------------------
const MdcHit* MdcSegData::hit( int hitno ) const {
  //------------------------------------------------------------------------
  return ( *_theHits )[hitno];
}
 
//------------------------------------------------------------------------
void MdcSegData::poisonHits( const MdcDetector* gm, int debug ) {
  //------------------------------------------------------------------------
  // Mark hits as poor candidates for segment-finding, if they are located
  //  within a string of consecutive hits in the same layer (end hits are
  //  not marked).
  for ( int ilayer = 0; ilayer < gm->nLayer(); ilayer++ )
  {
    int nwire = gm->Layer( ilayer )->nWires();
    for ( int iwire = 0; iwire < nwire; ++iwire )
    {
      int wireNext = mdcWrapWire( iwire + 1, nwire );
      if ( 0 == hitMap()->hitWire( ilayer, wireNext ) )
      {
        //                        0 1 2 3 4
        // next wire not there:   ? o ? ? ?
        // so the next-to-next wire shouldn't be poisened
        // (it's left neighbour is empty)
        // and the next candidate to check is + 3
        iwire += 2;
        continue;
      }
 
      MdcHit* theHit = hitMap()->hitWire( ilayer, iwire );
      if ( theHit == 0 )
      {
        //                                        0 1 2 3 4
        // next wire is there, this one isn't :   o x ? ? ?
        // and the next candidate to check is + 2
        ++iwire;
        continue;
      }
 
      int wirePrev = mdcWrapWire( iwire - 1, nwire );
      if ( 0 != hitMap()->hitWire( ilayer, wirePrev ) )
      {
        segUsage().get( theHit )->killHit();
        if ( debug > 1 )
        {
          theHit->print( std::cout );
          std::cout << " killed " << std::endl;
        }
      }
      else
      {
        if ( debug > 1 )
        {
          theHit->print( std::cout );
          std::cout << " ok " << std::endl;
        }
      }
    } // end wire loop
  }   // end layer loop
 
  return;
}