TileRawChannelBuilderWienerFilter.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef TILERECUTILS_TILERAWCHANNELBUILDERWIENERFILTER_H
#define TILERECUTILS_TILERAWCHANNELBUILDERWIENERFILTER_H
 
//
//     AUTHOR :  G. Goncalves <ginaciog@cern.ch>
//     CREATED:  May 2019
//
//     TileRawChannelBuilderWienerFilter.h
//
//     implementation of the Wiener Filter for
//          energy/time reconstruction in TileCal
//
 
// Tile includes
#include "TileRecUtils/TileRawChannelBuilder.h"
#include "TileConditions/ITileCondToolOfc.h"
#include "TileConditions/TileCondToolTiming.h"
#include "TileConditions/TileCondToolNoiseSample.h"
#include "TileConditions/TileWienerFilterWeights.h"
 
// Atlas includes
#include "StoreGate/ReadCondHandleKey.h"
#include "LumiBlockData/BunchCrossingCondData.h"
 
#include <vector>
#include <string>
 
class TileRawChannelBuilderWienerFilter: public TileRawChannelBuilder {
  public:
 
    TileRawChannelBuilderWienerFilter(const std::string& type, const std::string& name,
        const IInterface *parent); 
    ~TileRawChannelBuilderWienerFilter(); 
 
    // virtual methods
    virtual StatusCode initialize() override; 
    virtual StatusCode finalize() override; 
 
    // Inherited from TileRawChannelBuilder
    virtual TileRawChannel* rawChannel(const TileDigits* digits, const EventContext& ctx) override;
 
    static const InterfaceID& interfaceID();
 
  private:
 
    ToolHandle<TileCondToolNoiseSample> m_tileToolNoiseSample{this,
        "TileCondToolNoiseSample", "TileCondToolNoiseSample", "Tile noise sample tool"};
 
    SG::ReadCondHandleKey<BunchCrossingCondData> m_bunchCrossingKey{this,
        "BunchCrossingCondDataKey", "BunchCrossingData" ,"SG Key of BunchCrossing CDO"}; 
 
    double filter(int ros, int drawer, int channel, int &gain, double &pedestal, double &amplitude, double &time, const EventContext &ctx);
    int findMaxDigitPosition();  
    float getPedestal(int ros, int drawer, int channel, int gain, const EventContext &ctx);
    double compute(int ros, int drawer, int channel, int gain, double &pedestal, double &amplitude, double &time, double& phase);
    int getBCIDIndex(const EventContext &ctx);
 
    int m_maxIterations; 
    int m_pedestalMode;  
    bool m_confTB;       
    double m_timeForConvergence; 
    bool m_isMC;    
    bool m_minus1Iter;   
    bool m_correctAmplitude; 
    bool m_correctTimeNI; 
 
    bool m_bestPhase; // if true, use best phase from COOL DB in "fixed phase" mode (i.e., no iterations)
    bool m_emulateDsp; // if true, emulate DSP reconstruction algorithm
    int m_nSignal; 
    int m_nNegative;  
    int m_nCenter; 
    int m_nConst;  
 
    int m_nSamples; 
    int m_t0SamplePosition;  
    double m_maxTime; 
    double m_minTime; 
 
    const TileWienerFilterWeightsStruct *m_weights; 
 
    std::vector<float> m_digits;
};
 
#endif