source: EEW/trunk/mcast2eew/GCBRTSPModule.cc @ 3772

#include "GCBRTSPModule.h"
#include "eewPropertiesST.h"
#include "Duration.h"
#include "EEW.h"
//#define FILEDUMP 0

GCBRTSPModule::GCBRTSPModule(){}

GCBRTSPModule::GCBRTSPModule(Channel chn): Module(chn){
  //Initialize per-channel datastructures here//
    LTA = 0;
    tflag = false;
    gain = chn.gain;
    
//     eewPropertiesST* prop = eewPropertiesST::getInstance();

//     if(string(channel.channel) == string(string(prop->getChannelType())+"Z")){
//       gain = prop->getgainZ();
//     }
//     else if(string(channel.channel) == string(string(prop->getChannelType())+"E")){
//       gain = prop->getgainE();
//     }
//     else if(string(channel.channel) == string(string(prop->getChannelType())+"N")){
//       gain = prop->getgainN();
//     }
    cout <<"Gain for "<<chn<<" is "<<gain<<endl;

#ifdef WRITE2SAC
    int samplerate = chn.samprate;//eewPropertiesST::getInstance()->getSampleRate();
    sacwrt.open(channel,"acc",1.0/samplerate);
#endif

}

GCBRTSPModule::~GCBRTSPModule(){
  //Release per-channel resources here//
#ifdef WRITE2SAC
    sacwrt.close();
#endif
}


void GCBRTSPModule::load(string config_filename) throw(Error){
    //initialize the module
  cout <<"Loading GCBR Module"<<endl;
}

void GCBRTSPModule::unload(){
  //release the resources
  cout <<"Unloading GCBR Module"<<endl;
}


void GCBRTSPModule::process(timeval timestamp,
                           unsigned int in_nsamp,
                           int* in,
                           unsigned int& out_nsamp,
                           float* out){

    //Write code to process the time series data//
    //cout <<"["<<channel<<"] "<<"Start time "<<TimeStamp(timestamp)<<" Delay :"<<(double)(TimeStamp::current_time() - TimeStamp(timestamp))<<" seconds, samples="<<in_nsamp<<endl;


#ifdef WRITE2SAC
    sacwrt.write(in,in_nsamp);
#endif

    double sum = 0.0;
    double avg = 0.0;
    float multfact = (float)eewPropertiesST::getInstance()->getSampleMultFactor();

    for(int i=0;i<in_nsamp;i++){
        in[i] = in[i]*multfact;
        sum += in[i];
    }

    avg = sum/((in_nsamp)*60.0);    
    ltalist.push_back(avg);
    LTA += avg;
    if(ltalist.size() > 60){
      LTA -= ltalist.front();
      ltalist.pop_front();
    }

    unsigned long tsec = timestamp.tv_sec;
    //cout <<channel<<" gain is "<<gain<<endl;
    for(int j=0;j<in_nsamp;j++){
      out[j] = (float)((float)(in[j] - LTA)/gain); 
        tsec++;
    }
    out_nsamp = in_nsamp;

    //   for(int i=0;i<in_nsamp;i++){
    //  cout <<out[i]<<",";
    // }
    //cout <<endl;
}



void GCBRTSPModule::restart(timeval timestamp){
    //Write code to handle gaps in time series data//
  //    cout <<channel<<" GAP Found."<<endl;

    LTA = 0;
    ltalist.clear();

#ifdef WRITE2SAC
    sacwrt.close();
    sacwrt.reopen();
#endif
}