RootTofCalibDataCnv.cxx

Go to the documentation of this file.

// $Header:
// /bes/bes/BossCvs/Calibration/CalibSvc/CalibROOTCnv/src/cnv/RootTofCalibDataCnv.cxx,v 1.20
// 2019/09/20 07:01:13 sunss Exp $
#include "RootTofCalibDataCnv.h"
#include "CalibData/Tof/TofCalibData.h"
#include "CalibData/Tof/bTofCalibBase.h"
#include "CalibData/Tof/eTofCalibBase.h"
#include "CalibData/Tof/etfBunchCalibBase.h"
#include "CalibData/Tof/etfCalibBase.h"
#include "CalibDataSvc/IInstrumentName.h"
#include "GaudiKernel/MsgStream.h"
 
#include "TDirectory.h"
#include "TFile.h"
#include "TObject.h"
#include "TTree.h"
 
#include "GaudiKernel/DataObject.h"
#include "GaudiKernel/GenericAddress.h"
#include "GaudiKernel/IAddressCreator.h"
#include "GaudiKernel/IConversionSvc.h"
#include "GaudiKernel/IDataProviderSvc.h"
#include "GaudiKernel/IOpaqueAddress.h"
 
#include "CalibDataSvc/ICalibMetaCnvSvc.h"
#include "CalibDataSvc/ICalibRootSvc.h" //maybe
 
// Temporary.  Hope to find a better way to do this
#include "CalibData/CalibModel.h"
using namespace CalibData;
 
const unsigned int nBarAtten  = 8;
const unsigned int nBarSpeed  = 2;
const unsigned int nBarPar    = 10;
const unsigned int nBarParOff = 20;
const unsigned int nBarSigma  = 5;
const unsigned int nBarSigCnt = 10;
 
const unsigned int nBarSigCor = 8;
const unsigned int nBarOffset = 2;
 
const unsigned int nEndAtten = 5;
const unsigned int nEndSpeed = 4;
const unsigned int nEndPar   = 8;
const unsigned int nEndSigma = 5;
 
const unsigned int nEtfSpeed = 4;
const unsigned int nEtfPar   = 20;
const unsigned int nEtfBunch = 4;
 
RootTofCalibDataCnv::RootTofCalibDataCnv( ISvcLocator* svc )
    : RootCalBaseCnv( svc, CLID_Calib_TofCal ) {}
 
const CLID& RootTofCalibDataCnv::objType() const { return CLID_Calib_TofCal; }
 
const CLID& RootTofCalibDataCnv::classID() { return CLID_Calib_TofCal; }
 
