TofRawDataProvider.cxx

Go to the documentation of this file.

#include "GaudiKernel/Kernel.h"
#include "GaudiKernel/IInterface.h"
#include "GaudiKernel/IIncidentSvc.h"
#include "GaudiKernel/Incident.h"
#include "GaudiKernel/IIncidentListener.h"
#include "GaudiKernel/StatusCode.h"
#include "GaudiKernel/SvcFactory.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/IDataProviderSvc.h"
#include "GaudiKernel/SmartDataPtr.h"
#include "GaudiKernel/DataSvc.h"
#include "EventModel/EventHeader.h"
#include "Identifier/Identifier.h"
#include "Identifier/TofID.h"
#include "TofRawEvent/TofDigi.h"
#include "TofCaliSvc/ITofCaliSvc.h"
#include "TofQCorrSvc/ITofQCorrSvc.h"
#include "TofQElecSvc/ITofQElecSvc.h"
#include "RawDataProviderSvc/TofRawDataProvider.h"
 
#include "RawEvent/RawDataUtil.h"
 
#include <cmath>
// tianhl for mt
#include "GaudiKernel/ThreadGaudi.h"
// tianhl for mt
 
using namespace std;
 
ITofCaliSvc* tofCaliSvc;
ITofQCorrSvc* tofQCorrSvc;
ITofQElecSvc* tofQElecSvc;
 
