Lvl1MuCTPIInputPhase1.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
*/
#include <iostream>
#include <iomanip>
 
#include "TrigT1Interfaces/Lvl1MuCTPIInputPhase1.h"
 
namespace LVL1MUONIF {
 
  Lvl1MuCTPIInputPhase1::Lvl1MuCTPIInputPhase1() {
    m_isFilledOutOfTimeCandidates.fill(false);
  }
 
  Lvl1MuCTPIInputPhase1::Lvl1MuCTPIInputPhase1( const Lvl1MuCTPIInputPhase1& right ) {
    *this = right;
  }
 
 
 
  Lvl1MuCTPIInputPhase1& Lvl1MuCTPIInputPhase1::operator=( const Lvl1MuCTPIInputPhase1& right ) {
    if ( this != &right ) {
      clearAll();
      for (size_t idSys=0; idSys<NumberOfMuonSystem; idSys++){
    for( size_t ip=0; ip<right.m_data[idSys].size(); ip++){
      int bc=((right.m_data[idSys]).at(ip)).first;
      Lvl1MuVect vSL(((right.m_data[idSys]).at(ip)).second);
      m_data[idSys].emplace_back(bc,std::move(vSL));
    }
    m_isFilledOutOfTimeCandidates[idSys] = right.m_isFilledOutOfTimeCandidates[idSys];
      } 
    }  
    return *this;
  }
 
 
 
  const Lvl1MuSectorLogicDataPhase1& Lvl1MuCTPIInputPhase1::getSectorLogicData( size_t systemAddress,
                                        size_t subSystemAddress,
                                        size_t sectorAddress,
                                        int    bcid             ) const {
    static const Lvl1MuBarrelSectorLogicDataPhase1 dummy;
    for( size_t ip=0; ip<m_data[systemAddress].size(); ip++){
      int bc=((m_data[systemAddress]).at(ip)).first;
      if (bc != bcid) continue;
      const Lvl1MuVect &vecSL( ((m_data[systemAddress]).at(ip)).second);
      return *(vecSL.at(getSystemIndex(systemAddress,subSystemAddress,sectorAddress)));
    }
    return dummy;
  }
 
  std::shared_ptr<Lvl1MuSectorLogicDataPhase1> Lvl1MuCTPIInputPhase1::getSectorLogicDataPtr( size_t systemAddress,
                                                 size_t subSystemAddress,
                                                 size_t sectorAddress,
                                                 int    bcid             ) const {
    for( size_t ip=0; ip<m_data[systemAddress].size(); ip++){
      int bc=((m_data[systemAddress]).at(ip)).first;
      if (bc != bcid) continue;
      return m_data[systemAddress].at(ip).second.at(getSystemIndex(systemAddress,subSystemAddress,sectorAddress));
    }
    return nullptr;
  }
 
  void Lvl1MuCTPIInputPhase1::setSectorLogicData( const Lvl1MuSectorLogicDataPhase1& data,
                        size_t systemAddress,
                        size_t subSystemAddress,
                        size_t sectorAddress,
                        int    bcid             ) {
    size_t ip=0;
    if (isEmpty(systemAddress,bcid)) {
      ip = reserve( systemAddress, bcid );
    } else {
      ip = getBcidIndex( systemAddress, bcid );
    }
 
    const Lvl1MuVect &vecSL((((m_data[systemAddress]).at(ip)).second));
    size_t idx= getSystemIndex(systemAddress,subSystemAddress,sectorAddress);
    if ( systemAddress == Barrel ) {
      *static_cast<Lvl1MuBarrelSectorLogicDataPhase1*>(vecSL.at(idx).get()) = data;
    } else if ( systemAddress == Endcap ) {
      *static_cast<Lvl1MuEndcapSectorLogicDataPhase1*>(vecSL.at(idx).get())= data;
    } else if ( systemAddress == Forward ) {
      *static_cast<Lvl1MuForwardSectorLogicDataPhase1*>(vecSL.at(idx).get())= data;
    }
  }
 
  void Lvl1MuCTPIInputPhase1::setSectorLogicData( const Lvl1MuSectorLogicDataPhase1::SectorLogicWord & sectorWord,
                        size_t systemAddress,
                        size_t subSystemAddress,
                        size_t sectorAddress,
                        int    bcid              )  {
    // convert to SectorLogicDataPhase1
    std::unique_ptr<Lvl1MuSectorLogicDataPhase1> pData;
    if ( systemAddress == Barrel ) {
      pData = std::make_unique<Lvl1MuBarrelSectorLogicDataPhase1>();
    } else if ( systemAddress == Endcap ) {
      pData = std::make_unique<Lvl1MuEndcapSectorLogicDataPhase1>();
    } else if ( systemAddress == Forward ) {
      pData = std::make_unique<Lvl1MuForwardSectorLogicDataPhase1>();
    } else {
      return ;
    }
    pData->convertFromWordFormat( sectorWord );
    setSectorLogicData( *pData, systemAddress, subSystemAddress, sectorAddress, bcid );
  }
 
  void Lvl1MuCTPIInputPhase1::merge( const Lvl1MuCTPIInputPhase1& right ) {
    if ( this == &right ) return;
    
    for (size_t idSys=0; idSys<NumberOfMuonSystem; idSys++){
      for( size_t ip=0; ip<right.m_data[idSys].size(); ip++){
    int bc=((right.m_data[idSys]).at(ip)).first;
    if (isEmpty( idSys, bc)){
      Lvl1MuVect vSL(((right.m_data[idSys]).at(ip)).second);
      m_data[idSys].emplace_back(bc,std::move(vSL));
      m_isFilledOutOfTimeCandidates[idSys] = right.m_isFilledOutOfTimeCandidates[idSys];
    }
      } 
    }  
 
  }
 
