MagnifySink.cxx File Reference

#include "WireCellRoot/MagnifySink.h"
#include "WireCellIface/ITrace.h"
#include "TH1F.h"
#include "TH2F.h"
#include "TH2I.h"
#include "TFile.h"
#include "TTree.h"
#include "WireCellUtil/NamedFactory.h"
#include "WireCellIface/FrameTools.h"
#include <vector>
#include <string>

Go to the source code of this file.

Typedefs
typedef std::unordered_set< std::string >	string_set_t

Functions
	WIRECELL_FACTORY (MagnifySink, WireCell::Root::MagnifySink, WireCell::IFrameFilter, WireCell::IConfigurable) using namespace WireCell
string_set_t	getset (const WireCell::Configuration &cfg)
std::vector< WireCell::Binning >	collate_byplane (const ITrace::vector &traces, const IAnodePlane::pointer anode, ITrace::vector byplane[])

Typedef Documentation

typedef std::unordered_set<std::string> string_set_t

Definition at line 113 of file MagnifySink.cxx.

Function Documentation

std::vector<WireCell::Binning> collate_byplane	(	const ITrace::vector &	traces,
		const IAnodePlane::pointer	anode,
		ITrace::vector	byplane[]
	)

Definition at line 143 of file MagnifySink.cxx.

{
    std::vector<int> uvwt[4];
    for (auto trace : traces) {
        const int chid = trace->channel();
        auto wpid = anode->resolve(chid);
        const int iplane = wpid.index();
        if (iplane<0 || iplane>=3) {
            THROW(RuntimeError() << errmsg{"Illegal wpid"});
        }
        uvwt[iplane].push_back(chid);
        byplane[iplane].push_back(trace);
        uvwt[3].push_back(trace->tbin());
        uvwt[3].push_back(trace->tbin() + trace->charge().size());
    }

    std::vector<Binning> binnings(4);
    for (int ind=0; ind<4; ++ind) {
        auto const& one = uvwt[ind];
        // std::cerr << "[wgu] get ind: " << ind << " size: " << one.size() << std::endl;
        if (one.empty()) {
            // THROW(ValueError() << errmsg{"MagnifySink: bogus bounds"});
            std::cerr << "[wgu] plane: " << ind << " has not traces. " << std::endl;
        }
        else{    
            auto mme = std::minmax_element(one.begin(), one.end());
            const int vmin = *mme.first;
            const int vmax = *mme.second;
            if (ind == 3) {
                const int n = vmax - vmin;
                // binnings.push_back(Binning(n, vmin, vmax););
                binnings.at(ind) = Binning(n, vmin, vmax);
            }
            else {
                // Channel-centered binning
                const double diff = vmax - vmin;
                // binnings.push_back(Binning(diff+1, vmin-0.5, vmax+0.5));
                binnings.at(ind) = Binning(diff+1, vmin-0.5, vmax+0.5);
            }
        }
    }
    return binnings;
}

string_set_t getset

(

const WireCell::Configuration &

cfg

)

Definition at line 114 of file MagnifySink.cxx.

{
    string_set_t ret;
    for (auto jone : cfg) {
        ret.insert(jone.asString());
    }
    return ret;
}

WIRECELL_FACTORY	(	MagnifySink	,
		WireCell::Root::MagnifySink	,
		WireCell::IFrameFilter	,
		WireCell::IConfigurable
	)