//*-- AUTHOR : Vladimir Pechenov
//*-- $Id: hmdcclusinf.cc,v 1.10 2009-10-08 16:48:53 pechenov Exp $
//*-- Modified : 04.04.2003 by Vladimir Pechenov
//*-- Modified : 09.07.2002 by Vladimir Pechenov

//_HADES_CLASS_DESCRIPTION 
/////////////////////////////////////////////////////////////////////////
//
//  HMdcClusInf
//
//  this class keep some information about cluster 
// * Short_t getClusIndex()
//      Returns the index of HMdcClus container
// * Short_t getClusFitIndex()
//      Returns the index of HMdcClusFit container
// * Char_t getFitStat()
//      Returns the fitStatus, 0 - hit filled by cluster, 1 - by fit
// * Short_t getClusSize()
//      Returns the number of bins in cluster (cluster size)
// * Short_t getNMergClus()
//      Returns the number of merged clusters
// * Short_t getNDrTimes()
//      Returns the number of drift times in cluster
// * Short_t getNDrTimesFitted()
//      Returns the number of drift times used by the fitter
// * Short_t getNLayers()
//      Returns the number of layers in cluster
// * Short_t getNLayersFitted()
//      Returns the number of layers used by the fitter
// * void getShape(Float_t& sigma1,Float_t& sigma2,Float_t& alpha)
// * Float_t getAlpha()
// * Float_t getSigma1()
// * Float_t getSigma2()
//      Returns the cluster shape (all in sector coor. sys.):
//      sigma1 - long axis, (mm)
//      sigma2 - short axis (mm)
//      alpha - angle (deg.) between long axis of cluster and axis X
// * Float_t getRatio()
//      Returns sigma2/sigma1
// * Char_t getLevelClFinding()
//      Returns the level of cl.finding (minimal num. of layers)
// * UChar_t getNModInCl() {return uch3&0xF;}
//      Returns the num. of modules in cluster
// * UChar_t getNModInFit()
//      Returns the number of modules in fit
//
//
//
//
//
/////////////////////////////////////////////////////////////////////////

#include "hmdcclusinf.h"
#include <cstdlib> 

ClassImp(HMdcClusInf)

void HMdcClusInf::print(void) const{
  printf("fit_status=%i clusSize=%i nMergedClus=%i nLayers=%i nLayersFitted=%i nDrTimes=%i nDrTimesFitted=%i sigma1=%f sigma2=%f alpha=%f dedx=%f dedxSigma=%f dedxNWire=%i dedxNWireCut=%i\n",
      fitStatus,clusSize,nMergClus,nLayers,nLayersFitted,nDrTimes,nDrTimesFitted,sigma1,sigma2,alpha,dedx,dedxSigma,dedxNWire,dedxNWireCut);
}

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

   UInt_t R__s, R__c;
   if (R__b.IsReading()) {
      Version_t R__v = R__b.ReadVersion(&R__s, &R__c); if (R__v) { }
      TObject::Streamer(R__b);
      R__b >> fitStatus;
      R__b >> levelClFinding;
      R__b >> nModsClus;
      R__b >> nModsFit;
      R__b >> clusSize;
      R__b >> nMergClus;
      if(R__v>2){
      R__b >> nLayers;
      R__b >> nLayersFitted;
      }
      R__b >> nDrTimes;
      R__b >> nDrTimesFitted;
      R__b >> indClus;
      R__b >> indClusFit;
      R__b >> sigma1;
      R__b >> sigma2;
      R__b >> alpha;
      if(R__v>3){
        R__b >> xHit;
        R__b >> yHit;
        R__b >> xClus;
        R__b >> yClus;
      }
      if(R__v>1){
	  R__b >> dedx;
	  R__b >> dedxSigma;
	  R__b >> dedxNWire;
	  R__b >> dedxNWireCut;
      }
      
      R__b.CheckByteCount(R__s, R__c, HMdcClusInf::IsA());
   } else {
      R__c = R__b.WriteVersion(HMdcClusInf::IsA(), kTRUE);
      TObject::Streamer(R__b);
      R__b << fitStatus;
      R__b << levelClFinding;
      R__b << nModsClus;
      R__b << nModsFit;
      R__b << clusSize;
      R__b << nMergClus;
      R__b << nLayers;
      R__b << nLayersFitted;
      R__b << nDrTimes;
      R__b << nDrTimesFitted;
      R__b << indClus;
      R__b << indClusFit;
      R__b << sigma1;
      R__b << sigma2;
      R__b << alpha;
      R__b << xHit;
      R__b << yHit;
      R__b << xClus;
      R__b << yClus;
      R__b << dedx;
      R__b << dedxSigma;
      R__b << dedxNWire;
      R__b << dedxNWireCut;
      R__b.SetByteCount(R__c, kTRUE);
   }
}