/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Reconstruction/EmcRec/EmcRec-01-02-57/src/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;
     }
   }

}

---