StEEmcPoint.cxx

#include "StEEmcPoint.h"
#include "StEEmcUtil/EEmcGeom/EEmcGeomDefs.h"
#include <iostream>

#include "StEvent/StEmcPoint.h"

ClassImp(StEEmcPoint);

// ----------------------------------------------------------------------------
StEEmcPoint::StEEmcPoint()
    : TObject()
{
  mEmcPoint=0;
  mRelatives=999;
  for ( Int_t i=0;i<4;i++ ) mEnergy[i]=0.; 
  mResidueU=0.;
  mResidueV=0.;
  StEEmcClusterVec_t temp;
  for ( Int_t i=0;i<4;i++ ) mTowerClusters.push_back(temp);
  mKey=-1; // default invalid key
}

// ----------------------------------------------------------------------------
StEEmcPoint::StEEmcPoint( const StEEmcPoint &p )
    : TObject(p)
{

  mPosition       = p.mPosition;
  for ( Int_t i = 0; i < 4; i++ ) 
  mEnergy[i]      = p.mEnergy[i]; 
  mFraction       = p.mFraction;
  mTowers         = p.mTowers;
  mWeights        = p.mWeights;
  mSmdClusters[0] = p.mSmdClusters[0];
  mSmdClusters[1] = p.mSmdClusters[1];
  mEmcPoint       = p.mEmcPoint;
  mRelatives      = p.mRelatives;
  mTowerClusters  = p.mTowerClusters;

  mSector = p.mSector;
  mSigma  = p.mSigma;
  mU      = p.mU;
  mV      = p.mV;

  mResidueU = p.mResidueU;
  mResidueV = p.mResidueV;

  mKey=p.mKey;

}

// ----------------------------------------------------------------------------
StEmcPoint *StEEmcPoint::stemc()
{

  StThreeVectorF position( mPosition.X(), mPosition.Y(), mPosition.Z() );

  mEmcPoint = new StEmcPoint();
  mEmcPoint->setEnergy( energy() );
  mEmcPoint->setPosition( position );

  StThreeVectorF error( 1.0, 1.0, 1.0 );
  mEmcPoint->setPositionError( error );
  
  mEmcPoint->addCluster( kEndcapEmcTowerId, 0 );

  mEmcPoint->addCluster( kEndcapSmdUStripId, mSmdClusters[0].stemc() );
  mEmcPoint->addCluster( kEndcapSmdVStripId, mSmdClusters[1].stemc() );

  return mEmcPoint;

}

// --------------------------------------------------------------------------
Bool_t StEEmcPoint::chiSquare( const StEEmcPoint &other ) const 
{


  Float_t myChi2 = ( mSmdClusters[0].energy() - mSmdClusters[1].energy() );
  myChi2 *= myChi2;

  Float_t otherChi2 = ( other.mSmdClusters[0].energy() - other.mSmdClusters[1].energy() );
  otherChi2 *= otherChi2;

  return myChi2 < otherChi2;

}

// ----------------------------------------------------------------------------
void StEEmcPoint::print() const
{

  std::cout << "---------------------------------" << std::endl;
  std::cout << " X=" << mPosition.X() 
            << " Y=" << mPosition.Y() 
            << " energy: T=" << mEnergy[0] << " P=" << mEnergy[1] << " Q=" << mEnergy[2] << " R=" << mEnergy[3]
            << " frac=" << mFraction
            << std::endl;

  for ( UInt_t i=0; i<mTowers.size(); i++ ) 
    mTowers[i].print();

  for (UInt_t l = 0;l <= 3;l++) {
    const StEEmcClusterVec_t &v = clusters(l);
    std::cout << "layer " << l << ": # clusters " << v.size() << std::endl;
    for (UInt_t i = 0;i < v.size();i++) {
        const StEEmcCluster &c = v[i];
        std::cout << "layer " << l << ", cluster " << i << std::endl;
        c.print();
    }
  }

  std::cout << "ucluster:" << std::endl;
  mSmdClusters[0].print();
  std::cout << "vcluster:" << std::endl;
  mSmdClusters[1].print();

}