BmnHistHgnd.h

Go to the documentation of this file.

#ifndef BMNHISTHGND_H
#define BMNHISTHGND_H
 
#include "BmnHist.h"
#include "BmnNdetAddress.h"
#include "BmnNdetGeo.h"
#include "PadInfo.h"
#include "TCanvas.h"
#include "TH1F.h"
#include "TH2F.h"
#include "THttpServer.h"
 
#include <limits>
#include <map>
#include <memory>
#include <set>
#include <string>
#include <vector>
 
class BmnHistHgnd : public BmnHist
{
  public:
    struct CanvasInfo
    {
        TCanvas* canvas{nullptr};
        std::vector<PadInfo*> pads;
 
        CanvasInfo() = default;
        CanvasInfo(TCanvas* c)
            : canvas(c)
        {}
        ~CanvasInfo() { delete canvas; }
 
        CanvasInfo(const CanvasInfo&) = delete;
        CanvasInfo& operator=(const CanvasInfo&) = delete;
        CanvasInfo(CanvasInfo&& other) noexcept
            : canvas(std::exchange(other.canvas, nullptr))
            , pads(std::move(other.pads))
        {}
        CanvasInfo& operator=(CanvasInfo&& other) noexcept
        {
            if (this != &other) {
                delete canvas;
                canvas = std::exchange(other.canvas, nullptr);
                pads = std::move(other.pads);
            }
            return *this;
        }
    };
 
    BmnHistHgnd(TString title, TString path, Int_t PeriodID, BmnSetup stp);
    virtual ~BmnHistHgnd();
 
    void CreateHistos();
    void Register(THttpServer* serv) override;
    void SetDir(TFile* outFile, TTree* recoTree) override;
    void DrawBoth() override;
    void FillFromDigi(DigiArrays* fDigiArrays) override;
    void Reset() override;
 
    BmnStatus SetRefRun(Int_t id);
    void ClearRefRun() override;
 
    void SetSelection(Int_t Trigger);
 
  private:
    void SetGeometry(const std::string& path);
    void InitializeAddresses();
    void setupAxis(TH1* hist, const std::string& xTitle, const std::string& yTitle);
 
    template<typename T, typename... Args>
    T* createH1(Args&&... args);
 
    template<typename T, typename... Args>
    T* createH2(Args&&... args);
 
    // Geometry
    std::unique_ptr<BmnNdetGeo> fGeoHandler;
 
    // All possible detector addresses
    std::set<uint32_t> fAddresses;
 
    // Address-indexed histograms
    std::map<uint32_t, TH1F*> hProfCell;      // profile across layers for each cell
    std::map<uint32_t, TH2F*> hCellAmpTime;   // time vs amplitude per cell
    std::map<uint32_t, TH2F*> hLayerAmps;     // 2d row:col for each layer
 
    // Individual histograms (for convenient access)
    TH1F* h1_ECells{nullptr};
    TH1F* h1_ProfileSlice{nullptr};
 
    TH2F* h2_grid{nullptr};
    TH2F* h2_Layer_Cell{nullptr};
 
    // All histograms for bulk operations
    std::vector<TH1*> m_allHistograms;
 
    // Canvases
    std::vector<std::unique_ptr<CanvasInfo>> m_canvases;
 
    // Selection
    bool fSelection{false};
    Int_t fSelectedTrigger{-1};
 
    UInt_t fEventCounter{1};
 
    static constexpr UInt_t kLAYRS = 2;   // number of layers
    static constexpr UInt_t kROWS = 1;    // rows per layer
    static constexpr UInt_t kCOLS = 1;    // columns per layer
 
    ClassDefOverride(BmnHistHgnd, 3);
};
 
// Template implementations must be in header
template<typename T, typename... Args>
T* BmnHistHgnd::createH1(Args&&... args)
{
    T* hist = new T(std::forward<Args>(args)...);
    m_allHistograms.push_back(hist);
    return hist;
}
 
template<typename T, typename... Args>
T* BmnHistHgnd::createH2(Args&&... args)
{
    T* hist = new T(std::forward<Args>(args)...);
    m_allHistograms.push_back(hist);
    return hist;
}
 
#endif