TGCCableInASD.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2026 CERN for the benefit of the ATLAS collaboration
*/
 
#include "MuonTGC_Cabling/TGCCableInASD.h"
 
#include "MuonTGC_Cabling/TGCChannelASDIn.h"
#include "MuonTGC_Cabling/TGCChannelASDOut.h"
#include "MuonTGC_Cabling/TGCDatabaseASDToPP.h"
 
namespace MuonTGC_Cabling {
 
// Constructor & Destructor
TGCCableInASD::TGCCableInASD(const std::string& filename)
    : TGCCable(TGCCable::InASD), m_database{{{nullptr}}} {
    m_database.at(+TGCId::RegionType::Endcap).at(+TGCId::ModuleType::WD) =
        std::make_unique<TGCDatabaseASDToPP>(filename, "EWD");
    m_database.at(+TGCId::RegionType::Endcap).at(+TGCId::ModuleType::WT) =
        std::make_unique<TGCDatabaseASDToPP>(filename, "EWT");
    m_database.at(+TGCId::RegionType::Endcap).at(+TGCId::ModuleType::SD) =
        std::make_unique<TGCDatabaseASDToPP>(filename, "ESD");
    m_database.at(+TGCId::RegionType::Endcap).at(+TGCId::ModuleType::ST) =
        std::make_unique<TGCDatabaseASDToPP>(filename, "EST");
    m_database.at(+TGCId::RegionType::Endcap).at(+TGCId::ModuleType::WI) =
        std::make_unique<TGCDatabaseASDToPP>(filename, "EWI");
    m_database.at(+TGCId::RegionType::Endcap).at(+TGCId::ModuleType::SI) =
        std::make_unique<TGCDatabaseASDToPP>(filename, "ESI");
    m_database.at(+TGCId::RegionType::Forward).at(+TGCId::ModuleType::WD) =
        std::make_unique<TGCDatabaseASDToPP>(filename, "FWD");
    m_database.at(+TGCId::RegionType::Forward).at(+TGCId::ModuleType::WT) =
        std::make_unique<TGCDatabaseASDToPP>(filename, "FWT");
    m_database.at(+TGCId::RegionType::Forward).at(+TGCId::ModuleType::SD) =
        std::make_unique<TGCDatabaseASDToPP>(filename, "FSD");
    m_database.at(+TGCId::RegionType::Forward).at(+TGCId::ModuleType::ST) =
        std::make_unique<TGCDatabaseASDToPP>(filename, "FST");
    m_database.at(+TGCId::RegionType::Forward).at(+TGCId::ModuleType::WI) =
        std::make_unique<TGCDatabaseASDToPP>(filename, "FWI");
    m_database.at(+TGCId::RegionType::Forward).at(+TGCId::ModuleType::SI) =
        std::make_unique<TGCDatabaseASDToPP>(filename, "FSI");
}
 
TGCCableInASD::~TGCCableInASD() = default;
 
std::unique_ptr<TGCChannelId> TGCCableInASD::getChannel(
    const TGCChannelId& channelId, bool orChannel) const {
    if (channelId.getChannelIdType() == TGCChannelId::ChannelIdType::ASDIn) {
        return getChannelOut(channelId, orChannel);
    }
    if (channelId.getChannelIdType() == TGCChannelId::ChannelIdType::ASDOut) {
        return getChannelIn(channelId, orChannel);
    }
 
    return nullptr;
}
 
std::unique_ptr<TGCChannelId> TGCCableInASD::getChannelIn(
    const TGCChannelId& asdout, bool orChannel) const {
    if (orChannel || !asdout.isValid()) {
        return nullptr;
    }
 
    TGCDatabase* databaseP = m_database.at(+asdout.getRegionType()).at(+asdout.getModuleType()).get();
 
    // sector ASDIn [1..48, 1..24], ASDOut [0..47, 0..23]
    int sector;
    if (asdout.isEndcap() && !asdout.isInner()) {
        sector = asdout.getSector() - 1;
        if (sector <= 0) {
            sector += TGCId::NUM_ENDCAP_SECTOR;
        }
    } else {
        sector = asdout.getSector();
        if (sector <= 0) {
            sector += TGCId::NUM_FORWARD_SECTOR;
        }
    }
 
    // chamber ASDIn [1(F),1,2,3,4,5(E)], ASDOut [0(F),4,3,2,1,0(E)]
    int chamber;
    if (asdout.isEndcap() && !asdout.isInner()) {
        chamber = 5 - asdout.getChamber();
    } else {
        chamber = asdout.getChamber() + 1;
    }
 
    int channel = -1;
    // channel ASDIn [1..32(S),1..n(W chamber)], ASDOut [0..31(S),n..0(W
    // sector)]
    if (asdout.isWire()) {
        // Endcap Triplet chamberId start from 1 in ASDOut
        int dbChamber = asdout.getChamber();
        if (asdout.isEndcap() && asdout.isTriplet()) {
            dbChamber = dbChamber - 1;
        }
        int indexIn[TGCDatabaseASDToPP::NIndexIn] = {
            asdout.getLayer(), dbChamber, asdout.getChannel()};
        int i = databaseP->getIndexDBIn(indexIn);
        if (i < 0) {
            return nullptr;
        }
        channel = databaseP->getEntry(i, 7) + 1;
    } else {
        if ((asdout.isBackward() && asdout.isAside()) ||
            (!asdout.isBackward() && asdout.isCside())) {
            channel = 32 - asdout.getChannel();
        } else {
            channel = asdout.getChannel() + 1;
        }
    }
    if (channel == -1) {
        return nullptr;
    }
 
    return std::make_unique<TGCChannelASDIn>(
        asdout.getSideType(), asdout.getSignalType(), asdout.getRegionType(),
        sector, asdout.getLayer(), chamber, channel);
}
 
std::unique_ptr<TGCChannelId> TGCCableInASD::getChannelOut(
    const TGCChannelId& asdin, bool orChannel) const {
    if (orChannel || asdin.isValid() == false) {
        return nullptr;
    }
 
    const bool asdinisEndcap = asdin.isEndcap();
    const bool asdinisTriplet = asdin.isTriplet();
    const int asdinLayer = asdin.getLayer();
    const int asdinChannel = asdin.getChannel();
 
    TGCDatabase* databaseP =
        m_database.at(+asdin.getRegionType()).at(+asdin.getModuleType()).get();
 
    if (!databaseP) {
        return nullptr;
    }
    // sector ASDIn [1..48, 1..24], ASDOut [2..47.0.1, 1..23.0]
    int sector{0};
    if (asdin.isEndcap()) {
        if (!asdin.isInner()) {
            // Endcap
            sector = (asdin.getSector() + 1) % TGCId::NUM_ENDCAP_SECTOR;
        } else {
            // EI
            sector = (asdin.getSector()) % TGCId::NUM_INNER_SECTOR;
        }
    } else {
        if (!asdin.isInner()) {
            // Forward
            sector = (asdin.getSector()) % TGCId::NUM_FORWARD_SECTOR;
        } else {
            // FI
            sector = (asdin.getSector()) % TGCId::NUM_INNER_SECTOR;
        }
    }
 
    // chamber ASDIn [1(F),1,2,3,4,5(E)], ASDOut [0(F),4,3,2,1,0(E)]
    int chamber{0};
    if (asdin.isEndcap() && !asdin.isInner()) {
        chamber = 5 - asdin.getChamber();
    } else {
        chamber = asdin.getChamber() - 1;
    }
 
    int channel = -1;
    // channel ASDIn [1..32(S),1..n(W chamber)], ASDOut [0..31(S),n..0(W
    // sector)]
    if (asdin.isWire()) {
        const int MaxEntry = databaseP->getMaxEntry();
        for (int i = 0; i < MaxEntry; i++) {
            // Endcap Triplet chamberId start from 1 in ASDOut
            int dbChamber = chamber;
            if (asdinisEndcap && asdinisTriplet) {
                dbChamber = dbChamber - 1;
            }
 
            int dbChannel = asdinChannel - 1;
            if (databaseP->getEntry(i, 7) == dbChannel &&
                databaseP->getEntry(i, 1) == dbChamber &&
                databaseP->getEntry(i, 0) == asdinLayer) {
                channel = databaseP->getEntry(i, 6);
                break;
            }
        }
    } else {
        // asdin.isBackward() can not be used because this method rely on
        // sector number for asdout
        bool is_Backward = false;
        if (asdin.isEndcap()) {
            if (!asdin.isInner()) {
                if (asdin.isAside()) {
                    is_Backward = (sector % 2 == 1);
                } else {
                    is_Backward = (sector % 2 == 0);
                }
            } else {
                // EI
                // Special case of EI11
                if (sector == 15) {
                    if (asdin.isAside()) {
                        is_Backward = false;
                    } else {
                        is_Backward = true;
                    }
                } else if (sector == 16) {
                    if (asdin.isAside()) {
                        is_Backward = true;
                    } else {
                        is_Backward = false;
                    }
                } else {
                    //  A-side phi0 F: phi1 F: phi2 B
                    //  C-side phi0 B: phi1 B: phi2 F
                    if (asdin.isAside()) {
                        is_Backward = (sector % 3 == 2);
                    } else {
                        is_Backward = (sector % 3 != 2);
                    }
                }
            }
        } else {
            if (asdin.isAside()) {
                is_Backward = true;  // All Backward for A-side
            }
        }
        if ((is_Backward && asdin.isAside()) ||
            (!is_Backward && asdin.isCside())) {
            channel = 32 - asdin.getChannel();
        } else {
            channel = asdin.getChannel() - 1;
        }
    }
    if (channel == -1) {
        return nullptr;
    }
 
    return std::make_unique<TGCChannelASDOut>(
        asdin.getSideType(), asdin.getSignalType(), asdin.getRegionType(),
        sector, asdin.getLayer(), chamber, channel);
}
 
}  // namespace MuonTGC_Cabling