EmcRecHit2Cluster Class Reference

#include <EmcRecHit2Cluster.h>

List of all members.

Public Member Functions
	EmcRecHit2Cluster ()
	~EmcRecHit2Cluster ()
void	Convert (const RecEmcHitMap &aHitMap, RecEmcClusterMap &aClusterMap)

---

Detailed Description

Definition at line 12 of file EmcRecHit2Cluster.h.

---

Constructor & Destructor Documentation

EmcRecHit2Cluster::EmcRecHit2Cluster

(

)

Definition at line 12 of file EmcRecHit2Cluster.cxx.

{}

EmcRecHit2Cluster::~EmcRecHit2Cluster

(

)

Definition at line 15 of file EmcRecHit2Cluster.cxx.

{}

---

Member Function Documentation

void EmcRecHit2Cluster::Convert	(	const RecEmcHitMap &	aHitMap,
		RecEmcClusterMap &	aClusterMap
	)

Definition at line 18 of file EmcRecHit2Cluster.cxx.

References RecEmcCluster::Clear(), RecEmcCluster::ClusterId(), EmcRecParameter::EThresholdCluster(), EmcRecParameter::GetInstance(), EmcRecNeighbor::GetNeighbors(), and RecEmcCluster::Insert().

Referenced by EmcRec::execute().

{
   RecEmcHitMap tHitMap,aBufferHitMap;
   RecEmcHitMap::iterator pHitMap,pBufferHitMap;
   tHitMap=aHitMap;
   
   RecEmcID CellId;
   RecEmcIDVector NearCell;
   i_RecEmcIDVector pNearCell;
   
   RecEmcCluster aCluster;

   EmcRecParameter& Para=EmcRecParameter::GetInstance();
   EmcRecNeighbor nhb;
   
   while(!tHitMap.empty())
   {
     double eCluster=0;   //total energy of a cluster
     double emax=0;   //seed's energy of a cluster
     RecEmcID idmax;    //seed's id of a cluster
     
           // Search a cluster.
           pHitMap=tHitMap.begin();
           aBufferHitMap[pHitMap->first]=pHitMap->second;
           tHitMap.erase(pHitMap);
           aCluster.Clear();
           while(!aBufferHitMap.empty())
           {
             pBufferHitMap=aBufferHitMap.begin();
             CellId=pBufferHitMap->first;
             NearCell=nhb.GetNeighbors(CellId);
             for(pNearCell=NearCell.begin();
                       pNearCell!=NearCell.end();
                       pNearCell++)
             {
                     pHitMap=tHitMap.find(*pNearCell);
                     if(pHitMap!=tHitMap.end())
                     {
                       aBufferHitMap[pHitMap->first]=pHitMap->second;
                       tHitMap.erase(pHitMap);
                     }
             }

       // Find seed
       if(pBufferHitMap->second.getEnergy()>emax) {
         emax=pBufferHitMap->second.getEnergy();
         idmax=pBufferHitMap->first;
       }
         
       eCluster+=pBufferHitMap->second.getEnergy();
             aCluster.Insert(pBufferHitMap->second);
             aBufferHitMap.erase(pBufferHitMap);
           }
         // Now, done and save it.
     if(eCluster>Para.EThresholdCluster())
     {
       aCluster.ClusterId(idmax);
             aClusterMap[idmax]=aCluster;
             //aClusterMap[aCluster.ClusterId()]=aCluster;
     }
   }

}

---