TreeTofQElecCnv.cxx

Go to the documentation of this file.

#include "GaudiKernel/MsgStream.h"
#include "TreeTofQElecCnv.h"
#include "CalibData/Tof/TofElecData.h"
#include "CalibDataSvc/IInstrumentName.h"
#include "CalibMySQLCnv/TreeAddress.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<TreeTofElecDataCnv> TofCal_factory;
//const  ICnvFactory& TreeTofElecDataCnvFactory = TofCal_factory;
 
 
 
TreeTofElecDataCnv::TreeTofElecDataCnv( ISvcLocator* svc) :
    TreeCalBaseCnv(svc, CLID_TofQ_Elec) { 
 
    }
 
 
const CLID& TreeTofElecDataCnv::objType() const {
    return CLID_TofQ_Elec;
}
 
const CLID& TreeTofElecDataCnv::classID() {
    return CLID_TofQ_Elec;
}
 
StatusCode TreeTofElecDataCnv::i_createObj(IOpaqueAddress* addr,
        DataObject*& refpObject) {
 
    MsgStream log(msgSvc(), "TreeTofElecDataCnv");
    log<<MSG::DEBUG<<"SetProperty"<<endreq;
 
    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);
    TBufferFile *buf4 = new TBufferFile(TBuffer::kRead);
    TBufferFile *buf5 = new TBufferFile(TBuffer::kRead);
    TBufferFile *buf6 = new TBufferFile(TBuffer::kRead);
    TBufferFile *buf7 = new TBufferFile(TBuffer::kRead);
    TBufferFile *buf8 = new TBufferFile(TBuffer::kRead);
 
    buf1->SetBuffer((*records)["BarBoardNum"],512000,kFALSE);
    buf2->SetBuffer((*records)["EndBoardNum"],512000,kFALSE);
    buf3->SetBuffer((*records)["QELecBarParEast"],512000,kFALSE);
    buf4->SetBuffer((*records)["QELecBarParWest"],512000,kFALSE);
    buf5->SetBuffer((*records)["QELecEndPar"],512000,kFALSE);
    buf6->SetBuffer((*records)["SimQELecBarParEast"],512000,kFALSE);
    buf7->SetBuffer((*records)["SimQELecBarParWest"],512000,kFALSE);
    buf8->SetBuffer((*records)["SimQELecEndPar"],512000,kFALSE);
 
    std::cout<<" SftVer is "<<(*records)["SftVer"];
    std::cout<<" CalVerSft is "<<(*records)["CalParVer"];
    std::cout<<"  File name  is "<<(*records)["FileName"]<<std::endl;
 
 
    TTree *BarBoardNum= new TTree();
    TTree *EndBoardNum= new TTree();
    TTree *QElecBarParEast = new TTree();
    TTree *QElecBarParWest = new TTree();
    TTree *QElecEndPar = new TTree();
    TTree *SimQElecBarParEast = new TTree();
    TTree *SimQElecBarParWest = new TTree();
    TTree *SimQElecEndPar = new TTree();
    BarBoardNum->Streamer(*buf1);
    EndBoardNum->Streamer(*buf2);
    QElecBarParEast->Streamer(*buf3);
    QElecBarParWest->Streamer(*buf4);
    QElecEndPar->Streamer(*buf5);
    SimQElecBarParEast->Streamer(*buf6);
    SimQElecBarParWest->Streamer(*buf7);
    SimQElecEndPar->Streamer(*buf8);
 
 
    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;
 
 
}