WireCell::RayGrid::Tiling Class Reference

#include <RayTiling.h>

Public Member Functions
	Tiling (const Coordinates &coords)
Activity	projection (const Blob &blob, const Activity &activity)
blobs_t	operator() (const Activity &activity)
blobs_t	operator() (const blobs_t &prior, const Activity &activity)

Private Attributes
const Coordinates &	m_coords

Detailed Description

Definition at line 136 of file RayTiling.h.

Constructor & Destructor Documentation

Tiling::Tiling

(

const Coordinates &

coords

)

Definition at line 209 of file RayTiling.cxx.

    : m_coords(coords)
{
}

Member Function Documentation

blobs_t Tiling::operator()

(

const Activity &

activity

)

Definition at line 216 of file RayTiling.cxx.

{
    auto strips = activity.make_strips();
    const size_t nstrips = strips.size();
    blobs_t ret(nstrips);
    for (size_t ind=0; ind<nstrips; ++ind) {
        ret[ind].add(m_coords, strips[ind]);
    }
    return ret;
}

blobs_t Tiling::operator()	(	const blobs_t &	prior,
		const Activity &	activity
	)

Definition at line 301 of file RayTiling.cxx.

{
    blobs_t ret;

    for (const auto& blob : prior_blobs) {
        Activity proj = projection(blob, activity);
        if (proj.empty()) {
            continue;
        }
        auto strips = proj.make_strips();
        for (auto strip : strips) {
            Blob newblob = blob; // copy
            newblob.add(m_coords, strip);
            if (newblob.corners().empty()) {
                continue;
            }
            ret.push_back(newblob);
        }
    }

    return ret;
}

Activity Tiling::projection	(	const Blob &	blob,
		const Activity &	activity
	)

Definition at line 227 of file RayTiling.cxx.

{
    // special case.  A single layer blob effectively extends to
    // infinity along the ray direction so any possible activity
    // projects.
    if (blob.strips().size() == 1) {
        return activity;
    }

    const double pitch_mag = m_coords.pitch_mags()[activity.layer()];

    // find extreme pitches
    std::vector<double> pitches;
    const auto corners = blob.corners();
    if (corners.empty()) {
        return Activity(activity.layer());
    }
    for (const auto& c : blob.corners()) {
        const double p = m_coords.pitch_location(c.first, c.second, activity.layer());
        pitches.push_back(p);
    }
    if (pitches.empty()) {
        return Activity(activity.layer());
    }

    auto pbeg = pitches.begin();
    auto pend = pitches.end();

    const auto mm = std::minmax_element(pbeg, pend);
    int pind1 = std::floor((*mm.first)/pitch_mag);
    int pind2 = std::ceil((*mm.second)/pitch_mag);

    const int apind1 = activity.pitch_index(activity.begin());
    const int apind2 = activity.pitch_index(activity.end());

    if (pind2 <= apind1 or pind1 >= apind2) {
        return Activity(activity.layer());
    }

    pind1 = std::max(pind1, activity.pitch_index(activity.begin()));
    pind2 = std::min(pind2, activity.pitch_index(activity.end()));
    
    Activity ret = activity.subspan(pind1, pind2);
    return ret;
}

Member Data Documentation

const Coordinates& WireCell::RayGrid::Tiling::m_coords

private

Definition at line 152 of file RayTiling.h.

---

The documentation for this class was generated from the following files:

• wire-cell-build/util/inc/WireCellUtil/RayTiling.h
• wire-cell-build/util/src/RayTiling.cxx