EmcRecCluster2Shower.cxx

Go to the documentation of this file.

//
//  Convert from Cluster Map to Shower Map
//
//  Created by Zhe Wang 2004, 3, 29
//
#include <iostream>
 
#include "EmcRec/EmcRecCluster2Shower.h"
//#include "EmcRec/EmcRecSeedLocalMax.h"
#include "EmcRec/EmcRecSeedEThreshold.h"
#include "EmcRec/EmcRecSplitWeighted.h"
 
// Constructors and destructors
EmcRecCluster2Shower::EmcRecCluster2Shower()
{   
  //cout<<"====== EmcRec: Offline Mode ======"<<endl;
  //   fSeedFinder=new EmcRecSeedLocalMax;
  fSeedFinder=new EmcRecSeedEThreshold;
  fSplitter=new EmcRecSplitWeighted;
}
 
EmcRecCluster2Shower:: ~EmcRecCluster2Shower()
{
  delete fSeedFinder;
  delete fSplitter;
}
 
void EmcRecCluster2Shower::Convert(RecEmcClusterMap& aClusterMap,
    RecEmcShowerMap& aShowerMap)
{
  RecEmcClusterMap::iterator ciClusterMap;
 
  RecEmcIDVector aMaxVec;
  for(ciClusterMap=aClusterMap.begin();
      ciClusterMap!=aClusterMap.end();
      ++ciClusterMap)
  {
    //++++++++++++++++++++++++++
    //get its local maximum list
    fSeedFinder->Seed(ciClusterMap->second,aMaxVec);
    //++++++++++++++++++++++++++++++++++++++++++++++
    //put seeds to cluster
    if(!aMaxVec.empty()) {
      ci_RecEmcIDVector ciMax;
      for(ciMax=aMaxVec.begin();
          ciMax!=aMaxVec.end();
          ++ciMax) {
        ciClusterMap->second.InsertSeed(ciClusterMap->second.Find(*ciMax)->second);
      }
    }
    //++++++++++++++++++++++++++++++++++++++++++++++
    //split it into showers and push into shower map
    fSplitter->Split(ciClusterMap->second,aMaxVec,aShowerMap);
  }
}