HitMerger2D.h

Go to the documentation of this file.

/* Copyright (C) 2024 Facility for Antiproton and Ion Research in Europe, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Dominik Smith [committer], Alexandru Bercuci */
 
#pragma once
 
#include "CbmTrdDigi.h"
#include "Math/Rotation3D.h"
#include "Math/Vector3Dfwd.h"
#include "algo/detectors/trd/DigiRec.h"
#include "algo/detectors/trd/Hit.h"
#include "algo/detectors/trd/HitFactory2D.h"
#include "algo/detectors/trd/HitFinder2DPars.h"
 
#include <cstdint>
#include <tuple>
#include <unordered_map>
#include <vector>
 
using std::vector;
class DigiRec;
 
namespace cbm::algo::trd
{

 
  class HitMerger2D {
   public:
    typedef std::pair<Hit, std::vector<DigiRec>> inputType;
    typedef std::pair<std::vector<inputType>, std::vector<inputType>> outputType;

    HitMerger2D(){};

    HitMerger2D(HitFinder2DModPar params);
 
    virtual ~HitMerger2D(){};

    outputType operator()(std::vector<inputType>& hitsRow1, std::vector<inputType>& hitsRow2);
 
    std::pair<int, HitFactory2D> ProjectDigis(std::vector<DigiRec>* cid, std::vector<DigiRec>* cjd);
    int CheckMerge(std::vector<DigiRec>* cid, std::vector<DigiRec>* cjd);
    bool MergeHits(Hit* h, int a0, HitFactory2D& hitF);
 

    void SetHitTimeOffset(int dt) { fHitTimeOff = dt; }
 
   private:
    HitMerger2D(const HitMerger2D& ref);
    const HitMerger2D& operator=(const HitMerger2D& ref);
 
    Hit MakeHit(int cId);
 
    const float fClk = CbmTrdDigi::Clk(CbmTrdDigi::eCbmTrdAsicType::kFASP);
 
    HitFinder2DModPar fParams;  
 
    void CalibrateHit(Hit* h, const double dx, const double dy, const double edx, const double edy, const double edt,
                      const double time, const double tdrift, const double eloss, const HitFactory2D& hitF);

    inline int GetPadRowCol(int address, int& c);
 
    int fHitTimeOff = 0;  
  };
 
 
}  // namespace cbm::algo::trd