checkEstTof.cxx

Go to the documentation of this file.

//$Header: /bes/bes/BossCvs/Calibration/CalibSvc/CalibTreeCnv/src/test/checkEstTof.cxx,v 1.2
//2010/10/14 05:33:40 huangb Exp $
/*#include <stdio.h>
#include "CalibData/CalibModel.h"
#include "CalibData/Tof/TofCalibData.h"
#include "CalibData/Tof/bTofCalibBase.h"
#include "CalibDataSvc/ICalibTreeSvc.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 checkTof : public Algorithm {
 
public:
  checkTof(const std::string& name, ISvcLocator* pSvcLocator);
 
  StatusCode initialize();
 
  StatusCode execute();
 
  StatusCode finalize();
 
private:
  IDataProviderSvc* m_pCalibDataSvc;
  ICalibTreeSvc*    m_pTreeSvc;
  // 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<checkTof> Factory;
//const IAlgFactory& UseCalibFactory = Factory;
//const IAlgFactory& checkTofFactory = Factory;
*/
#include "checkEstTof.h"
checkEstTof::checkEstTof( const std::string& name, ISvcLocator* pSvcLocator )
    : Algorithm( name, pSvcLocator ), m_pCalibDataSvc( 0 ) {
  // Declare properties here.
}
 
StatusCode checkEstTof::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( "CalibTreeCnvSvc", m_pTreeSvc, true );
  if ( !sc.isSuccess() )
  {
    log << MSG::ERROR << "Could not get ICalibTreeSvc interface of CalibTreeCnvSvc" << endmsg;
    return sc;
  }
  // Get properties from the JobOptionsSvc
 
  sc = setProperties();
 
  return StatusCode::SUCCESS;
}
 
StatusCode checkEstTof::execute() {
 
  MsgStream log( msgSvc(), name() );
 
  // Cheat for now since Windows is having trouble finding definition
  // of Calibdata::Test_t
  std::string fullPath = "/Calib/EstTofCal";
  log << MSG::INFO << "execute() fullPath = " << fullPath << endmsg;
 
  SmartDataPtr<CalibData::TofCalibData> btof( m_pCalibDataSvc, fullPath );
 
  for ( int i = 0; i < 176; i++ )
  {
    double bAtten0 = btof->getBTofAtten( i, 0 );
    double bAtten1 = btof->getBTofAtten( i, 1 );
    double bP0     = btof->getBTofPleft( i, 0 );
    double bP1     = btof->getBTofPleft( i, 1 );
    double bP2     = btof->getBTofPleft( i, 2 );
    double bP3     = btof->getBTofPleft( i, 3 );
    double bP4     = btof->getBTofPleft( i, 4 );
    double bP5     = btof->getBTofPleft( i, 5 );
    double bP10    = btof->getBTofPright( i, 0 );
 
    //      double bQ = btof->getBTofQ(i);
    double bSpeed0 = btof->getBTofSpeed( i, 0 );
    double bSpeed1 = btof->getBTofSpeed( i, 1 );
    //      double bW1 = btof->getBTofW(i,0);
    //      double bW2 = btof->getBTofW(i,1);
    //      double bW3 = btof->getBTofW(i,2);
    //      double bW4 = btof->getBTofW(i,3);
    std::cout << "=========================================="
              << "\n";
    std::cout << "cnt =" << i << "\n";
    std::cout << "Atten0=" << bAtten0 << "Atten1=" << bAtten1 << "Q="
              << "bQ "
              << "Speed0=" << bSpeed0 << "Speed1=" << bSpeed1 << "P0=" << bP0 << "P1=" << bP1
              << "P2=" << bP2 << "P3=" << bP3 << "P4=" << bP4 << "P10=" << bP10;
    std::cout << "\n";
    std::cout << "P5=" << bP5 << "Q="
              << "bQ";
    //      std::cout<<"W1="<<bW1<<"W2="<<bW2<<"W3="<<bW3<<"W4="<<bW4;
    std::cout << "\n";
 
    /*
    std::cout<<"Atten[1]="<<*(bAtten+1)<<"P0[1]="<<*(bP0+1)<<"P1[1]="<<*(bP1+1)<<"P2[1]="<<*(bP2+1)<<"P3[1]="<<*(bP3+1)<<"P4[1]="<<*(bP4+1);
    std::cout<<"\n";
    std::cout<<"P5[1]="<<*(bP5+1)<<"S[1]="<<*(bS+1)<<"SS[1]="<<*(bSS+1)<<"Veff[1]="<<*(bVeff+1);
    */
  }
  for ( int i = 0; i < 96; i++ )
  {
    double eAtten0 = btof->getETofAtten( i, 0 );
    double eAtten1 = btof->getETofAtten( i, 1 );
    double eP0     = btof->getETofP( i, 0 );
    double eP1     = btof->getETofP( i, 1 );
    double eP2     = btof->getETofP( i, 2 );
    double eP3     = btof->getETofP( i, 3 );
    double eP4     = btof->getETofP( i, 4 );
    double eP5     = btof->getETofP( i, 5 );
    double eP6     = btof->getETofP( i, 6 );
 
    double eSpeed0 = btof->getETofSpeed( i, 0 );
    double eSpeed1 = btof->getETofSpeed( i, 1 );
 
    std::cout << "=========================================="
              << "\n";
    std::cout << "cnt =" << i << "\n";
 
    std::cout << "Atten0=" << eAtten0 << "Atten1=" << eAtten1 << "Speed0=" << eSpeed0
              << "Speed1=" << eSpeed1 << "P0=" << eP0 << "P1=" << eP1 << "P2=" << eP2
              << "P3=" << eP3 << "P4=" << eP4 << "P6=" << eP6;
    std::cout << "\n";
    std::cout << "P5=" << eP5;
    std::cout << "\n";
  }
  // maqm comment because cann't write
  // m_pTreeSvc->writeToTree("./Tof_test.root", fullPath);
 
  return StatusCode::SUCCESS;
}
/*
DataSvc::retrieveObject(fullPath,CalibData::CalibTest1);
*/
StatusCode checkEstTof::finalize() {
 
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "        checkEstTof FINALIZE!! " << endmsg;
 
  return StatusCode::SUCCESS;
}