mEmcGeaPatchSimulatedCluster.C

Go to the documentation of this file.

#include "mEmcGeaPatchSimulatedCluster.h"
#include "emcNodeHelper.h"
#include "emcClusterContainer.h"
#include "emcClusterContent.h"
#include <iostream>
#include <cassert>
#include "phool.h"

using std::cerr;
using std::cout;
using std::endl;

//_____________________________________________________________________________
mEmcGeaPatchSimulatedCluster::mEmcGeaPatchSimulatedCluster()
{
  name = "mEmcGeaPatchSimulatedCluster";
}

//_____________________________________________________________________________
mEmcGeaPatchSimulatedCluster::~mEmcGeaPatchSimulatedCluster()
{
}

//_____________________________________________________________________________
PHBoolean
mEmcGeaPatchSimulatedCluster::event(PHCompositeNode* topNode)
{
  emcNodeHelper nh;

  PHCompositeNode* dstNode = emcNodeHelper::findCompositeNode(topNode,"DST");

  if (!dstNode)
    {
      cerr << PHWHERE << " Could not find DST node !" << endl;
      return False;
    }

  emcClusterContainer* clusters = 
    nh.getObject<emcClusterContainer>("emcClusterContainer",dstNode);

  if (!clusters)
    {
      cerr << PHWHERE << " Could not find emcClusterContainer !" << endl;
      return False;
    }
  
  for ( size_t i = 0; i < clusters->size(); ++i )
    {
      emcClusterContent* c = clusters->getCluster(i);
      assert(c!=0);
      // indicate that clusters are simulated, i.e. Fraction of simulated
      // energy = 1 (100%)
      c->set_simfrac(1);    
    }

  return True;
}