DumpRawTowerContainer.cc

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

#include "DumpRawTowerContainer.h"

#include <phool/PHIODataNode.h>

#include <calobase/RawTower.h>
#include <calobase/RawTowerContainer.h>

#include <map>
#include <ostream>
#include <string>
#include <utility>

using MyNode_t = PHIODataNode<RawTowerContainer>;

DumpRawTowerContainer::DumpRawTowerContainer(const std::string &NodeName)
  : DumpObject(NodeName)
{
  return;
}

int DumpRawTowerContainer::process_Node(PHNode *myNode)
{
  RawTowerContainer *rawtowercontainer = nullptr;
  MyNode_t *thisNode = static_cast<MyNode_t *>(myNode);
  if (thisNode)
  {
    rawtowercontainer = thisNode->getData();
  }
  if (rawtowercontainer)
  {
    RawTowerContainer::ConstIterator hiter;
    RawTowerContainer::ConstRange begin_end = rawtowercontainer->getTowers();
    *fout << "size: " << rawtowercontainer->size() << std::endl;
    for (hiter = begin_end.first; hiter != begin_end.second; ++hiter)
    {
      RawTower *rawtwr = hiter->second;
      *fout << "bineta: " << rawtwr->get_bineta() << std::endl;
      *fout << "binphi: " << rawtwr->get_binphi() << std::endl;
      *fout << "energy: " << rawtwr->get_energy() << std::endl;
      RawTower::CellConstRange cbegin_end = rawtwr->get_g4cells();
      for (RawTower::CellConstIterator iter = cbegin_end.first; iter != cbegin_end.second; ++iter)
      {
        *fout << "cell key: 0x" << std::hex << iter->first
              << std::dec << " edep: " << iter->second << std::endl;
      }
      RawTower::ShowerConstRange sbegin_end = rawtwr->get_g4showers();
      for (RawTower::ShowerConstIterator iter = sbegin_end.first; iter != sbegin_end.second; ++iter)
      {
        *fout << "shower id: " << iter->first
              << " edep: " << iter->second << std::endl;
      }
    }
  }
  return 0;
}