LAr2DWaveBase.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2019 CERN for the benefit of the ATLAS collaboration
*/
 
// $Id$
#include "LArRawConditions/LAr2DWaveBase.h"
#include "LArRawConditions/LArCompactSubset.h"
#include <cassert>
 
 
LAr2DWaveBase::LAr2DWaveBase(float timeOffset,
                             float timeBinWidth,
                             const std::vector<std::vector<float> >& wave0,
                             const std::vector<std::vector<float> >& wave1)
  : m_timeOffset(timeOffset),
    m_timeBinWidth(timeBinWidth)
{
  m_waves[0] = wave0;
  m_waves[1] = wave1;
} 
 
 
LAr2DWaveBase::LAr2DWaveBase(float timeOffset,
                             float timeBinWidth,
                             unsigned int nPhases,
                             unsigned int nSamples,
                             const std::vector<float>& wave0,
                             const std::vector<float>& wave1,
                             unsigned int index)
 : m_timeOffset(timeOffset),
   m_timeBinWidth(timeBinWidth)
{
  assert (index + nPhases*nSamples <= wave0.size());
  assert (index + nPhases*nSamples <= wave1.size());
 
  m_waves[0].resize (nPhases);
  m_waves[1].resize (nPhases);
 
  for (unsigned int i = 0; i < nPhases; i++) {
    m_waves[0][i].assign (wave0.begin() + index,
                          wave0.begin() + index + nSamples);
    m_waves[1][i].assign (wave1.begin() + index,
                          wave1.begin() + index + nSamples);
    index += nSamples;
  }
}
 
 
LAr2DWaveBase::LAr2DWaveBase (const LArCompactSubsetConstChannelProxy& other)
{
  if (!other.isValid()) {
    // Invalid source --- do the same as the default ctor.
    m_timeOffset = 0;
    m_timeBinWidth = 25./24.;
  }
  else {
    m_timeOffset = other.timeOffset();
    m_timeBinWidth = other.timeBinWidth();
    size_t chanSize = other.getSize();
    for (unsigned int which = 0; which < nWaves; ++which) {
      m_waves[which].resize (chanSize);
      for (unsigned int i = 0; i < chanSize; i++) {
        LArVectorProxy wave = other.getData (which, i);
        m_waves[which][i].assign (wave.begin(), wave.end());
      }
    }
  }
}
 
 
void LAr2DWaveBase::setFrom (LAr2DWaveBase& other)
{
  m_timeOffset = other.m_timeOffset;
  m_timeBinWidth = other.m_timeBinWidth;
  for (unsigned int which = 0; which < nWaves; ++which)
    m_waves[which].swap (other.m_waves[which]);
}