TRT_AlignDbSvc.cxx

Go to the documentation of this file.

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

 
#include "TRT_AlignDbSvc.h"
 
 
 
#include "RegistrationServices/IIOVRegistrationSvc.h"
#include "AthenaKernel/IAthenaOutputStreamTool.h"
#include "InDetIdentifier/TRT_ID.h"
 
#include "TRT_ReadoutGeometry/TRT_DetElementCollection.h"
#include "TRT_ReadoutGeometry/TRT_BaseElement.h"
#include "ReadoutGeometryBase/InDetDD_Defs.h"
#include "TRT_ReadoutGeometry/TRT_DetectorManager.h"
 
#include "TRT_ConditionsData/FloatArrayStore.h"
 
// Amg
#include "EventPrimitives/EventPrimitives.h"
#include "GeoPrimitives/GeoPrimitives.h"
#include "GeoPrimitives/CLHEPtoEigenConverter.h"
#include "GeoPrimitives/CLHEPtoEigenEulerAnglesConverters.h"
#include <fstream>
#include <iostream>
#include <sstream>
 
ATLAS_NO_CHECK_FILE_THREAD_SAFETY; // This class uses const_cast and regFcn (callback). Legacy code
 
TRT_AlignDbSvc::TRT_AlignDbSvc( const std::string& name, ISvcLocator* pSvcLocator )
  : base_class(name,pSvcLocator),
    m_detStore("DetectorStore",name),
    m_trtStrawAlignDbSvc("TRT_StrawAlignDbSvc",name),
    m_trtid(nullptr),
    m_trtman(nullptr),
    m_alignroot("/TRT/Align"),
    m_alignString("AL"),
    m_par_alitextfile(""),
    m_streamer("AthenaOutputStreamTool/CondStream1"),
    m_alignDBprefix("/TRT/Align/"),
    m_dynamicDB(false),
    m_forceUserDBConfig(false)
{
  declareProperty("StreamTool",m_streamer);
  declareProperty("alignroot",m_alignroot);
  declareProperty("alignTextFile",m_par_alitextfile);
  declareProperty("DetectorStore",m_detStore);
  declareProperty("TrtStrawAlignDbSvc",   m_trtStrawAlignDbSvc,"Service for interaction with the TRT straw alignment DB");
  declareProperty("alignString",m_alignString);
  declareProperty("alignDBprefix",m_alignDBprefix);
  declareProperty("forceUserDBConfig",m_forceUserDBConfig, "Set to true to override any DB auto-configuration");
}
 
//-------------------------------------------------------------------------
 
TRT_AlignDbSvc::~TRT_AlignDbSvc()
{}
 
//
//-----------------------------------------------------------------------------
StatusCode TRT_AlignDbSvc::initialize() 
{
  ATH_MSG_INFO( " in initialize "  );
  
  ATH_CHECK( m_detStore.retrieve() );
  ATH_CHECK (m_trtStrawAlignDbSvc.retrieve() );
  ATH_MSG_INFO ("retrieved " << m_trtStrawAlignDbSvc);
 
  // Get the geometry.
  InDetDD::TRT_DetElementCollection::const_iterator iter,itermin,itermax;
  ATH_CHECK (m_detStore->retrieve(m_trtman,"TRT") );
 
  // Get TRT manager and ID helper
  itermin=m_trtman->getDetectorElementBegin();
  itermax=m_trtman->getDetectorElementEnd();
  ATH_CHECK( m_detStore->retrieve(m_trtid,"TRT_ID") );
  ATH_MSG_INFO( "found all TRT services "  );
 
  if (m_trtman->m_alignfoldertype == InDetDD::timedependent_run2 && !m_forceUserDBConfig){
    m_dynamicDB = true;
    m_alignroot = "/TRT/AlignL2";                                                         
    m_alignDBprefix = "/TRT/AlignL2/";
  }
  if (m_trtman->m_alignfoldertype == InDetDD::static_run1 && !m_forceUserDBConfig){
    m_dynamicDB = false;
    m_alignroot = "/TRT/Align";
    m_alignDBprefix = "/TRT/Align/";
  }
 
 
  // setup list of TDS objects from geometry description
  m_alignobjs.clear();
  m_alignchans.clear();
 
  // Include a level 1 alignment
  int ichan(0) ;
  if (!m_dynamicDB){
    m_alignobjs.emplace_back("/TRT/Align/TRT");
    ATH_MSG_INFO( "Adding key: /TRT/Align/TRT --> We are using static DB folder scheme"  );
    m_alignchans.push_back(ichan++);
  }
  
  //Check all detector elements in the present geometry setup
  for (iter=itermin;iter!=itermax;++iter) {
    const InDetDD::TRT_BaseElement* element=*iter;
    if (element!=nullptr) {
      const Identifier ident=element->identify();
      if (m_trtid->is_trt(ident)) {           //OK this Element is included
    std::string key = findkey(ident,m_alignString); 
    std::vector<std::string>::const_iterator ix = find(m_alignobjs.begin(),m_alignobjs.end(),key);
    if (ix==m_alignobjs.end()) {
      m_alignobjs.push_back(key);
      ATH_MSG_INFO( "Adding key: " << key  );
      m_alignchans.push_back(ichan++);
    }
      }
    }
  }
  

  bool alignFolderExists = m_detStore->contains<AlignableTransformContainer>(m_alignroot) ;
  if (m_dynamicDB) alignFolderExists = ( alignFolderExists && m_detStore->contains<CondAttrListCollection>("/TRT/AlignL1/TRT"));

  bool alignTextFileExists = !m_par_alitextfile.empty();
 
  if( alignFolderExists ) {    
    
    if( alignTextFileExists ) 
      ATH_MSG_INFO( "AlignableTransformContainer with name " << m_alignroot
                    << " exists. Will read text file from callback function."  );
  
  } else if ( alignTextFileExists ) {
    
    ATH_CHECK( this->createAlignObjects() );
    
    ATH_CHECK( this->readAlignTextFile(m_par_alitextfile) );
  
  } else {
    ATH_MSG_ERROR( "AlignableTransformContainer not in IOV service and no textfile specified."  );
  }
  
  ATH_MSG_DEBUG ("The TRT/Align keys and channels are:");
  for (unsigned int i=0;i<m_alignobjs.size();++i)
    ATH_MSG_DEBUG (" key " << m_alignobjs[i] << " chan " << m_alignchans[i]);
  
  return StatusCode::SUCCESS;
}
//
//-----------------------------------------------------------------------------
 
StatusCode TRT_AlignDbSvc::finalize()
{
  ATH_MSG_INFO( "TRT_AlignDbSvc finalize method called"  );
  for(std::vector<Amg::Transform3D*>::iterator it = m_amgTransformCache.begin(); it != m_amgTransformCache.end(); ++it){
      delete *it;
  }
  return StatusCode::SUCCESS;
}

StatusCode TRT_AlignDbSvc::IOVCallBack()
{
  ATH_MSG_DEBUG( "In IOVCallBack"  );
  
  ATH_CHECK( createAlignObjectsWhichDoNotAlreadyExist() );
  
//  for (std::list<std::string>::const_iterator itr=keys.begin(); itr!=keys.end(); ++itr)
//    ATH_MSG_INFO( "IOVCALLBACK for key " << *itr<< " number " << I  );
  
  
  if(!m_par_alitextfile.empty()){
    ATH_CHECK( this->readAlignTextFile(m_par_alitextfile) );
  }else{
    ATH_MSG_INFO( "Text file " << m_par_alitextfile << " is empty"  );
  }

  printCondObjects() ;
 
  ATH_MSG_DEBUG( "Leaving IOVCallBack"  );
  return StatusCode::SUCCESS;
}

