/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/MagneticField/MagneticField-00-01-38/src/ConnectionDB.cxx

Go to the documentation of this file.

#ifndef BEAN
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/SvcFactory.h"
#include "GaudiKernel/ISvcLocator.h"
#include "GaudiKernel/SmartDataPtr.h"
#include "GaudiKernel/Bootstrap.h"

#include "facilities/Util.h"
#else
#include "DatabaseSvc/DatabaseSvc.h"
#endif

#include <iostream>
#include <cstdio>
#include <cstdlib>

#include "MagneticField/ConnectionDB.h"

using namespace std;

namespace FieldDBUtil {

#ifndef BEAN
  ConnectionDB::ConnectionDB() {
    StatusCode sc = Gaudi::svcLocator()->service("DatabaseSvc", m_dbsvc, true);
    if (sc .isFailure() ) {
      std::cout << "ERROR: In ConnectionDB()--> Unable to find DatabaseSvc " << std::endl;
    }
  }
#endif

  ConnectionDB::eRet ConnectionDB::getReadSC_MagnetInfo( std::vector<double>& current, int runNo) {

#ifdef BEAN
    static DatabaseSvc* m_dbsvc = 0;          
                                         
    if( !m_dbsvc ) {                          
      m_dbsvc = DatabaseSvc::instance();      
    }                                         
#endif

    //Read magnetic field map                                                                                             
    char stmt1[200];
    int run_No = std::abs(runNo);                                                                                             

    sprintf(stmt1,"select Magnet_Current,SCQL,SCQR from SC_magnet where run_number = %d ",run_No);

    DatabaseRecordVector results;                                                                                         
    results.clear();                                                                                                      
    int status = m_dbsvc->query("run",stmt1,results);                                                                      

    if(status<0){                                                                                                         
      std::cout << "ERROR Read the SSM and SCQ current from the Database" << endl;                                        
      exit(1);                                                                                                            
    }                                                                                                                     
                                                                                                                          
    int RowNumber = results.size();                                                                                       
    if(RowNumber!=1){                                                                                                     
      std::cout<<"ERROR:error searching SC_Magnet Data in the database, check your selection criterions"<<std::endl;      
      return RETMySQLError;                                                                                               
    }                                                                                                                     

    DatabaseRecord& rec = *results[0];                                                                                    
    double ssm_curr = rec.GetDouble("Magnet_Current");                                                                    
    double scql_curr = rec.GetDouble("SCQL");                                                                             
    double scqr_curr = rec.GetDouble("SCQR");                                                                             

    // save results in vector                                                                                             
    current.resize(3);                                                                                                    
    current[0] = ssm_curr;                                                                                                
    current[1] = scql_curr;                                                                                               
    current[2] = scqr_curr;                                                                                               

    return RETOk;                                                                                                         
  }

  ConnectionDB::eRet ConnectionDB::getBeamEnergy( std::vector<double>& beamE, int runNo) {
    //Read magnetic field map
     char stmt1[200];
     int run_No = std::abs(runNo);

     sprintf(stmt1,"select BPR_PRB,BER_PRB from RunParams where run_number = %d ",run_No);

     DatabaseRecordVector results;
     results.clear();
     int status = m_dbsvc->query("run",stmt1,results);
     if(status<0){
       std::cout << "ERROR Read the beam energy from the Database" << endreq;
       exit(1);
     }

     int RowNumber = results.size();

     if(RowNumber == 0) {
       beamE.push_back(1.843); // for positron
       beamE.push_back(1.843); // for electron

       return RETOk;
     }

     if(RowNumber!=1){
         std::cout<<"ERROR:error searching beam energy in the database, check your selection criterions"<<std::endl;
        return RETMySQLError;
     }

     beamE.push_back(atof((*results[0])["BPR_PRB"])); // for positron
     beamE.push_back(atof((*results[0])["BER_PRB"])); // for electron

     return RETOk;
  }                                                                                                                       

}

---