UseDedx.cxx

Go to the documentation of this file.

//$Header: /bes/bes/BossCvs/Calibration/CalibSvc/CalibROOTCnv/src/test/UseDedx.cxx,v 1.7
//2010/03/05 02:17:54 huangb Exp $
/*#include <stdio.h>
#include "CalibData/CalibModel.h"
#include "CalibData/Dedx/DedxCalibData.h"
#include "CalibDataSvc/ICalibRootSvc.h"
#include "GaudiKernel/AlgFactory.h"
#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/DataSvc.h"
#include "GaudiKernel/IDataProviderSvc.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/Service.h"
#include "GaudiKernel/SmartDataPtr.h"
using namespace std;
/// Simple algorithm to test functioning of "the other" TDS
class UseDedx : public Algorithm {
 
public:
  UseDedx(const std::string& name, ISvcLocator* pSvcLocator);
 
  StatusCode initialize();
 
  StatusCode execute();
 
  StatusCode finalize();
 
private:
  IDataProviderSvc* m_pCalibDataSvc;
  ICalibRootSvc*    m_pRootSvc;
  // Maybe something to say which kind of data to look up?
 
};
 
/// Instantiation of a static factory to create instances of this algorithm
//static const AlgFactory<UseDedx> Factory;
//const IAlgFactory& UseCalibFactory = Factory;
//const IAlgFactory& UseDedxFactory = Factory;
*/
#include "UseDedx.h"
UseDedx::UseDedx( const std::string& name, ISvcLocator* pSvcLocator )
    : Algorithm( name, pSvcLocator ), m_pCalibDataSvc( 0 ) {
  // Declare properties here.
}
 
StatusCode UseDedx::initialize() {
  StatusCode sc;
  MsgStream  log( msgSvc(), name() );
  log << MSG::INFO << "Initialize()" << endmsg;
 
  // So far don't have any properties, but in case we do some day..
  //  setProperties();
 
  log << MSG::INFO << "setProperties()" << endmsg;
 
  sc = service( "CalibDataSvc", m_pCalibDataSvc, true );
 
  if ( !sc.isSuccess() )
  {
    log << MSG::ERROR << "Could not get IDataProviderSvc interface of CalibXmlCnvSvc"
        << endmsg;
    return sc;
  }
  else
  { log << MSG::DEBUG << "Retrieved IDataProviderSvc interface of CalibXmlCnvSvc" << endmsg; }
 
  sc = service( "CalibRootCnvSvc", m_pRootSvc, true );
  if ( !sc.isSuccess() )
  {
    log << MSG::ERROR << "Could not get ICalibRootSvc interface of CalibRootCnvSvc" << endmsg;
    return sc;
  }
 
  // Get properties from the JobOptionsSvc
 
  sc = setProperties();
 
  return StatusCode::SUCCESS;
}
 
StatusCode UseDedx::execute() {
 
  MsgStream log( msgSvc(), name() );
 
  std::string fullPath = "/Calib/DedxCal";
 
  log << MSG::INFO << "execute() fullPath = " << fullPath << endmsg;
 
  SmartDataPtr<CalibData::DedxCalibData> test( m_pCalibDataSvc, fullPath );
  int                                    i;
  int                                    j;
  double                                 ddg;
  double                                 ggs;
  double                                 wireg;
  double                                 gain;
  double                                 zdep;
  double                                 resol;
  double                                 layerg;
  double                                 m_id_doca;
  double                                 m_iner_chi;
  double                                 m_iner_gain;
  double                                 m_iner_hits;
  double                                 m_ip_eangle;
  double                                 m_out_chi;
  double                                 m_out_gain;
  double                                 m_out_hits;
 
  for ( i = 0; i < 4; i++ )
  {
    for ( j = 0; j < 43; j++ )
    {
      ddg  = test->getddg( i, j );
      ggs  = test->getggs( i, j );
      zdep = test->getzdep( i, j );
      //    std::cout<<"ddg["<<i<<"]["<<j<<"]="<<ddg;
      //        std::cout<<"         ggs["<<i<<"]["<<j<<"]="<<ggs;
      //        std::cout<<"        zdep["<<i<<"]["<<j<<"]="<<zdep;
      //    std::cout<<"\n";
    }
  }
 
  for ( int i = 0; i < 1600; i++ )
  {
    m_id_doca   = test->get_id_doca( i );
    m_iner_chi  = test->get_iner_chi( i );
    m_iner_gain = test->get_iner_gain( i );
    m_iner_hits = test->get_iner_hits( i );
    m_ip_eangle = test->get_ip_eangle( i );
    m_out_chi   = test->get_out_chi( i );
    m_out_gain  = test->get_out_gain( i );
    m_out_hits  = test->get_out_hits( i );
    std::cout << "m_id_doca:  " << m_id_doca << " m_iner_chi: " << m_iner_chi
              << " m_iner_gain: " << m_iner_gain << " m_iner_hits:" << m_iner_hits
              << " m_ip_eangle:" << m_ip_eangle << "  m_out_chi=" << m_out_chi
              << " m_out_gain:" << m_out_gain << " m_out_hits=" << m_out_hits << std::endl;
  }
  gain = test->getgain();
  //    std::cout<<"gain="<<gain<<"\n";
  double m_costheta;
  for ( int i = 0; i < 80; i++ )
  {
    m_costheta = test->get_costheta( i );
    std::cout << " m_costheta=" << m_costheta;
    if ( i % 5 == 0 ) std::cout << "\n" << std::endl;
  }
  double t0, dedx;
  for ( int i = 0; i < 35; i++ )
  {
    t0   = test->get_t0( i );
    dedx = test->get_dedx( i );
    std::cout << "t0  =" << t0 << "  dedx=" << dedx;
    if ( i % 5 == 0 ) std::cout << "\n" << std::endl;
  }
 
  resol = test->getresol();
  //  std::cout<<"resol="<<resol<<"\n";
 
  for ( i = 0; i < 43; i++ )
  {
    layerg = test->getlayerg( i );
    //   std::cout<<"layerg["<<i<<"]="<<layerg;
  }
 
  for ( i = 0; i < 60; i++ )
  {
    wireg = test->getwireg( i );
    // std::cout<<"wireg["<<i<<"]="<<wireg<<"\n";
  }
 
  for ( i = 0; i < 20; i++ ) { std::cout << "hadron " << i << " = " << test->get_hadron( i ); }
 
  // for pah is error
  // m_pRootSvc->writeToRoot("/home/huangb/2.root", fullPath);
  return StatusCode::SUCCESS;
}
 
StatusCode UseDedx::finalize() {
 
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "        UseDedx FINALIZE!! " << endmsg;
 
  return StatusCode::SUCCESS;
}