wct::sigproc Namespace Reference

Functions
void	dump_frame (WireCell::IFrame::pointer frame)
void	raster (WireCell::Array::array_xxf &block, WireCell::ITrace::vector traces, const std::vector< int > &channels)
int	maxcount_baseline (const WireCell::ITrace::vector &traces, const WireCell::Binning &binning=WireCell::Binning(4096, 0, 4096))
WireCell::ITrace::vector	tagged_traces (WireCell::IFrame::pointer frame, WireCell::IFrame::tag_t tag)

Function Documentation

void wct::sigproc::dump_frame

(

WireCell::IFrame::pointer

frame

)

Definition at line 10 of file FrameUtils.cxx.

{
    auto traces = frame->traces();
    const size_t ntraces = traces->size();
    std::vector<double> means, rmses, lengths, tbins;
    for (auto trace : *traces) {
        auto const& charge = trace->charge();

        auto mr = Waveform::mean_rms(charge);
        means.push_back(mr.first);
        rmses.push_back(mr.second*mr.second);
        const int nsamps = charge.size();
        lengths.push_back(nsamps);
        tbins.push_back(trace->tbin());
        
        if (std::isnan(mr.second)) {
            std::cerr << "Frame: channel " << trace->channel() << " rms is NaN\n";
        }

        for (int ind=0; ind<nsamps; ++ind) {
            float val = charge[ind];
            if (std::isnan(val)) {
                std::cerr << "Frame: channel " << trace->channel() << " sample " << ind << " is NaN\n";
            }
            if (std::isinf(val)) {
                std::cerr << "Frame: channel " << trace->channel() << " sample " << ind << " is INF\n";
            }
        }
    }
    double meanmean = Waveform::sum(means)/ntraces;
    double totrms = sqrt(Waveform::sum(rmses));
    double meanlen = Waveform::sum(lengths)/ntraces;
    double meantbin = Waveform::sum(tbins)/ntraces;
    std::cerr << "Frame: " << ntraces << " traces,"
              << " <mean>=" << meanmean
              << " TotRMS=" << totrms
              << " <tbin>=" << meantbin
              << " <len>=" << meanlen
              << std::endl;
    for (auto it : frame->masks()) {
        std::cerr << "\t" << it.first << " : " << it.second.size() << std::endl;
    }
}

int wct::sigproc::maxcount_baseline	(	const WireCell::ITrace::vector &	traces,
		const WireCell::Binning &	binning = WireCell::Binning(4096, 0, 4096)
	)

void wct::sigproc::raster	(	WireCell::Array::array_xxf &	block,
		WireCell::ITrace::vector	traces,
		const std::vector< int > &	channels
	)

Definition at line 55 of file FrameUtils.cxx.

{
    const size_t nchannels = channels.size();
    std::unordered_map<int, size_t> ch2col;
    for (size_t ind=0; ind<nchannels; ++ind) {
        ch2col[channels[ind]] = ind;
    }

    const size_t ncols = block.cols();
    for (auto trace : traces) {
        const auto& samples = trace->charge();
        const size_t nsamples = samples.size();
        const size_t tbin = trace->tbin();        

        if (tbin >= ncols) {    // underflow impossible as they are unsigned.
            continue;           // trace is off screen
        }

        // find the row to fill
        const int ch = trace->channel();
        const auto chit = ch2col.find(ch);
        if (chit == ch2col.end()) {
            continue;
        }
        const size_t irow = chit->second;
        WireCell::Array::array_xf tofill(ncols-tbin);
        const size_t maxbin = std::min(nsamples, ncols-tbin);
        for (size_t bin=0; bin<maxbin; ++bin) {
            tofill[bin] = samples[bin];
        }
        block.block(irow, tbin, 1, ncols-tbin).row(0) += tofill;
    }
}

WireCell::ITrace::vector wct::sigproc::tagged_traces	(	WireCell::IFrame::pointer	frame,
		WireCell::IFrame::tag_t	tag
	)