TileRawChannelBuilderFlatFilter.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef TILERECUTILS_TILERAWCHANNELBUILDERFLATFILTER_H
#define TILERECUTILS_TILERAWCHANNELBUILDERFLATFILTER_H
 
/********************************************************************
 *
 * NAME:     TileRawChannelBuilderFlatFilter.h 
 * PACKAGE:  offline/TileCalorimeter/TileRecUtils
 *
 * AUTHOR :  K. Gellerstedt
 * CREATED:  Oct 17 2003
 *
 * PURPOSE:  Build TileRawChannels from digits using flat filter
 *
 *  Input: TileDigitsContainer
 *  Output: TileRawChannelContainer
 *  Parameters: TileRawChannelCont - Name of output container in SG
 *   PedStart
 *   PedLength
 *   PedOffset
 *   SignalStart
 *   SignalLength
 *   FilterLength
 *   FrameLength
 *   DeltaCutLo
 *   DeltaCutHi
 *   RMSCutLo
 *   RMSCutHi
 ********************************************************************/
 
// Tile includes
#include "TileRecUtils/TileRawChannelBuilder.h"
 
#include <inttypes.h>
 
class TileRawChannelBuilderFlatFilter: public TileRawChannelBuilder {
  public:
 
    // constructor
    TileRawChannelBuilderFlatFilter(const std::string& type, const std::string& name,
        const IInterface *parent);
    // destructor
    ~TileRawChannelBuilderFlatFilter();
 
    // virtual methods
    virtual StatusCode initialize();
    virtual StatusCode finalize();
 
    // Inherited from TileRawChannelBuilder
    virtual TileRawChannel * rawChannel(const TileDigits* digits, const EventContext& ctx);
 
    void flatFilter(const std::vector<uint32_t> &digits, const int gain
                    , double& amplitude, double& time) const;
 
    void flatFilter(const std::vector<float> &digits, const int gain
                    , double& amplitude, double& time) const;
 
    static const InterfaceID& interfaceID();
 
  private:
 
    bool isSignalInFrame(const std::vector<float> &digits, double deltaCut, double rmsCut,
                          int &frameMax, int &frameMin, double &frameMean, double &frameRMS) const;
    double getPedestal(const std::vector<float> &digits, int pedStart, int pedLgt) const;
 
    int getMaxAdder(const std::vector<float> &digits, int filterLength, int signalStart,
                    int signalLength, int &tMax, int &tMaxFrame, int &adderFrame) const;
 
    int getMaxSample(const std::vector<float> &digits, int signalStart, int signalLength
                    , int &tMax, int &tMaxFrame, int &sampleFrame) const;
 
    double calculatePeak(const std::vector<float> &digits, int peakPos, double ped,
                         double &position) const;
 
    double calculateFlatFilter(const std::vector<float> &digits, int filterStart
                              , int filterLength, double ped, double &position) const;
 
    double getTime(const std::vector<float> &digits, int signalStart
                  , int signalLength, double ped) const;
 
    // parameters
    int m_pedStart; // start of pedestal window
    int m_pedLength; // length of pedestal window
    int m_pedOffset; // if > 0 ped can be taken after signal as well
    int m_signalStart; // start of signal window
    int m_signalLength; // length of signal window
    int m_filterLength;
    int m_frameLength;
 
    double m_deltaCut[2]; // Cut for signalInFrame, lo/hi gain
    double m_rmsCut[2]; // RMS cut for signalInFram, lo/hi gain
 
};
 
#endif