SCT_TestCablingAlg.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/

 
//Package
#include "SCT_TestCablingAlg.h"
 
#include "SCT_CablingUtilities.h"
#include "SCT_CablingXmlTags.h"
 
//Athena
#include "Identifier/IdentifierHash.h"
#include "InDetIdentifier/SCT_ID.h"
 
//Gaudi
#include "GaudiKernel/StatusCode.h"
 
//STL
#include <set>
#include <cstdint>
//for o/p to file
#include <iostream>
#include <fstream>
//to access environment
#include <cstdlib>
#include <algorithm>
#include <map>
 
using namespace SCT_Cabling;
 
SCT_TestCablingAlg::SCT_TestCablingAlg(const std::string& name, ISvcLocator* pSvcLocator):
  AthReentrantAlgorithm(name, pSvcLocator) {
  //nop
}
 
StatusCode
SCT_TestCablingAlg::initialize() {
  ATH_CHECK(detStore()->retrieve(m_idHelper, "SCT_ID"));
  ATH_CHECK(m_cablingTool.retrieve());
  ATH_MSG_INFO("Test algorithm for SCT_Cabling");
  return StatusCode::SUCCESS;
}
 
std::string
SCT_TestCablingAlg::coordString(const Identifier& offlineId) const {
  using std::to_string;
  const std::string sep{", "};
  std::string result{std::string("[") + to_string(m_idHelper->barrel_ec(offlineId)) + sep};
  result += to_string(m_idHelper->layer_disk(offlineId)) + sep;
  result += to_string(m_idHelper->phi_module(offlineId)) + sep;
  result += to_string(m_idHelper->eta_module(offlineId)) + sep;
  result += to_string(m_idHelper->side(offlineId)) + "]";
  return result;
}
 
StatusCode
SCT_TestCablingAlg::execute(const EventContext& ctx) const {
  // const string testAreaPath{CoveritySafe::getenv("TestArea")};
  const std::string testAreaPath = ".";
  std::string filename{testAreaPath+"/cabling.txt"};
  ATH_MSG_INFO("Filename: " << filename << " will be written to your $TestArea.");
  std::ofstream opFile1{filename.c_str(), std::ios::out};
  ATH_MSG_INFO("Executing...");
  ATH_MSG_INFO("hash, offline Id, online Id(hex), serial number");
  const unsigned int nHashesInCabling{2*m_cablingTool->size(ctx)};
  for (unsigned int i{0}; i!=nHashesInCabling; ++i) {
    IdentifierHash hash{i};
    Identifier offlineId{m_idHelper->wafer_id(hash)};
    SCT_OnlineId onlineId{m_cablingTool->getOnlineIdFromHash(hash, ctx)};
    SCT_SerialNumber sn{m_cablingTool->getSerialNumberFromHash(hash, ctx)};
    ATH_MSG_INFO(i << " " << offlineId << " " << std::hex << onlineId << std::dec << " " << sn << " " << coordString(offlineId));
    opFile1 << i << " " << offlineId << " " << std::hex << onlineId << std::dec << " " << sn << " " << coordString(offlineId) << std::endl;
    if (m_cablingTool->getHashFromOnlineId(onlineId, ctx) != hash){
      ATH_MSG_ERROR("?? " << m_cablingTool->getHashFromOnlineId(onlineId, ctx));
    }
  }
  opFile1.close();
  ATH_MSG_INFO("Size: " << m_cablingTool->size(ctx));
  std::vector<unsigned int> rods;
  m_cablingTool->getAllRods(rods, ctx);
  ATH_MSG_INFO("Num. of rods= " << rods.size());
  ATH_MSG_INFO("First rod id " << std::hex << rods[0] << std::dec);
  std::string sn{"20220130000299"};
  ATH_MSG_INFO("Hash from serial number " << m_cablingTool->getHashFromSerialNumber(sn, ctx));
  int tsn{130000299};
  ATH_MSG_INFO("Hash from truncated serial number " << m_cablingTool->getHashFromSerialNumber(tsn, ctx));
  unsigned long long snll{20220130000299LL};
  ATH_MSG_INFO("Hash from truncated serial number " << m_cablingTool->getHashFromSerialNumber(snll, ctx));
  ATH_MSG_INFO("Hash from onlineid " << m_cablingTool->getHashFromOnlineId(0x3d230006, ctx));
  ATH_MSG_INFO("Hash from invalid onlineid " << m_cablingTool->getHashFromOnlineId(0x3d250006, ctx));
  ATH_MSG_INFO("Hash from textfile onlineid " << m_cablingTool->getHashFromOnlineId(0x3d220010, ctx));
  ATH_MSG_INFO("Hash from db onlineid " << m_cablingTool->getHashFromOnlineId(0x3d220105, ctx));
  std::vector<IdentifierHash> hashVec;
  m_cablingTool->getHashesForRod(hashVec, 0x220005, ctx);
  ATH_MSG_INFO("number of hashes for rod 0x220005: " << hashVec.size());
  //new section December 2014
  
  //make a 'rodHash'
  std::sort(rods.begin(), rods.end());
  std::map<unsigned int, unsigned int> rodHashMap;
  for (unsigned int i{0}; i!=rods.size(); ++i) {
    rodHashMap[rods[i]]=i;
  }
  //open a file for writing
  const std::string fullFileName{testAreaPath + "/SCT_Cabling_svc" + (m_POSIXtime ? ("_" +std::to_string(m_POSIXtime)) : "") + ".dat"};
  ATH_MSG_INFO("Open file for write " << fullFileName);
  std::ofstream opFile{fullFileName.c_str(), std::ios::out};
  opFile << XmlHeader << std::endl;
  opFile << OpenRootTag << std::endl;
  opFile << OpenRodMappingTag << std::endl;
  for (auto r: rods) {
    opFile << formatRodOutput(r) << std::endl;
  }
  opFile << CloseRodMappingTag << std::endl;
  opFile << OpenModuleMappingTag << std::endl;
  for (unsigned int i{0}; i!=nHashesInCabling; ++i) {
    IdentifierHash hash{i};
    Identifier offlineId{m_idHelper->wafer_id(hash)};
    SCT_OnlineId onlineId{m_cablingTool->getOnlineIdFromHash(hash, ctx)};
    SCT_SerialNumber sn{m_cablingTool->getSerialNumberFromHash(hash, ctx)};
    //rod, fibre, bec, layerDisk, eta,  phi, side,  robId, sn
    const int bec{m_idHelper->barrel_ec(offlineId)};
    const int layer{m_idHelper->layer_disk(offlineId)};
    const int eta{m_idHelper->eta_module(offlineId)};
    const int phi{m_idHelper->phi_module(offlineId)};
    const int side{m_idHelper->side(offlineId)};
    const unsigned int rodIndex{rodHashMap[onlineId.rod()]};
    opFile << formatModuleOutput(rodIndex, onlineId.fibre(), bec, layer, eta, phi, side, onlineId.rod() ,sn.str()) << std::endl;
  }
  opFile << CloseModuleMappingTag << std::endl;
  opFile << CloseRootTag << std::endl;
  opFile.close();
  return StatusCode::SUCCESS;
}
 
 
StatusCode
SCT_TestCablingAlg::finalize() {
  ATH_MSG_INFO("finalize()");
  return StatusCode::SUCCESS;
}