/**
    @file DarkFunctor.cpp
    @author Brian Magill
    @creation date 1/01/2007

	$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 Function class for correcting for dark current
 
*/

#include <cmath>
#include <iostream>
#include <vector>
#include <stdexcept>
#include "SOFIE_namespace.h"
#include "GATS_Utilities.hpp"
#include "DarkFunctor.h"
#include "ChannelSignals.h"
#include "DarkSigParam.h"
#include "CalcDark.h"
#include "ApplyDark.h"

using namespace std;
using namespace GATS_Utilities;

        
MultiSignal DarkFunctor::operator() (MultiSignal const &inSignal) const
    {
    
    valarray<double> timeArray;
    vector<ChannelSignals> signalArray;
    double point_150km = inSignal.get150kmTime(); 
    bool mode          = inSignal.getModeFlag();
    Journal auditJournal;
    
    string processName = "Dark Signal Correction";
    ChannelSignals inChannel;
    ChannelSignals outChannel;
    string comments;

    DarkSigParam m_param;    
    CalcDark detDark(checkDetect, checkDiff);            
    ApplyDark applyCorrection;  
  
    inSignal.getTime(timeArray);    

    auditJournal = inSignal.getAudit();

    auditJournal.append(processName);

    for(int indx = 0; indx < inSignal.numberOfSignals(); indx++)
        try {
        
            inChannel = inSignal.getChannel(indx);        
            detDark(timeArray, inChannel, m_param, comments); 
            
            if (comments.size() > 0)
                auditJournal.append(comments);
                    
            outChannel = applyCorrection(inChannel, m_param);       
            signalArray.push_back(outChannel);   

            darkParam.push_back(m_param);    
             
            } catch (exception &ex) {
            
                cerr << "Exception thrown in DarkFunctor :" << 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;
            }
                
    MultiSignal outSignal(inSignal.getEventNumber(), point_150km, mode, timeArray, signalArray, auditJournal);
                
    return outSignal;
    };