#include <EmcRecHit2Cluster.h>

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::EmcRecHit2Cluster ( )

Definition at line 12 of file EmcRecHit2Cluster.cxx.

13{}

◆ ~EmcRecHit2Cluster()

EmcRecHit2Cluster::~EmcRecHit2Cluster ( )

Definition at line 15 of file EmcRecHit2Cluster.cxx.

16{}

Member Function Documentation

◆ Convert()

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

Definition at line 18 of file EmcRecHit2Cluster.cxx.

{
   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;
     }
   }
 
}

Referenced by EmcRec::execute().