LArVoxelData.cxx

Go to the documentation of this file.

#include "larsim/Simulation/LArVoxelData.h"

#include <ostream>

namespace sim {

  //----------------------------------------------------------------------------
  // Constructor; take care of any initializations.
  LArVoxelData::LArVoxelData() : fenergy(0) {}

  //----------------------------------------------------------------------------
  // Destructor.
  LArVoxelData::~LArVoxelData() {}

  //----------------------------------------------------------------------------
  const LArVoxelData::key_type& LArVoxelData::TrackID(const size_type index) const
  {
    const_iterator i = ftrackEnergy.begin();
    std::advance(i, index);
    return (*i).first;
  }

  //----------------------------------------------------------------------------
  const LArVoxelData::mapped_type& LArVoxelData::Energy(const size_type index) const
  {
    const_iterator i = ftrackEnergy.begin();
    std::advance(i, index);
    return (*i).second;
  }

  //----------------------------------------------------------------------------
  void LArVoxelData::Add(const LArVoxelData& other)
  {
    // When we add one voxel's data to another, it becomes impossible
    // to keep the particle<->energy assignments anymore; the most
    // likely reason to add two voxels is because we're adding events
    // to make overlays, and so the particles' track numbers change.

    // So if we're adding another LArVoxelData to this one, move all
    // the energies to "unassigned" in the sum.

    this->RemoveAllTracks();
    fenergy += other.Energy();
  }

  //----------------------------------------------------------------------------
  LArVoxelData& LArVoxelData::operator*=(const double& value)
  {
    // Multiply all energies by the value.
    for (iterator i = ftrackEnergy.begin(), end = ftrackEnergy.end(); i != end; ++i) {
      (*i).second *= value;
    }
    fenergy *= value;

    return (*this);
  }

  //----------------------------------------------------------------------------
  const LArVoxelData operator*(const double& value, const LArVoxelData& data)
  {
    return LArVoxelData(data) *= value;
  }

  //----------------------------------------------------------------------------
  std::ostream& operator<<(std::ostream& output, const LArVoxelData& data)
  {
    output << "Voxel: " << data.VoxelID() << std::endl;

    double unassigned = data.UnassignedEnergy();
    // Display the total energy then the breakdown of
    // the sum.
    output << data.Energy() << " = <ID,E>=";
    for (LArVoxelData::const_iterator i = data.begin(); i != data.end(); ++i) {
      if (i != data.begin()) output << ",";

      output << "<" << (*i).first << "," << (*i).second << ">";
    }
    if (unassigned > 0) output << ",<*," << unassigned << ">";

    return output;
  }

} // namespace sim