BmnRun8cgHalfDM.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
 
    BmnRun8cgHalfDM class declaration
    Run8 coarse-grained representation, single half of each station only
*/
 
#ifndef BMNRUN8CGHALFDM_H
#define BMNRUN8CGHALFDM_H
 
#include "BmnRun8cgDM.h"
 
/*
 *  Common base for the coarse-grained half-station models. Same scheme as
 *  BmnRun8cgDM, but only one half of every station is alignable (one element
 *  per station). Top elements carry the first half of module numbers, bottom
 *  elements the second half (see BmnRun8DM). Derived classes pick the kept
 *  half via KeepBottomHalf(). Intended for data with half-station hits only.
 */
 
class BmnRun8cgHalfDM : public BmnRun8cgDM
{
  protected:
    // selects which half of every station is alignable
    virtual Bool_t KeepBottomHalf() const noexcept = 0;
 
    // only the kept half of the station survives
    virtual Int_t StationAlignables(Int_t /*station*/) const noexcept override { return 1; }
 
    virtual void Reset() noexcept override
    {
        fSystem = SIL_SYSTEM_ID;
        fStation = SIL_STATION_IDX[kFirst] - 1;   // set iterator BEFORE first station
        fModule = MOD_NUM_START;                  // recomputed per station in NextID
    }
 
    virtual Int_t NextID() noexcept override
    {
        if (-1 == fSystem)
            return -1;   // EOF marker must be negative
 
        fStation++;   // one alignable element per half-station
 
        // EOF check
        if (fStation > GEM_STATION_IDX[kLast])
            return fSystem = -1;
 
        // next system check
        if (fStation > SIL_STATION_IDX[kLast])
            fSystem = GEM_SYSTEM_ID;
 
        // top halfs allways has fModule = MOD_NUM_START
        // but bottom halfs depends on modules count
        if (KeepBottomHalf())
            fModule = MOD_NUM_START + StationModules(fStation - STA_NUM_START) / 2;
 
        return (fSystem << 24) | (fStation << 16) | (fModule << 4) | 0b0010;
    }
 
  public:
    virtual Bool_t IgnoreHit(Int_t CbmDetectorID) const noexcept override
    {
        const Int_t module{((CbmDetectorID >> 4) & 0xFFF) - MOD_NUM_START},
            station{((CbmDetectorID >> 16) & 0xFF) - STA_NUM_START}, half{StationModules(station) / 2};
        return KeepBottomHalf() ? module < half : module >= half;
    }
    virtual Bool_t IgnoresSome() const noexcept override { return kTRUE; }
};
 
#endif   // BMNRUN8CGHALFDM_H