PeakFindProcessor.hpp

// PeakFindProcessor.h
//
// Author : Jeromy Tompkins
// Date   : 8/14/2013


#ifndef PEAKFINDPROCESSOR_H
#define PEAKFINDPROCESSOR_H
    
#include "Trace.hpp"
#include "AlgoIterator.hpp"
#include "TObject.h"

namespace TrAnal 
{

template<class T>
class PeakFindProcResult 
{
    public:
    double max; 
    double umax; 
    
    TrIterator<T> index;  

    // ROOT dictionary generation 
    ClassDef(PeakFindProcResult,1);

}; // end of PeakFindProcResult class


template<class T>
PeakFindProcResult<T> FindPeak(const TrIterator<T>& begin, const TrIterator<T>& end)
{
    
    double max=0;
    double current=0;
    TrIterator<T> index=begin; 

    TrIterator<T> it=begin;
    while ( it<end ) {
        current = *it;
        if (current>max) {
            max = current;
            index = it;
        }
        ++it;
    }

    // We successfully completed the algorithm...store the result
    // and pass it to the caller
    PeakFindProcResult<T> res;
    res.max = max;
    res.umax = 0;
    res.index = index;

    return res;
} // end of FindPeak


template<class T>
PeakFindProcResult<T> FindPeak(AlgoIterator<T>& begin, const TrIterator<T>& end)
{
    
    double max=0;
    double current=0;
    TrIterator<T> index=begin.min_extent(); 

    AlgoIterator<T>& it=begin;
    while ( it<end ) {
        current = *it;
        if (current>max) {
            max = current;
            index = it.min_extent();
        }
        ++it;
    }

    // We successfully completed the algorithm...store the result
    // and pass it to the caller
    PeakFindProcResult<T> res;
    res.max = max;
    res.umax = 0;
    res.index = index;

    return res;
} // end of FindPeak

} // end namespace
#endif