MdcCalibData.cxx

Go to the documentation of this file.

#include "CalibData/Mdc/MdcCalibData.h"
//typedef std::map<int, double>::value_type valType;
#include "GaudiKernel/MsgStream.h"
namespace CalibData {
  MdcCalibData::MdcCalibData(){
     newXtTreesBuf = NULL;
     r2tTreesBuf = NULL;
     newXtTreesArray = NULL;
     r2tTreesArray = NULL;
     for(int layid=0; layid<43; layid++){
          m_r2t[layid] = NULL;
          for(int entr=0; entr<18; entr++){
          for(int lr=0; lr<2; lr++){
               m_newxt[layid][entr][lr] = NULL;
          }//end lr
      }//end entr
     }//end for layid
  } 
 
  MdcCalibData::~MdcCalibData(){
     for(int layid=0; layid<43; layid++){
      delete m_r2t[layid];
      for(int entr=0; entr<18; entr++){
           for(int lr=0; lr<2; lr++){
         if(m_newxt[layid][entr][lr]!=NULL) delete m_newxt[layid][entr][lr];
           }
      }
     }
     if(newXtTreesBuf!=NULL) newXtTreesBuf->Clear();
     if(r2tTreesBuf!=NULL) r2tTreesBuf->Clear();
     if(newXtTreesArray!=NULL) delete newXtTreesArray;
     if(r2tTreesArray!=NULL) delete r2tTreesArray;
  }
 
  StatusCode MdcCalibData::update(CalibBase1& other, MsgStream* log)
  {
    MdcCalibData& other1 = dynamic_cast<MdcCalibData& >(other);
    std::cout<<"\n"<<"here is the update in the MdcCalConst"<<std::endl;
    CalibBase1::update(other, log);
    clear();
    m_xtmap = other1.m_xtmap;
    m_t0 = other1.m_t0;
    m_delt0 = other1.m_delt0;
    m_qtpar0 = other1.m_qtpar0;
    m_qtpar1 = other1.m_qtpar1;
    m_sdmap = other1.m_sdmap;
    return StatusCode::SUCCESS;  
}   
                
void MdcCalibData::setXtpar(int xtkey,double val){
     m_xtmap.insert( valType( int(xtkey), val) );
}
 
void MdcCalibData::setNewXtpar(TObjArray* newXtTrees){
     //put TObjArray to temp buffer
     TBufferFile* tempBuf = new TBufferFile(TBuffer::kWrite);
     newXtTrees->Streamer(*tempBuf);
 
     //put content of buffer to array
     newXtTreesBuf = new TBufferFile(TBuffer::kRead);
     newXtTreesBuf->SetBuffer(tempBuf->Buffer(), tempBuf->BufferSize());
 
     newXtTreesArray = new TObjArray();
     newXtTreesArray->Streamer(*newXtTreesBuf);
     for(int layid=0; layid<43; layid++){
      for(int entr=0; entr<18; entr++){
           for(int lr=0; lr<2; lr++){
            char newXtTreeName[20];
            sprintf(newXtTreeName,"trNewXt%02d_%02d_%d",layid,entr,lr);
            m_newxt[layid][entr][lr] = (TTree*)newXtTreesArray->FindObject(newXtTreeName);; 
           }
      }
     }
     delete tempBuf;
}
 
void MdcCalibData::setR2tpar(TObjArray* r2tTrees){
     //put TObjArray to temp buffer
     TBufferFile* tempBuf = new TBufferFile(TBuffer::kWrite);
     r2tTrees->Streamer(*tempBuf);
 
     //put content of buffer to array
     r2tTreesBuf = new TBufferFile(TBuffer::kRead);
     r2tTreesBuf->SetBuffer(tempBuf->Buffer(), tempBuf->BufferSize());
 
     r2tTreesArray = new TObjArray();
     r2tTreesArray->Streamer(*r2tTreesBuf);
     for(int layid=0; layid<43; layid++){
      char r2tTreeName[20];
      sprintf(r2tTreeName,"r2t%02d",layid);
      m_r2t[layid] = (TTree*)r2tTreesArray->FindObject(r2tTreeName);; 
     }
     delete tempBuf;
}
 
void MdcCalibData::setT0(double val){
     m_t0.push_back( val );
}
void MdcCalibData::setDelT0( double val){
     m_delt0.push_back( val );
}
 
void MdcCalibData::setQtpar0(double val){
     m_qtpar0.push_back( val );
}
 
void MdcCalibData::setQtpar1(double val){
     m_qtpar1.push_back( val );
}
 
void MdcCalibData::setSdpar(int sdkey, double val){
     m_sdmap.insert( valType( int(sdkey), val) );
}
 
// wulh updated on 08-3-24
double MdcCalibData::getXtpar(int layid, int entr, int lr, int order){
     double xtpar = -999.0;
     int key = getXtKey(layid, entr, lr, order);
     if( 1 == m_xtmap.count(key) )
      xtpar = m_xtmap[key];
 
     return xtpar;
}
 
// yzhang add 2011-12-12 
TTree* MdcCalibData::getNewXtpar(int layid, int entr, int lr){
     return m_newxt[layid][entr][lr];
}
 
// yzhang add 2011-12-15
TTree* MdcCalibData::getR2tpar(int layid){
     return m_r2t[layid];
}
 
// wulh updated on 08-3-24
double MdcCalibData::getSdpar(int layid, int entr, int lr, int bin){
     double sdpar = -999.0;
     int key = getSdKey(layid, entr, lr, bin);
     if( 1 == m_sdmap.count(key) ){
      return m_sdmap[key];
     }
 
     return sdpar;
}
 
// wulh updated on 08-3-24
int MdcCalibData::getXtKey(int layid, int entr, int lr, int order){
     int key;
 
     key = ( (layid << XTLAYER_INDEX) & XTLAYER_MASK ) |
      ( (entr << XTENTRA_INDEX) & XTENTRA_MASK ) |
      ( (lr << XTLR_INDEX) & XTLR_MASK ) |
      ( (order << XTORDER_INDEX) & XTORDER_MASK );
 
     return key;
}
 
// wulh updated on 08-3-24
int MdcCalibData::getSdKey(int layid, int entr, int lr, int bin){
     int key;
 
     key = ( (layid << SDLAYER_INDEX) & SDLAYER_MASK ) |
      ( (entr << SDENTRA_INDEX) & SDENTRA_MASK ) |
      ( (lr << SDLR_INDEX) & SDLR_MASK ) |
      ( (bin << SDBIN_INDEX) & SDBIN_MASK );
 
     return key;
}
 
void MdcCalibData::setXtBegin(){
     m_xtiter = m_xtmap.begin();
}
 
int MdcCalibData::getNextXtpar(int& key, double& xtpar){
     if( m_xtiter != m_xtmap.end() ){
      key = (*m_xtiter).first;
      xtpar = (*m_xtiter).second;
      m_xtiter++;
      return 1;
     }
     else return 0;
}
 
void MdcCalibData::setSdBegin(){
     m_sditer = m_sdmap.begin();
}
 
int MdcCalibData::getNextSdpar(int& key, double& sdpar){
     if( m_sditer != m_sdmap.end() ){
      key = (*m_sditer).first;
      sdpar = (*m_sditer).second;
      m_sditer++;
      return 1;
     }
     else return 0;
}
void MdcCalibData::clear(){
     m_xtmap.clear();
     m_t0.clear();
     m_delt0.clear();
     m_qtpar0.clear();
     m_qtpar1.clear();
     m_sdmap.clear();
}
 
}