TreeTofSimCnv.cxx

Go to the documentation of this file.

#include "TreeTofSimCnv.h"
#include "CalibData/Tof/TofSimData.h"
#include "CalibDataSvc/IInstrumentName.h"
#include "CalibMySQLCnvSvc/TreeAddress.h"
#include "GaudiKernel/MsgStream.h"
#include "TBuffer.h"
#include "TDirectory.h"
#include "TFile.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<TreeTofSimDataCnv> TofCal_factory;
// const  ICnvFactory& TreeTofSimDataCnvFactory = TofCal_factory;
 
TreeTofSimDataCnv::TreeTofSimDataCnv( ISvcLocator* svc )
    : TreeCalBaseCnv( svc, CLID_Calib_TofSim ) {}
 
const CLID& TreeTofSimDataCnv::objType() const { return CLID_Calib_TofSim; }
 
const CLID& TreeTofSimDataCnv::classID() { return CLID_Calib_TofSim; }
 
StatusCode TreeTofSimDataCnv::i_createObj( IOpaqueAddress* addr, DataObject*& refpObject ) {
 
  MsgStream log( msgSvc(), "TreeTofSimDataCnv" );
  log << MSG::DEBUG << "SetProperty" << endmsg;
 
  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 )["BTofSim"], 512000, kFALSE );
  buf2->SetBuffer( ( *records )["ETofSim"], 512000, kFALSE );
  buf3->SetBuffer( ( *records )["SimConstants"], 512000, kFALSE );
 
  std::cout << " SftVer is " << ( *records )["SftVer"];
  std::cout << " CalVerSft is " << ( *records )["TofSimParVer"];
  std::cout << "  File name  is " << ( *records )["FileName"] << std::endl;
 
  TTree* btoftree       = new TTree();
  TTree* etoftree       = new TTree();
  TTree* btofcommontree = new TTree();
  btoftree->Streamer( *buf1 );
  etoftree->Streamer( *buf2 );
  btofcommontree->Streamer( *buf3 );
 
  CalibData::BTofSimBase                  bTof;
  CalibData::ETofSimBase                  eTof;
  CalibData::TofSimConstBase              tofbase;
  std::vector<CalibData::BTofSimBase>     tmpbTof; //; = new vector<CalibData::bTofCalibBase>;
  std::vector<CalibData::ETofSimBase>     tmpeTof;
  std::vector<CalibData::TofSimConstBase> tofbaseCol;
  // Read in the object
  int cnt;
  // read btoftree ------------------------------------------------------------
  double AttenLength;
  double Gain;
  double Ratio;
  double NoiseSmear;
  btoftree->SetBranchAddress( "AttenLength", &AttenLength );
  btoftree->SetBranchAddress( "Gain", &Gain );
  btoftree->SetBranchAddress( "Ratio", &Ratio );
  int entries = btoftree->GetEntries();
  for ( cnt = 0; cnt < entries; cnt++ )
  {
    btoftree->GetEntry( cnt );
    bTof.setGain( Gain );
    bTof.setRatio( Ratio );
    bTof.setAttenLength( AttenLength );
    tmpbTof.push_back( bTof );
  }
  // read etoftree
  etoftree->SetBranchAddress( "Gain", &Gain );
  etoftree->SetBranchAddress( "AttenLength", &AttenLength );
  etoftree->SetBranchAddress( "NoiseSmear", &NoiseSmear );
  entries = etoftree->GetEntries();
  for ( cnt = 0; cnt < entries; cnt++ )
  {
    etoftree->GetEntry( cnt );
    eTof.setGain( Gain );
    eTof.setAttenLength( AttenLength );
    eTof.setNoiseSmear( NoiseSmear );
    tmpeTof.push_back( eTof );
  }
  // read SimConstants
  double BarConstant, BarPMTGain, BarHighThres, BarLowThres, EndConstant, EndPMTGain,
      EndHighThres, EndLowThres, EndNoiseSwitch;
  btofcommontree->SetBranchAddress( "BarConstant", &BarConstant );
  btofcommontree->SetBranchAddress( "BarPMTGain", &BarPMTGain );
  btofcommontree->SetBranchAddress( "BarHighThres", &BarHighThres );
  btofcommontree->SetBranchAddress( "BarLowThres", &BarLowThres );
  btofcommontree->SetBranchAddress( "EndConstant", &EndConstant );
  btofcommontree->SetBranchAddress( "EndPMTGain", &EndPMTGain );
  btofcommontree->SetBranchAddress( "EndHighThres", &EndHighThres );
  btofcommontree->SetBranchAddress( "EndLowThres", &EndLowThres );
  btofcommontree->SetBranchAddress( "EndNoiseSwitch", &EndNoiseSwitch );
  entries = btofcommontree->GetEntries();
  for ( cnt = 0; cnt < entries; cnt++ )
  {
    btofcommontree->GetEntry( cnt );
    tofbase.setBarLowThres( BarLowThres );
    tofbase.setBarHighThres( BarHighThres );
    tofbase.setEndLowThres( EndLowThres );
    tofbase.setEndHighThres( EndHighThres );
    tofbase.setBarPMTGain( BarPMTGain );
    tofbase.setEndPMTGain( EndPMTGain );
    tofbase.setBarConstant( BarConstant );
    tofbase.setEndConstant( EndConstant );
    tofbase.setEndNoiseSwitch( EndNoiseSwitch );
    tofbaseCol.push_back( tofbase );
  }
 
  CalibData::TofSimData* tmpObject =
      new CalibData::TofSimData( &tofbaseCol, &tmpbTof, &tmpeTof );
  refpObject = tmpObject;
  delete btoftree;
  delete etoftree;
  delete btofcommontree;
 
  return StatusCode::SUCCESS;
}