/**
    @class CalcDarkValue
    @author Brian Magill
    @creation date 4/4/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 Calculates the dark signal for a given band
*/


#include <stdexcept>
#include <cmath>
#include <cassert>
#include "CalcDarkValue.h"
#include "ValReversi.hpp"

using namespace std;

double CalcDarkValue::operator()(valarray<double> const &signal) const
{
return signal.min();
}

double CalcDarkValue::timeAtMin(valarray<double> const &signal, 
                                  valarray<double> const &time) const
{
double minValue = this->operator()(signal);
unsigned long i;
double timeValue = 0.0;
assert(signal.size() == time.size());

for (i = 0; i < signal.size(); i++)
    {
    if( signal[i] <= minValue)
        {
        timeValue = time[i];
        break;
        }
    
    }
    
return timeValue;
}