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;
}