TwoD_USB.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_usb_h_
#define _nidas_dynld_raf_2d_usb_h_

#include <nidas/core/DSMSensor.h>
#include <nidas/core/DerivedDataClient.h>

#include <nidas/util/EndianConverter.h>
#include <nidas/util/InvalidParameterException.h>

#include <nidas/linux/usbtwod/usbtwod.h>

namespace nidas { namespace dynld { namespace raf {

using namespace nidas::core;

class TwoD_USB : public DSMSensor, public DerivedDataClient
{

public:
    TwoD_USB();
    ~TwoD_USB();

    IODevice *buildIODevice() throw(nidas::util::IOException);

    SampleScanner *buildSampleScanner()
        throw(nidas::util::InvalidParameterException);

    void open(int flags)
        throw(nidas::util::IOException,nidas::util::InvalidParameterException);

    void close() throw(nidas::util::IOException);

    int getTASRate() const { return _tasRate; }

    void setTASRate(int val) { _tasRate = val; }

    float getResolution() const { return _resolutionMeters; }

    unsigned int getResolutionMicron() const { return _resolutionMicron; }

    virtual int NumberOfDiodes() const = 0;

    void init() throw(nidas::util::InvalidParameterException);

    virtual void
    derivedDataNotify(const nidas::core:: DerivedDataReader * s)
        throw();

    void printStatus(std::ostream& ostr) throw();

   // virtual int TASToTap2D(Tap2D * t2d, float tas);
    virtual int TASToTap2D(void * t2d, float tas);

    virtual float
    Tap2DToTAS(const Tap2D * t2d) const;
    virtual float
    Tap2DToTAS(const Tap2Dv1 * t2d) const;


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;
    } ;

    // Probe produces Big Endian.
    static const nidas::util::EndianConverter * bigEndian;

    // Tap2D value sent back from driver has little endian ntap value
    static const nidas::util::EndianConverter * littleEndian;

    virtual void init_parameters()
        throw(nidas::util::InvalidParameterException);

    virtual void sendTrueAirspeed(float tas) throw(nidas::util::IOException);

    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();

    int _tasRate;

    unsigned int _numImages;

    long long _lastStatusTime;


    float _resolutionMeters;

    unsigned int _resolutionMicron;


    dsm_sample_id_t _sorID;

    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 _tasOutOfRange;
    unsigned int _misAligned;
    unsigned int _suspectSlices;
    unsigned int _recordsPerSecond;

    unsigned int _totalPixelsShadowed;

    dsm_time_t _prevTime;

    long long _histoEndTime;

    float _twoDAreaRejectRatio;

    Particle _particle;

    float _trueAirSpeed;

    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;

    static const float DefaultTrueAirspeed;

private:

    TwoD_USB(const TwoD_USB&);

    TwoD_USB& operator=(const TwoD_USB&);
};


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