AmplitudeProcessor.hpp

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



#ifndef AMPLITUDEPROCESSOR_H
#define AMPLITUDEPROCESSOR_H
    
#include "BaseLineProcessor.hpp"
#include "PeakFindProcessor.hpp"
#include "Exceptions.h"

#include <TObject.h>
namespace TrAnal 
{

class AmplitudeProcResult 
{
    public:
    double amp; 
    double uamp; 

    double baseline; 
    double ubaseline; 

    // ROOT dictionary generation
    ClassDef(AmplitudeProcResult,1);
};


template<class T>
AmplitudeProcResult ComputeAmplitude( const TrRange<T>& bline_range, const TrRange<T>& pf_range)
{ 
    // Call the primitive processors
    // There results will be deleted automatically when they go out of scope
    BaseLineProcResult bl_res 
        = ComputeBaseLine(bline_range.begin(), bline_range.end());

    PeakFindProcResult<uint16_t> pf_res 
        = FindPeak(pf_range.begin(), pf_range.end());

    if (pf_res.max < bl_res.mean) {
        throw InvalidResultException("amplitude less than baseline");
    }

    // No exceptions were thrown, create the result and fill it with values
    AmplitudeProcResult res;
    res.amp = pf_res.max - bl_res.mean;
    res.uamp = ::sqrt(::pow(bl_res.stdev,2.0) );

    res.baseline = bl_res.mean;
    res.ubaseline = bl_res.stdev;

    return res;
}

} // end namespace
#endif