TreeMucCalibDataCnv.cxx

Go to the documentation of this file.

#include "TreeMucCalibDataCnv.h"
#include "CalibData/Muc/MucCalibData.h"
#include "CalibDataSvc/IInstrumentName.h"
#include "CalibMySQLCnvSvc/TreeAddress.h"
#include "GaudiKernel/MsgStream.h"
#include "TBuffer.h"
#include "TDirectory.h"
#include "TFile.h"
#include "TH1F.h"
#include "TObject.h"
#include "TTree.h"
 
#include "GaudiKernel/DataObject.h"
#include "GaudiKernel/GenericAddress.h"
#include "GaudiKernel/IAddressCreator.h"
#include "GaudiKernel/IConversionSvc.h"
#include "GaudiKernel/IDataProviderSvc.h"
#include "GaudiKernel/IOpaqueAddress.h"
 
#include "CalibDataSvc/ICalibMetaCnvSvc.h"
#include "CalibDataSvc/ICalibTreeSvc.h" //maybe
 
// Temporary.  Hope to find a better way to do this
#include "CalibData/CalibModel.h"
// using namespace CalibData;
// static CnvFactory<TreeMucCalibDataCnv> MucCal_factory;
// const  ICnvFactory& TreeMucCalibDataCnvFactory = MucCal_factory;
 
using CalibData::MucCalibData;
using CalibData::MucIdTransform;
 
TreeMucCalibDataCnv::TreeMucCalibDataCnv( ISvcLocator* svc )
    : TreeCalBaseCnv( svc, CLID_Calib_MucCal ) {
  m_ptrIdTr = new MucIdTransform();
}
 
const CLID& TreeMucCalibDataCnv::objType() const { return CLID_Calib_MucCal; }
 
const CLID& TreeMucCalibDataCnv::classID() { return CLID_Calib_MucCal; }
 