StatusCode TRT_AlignDbSvc::writeAlignTextFile(const std::string &file) const
{
  ATH_MSG_DEBUG(" in writeAlignTextFile ");
  ATH_MSG_INFO( " Write AlignableTransforms to text file: "<< file );
  std::ofstream  outfile(file);
 
  int ntrans=0;
  int nobj=0;

  std::vector<std::string>::const_iterator iobj=m_alignobjs.begin();
  std::vector<std::string>::const_iterator iobjE=m_alignobjs.end();
  for (;iobj!=iobjE; ++iobj) {
    const AlignableTransform* pat=cgetTransPtr(*iobj);
    std::string key=*iobj;
    
    if(isOldKey(key) && m_alignString != "ALold")
      continue;

    if (!pat){
      ATH_MSG_ERROR( "Cannot find AlignableTransform for key "<< key  );
      ATH_MSG_ERROR( "Skipping it. "  );
      outfile << key << std::endl;
      continue;
    }
    
    // first record is the name
    ATH_MSG_INFO( " Write folder: " << key  );
    outfile << key << std::endl;
    ++nobj;
    
    // loop over all transforms in the AlignableTransform object
    AlignableTransform::AlignTransMem_citr cit= pat->begin();
    AlignableTransform::AlignTransMem_citr citE= pat->end();
    for (;cit!=citE;++cit) {
      
      int bec(0);
      int layer(0);
      int sector(0);
      int strawLayer(0);
      
      ++ntrans;
      
      const Identifier& ident=cit->identify();
      bec=m_trtid->barrel_ec(ident);
      layer=m_trtid->layer_or_wheel(ident);
      sector=m_trtid->phi_module(ident);
      strawLayer=m_trtid->straw_layer(ident);
       
      const Amg::Transform3D& trans = Amg::CLHEPTransformToEigen(cit->transform());
 
      if(msgLvl(MSG::DEBUG)){
    std::string thisMess = key+" for "+file;
    printTransform(thisMess,trans);
      }
 
 
      HepGeom::Transform3D transCLHEP  =  Amg::EigenTransformToCLHEP(trans);
      CLHEP::HepRotation rot=transCLHEP.getRotation(); 
      Amg::Vector3D eulerangles;
      eulerangles[0] = rot.getPhi();
      eulerangles[1] = rot.getTheta();
      eulerangles[2] = rot.getPsi();
 
      Amg::Vector3D shift=trans.translation();
      //Matthias: was initially 2,1,2 --> this should be equivalent to complicated method above
      //Amg::Vector3D eulerangles = trans.rotation().eulerAngles(2,0,2) ; 
 
      if(key == (m_alignDBprefix+"L2A") || key == (m_alignDBprefix+"L2C") ){
 
    ATH_MSG_INFO( "   Write member: "
                      << ident << " " << bec << " " << layer << " " << strawLayer << " "  << key  );
    
    outfile << bec << " " << layer << " " << strawLayer << " "
         << std::setprecision(10) << " "
         << shift.x() << " " << shift.y() << " " << shift.z() << " "
         << eulerangles[0] << " " << eulerangles[1] << " " << eulerangles[2]
         << std::endl;
 
    ATH_MSG_DEBUG( "Writing transform to log file "  );
    ATH_MSG_DEBUG( bec << " " << layer << " " << strawLayer << " "
            << std::setprecision(10) << " "
            << shift.x() << " " << shift.y() << " " << shift.z() << " "
                << eulerangles[0] << " " << eulerangles[1] << " " << eulerangles[2]
                       );
 
      }else{
    //Need to validate this for writing out L1 Endcap alignments 
    ATH_MSG_INFO( "   Write member: "
                      << ident << " " << bec << " " << layer << " " << sector << " "  << key  );
    
    outfile << bec << " " << layer << " " << sector << " "
         << std::setprecision(10) << " "
         << shift.x() << " " << shift.y() << " " << shift.z() << " "
         << eulerangles[0] << " " << eulerangles[1] << " " << eulerangles[2]
         << std::endl;
 
    ATH_MSG_DEBUG( "Writing transform to log file "  );
    ATH_MSG_DEBUG( bec << " " << layer << " " << sector << " "
            << std::setprecision(10) << " "
            << shift.x() << " " << shift.y() << " " << shift.z() << " "
                << eulerangles[0] << " " << eulerangles[1] << " " << eulerangles[2]
                       );
      }
      
    }
  }
  
  outfile.close();
  
  ATH_MSG_INFO( "Written " << nobj << " AlignableTransform objects" << " with " << ntrans 
         << " transforms to text file" );
  
  ATH_MSG_DEBUG( " leaving writeAlignTextFile " );
  return StatusCode::SUCCESS;
}
 
StatusCode TRT_AlignDbSvc::writeStrawAlignTextFile(const std::string & file) const
{
  ATH_MSG_DEBUG( " in writeAlignTextFile "  );
  
  ATH_MSG_INFO( "Writing the straw  alignments "  );
  if( m_trtStrawAlignDbSvc -> writeTextFile(file) != StatusCode::SUCCESS )
    ATH_MSG_ERROR("Cannot write to file "<< file  );
  
  return StatusCode::SUCCESS;
}

StatusCode TRT_AlignDbSvc::readAlignTextFile(const std::string & file) {
  ATH_MSG_DEBUG( " In  readAlignTextFile " );
  ATH_MSG_INFO( "Read alignment constants from text file: " << file );
  std::ifstream infile;
  infile.open(file.c_str());

  const AlignableTransformContainer* container= getContainer();
  if(container) {
    AlignableTransformContainer* nccontainer = const_cast<AlignableTransformContainer*>(container) ;
 
    DataVector<AlignableTransform>::iterator it=nccontainer->begin();
    DataVector<AlignableTransform>::iterator itE=nccontainer->end();
    for(;it!=itE;++it)
      // there is no 'clear' function in alignable transform, so we do it with this hack
      **it = AlignableTransform( (*it)->tag() ) ;
  } else {
    return StatusCode::FAILURE ;
  }
  
  // loop over file
  int nobj=0;
  int ntrans=0;
  std::string atname;
  AlignableTransform* pat(nullptr) ;
  char line[512] ;
  
  while( infile.getline(line,512) ) {
    if(line[0] != '#') {
      
      std::string linestring(line) ;
      if( linestring.find('/') != std::string::npos) {
    // this must be a line with a container name
    atname = linestring ;
    ATH_MSG_INFO( "now reading container: " << atname  );
    pat=getTransPtr(atname);
    if (!pat) {
      ATH_MSG_ERROR( "Cannot find AlignableTransform object for key " << atname  );
      return StatusCode::FAILURE;
    } else {
      nobj++;
    } 
      } else if(pat!=nullptr) {
    // this must be a line with constants
    std::istringstream is(line) ;
    
    int bec(0),layer(0),sector(0),strawLayer(0);
    double dx,dy,dz,phi,theta,psi;
    
    if(atname == (m_alignDBprefix+"L2A") || atname == (m_alignDBprefix+"L2C") ){
      if( is >> bec >> layer >> strawLayer >> dx >> dy >> dz >> phi ) {
         CLHEP::Hep3Vector translation(dx, dy, dz);
         CLHEP::HepRotation rotation;
         if (is >> theta >> psi) {
           ATH_MSG_DEBUG ("read a line with euler angles!" << phi
                          << " " << theta << " " << psi);
            rotation.set(phi, theta, psi);
         }
        //We have to give a phi sector, so we'll give 0.
        Identifier ident=m_trtid->layer_id(bec,0,layer,strawLayer);
        ATH_MSG_DEBUG( "The identifier is for bec " << bec 
                           << " layer " << layer << " strawLayer " << strawLayer  );
        ATH_MSG_DEBUG( "The code is  " << m_trtid->show_to_string(ident)  );
        if( pat->findIdent(ident)!=pat->end() ) {
          ATH_MSG_WARNING ("WARNING: read module from file which was already in AlignableTransform. Will skip it."
                               << " bec,lay,sec,strawlay = " << bec << " " << layer << " " << sector << " " );
        } else {
          pat->add(ident,HepGeom::Transform3D(rotation, translation));
          HepGeom::Transform3D testtrans = HepGeom::Transform3D(rotation, translation);
          if(msgLvl(MSG::DEBUG)){
        std::string thisMess = atname+" for "+file;
        printTransform(thisMess,Amg::CLHEPTransformToEigen(testtrans));
          }
 
 
          ++ntrans;
          ATH_MSG_DEBUG( "Set transform: "
                             << " [" << bec << "," << layer << ","  << strawLayer 
                             << "] key " << atname 
                             << " rotation=(" << phi << "," << theta << "," << psi
                             << "),  translation=(" << dx << "," << dy << "," << dz << "])" 
                             );
        }
      } else if(!linestring.empty()) {
        ATH_MSG_WARNING( "Read invalid line from textfile. Line=\'" << line << "\'"  );
      }
    }//if endcap
    else{
      if( is >> bec >> layer >> sector >> dx >> dy >> dz >> phi ) {
        CLHEP::Hep3Vector translation(dx,dy,dz)  ;
        CLHEP::HepRotation rotation;
        if( is >> theta >> psi ) {
          ATH_MSG_DEBUG( "read a line with euler angles!" << phi << " " << theta << " " << psi );
            rotation.set(phi,theta,psi) ;
 
        } else {
          ATH_MSG_DEBUG( "read a line with angle in phi!" << phi  );
          rotation = CLHEP::HepRotationZ(phi) ;
        }
        
        Identifier ident=m_trtid->module_id(bec,sector,layer);
        if( pat->findIdent(ident)!=pat->end() ) {
          ATH_MSG_WARNING ("WARNING: read module from file which was already in AlignableTransform. Will skip it."
                               << " bec,lay,sec = " << bec << " " << layer << " " << sector);
        } else {
          pat->add(ident,HepGeom::Transform3D(rotation, translation));
          ++ntrans ;
          ATH_MSG_DEBUG( "Set transform: "
                             << " [" << bec << "," << layer << "," << sector 
                             << "] key " << atname << " rotation=(" << phi << "," << theta << "," << psi
                             << "),  translation=(" << dx << "," << dy << "," << dz << "])"  );
        }
      } else if(!linestring.empty()) {
        ATH_MSG_WARNING( "Read invalid line from textfile. Line=\'" << line << "\'"  );
      }
    }//if not the endcap
 
    
      }//if pat!=0
    }//if line != "#"
  }//while 
 
  infile.close() ;
  ATH_MSG_INFO( "Read " << nobj << " objects from file with " << ntrans << " transforms." 
                << " Now forcing callback in detector manager."  );

  if((const_cast<InDetDD::TRT_DetectorManager*>(m_trtman))->align().isFailure()){
    ATH_MSG_ERROR("Failed to force the alignment callback!" );
  }
  
  ATH_MSG_DEBUG( " Leaving  readAlignTextFile "  );
  return StatusCode::SUCCESS;
}  

