EmcRecDigit2Hit Class Reference

#include <EmcRecDigit2Hit.h>

Public Member Functions
	EmcRecDigit2Hit ()
	~EmcRecDigit2Hit ()
void	Convert (const RecEmcDigitMap &aDigitMap, RecEmcHitMap &aHitMap)
void	Output (const RecEmcHitMap &aHitMap) const
void	OutputEndcap (const RecEmcHitMap &aHitMap, const unsigned int module_ew) const
void	SetAlgName (const string &name)
	EmcRecDigit2Hit ()
	~EmcRecDigit2Hit ()
void	Convert (const RecEmcDigitMap &aDigitMap, RecEmcHitMap &aHitMap)
void	Output (const RecEmcHitMap &aHitMap) const
void	OutputEndcap (const RecEmcHitMap &aHitMap, const unsigned int module_ew) const
void	SetAlgName (const string &name)
	EmcRecDigit2Hit ()
	~EmcRecDigit2Hit ()
void	Convert (const RecEmcDigitMap &aDigitMap, RecEmcHitMap &aHitMap)
void	Output (const RecEmcHitMap &aHitMap) const
void	OutputEndcap (const RecEmcHitMap &aHitMap, const unsigned int module_ew) const
void	SetAlgName (const string &name)

Detailed Description

Definition at line 12 of file InstallArea/x86_64-el9-gcc13-dbg/include/EmcRec/EmcRecDigit2Hit.h.

Constructor & Destructor Documentation

EmcRecDigit2Hit() [1/3]

EmcRecDigit2Hit::EmcRecDigit2Hit

(

)

Definition at line 22 of file EmcRecDigit2Hit.cxx.

{}

~EmcRecDigit2Hit() [1/3]

EmcRecDigit2Hit::~EmcRecDigit2Hit

(

)

Definition at line 24 of file EmcRecDigit2Hit.cxx.

{}

EmcRecDigit2Hit() [2/3]

EmcRecDigit2Hit::EmcRecDigit2Hit

(

)

~EmcRecDigit2Hit() [2/3]

EmcRecDigit2Hit::~EmcRecDigit2Hit

(

)

EmcRecDigit2Hit() [3/3]

EmcRecDigit2Hit::EmcRecDigit2Hit

(

)

~EmcRecDigit2Hit() [3/3]

EmcRecDigit2Hit::~EmcRecDigit2Hit

(

)

Member Function Documentation

Convert() [1/3]

void EmcRecDigit2Hit::Convert	(	const RecEmcDigitMap &	aDigitMap,
		RecEmcHitMap &	aHitMap )

Definition at line 26 of file EmcRecDigit2Hit.cxx.

                                                                                      {
  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;
    }
  }
}

Convert() [2/3]

void EmcRecDigit2Hit::Convert	(	const RecEmcDigitMap &	aDigitMap,
		RecEmcHitMap &	aHitMap )

Convert() [3/3]

void EmcRecDigit2Hit::Convert	(	const RecEmcDigitMap &	aDigitMap,
		RecEmcHitMap &	aHitMap )

Output() [1/3]

void EmcRecDigit2Hit::Output

(

const RecEmcHitMap &

aHitMap

)

const

Definition at line 117 of file EmcRecDigit2Hit.cxx.

                                                                {
  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 );
}

Output() [2/3]

void EmcRecDigit2Hit::Output

(

const RecEmcHitMap &

aHitMap

)

const

Output() [3/3]

void EmcRecDigit2Hit::Output

(

const RecEmcHitMap &

aHitMap

)

const

OutputEndcap() [1/3]

void EmcRecDigit2Hit::OutputEndcap	(	const RecEmcHitMap &	aHitMap,
		const unsigned int	module_ew ) const

Definition at line 160 of file EmcRecDigit2Hit.cxx.

                                                                          {
  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;
}

Referenced by Output().

OutputEndcap() [2/3]

void EmcRecDigit2Hit::OutputEndcap	(	const RecEmcHitMap &	aHitMap,
		const unsigned int	module_ew ) const

OutputEndcap() [3/3]

void EmcRecDigit2Hit::OutputEndcap	(	const RecEmcHitMap &	aHitMap,
		const unsigned int	module_ew ) const

SetAlgName() [1/3]

void EmcRecDigit2Hit::SetAlgName ( const string & name )

inline

Definition at line 22 of file InstallArea/x86_64-el9-gcc13-dbg/include/EmcRec/EmcRecDigit2Hit.h.

{ m_algName = name; }

SetAlgName() [2/3]

void EmcRecDigit2Hit::SetAlgName ( const string & name )

inline

Definition at line 22 of file InstallArea/x86_64-el9-gcc13-opt/include/EmcRec/EmcRecDigit2Hit.h.

{ m_algName = name; }

SetAlgName() [3/3]

void EmcRecDigit2Hit::SetAlgName ( const string & name )

inline

Definition at line 22 of file Reconstruction/EmcRec/include/EmcRec/EmcRecDigit2Hit.h.

{ m_algName = name; }

---

The documentation for this class was generated from the following files:

• InstallArea/x86_64-el9-gcc13-dbg/include/EmcRec/EmcRecDigit2Hit.h
• InstallArea/x86_64-el9-gcc13-opt/include/EmcRec/EmcRecDigit2Hit.h
• Reconstruction/EmcRec/include/EmcRec/EmcRecDigit2Hit.h
• EmcRecDigit2Hit.cxx