BmnGemRaw2Digit.h

Go to the documentation of this file.

#ifndef BMNGEMRAW2DIGIT_H
#define BMNGEMRAW2DIGIT_H
 
#include <map>
#include <vector>
// ROOT
#include "BmnAdcProcessor.h"
#include "TClonesArray.h"
#include "TString.h"
 
#define N_CH_BUF 4096
 
using std::map;
using std::vector;
 
class BmnGemRaw2Digit : public BmnAdcProcessor
{
 
  public:
    BmnGemRaw2Digit(Int_t period, Int_t run, TString mapFileName, BmnSetup bmnSetup = kBMNSETUP);
    ~BmnGemRaw2Digit();
 
    BmnStatus FillEvent(TClonesArray* adc, TClonesArray* gem);
    BmnStatus FillProfiles(TClonesArray* adc);
    BmnStatus FillNoisyChannels();
 
    struct BmnGemMap
    {
        Int_t strip;
        Int_t lay;   // strip type: 0 - x, 1 - y
        Int_t mod;   // hot zones: 1 - inner zone, 0 - outer zone
 
        BmnGemMap(Int_t s, Int_t l, Int_t m)
            : strip(s)
            , lay(l)
            , mod(m)
        {}
 
        BmnGemMap()
            : strip(-1)
            , lay(-1)
            , mod(-1)
        {}
    };
    typedef std::array<BmnGemMap, N_CH_BUF> StripMapArray;
 
    struct GemMapLine
    {
        Int_t Serial;
        Int_t Ch_lo;
        Int_t Ch_hi;
        Int_t GEM_id;
        Int_t Side;
        Int_t Type;
        Int_t Station;
        Int_t Module;
        Int_t Zone;
        int32_t Channel_shift;
        StripMapArray* Map_arr_ptr;
 
        void Print()
        {
            cout << std::hex << Serial << std::dec << " " << Ch_lo << " " << Ch_hi << " " << GEM_id << " " << Side
                 << " " << Type << " " << Station << " " << Module << " " << Zone << endl;
        }
 
        GemMapLine()
            : Serial(-1)
            , Ch_lo(-1)
            , Ch_hi(-1)
            , GEM_id(-1)
            , Side(-1)
            , Type(-1)
            , Station(-1)
            , Module(-1)
            , Zone(-1)
            , Channel_shift(0)
            , Map_arr_ptr(nullptr)
        {}
    };
    typedef map<int32_t, GemMapLine*> InChanMap;
 
  private:
    unique_ptr<BmnGemStripStationSet> fGemStationSet;
    // starting thresholds, number of iterations
    Int_t niterped = 3;
    Float_t thrped = 35;
    Float_t thrpedcsc = 80;
 
    unordered_map<uint32_t, InChanMap> fOuterMap;   // serial to inner (channel range to local strip map) map
 
    StripMapArray fSmall;
    StripMapArray fMid;
    vector<StripMapArray> fBigHot;
    vector<StripMapArray> fBig;
    TString fMapFileName;
 
    typedef map<uint8_t, int32_t> MapLayer2MaxStrip;
    map<KeySM, MapLayer2MaxStrip> fSM2MapMaxStrip;   // {station, module} to max strip
    inline void MapStrip(GemMapLine* gemM, UInt_t ch, Int_t iSmpl, Int_t& lay, Int_t& strip);
    /**
     * Uses raw ADC data to either fill profile histograms or produce output digits
     * @param silicon raw ADC signals
     * @param doFill true - fill profile histograms, false - produce digits
     */
    void ProcessAdc(TClonesArray* silicon, Bool_t doFill);
    BmnStatus ReadMap(TString parName, StripMapArray& m, Int_t lay, Int_t mod);
    BmnStatus ReadLocalMap(TString parName, StripMapArray& m, Int_t lay, Int_t mod);
    BmnStatus ReadGlobalMap(TString FileName);
    void InitStripRanges(unique_ptr<BmnGemStripStationSet>& ss);
 
    ClassDef(BmnGemRaw2Digit, 1);
};
 
#endif /* BMNGEMRAW2DIGIT_H */