StatusCode TRT_AlignDbSvc::streamOutAlignObjects() const{
  ATH_MSG_DEBUG( "In streamOutAlignObjects "   );
 
  // Get Output Stream tool for writing
  ATH_CHECK( m_streamer.retrieve() );
  
  IAthenaOutputStreamTool*  streamer=const_cast<IAthenaOutputStreamTool*>(&(*m_streamer));
  
  ATH_CHECK( streamer->connectOutput() );
  
  IAthenaOutputStreamTool::TypeKeyPairs  typeKeys;
  IAthenaOutputStreamTool::TypeKeyPair arraypair("AlignableTransformContainer",m_alignroot);
  typeKeys.push_back(arraypair);
  ATH_CHECK( m_detStore->contains<AlignableTransformContainer>(m_alignroot) );
  
  ATH_CHECK( streamer->streamObjects(typeKeys) );
  ATH_CHECK( streamer->commitOutput() );
  
  ATH_MSG_INFO( "   Streamed out and committed AlignableTransformContainer"  );
  printCondObjects() ;
 
  ATH_MSG_DEBUG( "Leaving streamOutAlignObjects "   );
  return StatusCode::SUCCESS;
}

StatusCode TRT_AlignDbSvc::registerAlignObjects(const std::string & tag, int run1, int event1, int run2, int event2) const{
  
  ATH_MSG_INFO( "registerAlignObjects with IOV "  );
  ATH_MSG_INFO( "Run/evt1 [" << run1 << "," << event1 << "]"  );
  ATH_MSG_INFO( "Run/evt2 [" << run2 << "," << event2 << "]"  );
 
  // get pointer to registration svc
  SmartIF<IIOVRegistrationSvc> regsvc(service("IOVRegistrationSvc"));
  ATH_CHECK( regsvc.isValid() );
  
  if (StatusCode::SUCCESS==regsvc->registerIOV("AlignableTransformContainer",m_alignroot,tag,run1,run2,event1,event2)){
    ATH_MSG_INFO( " Register AlignableTransformContainer object " 
                  << m_alignroot  );
  } else {
    ATH_MSG_ERROR( " Failed to register AlignableTranformContainer " 
                   << m_alignroot  );
    return StatusCode::FAILURE;
  }
  
  return StatusCode::SUCCESS;
}

const Amg::Transform3D TRT_AlignDbSvc::getAlignmentTransform(const Identifier& ident, unsigned int level) const {
  ATH_MSG_DEBUG( "In getAlignmentTransform"  );
  
  if(level != 1 && level != 2){
    ATH_MSG_ERROR( "Call to trans wrong!"  );
    ATH_MSG_ERROR( "Level must be 1 or 2!"  );
    ATH_MSG_ERROR( "Returning NULL"  );
    return Amg::Transform3D();
  }
 
  return (level == 1) ? getAlignmentTransformL1(ident) : getAlignmentTransformL2(ident);
}

const Amg::Transform3D  TRT_AlignDbSvc::getAlignmentTransformL1(Identifier const& ident) const{
  ATH_MSG_DEBUG( "In getAlignmentTransformL1"  );
  // return level 1 AlignableTransform for the TRT subdetector bec (+-1,+-2)
  // or return null transform.
  
  std::string key="/TRT/Align/TRT";
  const AlignableTransform* pat=cgetTransPtr(key);
  if (!pat) return Amg::Transform3D();
 
  // get Amg::Transform3d via AlignTransMember
  AlignableTransform::AlignTransMem_citr itr = pat->findIdent(ident);
  return ( itr!=pat->end() ? Amg::CLHEPTransformToEigen(itr->transform()) : Amg::Transform3D() );
}

const Amg::Transform3D TRT_AlignDbSvc::getAlignmentTransformL2(Identifier const& ident) const{
  ATH_MSG_DEBUG( "In getAlignmentTransformL2"  );
  // return level 2 AlignableTransform for the module containing ident
  // or return null transform.
  
  std::string key=findkey(ident,m_alignString);
  //does a key exist corresponding to this identifier?
  
  const AlignableTransform* pat=cgetTransPtr(key);
  if (!pat) return Amg::Transform3D();
 
  // OK it exists
  // get "module identifier" from identifier
  int bec=m_trtid->barrel_ec(ident);
  int layer=m_trtid->layer_or_wheel(ident);
  int sector=m_trtid->phi_module(ident);
  int strawLayer=m_trtid->straw_layer(ident);
  int ring = getRingForStrawLayer(strawLayer);
 
  Identifier mid;
  if(key == (m_alignDBprefix+"L2A") || key == (m_alignDBprefix+"L2C") )
    mid=m_trtid->layer_id(bec,0,layer,ring);
  else
    mid=m_trtid->module_id(bec,sector,layer);
  
  // get Amg::Transform3d via AlignTransMember
  AlignableTransform::AlignTransMem_citr itr = pat->findIdent(mid);
  return ( itr!=pat->end() ? Amg::CLHEPTransformToEigen(itr->transform()) : Amg::Transform3D() );
}

const Amg::Transform3D* TRT_AlignDbSvc::getAlignmentTransformPtr(const Identifier& ident, unsigned int level) const {
  ATH_MSG_DEBUG( "In getAlignmentTransformPtr"  );
  
  if(level != 1 && level != 2){
    ATH_MSG_ERROR( "Call to trans wrong!"  );
    ATH_MSG_ERROR( "Level must be 1 or 2!"  );
    ATH_MSG_ERROR( "Returning NULL"  );
    return nullptr;
  }
 
  return (level == 1) ? getAlignmentTransformL1Ptr(ident) : getAlignmentTransformL2Ptr(ident);
}

