/**
     @file NonLinFunctor.cpp
     @author Brian Magill
     @date 4/08/2008

	$Date:$
    $Revision:$

    @copyright
     (©) Copyright 2006 by GATS Inc.
     11864 Canon Blvd., Suite 101, Newport News, VA 23606

     All Rights Reserved. No part of this software or publication may be
     reproduced, stored in a retrieval system, or transmitted, in any form
     or by any means, electronic, mechanical, photocopying, recording, or
     otherwise without the prior written permission of GATS Inc.
        
    @brief Corrects the nonlinearity in the signals
     
*/

#include <cmath>
#include <iostream>
#include <vector>
#include <stdexcept>
#include "SOFIE_namespace.h"
#include "GATS_Utilities.hpp"
#include "NonLinFunctor.h"
#include "ChannelSignals.h"
#include "RemoveNonLinAlgo.h"

using namespace std;
using namespace GATS_Utilities;

        
MultiSignal NonLinFunctor::operator() (MultiSignal const &inSignal) const
    {
    
    valarray<double> timeArray;
    vector<ChannelSignals> signalArray;

    Journal auditJournal;
      
    string processName = "Nonlinearity Correction";
    ChannelSignals inChannel;
    ChannelSignals outChannel;
    string comments;

    RemoveNonLinAlgo correction(labGain);
             
    inSignal.getTime(timeArray);    

    auditJournal = inSignal.getAudit();

    for(int indx = 0; indx < inSignal.numberOfSignals(); indx++)
        try {
                
            inChannel = inSignal.getChannel(indx);        
            outChannel = correction(timeArray, inChannel, 
                                    setup.getStrong(indx), setup.getWeak(indx) );       
            signalArray.push_back(outChannel);   

            } catch (exception &ex) {
            
            cerr << "Exception thrown in NonLinFunctor :" << endl;
            cerr << "Channel Gas ID: " << inChannel.getID() << endl;
            cerr << "Exception Message: " << ex.what() << endl;
            cerr << "Additional comments: " << comments << endl;
            cerr << "Current Audit Trail: ";
            inSignal.getAudit().print(cerr);
            cerr << endl;
            
            throw;
            }

    auditJournal.append(processName);
    
    MultiSignal outSignal(inSignal.getEventNumber(), inSignal.get150kmTime(), 
                            inSignal.getModeFlag(), timeArray, signalArray, auditJournal) ;
                
    return outSignal;
    };