TileRawChannelBuilderOpt2Filter.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef TILERECUTILS_TILERAWCHANNELBUILDEROPT2FILTER_H
#define TILERECUTILS_TILERAWCHANNELBUILDEROPT2FILTER_H
 
//
//     Based on the code of Ximo Poveda@cern.ch. June 2007
//     Andrei.Artamonov@cern.ch, July 2008
//
//     TileRawChannelBuilderOpt2Filter.h
//
//     implementation of the Optimal Filtering based on Lagrange multipliers
//       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 <vector>
#include <string>
 
class TileRawChannelBuilderOpt2Filter: public TileRawChannelBuilder {
  public:
 
    TileRawChannelBuilderOpt2Filter(const std::string& type, const std::string& name,
        const IInterface *parent); 
    ~TileRawChannelBuilderOpt2Filter(); 
 
    // virtual methods
    virtual StatusCode initialize(); 
    //virtual StatusCode execute();
    virtual StatusCode finalize(); 
 
    // Inherited from TileRawChannelBuilder
    virtual TileRawChannel* rawChannel(const TileDigits* digits, const EventContext& ctx);
 
    static const InterfaceID& interfaceID();
 
  private:
 
    ToolHandle<ITileCondToolOfc> m_tileCondToolOfc{this,
        "TileCondToolOfc", "TileCondToolOfc", "Tile OFC tool"};
 
    ToolHandle<TileCondToolNoiseSample> m_tileToolNoiseSample{this,
        "TileCondToolNoiseSample", "TileCondToolNoiseSample", "Tile noise sample tool"};
 
    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);
    int iterate(int ros, int drawer, int channel, int gain, double &pedestal, double &amplitude, double &time, double &chi2, const EventContext &ctx);
    double compute(int ros, int drawer, int channel, int gain, double &pedestal, double &amplitude, double &time, double& phase, const EventContext &ctx);
 
    void ofc2int(int nDigits, double* w_off, short* w_int, short& scale); // convert weights to dsp short int format
 
    int m_maxIterations; 
    int m_pedestalMode;  
    bool m_confTB;       
    double m_timeForConvergence; 
    bool m_of2;    
    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; 
 
 
    std::vector<float> m_digits;
 
    int m_noiseThresholdHG;
    int m_noiseThresholdLG;
};
 
#endif