WireCell::Root::HistFrameSink Class Reference

#include <HistFrameSink.h>

Inheritance diagram for WireCell::Root::HistFrameSink:

Public Member Functions
	HistFrameSink ()
virtual	~HistFrameSink ()
virtual bool	operator() (const IFrame::pointer &frame)
	Frame sink interface. More...
virtual void	configure (const WireCell::Configuration &config)
	Configurable interface. More...
virtual WireCell::Configuration	default_configuration () const
	Optional, override to return a hard-coded default configuration. More...
Public Member Functions inherited from WireCell::IFrameSink
virtual	~IFrameSink ()
virtual std::string	signature ()
Public Member Functions inherited from WireCell::ISinkNode< IFrame >
virtual	~ISinkNode ()
virtual bool	operator() (const boost::any &anyin)
virtual std::vector< std::string >	input_types ()
Public Member Functions inherited from WireCell::ISinkNodeBase
virtual	~ISinkNodeBase ()
virtual NodeCategory	category ()
	Return the behavior category type. More...
Public Member Functions inherited from WireCell::INode
virtual	~INode ()
virtual int	concurrency ()
virtual std::vector< std::string >	output_types ()
virtual void	reset ()
Public Member Functions inherited from WireCell::IComponent< INode >
virtual	~IComponent ()
Public Member Functions inherited from WireCell::Interface
virtual	~Interface ()
Public Member Functions inherited from WireCell::IConfigurable
virtual	~IConfigurable ()
Public Member Functions inherited from WireCell::IComponent< IConfigurable >
virtual	~IComponent ()

Private Attributes
std::string	m_filepat
std::string	m_anode_tn
IAnodePlane::pointer	m_anode
double	m_units

Additional Inherited Members
Public Types inherited from WireCell::IFrameSink
typedef std::shared_ptr< IFrameSink >	pointer
Public Types inherited from WireCell::ISinkNode< IFrame >
typedef IFrame	input_type
typedef std::shared_ptr< const IFrame >	input_pointer
Public Types inherited from WireCell::ISinkNodeBase
typedef std::shared_ptr< ISinkNodeBase >	pointer
Public Types inherited from WireCell::INode
enum	NodeCategory {   unknown, sourceNode, sinkNode, functionNode,   queuedoutNode, joinNode, splitNode, faninNode,   fanoutNode, multioutNode, hydraNode }
Public Types inherited from WireCell::IComponent< INode >
typedef std::shared_ptr< INode >	pointer
	Access subclass facet by pointer. More...
typedef std::vector< pointer >	vector
	Vector of shared pointers. More...
Public Types inherited from WireCell::Interface
typedef std::shared_ptr< Interface >	pointer
Public Types inherited from WireCell::IComponent< IConfigurable >
typedef std::shared_ptr< IConfigurable >	pointer
	Access subclass facet by pointer. More...
typedef std::vector< pointer >	vector
	Vector of shared pointers. More...

Detailed Description

Definition at line 19 of file HistFrameSink.h.

Constructor & Destructor Documentation

Root::HistFrameSink::HistFrameSink

(

)

Definition at line 20 of file HistFrameSink.cxx.

    : m_filepat("histframe-%02d.root")
    , m_anode_tn("AnodePlane")
    , m_units(1.0)
{
}

Root::HistFrameSink::~HistFrameSink ( )

virtual

Definition at line 27 of file HistFrameSink.cxx.

{
}

Member Function Documentation

void Root::HistFrameSink::configure ( const WireCell::Configuration &  config )

virtual

Configurable interface.

Implements WireCell::IConfigurable.

Definition at line 41 of file HistFrameSink.cxx.

{
    m_filepat = get<std::string>(cfg, "filename", m_filepat);
    m_units = get<double>(cfg, "units", m_units);
    m_anode_tn = get<std::string>(cfg, "anode", m_anode_tn);
    m_anode = Factory::lookup_tn<IAnodePlane>(m_anode_tn);
    if (!m_anode) {
        cerr << "Root::HistFrameSink: failed to get anode: \"" << m_anode_tn << "\"\n";
        return;
    }
    cerr << "Root::HistFrameSink: configured with: "
         << "file:" << m_filepat << ", "
         << "units:" << m_units << ", "
         << "anode:" << m_anode_tn << endl;
}

WireCell::Configuration Root::HistFrameSink::default_configuration ( ) const

virtual

Optional, override to return a hard-coded default configuration.

Reimplemented from WireCell::IConfigurable.

Definition at line 32 of file HistFrameSink.cxx.

