StEEmcTower.cxx

#include "StEEmcTower.h"
#include <iostream>
ClassImp(StEEmcTower);

// ----------------------------------------------------------------------------
StEEmcTower::StEEmcTower()
{
  mIndex=-1;
  mLayer=-1;
  mSector=-1;
  mPhibin=-1;
  mSubsector=-1;
  mEtabin=-1;
}

// ----------------------------------------------------------------------------
StEEmcTower::StEEmcTower( const StEEmcTower &other )
{
  mIndex=other.mIndex;
  mLayer=other.mLayer;
  mSector=other.mSector;
  mPhibin=other.mPhibin;
  mSubsector=other.mSubsector;
  mEtabin=other.mEtabin;
  mET=other.mET;
  mNeighbors=other.mNeighbors;

  mRaw=other.mRaw;
  mAdc=other.mAdc;
  mEnergy=other.mEnergy;
  mStat=other.mStat;
  mFail=other.mFail;
  mName=other.mName;
  mstRawHit=other.mstRawHit;
}

// ----------------------------------------------------------------------------
void StEEmcTower::print()
{

    std::cout << "--------------------------- StEEmcTower::print()" << std::endl; 
    std::cout << "name   = " << mName << std::endl;
    std::cout << "index  = " << mIndex << std::endl; 
    std::cout << "energy = " << energy() << std::endl;
    std::cout << "raw    = " << raw() << std::endl;
    std::cout << "adc    = " << adc() << std::endl;

}

void StEEmcTower::printLine()
{
  std::cout << mName << " ADC=" << raw() << " " << adc() << " E=" << energy();
}

// ----------------------------------------------------------------------------
void StEEmcTower::index( Int_t i )
{

  mIndex  = i;
  mPhibin = i / 12;
  mEtabin = i % 12;
  mSector    = mPhibin / 5;
  mSubsector = mPhibin % 5;

  const Char_t *secs[]={"01","02","03","04","05","06","07","08","09","10","11","12"};
  const Char_t *subs[]={"A","B","C","D","E"};
  const Char_t *lays[]={"T","P","Q","R","S"};

  mName="";
  mName += secs[mSector];
  mName += lays[mLayer];
  mName += subs[mSubsector];
  mName += secs[mEtabin];
  
}

// ----------------------------------------------------------------------------
Bool_t StEEmcTower::isNeighbor( StEEmcTower t )
{

  if ( t.index() == index() ) return true;

  if ( TMath::Abs(t.etabin()-etabin()) > 1 ) return false;

  Int_t phi1 = t.phibin();
  Int_t phi2 = phibin();
  
  Int_t dphi = phi1 - phi2;

  if ( TMath::Abs(dphi) <= 1 ) return true;

  if ( phi1==59 ) phi1=-1;
  if ( phi2==59 ) phi2=-1;
  dphi = phi1 - phi2;

  if ( TMath::Abs(dphi) <= 1 ) return true;

  return false;

}