PHG4CylinderCellv1.cc

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file PHG4CylinderCellv1.cc

#include "PHG4CylinderCellv1.h"

#include <g4main/PHG4HitDefs.h>  // for keytype

using namespace std;

PHG4CylinderCellv1::PHG4CylinderCellv1():
  layer(0xFFFFFFFF),
  cellid(0xFFFFFFFF),
  binz(-1),
  binphi(-1),
  edeps(), showeredeps(), light_yield(0)
{}

void
PHG4CylinderCellv1::add_edep(const PHG4HitDefs::keytype g4hitid, const float edep)
{
  if (edeps.find(g4hitid) == edeps.end())
    {
      edeps[g4hitid] = edep;
    }
  else
    {
      edeps[g4hitid]+= edep;
    }
}

void
PHG4CylinderCellv1::add_edep(const PHG4HitDefs::keytype g4hitid, const float edep, const float ly)
{
  add_edep(g4hitid, edep);
  light_yield += ly;
}

void
PHG4CylinderCellv1::add_shower_edep(const int g4showerid, const float edep)
{
  if (showeredeps.find(g4showerid) == showeredeps.end())
    {
      showeredeps[g4showerid] = edep;
    }
  else
    {
      showeredeps[g4showerid]+= edep;
    }
}

double PHG4CylinderCellv1::get_edep() const {
  
  double esum = 0.0;
  EdepConstIterator iter;
  for (iter = edeps.begin(); iter != edeps.end(); ++iter) {
    esum += iter->second;
  }
  return esum;
}

void
PHG4CylinderCellv1::identify(std::ostream& os) const
{
  os << "PHG4CylinderCellv1: #" << cellid << " ";
  os << "(layer,binz,binphi,e) = (";
  os << layer << ","; 
  os << binz << ",";
  os << binphi << ",";
  os << get_edep();
  os << ")";
  os << endl;
}