const Amg::Transform3D*  TRT_AlignDbSvc::getAlignmentTransformL1Ptr(Identifier const& ident) const{
  ATH_MSG_DEBUG( "In getAlignmentTransformL1Ptr"  );
  // get level 1 (despite name, will change soon) AlignableTransform for the subdetector containing ident
  const Amg::Transform3D* rc(nullptr) ;
  if( m_trtid->is_trt(ident) ) {
    const AlignableTransform* pat = cgetTransPtr( "/TRT/Align/TRT" ) ;
    if( pat ) {
      int bec=m_trtid->barrel_ec(ident);
      Identifier mid=m_trtid->module_id(bec,0,0);
      AlignableTransform::AlignTransMem_citr iter = pat->findIdent(mid) ;
      if(iter != pat->end() ) {
          Amg::Transform3D* amgTransform = new Amg::Transform3D((Amg::CLHEPTransformToEigen(iter->transform())));
          m_amgTransformCache.push_back(amgTransform);
          rc = amgTransform;
      }
    }
  }
  
  ATH_MSG_DEBUG( "Leaving getAlignmentTransformL1Ptr"  );
  return rc ;
}

const Amg::Transform3D* TRT_AlignDbSvc::getAlignmentTransformL2Ptr(Identifier const& ident) const{
  ATH_MSG_DEBUG( "In getAlignmentTransformL2Ptr"  );
  // set level 2 (despite name, will change soon) AlignableTransform for the module containing ident
  // or add a new one.
  const Amg::Transform3D* rc(nullptr) ;
  //does a folder exist corresponding to this identifier?
  if( m_trtid->is_trt(ident) ) {
    std::string key=findkey(ident,m_alignString);
    const AlignableTransform* pat=cgetTransPtr(key);
    if(pat) {
      // make sure the identifier is a "module identifier"
      int bec=m_trtid->barrel_ec(ident);
      int layer=m_trtid->layer_or_wheel(ident);
      int sector=m_trtid->phi_module(ident);
      int strawLayer=m_trtid->straw_layer(ident);
      int ring = getRingForStrawLayer(strawLayer);
      Identifier mid;
      if(key == (m_alignDBprefix+"L2A") || key == (m_alignDBprefix+"L2C") ){
    mid =m_trtid->layer_id(bec,0,layer,ring);
      }else
    mid =m_trtid->module_id(bec,sector,layer);
      
      AlignableTransform::AlignTransMem_citr iter = pat->findIdent(mid) ;
      if(iter != pat->end() ){
          AlignableTransform::AlignTransMem_citr iter = pat->findIdent(mid) ;
          Amg::Transform3D* amgTransform = new Amg::Transform3D((Amg::CLHEPTransformToEigen(iter->transform())));
          m_amgTransformCache.push_back(amgTransform);
          rc = amgTransform;
      }
 
      else
    ATH_MSG_WARNING( "Did not find the transform for " << m_trtid->show_to_string(mid)  );
    }
  }
  
  ATH_MSG_DEBUG( "Leaving getAlignmentTransformL2Ptr"  );
  return rc ;
}

StatusCode TRT_AlignDbSvc::setAlignTransform(Identifier ident, Amg::Transform3D trans, unsigned int level) {
  ATH_MSG_DEBUG( "In getAlignmentTransform"  );
  
  if(level != 1 && level != 2 && level != 3){
    ATH_MSG_FATAL( "Call to setAlignTransform wrong!"  );
    ATH_MSG_FATAL( "Level must be 1,2 or 3!"  );
    return StatusCode::FAILURE;
  }
 
  switch (level) {
    
  case 1:
    ATH_CHECK( setAlignTransformL1(ident,trans) );
    break;
  case 2:
    ATH_CHECK( setAlignTransformL2(ident,trans) );
    break;
  case 3:
    ATH_CHECK( setAlignTransformL3(ident,trans) );
    break;
  }
 
  return StatusCode::SUCCESS;
}

StatusCode TRT_AlignDbSvc::setAlignTransformL1(Identifier ident, const Amg::Transform3D& trans) {
  ATH_MSG_DEBUG( "In setAlignTransformL1"  );
  if(msgLvl(MSG::DEBUG)){
    printTransform("Transform in setAlignTransformL1",trans);
  }
 
  // New additions for new global folder structure
  // No ATs exist for levels 1 & 2 --> need alternative
  if (m_dynamicDB){
    if(!tweakGlobalFolder(ident, trans)) {
      ATH_MSG_ERROR( "Attempt tweak GlobalDB folder failed" );
      return StatusCode::FAILURE;
    }
  }
  else {
  
    //does a folder exist corresponding to this identifier?
    AlignableTransform* pat=getTransPtr("/TRT/Align/TRT");
    if (!pat){
      ATH_MSG_FATAL( "The AlignableTransform for key " << "/TRT/Align/TRT" << " does not exist "  );
      ATH_MSG_FATAL( "Failing ... "  );
      return StatusCode::FAILURE;
    }
    
    // make sure the identifier is a "subdetector identifier"
    if( !(m_trtid->is_trt(ident)) ){
      ATH_MSG_FATAL( "The identifier " << ident << " is not from the TRT "  );
      ATH_MSG_FATAL( "Failing ... "  );
      return StatusCode::FAILURE;
    }
    
    int bec=m_trtid->barrel_ec(ident);
    Identifier mid=m_trtid->module_id(bec,0,0);
    // update or add (mid,trans) pair
    HepGeom::Transform3D clhepTrans = Amg::EigenTransformToCLHEP(trans);
    if( !(pat->update(mid,clhepTrans)) ) {
      pat->add(mid,clhepTrans);
      pat->sortv() ;
    }
  }
 
  ATH_MSG_DEBUG( "Leaving setAlignTransformL1"  );
  return StatusCode::SUCCESS;
}

StatusCode TRT_AlignDbSvc::setAlignTransformL2 (Identifier ident, Amg::Transform3D trans) {
  ATH_MSG_DEBUG( "In setAlignTransformL2"   );
  
  //does a folder exist corresponding to this identifier?
  std::string key=findkey(ident,m_alignString);
  AlignableTransform* pat=getTransPtr(key);
  if (!pat) return StatusCode::FAILURE;
  
  // make sure the identifier is a "module identifier"
  if( !(m_trtid->is_trt(ident)) ) return StatusCode::FAILURE;
  int bec=m_trtid->barrel_ec(ident);
  int layer=m_trtid->layer_or_wheel(ident);
  int sector=m_trtid->phi_module(ident);
  int strawLayer=m_trtid->straw_layer(ident);
  int ring = getRingForStrawLayer(strawLayer);
  Identifier ident2;
 
  ATH_MSG_DEBUG( "Setting the L2 Alignment for: " 
                 << " Bec= " << bec << " layer= " << layer << " sector= " << sector 
                 << " strawLayer= " << strawLayer  );
  ATH_MSG_DEBUG( "The translations are: x= " << trans(0,3) << " y= "<<trans(1,3) << " z= "<<trans(2,3)  );
 
  if(key == (m_alignDBprefix+"L2A") || key == (m_alignDBprefix+"L2C"))
    ident2=m_trtid->layer_id(bec,0,layer,ring);    
  else
    ident2=m_trtid->module_id(bec,sector,layer);
 
  // make sure the identifier is a TRT identifier
  if( !(m_trtid->is_trt(ident2)) ){
    ATH_MSG_FATAL( "The identifier " << ident2 << " is not from the TRT "  );
    ATH_MSG_FATAL( "Failing ... "  );
    return StatusCode::FAILURE;
  }
  
  // update or add (mid,trans) pair
  HepGeom::Transform3D clhepTrans = Amg::EigenTransformToCLHEP(trans);
  if( !(pat->update(ident2,clhepTrans)) ) {
    pat->add(ident2,clhepTrans);
    pat->sortv() ;
  }
  
  ATH_MSG_DEBUG( "Leaving setAlignTransformL2 "   );
  return StatusCode::SUCCESS;
}

