source: EEW/trunk/mcast2eew/EEW.h @ 3716

#ifndef __EEW_h
#define __EEW_h

#include <sys/time.h>
#include "GCBRTSPModule.h"
#include "PPickerTSPModule.h"
#include "ButhpinFilterTSPModule.h"
#include "Channel.h"
#include "DataTypes.h"
#include <queue>
#include <map>

using namespace std;

const int MAXQ = 10000;
const double PI = 3.1415926;

typedef map<TimeStamp,float> GCBRMap;
typedef deque< pair<TimeStamp,float> > GCBRQueue;
typedef deque< pair<TimeStamp,struct vd> > VDQueue;

class EEW{
 private:
  Channel chanZ,chanE,chanN;
  int samplerate;
  static unsigned int serverip;
  static bool udpsetup;

  GCBRTSPModule *gcbrZ;
  GCBRTSPModule *gcbrN;
  GCBRTSPModule *gcbrE;

  PPickerTSPModule *picker;
  ButhpinFilterTSPModule *buthfilter;

  float* buf;
  float* pickbuf;
  struct vd* vdbuf;
  GCBRQueue zq,eq,nq;
  GCBRMap   zmap,emap,nmap;
  VDQueue   vdq;
  bool trigger;
  TimeStamp trigtime;
  TimeStamp lasttimeE,lasttimeN,lasttimeZ;
  bool firsttimeE,firsttimeN,firsttimeZ;


  //TauC/PD Calculation Vars
  TimeStamp clipbegin;
  double TC;
  double Pd3;
  double Q,mag,pgv; 
  int clipped;

  void maintainGCBRData(GCBRQueue&,GCBRMap&,timeval stime,int size);
  void maintainVDQueue(VDQueue&,timeval stime,int size);
  void calcTCPD();
  double calTauC(struct vd* data,unsigned int size);
  double calPd3(struct vd* data,unsigned int size);
  bool isClipping(struct vd* data, unsigned int size);
  int tauc_pd_trigger(double tau_c,double pd);

  void generateRandomTrigger(TimeStamp stime);//for testing
  void checkgap(TimeStamp firstsamptime,TimeStamp& lastsamptime,int nsamps,bool& firsttime);
  void reset();
  float findSample(GCBRQueue& q,TimeStamp timestamp);
  void sendToServer();
  const char* gethostnametoip(const char* hostname);
 public:

  EEW(Channel Z,Channel E,Channel N);
  void init();
  void process(Channel ch,timeval stime,int*,int size);
  void clear();
  void setTriggerOn(TimeStamp timestamp);
};

#endif