Fourdee.cxx File Reference

#include "WireCellGen/Fourdee.h"
#include "WireCellGen/FrameUtil.h"
#include "WireCellUtil/NamedFactory.h"
#include "WireCellUtil/ConfigManager.h"
#include "WireCellUtil/ExecMon.h"
#include "WireCellGen/GenPipeline.h"

Go to the source code of this file.

Classes
class	NoiseAdderProc

Functions
	WIRECELL_FACTORY (FourDee, WireCell::Gen::Fourdee, WireCell::IApplication, WireCell::IConfigurable) using namespace std
void	dump (const IFrame::pointer frame)
template<typename DEPOS >
void	dump (DEPOS &depos)

Function Documentation

void dump

(

const IFrame::pointer

frame

)

Definition at line 120 of file Fourdee.cxx.

{
    if (!frame) {
        cerr << "Fourdee: dump: empty frame\n";
        return;
    }

    for (auto tag: frame->frame_tags()) {
        const auto& tlist = frame->tagged_traces(tag);
        cerr << "Fourdee: frame tag: " << tag << " with " << tlist.size() << " traces\n";
    }
    for (auto tag: frame->trace_tags()) {
        const auto& tlist = frame->tagged_traces(tag);
        cerr << "Fourdee: trace tag: " << tag << " with " << tlist.size() << " traces\n";
    }

    auto traces = frame->traces();
    const int ntraces = traces->size();

    if (ntraces <= 0) {
        cerr << "Fourdee: dump: no traces\n";
        return;
    }

    std::vector<int> tbins, tlens;
    for (auto trace : *traces) {
        const int tbin = trace->tbin();
        tbins.push_back(tbin);
        tlens.push_back(tbin+trace->charge().size());
    }

    int tmin = *(std::minmax_element(tbins.begin(), tbins.end()).first);
    int tmax = *(std::minmax_element(tlens.begin(), tlens.end()).second);

    cerr << "frame: #" << frame->ident()
         << " @" << frame->time()/units::ms
         << "ms with " << ntraces << " traces, tbins in: "
         << "[" << tmin << "," << tmax << "]"
         << endl;
}

template<typename DEPOS >

void dump

(

DEPOS &

depos

)

Definition at line 162 of file Fourdee.cxx.

{
    if (depos.empty() or (depos.size() == 1 and !depos[0])) {
        std::cerr << "Fourdee::dump: empty depos set\n";
        return;
    }

    std::vector<double> t,x,y,z;
    double qtot = 0.0;
    double qorig = 0.0;

    for (auto depo : depos) {
        if (!depo) {
            cerr << "Gen::Fourdee: null depo" << endl;
            break;
        }
        auto prior = depo->prior();
        if (!prior) {
            cerr << "Gen::Fourdee: null prior depo" << endl;
        }
        else {
            qorig += prior->charge();
        }
        t.push_back(depo->time());
        x.push_back(depo->pos().x());
        y.push_back(depo->pos().y());
        z.push_back(depo->pos().z());
        qtot += depo->charge();
    }

    auto tmm = std::minmax_element(t.begin(), t.end());
    auto xmm = std::minmax_element(x.begin(), x.end());
    auto ymm = std::minmax_element(y.begin(), y.end());
    auto zmm = std::minmax_element(z.begin(), z.end());
    const int ndepos = depos.size();

    std::cerr << "Gen::FourDee: drifted " << ndepos << ", extent:\n"
              << "\tt in [ " << (*tmm.first)/units::us << "," << (*tmm.second)/units::us << "]us,\n"
              << "\tx in [" << (*xmm.first)/units::mm << ","<<(*xmm.second)/units::mm<<"]mm,\n"
              << "\ty in [" << (*ymm.first)/units::mm << ","<<(*ymm.second)/units::mm<<"]mm,\n"
              << "\tz in [" << (*zmm.first)/units::mm << ","<<(*zmm.second)/units::mm<<"]mm,\n"
              << "\tQtot=" << qtot/units::eplus << ", <Qtot>=" << qtot/ndepos/units::eplus << " "
              << " Qorig=" << qorig/units::eplus
              << ", <Qorig>=" << qorig/ndepos/units::eplus << " electrons" << std::endl;

        
}

WIRECELL_FACTORY	(	FourDee	,
		WireCell::Gen::Fourdee	,
		WireCell::IApplication	,
		WireCell::IConfigurable
	)