StatusCode TRT_AlignDbSvc::setAlignTransformL3 (Identifier ident, Amg::Transform3D trans) {
  ATH_MSG_DEBUG( "In setAlignTransformL3"   );
    
  // make sure the identifier is a "module identifier"
  if( !(m_trtid->is_trt(ident)) ) return StatusCode::FAILURE;
  int bec=m_trtid->barrel_ec(ident);
  
  //Only L3 in the barrel is currently supported 
  if(abs(bec) == 1){
    
    // Currently at L3 only displacements in the plane of the straw layer can be stored.
    // For tracks coming from the interaction point this coorisponds to the sensitive coordinate
    ATH_MSG_INFO("Storing L3 Barrel constants, Only dx and rotz will be written to DB");
    float dx = trans.translation().x();
    float rotz = std::atan2(trans.rotation()(0,1),trans.rotation()(0,0));
 
    // Need the length of the wires to translate the rotation to a translation.
    const InDetDD::TRT_BaseElement* strawElement =  m_trtman->getElement( ident );
    double strawLenthOver2 = 0.5* strawElement->strawLength();
 
    // The alignment frames are the same for straws on side A and side C (verrify this!!)
    // but variables stored in the db are not. We'll calculate the dispacement of each end 
    // in the alignment frame and in a second step, convert to the values in the DB. 
    double delta_dx_atLargerZ = dx + strawLenthOver2 * std::sin(rotz);
    double delta_dx_atSmallerZ = dx - strawLenthOver2 * std::sin(rotz);
    
    // For the definition of dx1 and dx2 see TRT_ConditionsData/StrawDx.h
    // Now we need to know the real meaning of dx1 and dx2. Alas, we
    // are in the frame of either of the two barrel sides, but not
    // both, so we need to know which one.
    bool sideA = m_trtid->barrel_ec(ident) == 1;
    
    // Straw position closest to the electronics 
    double dx1_new = sideA ? delta_dx_atLargerZ : delta_dx_atSmallerZ;
    
    // Straw position away from the electronics 
    double dx2_new = sideA ? delta_dx_atSmallerZ : delta_dx_atLargerZ;
 
    // Uncertianty on straw positions (Arbitrary for now)
    double dxErr = 0.001;
    
    m_trtStrawAlignDbSvc->setDx(ident, dx1_new, dx2_new, dxErr);
  }else{
    // Currently at L3 only displacements in the plane of the straw layer can be stored.
    // For tracks coming from the interaction point this coorisponds to the sensitive coordinate
    ATH_MSG_INFO("Storing L3 Endcap constants, Only dy and rotz will be written to DB");
    
    // The expected transformation is in the local frame of the Endcap Straws
    // In the local straw frame:
    //   - the z-axis is along the straw and points toward the beampipe
    //   - the x-axis is along global-z and away from the interaction point 
    //          (locX = globZ A-side / locX = -1 *gobZ C-side) 
    //   - the y-axis is along global phi_hat direction determined by the other 2. 
    //          (clockwise C-side, counter clockwise A-Side)
    float dy = trans.translation().y();
    float rotx = std::atan2(trans.rotation()(1,2),trans.rotation()(2,2));
 
    // Need the length of the wires to translate the rotation to a translation.
    const InDetDD::TRT_BaseElement* strawElement =  m_trtman->getElement( ident );
    double strawLenthOver2 = 0.5* strawElement->strawLength();
 
    // In the global frame, 'dx1' corresponds to the readout side and 'dx2'
    // to the side closest the beampipe.
    double delta_dx_nearBeamPipe = dy + strawLenthOver2 * std::sin(rotx);
    double delta_dx_nearReadOut = dy - strawLenthOver2 * std::sin(rotx);
    
    // Uncertianty on straw positions (Arbitrary for now)
    double dxErr = 0.001;
    
    m_trtStrawAlignDbSvc->setDx(ident, delta_dx_nearReadOut, delta_dx_nearBeamPipe, dxErr);
  }
  
  ATH_MSG_DEBUG( "Leaving setAlignTransformL3 "   );
  return StatusCode::SUCCESS;
}

StatusCode TRT_AlignDbSvc::tweakAlignTransform(Identifier ident, Amg::Transform3D trans, unsigned int level) {
  
  ATH_MSG_DEBUG( "In tweakAlignTransform"  );
  
  if(level != 1 && level != 2 && level != 3){
    ATH_MSG_FATAL( "Incorrect call to tweakTrans"  );
    ATH_MSG_FATAL( "level must be 1,2 or 3"  );
    return StatusCode::FAILURE;
  }
  
  if(level == 1){
    if(tweakAlignTransformL1(ident,trans).isFailure()){
      
      ATH_MSG_WARNING( "tweak failed...just setting it."  );
      
      if(setAlignTransformL1(ident,trans).isFailure()){
    ATH_MSG_FATAL( "Set also failed!"  );
    ATH_MSG_FATAL( "Fail for real."  );
    return StatusCode::FAILURE;
      }
    }
  }
  
  else if(level == 2){
    if(tweakAlignTransformL2(ident,trans).isFailure()){
      
      ATH_MSG_WARNING( "tweak failed...just setting it."  );
      
      if(setAlignTransformL2(ident,trans).isFailure()){
    ATH_MSG_FATAL( "Set also failed!"  );
    ATH_MSG_FATAL( "Fail for real."  );
    return StatusCode::FAILURE;
      }
    }    
  }
  
  else if(level == 3){
    if(tweakAlignTransformL3(ident,trans).isFailure()){
      
      ATH_MSG_WARNING( "tweak failed...just setting it."  );
      
      if(setAlignTransformL3(ident,trans).isFailure()){
    ATH_MSG_FATAL( "Set also failed!"  );
    ATH_MSG_FATAL( "Fail for real."  );
    return StatusCode::FAILURE;
      }
    }    
  }
    
  ATH_MSG_DEBUG( "Leaving tweakAlignTransform"  );
  return StatusCode::SUCCESS;
}

StatusCode TRT_AlignDbSvc::tweakAlignTransformL1(Identifier ident, const Amg::Transform3D& trans) {
  ATH_MSG_DEBUG( "In tweakAlignTransformL1"  );
 
  // New additions for new global folder structure                                                                                                                    
  // No ATs exist for levels 1 & 2 --> need alternative                                                                                                               
  if (m_dynamicDB){
    if(!tweakGlobalFolder(ident, trans)) {
      ATH_MSG_ERROR( "Attempt tweak GlobalDB folder failed" );
      return StatusCode::FAILURE;
    }
  }
  else {
    
    AlignableTransform* pat=getTransPtr("/TRT/Align/TRT");
    if (!pat){
      ATH_MSG_WARNING( "The AlignableTransform for key " << "/TRT/Align/TRT" << " does not exist "  );
      ATH_MSG_WARNING( "Failing ... "  );
      return StatusCode::FAILURE;
    }
    
    if( !(m_trtid->is_trt(ident)) ){
      ATH_MSG_WARNING( "The identifier " << ident << " is not from the TRT "  );
      ATH_MSG_WARNING( "Failing ... "  );
      return StatusCode::FAILURE;
    }
    
    int bec=m_trtid->barrel_ec(ident);
    Identifier mid=m_trtid->module_id(bec,0,0);
    
    if( !(pat->tweak(mid,Amg::EigenTransformToCLHEP(trans))) ){
      ATH_MSG_WARNING( "The Alignable transfor tweek failed for " << mid  );
      ATH_MSG_WARNING( "Failing ... "  );
      return StatusCode::FAILURE;
    }
  }
 
  ATH_MSG_DEBUG( "Leaving tweakAlignTransformL1"  );
  return StatusCode::SUCCESS;
}

