/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Reconstruction/MdcPatRec/MdcTrkRecon/MdcTrkRecon-00-03-45/src/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 "MdcTrkRecon/MdcSegUsage.h"
#include "MdcGeom/MdcDetector.h"
#include "MdcData/MdcHitMap.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;
}

---