BmnSiliconRaw2Digit.h

Go to the documentation of this file.

#ifndef BMNSILICONRAW2DIGIT_H
#define BMNSILICONRAW2DIGIT_H
 
#include "BmnADCDigit.h"
#include "BmnAdcProcessor.h"
#include "BmnEnums.h"
#include "BmnSiliconDigit.h"
#include "TClonesArray.h"
#include "TColor.h"
#include "TMath.h"
#include "TString.h"
#include "TTree.h"
 
#include <BmnEventHeader.h>
#include <BmnFunctionSet.h>
#include <TSystem.h>
#include <fstream>
#include <iostream>
#include <list>
#include <map>
#include <tuple>
#include <vector>
 
using std::map;
using std::string;
using std::unordered_map;
using std::vector;
 
struct BmnSiliconMapInfo
{
    uint32_t serial;
    uint16_t station;
    uint16_t module;
    uint8_t layer;   // X = 0 or X' = 1
    bool inverted;
    uint16_t channel_low;
    uint16_t channel_high;
    int32_t strip_min;
    int32_t strip_max;
};
/**
 *  upper channel to map info
 * (range map (C++20 is still unavailable, so we have to use crutches))
 */
typedef map<int32_t, BmnSiliconMapInfo*> InChanMapSil;
 
class BmnSiliconRaw2Digit : public BmnAdcProcessor
{
  public:
    BmnSiliconRaw2Digit(uint32_t period, uint32_t run, TString MapFileName, BmnSetup bmnSetup = kBMNSETUP);
    virtual ~BmnSiliconRaw2Digit();
 
    BmnStatus FillEvent(TClonesArray* adc, TClonesArray* sts);
    BmnStatus FillProfiles(TClonesArray* adc);
    BmnStatus FillNoisyChannels();
    void CreatePedestalReport();
 
  private:
    unique_ptr<BmnSiliconStationSet> fSilStationSet;
 
    unordered_map<uint32_t, InChanMapSil> fOuterMap;   // serial to inner (channel range to ) map
    map<KeySML, BmnSiliconMapInfo*> fSML2MapInfo;      // {station, module, layer} to map info
    TString fMapFileName;
 
    TCanvas* canStrip = nullptr;
    TH2F* hraw = nullptr;
    TH2F* hrms = nullptr;
    TH2F* hcorrp = nullptr;
    TH2F* hcorr = nullptr;
    TH2F* hfilter = nullptr;
    TH2F* hped = nullptr;
    TH2F* hcms = nullptr;
    TH2F* hscms = nullptr;
    TH1F* hscms1 = nullptr;
    TH1F* hscms1full = nullptr;
    TH1F* hped1 = nullptr;
    TH1F* hsig = nullptr;
    int32_t nev = -1;
    int32_t niterped;
    int32_t nchip;
    int32_t nchmin;
    int32_t npevents;
 
    int32_t zeroCnt = 0;
    int32_t allCnt = 0;
 
    Double_t thresh;
 
    Double_t thrnoise;
 
    BmnStatus ReadMapFile();
    inline int32_t MapStrip(BmnSiliconMapInfo* v, uint16_t& iCh, uint16_t& iSmpl)
    {
        if (v->inverted)
            return (v->channel_high - iCh) * GetNSamples() + iSmpl;
        else
            return (iCh - v->channel_low) * GetNSamples() + iSmpl;
    }
    void ProcessAdc(TClonesArray* silicon, Bool_t doFill);
 
    ClassDef(BmnSiliconRaw2Digit, 1);
};
 
#endif /* BMNSILICONRAW2DIGIT_H */