StatusCode RootTofCalibDataCnv::i_createObj( const std::string& fname,
                                             DataObject*&       refpObject ) {
 
  MsgStream log( msgSvc(), "RootTofCalibDataCnv" );
  log << MSG::DEBUG << "SetProperty" << endmsg;
 
  // open the file
  StatusCode sc = openRead( fname );
  if ( !sc ) { log << MSG::ERROR << "unable to open files" << endmsg; }
 
  CalibData::bTofCalibBase       bTof;
  CalibData::eTofCalibBase       eTof;
  CalibData::etfCalibBase        etf;
  CalibData::etfBunchCalibBase   etfBunch;
  CalibData::bTofCommonCalibBase bTofCommon;
  CalibData::tofCalibInfoBase    tofinfo;
 
  std::vector<CalibData::bTofCalibBase>       tmpbTof;
  std::vector<CalibData::eTofCalibBase>       tmpeTof;
  std::vector<CalibData::etfCalibBase>        tmpetf;
  std::vector<CalibData::etfBunchCalibBase>   tmpetfBunch;
  std::vector<CalibData::bTofCommonCalibBase> tmpbTofCommon;
  std::vector<CalibData::tofCalibInfoBase>    tofinfoCol;
  // Read in the object
  int cnt;
  // read btoftree ------------------------------------------------------------
 
  double cnvBarAtten[nBarAtten];
  double cnvBarSpeed[nBarSpeed];
  double cnvBarPar1[nBarPar];
  double cnvBarPar2[nBarPar];
  double cnvBarParOff1_bunch0[nBarParOff];
  double cnvBarParOff2_bunch0[nBarParOff];
  double cnvBarParOff1_bunch1[nBarParOff];
  double cnvBarParOff2_bunch1[nBarParOff];
  double cnvBarParOff1_bunch2[nBarParOff];
  double cnvBarParOff2_bunch2[nBarParOff];
  double cnvBarParOff1_bunch3[nBarParOff];
  double cnvBarParOff2_bunch3[nBarParOff];
  double cnvFLeft[nBarSigma];
  double cnvFRight[nBarSigma];
  double cnvFCount[nBarSigCnt];
 
  TTree* btoftree = (TTree*)m_inFile->Get( "BarTofPar" );
 
  char brname[10];
  for ( unsigned int i = 0; i < nBarAtten; i++ )
  {
    sprintf( brname, "Atten%i", i );
    btoftree->SetBranchAddress( brname, &cnvBarAtten[i] );
  }
  for ( unsigned int i = 0; i < nBarSpeed; i++ )
  {
    sprintf( brname, "Speed%i", i );
    btoftree->SetBranchAddress( brname, &cnvBarSpeed[i] );
  }
  for ( unsigned int i = 0; i < nBarPar; i++ )
  {
    sprintf( brname, "P%i", i );
    btoftree->SetBranchAddress( brname, &cnvBarPar1[i] );
  }
  for ( unsigned int i = 0; i < nBarPar; i++ )
  {
    sprintf( brname, "P%i", i + nBarPar );
    btoftree->SetBranchAddress( brname, &cnvBarPar2[i] );
  }
  for ( unsigned int i = 0; i < nBarParOff; i++ )
  {
    sprintf( brname, "Bunch0_Poff%i", i );
    btoftree->SetBranchAddress( brname, &cnvBarParOff1_bunch0[i] );
  }
  for ( unsigned int i = 0; i < nBarParOff; i++ )
  {
    sprintf( brname, "Bunch0_Poff%i", i + nBarParOff );
    btoftree->SetBranchAddress( brname, &cnvBarParOff2_bunch0[i] );
  }
  for ( unsigned int i = 0; i < nBarParOff; i++ )
  {
    sprintf( brname, "Bunch1_Poff%i", i );
    btoftree->SetBranchAddress( brname, &cnvBarParOff1_bunch1[i] );
  }
  for ( unsigned int i = 0; i < nBarParOff; i++ )
  {
    sprintf( brname, "Bunch1_Poff%i", i + nBarParOff );
    btoftree->SetBranchAddress( brname, &cnvBarParOff2_bunch1[i] );
  }
  for ( unsigned int i = 0; i < nBarParOff; i++ )
  {
    sprintf( brname, "Bunch2_Poff%i", i );
    btoftree->SetBranchAddress( brname, &cnvBarParOff1_bunch2[i] );
  }
  for ( unsigned int i = 0; i < nBarParOff; i++ )
  {
    sprintf( brname, "Bunch2_Poff%i", i + nBarParOff );
    btoftree->SetBranchAddress( brname, &cnvBarParOff2_bunch2[i] );
  }
  for ( unsigned int i = 0; i < nBarParOff; i++ )
  {
    sprintf( brname, "Bunch3_Poff%i", i );
    btoftree->SetBranchAddress( brname, &cnvBarParOff1_bunch3[i] );
  }
  for ( unsigned int i = 0; i < nBarParOff; i++ )
  {
    sprintf( brname, "Bunch3_Poff%i", i + nBarParOff );
    btoftree->SetBranchAddress( brname, &cnvBarParOff2_bunch3[i] );
  }
  for ( unsigned int i = 0; i < nBarSigma; i++ )
  {
    sprintf( brname, "FLeft%i", i );
    btoftree->SetBranchAddress( brname, &cnvFLeft[i] );
  }
  for ( unsigned int i = 0; i < nBarSigma; i++ )
  {
    sprintf( brname, "FRight%i", i );
    btoftree->SetBranchAddress( brname, &cnvFRight[i] );
  }
  for ( unsigned int i = 0; i < nBarSigCnt; i++ )
  {
    sprintf( brname, "FCounter%i", i );
    btoftree->SetBranchAddress( brname, &cnvFCount[i] );
  }
 
  for ( cnt = 0; cnt < btoftree->GetEntries(); cnt++ )
  {
    btoftree->GetEntry( cnt );
    bTof.setAtten( cnvBarAtten );
    bTof.setSpeed( cnvBarSpeed );
    bTof.setP1( cnvBarPar1 );
    bTof.setP2( cnvBarPar2 );
    bTof.setPoff1_bunch0( cnvBarParOff1_bunch0 );
    bTof.setPoff2_bunch0( cnvBarParOff2_bunch0 );
    bTof.setPoff1_bunch1( cnvBarParOff1_bunch1 );
    bTof.setPoff2_bunch1( cnvBarParOff2_bunch1 );
    bTof.setPoff1_bunch2( cnvBarParOff1_bunch2 );
    bTof.setPoff2_bunch2( cnvBarParOff2_bunch2 );
    bTof.setPoff1_bunch3( cnvBarParOff1_bunch3 );
    bTof.setPoff2_bunch3( cnvBarParOff2_bunch3 );
    bTof.setFPleft( cnvFLeft );
    bTof.setFPright( cnvFRight );
    bTof.setFPcounter( cnvFCount );
    tmpbTof.push_back( bTof );
  }
 
  // read etoftree
  double cnvEndAtten[nEndAtten];
  double cnvEndSpeed[nEndSpeed];
  double cnvEndPar[nEndPar];
  double cnvEndFPar[nEndSigma];
 
  TTree* etoftree = (TTree*)m_inFile->Get( "EndTofPar" );
 
  char ecname[10];
  for ( unsigned int i = 0; i < nEndAtten; i++ )
  {
    sprintf( ecname, "Atten%i", i );
    etoftree->SetBranchAddress( ecname, &cnvEndAtten[i] );
  }
  for ( unsigned int i = 0; i < nEndSpeed; i++ )
  {
    sprintf( ecname, "Speed%i", i );
    etoftree->SetBranchAddress( ecname, &cnvEndSpeed[i] );
  }
  for ( unsigned int i = 0; i < nEndPar; i++ )
  {
    sprintf( ecname, "P%i", i );
    etoftree->SetBranchAddress( ecname, &cnvEndPar[i] );
  }
  for ( unsigned int i = 0; i < nEndSigma; i++ )
  {
    sprintf( ecname, "FCounter%i", i );
    etoftree->SetBranchAddress( ecname, &cnvEndFPar[i] );
  }
 
  for ( cnt = 0; cnt < etoftree->GetEntries(); cnt++ )
  {
    etoftree->GetEntry( cnt );
    eTof.setAtten( cnvEndAtten );
    eTof.setSpeed( cnvEndSpeed );
    eTof.setP( cnvEndPar );
    eTof.setFP( cnvEndFPar );
    tmpeTof.push_back( eTof );
  }
 
  // read etftree
  double cnvEtfSpeed[nEtfSpeed];
  double cnvEtfPar[nEtfPar];
  double cnvEtfPar1[nEtfPar];
  double cnvEtfPar2[nEtfPar];
 
  if ( NULL != m_inFile->Get( "EtfTofPar" ) )
  {
    TTree* etftree = (TTree*)m_inFile->Get( "EtfTofPar" );
 
    char etfname[10];
    for ( unsigned int i = 0; i < nEtfSpeed; i++ )
    {
      sprintf( etfname, "Speed%i", i );
      etftree->SetBranchAddress( etfname, &cnvEtfSpeed[i] );
    }
    for ( unsigned int i = 0; i < nEtfPar; i++ )
    {
      sprintf( etfname, "P%i", i );
      etftree->SetBranchAddress( etfname, &cnvEtfPar[i] );
    }
    for ( unsigned int i = 0; i < nEtfPar; i++ )
    {
      sprintf( etfname, "P%i", i + nEtfPar );
      etftree->SetBranchAddress( etfname, &cnvEtfPar1[i] );
    }
    for ( unsigned int i = 0; i < nEtfPar; i++ )
    {
      sprintf( etfname, "P%i", i + 2 * nEtfPar );
      etftree->SetBranchAddress( etfname, &cnvEtfPar2[i] );
    }
 
    for ( cnt = 0; cnt < etftree->GetEntries(); cnt++ )
    {
      etftree->GetEntry( cnt );
      etf.setSpeed( cnvEtfSpeed );
      etf.setP( cnvEtfPar );
      etf.setP1( cnvEtfPar1 );
      etf.setP2( cnvEtfPar2 );
      tmpetf.push_back( etf );
    }
  }
 
  // read etfbunchtree
  double cnvEtfBunchP[nEtfBunch];
 
  if ( NULL != m_inFile->Get( "EtfTofBunch" ) )
  {
    TTree* etfbunchtree = (TTree*)m_inFile->Get( "EtfTofBunch" );
 
    char etfbunchname[10];
    for ( unsigned int i = 0; i < nEtfBunch; i++ )
    {
      sprintf( etfbunchname, "pbunch%i", i );
      etfbunchtree->SetBranchAddress( etfbunchname, &cnvEtfBunchP[i] );
    }
 
    int entries = etfbunchtree->GetEntries();
    for ( cnt = 0; cnt < entries; cnt++ )
    {
      etfbunchtree->GetEntry( cnt );
      etfBunch.setPBunch( cnvEtfBunchP );
      tmpetfBunch.push_back( etfBunch );
    }
  }
 
  // read bTofCommonCalibBase
  double cnvBarSigCor[nBarSigCor];
  double cnvBarOffset[nBarOffset];
  TTree* btofcommontree = (TTree*)m_inFile->Get( "BarTofParCommon" );
  for ( unsigned int i = 0; i < nBarSigCor; i++ )
  {
    sprintf( brname, "sigmaCorr%i", i );
    btofcommontree->SetBranchAddress( brname, &cnvBarSigCor[i] );
  }
  for ( unsigned int i = 0; i < nBarOffset; i++ )
  {
    sprintf( brname, "t0offset%i", i );
    btofcommontree->SetBranchAddress( brname, &cnvBarOffset[i] );
  }
 
  int entries = btofcommontree->GetEntries();
  for ( cnt = 0; cnt < entries; cnt++ )
  {
    btofcommontree->GetEntry( cnt );
    bTofCommon.setSigmaCorr( cnvBarSigCor );
    bTofCommon.setOffset( cnvBarOffset );
    tmpbTofCommon.push_back( bTofCommon );
  }
 
  int m_run1, m_run2, m_version;
  int m_qCorr, m_qElec, m_misLable;
  int m_tofidEast[5], m_tofidWest[5], m_tofidEndcap[5];
  int m_runFrom, m_runTo, m_eventFrom, m_eventTo;
 
  TTree* CalibInfo = (TTree*)m_inFile->Get( "CalibInfo" );
  CalibInfo->SetBranchAddress( "Run1", &m_run1 );
  CalibInfo->SetBranchAddress( "Run2", &m_run2 );
  CalibInfo->SetBranchAddress( "Version", &m_version );
  CalibInfo->SetBranchAddress( "ebrId0", &m_tofidEast[0] );
  CalibInfo->SetBranchAddress( "ebrId1", &m_tofidEast[1] );
  CalibInfo->SetBranchAddress( "ebrId2", &m_tofidEast[2] );
  CalibInfo->SetBranchAddress( "ebrId3", &m_tofidEast[3] );
  CalibInfo->SetBranchAddress( "ebrId4", &m_tofidEast[4] );
  CalibInfo->SetBranchAddress( "ecId0", &m_tofidEndcap[0] );
  CalibInfo->SetBranchAddress( "ecId1", &m_tofidEndcap[1] );
  CalibInfo->SetBranchAddress( "ecId2", &m_tofidEndcap[2] );
  CalibInfo->SetBranchAddress( "ecId3", &m_tofidEndcap[3] );
  CalibInfo->SetBranchAddress( "ecId4", &m_tofidEndcap[4] );
  CalibInfo->SetBranchAddress( "wbrId0", &m_tofidWest[0] );
  CalibInfo->SetBranchAddress( "wbrId1", &m_tofidWest[1] );
  CalibInfo->SetBranchAddress( "wbrId2", &m_tofidWest[2] );
  CalibInfo->SetBranchAddress( "wbrId3", &m_tofidWest[3] );
  CalibInfo->SetBranchAddress( "wbrId4", &m_tofidWest[4] );
  CalibInfo->SetBranchAddress( "misLable", &m_misLable );
  CalibInfo->SetBranchAddress( "qCorr", &m_qCorr );
  CalibInfo->SetBranchAddress( "qElec", &m_qElec );
  if ( CalibInfo->GetBranchStatus( "runFrom" ) )
  {
    CalibInfo->SetBranchAddress( "runFrom", &m_runFrom );
    CalibInfo->SetBranchAddress( "runTo", &m_runTo );
    CalibInfo->SetBranchAddress( "eventFrom", &m_eventFrom );
    CalibInfo->SetBranchAddress( "eventTo", &m_eventTo );
  }
  else
  {
    m_runFrom   = -1;
    m_runTo     = -1;
    m_eventFrom = -1;
    m_eventTo   = -1;
  }
 
  entries = CalibInfo->GetEntries();
  for ( cnt = 0; cnt < entries; cnt++ )
  {
    CalibInfo->GetEntry( cnt );
    tofinfo.setRunBegin( m_run1 );
    tofinfo.setRunEnd( m_run2 );
    tofinfo.setVersion( m_version );
    tofinfo.setQCorr( m_qCorr );
    tofinfo.setQElec( m_qElec );
    tofinfo.setMisLable( m_misLable );
    tofinfo.setBrEast( m_tofidEast );
    tofinfo.setBrWest( m_tofidWest );
    tofinfo.setEndcap( m_tofidEndcap );
    tofinfo.setRunFrom( m_runFrom );
    tofinfo.setRunTo( m_runTo );
    tofinfo.setEventFrom( m_eventFrom );
    tofinfo.setEventTo( m_eventTo );
    tofinfoCol.push_back( tofinfo );
  }
 
  CalibData::TofCalibData* tmpObject = new CalibData::TofCalibData(
      tmpbTof, tmpbTofCommon, tmpeTof, tmpetf, tmpetfBunch, tofinfoCol );
 
  refpObject = tmpObject;
 
  return StatusCode::SUCCESS;
}
 
