/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Emc/EmcRecEventModel/EmcRecEventModel-01-01-18/src/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;
}

---