wcls::CookedFrameSink Class Reference

#include <CookedFrameSink.h>

Inheritance diagram for wcls::CookedFrameSink:

Public Member Functions
	CookedFrameSink ()
virtual	~CookedFrameSink ()
virtual void	produces (art::ProducesCollector &collector)
	IArtEventVisitor. More...
virtual void	visit (art::Event &event)
	Implement to visit an Art event. More...
virtual bool	operator() (const WireCell::IFrame::pointer &frame)
	IFrameSink. More...
virtual WireCell::Configuration	default_configuration () const
	IConfigurable. More...
virtual void	configure (const WireCell::Configuration &config)
Public Member Functions inherited from wcls::IArtEventVisitor
virtual	~IArtEventVisitor ()

Private Attributes
WireCell::IFrame::pointer	m_frame
WireCell::IAnodePlane::pointer	m_anode
std::vector< std::string >	m_frame_tags
int	m_nticks

Detailed Description

Definition at line 20 of file CookedFrameSink.h.

Constructor & Destructor Documentation

CookedFrameSink::CookedFrameSink

(

)

Definition at line 21 of file CookedFrameSink.cxx.

    : m_frame(nullptr)
    , m_nticks(0)
{
}

CookedFrameSink::~CookedFrameSink ( )

virtual

Definition at line 27 of file CookedFrameSink.cxx.

{
}

Member Function Documentation

void CookedFrameSink::configure ( const WireCell::Configuration &  config )

virtual

Definition at line 43 of file CookedFrameSink.cxx.

{
    const std::string anode_tn = cfg["anode"].asString();
    if (anode_tn.empty()) {
        THROW(ValueError() << errmsg{"CookedFrameSink requires an anode plane"});
    }

    // this throws if not found
    m_anode = Factory::find_tn<IAnodePlane>(anode_tn);

    auto jtags = cfg["frame_tags"];
    std::cerr << "CookedFrameSink: saving " << jtags.size() << " tags\n";
    for (auto jtag : jtags) {
        std::string tag = jtag.asString();
        std::cerr << "\t" << tag << "\n";
        m_frame_tags.push_back(tag);
    }
    m_nticks = get(cfg, "nticks", m_nticks);
}

WireCell::Configuration CookedFrameSink::default_configuration ( ) const

virtual

IConfigurable.

Definition at line 32 of file CookedFrameSink.cxx.

{
    Configuration cfg;
    cfg["anode"] = "AnodePlane";
    // frames to output
    cfg["frame_tags"][0] = "gauss";
    cfg["frame_tags"][1] = "wiener";
    cfg["nticks"] = m_nticks; // if nonzero, force number of ticks in output waveforms.
    return cfg;
}

bool CookedFrameSink::operator() ( const WireCell::IFrame::pointer &  frame )

virtual

IFrameSink.

Definition at line 152 of file CookedFrameSink.cxx.

{
    // set an IFrame based on last visited event.
    m_frame = frame;
    return true;
}

void CookedFrameSink::produces ( art::ProducesCollector &  collector )

virtual

IArtEventVisitor.

Reimplemented from wcls::IArtEventVisitor.

Definition at line 63 of file CookedFrameSink.cxx.

{
    for (auto tag : m_frame_tags) {
        std::cerr << "CookedFrameSink: promising to produce recob::Wires named \"" << tag << "\"\n";
        collector.produces< std::vector<recob::Wire> >(tag);
    }
}

void CookedFrameSink::visit ( art::Event &  event )

virtual

Implement to visit an Art event.

Implements wcls::IArtEventVisitor.

Definition at line 91 of file CookedFrameSink.cxx.

{
    if (!m_frame) {
        std::cerr << "CookedFrameSink: I have no frame to save to art::Event\n";
        return;
    }

    std::cerr << "CookedFrameSink: got " << m_frame->traces()->size() << " total traces\n";

    for (auto tag : m_frame_tags) {

        auto traces = tagged_traces(m_frame, tag);
        if (traces.empty()) {
            std::cerr << "CookedFrameSink: no traces for tag \"" << tag << "\"\n";
            continue;
        }


        std::unique_ptr<std::vector<recob::Wire> > outwires(new std::vector<recob::Wire>);

        // what about the frame's time() and ident()?

        for (const auto& trace : traces) {

            const int tbin = trace->tbin();
            const int chid = trace->channel();
            const auto& charge = trace->charge();

            //std::cerr << tag << ": chid=" << chid << " tbin=" << tbin << " nq=" << charge.size() << std::endl;

            // enforce number of ticks if we are so configured.
            size_t ncharge = charge.size();
            int nticks = tbin + ncharge;
            if (m_nticks) {     // force output waveform size
                if (m_nticks < nticks) {
                    ncharge  = m_nticks - tbin;
                }
                nticks = m_nticks;
            }
            recob::Wire::RegionsOfInterest_t roi(nticks);
            roi.add_range(tbin, charge.begin(), charge.begin() + ncharge);

                // FIXME: the current assumption in this code is that LS channel
                // numbers are identified with WCT channel IDs.
                // Fact: the plane view for the ICARUS induction-1 is "geo::kY",
                // instead of "geo::kU"
                auto const& gc = *lar::providerFrom<geo::Geometry>();
                auto view = gc.View(chid);

            // what about those pesky channel map masks?
            // they are dropped for now.

            outwires->emplace_back(recob::Wire(roi, chid, view));
        }
        std::cerr << "CookedFrameSink saving " << outwires->size() << " recob::Wires named \""<<tag<<"\"\n";
        event.put(std::move(outwires), tag);
    }

    m_frame = nullptr;
}

Member Data Documentation

WireCell::IAnodePlane::pointer wcls::CookedFrameSink::m_anode

private

Definition at line 40 of file CookedFrameSink.h.

WireCell::IFrame::pointer wcls::CookedFrameSink::m_frame

private

Definition at line 39 of file CookedFrameSink.h.

std::vector<std::string> wcls::CookedFrameSink::m_frame_tags

private

Definition at line 41 of file CookedFrameSink.h.

int wcls::CookedFrameSink::m_nticks

private

Definition at line 42 of file CookedFrameSink.h.

---

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

• CookedFrameSink.h
• CookedFrameSink.cxx