StatusCode TRT_AlignDbSvc::tweakAlignTransformL2(Identifier ident, const Amg::Transform3D& trans) {
  ATH_MSG_DEBUG( "In tweakAlignTransformL2"  );
  // multiply level 2 AlignableTransform for the module containing ident
  // by an additional transform.
 
  //does a folder exist corresponding to this identifier?
  std::string key=findkey(ident,m_alignString);
  AlignableTransform* pat=getTransPtr(key);
  if (!pat){
    ATH_MSG_WARNING( "The transfor for key: " << key << "Does not exist"  );
    ATH_MSG_WARNING( " Failing ... "  );
    return StatusCode::FAILURE;
  }
 
  //OK the key exists
  // make sure the identifier is a "module identifier"
  if( !(m_trtid->is_trt(ident)) ) return StatusCode::FAILURE;
  int bec=m_trtid->barrel_ec(ident);
  int layer=m_trtid->layer_or_wheel(ident);
  int sector=m_trtid->phi_module(ident);
  int strawLayer=m_trtid->straw_layer(ident);
  int ring = getRingForStrawLayer(strawLayer);
  Identifier mid;
  if(key == (m_alignDBprefix+"L2A") || key == (m_alignDBprefix+"L2C"))
    mid = m_trtid->layer_id(bec,0,layer,ring);
  else
    mid = m_trtid->module_id(bec,sector,layer);
 
  // multiply the additional transformation
  if( !(pat->tweak(mid,Amg::EigenTransformToCLHEP(trans))) ){
    ATH_MSG_WARNING( "Leaving tweakAlignTransformL2.  TWEAK FAILED!!"  );
    return StatusCode::FAILURE;
  }
 
  ATH_MSG_DEBUG( "Leaving tweakAlignTransformL2"  );
  return StatusCode::SUCCESS;
}

StatusCode TRT_AlignDbSvc::tweakAlignTransformL3 (Identifier ident, Amg::Transform3D trans) {
  ATH_MSG_DEBUG( "In tweakAlignTransformL3"   );
    
  // make sure the identifier is a "module identifier"
  if( !(m_trtid->is_trt(ident)) ) return StatusCode::FAILURE;
  int bec=m_trtid->barrel_ec(ident);
  
  //Only L3 in the barrel is currently supported 
  if(abs(bec) == 1){
    
    // Currently at L3 only displacements in the plane of the straw layer can be stored.
    // For tracks coming from the interaction point this coorisponds to the sensitive coordinate
    ATH_MSG_INFO("Storing L3 constants, Only dx and rotz will be written to DB");
    float dx = trans.translation().x();
 
    ATH_MSG_DEBUG("xx: " << trans.rotation()(0,0));    
    ATH_MSG_DEBUG("xy: " << trans.rotation()(0,1));    
    ATH_MSG_DEBUG("xz: " << trans.rotation()(0,2));    
    ATH_MSG_DEBUG("yx: " << trans.rotation()(1,0));    
    ATH_MSG_DEBUG("yy: " << trans.rotation()(1,1));    
    ATH_MSG_DEBUG("yz: " << trans.rotation()(1,2));    
    ATH_MSG_DEBUG("zx: " << trans.rotation()(2,0));    
    ATH_MSG_DEBUG("zy: " << trans.rotation()(2,1));    
    ATH_MSG_DEBUG("zz: " << trans.rotation()(2,2));    
    
    // For the definition of dx1 and dx2 see TRT_ConditionsData/StrawDx.h
    // Now we need to know the real meaning of dx1 and dx2. Alas, we
    // are in the frame of either of the two barrel sides, but not
    // both, so we need to know which one.
    bool sideA = m_trtid->barrel_ec(ident) == 1;
 
    float rotz = std::atan2(trans.rotation()(0,2),trans.rotation()(0,0));
 
    // Need the length of the wires to translate the rotation to a translation.
    const InDetDD::TRT_BaseElement* strawElement =  m_trtman->getElement( ident );
    double strawLenthOver2 = 0.5* strawElement->strawLength();
 
    // Old way (buggy)
    //double delta_dx_atLargerZ = dx + strawLenthOver2 * std::sin(rotz);
    //double delta_dx_atSmallerZ = dx - strawLenthOver2 * std::sin(rotz);
 
    // New way - The rotation is opposite for side A as compared to side C
    double delta_dx_atLargerZ;
    double delta_dx_atSmallerZ;
    if (sideA){
      delta_dx_atLargerZ = dx - strawLenthOver2 * std::sin(rotz);
      delta_dx_atSmallerZ = dx + strawLenthOver2 * std::sin(rotz);
    }
    else{
      delta_dx_atLargerZ = dx + strawLenthOver2 * std::sin(rotz);
      delta_dx_atSmallerZ = dx - strawLenthOver2 * std::sin(rotz);
    }
 
    // Straw position closest to the electronics 
    double delta_dx1 = sideA ? delta_dx_atLargerZ : delta_dx_atSmallerZ;
    
    // Straw position away from the electronics 
    double delta_dx2 = sideA ? delta_dx_atSmallerZ : delta_dx_atLargerZ;
 
    // Uncertianty on straw positions
    double dxErr = 0.001;
    
    // now we do some gymnastics: because the strawdisplacement is
    // binned, we would rather not just 'add' the current change,
    // but rather 'set' the total change. note that this goes
    // wrong if there were other straw corrections
    // applied. We might have to worry about that at
    // some point. 
    double dx1_current = m_trtStrawAlignDbSvc->getDx1(ident);
    double dx2_current = m_trtStrawAlignDbSvc->getDx2(ident) ;
    ATH_MSG_DEBUG("Current dx1 is: " << dx1_current);
    ATH_MSG_DEBUG("Current dx2 is: " << dx2_current);
    
    double dx1_new = dx1_current + delta_dx1;
    double dx2_new = dx2_current + delta_dx2;
 
    //ATH_MSG_INFO("New dx1 is: " << dx1_new);
    //ATH_MSG_INFO("New dx2 is: " << dx2_new);
 
    m_trtStrawAlignDbSvc->setDx(ident, dx1_new, dx2_new, dxErr);
  }else{
    
    // Currently at L3 only displacements in the plane of the straw layer can be stored.
    // For tracks coming from the interaction point this coorisponds to the sensitive coordinate
    ATH_MSG_INFO("Storing L3 Endcap constants, Only dy and rotz will be written to DB");
    
    // The expected transformation is in the local frame of the Endcap Straws
    // In the local straw frame:
    //   - the z-axis is along the straw and points toward the beampipe
    //   - the x-axis is along global-z and away from the interaction point 
    //          (locX = globZ A-side / locX = -1 *gobZ C-side) 
    //   - the y-axis is along global phi_hat direction determined by the other 2. 
    //          (clockwise C-side, counter clockwise A-Side)
    float dy = trans.translation().y();
    float rotx = std::atan2(trans.rotation()(1,2),trans.rotation()(2,2));
 
    // Need the length of the wires to translate the rotation to a translation.
    const InDetDD::TRT_BaseElement* strawElement =  m_trtman->getElement( ident );
    double strawLenthOver2 = 0.5* strawElement->strawLength();
 
    // In the global frame, 'dx1' corresponds to the readout side and 'dx2'
    // to the side closest the beampipe.
    double delta_dx_nearBeamPipe = dy + strawLenthOver2 * std::sin(rotx);
    double delta_dx_nearReadOut = dy - strawLenthOver2 * std::sin(rotx);
 
    // Uncertianty on straw positions
    double dxErr = 0.001;
 
    // now we do some gymnastics: because the strawdisplacement is
    // binned, we would rather not just 'add' the current change,
    // but rather 'set' the total change. note that this goes
    // wrong if there were other straw corrections
    // applied. We might have to worry about that at
    // some point. 
    double dx1_current = m_trtStrawAlignDbSvc->getDx1(ident);
    double dx2_current = m_trtStrawAlignDbSvc->getDx2(ident) ;
    ATH_MSG_DEBUG("Current dx1 is: " << dx1_current);
    ATH_MSG_DEBUG("Current dx2 is: " << dx2_current);
    
    double dx1_new = dx1_current + delta_dx_nearReadOut;
    double dx2_new = dx2_current + delta_dx_nearBeamPipe;
 
    m_trtStrawAlignDbSvc->setDx(ident, dx1_new, dx2_new, dxErr);
  }
  
  ATH_MSG_DEBUG( "Leaving tweakAlignTransformL3 "   );
  return StatusCode::SUCCESS;
}

