ElecResponse.cxx

Go to the documentation of this file.

#include "WireCellGen/ElecResponse.h"

#include "WireCellUtil/NamedFactory.h"
#include "WireCellUtil/Response.h"

#include <iostream>

WIRECELL_FACTORY(ElecResponse, WireCell::Gen::ElecResponse,
                 WireCell::IWaveform, WireCell::IConfigurable)

using namespace std;
using namespace WireCell;

Gen::ElecResponse::ElecResponse(int nticks, double t0, double gain, double shaping, double postgain, double tick)
{
    m_cfg["gain"] = gain;
    m_cfg["shaping"] = shaping;
    m_cfg["postgain"] = postgain;
    m_cfg["start"] = t0;
    m_cfg["tick"] = tick;
    m_cfg["nticks"] = nticks;
}
WireCell::Configuration Gen::ElecResponse::default_configuration() const
{
    return m_cfg;
}
void Gen::ElecResponse::configure(const WireCell::Configuration& cfg)
{
    m_cfg = cfg;

    Response::ColdElec ce(m_cfg["gain"].asDouble(), m_cfg["shaping"].asDouble());
    const int nbins = m_cfg["nticks"].asInt();
    const double t0 = waveform_start();
    const double tick = waveform_period();
    Binning bins(nbins, t0, t0+nbins*tick);
    m_wave = ce.generate(bins);
    Waveform::scale(m_wave, m_cfg["postgain"].asDouble());
}

double Gen::ElecResponse::waveform_start() const
{
    return m_cfg["start"].asDouble();
}

double Gen::ElecResponse::waveform_period() const
{
    return m_cfg["tick"].asDouble();
}

const IWaveform::sequence_type& Gen::ElecResponse::waveform_samples() const
{
    return m_wave;
}