BlobGrouping.cxx File Reference

#include "WireCellImg/BlobGrouping.h"
#include "WireCellIface/SimpleCluster.h"
#include "WireCellUtil/NamedFactory.h"
#include <boost/graph/connected_components.hpp>

Go to the source code of this file.

Typedefs
typedef std::unordered_map< WirePlaneLayer_t, cluster_indexed_graph_t >	layer_graphs_t

Functions
	WIRECELL_FACTORY (BlobGrouping, WireCell::Img::BlobGrouping, WireCell::IClusterFilter, WireCell::IConfigurable) using namespace WireCell
static void	fill_blob (layer_graphs_t &lgs, const cluster_indexed_graph_t &grind, IBlob::pointer iblob)
static void	fill_slice (cluster_indexed_graph_t &grind, ISlice::pointer islice)

Typedef Documentation

typedef std::unordered_map<WirePlaneLayer_t, cluster_indexed_graph_t> layer_graphs_t

Definition at line 14 of file BlobGrouping.cxx.

Function Documentation

static void fill_blob ( layer_graphs_t &  lgs, const cluster_indexed_graph_t &  grind, IBlob::pointer  iblob  )

static

Definition at line 37 of file BlobGrouping.cxx.

{
    cluster_node_t nblob{iblob};

    for (auto wvtx : grind.neighbors(nblob)) {
        if (wvtx.code() != 'w') {
            continue;
        }
        auto iwire = std::get<IWire::pointer>(wvtx.ptr);
        auto layer = iwire->planeid().layer();
        auto& lg = lgs[layer];

        for (auto cvtx : grind.neighbors(wvtx)) {
            if (cvtx.code() != 'c') {
                continue;
            }
            //auto ich = std::get<IChannel::pointer>(cvtx.ptr);
            lg.edge(nblob, cvtx);
        }
    }
}

static void fill_slice ( cluster_indexed_graph_t &  grind, ISlice::pointer  islice  )

static

Definition at line 63 of file BlobGrouping.cxx.

{
    layer_graphs_t lgs;

    for (auto other : grind.neighbors(islice)) {
        if (other.code() != 'b') {
            continue;
        }
        IBlob::pointer iblob = std::get<IBlob::pointer>(other.ptr);
        fill_blob(lgs, grind, iblob);
    }
    
    for (auto lgit : lgs) {
        auto& lgrind = lgit.second;
        auto groups = lgrind.groups();
        for (auto& group : groups) {
            // add a "measurement" to the graph
            IChannel::vector* chans = new IChannel::vector;
            IChannel::shared_vector imeas = IChannel::shared_vector(chans);

            for (auto& v : group.second) {
                if (v.code() == 'b') {
                    // (b-m)
                    grind.edge(v.ptr, imeas);
                    continue;
                }
                if (v.code() == 'c') {
                    // (c-m)
                    grind.edge(v.ptr, imeas);
                    chans->push_back(std::get<IChannel::pointer>(v.ptr));
                    continue;
                }
            }
        }
    }
}

WIRECELL_FACTORY	(	BlobGrouping	,
		WireCell::Img::BlobGrouping	,
		WireCell::IClusterFilter	,
		WireCell::IConfigurable
	)