std::string TRT_AlignDbSvc::findkey(const Identifier& ident, const std::string& type) const {
  std::ostringstream result;
  if (!m_trtid->is_trt(ident)) return result.str();
  
  if(type=="AL" || type=="ALold") {
    //    result << "/TRT/Align/";
    result << m_alignDBprefix;
  } else if(type=="DF") {
    result << "/TRT/Calib/DF";
  } else return result.str();
 
  int bec=m_trtid->barrel_ec(ident);
 
  if(type!="AL" && type!="ALold") result << "_";
 
  if(type=="AL"){
    if(bec==2) {
      result << "L2A";
    } else if(bec==-2) {
      result << "L2C";
    } else {
      result << "B";
      int layer = m_trtid->layer_or_wheel(ident);
      result << layer;
    }
  }else if(type=="ALold"){
    if(bec==2) {
      result << "A";
    } else if(bec==-2) {
      result << "C";
    } else {
      result << "B";
    }
    int layer = m_trtid->layer_or_wheel(ident);
    result << layer;
  }
  
  
  return result.str();
}     

std::string TRT_AlignDbSvc::prefixtag(const std::string& key) {
 
  std::ostringstream result;
  std::string  detector=key.substr(1,3);
 
  std::string type,quantity;
  int ind,ind1;
  if(detector=="TRT") {
    type=key.substr(5,5);
    ind=11;
    ind1=11;
    quantity=key.substr(ind,2);
  } else if(detector=="Indet") {
    type=key.substr(7,5);
    ind=17;
    ind1=13;
    quantity="";
    if(type=="Align" && key.substr(ind1,2)!="ID") {
      if(key.substr(ind1+3,1)=="_") quantity="DF";
    }
  } else {
    return "";
  }
  if(type=="Calib") {
 
    if(quantity=="DF") {
       result << "TrtDf";
    } else if (quantity=="RT") {
       result <<  "TrtRt";
    } else if (quantity=="T0") {
       result <<  "TrtT0";
    } else {
       return "TrtDe_";
    }
    std::string module=key.substr(ind+3,2);
    result << module;
  } else if(type=="Align") {
    std::string  module=key.substr(ind1,2);
    if(module=="ID" || module=="TR") {
      return "TrtAl_";
    } else {
      result << "TrtAl";
    }
    module=key.substr(ind+3,2);
    result << module;
  } else {
    return "";
  }
  result << "_";
  return result.str();
}

StatusCode TRT_AlignDbSvc::createAlignObjects() const{

  ATH_MSG_DEBUG( "createAlignObjects method called"  );

  if (m_detStore->contains<AlignableTransformContainer>(m_alignroot)) {
    ATH_MSG_WARNING( "createAlignObjects: AlignableTransformContainer with name " << m_alignroot
                     << " already exists."  );
    return StatusCode::FAILURE;
  }
  
  auto patc = std::make_unique<AlignableTransformContainer>();
 
  for (unsigned int i=0;i<m_alignobjs.size();++i) {
    patc->push_back(std::make_unique<AlignableTransform>(m_alignobjs[i]));
    patc->add(m_alignchans[i]);
    ATH_MSG_INFO( " added empty AlignableTransform for key " << m_alignobjs[i]  );
  }

  ATH_CHECK( m_detStore->record(std::move(patc),m_alignroot) );
 
  ATH_MSG_DEBUG( "Leaving createAlignObjects"  );
  return StatusCode::SUCCESS;
}

StatusCode TRT_AlignDbSvc::createAlignObjectsWhichDoNotAlreadyExist(){

  ATH_MSG_DEBUG( "createAlignObjectsWhichDoNotAlreadyExist method called"  );

  ATH_CHECK( m_detStore->contains<AlignableTransformContainer>(m_alignroot) );
  
  const AlignableTransformContainer* patc= getContainer();

  AlignableTransformContainer* patc_NonConst = const_cast<AlignableTransformContainer*>(patc);
  
  for (unsigned int i=0;i<m_alignobjs.size();++i) {
    
    if(m_alignString != "Alold" && !isOldKey(m_alignobjs[i] )){
      
      if(!cgetTransPtr(m_alignobjs[i])){
    
    AlignableTransform* pat=new AlignableTransform(m_alignobjs[i]);
    patc_NonConst->push_back(pat);
    patc_NonConst->add(m_alignchans[i]);
    ATH_MSG_INFO( " added empty AlignableTransform for key " << m_alignobjs[i]  );
 
      }
    }
  }
 
  ATH_MSG_DEBUG( "Leaving createAlignObjectsWhichDoNotAlreadyExist"  );
  return StatusCode::SUCCESS;
}

void TRT_AlignDbSvc::printCondObjects() const {
  ATH_MSG_DEBUG( " In printCondObjects "  );
  ATH_MSG_INFO( " TRT Conditions DB contains " << m_alignobjs.size() <<" AlignableTransforms: "  );
  
  const AlignableTransform *pat;
  std::vector<std::string>::const_iterator iobj=m_alignobjs.begin();
  std::vector<std::string>::const_iterator iobjE=m_alignobjs.end();
  for (;iobj!=iobjE;++iobj) {
    ATH_MSG_INFO( " " << *iobj  );
    pat=cgetTransPtr(*iobj);
    if (pat) {
      pat->print2();
    } else {
      ATH_MSG_ERROR( " Could not find key: " << *iobj  );
    }
  }
 
  ATH_MSG_DEBUG( " Leaving printCondObjects"  );
  return;
}

AlignableTransform* TRT_AlignDbSvc::getTransPtr(const std::string& key) const {
  ATH_MSG_DEBUG( "In (and leaving) getTransPtr"  );
  return const_cast<AlignableTransform*>(cgetTransPtr(key)) ;
}

const AlignableTransform* TRT_AlignDbSvc::cgetTransPtr(const std::string& key) const {
  ATH_MSG_DEBUG( "In cgetTransPtr"  );
  ATH_MSG_DEBUG( "Getting the poiter for key " << key  );
 
  // Retrieve AlignableTransform pointer for a given key - const version
  const AlignableTransform* pat=nullptr;
  const AlignableTransformContainer* patc= getContainer();
  // the retrieve is unnecessary. in fact, if patc==0, retrieve should
  // fail as well.
 
  if( patc || StatusCode::SUCCESS==m_detStore->retrieve(patc,m_alignroot)) {
    for(DataVector<AlignableTransform>::const_iterator dva=patc->begin();dva!=patc->end();++dva) {
      if( (*dva)->tag()==key ) {
    pat = *dva;
    break;
      }
    }
  }
  
  if(!pat){
    ATH_MSG_WARNING( "AlignableTransform poiter is NULL!!!!!"  );
    ATH_MSG_WARNING( "Failed to get the AlignableTransform for key " << key  );
  }
  
  ATH_MSG_DEBUG( "leaving cgetTransPtr "  );
  return pat;
}

int TRT_AlignDbSvc::getRingForStrawLayer(int strawlayer) {
  return strawlayer / 4 * 4;
}
 

bool TRT_AlignDbSvc::isOldKey(const std::string& input) const{
  
  for(unsigned int i=0; i<14; ++i){
    std::string testA = "A"+intToString(i);
    ATH_MSG_VERBOSE( " The testA is " << testA << " " 
                     << bool(input.find(testA)!=std::string::npos)  );
    if(input.find(testA)!=std::string::npos){
      return true;
    }
 
    std::string testC = "C"+intToString(i);
    ATH_MSG_VERBOSE( " The testC is " << testC << " " 
                     << bool(input.find(testC)!=std::string::npos)  );
    if(input.find(testC)!=std::string::npos){
      return true;
    }
  
  }
  
  return false;
}

