#pragma implementation
#include "hshowerhit.h"
#include <iostream>
ClassImp(HShowerHit)

//_HADES_CLASS_DESCRIPTION 
//////////////////////////////////////////////////////////////////////
//HShowerHit
//
// Structure for Hit Level
// One object describes one hit in detector
//
////////// **********ShowerHit defintion ****************/////////////77

Int_t HShowerHit::clear()
{
   m_fCharge=0.0F;
   m_fX = 0.0;
   m_fY = 0.0;
   m_fSigmaX = 0.0;
   m_fSigmaY = 0.0;


   m_fSum0 = 0.0;  //sum of charge
   m_fSum1 = 0.0;  //sum of charge
   m_fSum2 = 0.0;  //sum of charge
   m_fSum2_25 = 0.0;  //sum of charge
   m_fSum2_N = 0.0;  //sum of charge
   m_fVar0 = 0.0;  //var. of sum of charge
   m_fVar1 = 0.0;  //var. of sum of charge
   m_fVar2 = 0.0;  //var. of sum of charge
   m_fLocalMax = 0.0; //charge in local max

   m_fSumDiff01 = 0.0;
   m_fSumDiff02 = 0.0;
   m_fSumDiff12 = 0.0;
   m_fShower = 0.0; //shower - diff m_fSum2 - m_fSum1 if > 0

   m_nSector = -1;
   m_nModule = -1;
   m_nRow = -1;
   m_nCol = -1;
   m_nAddress = -1;

   m_fLabX =0.0;  //spatial coord in lab
   m_fLabY =0.0;
   m_fLabZ =0.0;

   m_fRad  =0.0; // sphere coord in lab
   m_fPhi  =0.0;
   m_fTheta=0.0;

   m_nClusterSize0=0;
   m_nClusterSize1=0;
   m_nClusterSize2=0;
 
   return 1;
}
Float_t HShowerHit::getZ()
{
  //Z coordinate of the hit in case of low
  //efficiency mode;

  if( m_fSum0 < 0.0000001 )
    { 
      return 67.6;
    }
  else
    {
      return -4.;
    } 
}

void HShowerHit::setSum(Int_t i, Float_t fSum) {
  switch(i) {
     case 0: m_fSum0 = fSum;
             break;
 
     case 1: m_fSum1 = fSum;
             break;
 
     case 2: m_fSum2 = fSum;
             break;
 
     default: break;
  }
}

void HShowerHit::setVar(Int_t i, Float_t fVar) {
  switch(i) {
     case 0: m_fVar0 = fVar;
             break;

     case 1: m_fVar1 = fVar;
             break;

     case 2: m_fVar2 = fVar;
             break;

     default: break;
  }
}

void HShowerHit::setClusterSize(Int_t i, Int_t cs) {
  switch(i) {
     case 0: m_nClusterSize0 = cs;
             break;

     case 1: m_nClusterSize1 = cs;
             break;

     case 2: m_nClusterSize2 = cs;
             break;

     default: break;
  }

}

void HShowerHit::setSum25(Float_t fSum) {
    m_fSum2_25 = fSum;
    m_fSum2_N = fSum * 0.36;
}

Float_t HShowerHit::getSum(Int_t i) {
  Float_t fSum = -1.0;
  switch(i) {
     case 0: fSum = m_fSum0;
             break;

     case 1: fSum = m_fSum1;
             break;

     case 2: fSum = m_fSum2;
             break;

     default: break;
  }

  return fSum;
}

Float_t HShowerHit::getVar(Int_t i) {
  Float_t fVar = -1.0;
  switch(i) {
     case 0: fVar = m_fVar0;
             break;

     case 1: fVar = m_fVar1;
             break;

     case 2: fVar = m_fVar2;
             break;

     default: break;
  }

  return fVar;
}

Int_t   HShowerHit::getClusterSize(Int_t i) {
  Int_t cs = -1;
  switch(i) {
     case 0: cs = m_nClusterSize0;
             break;

     case 1: cs = m_nClusterSize1;
             break;

     case 2: cs = m_nClusterSize2;
             break;

     default: break;
  }

  return cs;
}

void HShowerHit::updateCalc() {
   m_fSumDiff01 = m_fSum1 - m_fSum0;
   m_fSumDiff02 = m_fSum2 - m_fSum0;
   m_fSumDiff12 = m_fSum2 - m_fSum1;
}

// -----------------------------------------------------------------------------

void HShowerHit::Streamer(TBuffer &R__b)
{
   // Stream an object of class HShowerHit.

   if (R__b.IsReading()) {
      Version_t R__v = R__b.ReadVersion(); 
      if (R__v == 1) 
      {
      Version_t R__v;
            R__v = R__b.ReadVersion();	// simulate HDataObject v 1 readout
      }
      TObject::Streamer(R__b);
      R__b >> m_fX;
      R__b >> m_fY;
      R__b >> m_fLabX;
      R__b >> m_fLabY;
      R__b >> m_fLabZ;
      R__b >> m_fRad;
      R__b >> m_fPhi;
      R__b >> m_fTheta;
      R__b >> m_fSum0;
      R__b >> m_fSum1;
      R__b >> m_fSum2;
      R__b >> m_fSum2_25;
      R__b >> m_fSum2_N;
      R__b >> m_fVar0;
      R__b >> m_fVar1;
      R__b >> m_fVar2;
      R__b >> m_fLocalMax;
      R__b >> m_fSumDiff01;
      R__b >> m_fSumDiff02;
      R__b >> m_fSumDiff12;
      R__b >> m_fShower;
      R__b >> m_nClusterSize0;
      R__b >> m_nClusterSize1;
      R__b >> m_nClusterSize2;
      R__b >> m_nSector;
      R__b >> m_nModule;
      R__b >> m_nRow;
      R__b >> m_nCol;
      R__b >> m_fCharge;
      R__b >> m_nAddress;

	  if(R__v <= 2)
	  {
	  Int_t iEventNr;
		  R__b >> iEventNr;
      }
      if(R__v >= 4)
      {
         R__b >>m_fSigmaX;
         R__b >>m_fSigmaY;
      } 

   } else {
      R__b.WriteVersion(HShowerHit::IsA());
      TObject::Streamer(R__b);
      R__b << m_fX;
      R__b << m_fY;
      R__b << m_fLabX;
      R__b << m_fLabY;
      R__b << m_fLabZ;
      R__b << m_fRad;
      R__b << m_fPhi;
      R__b << m_fTheta;
      R__b << m_fSum0;
      R__b << m_fSum1;
      R__b << m_fSum2;
      R__b << m_fSum2_25;
      R__b << m_fSum2_N;
      R__b << m_fVar0;
      R__b << m_fVar1;
      R__b << m_fVar2;
      R__b << m_fLocalMax;
      R__b << m_fSumDiff01;
      R__b << m_fSumDiff02;
      R__b << m_fSumDiff12;
      R__b << m_fShower;
      R__b << m_nClusterSize0;
      R__b << m_nClusterSize1;
      R__b << m_nClusterSize2;
      R__b << m_nSector;
      R__b << m_nModule;
      R__b << m_nRow;
      R__b << m_nCol;
      R__b << m_fCharge;
      R__b << m_nAddress;
      R__b << m_fSigmaX;
      R__b << m_fSigmaY;
   }
}