DstEmcShower.cxx

Go to the documentation of this file.

#include "DstEvent/DstEmcShower.h"
 
double DstEmcShower::dx() const
{
  double xx=99.;
  if(m_errorMatrix[0][0]>0) xx=sqrt(m_errorMatrix[0][0]);
  else if(m_errorMatrix[0][0]==0) xx=0.;
  return xx;
}
 
double DstEmcShower::dy() const
{
  double yy=99.;
  if(m_errorMatrix[1][1]>0) yy=sqrt(m_errorMatrix[1][1]);
  else if(m_errorMatrix[1][1]==0) yy=0.;
  return yy;
}
 
double DstEmcShower::dz() const
{
  double zz=99.;
  if(m_errorMatrix[2][2]>0) zz=sqrt(m_errorMatrix[2][2]);
  else if(m_errorMatrix[2][2]==0) zz=0.;
  return zz;
}
 
ostream& operator<<(ostream & os, const DstEmcShower& aShower)
{
  os<<"------------------DstEmcShower:"<<endl;
  os<<"Track id: "<<aShower.trackId()<<", number of hits: "<<aShower.numHits()
    <<", status: "<<aShower.status()<<", cell id: "<<aShower.cellId()
    <<", module: "<<aShower.module()<<", time: "<<aShower.time()<<endl;
  os<<"Energy: "<<aShower.energy()<<", de: "<<aShower.dE()
    <<", eseed: "<<aShower.eSeed()
    <<", e3x3: "<<aShower.e3x3()<<", e5x5: "<<aShower.e5x5()<<endl;
  os<<"Position: "<<aShower.position()
    <<", theta: "<<aShower.theta()<<", phi: "<<aShower.phi()
    <<"\ndx: "<<aShower.dx()<<", dy: "<<aShower.dy()<<", dz: "<<aShower.dz()
    <<", dtheta: "<<aShower.dtheta()<<", dphi: "<<aShower.dphi()<<endl;
  os<<"Second moment: "<<aShower.secondMoment()
    <<", lat moment: "<<aShower.latMoment()
    <<", a20 moment: "<<aShower.a20Moment()
    <<", a42 moment: "<<aShower.a42Moment()<<endl;
  os<<"Error matrix: \n"<<aShower.errorMatrix()<<endl;
 
  return os;
}