/*
 * To change this license header, choose License Headers in Project Properties.
 * To change this template file, choose Tools | Templates
 * and open the template in the editor.
 */

/* 
 * File:   SpectrumProcess.cpp
 * Author: hacene
 * 
 * Created on December 31, 2020, 3:42 PM
 */


#include <stdlib.h>
#include <string>

#include <boost/lexical_cast.hpp>

#include "DicError.hh"
#include "AMDA_exception.hh"

#include "Operation.hh"
#include "ParameterManager.hh"
#include "ParamMgr.hh"
#include "ParamTable.hh"
#include "ParamInfo.hh"
#include "ParamMgr.hh"
#include "ParameterCreatorFromExpression.hh"

#include "SpectrumProcess.hh"


using namespace std;
using namespace boost;
using namespace log4cxx;

namespace AMDA {
    namespace Parameters {

        SpectrumProcess::SpectrumProcess(Parameter &parameter) : SingleParamProcess_CRTP(parameter) {

            _type = SpectrumCreator::SpectrumType::FOURIER;
            _paramInfoId = parameter.getId();
            parameter.setInfoId(_paramInfoId);
        }

        SpectrumProcess::SpectrumProcess(const SpectrumProcess& pProcess, Parameter &parameter) : SingleParamProcess_CRTP(pProcess, parameter) {
            _type = SpectrumCreator::SpectrumType::FOURIER;
             _paramInfoId = parameter.getId();
            parameter.setInfoId(_paramInfoId);
        }

        SpectrumProcess::~SpectrumProcess() {
        }

        TimeStamp SpectrumProcess::init() {


            if (_attributList.size() < 1) {
                BOOST_THROW_EXCEPTION(AMDA::AMDA_exception() << AMDA::errno_code(AMDA_PROCESS_ERR) << AMDA::ex_msg(std::string("ProcessSpectrum::parse require 1 attribute'")));
            }
            _windowtime = atof(_attributList[0].c_str());

            _sampling = "center";
            if (_attributList.size() > 1)
                _sampling = _attributList[1];

            /*****/


            TimeStamp time = _parameterInput->init(this, _timeIntervalList);

            _paramInput = _parameterInput->getParamData(this).get();


            /**  in case we need to sample 
             ParameterSPtr sampledParamInput = _parameterInput->getParameterManager().getSampledParameter(_parameterInput->getId(), "classic", _parameterInput->getGapThreshold(), false);
                
             sampledParamInput->openConnection(this);
                
             time = sampledParamInput->init(this, _timeIntervalList);
                
             ParamData* sampledParamData = sampledParamInput->getParamData(this).get();
                
             SpectrumCreator lCreator(*this, _timeIntervalList, *sampledParamData, _type, _windowtime, _sampling);
             _operation = lCreator.getOperation();
             */

            SpectrumCreator lCreator(*this, _timeIntervalList, *_paramInput, _type, _windowtime, _sampling);
            _operation = lCreator.getOperation();

            _paramData = ParamDataSPtr(_operation->getParamOutput());

            _paramData->setMinSampling(_paramInput->getMinSampling());

           
            int size_ = std::ceil(_windowtime/_paramInput->getMinSampling()) +1;
             std::vector<double> frequencies(size_); 
             for(int i =0; i<size_; i++ )
                        frequencies[i] = (double) 1./_paramInput->getMinSampling()/(double) (i+1);

            AMDA::Info::ParamInfoSPtr paramInfo = AMDA::Info::ParamMgr::getInstance()->getParamInfoFromId(_paramInfoId);

            boost::shared_ptr<AMDA::Info::ParamTable> table(new AMDA::Info::EmbeddedCenterTable(
                    _paramInfoId.c_str(), 0, frequencies,  false));

            table->addTableParam(AMDA::Info::EmbeddedCenterTable::_paramKey,  _paramInfoId );
            table->setName("Frequencies");
            table->setUnits("Hz");
            table->setIsVariable(false);
            table->setIsFullVariable(false);

            paramInfo->addTable(0, table);

            return time;
        }

    }
}