EmcRecDigit2Hit.cxx

Go to the documentation of this file.

//
//  Convert from Digit Map to Hit Map
//
//  Wang.Zhe 2004, 3, 29
//  He.Miao 2005
//
#include <iostream>
 
#include "EmcRec/EmcRecDigit2Hit.h"
#include "EmcRec/EmcRecParameter.h"
 
#include "EmcCalibConstSvc/IEmcCalibConstSvc.h"
#include "GaudiKernel/Bootstrap.h"
#include "GaudiKernel/ISvcLocator.h"
#include "RawDataProviderSvc/EmcRawDataProvider.h"
#include "RawDataProviderSvc/IRawDataProviderSvc.h"
 
// tianhl for mt
// #include "GaudiKernel/ThreadGaudi.h"
// tianhl for mt
// Constructors and destructors
EmcRecDigit2Hit::EmcRecDigit2Hit() {}
 
EmcRecDigit2Hit::~EmcRecDigit2Hit() {}
 
void EmcRecDigit2Hit::Convert( const RecEmcDigitMap& aDigitMap, RecEmcHitMap& aHitMap ) {
  RecEmcDigitMap::const_iterator ciDigitMap;
  RecEmcHit                      aHit;
 
  EmcRecParameter& Para = EmcRecParameter::GetInstance();
 
  StatusCode sc;
  // Get EmcCalibConstSvc.
  IEmcCalibConstSvc* emcCalibConstSvc;
  sc = Gaudi::svcLocator()->service( "EmcCalibConstSvc", emcCalibConstSvc );
  if ( sc != StatusCode::SUCCESS )
  {
    // cout << "EmcRecDigit2Hit Error: Can't get EmcCalibConstSvc." << endl;
  }
 
  // Get RawDataProviderSvc
  IRawDataProviderSvc* rawDataProviderSvc;
  // tianhl for mt
  std::string rawDataProviderSvc_name( "RawDataProviderSvc" );
  // if ( isGaudiThreaded( m_algName ) ) { rawDataProviderSvc_name += getGaudiThreadIDfromName(
  // m_algName ); } tianhl for mt
  sc = Gaudi::svcLocator()->service( rawDataProviderSvc_name, rawDataProviderSvc );
  if ( sc != StatusCode::SUCCESS )
  {
    // cout << "EmcRecDigit2Hit Error: Can't get RawDataProviderSvc." << endl;
  }
 
  for ( ciDigitMap = aDigitMap.begin(); ciDigitMap != aDigitMap.end(); ++ciDigitMap )
  {
 
    const RecEmcDigit& aDigit = ciDigitMap->second;
 
    aHit.CellId( aDigit.CellId() );
 
    double ein, ecorr, eout;
    ein = aDigit.ADC();
 
    // Times the calibration constant in offline mode
    if ( !( rawDataProviderSvc->isOnlineMode() ) && Para.DigiCalib() )
    {
 
      unsigned int partId = EmcID::barrel_ec( aDigit.CellId() );
      unsigned int theta  = EmcID::theta_module( aDigit.CellId() );
      unsigned int phi    = EmcID::phi_module( aDigit.CellId() );
 
      int index = emcCalibConstSvc->getIndex( partId, theta, phi );
 
      // liucx
 
      double adc2e = emcCalibConstSvc->getDigiCalibConst( index );
 
      if ( Para.DataMode() == 0 )
      {
        if ( Para.ElecSaturation() == 1 )
        {
          double emaxData = emcCalibConstSvc->getCrystalEmaxData( index );
          if ( emaxData > 0 )
          {
            //  cout<<adc2e<<"\t"<<emaxData/2.5<<endl;
            adc2e = emaxData / 2.5;
          }
        }
      }
 
      ein *= adc2e;
      // liucx
 
      // liucx ein *= emcCalibConstSvc->getDigiCalibConst(index);
    }
 
    // Correct electronic noise bias. Not used now.
    if ( ein < Para.ElecBias( 4 ) )
    {
      ecorr = Para.ElecBias( 0 ) * log( ein / Para.ElecBias( 1 ) ); // noise=0.5
    }
    else { ecorr = Para.ElecBias( 2 ) + Para.ElecBias( 3 ) * ein; }
    eout = ein;
    aHit.Energy( eout / GeV );
 
    aHit.Time( aDigit.TDC() );
    // For some reason, there are 0 energy hit.
    // Avoid this to happen, dangerous!
    if ( aHit.getEnergy() >= Para.ElectronicsNoiseLevel() )
    { // &&
      //(aHit.time()>=Para.TimeMin()) &&
      //(aHit.time()<=Para.TimeMax()))
      aHitMap[aHit.getCellId()] = aHit;
    }
  }
}
 
