In This Package:

PerDaqPmtCrate.cc

Go to the documentation of this file.

/*
 *  PerDaqPmtCrate.cc
 *  PerReadoutEvent
 *
 *  Created by Simon Patton on 2/22/11.
 *  Copyright 2011 DayaBay Collaboration. All rights reserved.
 *
 */
#include "PerReadoutEvent/PerDaqPmtCrate.h"

#include "Conventions/Electronics.h"
#include "PerReadoutEvent/PerUintValue.h"

using DayaBay::DaqPmtChannel;
using DayaBay::FeeChannelId;
using DayaBay::DaqPmtCrate;

DaqPmtCrate::PmtChannelPtrList* PerDaqPmtCrate::m_channels = new DaqPmtCrate::PmtChannelPtrList();

PerDaqPmtCrate::PerDaqPmtCrate(const DybDaq::EventReadout* eventReadout) :
  DayaBay::DaqPmtCrate(eventReadout),
  m_returnedPeaks(0),
  m_returnedAdcs(0),
  m_returnedPreAdcs(0),
  m_returnedGains(0),
  m_returnedTdcs(0) {
}
  
PerDaqPmtCrate::PerDaqPmtCrate(const PerDaqPmtCrate& rhs) :
  DayaBay::DaqPmtCrate(rhs),
  m_returnedPeaks(0),
  m_returnedAdcs(0),
  m_returnedPreAdcs(0),
  m_returnedGains(0),
  m_returnedTdcs(0) {
}

PerDaqPmtCrate::~PerDaqPmtCrate() {
  if (0 != m_returnedTdcs) {
    const std::vector<const PerUintValue*>::const_iterator finished = m_returnedTdcs->end();
    for (std::vector<const PerUintValue*>::const_iterator iterator = m_returnedTdcs->begin();
         finished != iterator;
         ++iterator) {
      delete (PerUintValue*)(*iterator);
    }
    delete m_returnedTdcs;
  }
  if (0 != m_returnedGains) {
    const std::vector<const PerUintValue*>::const_iterator finished = m_returnedGains->end();
    for (std::vector<const PerUintValue*>::const_iterator iterator = m_returnedGains->begin();
         finished != iterator;
         ++iterator) {
      delete (PerUintValue*)(*iterator);
    }
    delete m_returnedGains;
  }
  if (0 != m_returnedAdcs) {
    const std::vector<const PerUintValue*>::const_iterator finished = m_returnedAdcs->end();
    for (std::vector<const PerUintValue*>::const_iterator iterator = m_returnedAdcs->begin();
         finished != iterator;
         ++iterator) {
      delete (PerUintValue*)(*iterator);
    }
    delete m_returnedAdcs;
  }
  if (0 != m_returnedPreAdcs) {
    const std::vector<const PerUintValue*>::const_iterator finished = m_returnedPreAdcs->end();
    for (std::vector<const PerUintValue*>::const_iterator iterator = m_returnedPreAdcs->begin();
         finished != iterator;
         ++iterator) {
      delete (PerUintValue*)(*iterator);
    }
    delete m_returnedPreAdcs;
  }
  if (0 != m_returnedPeaks) {
    const std::vector<const PerUintValue*>::const_iterator finished = m_returnedPeaks->end();
    for (std::vector<const PerUintValue*>::const_iterator iterator = m_returnedPeaks->begin();
         finished != iterator;
         ++iterator) {
      delete (PerUintValue*)(*iterator);
    }
    delete m_returnedPeaks;
  }
}

const DaqPmtChannel& PerDaqPmtCrate::channel(const DayaBay::FeeChannelId& channelId) const {
  return DaqPmtCrate::channel(channelId);
}

const DaqPmtCrate::PmtChannelPtrList& PerDaqPmtCrate::channel(int board, int chan) const {
  m_channels->clear();
  m_channels->push_back(const_cast<DayaBay::DaqPmtChannel*>(&(channel(FeeChannelId(board, chan, (Site::Site_t)0, (DetectorId::DetectorId_t)0)))));
  return *m_channels;
}

