///////////////////////////////////////////////////////////////////////////////
// Apr 15 2001 P.Murat: data part of the muon stub (should fit everything -
// CMU/CMP/CMX)
///////////////////////////////////////////////////////////////////////////////
#include "TMath.h"
#include <Stntuple/obj/TBaseMuonStub.hh>


ClassImp(TBaseMuonStub)

//_____________________________________________________________________________
TBaseMuonStub::TBaseMuonStub() {
  memset(&fGeomCode,0,((char*)&fEOR)-((char*)&fGeomCode));
}

//_____________________________________________________________________________
TBaseMuonStub::TBaseMuonStub(const TBaseMuonStub& stub) {
  // trying to make the copy constructor as efficient as possible

  memcpy(this,&stub,((char*)&fEOR)-((char*)this));
}

//_____________________________________________________________________________
TBaseMuonStub::TBaseMuonStub(const TBaseMuonStub& stub, Int_t define_all) {
  // trying to make the copy constructor as efficient as possible

  memcpy(this,&stub,((char*)&fEOR)-((char*)this));

  if (define_all) {
					// calculate missing parameters, such
					// as theta and phi angles of the stub
					// phi with respect to the wedge
    Double_t dphi = atan(fDyDx);
					// wedge phi0 refers to the angle of 
					// X axis

    //    fPhi = dphi+fChamber->Phi0()+TMath::Pi()/2;
  }
}

//_____________________________________________________________________________
TBaseMuonStub::~TBaseMuonStub() {
}

//_____________________________________________________________________________
void TBaseMuonStub::Print(Option_t* option) const {
  printf(" *** dummy TBaseMuonStub::Print called\n");
}

//_____________________________________________________________________________
void TBaseMuonStub::Clear(Option_t* option) {
  for (int i=0; i<fMaxLayer; i++) {
    fHit[i] = NULL;
  }
  fNHits    = 0;
  fChi2XY   = 1.e6;
  fChi2XZ   = 1.e6;
  fMaxLayer = 0;
}


//_____________________________________________________________________________
void TBaseMuonStub::Copy(TObject& o) const {
  // copy contents of this into 'o'

  TBaseMuonStub& stub = (TBaseMuonStub&) o;

					// we're using fUniqueID to store
					// detector code
  stub.SetUniqueID(GetUniqueID());

  stub.fFitType   = fFitType;
  stub.fNHits     = fNHits;
  stub.fTrigTower = fTrigTower;
  stub.fTheta     = fTheta;
  stub.fChi2XY    = fChi2XY;
  stub.fChi2XZ    = fChi2XZ;
  stub.fY0        = fY0;
  stub.fDyDx      = fDyDx;
  stub.fZ0        = fZ0;
  stub.fDzDx      = fDzDx;
  stub.fChamber   = fChamber;
  for (int i=0; i<8; i++) {
    stub.fDriftSign[i] = fDriftSign[i];
    stub.fHit[i]       = fHit[i];
  }
  stub.fMaxLayer = fMaxLayer;
}



//____________________________________________________________________________
void TBaseMuonStub::Streamer(TBuffer &R__b)
{

  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 >> fGeomCode;
    R__b >> fDigiCode;
    R__b >> fFitType;
    R__b >> fNHits;
    if(R__v>=2) {
      R__b >> fNcHits;
    } else {
      fNcHits = 0;
    }
    R__b >> fTrigTower;
    R__b >> fPhi;
    R__b >> fTheta;
    R__b >> fChi2XY;
    R__b >> fChi2XZ;
    R__b >> fY0;
    R__b >> fDyDx;
    R__b >> fZ0;
    R__b >> fDzDx;
    R__b.CheckByteCount(R__s, R__c, TBaseMuonStub::IsA());
  } else {
    R__c = R__b.WriteVersion(TBaseMuonStub::IsA(), kTRUE);
    TObject::Streamer(R__b);
    R__b << fGeomCode;
    R__b << fDigiCode;
    R__b << fFitType;
    R__b << fNHits;
    R__b << fNcHits;
    R__b << fTrigTower;
    R__b << fPhi;
    R__b << fTheta;
    R__b << fChi2XY;
    R__b << fChi2XZ;
    R__b << fY0;
    R__b << fDyDx;
    R__b << fZ0;
    R__b << fDzDx;
    R__b.SetByteCount(R__c, kTRUE);
  }
}