StatusCode RootTofCalibDataCnv::createRoot( const std::string&     fname,
                                            CalibData::CalibBase1* pTDSObj ) {
  MsgStream log( msgSvc(), "RootTofCalibDataCnv" );
 
  // Open the file, create the branch
  StatusCode sc = openWrite( fname );
  if ( !sc ) { log << MSG::ERROR << "unable to open files" << endmsg; }
  // write the Data in the TCDS to RootFile
  int                      j;
  CalibData::TofCalibData* btof = dynamic_cast<CalibData::TofCalibData*>( pTDSObj );
 
  // write  btoftree----------------------------------------------------------------
  double cnvBarAtten[nBarAtten];
  double cnvBarSpeed[nBarSpeed];
  double cnvBarPar1[nBarPar];
  double cnvBarPar2[nBarPar];
  double cnvBarParOff1_bunch0[nBarParOff];
  double cnvBarParOff2_bunch0[nBarParOff];
  double cnvBarParOff1_bunch1[nBarParOff];
  double cnvBarParOff2_bunch1[nBarParOff];
  double cnvBarParOff1_bunch2[nBarParOff];
  double cnvBarParOff2_bunch2[nBarParOff];
  double cnvBarParOff1_bunch3[nBarParOff];
  double cnvBarParOff2_bunch3[nBarParOff];
  double cnvFLeft[nBarSigma];
  double cnvFRight[nBarSigma];
  double cnvFCount[nBarSigCnt];
 
  char   brname[8], ibrname[8];
  TTree* btoftree = new TTree( "BarTofPar", "BarTofPar" );
  for ( unsigned int i = 0; i < nBarAtten; i++ )
  {
    sprintf( brname, "Atten%i", i );
    sprintf( ibrname, "Atten%i/D", i );
    btoftree->Branch( brname, &cnvBarAtten[i], ibrname );
  }
  for ( unsigned int i = 0; i < nBarSpeed; i++ )
  {
    sprintf( brname, "Speed%i", i );
    sprintf( ibrname, "Speed%i/D", i );
    btoftree->Branch( brname, &cnvBarSpeed[i], ibrname );
  }
  for ( unsigned int i = 0; i < nBarPar; i++ )
  {
    sprintf( brname, "P%i", i );
    sprintf( ibrname, "P%i/D", i );
    btoftree->Branch( brname, &cnvBarPar1[i], ibrname );
  }
  for ( unsigned int i = 0; i < nBarPar; i++ )
  {
    sprintf( brname, "P%i", i + nBarPar );
    sprintf( ibrname, "P%i/D", i + nBarPar );
    btoftree->Branch( brname, &cnvBarPar2[i], ibrname );
  }
  for ( unsigned int i = 0; i < nBarParOff; i++ )
  {
    sprintf( brname, "Bunch0_Poff%i", i );
    sprintf( ibrname, "Bunch0_Poff%i/D", i );
    btoftree->Branch( brname, &cnvBarParOff1_bunch0[i], ibrname );
  }
  for ( unsigned int i = 0; i < nBarParOff; i++ )
  {
    sprintf( brname, "Bunch0_Poff%i", i + nBarParOff );
    sprintf( ibrname, "Bunch0_Poff%i/D", i + nBarParOff );
    btoftree->Branch( brname, &cnvBarParOff2_bunch0[i], ibrname );
  }
  for ( unsigned int i = 0; i < nBarParOff; i++ )
  {
    sprintf( brname, "Bunch1_Poff%i", i );
    sprintf( ibrname, "Bunch1_Poff%i/D", i );
    btoftree->Branch( brname, &cnvBarParOff1_bunch1[i], ibrname );
  }
  for ( unsigned int i = 0; i < nBarParOff; i++ )
  {
    sprintf( brname, "Bunch1_Poff%i", i + nBarParOff );
    sprintf( ibrname, "Bunch1_Poff%i/D", i + nBarParOff );
    btoftree->Branch( brname, &cnvBarParOff2_bunch1[i], ibrname );
  }
  for ( unsigned int i = 0; i < nBarParOff; i++ )
  {
    sprintf( brname, "Bunch2_Poff%i", i );
    sprintf( ibrname, "Bunch2_Poff%i/D", i );
    btoftree->Branch( brname, &cnvBarParOff1_bunch2[i], ibrname );
  }
  for ( unsigned int i = 0; i < nBarParOff; i++ )
  {
    sprintf( brname, "Bunch2_Poff%i", i + nBarParOff );
    sprintf( ibrname, "Bunch2_Poff%i/D", i + nBarParOff );
    btoftree->Branch( brname, &cnvBarParOff2_bunch2[i], ibrname );
  }
  for ( unsigned int i = 0; i < nBarParOff; i++ )
  {
    sprintf( brname, "Bunch3_Poff%i", i );
    sprintf( ibrname, "Bunch3_Poff%i/D", i );
    btoftree->Branch( brname, &cnvBarParOff1_bunch3[i], ibrname );
  }
  for ( unsigned int i = 0; i < nBarParOff; i++ )
  {
    sprintf( brname, "Bunch3_Poff%i", i + nBarParOff );
    sprintf( ibrname, "Bunch3_Poff%i/D", i + nBarParOff );
    btoftree->Branch( brname, &cnvBarParOff2_bunch3[i], ibrname );
  }
  for ( unsigned int i = 0; i < nBarSigma; i++ )
  {
    sprintf( brname, "FLeft%i", i );
    sprintf( ibrname, "FLeft%i/D", i );
    btoftree->Branch( brname, &cnvFLeft[i], ibrname );
  }
  for ( unsigned int i = 0; i < nBarSigma; i++ )
  {
    sprintf( brname, "FRight%i", i );
    sprintf( ibrname, "FRight%i/D", i );
    btoftree->Branch( brname, &cnvFRight[i], ibrname );
  }
  for ( unsigned int i = 0; i < nBarSigCnt; i++ )
  {
    sprintf( brname, "FCounter%i", i );
    sprintf( ibrname, "FCounter%i/D", i );
    btoftree->Branch( brname, &cnvFCount[i], ibrname );
  }
 
  for ( int i = 0; i < 176; i++ )
  {
    for ( j = 0; j < static_cast<int>( nBarAtten ); j++ )
    { cnvBarAtten[j] = btof->getBTofAtten( i, j ); }
    for ( j = 0; j < static_cast<int>( nBarSpeed ); j++ )
    { cnvBarSpeed[j] = btof->getBTofSpeed( i, j ); }
    for ( j = 0; j < static_cast<int>( nBarPar ); j++ )
    {
      cnvBarPar1[j] = btof->getBTofPleft( i, j );
      cnvBarPar2[j] = btof->getBTofPright( i, j );
    }
    for ( j = 0; j < static_cast<int>( nBarParOff ); j++ )
    {
      cnvBarParOff1_bunch0[j] = btof->getBTofPoffleft_bunch0( i, j );
      cnvBarParOff2_bunch0[j] = btof->getBTofPoffright_bunch0( i, j );
    }
    for ( j = 0; j < static_cast<int>( nBarParOff ); j++ )
    {
      cnvBarParOff1_bunch1[j] = btof->getBTofPoffleft_bunch1( i, j );
      cnvBarParOff2_bunch1[j] = btof->getBTofPoffright_bunch1( i, j );
    }
    for ( j = 0; j < static_cast<int>( nBarParOff ); j++ )
    {
      cnvBarParOff1_bunch2[j] = btof->getBTofPoffleft_bunch2( i, j );
      cnvBarParOff2_bunch2[j] = btof->getBTofPoffright_bunch2( i, j );
    }
    for ( j = 0; j < static_cast<int>( nBarParOff ); j++ )
    {
      cnvBarParOff1_bunch3[j] = btof->getBTofPoffleft_bunch3( i, j );
      cnvBarParOff2_bunch3[j] = btof->getBTofPoffright_bunch3( i, j );
    }
    for ( j = 0; j < static_cast<int>( nBarSigma ); j++ )
    {
      cnvFLeft[j]  = btof->getBTofFleft( i, j );
      cnvFRight[j] = btof->getBTofFright( i, j );
    }
    for ( j = 0; j < static_cast<int>( nBarSigCnt ); j++ )
    { cnvFCount[j] = btof->getBTofFcounter( i, j ); }
    btoftree->Fill();
  }
 
  // write  etoftree----------------------------------------------------------------
  double cnvEndAtten[nEndAtten];
  double cnvEndSpeed[nEndSpeed];
  double cnvEndPar[nEndPar];
  double cnvEndFPar[nEndSigma];
 
  char   ecname[8], iecname[8];
  TTree* etoftree = new TTree( "EndTofPar", "EndTofPar" );
  for ( unsigned int i = 0; i < nEndAtten; i++ )
  {
    sprintf( ecname, "Atten%i", i );
    sprintf( iecname, "Atten%i/D", i );
    etoftree->Branch( ecname, &cnvEndAtten[i], iecname );
  }
  for ( unsigned int i = 0; i < nEndSpeed; i++ )
  {
    sprintf( ecname, "Speed%i", i );
    sprintf( iecname, "Speed%i/D", i );
    etoftree->Branch( ecname, &cnvEndSpeed[i], iecname );
  }
  for ( unsigned int i = 0; i < nEndPar; i++ )
  {
    sprintf( ecname, "P%i", i );
    sprintf( iecname, "P%i/D", i );
    etoftree->Branch( ecname, &cnvEndPar[i], iecname );
  }
  for ( unsigned int i = 0; i < nEndSigma; i++ )
  {
    sprintf( ecname, "FCounter%i", i );
    sprintf( iecname, "FCounter%i/D", i );
    etoftree->Branch( ecname, &cnvEndFPar[i], iecname );
  }
 
  for ( int i = 0; i < 96; i++ )
  {
    for ( j = 0; j < static_cast<int>( nEndAtten ); j++ )
    { cnvEndAtten[j] = btof->getETofAtten( i, j ); }
    for ( j = 0; j < static_cast<int>( nEndSpeed ); j++ )
    { cnvEndSpeed[j] = btof->getETofSpeed( i, j ); }
    for ( j = 0; j < static_cast<int>( nEndPar ); j++ )
    { cnvEndPar[j] = btof->getETofP( i, j ); }
    for ( j = 0; j < static_cast<int>( nEndSigma ); j++ )
    { cnvEndFPar[j] = btof->getETofFP( i, j ); }
    etoftree->Fill();
  }
 
  // write  etftree----------------------------------------------------------------
  double cnvEtfSpeed[nEtfSpeed];
  double cnvEtfPar[nEtfPar];
  double cnvEtfPar1[nEtfPar];
  double cnvEtfPar2[nEtfPar];
 
  char   etfname[8], ietfname[8];
  TTree* etftree = new TTree( "EtfTofPar", "EtfTofPar" );
  for ( unsigned int i = 0; i < nEtfSpeed; i++ )
  {
    sprintf( etfname, "Speed%i", i );
    sprintf( ietfname, "Speed%i/D", i );
    etftree->Branch( etfname, &cnvEtfSpeed[i], ietfname );
  }
  for ( unsigned int i = 0; i < nEtfPar; i++ )
  {
    sprintf( etfname, "P%i", i );
    sprintf( ietfname, "P%i/D", i );
    etftree->Branch( etfname, &cnvEtfPar[i], ietfname );
  }
  for ( unsigned int i = 0; i < nEtfPar; i++ )
  {
    sprintf( etfname, "P%i", i + nEtfPar );
    sprintf( ietfname, "P%i/D", i + nEtfPar );
    etftree->Branch( etfname, &cnvEtfPar1[i], ietfname );
  }
  for ( unsigned int i = 0; i < nEtfPar; i++ )
  {
    sprintf( etfname, "P%i", i + 2 * nEtfPar );
    sprintf( ietfname, "P%i/D", i + 2 * nEtfPar );
    etftree->Branch( etfname, &cnvEtfPar2[i], ietfname );
  }
 
  for ( int i = 0; i < 72; i++ )
  {
    for ( int k = 0; k < 12; k++ )
    {
      for ( j = 0; j < static_cast<int>( nEtfSpeed ); j++ )
      { cnvEtfSpeed[j] = btof->getEtfSpeed( i, k, j ); }
      for ( j = 0; j < static_cast<int>( nEtfPar ); j++ )
      {
        cnvEtfPar[j]  = btof->getEtfPcombine( i, k, j );
        cnvEtfPar1[j] = btof->getEtfPleft( i, k, j );
        cnvEtfPar2[j] = btof->getEtfPright( i, k, j );
      }
      etftree->Fill();
    }
  }
 
  // write  etftree----------------------------------------------------------------
  double cnvEtfBunchP[nEtfBunch];
 
  char   etfbunchname[8], ietfbunchname[8];
  TTree* etfbunchtree = new TTree( "EtfTofBunch", "EtfTofBunch" );
  for ( unsigned int i = 0; i < nEtfBunch; i++ )
  {
    sprintf( etfbunchname, "pbunch%i", i );
    sprintf( ietfbunchname, "pbunch%i/D", i );
    etfbunchtree->Branch( etfbunchname, &cnvEtfBunchP[i], ietfbunchname );
  }
 
  for ( j = 0; j < static_cast<int>( nEtfBunch ); j++ )
  { cnvEtfBunchP[j] = btof->getEtfPBunch( 0, j ); }
  etfbunchtree->Fill();
 
  //  write all the trees
  btoftree->Write();
  etoftree->Write();
  etftree->Write();
  etfbunchtree->Write();
  delete btoftree;
  delete etoftree;
  delete etftree;
  delete etfbunchtree;
  closeWrite();
  log << MSG::INFO << "successfully create RootFile" << endmsg;
 
  return sc;
}