  size_t Lvl1MuCTPIInputPhase1::reserve( size_t systemAddress ,
                   int    bcid)          {
    
    Lvl1MuVect vecSL;
    
    if ( systemAddress == Barrel ) {
      vecSL.reserve(NumberOfBarrelSector*NumberOfMuonSubSystem);
      for ( size_t id = 0; id < NumberOfBarrelSector*NumberOfMuonSubSystem; id++ ) {
    vecSL.emplace_back(std::make_shared<Lvl1MuBarrelSectorLogicDataPhase1>());
      }
 
    } else if ( systemAddress == Endcap ) {
      vecSL.reserve(NumberOfEndcapSector*NumberOfMuonSubSystem);
      for ( size_t id = 0; id < NumberOfEndcapSector*NumberOfMuonSubSystem; id++ ) {
    vecSL.emplace_back(std::make_shared<Lvl1MuEndcapSectorLogicDataPhase1>());
      }
 
    } else if ( systemAddress == Forward ) {
      vecSL.reserve(NumberOfForwardSector*NumberOfMuonSubSystem);
      for ( size_t id = 0; id < NumberOfForwardSector*NumberOfMuonSubSystem; id++ ) {
    vecSL.emplace_back(std::make_shared<Lvl1MuForwardSectorLogicDataPhase1>());
      }
    }
 
    size_t ip = m_data[systemAddress].size();
    m_data[systemAddress].emplace_back( bcid, std::move(vecSL));
 
    if (bcid!=0) m_isFilledOutOfTimeCandidates[systemAddress] = true; 
 
    return ip;
  }
 
  bool Lvl1MuCTPIInputPhase1::isEmpty( size_t systemAddress,
                 int    bcid           ) const {
    for( size_t ip=0; ip<m_data[systemAddress].size(); ip++){
      int bc=((m_data[systemAddress]).at(ip)).first;
      if (bc == bcid) return false;;
    }
    return true;
  }
 
  bool Lvl1MuCTPIInputPhase1::isEmptyAll(int bcid) const {
    if ( isEmpty(Barrel,bcid) && isEmpty(Endcap,bcid) && isEmpty(Forward,bcid) ) {
    return true;
      } else {
    return false;
      }
  }
 
  void Lvl1MuCTPIInputPhase1::clear( size_t systemAddress ) {
    m_data[systemAddress].clear();
  }
 
  void Lvl1MuCTPIInputPhase1::clearAll() {
    for ( size_t id = 0; id < NumberOfMuonSystem; id++ ) {
      m_data[ id ].clear();
    }
 
 
  }
 
  std::ostream& operator<<( std::ostream& out, const Lvl1MuCTPIInputPhase1& right ) {
 
    size_t systemAddress;
    size_t subSystemAddress;
    size_t side, id;
 
    out << "--------- Barrel Sector ---------------------------" << std::endl;
    systemAddress = right.idBarrelSystem();
    for ( side = 0; side < right.NumberOfMuonSubSystem; ++side ) {
      if ( side == 0 ) {
        out << "   ------ A Side (Z>0)  ---------------------------" << std::endl;
        subSystemAddress = right.idSideA();
      } else {
        out << "   ------ C Side (Z<0)  ---------------------------" << std::endl;
        subSystemAddress = right.idSideC();
      }
      for ( id = 0; id < right.numberOfBarrelSector() ;++id ) {
        out << "      --- Sector ID  :  " << id << "   ----------" << std::endl;
    for( size_t ip=0; ip < right.m_data[systemAddress].size(); ip++){
      int bc=((right.m_data[systemAddress]).at(ip)).first;
      out << right.getSectorLogicData( systemAddress, subSystemAddress, id, bc );
      out << " BC: " << bc;
      out << std::endl;
    }
      }
    }
 
    out << "--------- Endcap Sector ---------------------------" << std::endl;
    systemAddress = right.idEndcapSystem();
    for ( side = 0; side < right.NumberOfMuonSubSystem; ++side ) {
      if ( side == 0 ) {
        out << "   ------ A Side (Z>0)  ---------------------------" << std::endl;
        subSystemAddress = right.idSideA();
      } else {
        out << "   ------ C Side (Z<0)  ---------------------------" << std::endl;
        subSystemAddress = right.idSideC();
      }
      for ( id = 0; id < right.numberOfEndcapSector() ;++id ) {
        out << "      --- Sector ID  :  " << id << "   ----------" << std::endl;
    for( size_t ip=0; ip < right.m_data[systemAddress].size(); ip++){
      int bc=((right.m_data[systemAddress]).at(ip)).first;
      out << right.getSectorLogicData( systemAddress, subSystemAddress, id, bc );
      out << " BC: " << bc;
      out << std::endl;
    }
      }
    }
 
    out << "--------- Forward Sector ---------------------------" << std::endl;
    systemAddress = right.idForwardSystem();
    for ( side = 0; side < right.NumberOfMuonSubSystem; ++side ) {
      if ( side == 0 ) {
        out << "   ------ A Side (Z>0)  ---------------------------" << std::endl;
        subSystemAddress = right.idSideA();
      } else {
        out << "   ------ C Side (Z<0)  ---------------------------" << std::endl;
        subSystemAddress = right.idSideC();
      }
      for ( id = 0; id < right.numberOfForwardSector() ;++id ) {
        out << "      --- Sector ID  :  " << id << "   ----------" << std::endl;
    for( size_t ip=0; ip < right.m_data[systemAddress].size(); ip++){
      int bc=((right.m_data[systemAddress]).at(ip)).first;
      out << right.getSectorLogicData( systemAddress, subSystemAddress, id, bc );
      out << " BC: " << bc;
      out << std::endl;
    }
      }
    }
 
    return out;
  }
}