QOptimumShapeFilter.cc

Go to the documentation of this file.

#include "QOptimumShapeFilter.hh"
#include "QVector.hh"
 
using namespace Diana;
 
QOptimumShapeFilter::QOptimumShapeFilter(const QVector& ap1, const QVector& an, int maxJitter, bool useDiff, bool debugOn, size_t step) : QOptimumFilter(ap1,an), fStep(step) 
{
    fMaxJitter = maxJitter;
    fDebugOn = debugOn;
    fUseDiff = useDiff;
 
 
    // Save AP before it gets modified by ManipulateInputs()
    fOriginalAveragePulse = NormalizeVector(ap1);
 
    // Build OF with new AP
    QError err = QERR_SUCCESS;
    if(fCheckForBuiltFilter != QERR_SUCCESS) {
        // note that BuildFilter calls ManipulateInputs()
        err = BuildFilter();
        if(err != QERR_SUCCESS) {DianaThrow(err);}
        else fCheckForBuiltFilter = QERR_SUCCESS;
    }
 
 
    //  filter original AP with OF built with modified shape to study the shape
    err = Filter(fOriginalAveragePulse);
    if(err != QERR_SUCCESS) DianaThrow(err);
 
    fOriginalAveragePulseFilteredShifted = GetFilteredShifted();
    fOriginalAveragePulseFilteredMaxPos = fOriginalAveragePulseFilteredShifted.GetMaxIndex();
    double base = fOriginalAveragePulseFilteredShifted.GetMean(fOriginalAveragePulseFilteredMaxPos*3./4);
    double gain = fOriginalAveragePulseFilteredShifted.GetMax()-fOriginalAveragePulseFilteredShifted.GetMin();
    fOriginalAveragePulseFilteredShifted -= base;
    fOriginalAveragePulseFilteredShifted /= gain;
    double max = fOriginalAveragePulseFilteredShifted.GetMax();
 
    fDeltat50 = fOriginalAveragePulseFilteredMaxPos;
    fDeltat80 = fOriginalAveragePulseFilteredMaxPos;
    fDeltat20 = fOriginalAveragePulseFilteredMaxPos;
    for(int i = fOriginalAveragePulseFilteredMaxPos; i > 0; i--) {
        if((max-fOriginalAveragePulseFilteredShifted[i])<0.8) {
            fDeltat20 = fOriginalAveragePulseFilteredMaxPos-i;
        }
        if((max-fOriginalAveragePulseFilteredShifted[i])<0.5) {
            fDeltat50 = fOriginalAveragePulseFilteredMaxPos-i;
        }
        if((max-fOriginalAveragePulseFilteredShifted[i])<0.2) {
            fDeltat80 = fOriginalAveragePulseFilteredMaxPos-i;
        }
 
    }
    fFilter1 /= gain;
    fOffset = -base/gain;
}
 
 
 
QError QOptimumShapeFilter::ManipulateInputs()
{
    // substitute the average pulse with a fast one
    QError err = QERR_SUCCESS;
    fAveragePulse1 = NormalizeVector(fAveragePulse1);
    fMaxPos = fAveragePulse1.GetMaxIndex();
    fStep = 2;
 
    QVector fastAp(fSize);
    fastAp.Initialize(0);
    for(int i = fMaxPos; i > 0; i--) {
        // sample the rising edge at step rate
        int j = fStep*(i-fMaxPos)+fMaxPos;
        if(j > 0) {  
            fastAp[i] = fAveragePulse1[j];
            // mirror to the right
            fastAp[2*fMaxPos-i] = fAveragePulse1[j];
        }
    }
    fAveragePulse1 = fastAp;
    err = this->QOptimumFilter::ManipulateInputs();
    return err;
}
 
 
QError QOptimumShapeFilter::GetChiLeft(double& chiLeft) const
{
    if(fCheckForFilteredSamples != QERR_SUCCESS) return fCheckForFilteredSamples;
    QError err = QERR_SUCCESS;
    return err;
}
 
 
QError QOptimumShapeFilter::GetChiRight(double& chiRight) const
{
    if(fCheckForFilteredSamples != QERR_SUCCESS) return fCheckForFilteredSamples;
    QError err = QERR_SUCCESS;
    return err;
}
 
QError QOptimumShapeFilter::GetRiseTime(double& riseTime) const
{
    if(fCheckForFilteredSamples != QERR_SUCCESS) return fCheckForFilteredSamples;
    QError err = QERR_SUCCESS;
    return err;
}
QError QOptimumShapeFilter::GetDecayTime(double& decayTime) const
{
    if(fCheckForFilteredSamples != QERR_SUCCESS) return fCheckForFilteredSamples;
    QError err = QERR_SUCCESS;
    return err;
}