RecEmcFraction.cxx

Go to the documentation of this file.

//
//  Bes EMC Rec Shower Fraction
//  (a fraction of a hit)
//
//  Created by Wang.Zhe  2003, 10, 1
//
 
#include <iostream>
 
#include "EmcRecEventModel/RecEmcFraction.h"
#include "Identifier/EmcID.h"
 
//Constructors and destructors
RecEmcFraction::RecEmcFraction():RecEmcHit()
{
   fFraction=-1;
}
 
RecEmcFraction::RecEmcFraction(const RecEmcHit& aHit):RecEmcHit(aHit)
{
   fFraction=-1;
}
 
RecEmcFraction::~RecEmcFraction()
{
   //nothing to do
}
 
//Copy and assignment
//RecEmcFraction::RecEmcFraction(const RecEmcFraction& aFraction)
//:ContainedObject(aFraction),RecEmcHit(aFraction)
//{
//   //fCellId=aFraction.CellId();
//   //fEnergy=aFraction.Energy();
//   //fTime=aFraction.Time();
//   //
//   fFraction=aFraction.Fraction();
//}
//
//RecEmcFraction& RecEmcFraction::operator=(const RecEmcFraction& aFraction)
//{
//   if(this!=&aFraction)
//     {
//  fCellId=aFraction.CellId();
//  fEnergy=aFraction.Energy();
//  fTime=aFraction.Time();
//  //
//  fFraction=aFraction.Fraction();
//     }
//   return *this;
//}
 
//Other methods(read or write)
RecEmcFrac RecEmcFraction::getFraction() const
{
   return fFraction;
}
 
RecEmcFrac RecEmcFraction::Fraction(const RecEmcFrac& Fraction)
{
   fFraction=Fraction;
   return fFraction;
}
 
//Distance to another hit
double RecEmcFraction::Distance(const RecEmcID id) const
{
  unsigned int tht,phi;
  unsigned int tht2,phi2;
  unsigned int dtht,dphi;
  double dist;
 
  tht=EmcID::theta_module(fCellId);
  phi=EmcID::phi_module(fCellId);
  tht2=EmcID::theta_module(id);
  phi2=EmcID::phi_module(id);
 
  dtht=tht>tht2 ? tht-tht2 : tht2-tht;
  dphi=phi>phi2 ? phi-phi2 : phi2-phi;
  if(dphi>(EmcID::getPHI_BARREL_MAX()+1)/2)
  {
    dphi=EmcID::getPHI_BARREL_MAX()+1-dphi;
  }
  dist=sqrt(double(dtht*dtht+dphi*dphi));
 
  return dist;
}
 
//others
void RecEmcFraction::Dump() const
{
   cout<<"Fraction: ";
 
   cout<<fCellId<<", ";
   
   cout.width(12);
   cout.setf(ios::right);
   cout<<fEnergy<<", ";
   
   cout.width(12);
   cout.setf(ios::right);
   cout<<fTime<<", ";
 
   cout.width(12);
   cout.setf(ios::right);
   cout<<fFraction;
 
   cout<<endl;
}
 
ostream& operator<<(ostream & os, const RecEmcFraction& aFraction)
{
   os<<"Fraction: ";
 
   os<<aFraction.getCellId()<<", ";
 
   os.width(12);
   os.setf(ios::right);
   os<<aFraction.getEnergy()<<", ";
 
   os.width(12);
   os.setf(ios::right);
   os<<aFraction.getTime()<<", ";
 
   os.width(12);
   os.setf(ios::right);
   os<<aFraction.getFraction()<<", ";
   
   os.width(12);
   os.setf(ios::right);
   os<<aFraction.getEnergy()*aFraction.getFraction();
 
   os<<endl;
 
   return os;
}