TofRawDataProvider::TofRawDataProvider() :
  RawDataProviderBase(),
  m_hasFilled(0)
{}
 
 
TofRawDataProvider::TofRawDataProvider( const char* name ) :
  RawDataProviderBase( name ),
  m_hasFilled(0)
{}
 
 
TofRawDataProvider::~TofRawDataProvider() {
  IterTofDataMap iter = m_tofDataMap.begin();
  for( ; iter != m_tofDataMap.end(); iter++ ) {
    delete (*iter).second;
  }
  m_tofDataMap.clear();
 
  m_tofDataMapOnline.clear();
  m_tofDataVectorOnline.clear();
  m_tofDataMapEstime.clear();
  m_tofDataVectorEstime.clear();
  m_tofDataMapTof.clear();
  m_tofDataVectorTof.clear();
  m_tofDataMapEmc.clear();
  m_tofDataVectorEmc.clear();
 
  return;
}
 
 
StatusCode TofRawDataProvider::initialize( bool mode, ISvcLocator* pSvcLoc, IMessageSvc* pMsg ) {
  MsgStream log( m_msgSvc, m_name );
  RawDataProviderBase::initialize( pSvcLoc,pMsg );
 
  //Get TOF Calibtration Service
  StatusCode scc = m_svcLocator->service("TofCaliSvc", tofCaliSvc);
  if (scc ==  StatusCode::SUCCESS) { 
    log << MSG::INFO << "TofRec Get Calibration Service Sucessfully!" << endreq;
  } else if(!mode) {
    log << MSG::ERROR << "TofRec Get Calibration Service Failed !" << endreq;
    return StatusCode::FAILURE;
  }
 
  //Get TOF Q Correction Service
  StatusCode scq = m_svcLocator->service("TofQCorrSvc", tofQCorrSvc);
  if (scq ==  StatusCode::SUCCESS) { 
    log << MSG::INFO << "TofRec Get Q Correction Service Sucessfully!" << endreq;
  } else if(!mode){
    log << MSG::ERROR << "TofRec Get Q Correction Service Failed !" << endreq;
    return StatusCode::FAILURE;
  }
 
  //Get TOF Q Electronics Saturation Service
  StatusCode sce = m_svcLocator->service("TofQElecSvc", tofQElecSvc);
  if (sce ==  StatusCode::SUCCESS) { 
    log << MSG::INFO << "TofRec Get Q Electronics Calibration Service Sucessfully!" << endreq;
  } else if(!mode){
    log << MSG::ERROR << "TofRec Get Q Electronics Calibration Service Failed !" << endreq;
    return StatusCode::FAILURE;
  }
 
  return StatusCode::SUCCESS; 
}
 
 
void TofRawDataProvider::handle(const Incident& inc){
  MsgStream log( m_msgSvc, m_name );
  log << MSG::DEBUG << "handle: " << inc.type() << endreq;
  if ( inc.type() == "BeginEvent" ){
    log << MSG::DEBUG << "Begin Event" << endreq;
    IterTofDataMap iter = m_tofDataMap.begin();
    for( ; iter != m_tofDataMap.end(); iter++ ) {
      delete (*iter).second;
    }
    IterTofDataMap iterOnline = m_tofDataMapOnline.begin();
    for( ; iterOnline != m_tofDataMapOnline.end(); iterOnline++ ) {
      delete (*iterOnline).second;
    }
    m_tofDataMap.clear();
    m_tofDataMapOnline.clear();
    m_tofDataVectorOnline.clear();
    m_tofDataMapEstime.clear();
    m_tofDataVectorEstime.clear();
    m_tofDataMapTof.clear();
    m_tofDataVectorTof.clear();
    m_tofDataMapEmc.clear();
    m_tofDataVectorEmc.clear();
 
    m_hasFilled = false;
  }
  return;
}
 
 
// tofDataMapFull includes all possible combinations.
void TofRawDataProvider::tofDataMapFull() {
 
  MsgStream log(m_msgSvc, m_name);
  log << MSG::INFO << "TofRawDataProvider::tofDataMapFull()" << endreq;
 
  if( !m_tofDataMap.empty() ) {
 
    log << MSG::WARNING << "TofRawDataProvider::tofDataMapFull: m_tofDataMap is NOT empty!" << endreq;
 
    IterTofDataMap iter = m_tofDataMap.begin();
    for( ; iter != m_tofDataMap.end(); iter++ ) {
      delete (*iter).second;
    }
    m_tofDataMap.clear();
  }
   // tianhl for mt
    std::string evtDataSvc_name("EventDataSvc");
    if(isGaudiThreaded(m_name)){
      evtDataSvc_name += getGaudiThreadIDfromName(m_name);
    }
    // tianhl for mt
 
  IDataProviderSvc* eventSvc;
  StatusCode sc = m_svcLocator->service( evtDataSvc_name.c_str(),eventSvc, true );
  if( !sc.isSuccess() ) {
    log << MSG::FATAL << "TofRawDataProvider::tofDataMapFull: ERROR Could not load EventDataSvc" << endreq;
    return;
  }
 
  // Event Header
  SmartDataPtr<Event::EventHeader> eventHeader(eventSvc,"/Event/EventHeader");
  if( !eventHeader ) {
    log << MSG::FATAL << "TofRawDataProvider could not find Event Header!" << endreq;
    return;
  }
  int run = eventHeader->runNumber();
 
  bool mc = ( run < 0 ) ? true:false;
  bool qCorr = bool( tofCaliSvc->QCorr() );
  qCorr = qCorr && (!mc);
  bool qElec = bool( tofCaliSvc->QElec() );
  bool misLable = bool( tofCaliSvc->MisLable() );
  misLable = ( !mc && misLable );
  vector<int> deadId;
 
  int identmp = -1;
  for( unsigned int i=0; i<5; i++ ) {
    identmp = tofCaliSvc->BrEast(i);
    //std::cout << "TofRawDataProvider      DeadID: identmp_breast   "  << identmp <<  std::endl; 
    if( identmp != 0x2fffffff ) {
      //std::cout << "TofRawDataProvider    deadid pushback "   << identmp <<  std::endl;
      deadId.push_back( identmp );
       }
    identmp = tofCaliSvc->BrWest(i);
    //std::cout << "TofRawDataProvider      DeadID: identmp_brwest   "  << identmp <<  std::endl;
 
    if( identmp != 0x2fffffff ) {
      //std::cout << "TofRawDataProvider    deadid pushback "   << identmp <<  std::endl;
      deadId.push_back( identmp );
    }
    identmp = tofCaliSvc->Endcap(i);
    //std::cout << "TofRawDataProvider      DeadID: identmp_endcap   "  << identmp <<  std::endl;
    if( identmp != 0x2fffffff ) {
      //std::cout << "TofRawDataProvider    deadid pushback "   << identmp <<  std::endl;
      deadId.push_back( identmp );
    }
  }
 
  // get TDS data in a common class
  SmartDataPtr<TofDigiCol> tofDigiCol(eventSvc,"/Event/Digi/TofDigiCol");
  if( !tofDigiCol ) {
    log << MSG::WARNING << "TofRawDataProvider::tofDataMapFull: Could not find Tof Digi Data!" << endreq;
    return;
  } 
 
  int qnumber = 0;
  int tnumber = 0;
 
  std::vector<Adc*> chargeVec;
  std::vector<Tdc*> timeVec;
 
 
 
 
  TofDigiCol::iterator iter = tofDigiCol->begin();
 
for( ; iter != tofDigiCol->end(); iter++ ) {
 
 
 
  Identifier iden_help = (*iter)->identify();
  int  barrel_ec_help  = TofID::barrel_ec(iden_help);
 
 
  if(barrel_ec_help == 0 || barrel_ec_help ==1 || barrel_ec_help==2) // barrel tof or old endcap
    {
 
    unsigned int overflow   = (*iter)->getOverflow();
    // electronics error message
 
 
    if( ( overflow & 0xfe000000 ) == 0xfe000000 ) 
      {
    continue;
      }
 
    // overflow 0x   1 1         1 1            1 1
    //            no Q T / multi Q T / overflow Q T
    
    bool multiQ = ( ( overflow & 0x8 ) != 0 );
    bool multiT = ( ( overflow & 0x4 ) != 0 );
 
    Identifier iden = (*iter)->identify();
 
    unsigned int iden_value = iden.get_value();
 
 
 
 
    bool sameId = false;
    vector<int>::iterator iter_dead = deadId.begin();
    for( ; iter_dead != deadId.end(); iter_dead++ ) {
      if( iden_value == (*iter_dead) ) {
    sameId = true;
    break;
      }
    }
    if( sameId ) {
      if( mc ) { continue; }
      else {
    log << MSG::ERROR << "TofRawDataProvider::tofDataMapFull: Dead Channel Number is not Correct, please check TOF Calibration Constants!" << endreq;
     }
    }
 
 
    int  barrel_ec  = TofID::barrel_ec(iden);
    int  id         = TofID::phi_module(iden);
    int  end        = TofID::end(iden);
 
 
 
    if( ( barrel_ec != 1 ) && ( end == 1 ) ) continue;
    // TofDigiCol 48 and 96 for luminocity
    if( ( barrel_ec != 1 ) && ( id == 48 ) ) continue;
 
    unsigned int tdcChannel = (*iter)->getTimeChannel();
    unsigned int adcChannel = (*iter)->getChargeChannel();
    
 
 
 
    if( ( overflow & 0x2 ) != 0 ) 
      { 
    adcChannel = ( adcChannel | 0x80000 ); 
    }
 
    if( (adcChannel&0x7fffffff) != 0x7fffffff ) {
      
      Adc* adc = new Adc;
      if( qCorr ) adc->setCorr();
      if( qElec ) adc->setElec();
    
      adc->setValue( adcChannel );
      qnumber++;
      adc->setNumber( qnumber );
      chargeVec.push_back( adc );
    }
 
 
 
    if( tdcChannel != 0x7fffffff ) {
      Tdc* tdc = new Tdc;
      tdc->setValue( tdcChannel );
      tnumber++;
      tdc->setNumber( tnumber );
      timeVec.push_back( tdc );
    }
 
 
 
    if( multiT || multiQ ) {
      TofDigiCol::iterator iter_next = iter + 1;
      if( iter_next != tofDigiCol->end() ) {
    unsigned int idenNext_value = ((*iter_next)->identify()).get_value();
    if( iden_value == idenNext_value ) continue;
      }
    }
 
 
 
    qnumber = 0;
    tnumber = 0;
 
    std::vector<PmtData*> pmtDataVec;
    std::vector<Tdc*>::iterator iter_t = timeVec.begin();
    for( ; iter_t != timeVec.end(); iter_t++ ) {
      int tclock = (*iter_t)->clock();
      PmtData* pmt = new PmtData;
      pmt->setTdc( (*iter_t) );
 
      std::vector<Adc*>::iterator iter_q = chargeVec.begin();
      for( ; iter_q != chargeVec.end(); iter_q++ ) {
    int qclock = (*iter_q)->clock();
    
    if( abs(tclock-qclock)<2 ) {
       if( ( pmt->quality() & 0x2 ) != 0 ) {         // pmt has Q
         if( (*iter_q)->value() > pmt->adc() ) {
          pmt->qtimesmm();
          pmt->setAdc( (*iter_q) );
 
          
        }
       }//close pmt has q
      else {                                        // pmt has no Q
        pmt->setAdc( (*iter_q) );
      }
    }
      }
      pmtDataVec.push_back( pmt );
    }
 
 
 
    std::vector<Adc*>::iterator iter_q = chargeVec.begin();
    for( ; iter_q != chargeVec.end(); iter_q++ ) {
      //cout << "(*iter_q)->times() " << (*iter_q)->times() <<endl;
      if( (*iter_q)->times() != -1 )
    {
     continue;
    }
 
      //cout << "Hallo Hier, d.h. times() von charge  !=-1 und barrel_ec ist  "<< barrel_ec_help  <<endl;
      PmtData* pmt = new PmtData;
      pmt->setAdc( (*iter_q) );
      pmtDataVec.push_back( pmt );
    }
 
 
    //std::cout << "TofRawDataProvider     old chargeVec.size()   " << chargeVec.size() << std::endl;
    //std::cout << "TofRawDataProvider     old timeVec.size()     " << timeVec.size() << std::endl;
    chargeVec.clear();
    timeVec.clear();
 
    unsigned int iden_value_key      = (iden_value & 0xfffffffe );
    if( ( barrel_ec == 0 ) || ( barrel_ec == 2 ) ) {
      std::vector<PmtData*>::iterator iter_pmt = pmtDataVec.begin();
      for( ; iter_pmt != pmtDataVec.end(); iter_pmt++ ) {
    TofData* tof = new TofData;
    if( misLable ) { 
      tof->setMisLable(); 
 
      if( barrel_ec==2 && id==42 ) {
        unsigned int iden91 = TofID::getIntID( 2, 0, 91-48, 0 ); 
 
        iden_value_key = iden91; 
      } 
 
      if( barrel_ec==2 && id==43 ) { 
        unsigned int iden90 = TofID::getIntID( 2, 0, 90-48, 0 );
        iden_value_key = iden90;
      }
 
    }
 
 
    tof->setIdentify( iden );
    tof->setForward( (*iter_pmt) );
    tof->setTMatched( true );
    m_tofDataMap.insert( make_pair( iden_value_key, tof ) );
      }
    }
    else if( barrel_ec == 1 ) { //super
      unsigned int count = m_tofDataMap.count( iden_value_key );
      if( count == 0 ) {
    std::vector<PmtData*>::iterator iter_pmt = pmtDataVec.begin();
    for( ; iter_pmt != pmtDataVec.end(); iter_pmt++ ) {
      TofData* tof = new TofData;
      tof->setIdentify( iden );
      if( end == 0 ) {
        tof->setForward( (*iter_pmt) );
      }
      else {
        tof->setBackward( (*iter_pmt) );
      }
      tof->setTMatched( true );
 
      m_tofDataMap.insert( make_pair( iden_value_key, tof ) );
    }
      }
      else {
    pair< IterTofDataMap, IterTofDataMap > range = m_tofDataMap.equal_range( iden_value_key );
    std::vector<PmtData*>::iterator iter_pmt = pmtDataVec.begin();
    for( ; iter_pmt != pmtDataVec.end(); iter_pmt++ ) {
      bool used = false;
      IterTofDataMap iter = range.first;
      for( unsigned int tofDataNumber = 0; tofDataNumber < count; tofDataNumber++, iter++ ) {
        bool matched = false;
        if( end == 0 ) { // east, forward
          if( ( ( (*iter_pmt)->quality() & 0x1 ) != 0 ) && ( ( ((*iter).second)->quality() & 0x1 ) != 0 ) ) { // has east and west time
        double time1 = (*iter_pmt)->tdc();
        double time2 = ((*iter).second)->tdc2();
        matched = ( abs(time1-time2)<timeDiff );
          }
          else if( ( ( (*iter_pmt)->quality() & 0x1 ) == 0 ) && ( ( (*iter_pmt)->quality() & 0x2 ) != 0 ) && ( ( ((*iter).second)->quality() & 0x1 ) != 0 ) ) { // has east charge and west time
        double time1 = (*iter_pmt)->qclock();
        double time2 = ((*iter).second)->tclock2();
        double time3 = ((*iter).second)->qclock2();
        matched = ( ( abs(time1-time2)<=tClockDiff ) || ( abs(time1-time3)<=tClockDiff ) );
          }
          else if( ( ( (*iter_pmt)->quality() & 0x1 ) != 0 ) && ( ( ((*iter).second)->quality() & 0x1 ) == 0 ) && ( ( ((*iter).second)->quality() & 0x2 ) != 0 ) ) { // has east time and west charge
        double time1 = (*iter_pmt)->tclock();
        double time2 = (*iter_pmt)->qclock();
        double time3 = ((*iter).second)->qclock2();
        matched = ( ( abs(time1-time3)<=tClockDiff ) || ( abs(time2-time3)<=tClockDiff ) );
          }
          else if( ( ( (*iter_pmt)->quality() & 0x1 ) == 0 ) && ( ( (*iter_pmt)->quality() & 0x2 ) != 0 ) && ( ( ((*iter).second)->quality() & 0x1 ) == 0 ) && ( ( ((*iter).second)->quality() & 0x2 ) != 0 ) ) { // has east and west charge
        double time1 = (*iter_pmt)->qclock();
        double time2 = ((*iter).second)->qclock2();
        matched = ( abs(time1-time2)<=tClockDiff );
          }
          
          if( matched ) {
        used = true;
        if( ( ( (*iter).second)->quality() & 0xc ) == 0 ) {
          ((*iter).second)->setForward( (*iter_pmt) );
        }
        else {
          TofData* tof = new TofData;
          tof->setIdentify( iden );
          tof->setForward( (*iter_pmt) );
          tof->setBackward( ((*iter).second)->backward() );
                  
          m_tofDataMap.insert( make_pair( iden_value_key, tof ) );
        }
          }
          
        }//close if(end==0)
        else { // west, backward
          if( ( ( (*iter_pmt)->quality() & 0x1 ) != 0 ) && ( ( ((*iter).second)->quality() & 0x4 ) != 0 ) ) { // has west and east time
        double time1 = (*iter_pmt)->tdc();
        double time2 = ((*iter).second)->tdc1();
        matched = ( abs(time1-time2)<timeDiff );
          }
          else if( ( ( (*iter_pmt)->quality() & 0x1 ) == 0 ) && ( ( (*iter_pmt)->quality() & 0x2 ) != 0 ) && ( ( ((*iter).second)->quality() & 0x4 ) != 0 ) ) { // has west charge and east time 
        double time1 = (*iter_pmt)->qclock();
        double time2 = ((*iter).second)->tclock1();
        double time3 = ((*iter).second)->qclock1();
        matched = ( ( abs(time1-time2)<=tClockDiff ) || ( abs(time1-time3)<=tClockDiff ) );
          }
          else if( ( ( (*iter_pmt)->quality() & 0x1 ) != 0 ) && ( ( ((*iter).second)->quality() & 0x4 ) == 0 ) && ( ( ((*iter).second)->quality() & 0x8 ) != 0 ) ) { // has east time and west charge
        double time1 = (*iter_pmt)->tclock();
        double time2 = (*iter_pmt)->qclock();
        double time3 = ((*iter).second)->qclock1();
        matched = ( ( abs(time1-time3)<=tClockDiff ) || ( abs(time2-time3)<=tClockDiff ) );
          }
          else if( ( ( (*iter_pmt)->quality() & 0x1 ) == 0 ) && ( ( (*iter_pmt)->quality() & 0x2 ) != 0 ) && ( ( ((*iter).second)->quality() & 0x4 ) == 0 ) && ( ( ((*iter).second)->quality() & 0x8 ) != 0 ) ) { // has east and west charge
        double time1 = (*iter_pmt)->qclock();
        double time2 = ((*iter).second)->qclock1();
        matched = ( abs(time1-time2)<=tClockDiff );
          }
          
          if( matched ) {
        used = true;
        if( ( ( (*iter).second)->quality() & 0x3 ) == 0 ) {
          ((*iter).second)->setBackward( (*iter_pmt) );
          if( ( ( ((*iter).second)->quality() & 0x5 )==0x5 ) || ( ( ( ((*iter).second)->quality() & 0x5 )==0x4 ) && ( ((*iter).second)->ttimes1()==1 ) ) || ( ( ( ((*iter).second)->quality() & 0x5 )==0x1 ) && ( ((*iter).second)->ttimes2()==1 ) ) ) {
            ((*iter).second)->setTMatched( true );
          }
          else {
            ((*iter).second)->setTMatched( false );
          }
        }
        else {
          TofData* tof = new TofData;
          tof->setIdentify( iden );
          tof->setForward( ((*iter).second)->forward() );
          tof->setBackward( (*iter_pmt) );
          if( ( ( tof->quality() & 0x5 )==0x5 ) || ( ( ( tof->quality() & 0x5 )==0x4 ) && ( tof->ttimes1()==1 ) ) || ( ( ( tof->quality() & 0x5 )==0x1 ) && ( tof->ttimes2()==1 ) ) ) {
            tof->setTMatched( true );
          }
          else {
            tof->setTMatched( false );
          }
          m_tofDataMap.insert( make_pair( iden_value_key, tof ) );
        }
          }
        }
      }
      if( ! used ) {
        TofData* tof = new TofData;
        tof->setIdentify( iden );
        if( end == 0 ) {
          tof->setForward( (*iter_pmt) );
        }
        else {
          tof->setBackward( (*iter_pmt) );
        } 
        tof->setTMatched( true );
        m_tofDataMap.insert( make_pair( iden_value_key, tof ) );
      }
    }
      }
    }
 
    pmtDataVec.clear();
 
  }//close if for barrel/old tof
 
 
  else if(barrel_ec_help==3 || barrel_ec_help==4 || barrel_ec_help==5 || barrel_ec_help==6)//We treat the new MRPC endcap here!
    {
 
 
 
 
      //////////////////////////////////////////////////////
      ////            IMPORTANT                          ///
      /// We use the same datastructure as the olf tof,  ///
      /// even if we do not have PMTs and so on          ///
      //////////////////////////////////////////////////////
 
 
      
    //#Matthias
    Identifier iden = (*iter)->identify();
 
    unsigned int iden_value = iden.get_value();
    //std::cout << "TofRawDataProvider       iden   " << iden << std::endl;
    //std::cout << "TofRawDataProvider       iden.getvalue   " <<iden_value << std::endl;
 
    int  barrel_ec  = TofID::barrel_ec(iden);
    int  id         = TofID::phi_module(iden);
    int  end        = TofID::end(iden);
    //std::cout << "TofRawDataProvider      tofDataMapFull()   partID  |  phi_module  |  end   " << barrel_ec  << "  |  " << id << "  |  "<< end<< std::endl;
 
 
 
    /* Alt
    unsigned int tdcChannel = (*iter)->getTimeChannel();
    unsigned int tdc2Channel = (*iter)->getChargeChannel();
    double multiplicationfactor=1000.; //Convert charge from fc to pc as the time-amplitude correlation is in pc and ps                                                             
 
 
    double timedifference=RawDataUtil::TofTime(tdc2Channel-tdcChannel);
 
    unsigned int adcChannel =0;
    if(timedifference>4.164) adcChannel = exp(-4.50565/(timedifference-16.6653))/0.0812208;   // 4.164 ns = 10 fC; Everything smaller is expected to be noise! ->Deposit charge =0;
 
 
    //With the calculated charge we correct the tdc-channel using the time amplitude correlation.                                                                                   
    double my_chargein_pc = (double)adcChannel/multiplicationfactor;
    double time_of_flight = RawDataUtil::TofTime(tdcChannel);
 
 
    double correction_ta=0;
    if(my_chargein_pc<=2.3)
      {
        correction_ta = 268.81-134.586*pow(my_chargein_pc,1)+40.8928*pow(my_chargein_pc,2)+46.9348*pow(my_chargein_pc,3)-31.754*pow(my_chargein_pc,4)
          +0.62064*pow(my_chargein_pc,5)+1.5672*pow(my_chargein_pc,6);
      }
    else
      {
        correction_ta= 209.332-31.3283*my_chargein_pc;
        if(correction_ta<0)correction_ta=0;
      }
 
    time_of_flight=time_of_flight - correction_ta/1000.; //Divide by 1000 to convert in ns                                                                                          
    tdcChannel = RawDataUtil::TofTimeChannel(time_of_flight); //This is the correceted time of flight. All further calcualtion are performed with the corrected version.      
 
    adcChannel = ( (tdcChannel & 0x07e000) | adcChannel);
    */
 
    unsigned int tdcChannel = (*iter)->getTimeChannel();
    unsigned int tdc2Channel = (*iter)->getChargeChannel();
 
 
    double time_of_flight = RawDataUtil::TofTime(tdcChannel);
    double pulse_length=RawDataUtil::TofTime(tdc2Channel-tdcChannel);
    unsigned int pulse_length_channel_without_correction = tdc2Channel-tdcChannel;
 
    if(pulse_length < 4.164){continue;}   // 4.164 ns = 10 fC; Everything smaller is expected to be noise! ->Deposit charge =0;
 
 
    //Next steps: Do the Time-Amplitude-Correlation: 
    //std::cout << "TofRawDataProvider     pulselength channe no_correction  "<< tdc2Channel-tdcChannel  <<std::endl;
    //std::cout << "TofRawDataProvider     pulselength                       "<< pulse_length  <<std::endl;
    //std::cout << "TofRawDataProvider     tdcChannel no correction   "<< tdcChannel  <<std::endl;
    //std::cout << "TofRawDataProvider     tof  no correction   "<< time_of_flight  <<std::endl;
 
    double correction_ta=0;
    if(pulse_length<=13.9)
      {
        correction_ta= 321.2;
      }
    else 
      {
        correction_ta = 72310.5  -   15162.2*pulse_length  +  1066.06*pulse_length*pulse_length  - 25.0252*pulse_length*pulse_length*pulse_length;
        if(correction_ta<0)correction_ta=0;
      }
    
 
        time_of_flight=time_of_flight - correction_ta/1000.; //Divide by 1000 to convert in ns
        tdcChannel = RawDataUtil::TofTimeChannel(time_of_flight); //This is the correceted time of flight. All further calcualtion are performed with the corrected version.
    //std::cout << "TofRawDataProvider     correction [ns]  "<< correction_ta/1000.  <<std::endl;
    //std::cout << "TofRawDataProvider     tdcChannel correct   "<< tdcChannel  <<std::endl;
    //std::cout << "TofRawDataProvider     tof  correct   "<< time_of_flight  <<std::endl;
 
        
    //std::cout << "TofRawDataProvider      qCorr    |    qElec  " << qCorr    << "    |   " << qElec  << std::endl;
 
 
 
    
 
    //Store HIt information: The time of the first hit is stored in the "TDC" the lenght of the pulse in the "ADC", as the pulse length correspond to charge
      
      Adc* adc = new Adc;
      if( qCorr ) adc->setCorr(); //qCorr =0 for MC
      if( qElec ) adc->setElec(); //qElec= 1 for MC
 
      //std::cout << "TofRawDataProvider      tofDataMapFull()  qCorr   | qElec   " << qCorr << "    |   "<< qElec  << std::endl;
 
 
 
      
      adc->setValue_mrpc( pulse_length_channel_without_correction,(tdcChannel & 0x07e000)  ); //We store the length of the pulse into the ADC information as it corresponds to the charge! The other part is the clock information
      //std::cout << "TofRawDataProvider     adc->number() after initialize   " << adc->number() << std::endl;
      //std::cout << "TofRawDataProvider     adc->getValue   " << adc->channel()  << std::endl;
      qnumber++;
      adc->setNumber(qnumber);//old stuff ?????
      //std::cout << "TofRawDataProvider     adc->number() after set   " << adc->number() << std::endl;
    
     
      Tdc* tdc = new Tdc;
      tdc->setValue( tdcChannel );
      tnumber++;
 
      //std::cout << "TofRawDataProvider     tdc->number() after initialize   " << tdc->number() << std::endl;
      tdc->setNumber( tnumber );//old stuff ?????
      //std::cout << "TofRawDataProvider     tdc->number() after set   " << tdc->number() << std::endl;
 
 
      qnumber = 0;
      tnumber = 0;
 
      //Store the ADC and TDC information into a PMT
      PmtData* pmt = new PmtData;
      pmt->setTdc(tdc);
      pmt->setAdc(adc);
      
      
      chargeVec.clear();//old stuff 
      timeVec.clear();//old stuff
 
 
 
      //Store this data into the TOFDATAMAP !
      unsigned int iden_value_key_mrpc = (iden_value & 0xffffffff ); //The MRPC needs the last bit! 
 
      TofData* tof = new TofData;
      if( misLable ) { tof->setMisLable(); }
      tof->setIdentify( iden );
      tof->setForward( pmt );
      tof->setTMatched( true );//Sof far only MC data -> Match is true
 
      //std::cout << "TofRawDataProvider      tofDataMapFull():  Set the tofdata Endcap! pmt->adcchanel() =  " << (*iter_pmt)->adcChannel() << std::endl;
      //std::cout << "TofRawDataProvider      tofDataMapFull():  Set the tofdata Endcap! tof->adc() =  " << tof->adc() << std::endl;
      //std::cout << "TofRawDataProvider      tofDataMapFull():  Set the tofdata Endcap! tof->tdc() =  " << tof->tdc() << std::endl;
      //std::cout << "TofRawDataProvider      tofDataMapFull():  Set the tofdata Endcap! tof->adc1() =  " << tof->adc1() << std::endl;
      //std::cout << "TofRawDataProvider      tofDataMapFull():  Set the tofdata Endcap! tof->tdc1() =  " << tof->tdc1() << std::endl;
      //std::cout << "TofRawDataProvider      tofDataMapFull():  Set the tofdata Endcap! tof->adc2() =  " << tof->adc2() << std::endl;
      //std::cout << "TofRawDataProvider      tofDataMapFull():  Set the tofdata Endcap! tof->tdc2() =  " << tof->tdc2() << std::endl;
      //std::cout << "TofRawDataProvider      tofDataMapFull():  Set the tofdata Endcap! tof->zadc() =  " << tof->zadc() << std::endl;
      //std::cout << "TofRawDataProvider      tofDataMapFull():  Set the tofdata Endcap! tof->ztdc() =  " << tof->ztdc() << std::endl;
      
 
      m_tofDataMap.insert( make_pair( iden_value_key_mrpc, tof ) );
 
 
 
 
 
    }//close else if new endcap.
 
 
 
}//close for loop
 
 
 
  return;
 
}//close tofDataMapFull
 
 
// onlineMode only T is combined, no Q reserved
TofDataMap& TofRawDataProvider::tofDataMapOnlineMode(uint32_t control) {
 
  MsgStream log(m_msgSvc, m_name);
  log << MSG::INFO << "TofRawDataProvider::tofDataMapOnlineMode()" << endreq;
 
  if( ! m_tofDataMapOnline.empty() ) {
    if(control){
      IterTofDataMap iter = m_tofDataMapOnline.begin();
      for( ; iter != m_tofDataMapOnline.end(); iter++ ) {
    delete (*iter).second;
      }
      m_tofDataMapOnline.clear();
    }
    else return m_tofDataMapOnline;
  }
    // tianhl for mt
    std::string evtDataSvc_name("EventDataSvc");
    if(isGaudiThreaded(m_name)){
      evtDataSvc_name += getGaudiThreadIDfromName(m_name);
    }
    // tianhl for mt
 
  IDataProviderSvc* eventSvc;
  StatusCode sc = m_svcLocator->service( evtDataSvc_name.c_str(), eventSvc, true );
  if( !sc.isSuccess() ) {
    log << MSG::FATAL << "TofRawDataProvider::tofDataMapOnlineMode(): ERROR Could not load EventDataSvc" << endreq;
    return m_tofDataMapOnline;
  }
 
  // get TDS data in a common class
  SmartDataPtr<TofDigiCol> tofDigiCol(eventSvc,"/Event/Digi/TofDigiCol");
  if( !tofDigiCol ) {
    log << MSG::WARNING << "TofRawDataProvider::tofDataMapOnlineMode(): Could not find Tof Digi Data!" << endreq;
    return m_tofDataMapOnline;
  } 
 
  TofDigiCol::iterator iter = tofDigiCol->begin();
  for( ; iter != tofDigiCol->end(); iter++ ) {
    unsigned int overflow   = (*iter)->getOverflow();
    if( ( overflow & 0xfe000000 ) == 0xfe000000 ) continue;
    // overflow 0x   1 1         1 1            1 1
    //            no Q T / multi Q T / overflow Q T
    
    bool noT = ( ( overflow & 0x10 ) != 0 );
    //    bool multiQ = ( ( overflow & 0x8 ) != 0 );
    //    bool multiT = ( ( overflow & 0x4 ) != 0 );
    if( noT ) continue;
    
    Identifier iden = (*iter)->identify();
    unsigned int iden_value = iden.get_value();
    int  barrel_ec  = TofID::barrel_ec(iden);
    int  id         = TofID::phi_module(iden);
    int  end        = TofID::end(iden);
    if( ( barrel_ec != 1 ) && ( end == 1 ) ) continue;
    if( ( barrel_ec != 1 ) && ( id == 48 ) ) continue;
 
    unsigned int tdcChannel = (*iter)->getTimeChannel();
    if( tdcChannel == 0x7fffffff ) continue;
 
    Tdc* tdc = new Tdc;
    tdc->setValue( tdcChannel );
    PmtData* pmt = new PmtData;
    pmt->setTdc( tdc );
 
    unsigned int iden_value_key = ( iden_value & 0xfffffffe );
 
    if( barrel_ec == 0 || barrel_ec == 2 ) {
      TofData* tof = new TofData;
      tof->setIdentify( iden );
      tof->setForward( pmt );
      m_tofDataMapOnline.insert( make_pair( iden_value_key, tof ) );
    }
    else if( barrel_ec == 1 ) {
      unsigned int count = m_tofDataMapOnline.count( iden_value_key );
 
      if( count == 0 ) {
    TofData* tof = new TofData;
    tof->setIdentify( iden );
    if( end == 0 ) {
      tof->setForward( pmt );
    }
    else {
      tof->setBackward( pmt );
    } 
    m_tofDataMapOnline.insert( make_pair( iden_value_key, tof ) );
      }
      else {
    bool used = false;
    pair< IterTofDataMap, IterTofDataMap > range = m_tofDataMapOnline.equal_range( iden_value_key );
    IterTofDataMap iter = range.first;
    for( unsigned int i=0; i<count; i++, iter++ ) {
      if( ( end == 0 ) && ( ( ( (*iter).second)->quality() & 0x1 ) != 0 ) ) {
        double time1 = pmt->tdc();
        double time2 = ((*iter).second)->tdc2();
        if( abs(time1-time2)<timeDiff ) {
          used = true;
          if( ( ( (*iter).second)->quality() & 0x4 ) == 0 ) {
        ((*iter).second)->setForward( pmt );
          }
          else {
        TofData* tof = new TofData;
        tof->setIdentify( iden );
        tof->setForward( pmt );
        tof->setBackward( ((*iter).second)->backward() );
        m_tofDataMapOnline.insert( make_pair( iden_value_key, tof ) );
          }
        }
      }
      else if( ( end == 1 ) && ( ( ((*iter).second)->quality() & 0x4 ) != 0 ) ) {
        double time1 = ((*iter).second)->tdc1();
        double time2 = pmt->tdc();
 
        if( abs(time1-time2)<timeDiff ) {
          used = true;
          if( ( ( (*iter).second)->quality() & 0x1 ) == 0 ) {
        ((*iter).second)->setBackward( pmt );
          }
          else {
        TofData* tof = new TofData;
        tof->setIdentify( iden );
        tof->setForward( ((*iter).second)->forward() );
        tof->setBackward( pmt );
        m_tofDataMapOnline.insert( make_pair( iden_value_key, tof ) );
          }
        }
      }
    }
    if( ! used ) {
      TofData* tof = new TofData;
      tof->setIdentify( iden );
      if( end == 0 ) {
        tof->setForward( pmt );
      }
      else {
        tof->setBackward( pmt );
      } 
      m_tofDataMapOnline.insert( make_pair( iden_value_key, tof ) );
    }
      }
    }
  }
 
  return m_tofDataMapOnline;
 
}
 
 
TofDataVector& TofRawDataProvider::tofDataVectorOnlineMode(uint32_t control) {
 
  MsgStream log(m_msgSvc, m_name);
  log << MSG::INFO << "TofRawDataProvider::tofDataVectorOnlineMode()" << endreq;
  
  if( m_tofDataVectorOnline.size() != 0 ) {
    if(!control) return m_tofDataVectorOnline;
    else         m_tofDataVectorOnline.clear();
  }
  TofDataMap tofRawDataMap = TofRawDataProvider::tofDataMapOnlineMode(control);
  IterTofDataMap iter = tofRawDataMap.begin();
  for( ; iter != tofRawDataMap.end(); iter++ ) {
    m_tofDataVectorOnline.push_back( iter->second );
  }
 
  return m_tofDataVectorOnline;
}
 
 
// Event Start Time, only Time is required.
TofDataMap& TofRawDataProvider::tofDataMapEstime() {
 
  if( !m_hasFilled ) {
    TofRawDataProvider::tofDataMapFull();
    m_hasFilled = true;
  }
 
  if( !m_tofDataMapEstime.empty() ) {
    m_tofDataMapEstime.clear();
  }
 
  IterTofDataMap iter = m_tofDataMap.begin();
  for( ; iter != m_tofDataMap.end(); iter++ ) {
    if( ( ( ((*iter).second)->quality() & 0x5 ) != 0 ) && ((*iter).second)->tmatched() ) {
      m_tofDataMapEstime.insert( make_pair( (*iter).first, (*iter).second ) );
    }
  }
 
  iter = m_tofDataMapEstime.begin();
  while( iter != m_tofDataMapEstime.end() ) {
    unsigned int iden_value = (*iter).first;
    unsigned int count = m_tofDataMapEstime.count( iden_value );
    for( unsigned int i=0; i != count; i++, iter++ ) {
      ((*iter).second)->setTimes( count );
    }
  }
 
  return m_tofDataMapEstime;
}
 
 
TofDataVector& TofRawDataProvider::tofDataVectorEstime() {
  if( m_tofDataVectorEstime.size() != 0 ) return m_tofDataVectorEstime;
 
  TofDataMap tofDataMapEstime = TofRawDataProvider::tofDataMapEstime();
  IterTofDataMap iter = tofDataMapEstime.begin();
  for( ; iter != tofDataMapEstime.end(); iter++ ) {
    m_tofDataVectorEstime.push_back( iter->second );
  }
  return m_tofDataVectorEstime;
}
 
 
// TOF reconstruction, both time and charge are required
TofDataMap& TofRawDataProvider::tofDataMapTof( double estime ) {
 
  if( !m_hasFilled ) {
    TofRawDataProvider::tofDataMapFull();
    m_hasFilled = true;
  }
 
  if( !m_tofDataMapTof.empty() ) {
    m_tofDataMapTof.clear();
  }
 
  IterTofDataMap iter = m_tofDataMap.begin();
  for( ; iter != m_tofDataMap.end(); iter++ ) {
    if( ( ( ((*iter).second)->quality() & 0xc ) == 0xc ) || ( ( ((*iter).second)->quality() & 0x3 ) == 0x3 ) ) {
 
      if( estime > 1.0e-6 ) {
    double tdc1 = ((*iter).second)->tdc1();
    double tdc2 = ((*iter).second)->tdc2();
    bool  forward = ( ( ( estime - tdc1 ) < timeBackward ) && ( ( tdc1 - estime ) < timeForward ) );
    bool backward = ( ( ( estime - tdc2 ) < timeBackward ) && ( ( tdc2 - estime ) < timeForward ) );
    
    //std::cout << "TofRawDataProvider       tdc1: "  << tdc1 << std::endl;
    //std::cout << "TofRawDataProvider       tdc2: " << tdc1<< std::endl;
    //std::cout << "TofRawDataProvider      (est - tdc1)  |   (tdc1 - est):  "  <<  estime - tdc1 << "  |  " <<  tdc1 - estime<< std::endl;
    //std::cout << "TofRawDataProvider      (est - tdc2)  |   (tdc2 - est):  "  <<  estime - tdc2 << "  |  " <<  tdc2 - estime<< std::endl;
    //std::cout << "TofRawDataProvider    bool  (est - tdc1)< 24  &&   (tdc1 - est)<120:  "  << forward << std::endl;
    //std::cout << "TofRawDataProvider    bool  (est - tdc2)< 24  &&   (tdc2 - est)<120:  "  << backward << std::endl;
    if( !forward && !backward ) continue;
    
    
      }
 
      m_tofDataMapTof.insert( make_pair( (*iter).first, (*iter).second ) );
    }
  }
 
  iter = m_tofDataMapTof.begin();
  while( iter != m_tofDataMapTof.end() ) {
    unsigned int iden_value = (*iter).first;
    unsigned int count = m_tofDataMapTof.count( iden_value );
    for( unsigned int i=0; i != count; i++, iter++ ) {
      ((*iter).second)->setTimes( count );
    }
  }
 
  return m_tofDataMapTof;
}
 
 
TofDataVector& TofRawDataProvider::tofDataVectorTof( double estime ) {
  if( m_tofDataVectorTof.size() != 0 ) {
    m_tofDataVectorTof.clear();
  }
 
  TofDataMap tofDataMapTof = TofRawDataProvider::tofDataMapTof( estime );
  IterTofDataMap iter = tofDataMapTof.begin();
  for( ; iter != tofDataMapTof.end(); iter++ ) {
    m_tofDataVectorTof.push_back( iter->second );
  }
  return m_tofDataVectorTof;
}
 
 
// EMC reconstruction, only charge are required
TofDataMap& TofRawDataProvider::tofDataMapEmc( double estime ) {
 
  if( !m_hasFilled ) {
    TofRawDataProvider::tofDataMapFull();
    m_hasFilled = true;
  }
 
  if( !m_tofDataMapEmc.empty() ) {
    m_tofDataMapEmc.clear();
  }
 
  IterTofDataMap iter = m_tofDataMap.begin();
  for( ; iter != m_tofDataMap.end(); iter++ ) {
    if( ( ((*iter).second)->quality() & 0xa ) != 0 ) {
 
      if( estime > 1.0e-6 ) {
    int t0clock = static_cast<int>(estime/24.0);
    int clock1 = ((*iter).second)->qclock1();
    int clock2 = ((*iter).second)->qclock2();
    bool  forward = ( ( ( t0clock - clock1 ) < tClockBackward ) && ( ( clock1 - t0clock ) < tClockForward ) );
    bool backward = ( ( ( t0clock - clock2 ) < tClockBackward ) && ( ( clock2 - t0clock ) < tClockForward ) );
    if( !forward && !backward ) continue;
      }
 
      m_tofDataMapEmc.insert( make_pair( (*iter).first, (*iter).second ) );
    }
  }
 
  iter = m_tofDataMapEmc.begin();
  while( iter != m_tofDataMapEmc.end() ) {
    unsigned int iden_value = (*iter).first;
    unsigned int count = m_tofDataMapEmc.count( iden_value );
    for( unsigned int i=0; i != count; i++, iter++ ) {
      ((*iter).second)->setTimes( count );
    }
  }
 
  return m_tofDataMapEmc;
}
 
 
TofDataVector& TofRawDataProvider::tofDataVectorEmc( double estime ) {
  if( m_tofDataVectorEmc.size() != 0 ) {
    m_tofDataVectorEmc.clear();
  }
 
  TofDataMap tofDataMapEmc = TofRawDataProvider::tofDataMapEmc( estime );
  IterTofDataMap iter = tofDataMapEmc.begin();
  for( ; iter != tofDataMapEmc.end(); iter++ ) {
    m_tofDataVectorEmc.push_back( iter->second );
  }
  return m_tofDataVectorEmc;
}