void EmcRecDigit2Hit::Output( const RecEmcHitMap& aHitMap ) const {
  unsigned int                 module, theta, phi;
  RecEmcID                     id;
  RecEmcHitMap::const_iterator ci_HitMap;
  RecEmcEnergy                 e;
  int                          ie;
 
  // cout<<"Output a hitmap:"<<endl;
 
  module = EmcID::EmcID::getBARREL();
 
  // cout<<"Barrel:"<<endl;
  // cout<<"+--------------------------------------------+"<<endl;
  for ( phi = EmcID::getPHI_BARREL_MAX(); phi != EmcID::getPHI_BARREL_MIN() - 1; --phi )
  {
    // cout<<"|";
    for ( theta = EmcID::getTHETA_BARREL_MIN(); theta <= EmcID::getTHETA_BARREL_MAX();
          ++theta )
    {
      id        = EmcID::crystal_id( module, theta, phi );
      ci_HitMap = aHitMap.find( id );
      if ( ci_HitMap != aHitMap.end() )
      {
        e  = ci_HitMap->second.getEnergy();
        ie = (int)( e * 10 );
        if ( ie > 9 ) { ie = 9; }
        // cout<<ie;
      }
      else
      {
        // cout<<" ";
      }
    }
    // cout<<"|"<<endl;
  }
  // cout<<"+--------------------------------------------+"<<endl;
 
  module = EmcID::EmcID::getENDCAP_EAST();
  OutputEndcap( aHitMap, module );
  module = EmcID::EmcID::getENDCAP_WEST();
  OutputEndcap( aHitMap, module );
}
 
