In This Package:

DaqPmtChannel.cc

Go to the documentation of this file.

/*
 *  DaqPmtChannel.cc
 *  DaqEvent
 *
 *  Created by Simon Patton on 9/2/2010.
 *  Copyright 2010 DayaBay Collaboration. All rights reserved.
 *
 */
#include "Event/DaqPmtChannel.h"

#include "FeeReadoutFormat/FeeHead.h"
#include "FeeReadoutFormat/FeeHit.h"
#include "FeeReadoutFormat/FeeReadout.h"
#include "FeeReadoutFormat/FeeTraits.h"

using DayaBay::DaqPmtChannel;
using DybDaq::FeeHead;
using DybDaq::FeeHit;
using DybDaq::FeeReadout;

DaqPmtChannel::DaqPmtChannel(const DayaBay::FeeChannelId& channelId,
                             FeeReadout& feeReadout) :
    m_feeHits(0),
    m_feeReadout(&feeReadout),
    m_channelId(channelId) {
}

DaqPmtChannel::DaqPmtChannel(const FeeHit& feeHit,
                             const DayaBay::FeeChannelId& channelId) :
    m_feeHits(0),
    m_feeReadout(0),
    m_channelId(channelId) {
    addHit(feeHit);
}

DaqPmtChannel::~DaqPmtChannel() {
    if (0 != m_feeHits) {
        delete m_feeHits;
    }
}

unsigned int DaqPmtChannel::hitCount() const {
    return feeHits().size();
}

unsigned int DaqPmtChannel::peakCycle(const unsigned int hit) const {
    return feeHits().at(hit)->peakCycle();
}

bool DaqPmtChannel::isHighGainAdc(const unsigned int hit) const {
    return !(feeHits().at(hit)->isHighRangeAdc());
}

unsigned int DaqPmtChannel::preAdcRaw(const unsigned int hit) const {
    return feeHits().at(hit)->preAdc();
}

float DaqPmtChannel::preAdcAvg(const unsigned int hit) const {
    if(feeHits().at(hit)->feeTraits().majorVersion() == 0
       && feeHits().at(hit)->feeTraits().minorVersion() < 2){
        return float(this->preAdcRaw(hit));
    }
    return float(this->preAdcRaw(hit)) / 4.0;
}

unsigned int DaqPmtChannel::adc(const unsigned int hit) const {
    return feeHits().at(hit)->adc();
}

float DaqPmtChannel::deltaAdc(const unsigned int hit) const {
    return this->adc(hit) - this->preAdcAvg(hit);
}

unsigned int DaqPmtChannel::tdc(const unsigned int hit) const {
    return feeHits().at(hit)->tdc();
}

const DaqPmtChannel::FeeHitPtrList& DaqPmtChannel::feeHits() const {
    if (0 == m_feeHits) {
        FeeHitPtrList** hits = const_cast<FeeHitPtrList**>(&m_feeHits);
        *(hits) = new FeeHitPtrList();
    }
    return *m_feeHits;
}

void DaqPmtChannel::addHit(const unsigned int hitNumber,
                           const unsigned int peakCycle,
                           const bool isHighGainAdc,
                           const unsigned int pedestal,
                           const unsigned int charge,
                           const unsigned int time) {
    const FeeHit& feeHit = m_feeReadout->addHit(m_channelId.connector() - 1,
                                                peakCycle,
                                                !isHighGainAdc,
                                                pedestal,
                                                charge,
                                                hitNumber,
                                                time);
    addHit(feeHit);
}

void DaqPmtChannel::addHit(const FeeHit& feeHit) {
    // TODO: Make sure the hits are in order! (eventually)
    feeHits();
    m_feeHits->push_back(&feeHit);
}

---