WireGenerator.cxx File Reference

#include "WireCellGen/WireGenerator.h"
#include "WireCellUtil/Intersection.h"
#include "WireCellUtil/NamedFactory.h"
#include <cmath>
#include <iostream>
#include <vector>

Go to the source code of this file.

Classes
class	WireCell::GenWire
struct	SortByIndex

Namespaces
	WireCell

Functions
	WIRECELL_FACTORY (WireGenerator, WireCell::WireGenerator, WireCell::IWireGenerator) using namespace WireCell
static GenWire *	make_wire (int index, const Point &point, const Point &proto, const Ray &bounds)
static void	make_one_plane (IWire::vector &returned_wires, WirePlaneId wpid, const Ray &bounds, const Ray &step)

Function Documentation

static void make_one_plane ( IWire::vector &  returned_wires, WirePlaneId  wpid, const Ray &  bounds, const Ray &  step  )

static

Definition at line 83 of file WireGenerator.cxx.

{
    const Vector xaxis(1,0,0);
    const Point starting_point = step.first;
    const Vector pitch = step.second - starting_point;
    const Vector proto = pitch.cross(xaxis).norm();
    
    std::vector<GenWire*> these_wires;

    int pos_index = 0;
    Point offset = starting_point;
    while (true) {              // go in positive pitch direction
        GenWire* wire = make_wire(pos_index, offset, proto, bounds);
        if (! wire) { break; }
        these_wires.push_back(wire);
        offset = wire->center() + pitch;
        pos_index += 1;
    }

    int neg_index = -1;         // now go in negative pitch direction
    const Vector neg_pitch = -1.0 * pitch;
    offset = these_wires[0]->center() + neg_pitch; // start one below first upward going one
    while (true) {              // go in negative pitch direction
        GenWire* wire = make_wire(neg_index, offset, proto, bounds);
        if (! wire) { break; }
        these_wires.push_back(wire);
        offset = wire->center() + neg_pitch;
        neg_index -= 1;
    }

    // order by index
    std::sort(these_wires.begin(), these_wires.end(), SortByIndex());

    // load in to store and fix up index and plane
    for (size_t ind=0; ind<these_wires.size(); ++ind) {
        GenWire* pwire = these_wires[ind];
        pwire->set_index(ind);
        pwire->set_planeid(wpid);
        returned_wires.push_back(IWire::pointer(pwire));
    }

    //std::cerr << "Made "<<store.size()<<" wires for plane " << plane << std::endl;
    //std::cerr << "step = " << step << std::endl;
    //std::cerr << "bounds = " << bounds << std::endl;
}

static GenWire* make_wire ( int  index, const Point &  point, const Point &  proto, const Ray &  bounds  )

static

Definition at line 57 of file WireGenerator.cxx.

{
    const Point pt1 = point;
    const Point pt2  = pt1 + proto;
    const Ray wireray(pt1, pt2);

    Ray hits;
    int hitmask = box_intersection(bounds, wireray, hits);
    if (3 != hitmask) {
        return 0;
    }
    // ray should point generally towards +Y
    if (hits.first.y() > hits.second.y()) {
        hits = Ray(hits.second, hits.first);
    }
    return new GenWire(WirePlaneId(kUnknownLayer), index, hits);
}

WIRECELL_FACTORY	(	WireGenerator	,
		WireCell::WireGenerator	,
		WireCell::IWireGenerator
	)