/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Calibration/CalibSvc/CalibTreeCnv/CalibTreeCnv-00-01-18/src/cnv/TreeDedxSimDataCnv.cxx

Go to the documentation of this file.

#include "GaudiKernel/MsgStream.h"
#include "TreeDedxSimDataCnv.h"
#include "CalibData/Dedx/DedxSimData.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  "TBufferFile.h"
#include  "TROOT.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<TreeDedxSimDataCnv> DedxCal_factory;
//const  ICnvFactory& TreeDedxSimDataCnvFactory = DedxCal_factory;



TreeDedxSimDataCnv::TreeDedxSimDataCnv( ISvcLocator* svc) :
        TreeCalBaseCnv(svc, CLID_Dedx_Sim) { 

        }


const CLID& TreeDedxSimDataCnv::objType() const {
        return CLID_Dedx_Sim;
}

const CLID& TreeDedxSimDataCnv::classID() {
        return CLID_Dedx_Sim;
}

StatusCode TreeDedxSimDataCnv::i_createObj(IOpaqueAddress* addr,
                DataObject*& refpObject) {

        MsgStream log(msgSvc(), "TreeDedxSimDataCnv");

        CalibData::DedxSimData *tmpObject = new CalibData::DedxSimData;
        TreeAddress* add = dynamic_cast<TreeAddress*>(addr);
        DatabaseRecord *records=add->pp();

        TBufferFile *buf1 = new TBufferFile(TBuffer::kRead);
        TBufferFile *buf2 = new TBufferFile(TBuffer::kRead);

        buf1->SetBuffer((*records)["TH1F_Col"],51200000,kFALSE);
        buf2->SetBuffer((*records)["bin"],512000,kFALSE);

        std::cout<<" SftVer is "<<(*records)["SftVer"];
        std::cout<<" CalVerSft is "<<(*records)["CalParVer"];
        std::cout<<"  File name  is "<<(*records)["FileName"]<<std::endl;


        //TClass *myclass=gROOT->GetClass("TTree");
        //TTree *tree =(TTree*)buf1->ReadObject(myclass);
        TTree* tree = new TTree();
        tree->Streamer(*buf1);

        TTree* bin = new TTree(); 
        bin->Streamer(*buf2);

        TH1F *h1=new TH1F();

        std::vector<TH1F> hist;
        std::vector<double> hRange;
        // Read in the object
        int cnt;
        tree -> SetBranchAddress("TH1F_Col", &h1);
        int entries=tree->GetEntries();
        for(cnt=0; cnt<entries; cnt++){
                tree -> GetEntry(cnt);
                gROOT->cd();
                TH1F *h2=new TH1F();
                h2=(TH1F*)h1->Clone();
                hist.push_back(*h2);
        }
    int ver;
        double betagamma[5000];
        int totalNum;
        int bgNum;
    int thetaNum;
    if (bin->FindBranch("version")) bin->SetBranchAddress("version", &ver);
    else ver = 0;
    if (bin->FindBranch("thetaNum")) bin->SetBranchAddress("thetaNum", &thetaNum);
    else thetaNum = 16;
        bin-> SetBranchAddress("totalNum", &totalNum);
        bin-> SetBranchAddress("betagamma",  betagamma);
        bin-> SetBranchAddress("betagammaBounds", &bgNum);
        bin->GetEntry(0);
        //int tot=totalNum*2;
        for(cnt=0; cnt<bgNum; cnt++){
                hRange.push_back(betagamma[cnt]);
        }
    tmpObject->setVersion(ver);
        tmpObject->setRangeNo(bgNum);
        tmpObject->setRange(&hRange);
        tmpObject->setHistNo(totalNum);
    tmpObject->setThetaNo(thetaNum);
        tmpObject->setHist(&hist);
        refpObject=tmpObject;

        delete tree;
        delete bin;

        return StatusCode::SUCCESS;


}

---