/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Calibration/CalibSvc/CalibROOTCnv/CalibROOTCnv-00-01-13/src/test/UseMuc.cxx

Go to the documentation of this file.

//$Header: /bes/bes/BossCvs/Calibration/CalibSvc/CalibROOTCnv/src/test/UseMuc.cxx,v 1.2 2008/04/09 01:02:04 huangb Exp $
/*#include <stdio.h>
#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/AlgFactory.h"
#include "GaudiKernel/IDataProviderSvc.h"
#include "GaudiKernel/Service.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/SmartDataPtr.h"
#include "CalibData/Muc/MucCalibData.h"
#include "CalibData/CalibModel.h"
#include "GaudiKernel/DataSvc.h"
#include "CalibDataSvc/ICalibRootSvc.h"
#include "GaudiKernel/MsgStream.h"
#include "CalibData/CalibModel.h"
#include "CalibData/CalibBase1.h"
#include "CalibData/Muc/MucCalibConst.h"
//#include "CalibROOTCnv/MucIdTransform.h"
//#include "MucCalibAlg/MucIdTransform.h"


//using namespace std;
class UseMuc : public Algorithm {

public:
  UseMuc(const std::string& name, ISvcLocator* pSvcLocator); 

  StatusCode initialize();

  StatusCode execute();

  StatusCode finalize();

private:

  IDataProviderSvc* m_pCalibDataSvc;
  ICalibRootSvc*    m_pRootSvc;
  // MucIdTransform*   m_pIdTr;
  // Maybe something to say which kind of data to look up?

};

//static const AlgFactory<UseMuc> Factory;
//const IAlgFactory& UseCalibFactory = Factory;
//const IAlgFactory& UseMucFactory = Factory;
*/
#include "UseMuc.h"
UseMuc::UseMuc( const std::string&  name, 
                    ISvcLocator*        pSvcLocator )
  : Algorithm     ( name, pSvcLocator ), m_pCalibDataSvc(0)
{
  // Declare properties here.

}


StatusCode UseMuc::initialize() {
  StatusCode sc;
  MsgStream log(msgSvc(), name());
  log << MSG::INFO << "Initialize()" << endreq;

  // So far don't have any properties, but in case we do some day..
//  setProperties();

   log<<MSG::INFO << "setProperties()" << endreq;
  
   sc = service("CalibDataSvc", m_pCalibDataSvc, true);

  if ( !sc.isSuccess() ) {
    log << MSG::ERROR 
        << "Could not get IDataProviderSvc interface of CalibXmlCnvSvc" 
        << endreq;
    return sc;
  } else {
    log << MSG::DEBUG 
        << "Retrieved IDataProviderSvc interface of CalibXmlCnvSvc" 
        << endreq;
  }
  
  sc = service("CalibRootCnvSvc", m_pRootSvc, true);
  if ( !sc.isSuccess() ) {
    log << MSG::ERROR 
        << "Could not get ICalibRootSvc interface of CalibRootCnvSvc" 
        << endreq;
    return sc;
  }
  // Get properties from the JobOptionsSvc
  
 sc = setProperties();

 // m_pIdTr = new MucIdTransform();

 return StatusCode::SUCCESS;

}


StatusCode UseMuc::execute( ) {
   
  MsgStream log(msgSvc(), name());
/*  std::vector<double> emccalib;
  CalibData::MucCalibData* emc = new CalibData::MucCalibData;
  //set emcCalibdata--------------------------------
   for(i=0;i<6000;i++)
    {  emccalib.push_back(i);
    }
  
   emc->setDigiCalibConst(&emccalib);
   emccalib.clear();
  //set EnCoeff
   for(i=0;i<43;i++)
    { emccalib.push_back(1);
    }
   emc->setEnCoeff(&emccalib);
   emccalib.clear();  
   //set PosCoeff
  for(i=0;i<15;i++)
    { emccalib.push_back(0);
    }
  emc->setPosCoeff(&emccalib);
  emc->setrunfrm(0);
  emc->setrunto(0);
 */
   //register the emcCalibData into the TCDS ; nelectable
   std::string fullPath = "/Calib/MucCal";
    log << MSG::INFO<<"execute() fullPath = "<<fullPath<< endreq;
  //  m_pCalibDataSvc->registerObject(fullPath,emc);

    //write the MucCalibData to the rootfile; 
//   m_pRootSvc->writeToRoot("./Muc_dat.root", fullPath);
   //another way,write the MucCalibData without TCDS
//   m_pRootSvc->writeToRoot("./Muc_dat1.root", emc);
  //retrieve DataObject from the TCDS-----------------
  SmartDataPtr<CalibData::MucCalibData> test(m_pCalibDataSvc, fullPath);
  if(!test);
   std::cout<<"m_uniformEff="<<test->getUniformEff()<<std::endl;
   std::cout<<"m_UniformNos="<<test->getUniformNos()<<std::endl;
   std::cout<<"m_UniformClst="<<test->getUniformClst()<<std::endl;
   for(int i=0;i<LAYER_MAX;i++){
     std::cout<<"m_layerEff["<<i<<"]="<< test->getLayerEff(i)<<"  m_layerNos["<<i<<"]="<<test->getLayerNos(i)<<std::endl;
     for(int j=0;j<CLST_MAX;j++){
      std::cout<<"m_layerClstPro["<<i<<"]["<<j<<"]="<<test->getLayerClstPro(i,j)<<std::endl;
     }
   }

   for(int i=0; i<PART_MAX; i++)
        {
          for(int j=0; j<((i==1)?B_SEG_NUM:E_SEG_NUM); j++)
          {
            for(int k=0; k<((i==1)?B_LAY_NUM:E_LAY_NUM); k++)
            {
                    std::cout << "Box: " << i << "\t" << j << "\t" << k << "\t"
                    << test->getBoxEff(i, j, k) << endl;
            }
          }     
        }
/*
        int strip_max;
        for(int i=0; i<PART_MAX; i++)
        {
          for(int j=0; j<((i==1)?B_SEG_NUM:E_SEG_NUM); j++)
          {
            for(int k=0; k<((i==1)?B_LAY_NUM:E_LAY_NUM); k++)
            {
                strip_max = m_pIdTr->GetStripMax(i, j, k);
                for(int n=0; n<strip_max; n++)
                {
                        std::cout << "Strip: " << i << "\t" << j << "\t" << k << "\t" 
                            << "\t" << n << "\t" << test->getStripEff(i, j, k, n) << endl;
                }
            }
          }                     
        }
*/

      return StatusCode::SUCCESS;
}

StatusCode UseMuc::finalize( ) {

  MsgStream log(msgSvc(), name());
  log << MSG::INFO 
      << "        UseMuc FINALIZE!! "
      << endreq;
  
  return StatusCode::SUCCESS;
}

---