BmnRun8cgDetModel.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 "BmnRun8cgDetModel.h"
 
BmnRun8cgDetModel::BmnRun8cgDetModel()
{
    for (Int_t i = 0; i < SIL_STATION_CNT + GEM_STATION_CNT; i++)
        if (i < SIL_STATION_CNT)
            fStationStep[i] = SIL_MODULES[i] / 2;
        else
            fStationStep[i] = GEM_MODULES[i - SIL_STATION_CNT] / 2;
}
 
// convertion from CbmStsHit DetectorID
// for detailed CbmStsHit DetectorID layout see BmnRun8DetModel
Int_t BmnRun8cgDetModel::IDfromHit(Int_t CbmDetectorID) const noexcept
{
    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);
}
 
Int_t BmnRun8cgDetModel::MaxModulesInStation() const noexcept
{
    return 2;
}
 
Int_t BmnRun8cgDetModel::Alignables() const noexcept
{
    return 2 * (SIL_STATION_CNT + GEM_STATION_CNT);
}
 
void BmnRun8cgDetModel::Reset() noexcept
{
    fSystem = SIL_SYSTEM_ID;
    fStation = SIL_STATION_IDX[kFirst];
    fModule = MOD_NUM_START - fStationStep[fStation - STA_NUM_START];
}
 
Int_t BmnRun8cgDetModel::NextID() noexcept
{
    // 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)
        / (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;
}