RecEmcTofHit.cxx

Go to the documentation of this file.

//
// Bes Tof Rec Hit
//
// Created by He Miao 2006,9,15
//
#include "EmcRecEventModel/RecEmcTofHit.h"
#include "Identifier/TofID.h"
#include <iostream>
 
// Constructors and destructors
RecEmcTofHit::RecEmcTofHit() { Clear(); }
 
RecEmcTofHit::RecEmcTofHit( const Identifier& id, const double e, const HepPoint3D pos ) {
  m_cellId = id;
  m_energy = e;
  m_pos    = pos;
}
 
RecEmcTofHit::~RecEmcTofHit() {}
 
// Copy and assignment
RecEmcTofHit::RecEmcTofHit( const RecEmcTofHit& tofHit ) {
  m_cellId = tofHit.CellId();
  m_energy = tofHit.Energy();
  m_pos    = tofHit.Position();
}
 
RecEmcTofHit& RecEmcTofHit::operator=( const RecEmcTofHit& tofHit ) {
  if ( this != &tofHit )
  {
    m_cellId = tofHit.CellId();
    m_energy = tofHit.Energy();
    m_pos    = tofHit.Position();
  }
  return *this;
}
 
// Other methods
int RecEmcTofHit::IntCellId() const {
  int layer = TofID::layer( m_cellId );
  int im    = TofID::phi_module( m_cellId );
  return im + layer * 88;
}
 
void RecEmcTofHit::CellId( const Identifier& id ) { m_cellId = id; }
 
void RecEmcTofHit::Energy( const double e ) { m_energy = e; }
 
void RecEmcTofHit::Position( const HepPoint3D pos ) { m_pos = pos; }
 
void RecEmcTofHit::Assign( const Identifier& id, const double e, const HepPoint3D pos ) {
  m_cellId = id;
  m_energy = e;
  m_pos    = pos;
}
 
void RecEmcTofHit::Clear() {
  m_cellId.clear();
  m_energy = 0;
  m_pos    = HepPoint3D( 0, 0, 0 );
}
 
double RecEmcTofHit::DistanceTo( const RecEmcTofHit& tofHit ) {
  // double dphi = m_pos.phi()-tofHit.m_pos.phi();
 
  // if(dphi>CLHEP::pi) dphi = CLHEP::twopi - dphi;
  // if(dphi<-CLHEP::pi) dphi = CLHEP::twopi + dphi;
 
  Hep3Vector pos1( m_pos );
  Hep3Vector pos2( tofHit.m_pos );
 
  return fabs( pos1.deltaPhi( pos2 ) );
 
  // cout<<"dphi="<<dphi<<"\tdeltaPhi="<<deltaPhi<<endl;
  // return fabs(dphi);
}
 
ostream& operator<<( ostream& os, const RecEmcTofHit& tofHit ) {
  os << "TofHit: ";
 
  os << tofHit.IntCellId() << ", ";
 
  os.width( 12 );
  os.setf( ios::right );
  os << tofHit.CellId() << ", ";
 
  os.width( 12 );
  os.setf( ios::right );
  os << tofHit.Energy() << ", ";
 
  os.width( 12 );
  os.setf( ios::right );
  os << tofHit.ZPosition() << endl;
 
  return os;
}