EmcRecHit2Cluster.cxx

Go to the documentation of this file.

//
//  Convert from Hit Map to Cluster Map
//
//  Wang.Zhe 2004, 3, 29
//
#include <iostream>
 
#include "EmcRec/EmcRecHit2Cluster.h" 
#include "EmcRec/EmcRecParameter.h"
 
// Constructors and destructors
EmcRecHit2Cluster::EmcRecHit2Cluster()
{}
 
EmcRecHit2Cluster:: ~EmcRecHit2Cluster()
{}
 
void EmcRecHit2Cluster::Convert(const RecEmcHitMap& aHitMap,
                   RecEmcClusterMap& aClusterMap)
{
   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;
     }
   }
 
}