std::string TRT_AlignDbSvc::intToString(int input) {
  std::ostringstream stm;
  stm << input;
  return stm.str();
}

void TRT_AlignDbSvc::printTransform(const std::string& thisName,  const Amg::Transform3D& transform) { 
  std::cout << thisName << " " << transform(0,3) << " " << transform(1,3) << "  " << transform(2,3) << std::endl;
  std::cout << thisName << " " << transform(0,0) << " " << transform(0,1) << "  " << transform(0,2) << std::endl;
  std::cout << thisName << " " << transform(1,0) << " " << transform(1,1) << "  " << transform(1,2) << std::endl;
  std::cout << thisName << " " << transform(2,1) << " " << transform(2,1) << "  " << transform(2,2) << std::endl;
  return;
}
 
// write global DB to file
StatusCode TRT_AlignDbSvc::writeGlobalFolderFile( const std::string &file)
  const {
 
  if (m_dynamicDB){
    ATH_MSG_DEBUG( "writeFile: Write GlobalFolder DB in text file: " << file );
    std::ofstream outfile(file);
    std::vector<std::string> folder_list = {"/TRT/AlignL1/TRT"};
 
    for (std::vector<std::string>::iterator it = folder_list.begin(); it != folder_list.end(); ++it){
 
      outfile << *it << std::endl;
      const CondAttrListCollection* atrlistcol=nullptr;
      if (StatusCode::SUCCESS==m_detStore->retrieve(atrlistcol,*it)) {
        // loop over objects in collection
        for (CondAttrListCollection::const_iterator citr=atrlistcol->begin(); citr!=atrlistcol->end();++citr) {
 
          const coral::AttributeList& atrlist=citr->second;
          outfile  << atrlist["bec"].data<int>()
                    << " "     << atrlist["layer"].data<int>()
                    << " "     << atrlist["sector"].data<int>()
                    << " "     << atrlist["Tx"].data<float>()
                    << " "     << atrlist["Ty"].data<float>()
                    << " "     << atrlist["Tz"].data<float>()
                    << " "     << atrlist["phi"].data<float>()
                    << " "     << atrlist["theta"].data<float>()
                    << " "     << atrlist["psi"].data<float>() << std::endl;
        }
      }
      else {
        if (msgLvl(MSG::INFO)){
          ATH_MSG_INFO( "Cannot find " << *it << " Container - cannot write DB in text file "  );
      return StatusCode::FAILURE;
    }
      }
    }
 
    return StatusCode::SUCCESS;
  }
  else {
    ATH_MSG_DEBUG( "writeFile: No dynamic Run2 DB structure is present --> skipping writing file " << file );
    return StatusCode::SUCCESS;
  }
}
 
 
 
 
 
 
StatusCode TRT_AlignDbSvc::tweakGlobalFolder(Identifier ident, const Amg::Transform3D& trans ) {
 
  // find transform key, then set appropriate transform
  const CondAttrListCollection* atrlistcol1=nullptr;
  CondAttrListCollection* atrlistcol2=nullptr;
  bool result = false;
  std::string key="/TRT/AlignL1/TRT";
  ATH_MSG_DEBUG( " Identifier is valid: "<< ident.is_valid()  );
  int bec=m_trtid->barrel_ec(ident);
  const unsigned int DBident=1000+bec*100;
  // so far not a very fancy DB identifier, but seems elaborate enough for this simple structure
 
  if (StatusCode::SUCCESS==m_detStore->retrieve(atrlistcol1,key)) {
    // loop over objects in collection
    atrlistcol2 = const_cast<CondAttrListCollection*>(atrlistcol1);
    if (atrlistcol1!=nullptr){
      for (CondAttrListCollection::const_iterator citr=atrlistcol2->begin(); citr!=atrlistcol2->end();++citr) {
 
        const coral::AttributeList& atrlist=citr->second;
    coral::AttributeList& atrlist2  = const_cast<coral::AttributeList&>(atrlist);
 
        if(citr->first!=DBident){
          ATH_MSG_DEBUG( "tweakGlobalFolder fails due to identifier mismatch"  );
          continue;
    }
        else {
          ATH_MSG_DEBUG( "Tweak Old global DB -- channel: " << citr->first
                         << " ,bec: "    << atrlist2["bec"].data<int>()
                         << " ,layer: "  << atrlist2["layer"].data<int>()
                         << " ,sector: " << atrlist2["sector"].data<int>()
                         << " ,Tx: "     << atrlist2["Tx"].data<float>()
                         << " ,Ty: "     << atrlist2["Ty"].data<float>()
                         << " ,Tz: "     << atrlist2["Tz"].data<float>()
                         << " ,phi: "    << atrlist2["phi"].data<float>()
                         << " ,theta: "  << atrlist2["theta"].data<float>()
                         << " ,psi: "    << atrlist2["psi"].data<float>()  );
 
 
          // Follow same definitions as in TRT_AlignDbSvc.cxx                                                                                               
      CLHEP::Hep3Vector  oldtranslation(atrlist2["Tx"].data<float>(),atrlist2["Ty"].data<float>(),atrlist2["Tz"].data<float>());
      CLHEP::HepRotation oldrotation;
      oldrotation.set(atrlist["phi"].data<float>(),atrlist["theta"].data<float>(),atrlist["psi"].data<float>());
      HepGeom::Transform3D oldtransform(oldrotation, oldtranslation);
          
          // get the new transform
      HepGeom::Transform3D newtrans = Amg::EigenTransformToCLHEP(trans)*oldtransform;
      Amg::Transform3D newtransAMG = trans*Amg::CLHEPTransformToEigen(oldtransform);
 
          // Extract the values we need to write to DB
      Amg::Vector3D shift=newtransAMG.translation();
      CLHEP::HepRotation rot=newtrans.getRotation();
      Amg::Vector3D eulerangles;
          eulerangles[0] = rot.getPhi();
          eulerangles[1] = rot.getTheta();
          eulerangles[2] = rot.getPsi();      
 
          atrlist2["Tx"].data<float>() = shift.x();
          atrlist2["Ty"].data<float>() = shift.y();
          atrlist2["Tz"].data<float>() = shift.z();
          atrlist2["phi"].data<float>()   = eulerangles[0] ;
          atrlist2["theta"].data<float>() = eulerangles[1] ;
          atrlist2["psi"].data<float>()   = eulerangles[2] ;
 
      result = true;
      ATH_MSG_DEBUG( "Tweak New global DB -- channel: " << citr->first
                         << " ,bec: "    << atrlist2["bec"].data<int>()
                         << " ,layer: "  << atrlist2["layer"].data<int>()
                         << " ,sector: " << atrlist2["sector"].data<int>()
                         << " ,Tx: "     << atrlist2["Tx"].data<float>()
                         << " ,Ty: "     << atrlist2["Ty"].data<float>()
                         << " ,Tz: "     << atrlist2["Tz"].data<float>()
                         << " ,phi: "    << atrlist2["phi"].data<float>()
                         << " ,theta: "  << atrlist2["theta"].data<float>()
                         << " ,psi: "    << atrlist2["psi"].data<float>()  );
 
        }
      }
    }
    else {
      ATH_MSG_ERROR("tweakGlobalFolder: cast fails for DBident " << DBident );
      return StatusCode::FAILURE;
    }
  }
  else {
    ATH_MSG_ERROR("tweakGlobalFolder: cannot retrieve CondAttrListCollection for key " << key );
    return StatusCode::FAILURE;
  }
 
  if (result) return StatusCode::SUCCESS;
  else return StatusCode::FAILURE;
}
 
 
const AlignableTransformContainer* TRT_AlignDbSvc::getContainer() const
{
  const AlignableTransformContainer* ptr = nullptr;
  if (m_detStore->retrieve (ptr, m_alignroot).isFailure()) {
    ATH_MSG_ERROR ("Cannot retrieve " << m_alignroot << " from detStore ");
  }
  return ptr;
}