TileDetectorTool.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
/*
 * Updates:
 * - 2021, Riccardo Maria BIANCHI <riccardo.maria.bianchi@cern.ch>
 *         * Added TileDetectorFactoryLite to load the Tile geomnetry from SQLite
 */
 
#include "TileGeoModel/TileDetectorTool.h"
#include "TileDetectorFactory.h"
#include "TileAtlasFactory.h"
#include "TileTBFactory.h"
#include "TileDetectorFactoryLite.h"
 
#include "TileDetDescr/TileDetDescrManager.h"
#include "TileDetDescr/TileDddbManager.h"
#include "TileConditions/TileCablingService.h"
 
#include "CaloIdentifier/TileID.h"
#include "TileIdentifier/TileHWID.h"
#include "CaloIdentifier/CaloCell_ID.h"
 
#include "GeoModelInterfaces/IGeoModelSvc.h"
#include "GeoModelUtilities/GeoModelExperiment.h"
#include "GaudiKernel/IService.h"
#include "GaudiKernel/MsgStream.h"
#include "StoreGate/StoreGateSvc.h"
#include "RDBAccessSvc/IRDBAccessSvc.h"
 
#include "GeoModelRead/ReadGeoModel.h"
 
#include "AthenaKernel/ClassID_traits.h"
#include "CxxUtils/checker_macros.h"
#include "SGTools/DataProxy.h"
 
 
TileDetectorTool::TileDetectorTool(const std::string& type,
                                   const std::string& name,
                                   const IInterface* parent)
        : GeoModelTool(type, name, parent)
        , m_switches(false,true)
        , m_not_locked(true)
        , m_useNewFactory(true)
        , m_geometryConfig("FULL")
        , m_manager(0)
{
  declareProperty("UseNewFactory", m_useNewFactory);
  declareProperty("GeometryConfig",m_geometryConfig);
  declareProperty("TestBeam", m_switches.testBeam);
  // declareProperty("AddPlatesToCellVolume", m_switches.addPlatesToCell); to make user aware that the value is taken from DB rather from JO
  declareProperty("Ushape",m_switches.uShape);
  declareProperty("Glue",m_switches.glue);
  declareProperty("PVT",m_switches.pvt);
  declareProperty("Steel",m_switches.steel);
  declareProperty("CsTube",m_switches.csTube);
  declareProperty("CrackOption",m_switches.crackOption);
}
 
 
TileDetectorTool::~TileDetectorTool()
{
}
 
 
void TileDetectorTool::setSwitch(int & param, int value, const char * name)
{
  if (param < 0) {
    param = value;
    ATH_MSG_INFO(name << " parameter from database is: " << param);
  } else {
    if (param != value) {
      ATH_MSG_WARNING("Overriding " << name << " value from DB by value from jobOptions, using "
                      << param << " instead of " << value);
    } else {
      ATH_MSG_INFO(name << " parameter from jobOptions is: " << param);
    }
  }
}
 
 
StatusCode TileDetectorTool::create()
{
  MsgStream log(msgSvc(), name());
  ATH_MSG_INFO(" Entering TileDetectorTool::create()");
 
  // Get the detector configuration
  IGeoModelSvc *geoModel;
  ATH_CHECK(service ("GeoModelSvc",geoModel));
 
  std::string atlasVersion = geoModel->atlasVersion();
  std::string tileVersion = geoModel->tileVersionOverride();
 
  std::string versionTag = tileVersion.empty()? atlasVersion : tileVersion;
  std::string versionNode = tileVersion.empty()? "ATLAS" : "TileCal";
 
  if (versionTag == "AUTO") {
    versionTag = "TileCal-00";
    versionNode = "TileCal";
  }
  if (atlasVersion.compare(0,9,"ATLAS-CTB") == 0 || tileVersion.compare(0,6,"TileTB") == 0) {
    ATH_MSG_INFO("CTB geometry detected: " << atlasVersion  << " " << tileVersion);
    m_switches.testBeam = true;
  }
 
  //Locate the top level experiment node
  GeoModelExperiment* theExpt = nullptr;
  CHECK( detStore()->retrieve(theExpt, "ATLAS") );
 
  if ( 0 == m_detector ) {
    // Get the detector configuration.
    ServiceHandle<IGeoDbTagSvc> geoDbTag("GeoDbTagSvc",name());
    ATH_CHECK(geoDbTag.retrieve());
 
    // Get the 'new' accessSvc to get parameters / DB data from the DD SQLite input file.
    ServiceHandle<IRDBAccessSvc> accessSvc(geoDbTag->getParamSvcName(),name());
    ATH_CHECK(accessSvc.retrieve());
 
    // Get the SQLite reader, if specified in the jobOption
    GeoModelIO::ReadGeoModel* sqliteReader = geoDbTag->getSqliteReader();
 
    // Get the Tile 'DDDB' and 'DetDescr' managers
    // Note: we need to pass the raw IRDBAccess* to the TildDddbManager constructor.
    bool sqliteInput = false;
    if (sqliteReader) sqliteInput = true;
    TileDddbManager_ptr  dbManager(new TileDddbManager(&*accessSvc,versionTag,versionNode, sqliteInput));
    m_manager = new TileDetDescrManager(dbManager);
 
    // check what factory can be used
    if (0==dbManager->GetNumberOfEnv() && m_useNewFactory) {
      ATH_MSG_WARNING("New TileAtlasFactory can not be used because TileGlobals do not exist in Database");
      ATH_MSG_WARNING("Use old TileDetectorFactory instead");
      m_useNewFactory = false;
    }
 
    CHECK( initIds() );
 
    m_switches.addPlatesToCell = dbManager->addPlatesToCell();
    setSwitch(m_switches.uShape, dbManager->uShape(), "Ushape");
    setSwitch(m_switches.glue,   dbManager->glue(),   "Glue");
    setSwitch(m_switches.pvt,    dbManager->PVT(),    "PVT");
    setSwitch(m_switches.steel,  dbManager->steel(),  "Steel");
    setSwitch(m_switches.csTube, dbManager->csTube(), "CsTube");
 
    m_not_locked = false;
 
    GeoPhysVol *world=&*theExpt->getPhysVol();
 
    // build the geometry from the standalone SQLite file
    if (sqliteReader) {
      TileDetectorFactoryLite theTileFactoryLite(detStore().operator->(),
                                                 m_manager,
                                                 sqliteReader,
                                                 accessSvc.operator->(),
                                                 m_switches,
                                                 &log,
                                                 true);
      theTileFactoryLite.create(world);
    } else {
      // build the geometry from the Oracle-based GeometryDB
      if (m_switches.testBeam) {
        // TileTBFactory is not thread-safe. But since this code should only be
        // ever called once (and this is "only" for test beam geometry), we use
        // this rather ugly hack to silence the thread-checker warnings:
        [[maybe_unused]] static const bool do_once ATLAS_THREAD_SAFE = [&]() ATLAS_NOT_THREAD_SAFE {
          TileCablingService::getInstance_nc()->setTestBeam(true);
          TileTBFactory theTileTBFactory = TileTBFactory(detStore().operator->(),m_manager,m_switches,&log);
          theTileTBFactory.create(world);
          return true;
        }();
 
      } else if (m_useNewFactory) {
        std::vector<GeoPhysVol *> volumePtrs;
        std::vector<double> volumePositions;
        TileAtlasFactory theTileFactory(detStore().operator->(),m_manager,m_switches,m_volumes,volumePtrs,volumePositions,&log,m_geometryConfig=="FULL");
        theTileFactory.create(world);
 
      } else {
        TileDetectorFactory theTileFactory(detStore().operator->(),m_manager,m_switches,&log);
        theTileFactory.create(world);
      }
    } // end of building the geometry from the GeometryDB
    ATH_MSG_DEBUG( "The Tile raw geometry has been built.");
 
    CHECK( createElements() );
 
    // Register the TileDetDescrManager instance with the Transient Detector Store
    CHECK( detStore()->record(m_manager, m_manager->getName()) );
    theExpt->addManager(m_manager);
 
    // For reco jobs: release DB manager. Cannot do it here for simulation jobs, they
    // release DB manager as part of global GeoModel release
    if (m_geometryConfig=="RECO")
      m_manager->releaseDbManager();
 
    return StatusCode::SUCCESS;
 
  }
 
  return StatusCode::FAILURE;
}
 
 
StatusCode TileDetectorTool::clear()
{
  SG::DataProxy* proxy = detStore()->proxy(ClassID_traits<TileDetDescrManager>::ID(),m_manager->getName());
  if (proxy) {
    proxy->reset();
    m_manager = 0;
  }
  return StatusCode::SUCCESS;
}
 
 
StatusCode TileDetectorTool::initIds()
{
  const TileID* tileID;
  const CaloCell_ID* cellID;
  const TileHWID* tileHWID;
 
  // Retrieve all helpers from the detector store
  CHECK( detStore()->retrieve(tileID, "TileID") );
  CHECK( detStore()->retrieve(tileHWID, "TileHWID") );
  CHECK( detStore()->retrieve(cellID, "CaloCell_ID") );
 
  m_manager->set_helper(tileID);
  m_manager->set_helper(cellID);
  m_manager->set_helper(tileHWID);
 
  // instantiate Cabling Svc to initialize pointers to helpers there
  const TileCablingService * cabling = TileCablingService::getInstance();
  if (cabling==0) {
    ATH_MSG_ERROR("Could not get instance of TileCablingService");
    return StatusCode::FAILURE;
  }
 
  return StatusCode::SUCCESS;
}
 
 
StatusCode TileDetectorTool::createElements()
{
  m_manager->create_elements(true);
  return StatusCode::SUCCESS;
}