BaseLineProcessor.hpp

//  BaseLineProcessor.hpp
//
//  Author : Jeromy Tompkins
//  Date   : 8/14/2013

#ifndef BASELINEPROCESSOR_H
#define BASELINEPROCESSOR_H
    
#include "Trace.hpp"
#include <cmath>
#include "Exceptions.h"
#include <TObject.h>

namespace TrAnal 
{

class BaseLineProcResult 
{
    public:
    double mean; 
    double stdev; 

    // ROOT dict generation
    ClassDef(BaseLineProcResult,1);
};



template<class T>
double ComputeMean(const TrIterator<T>& begin, const TrIterator<T>& end)
{
   
    double sum=0;
    double n=0;
    TrIterator<T> it=begin;
    while (it<end) 
    {
        // access element
        sum += *it;
        ++n;

        // increment iterator
        ++it;
    }
    
    return sum/n;
}


template<class T>
double ComputeStDev(const TrIterator<T>& begin, const TrIterator<T>& end, double mean)
{

    double stdev2 =0;
    double tmpdiff = 0;

    double n=0;

    TrIterator<T> it=begin;
    while ( it<end ) 
    {
        tmpdiff = (*it-mean);
        stdev2 += ::pow(tmpdiff,2.0); 
        ++n;

        ++it;
    }

    return ::sqrt(stdev2 / (n - 1.0));
}


template<class T>
BaseLineProcResult ComputeBaseLine(const TrIterator<T>& begin, const TrIterator<T>& end)
{
    double mean = ComputeMean<T>(begin, end);
    double stdev = ComputeStDev<T>(begin, end, mean);

    BaseLineProcResult res;
    res.mean = mean;
    res.stdev = stdev;
    return res;
}

} // end namespace
#endif