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

Go to the documentation of this file.

#include "GaudiKernel/MsgStream.h"
#include "TreeMucCalibDataCnv.h"
#include "CalibData/Muc/MucCalibData.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 "TH1F.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<TreeMucCalibDataCnv> MucCal_factory;
//const  ICnvFactory& TreeMucCalibDataCnvFactory = MucCal_factory;

using CalibData::MucCalibData;
using CalibData::MucIdTransform;

TreeMucCalibDataCnv::TreeMucCalibDataCnv( ISvcLocator* svc) :
        TreeCalBaseCnv(svc, CLID_Calib_MucCal) { 
                m_ptrIdTr = new MucIdTransform();
        }


const CLID& TreeMucCalibDataCnv::objType() const {
        return CLID_Calib_MucCal;
}

const CLID& TreeMucCalibDataCnv::classID() {
        return CLID_Calib_MucCal;
}

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

        MsgStream log(msgSvc(), "TreeMucCalibDataCnv");
        log<<MSG::DEBUG<<"SetProperty"<<endreq;
        CalibData::MucCalibData *tmpObject = new CalibData::MucCalibData ;  
        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)["LayTree"],512000,kFALSE);
        buf2->SetBuffer((*records)["BoxTree"],512000,kFALSE);
        buf3->SetBuffer((*records)["StrTree"],512000,kFALSE);
        std::cout<<" SftVer is "<<(*records)["SftVer"];
        std::cout<<" CalVerSft is "<<(*records)["CalParVer"];
        std::cout<<"  File name  is "<<(*records)["FileName"]<<std::endl;


        TTree* laytree = new TTree(); 
        laytree->Streamer(*buf1); 
        log<<MSG::DEBUG<< "laytree name:\t" << laytree->GetName() << endreq;
        TTree* boxtree = new TTree(); 
        boxtree->Streamer(*buf2);
        log<<MSG::DEBUG<< "boxtree name:\t" << boxtree->GetName() << endreq;
        TTree* strtree= new TTree(); 
        strtree->Streamer(*buf3); 
        log<<MSG::DEBUG<< "strtree name:\t" << strtree->GetName() << endreq;

        // Read in the object
        //read layerConst-----------------------------
        Double_t lay_eff,box_eff,str_eff;
        Double_t lay_cnt, box_cnt, str_cnt;
        Double_t lay_nos, box_nos, str_nos;
        Double_t lay_nos_ratio, box_nos_ratio, str_nos_ratio;
        lay_eff = box_eff = str_eff = 0.0;
        lay_cnt = box_cnt = str_cnt = 0.0;
        lay_nos = box_nos = str_nos = 0.0;
        lay_nos_ratio = box_nos_ratio = str_nos_ratio = 0.0;

        char name[60];
        laytree->SetBranchAddress("layer_eff", &lay_eff);
        laytree->SetBranchAddress("layer_cnt", &lay_cnt);
        laytree->SetBranchAddress("layer_noise", &lay_nos);
        laytree->SetBranchAddress("layer_nosratio",&lay_nos_ratio);

        boxtree->SetBranchAddress("box_eff", &box_eff);
        boxtree->SetBranchAddress("box_cnt", &box_cnt);
        boxtree->SetBranchAddress("box_noise", &box_nos);
        boxtree->SetBranchAddress("box_nosratio", &box_nos_ratio); 

        strtree->SetBranchAddress("strip_eff", &str_eff);
        strtree->SetBranchAddress("strip_cnt", &str_cnt);
        strtree->SetBranchAddress("strip_noise", &str_nos);
        strtree->SetBranchAddress("strip_nosratio", &str_nos_ratio);

        int part, segment, layer, strip;
        part = segment = layer = strip = 0;
        for(int i=0; i<LAYER_MAX; i++)
        {
                laytree->GetEntry(i);
                tmpObject->setLayerEff(lay_eff,i);
                tmpObject->setLayerCnt(lay_cnt, i);
                tmpObject->setLayerNos(lay_nos, i);
                tmpObject->setLayerNosRatio(lay_nos_ratio, i);
                sprintf(name,"LayClstPro");
                for(int j=0; j<CLST_MAX; j++) {
                        //                if( ht_Pro[0] != NULL ) tmpObject->setLayerClstPro(ht_Pro[0]->GetBinContent(j),i,j);
                        //          else tmpObject->setLayerClstPro(DEFAULT_CLST_PRO[j],i,j);
                        tmpObject->setLayerClstPro(DEFAULT_CLST_PRO[j],i,j);
                }
                log<<MSG::DEBUG<<"layer: " << i << "\t" << lay_eff <<endreq;            
        }

        for(int i=0; i<BOX_MAX; i++)
        {
                m_ptrIdTr->SetBoxPos( i, &part, &segment, &layer );
                boxtree->GetEntry(i);           
                tmpObject->setBoxEff(box_eff, part, segment, layer);
                tmpObject->setBoxCnt(box_cnt, part, segment, layer);
                tmpObject->setBoxNos(box_nos, part, segment, layer);
                tmpObject->setBoxNosRatio(box_nos_ratio, part, segment, layer);

                sprintf(name,"BoxClstPro_B%d",i);
                for(int j=0; j<CLST_MAX; j++) {
                        //                if( ht_Pro[1] != NULL ) tmpObject->setBoxClstPro(ht_Pro[1]->GetBinContent(j),part,segment,layer,j);
                        //        else
                        tmpObject->setBoxClstPro(DEFAULT_CLST_PRO[j],part,segment,layer,j);
                }
                log<<MSG::DEBUG<<"box: " << i << "\t" << box_eff <<endreq;               
        }

        log<<MSG::DEBUG<<"strtree entries:\t" << strtree->GetEntries() << endreq;
        for(int i=0; i<STRIP_MAX; i++)
        {       
                m_ptrIdTr->SetStripPos( i, &part, &segment, &layer, &strip );
                strtree->GetEntry(i);
                tmpObject->setStripEff(str_eff, part, segment, layer, strip);
                tmpObject->setStripCnt(str_cnt, part, segment, layer, strip);
                tmpObject->setStripNos(str_nos, part, segment, layer, strip);
                tmpObject->setStripNosRatio(str_nos_ratio, part, segment, layer, strip);
                log<<MSG::DEBUG<<"strip: " << i << "\t" << str_eff <<endreq;               

        }  

        refpObject=tmpObject;
        delete laytree;
        delete boxtree;
        delete strtree;
        return StatusCode::SUCCESS;


}

---