BmnRun8cgDM.h

Go to the documentation of this file.

/*
    BM@N alignment routine
    BM@N experiment at NICA complex, JINR, 2025
 
    Department: Math & Soft Group of HEP lab
    Author:     Igor Polev, polev@jinr.int
 
    BmnRun8cgDM class declaration
    Run8 coarce-grained representation of detector elements via encoding
*/
 
#ifndef BMNRUN8CGDM_H
#define BMNRUN8CGDM_H
 
#include "BmnRun8DM.h"
 
/*
 *  Detector elements enumeration scheme
 *
 *  SystemID  - fixed constatns defined in BmnDetectorList.h
 *  StationID - sequence across all systems starting from 0 or 1 (currently 1)
 *  ModuleID  - separate sequence for each station starting from 0 or 1 (currently 1)
 */
 
class BmnRun8cgDM : public BmnRun8DM
{
  private:
    Int_t fStationStep[SIL_STATION_CNT + GEM_STATION_CNT];
 
  protected:
    // one element per station half
    virtual Int_t StationAlignables(Int_t /*station*/) const noexcept override { return 2; }
 
    virtual void Reset() noexcept override
    {
        fSystem = SIL_SYSTEM_ID;
        fStation = SIL_STATION_IDX[kFirst];
        fModule = MOD_NUM_START - fStationStep[fStation - STA_NUM_START];
    }
 
    virtual Int_t NextID() noexcept override
    {
        // out of range check
        if (-1 == fSystem)
            return -1;   // EOF marker must be negative
 
        // next module
        fModule += fStationStep[fStation - STA_NUM_START];
 
        // next station check
        if ((fModule - MOD_NUM_START) / StationModules(fStation - STA_NUM_START)) {
            fModule = MOD_NUM_START;
            fStation++;
            // remember: stations has common enumeration across all systems
            // currently only two systems are considered - silicon, then GEM
 
            // EOF check
            if (fStation > GEM_STATION_IDX[kLast])
                return fSystem = -1;
 
            // next system check
            if (fStation > SIL_STATION_IDX[kLast])
                fSystem = GEM_SYSTEM_ID;
        }
 
        return (fSystem << 24) | (fStation << 16) | (fModule << 4) | 0b0010;
    }
 
  public:
    // convertion from CbmStsHit DetectorID
    // for detailed CbmStsHit DetectorID layout see BmnRun8DM
    virtual Int_t IDfromHit(Int_t CbmDetectorID) const noexcept override
    {
        Int_t module{((CbmDetectorID >> 4) & 0xFFF) - MOD_NUM_START},
            station{((CbmDetectorID >> 16) & 0xFF) - STA_NUM_START};
 
        module /= fStationStep[station];
        module *= fStationStep[station];
        module += MOD_NUM_START;
 
        return (CbmDetectorID & 0xFFF000F) | (module << 4);
    }
 
    BmnRun8cgDM()
    {
        fDescription = "Run8 GEM and FSD hard parts (half stations): 22 elements";
 
        // StationModules() is not virtual, so calling it from the constructor is safe
        for (Int_t i = 0; i < SIL_STATION_CNT + GEM_STATION_CNT; i++)
            fStationStep[i] = StationModules(i) / 2;
    }
    virtual ~BmnRun8cgDM() override;   // out-of-line vtable anchor
};
 
#endif   // BMNRUN8CGDM_H