{
    Configuration cfg;
    cfg["filename"] = m_filepat;
    cfg["anode"] = m_anode_tn;
    cfg["units"] = units::mV;
    return cfg;
}

bool Root::HistFrameSink::operator() ( const IFrame::pointer &  frame )

virtual

Frame sink interface.

Implements WireCell::ISinkNode< IFrame >.

Definition at line 59 of file HistFrameSink.cxx.

{
    if (!frame) {
        cerr << "Root::HistFrameSink: no frame\n";
        return true;
    }

    std::string fname = Form(m_filepat.c_str(), frame->ident());
    TFile* file = TFile::Open(fname.c_str(), "recreate");
    

    typedef std::tuple< ITrace::vector, std::vector<int>, std::vector<int> > tct_tuple;
    std::unordered_map<int, tct_tuple> perplane; 

    // collate traces into per plane and also calculate bounds 
    ITrace::shared_vector traces = frame->traces();
    for (ITrace::pointer trace : *traces) {
        int ch = trace->channel();
        auto wpid = m_anode->resolve(ch);
        int wpident = wpid.ident(); 
        double tmin = trace->tbin();
        double tlen = trace->charge().size();

        auto& tct = perplane[wpident];
        get<0>(tct).push_back(trace);
        get<1>(tct).push_back(ch);
        get<2>(tct).push_back(tmin);
        get<2>(tct).push_back(tmin+tlen);

        if (wpident < 0) {
            cerr << "Channel "<<ch<<" has illegal wire plane ident: " << wpid << endl;
        }
    }

    const double t0 = frame->time();
    const double tick = frame->tick();

    for (auto& thisplane : perplane) {
        int wpident = thisplane.first;
        auto& tct = thisplane.second;
        auto& traces = get<0>(tct);
        auto& chans = get<1>(tct);
        auto chmm = std::minmax_element(chans.begin(), chans.end());
        auto& tbins = get<2>(tct);
        auto tbmm = std::minmax_element(tbins.begin(), tbins.end());


        const double tmin = t0 + tick*(*tbmm.first);
        const double tmax = t0 + tick*(*tbmm.second);
        const int ntbins = (*tbmm.second)-(*tbmm.first);

        const int chmin = round(*chmm.first);
        const int chmax = round(*chmm.second + 1);
        const int nchbins = chmax - chmin;
        
        TH2F* hist = new TH2F(Form("plane%d", wpident),
                              Form("Plane %d", wpident),
                              ntbins, tmin/units::us, tmax/units::us,
                              nchbins, chmin, chmax);
        hist->SetDirectory(file); // needed?
        hist->SetXTitle("time [us]");
        hist->SetYTitle("channel");

        double qtot = 0;
        int nbins_tot = 0;
        for (auto& trace : traces) {
            double fch = trace->channel() + 0.5; // make sure we land in bin-center.
            int tbin = trace->tbin();
            auto& charge = trace->charge();
            int nbins = charge.size();
            nbins_tot += nbins;
            for (int ibin=0; ibin<nbins; ++ibin) {
                const double t = t0 + (tick)*(tbin+ibin+0.5); // 0.5 to land in bin-center
                hist->Fill(t/units::us, fch, charge[ibin]/m_units);
                qtot += charge[ibin];
            }
        }

        cerr << wpident
             << " ntraces:" << traces.size() << " "
             << " nsamples:" << nbins_tot << " "
             << " qtot:" << qtot/m_units << " "
             << " qunit:" << m_units << " "
             << " integ:" << hist->Integral()
             << " min:" << hist->GetMinimum()
             << " max:" << hist->GetMaximum()
             << " chan:"<<nchbins<<"[" << chmin << "," << chmax << "] "
             << " time:"<<ntbins<<"[" << tmin/units::us << "," << tmax/units::us <<"]us\n";


        hist->Write();
    }
    file->Close();
    delete file;
    file = nullptr;
    return true;
}

Member Data Documentation

IAnodePlane::pointer WireCell::Root::HistFrameSink::m_anode

private

Definition at line 35 of file HistFrameSink.h.

std::string WireCell::Root::HistFrameSink::m_anode_tn

private

Definition at line 34 of file HistFrameSink.h.

std::string WireCell::Root::HistFrameSink::m_filepat

private

Definition at line 34 of file HistFrameSink.h.

double WireCell::Root::HistFrameSink::m_units

private

Definition at line 36 of file HistFrameSink.h.

---

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

• wire-cell-build/root/inc/WireCellRoot/HistFrameSink.h
• wire-cell-build/root/src/HistFrameSink.cxx