/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Reconstruction/EmcRec/EmcRec-01-02-57/src/EmcRecSeedLocalMax.cxx

Go to the documentation of this file.

//
//  Seed finder abstract base class
//  Use local max to determine a seed
//
//  Zhe Wang 2004, 8, 31
//
#include <iostream>

#include "EmcRec/EmcRecSeedLocalMax.h"

EmcRecSeedLocalMax::EmcRecSeedLocalMax()
{}

EmcRecSeedLocalMax::~EmcRecSeedLocalMax()
{}

void EmcRecSeedLocalMax::Seed(const RecEmcCluster& aCluster,
                              RecEmcIDVector& aMaxVec)
{
   RecEmcHitMap::const_iterator ciHitMap;
   RecEmcHitMap::const_iterator ciNeigh;
   RecEmcHit aHit;
   
   bool max;
   
   RecEmcIDVector NeighVec;
   RecEmcIDVector::const_iterator ciNeighVec;
      
   aMaxVec.clear();
   
   EmcRecNeighbor nhb;

   for(ciHitMap=aCluster.Begin();
       ciHitMap!=aCluster.End();
       ++ciHitMap) {
      //++++++++++++++++++++++++++++++++++++++++++++++++++
      //prepare a hit
      max=true;
      aHit=ciHitMap->second;
      //prepare its neighbours
      NeighVec.clear();
      NeighVec=nhb.GetNeighbors(aHit.getCellId());
      
      //to see if it is energy max
      for(ciNeighVec=NeighVec.begin();
          ciNeighVec!=NeighVec.end();
          ++ciNeighVec) {
         //find it
         ciNeigh=aCluster.Find(*ciNeighVec);
         if(ciNeigh!=aCluster.End()) {
            //largest!!!
            if(ciNeigh->second.getEnergy()>=aHit.getEnergy()) {
               max=false;
            }
         }
      }
         
      //once max hit is found
      if(max==true) {
        aMaxVec.push_back(aHit.getCellId());
      }
      //++++++++++++++++++++++++++++++++++++++++++++++++++
   }
}

---