BmnRun8DetModel.cxx

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.ru
 
    BmnRun8DetModel class implementation
*/
 
#include "BmnRun8DetModel.h"
 
/*
 * CbmStsHit DetectorID layout:
 *
 * 27:24 ( 4 bit) System ID
 * 23:16 ( 8 bit) Station number
 * 15:4  (12 bit) Sector (module) number
 *  3:1  ( 3 bit) Sensor number (not used, allways 0b001)
 *  0:0  ( 1 bit) Unknown (not used, allways 0b0)
 * --------------
 * 28 bit in total
 *
 * See CbmStsHit.h and BmnToCbmHitConverter::Exec()
 */
 
// convertion from CbmStsHit DetectorID
Int_t BmnRun8DetModel::IDfromHit(Int_t CbmDetectorID) const noexcept
{
    constexpr static const Int_t MOD_NUM_START4{MOD_NUM_START << 4}, ZERO_BIT4{0xFFFFFEF};
 
    if (GEM_SYSTEM_ID == (CbmDetectorID >> 24))
        // to skip GEM hot zone modules we need to erase the first bit
        // of module number if modules are enumerated from zero
        return ((CbmDetectorID - MOD_NUM_START4) & ZERO_BIT4) + MOD_NUM_START4;
    else
        return CbmDetectorID;
}
 
Int_t BmnRun8DetModel::Alignables() const noexcept
{
    Int_t al{0}, i;
    for (i = 0; i < SIL_STATION_CNT; i++)
        al += SIL_MODULES[i];
    for (i = 0; i < GEM_STATION_CNT; i++)
        al += GEM_MODULES[i] / 2;   // skip hot zones
    return al;
};
 
Int_t BmnRun8DetModel::MaxModulesInStation() const noexcept
{
    Int_t cnt{0}, i;
    for (i = 0; i < SIL_STATION_CNT; i++)
        if (SIL_MODULES[i] > cnt)
            cnt = SIL_MODULES[i];
    for (i = 0; i < GEM_STATION_CNT; i++)
        if (GEM_MODULES[i] / 2 > cnt)
            cnt = GEM_MODULES[i] / 2;
    return cnt;
}
 
void BmnRun8DetModel::Reset() noexcept
{
    fSystem = SIL_SYSTEM_ID;
    fStation = SIL_STATION_IDX[kFirst];
    fModule = MOD_NUM_START - 1;   // must set iterator BEFORE first module
                                   // according to base class definition
}
 
Int_t BmnRun8DetModel::NextID() noexcept
{
    // out of range check
    if (-1 == fSystem)
        return -1;   // EOF marker must be negative
 
    // next module
    if (fSystem == GEM_SYSTEM_ID)
        fModule += 2;   // skip hot zones in GEM stations
    else
        fModule++;
 
    // next station check
    if ((fModule - MOD_NUM_START)
        / (fSystem == GEM_SYSTEM_ID ? GEM_MODULES[fStation - GEM_STATION_IDX[kFirst]]
                                    : SIL_MODULES[fStation - SIL_STATION_IDX[kFirst]]))
    {
        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;
    // constant term 0b0010 is used to match CbmStsHit DetectorID encoding
}