void EmcRecDigit2Hit::OutputEndcap( const RecEmcHitMap& aHitMap,
                                    const unsigned int  module_ew ) const {
  unsigned int                 module, theta, phi;
  RecEmcID                     id;
  RecEmcHitMap::const_iterator ci_HitMap;
  RecEmcEnergy                 e;
  int                          ie;
 
  module = module_ew;
  if ( module == EmcID::EmcID::getENDCAP_EAST() )
  {
    // cout<<"East endcap:"<<endl;
  }
  else if ( module == EmcID::EmcID::getENDCAP_WEST() )
  {
    // cout<<"West endcap:"<<endl;
  }
 
  // cout<<"+--------------------------------------------------------+"<<endl;
  for ( theta = EmcID::getTHETA_ENDCAP_MAX(); theta != EmcID::getTHETA_ENDCAP_MIN() - 1;
        theta-- )
  {
    // cout<<"|";
    for ( phi = ( EmcID::getPHI_ENDCAP_MAX( theta ) - 1 ) / 2;
          phi != EmcID::getPHI_ENDCAP_MIN() - 1; phi-- )
    {
      id        = EmcID::crystal_id( module, theta, phi );
      ci_HitMap = aHitMap.find( id );
      if ( ci_HitMap != aHitMap.end() )
      {
        e  = ci_HitMap->second.getEnergy();
        ie = (int)( e * 10 );
        if ( ie > 9 ) { ie = 9; }
        if ( ( theta < 4 ) && ( phi % ( ( EmcID::getPHI_ENDCAP_MAX( theta ) + 1 ) / 16 ) ==
                                ( EmcID::getPHI_ENDCAP_MAX( theta ) + 1 ) / 16 - 1 ) )
        {
          // cout<<"*";
        }
        // cout<<ie;
        if ( ( theta == 0 || theta == 1 ) &&
             ( phi % ( ( EmcID::getPHI_ENDCAP_MAX( theta ) + 1 ) / 16 ) == 0 ) )
        {
          // cout<<"*";
        }
      }
      else
      {
        if ( ( theta < 4 ) && ( phi % ( ( EmcID::getPHI_ENDCAP_MAX( theta ) + 1 ) / 16 ) ==
                                ( EmcID::getPHI_ENDCAP_MAX( theta ) + 1 ) / 16 - 1 ) )
        {
          // cout<<" ";
        }
        // cout<<" ";
        if ( ( theta == 0 || theta == 1 ) &&
             ( phi % ( ( EmcID::getPHI_ENDCAP_MAX( theta ) + 1 ) / 16 ) == 0 ) )
        {
          // cout<<" ";
        }
      }
      if ( phi % ( ( EmcID::getPHI_ENDCAP_MAX( theta ) + 1 ) / 16 ) == 0 )
      {
        // cout<<"|";
      }
      if ( phi == ( EmcID::getPHI_ENDCAP_MAX( theta ) + 1 ) / 4 )
      {
        // cout<<"|";
      }
    }
    // cout<<endl;
  }
  // cout<<"+--------------------------------------------------------+"<<endl;
 
  for ( theta = EmcID::getTHETA_ENDCAP_MIN(); theta != EmcID::getTHETA_ENDCAP_MAX() + 1;
        theta++ )
  {
    // cout<<"|";
    for ( phi = ( EmcID::getPHI_ENDCAP_MAX( theta ) + 1 ) / 2;
          phi != EmcID::getPHI_ENDCAP_MAX( theta ) + 1; phi++ )
    {
      id        = EmcID::crystal_id( module, theta, phi );
      ci_HitMap = aHitMap.find( id );
      if ( ci_HitMap != aHitMap.end() )
      {
        e  = ci_HitMap->second.getEnergy();
        ie = (int)( e * 10 );
        if ( ie > 9 ) { ie = 9; }
        if ( ( theta < 4 ) &&
             ( phi % ( ( EmcID::getPHI_ENDCAP_MAX( theta ) + 1 ) / 16 ) == 0 ) )
        {
          // cout<<"*";
        }
        // cout<<ie;
        if ( ( theta == 0 || theta == 1 ) &&
             ( ( phi + 1 ) % ( ( EmcID::getPHI_ENDCAP_MAX( theta ) + 1 ) / 16 ) == 0 ) )
        {
          // cout<<"*";
        }
      }
      else
      {
        if ( ( theta < 4 ) &&
             ( phi % ( ( EmcID::getPHI_ENDCAP_MAX( theta ) + 1 ) / 16 ) == 0 ) )
        {
          // cout<<" ";
        }
        // cout<<" ";
        if ( ( theta == 0 || theta == 1 ) &&
             ( ( phi + 1 ) % ( ( EmcID::getPHI_ENDCAP_MAX( theta ) + 1 ) / 16 ) == 0 ) )
        {
          // cout<<" ";
        }
      }
      if ( ( phi + 1 ) % ( ( EmcID::getPHI_ENDCAP_MAX( theta ) + 1 ) / 16 ) == 0 )
      {
        // cout<<"|";
      }
      if ( phi == 3 * ( EmcID::getPHI_ENDCAP_MAX( theta ) + 1 ) / 4 - 1 )
      {
        // cout<<"|";
      }
    }
    // cout<<endl;
  }
  // cout<<"+--------------------------------------------------------+"<<endl;
}