/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Mdc/MdcAlignAlg/MdcAlignAlg-00-01-04/src/MdcAlignPar.cxx

Go to the documentation of this file.

#include "MdcAlignAlg/MdcAlignPar.h"

#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/IMessageSvc.h"
#include "GaudiKernel/StatusCode.h"
#include "GaudiKernel/ISvcLocator.h"
#include "GaudiKernel/Bootstrap.h"
#include "GaudiKernel/SmartDataPtr.h"
#include "GaudiKernel/IDataProviderSvc.h"
#include "GaudiKernel/PropertyMgr.h"

#include <iostream>
#include <fstream>
#include <string>
#include <cstring>

using namespace std;
using namespace Alignment;

MdcAlignPar::MdcAlignPar(){
     for(int iEP=0; iEP<NEP; iEP++){
          m_delDx[iEP] = 0.0;
          m_delDy[iEP] = 0.0;
          m_delDz[iEP] = 0.0;
          m_delRx[iEP] = 0.0;
          m_delRy[iEP] = 0.0;
          m_delRz[iEP] = 0.0;
     }
}

MdcAlignPar::~MdcAlignPar(){
}
void MdcAlignPar::initAlignPar(){
     for(int iEP=0; iEP<NEP; iEP++){
          m_dx[iEP] = 0.0;
          m_dy[iEP] = 0.0;
          m_dz[iEP] = 0.0;
          m_rx[iEP] = 0.0;
          m_ry[iEP] = 0.0;
          m_rz[iEP] = 0.0;

          m_errDx[iEP] = 0.0;
          m_errDy[iEP] = 0.0;
          m_errDz[iEP] = 0.0;
          m_errRx[iEP] = 0.0;
          m_errRy[iEP] = 0.0;
          m_errRz[iEP] = 0.0;
     }
}

bool MdcAlignPar::rdAlignPar(std::string alignFile){
     IMessageSvc* msgSvc;
     Gaudi::svcLocator() -> service("MessageSvc", msgSvc);
     MsgStream log(msgSvc, "MdcAlignPar");
     log << MSG::WARNING << "read MdcAlignPar data directly from run directory"
         << endreq;

     ifstream fpar(alignFile.c_str());
     log << MSG::INFO << "open file" << endreq;
     if( ! fpar.is_open() ){
          log << MSG::FATAL << "can not open alignment file " << alignFile.c_str() << endreq;
          return false;
     }

     int i;
     string strtmp;
     for(i=0; i<7; i++) fpar >> strtmp;
     for(int iEP=0; iEP<NEP; iEP++){
          fpar >> strtmp >> m_dx[iEP] >> m_dy[iEP] >> m_dz[iEP]
               >> m_rx[iEP] >> m_ry[iEP] >> m_rz[iEP];
     }
     fpar.close();
     return true;
}

void MdcAlignPar::wrtAlignPar(){
     IMessageSvc* msgSvc;
     Gaudi::svcLocator() -> service("MessageSvc", msgSvc);
     MsgStream log(msgSvc, "MdcAlignPar");
     log << MSG::INFO << "MdcAlignPar::wrtAlignPar()" << endreq;

     int i;
     int iEnd;
     string str[] = {"Inner_east", "Step0_east", "Step1_east",
                     "Step2_east", "Step3_east", "Step4_east",
                     "Step5_east", "Outer_east", "Inner_west",
                     "Step0_west", "Step1_west", "Step2_west",
                     "Step3_west", "Step4_west", "Step5_west",
                     "Outer_west"};

     string name[] = {"Elements", "DeltaX(mm)", "DeltaY(mm)",
                      "DeltaZ(mm)", "RX(rad)", "RY(rad)", "RZ(rad)"};

     ofstream fout("alignPar_new.txt");
     fout << setw(14) << name[0];
     for(i=1; i<7; i++) fout << setw(13) << name[i];
     fout << endl;
     for (iEnd=0; iEnd<NEP; iEnd++){
          fout << setw(14) << str[iEnd]
               << setw(13) << m_dx[iEnd] + m_delDx[iEnd]
               << setw(13) << m_dy[iEnd] + m_delDy[iEnd]
               << setw(13) << m_dz[iEnd] + m_delDz[iEnd]
               << setw(13) << m_rx[iEnd] + m_delRx[iEnd]
               << setw(13) << m_ry[iEnd] + m_delRy[iEnd]
               << setw(13) << m_rz[iEnd] + m_delRz[iEnd] << endl;
          if(7 == iEnd) fout << endl;
     }
     fout.close();

     ofstream fdel("delAlign_new.txt");
     fdel << setw(14) << name[0];
     for(i=1; i<7; i++) fdel << setw(13) << name[i];
     fdel << endl;
     for (iEnd=0; iEnd<NEP; iEnd++){
          fdel << setw(14) << str[iEnd] << setw(13) << m_delDx[iEnd]
               << setw(13) << m_delDy[iEnd] << setw(13) << m_delDz[iEnd]
               << setw(13) << m_delRx[iEnd] << setw(13) << m_delRy[iEnd]
               << setw(13) << m_delRz[iEnd] << endl;
          if(7 == iEnd) fdel << endl;
     }

     fdel << endl << "Fit error:" << endl;
     for (iEnd=0; iEnd<NEP; iEnd++){
          fdel << setw(14) << str[iEnd] << setw(13) << m_errDx[iEnd]
               << setw(13) << m_errDy[iEnd] << setw(13) << m_errDz[iEnd]
               << setw(13) << m_errRx[iEnd] << setw(13) << m_errRy[iEnd]
               << setw(13) << m_errRz[iEnd] << endl;
          if(7 == iEnd) fdel << endl;
     }
     fdel.close();
}

---