RootTofQElecDataCnv.cxx

Go to the documentation of this file.

#include "RootTofQElecDataCnv.h"
#include "CalibData/Tof/TofElecData.h"
#include "CalibDataSvc/IInstrumentName.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<RootTofElecDataCnv> TofCal_factory;
// const  ICnvFactory& RootTofElecDataCnvFactory = TofCal_factory;
 
RootTofElecDataCnv::RootTofElecDataCnv( ISvcLocator* svc )
    : RootCalBaseCnv( svc, CLID_TofQ_Elec ) {}
 
const CLID& RootTofElecDataCnv::objType() const { return CLID_TofQ_Elec; }
 
const CLID& RootTofElecDataCnv::classID() { return CLID_TofQ_Elec; }
 
StatusCode RootTofElecDataCnv::i_createObj( const std::string& fname,
                                            DataObject*&       refpObject ) {
 
  MsgStream log( msgSvc(), "RootTofElecDataCnv" );
  log << MSG::DEBUG << "SetProperty" << endmsg;
  StatusCode sc = openRead( fname );
  if ( !sc ) { log << MSG::ERROR << "unable to open files" << endmsg; }
 
  TTree* BarBoardNum        = (TTree*)m_inFile->Get( "BarBoardNum" );
  TTree* EndBoardNum        = (TTree*)m_inFile->Get( "EndBoardNum" );
  TTree* QElecBarParEast    = (TTree*)m_inFile->Get( "QElecBarParEast" );
  TTree* QElecBarParWest    = (TTree*)m_inFile->Get( "QElecBarParWest" );
  TTree* QElecEndPar        = (TTree*)m_inFile->Get( "QElecEndPar" );
  TTree* SimQElecBarParEast = (TTree*)m_inFile->Get( "SimQElecBarParEast" );
  TTree* SimQElecBarParWest = (TTree*)m_inFile->Get( "SimQElecBarParWest" );
  TTree* SimQElecEndPar     = (TTree*)m_inFile->Get( "SimQElecEndPar" );
 
  CalibData::BarrelElec bTof;
  CalibData::EndcapElec eTof;
 
  std::vector<CalibData::BarrelElec> tmpbTof; //; = new vector<CalibData::bTofCalibBase>;
  std::vector<CalibData::EndcapElec> tmpeTof;
 
  // Read in the object
  int cnt;
  // read btoftree ------------------------------------------------------------
  unsigned int num[4];
  double       num1[4];
  BarBoardNum->SetBranchAddress( "Board", &num1[0] );
  BarBoardNum->SetBranchAddress( "Channel", &num1[1] );
  BarBoardNum->SetBranchAddress( "Crate", &num1[2] );
  BarBoardNum->SetBranchAddress( "Fee", &num1[3] );
  double p[11];
  QElecBarParEast->SetBranchAddress( "P0", &p[0] );
  QElecBarParEast->SetBranchAddress( "P1", &p[1] );
  QElecBarParEast->SetBranchAddress( "P2", &p[2] );
  QElecBarParEast->SetBranchAddress( "P3", &p[3] );
  QElecBarParEast->SetBranchAddress( "P4", &p[4] );
  QElecBarParEast->SetBranchAddress( "P5", &p[5] );
  QElecBarParEast->SetBranchAddress( "P6", &p[6] );
  QElecBarParEast->SetBranchAddress( "P7", &p[7] );
  QElecBarParEast->SetBranchAddress( "P8", &p[8] );
  QElecBarParEast->SetBranchAddress( "P9", &p[9] );
  QElecBarParEast->SetBranchAddress( "P10", &p[10] );
  double p_w[11];
  QElecBarParWest->SetBranchAddress( "P0", &p_w[0] );
  QElecBarParWest->SetBranchAddress( "P1", &p_w[1] );
  QElecBarParWest->SetBranchAddress( "P2", &p_w[2] );
  QElecBarParWest->SetBranchAddress( "P3", &p_w[3] );
  QElecBarParWest->SetBranchAddress( "P4", &p_w[4] );
  QElecBarParWest->SetBranchAddress( "P5", &p_w[5] );
  QElecBarParWest->SetBranchAddress( "P6", &p_w[6] );
  QElecBarParWest->SetBranchAddress( "P7", &p_w[7] );
  QElecBarParWest->SetBranchAddress( "P8", &p_w[8] );
  QElecBarParWest->SetBranchAddress( "P9", &p_w[9] );
  QElecBarParWest->SetBranchAddress( "P10", &p_w[10] );
  double p_se[11];
  SimQElecBarParEast->SetBranchAddress( "P0", &p_se[0] );
  SimQElecBarParEast->SetBranchAddress( "P1", &p_se[1] );
  SimQElecBarParEast->SetBranchAddress( "P2", &p_se[2] );
  SimQElecBarParEast->SetBranchAddress( "P3", &p_se[3] );
  SimQElecBarParEast->SetBranchAddress( "P4", &p_se[4] );
  SimQElecBarParEast->SetBranchAddress( "P5", &p_se[5] );
  SimQElecBarParEast->SetBranchAddress( "P6", &p_se[6] );
  SimQElecBarParEast->SetBranchAddress( "P7", &p_se[7] );
  SimQElecBarParEast->SetBranchAddress( "P8", &p_se[8] );
  SimQElecBarParEast->SetBranchAddress( "P9", &p_se[9] );
  SimQElecBarParEast->SetBranchAddress( "P10", &p_se[10] );
  double p_sw[11];
  SimQElecBarParWest->SetBranchAddress( "P0", &p_sw[0] );
  SimQElecBarParWest->SetBranchAddress( "P1", &p_sw[1] );
  SimQElecBarParWest->SetBranchAddress( "P2", &p_sw[2] );
  SimQElecBarParWest->SetBranchAddress( "P3", &p_sw[3] );
  SimQElecBarParWest->SetBranchAddress( "P4", &p_sw[4] );
  SimQElecBarParWest->SetBranchAddress( "P5", &p_sw[5] );
  SimQElecBarParWest->SetBranchAddress( "P6", &p_sw[6] );
  SimQElecBarParWest->SetBranchAddress( "P7", &p_sw[7] );
  SimQElecBarParWest->SetBranchAddress( "P8", &p_sw[8] );
  SimQElecBarParWest->SetBranchAddress( "P9", &p_sw[9] );
  SimQElecBarParWest->SetBranchAddress( "P10", &p_sw[10] );
 
  for ( cnt = 0; cnt < 352; cnt++ )
  {
    BarBoardNum->GetEntry( cnt );
    num[0] = num1[0];
    num[1] = num1[1];
    num[2] = num1[2];
    num[3] = num1[3];
    if ( cnt % 2 != 0 ) { bTof.setNumWest( num ); }
    if ( cnt % 2 == 0 )
    {
      bTof.setNumEast( num );
      int kkk = cnt / 2;
      QElecBarParEast->GetEntry( kkk );
      QElecBarParWest->GetEntry( kkk );
      SimQElecBarParEast->GetEntry( kkk );
      SimQElecBarParWest->GetEntry( kkk );
      bTof.setP1( p );
      bTof.setP2( p_w );
      bTof.setSimP1( p_se );
      bTof.setSimP2( p_sw );
    }
    if ( cnt % 2 != 0 ) tmpbTof.push_back( bTof );
  }
 
  // read etoftree
  EndBoardNum->SetBranchAddress( "Board", &num1[0] );
  EndBoardNum->SetBranchAddress( "Channel", &num1[1] );
  EndBoardNum->SetBranchAddress( "Crate", &num1[2] );
  EndBoardNum->SetBranchAddress( "Fee", &num1[3] );
  QElecEndPar->SetBranchAddress( "P0", &p[0] );
  QElecEndPar->SetBranchAddress( "P1", &p[1] );
  QElecEndPar->SetBranchAddress( "P2", &p[2] );
  QElecEndPar->SetBranchAddress( "P3", &p[3] );
  QElecEndPar->SetBranchAddress( "P4", &p[4] );
  QElecEndPar->SetBranchAddress( "P5", &p[5] );
  QElecEndPar->SetBranchAddress( "P6", &p[6] );
  QElecEndPar->SetBranchAddress( "P7", &p[7] );
  QElecEndPar->SetBranchAddress( "P8", &p[8] );
  QElecEndPar->SetBranchAddress( "P9", &p[9] );
  QElecEndPar->SetBranchAddress( "P10", &p[10] );
  SimQElecEndPar->SetBranchAddress( "P0", &p_w[0] );
  SimQElecEndPar->SetBranchAddress( "P1", &p_w[1] );
  SimQElecEndPar->SetBranchAddress( "P2", &p_w[2] );
  SimQElecEndPar->SetBranchAddress( "P3", &p_w[3] );
  SimQElecEndPar->SetBranchAddress( "P4", &p_w[4] );
  SimQElecEndPar->SetBranchAddress( "P5", &p_w[5] );
  SimQElecEndPar->SetBranchAddress( "P6", &p_w[6] );
  SimQElecEndPar->SetBranchAddress( "P7", &p_w[7] );
  SimQElecEndPar->SetBranchAddress( "P8", &p_w[8] );
  SimQElecEndPar->SetBranchAddress( "P9", &p_w[9] );
  SimQElecEndPar->SetBranchAddress( "P10", &p_w[10] );
  for ( cnt = 0; cnt < 96; cnt++ )
  {
    EndBoardNum->GetEntry( cnt );
    QElecEndPar->GetEntry( cnt );
    SimQElecEndPar->GetEntry( cnt );
    num[0] = num1[0];
    num[1] = num1[1];
    num[2] = num1[2];
    num[3] = num1[3];
    eTof.setNum( num );
    eTof.setP( p );
    eTof.setSimP( p_w );
    tmpeTof.push_back( eTof );
  }
 
  CalibData::TofElecData* tmpObject = new CalibData::TofElecData( &tmpbTof, &tmpeTof );
 
  refpObject = tmpObject;
  delete BarBoardNum;
  delete EndBoardNum;
  delete QElecBarParEast;
  delete QElecBarParWest;
  delete QElecEndPar;
  delete SimQElecBarParEast;
  delete SimQElecBarParWest;
  delete SimQElecEndPar;
 
  return StatusCode::SUCCESS;
}