const PerDaqPmtCrate::AdcPeakValues& PerDaqPmtCrate::peaks(int board, int chan, int highGain) const {
  if (0 != m_returnedPeaks) {
    delete m_returnedPeaks;
  }
  m_returnedPeaks = new std::vector<const PerUintValue*>();

  const DaqPmtChannel& channelToUse = channel(FeeChannelId(board, chan, (Site::Site_t)0, (DetectorId::DetectorId_t)0));
  if (0 != &channelToUse) {
    const unsigned int finished = channelToUse.hitCount();
    for(unsigned int hit = 0; finished != hit; ++hit) {
      if (0 == highGain && !channelToUse.isHighGainAdc(hit)) {
        m_returnedPeaks->push_back(new PerUintValue(channelToUse.peakCycle(hit)));
      } else if (1 == highGain && channelToUse.isHighGainAdc(hit)) {
        m_returnedPeaks->push_back(new PerUintValue(channelToUse.peakCycle(hit)));
      }
    }
  }
  return *m_returnedPeaks;
}

const PerDaqPmtCrate::AdcValues& PerDaqPmtCrate::adcs(int board, int chan, int highGain) const {
  if (0 != m_returnedAdcs) {
    delete m_returnedAdcs;
  }
  m_returnedAdcs = new std::vector<const PerUintValue*>();

  const DaqPmtChannel& channelToUse = channel(FeeChannelId(board, chan, (Site::Site_t)0, (DetectorId::DetectorId_t)0));
  if (0 != &channelToUse) {
    const unsigned int finished = channelToUse.hitCount();
    for(unsigned int hit = 0; finished != hit; ++hit) {
      if (0 == highGain && !channelToUse.isHighGainAdc(hit)) {
        m_returnedAdcs->push_back(new PerUintValue(channelToUse.adc(hit)));
      } else if (1 == highGain && channelToUse.isHighGainAdc(hit)) {
        m_returnedAdcs->push_back(new PerUintValue(channelToUse.adc(hit)));
      }
    }
  }
  return *m_returnedAdcs;
}

const PerDaqPmtCrate::PreAdcValues& PerDaqPmtCrate::preAdcRaws(int board, int chan, int highGain) const {
  if (0 != m_returnedPreAdcs) {
    delete m_returnedPreAdcs;
  }
  m_returnedPreAdcs = new std::vector<const PerUintValue*>();

  const DaqPmtChannel& channelToUse = channel(FeeChannelId(board, chan, (Site::Site_t)0, (DetectorId::DetectorId_t)0));
  if (0 != &channelToUse) {
    const unsigned int finished = channelToUse.hitCount();
    for(unsigned int hit = 0; finished != hit; ++hit) {
      if (0 == highGain && !channelToUse.isHighGainAdc(hit)) {
        m_returnedPreAdcs->push_back(new PerUintValue(channelToUse.preAdcRaw(hit)));
      } else if (1 == highGain && channelToUse.isHighGainAdc(hit)) {
        m_returnedPreAdcs->push_back(new PerUintValue(channelToUse.preAdcRaw(hit)));
      }
    }
  }
  return *m_returnedPreAdcs;
}

const PerDaqPmtCrate::GainValues& PerDaqPmtCrate::gains(int board, int chan) const {
  if (0 != m_returnedGains) {
    delete m_returnedGains;
  }
  m_returnedGains = new std::vector<const PerUintValue*>();

  const DaqPmtChannel& channelToUse = channel(FeeChannelId(board, chan, (Site::Site_t)0, (DetectorId::DetectorId_t)0));
  if (0 != &channelToUse) {
    const unsigned int finished = channelToUse.hitCount();
    for(unsigned int hit = 0; finished != hit; ++hit) {
      if (!channelToUse.isHighGainAdc(hit)) {
        m_returnedGains->push_back(new PerUintValue(0));
      } else {
        m_returnedGains->push_back(new PerUintValue(1));
      }
    }
  }
  return *m_returnedGains;
}

const PerDaqPmtCrate::TdcValues& PerDaqPmtCrate::tdcs(int board, int chan, int highGain) const {
  if (0 != m_returnedTdcs) {
    delete m_returnedTdcs;
  }
  m_returnedTdcs = new std::vector<const PerUintValue*>();

  const DaqPmtChannel& channelToUse = channel(FeeChannelId(board, chan, (Site::Site_t)0, (DetectorId::DetectorId_t)0));
  if (0 != &channelToUse) {
    const unsigned int finished = channelToUse.hitCount();
    for(unsigned int hit = 0; finished != hit; ++hit) {
      if (0 == highGain && !channelToUse.isHighGainAdc(hit)) {
        m_returnedTdcs->push_back(new PerUintValue(channelToUse.tdc(hit)));
      } else if (1 == highGain && channelToUse.isHighGainAdc(hit)) {
        m_returnedTdcs->push_back(new PerUintValue(channelToUse.tdc(hit)));
      }
    }
  }
  return *m_returnedTdcs;
}

---