/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Calibration/CalibSvc/CalibROOTCnv/CalibROOTCnv-00-01-13/src/cnv/RootTofQElecDataCnv.cxx

Go to the documentation of this file.

#include "GaudiKernel/MsgStream.h"
#include "RootTofQElecDataCnv.h"
#include "CalibData/Tof/TofElecData.h"
#include "CalibDataSvc/IInstrumentName.h"
#include "TFile.h"
#include "TTree.h"
#include "TDirectory.h"
#include "TObject.h"
#include "TBuffer.h"

#include "GaudiKernel/CnvFactory.h"
#include "GaudiKernel/IOpaqueAddress.h"
#include "GaudiKernel/DataObject.h"
#include "GaudiKernel/IAddressCreator.h"
#include "GaudiKernel/IDataProviderSvc.h"
#include "GaudiKernel/IConversionSvc.h"
#include "GaudiKernel/GenericAddress.h"

#include "CalibDataSvc/ICalibTreeSvc.h"    //maybe
#include "CalibDataSvc/ICalibMetaCnvSvc.h"

// 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"<<endreq;
  StatusCode sc = openRead(fname);
  if(!sc)
    { log<<MSG::ERROR<<"unable to open files"<<endreq;
    }

  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;


}

---