TwoD_Processing.h

Go to the documentation of this file.

// -*- mode: C++; indent-tabs-mode: nil; c-basic-offset: 4; tab-width: 4; -*-
// vim: set shiftwidth=4 softtabstop=4 expandtab:
/*
 ********************************************************************
 ** NIDAS: NCAR In-situ Data Acquistion Software
 **
 ** 2007, Copyright University Corporation for Atmospheric Research
 **
 ** This program is free software; you can redistribute it and/or modify
 ** it under the terms of the GNU General Public License as published by
 ** the Free Software Foundation; either version 2 of the License, or
 ** (at your option) any later version.
 **
 ** This program is distributed in the hope that it will be useful,
 ** but WITHOUT ANY WARRANTY; without even the implied warranty of
 ** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 ** GNU General Public License for more details.
 **
 ** The LICENSE.txt file accompanying this software contains
 ** a copy of the GNU General Public License. If it is not found,
 ** write to the Free Software Foundation, Inc.,
 ** 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
 **
 ********************************************************************
*/
 
#ifndef _nidas_dynld_raf_2d_processing_h_
#define _nidas_dynld_raf_2d_processing_h_
 
#include <nidas/core/Sample.h>
 
 
namespace nidas { namespace dynld { namespace raf {
 
using namespace nidas::core;
 
class TwoD_Processing
{
 
public:
    TwoD_Processing();
    virtual ~TwoD_Processing();
 
    virtual int NumberOfDiodes() const = 0;
 
    float getResolution() const { return _resolutionMeters; }
 
    unsigned int getResolutionMicron() const { return _resolutionMicron; }
 
 
protected:
    class Particle
    {
    public:
        Particle() : height(0), width(0), area(0), edgeTouch(0), liveTime(0), dofReject(false) { } ;
        void zero() { height = width = area = liveTime = 0; edgeTouch = 0; dofReject = false; }
 
        unsigned int height;
        unsigned int width;
        unsigned int area;
        unsigned char edgeTouch;
        unsigned int liveTime;
 
        bool dofReject;
    } ;
 
    virtual void processParticleSlice(Particle& p, const unsigned char * slice);
 
    virtual void countParticle(const Particle& p, float resolutionUsec);
 
 
    virtual bool acceptThisParticle1D(const Particle& p) const;
    virtual bool acceptThisParticle2D(const Particle& p) const;
 
 
    virtual void createSamples(dsm_time_t nextTimeTag,std::list<const Sample *>&results) throw();
 
    virtual void clearData();
 
    unsigned int _numImages;
 
    long long _lastStatusTime;
 
 
    float _resolutionMeters;
 
    unsigned int _resolutionMicron;
 
    dsm_sample_id_t _1dcID, _2dcID;
 
 
    unsigned int * _size_dist_1D;
    unsigned int * _size_dist_2D;
 
    float _dead_time;
 
 
    unsigned int _totalRecords;
    unsigned int _totalParticles;
    unsigned int _rejected1D_Cntr, _rejected2D_Cntr;
    unsigned int _overLoadSliceCount;
    unsigned int _overSizeCount_2D;
    unsigned int _misAligned;
    unsigned int _suspectSlices;
    unsigned int _recordsPerSecond;
 
    unsigned int _totalPixelsShadowed;
 
    dsm_time_t _prevTime;
 
    long long _histoEndTime;
 
    float _twoDAreaRejectRatio;
 
    Particle _particle;
 
    int _nextraValues;
 
    void saveBuffer(const unsigned char * cp, const unsigned char * eod);
 
    void setupBuffer(const unsigned char ** cp, const unsigned char ** eod);
 
    unsigned char * _saveBuffer;
 
    int _savedBytes;
 
    int _savedAlloc;
 
 
private:
 
    TwoD_Processing(const TwoD_Processing&);
 
    TwoD_Processing& operator=(const TwoD_Processing&);
};
 
 
template <typename T, typename V>
inline
void
stream_histogram(T& out, V* sizedist, unsigned int nbins)
{
    bool zeros = true;
    for (unsigned int i = 0; i < nbins; ++i)
    {
        if (sizedist[i] == 0)
        {
            out << ".";
            zeros = true;
        }
        else if (zeros)
        {
            out << sizedist[i];
            zeros = false;
        }
        else
        {
            out << ",";
            out << sizedist[i];
        }
    }
}
 
}}}                     // namespace nidas namespace dynld namespace raf
 
#endif