StatusCode TreeMucCalibDataCnv::i_createObj( IOpaqueAddress* addr, DataObject*& refpObject ) {
 
  MsgStream log( msgSvc(), "TreeMucCalibDataCnv" );
  log << MSG::DEBUG << "SetProperty" << endmsg;
  CalibData::MucCalibData* tmpObject = new CalibData::MucCalibData;
  TreeAddress*             add       = dynamic_cast<TreeAddress*>( addr );
  DatabaseRecord*          records   = add->pp();
 
  TBufferFile* buf1 = new TBufferFile( TBuffer::kRead );
  TBufferFile* buf2 = new TBufferFile( TBuffer::kRead );
  TBufferFile* buf3 = new TBufferFile( TBuffer::kRead );
  buf1->SetBuffer( ( *records )["LayTree"], 512000, kFALSE );
  buf2->SetBuffer( ( *records )["BoxTree"], 512000, kFALSE );
  buf3->SetBuffer( ( *records )["StrTree"], 512000, kFALSE );
  std::cout << " SftVer is " << ( *records )["SftVer"];
  std::cout << " CalVerSft is " << ( *records )["CalParVer"];
  std::cout << "  File name  is " << ( *records )["FileName"] << std::endl;
 
  TTree* laytree = new TTree();
  laytree->Streamer( *buf1 );
  log << MSG::DEBUG << "laytree name:\t" << laytree->GetName() << endmsg;
  TTree* boxtree = new TTree();
  boxtree->Streamer( *buf2 );
  log << MSG::DEBUG << "boxtree name:\t" << boxtree->GetName() << endmsg;
  TTree* strtree = new TTree();
  strtree->Streamer( *buf3 );
  log << MSG::DEBUG << "strtree name:\t" << strtree->GetName() << endmsg;
 
  // Read in the object
  // read layerConst-----------------------------
  Double_t lay_eff, box_eff, str_eff;
  Double_t lay_cnt, box_cnt, str_cnt;
  Double_t lay_nos, box_nos, str_nos;
  Double_t lay_nos_ratio, box_nos_ratio, str_nos_ratio;
  lay_eff = box_eff = str_eff = 0.0;
  lay_cnt = box_cnt = str_cnt = 0.0;
  lay_nos = box_nos = str_nos = 0.0;
  lay_nos_ratio = box_nos_ratio = str_nos_ratio = 0.0;
 
  char name[60];
  laytree->SetBranchAddress( "layer_eff", &lay_eff );
  laytree->SetBranchAddress( "layer_cnt", &lay_cnt );
  laytree->SetBranchAddress( "layer_noise", &lay_nos );
  laytree->SetBranchAddress( "layer_nosratio", &lay_nos_ratio );
 
  boxtree->SetBranchAddress( "box_eff", &box_eff );
  boxtree->SetBranchAddress( "box_cnt", &box_cnt );
  boxtree->SetBranchAddress( "box_noise", &box_nos );
  boxtree->SetBranchAddress( "box_nosratio", &box_nos_ratio );
 
  strtree->SetBranchAddress( "strip_eff", &str_eff );
  strtree->SetBranchAddress( "strip_cnt", &str_cnt );
  strtree->SetBranchAddress( "strip_noise", &str_nos );
  strtree->SetBranchAddress( "strip_nosratio", &str_nos_ratio );
 
  int part, segment, layer, strip;
  part = segment = layer = strip = 0;
  for ( int i = 0; i < LAYER_MAX; i++ )
  {
    laytree->GetEntry( i );
    tmpObject->setLayerEff( lay_eff, i );
    tmpObject->setLayerCnt( lay_cnt, i );
    tmpObject->setLayerNos( lay_nos, i );
    tmpObject->setLayerNosRatio( lay_nos_ratio, i );
    sprintf( name, "LayClstPro" );
    for ( int j = 0; j < CLST_MAX; j++ )
    {
      //          if( ht_Pro[0] != NULL )
      // tmpObject->setLayerClstPro(ht_Pro[0]->GetBinContent(j),i,j);
      //          else tmpObject->setLayerClstPro(DEFAULT_CLST_PRO[j],i,j);
      tmpObject->setLayerClstPro( DEFAULT_CLST_PRO[j], i, j );
    }
    log << MSG::DEBUG << "layer: " << i << "\t" << lay_eff << endmsg;
  }
 
  for ( int i = 0; i < BOX_MAX; i++ )
  {
    m_ptrIdTr->SetBoxPos( i, &part, &segment, &layer );
    boxtree->GetEntry( i );
    tmpObject->setBoxEff( box_eff, part, segment, layer );
    tmpObject->setBoxCnt( box_cnt, part, segment, layer );
    tmpObject->setBoxNos( box_nos, part, segment, layer );
    tmpObject->setBoxNosRatio( box_nos_ratio, part, segment, layer );
 
    sprintf( name, "BoxClstPro_B%d", i );
    for ( int j = 0; j < CLST_MAX; j++ )
    {
      //          if( ht_Pro[1] != NULL )
      // tmpObject->setBoxClstPro(ht_Pro[1]->GetBinContent(j),part,segment,layer,j);      else
      tmpObject->setBoxClstPro( DEFAULT_CLST_PRO[j], part, segment, layer, j );
    }
    log << MSG::DEBUG << "box: " << i << "\t" << box_eff << endmsg;
  }
 
  log << MSG::DEBUG << "strtree entries:\t" << strtree->GetEntries() << endmsg;
  for ( int i = 0; i < STRIP_MAX; i++ )
  {
    m_ptrIdTr->SetStripPos( i, &part, &segment, &layer, &strip );
    strtree->GetEntry( i );
    tmpObject->setStripEff( str_eff, part, segment, layer, strip );
    tmpObject->setStripCnt( str_cnt, part, segment, layer, strip );
    tmpObject->setStripNos( str_nos, part, segment, layer, strip );
    tmpObject->setStripNosRatio( str_nos_ratio, part, segment, layer, strip );
    log << MSG::DEBUG << "strip: " << i << "\t" << str_eff << endmsg;
  }
 
  refpObject = tmpObject;
  delete laytree;
  delete boxtree;
  delete strtree;
  return StatusCode::SUCCESS;
}