/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

/* $Id: AliTRDtrackletMCM.cxx 28397 2008-09-02 09:33:00Z cblume $ */

////////////////////////////////////////////////////////////////////////////
//                                                                        //
//  MCM tracklet                                                          //
//                                                                        //
//  Author: J. Klein (Jochen.Klein@cern.ch)                               //
//                                                                        //
////////////////////////////////////////////////////////////////////////////

#include "AliTRDtrackletMCM.h"
#include "AliLog.h"

ClassImp(AliTRDtrackletMCM)

AliTRDtrackletMCM::AliTRDtrackletMCM(UInt_t trackletWord) :
  AliTRDtrackletBase(),
  fGeo(0x0),
  fHCId(-1),
  fTrackletWord(trackletWord),
  fMCM(-1),
  fROB(-1),
  fQ0(0),
  fQ1(0),
  fNHits(0),
  fNHits0(0),
  fNHits1(0),
  fSlope(0.),
  fOffset(0.),
  fError(0.),
  fNClusters(0),
  fResiduals(0x0),
  fClsCharges(0x0)
{
  // constructor

    fGeo = new AliTRDgeometry();
    fLabel[0] = -1;
    fLabel[1] = -1;
    fLabel[2] = -1;
}

AliTRDtrackletMCM::AliTRDtrackletMCM(UInt_t trackletWord, Int_t hcid) :
  AliTRDtrackletBase(),
  fGeo(0x0),
  fHCId(hcid),
  fTrackletWord(trackletWord),
  fMCM(-1),
  fROB(-1),
  fQ0(0),
  fQ1(0),
  fNHits(0),
  fNHits0(0),
  fNHits1(0),
  fSlope(0.),
  fOffset(0.),
  fError(0.),
  fNClusters(0),
  fResiduals(0x0),
  fClsCharges(0x0)
{
  // constructor

    fGeo = new AliTRDgeometry();
    fLabel[0] = -1;
    fLabel[1] = -1;
    fLabel[2] = -1;
}

AliTRDtrackletMCM::AliTRDtrackletMCM(UInt_t trackletWord, Int_t hcid, Int_t rob, Int_t mcm) :
  AliTRDtrackletBase(),
  fGeo(0x0),
  fHCId(hcid),
  fTrackletWord(trackletWord),
  fMCM(mcm),
  fROB(rob),
  fQ0(0),
  fQ1(0),
  fNHits(0),
  fNHits0(0),
  fNHits1(0),
  fSlope(0.),
  fOffset(0.),
  fError(0.),
  fNClusters(0),
  fResiduals(0x0),
  fClsCharges(0x0)
{
  // constructor

    fGeo = new AliTRDgeometry();
    fLabel[0] = -1;
    fLabel[1] = -1;
    fLabel[2] = -1;
}

AliTRDtrackletMCM::AliTRDtrackletMCM(const AliTRDtrackletMCM &rhs) :
  AliTRDtrackletBase(rhs),
  fGeo(0x0),
  fHCId(rhs.fHCId),
  fTrackletWord(rhs.fTrackletWord),
  fMCM(rhs.fMCM),
  fROB(rhs.fROB),
  fQ0(rhs.fQ0),
  fQ1(rhs.fQ1),
  fNHits(rhs.fNHits),
  fNHits0(rhs.fNHits0),
  fNHits1(rhs.fNHits1),
  fSlope(rhs.fSlope),
  fOffset(rhs.fOffset),
  fError(rhs.fError),
  fNClusters(rhs.fNClusters),
  fResiduals(0x0),
  fClsCharges(0x0)
{
  // copy constructor

    fGeo = new AliTRDgeometry();
    fResiduals = new Float_t[fNClusters];
    fClsCharges = new Float_t[fNClusters];
    for (Int_t iCls = 0; iCls < fNClusters; iCls++) {
      fResiduals[iCls] = rhs.fResiduals[iCls];
      fClsCharges[iCls] = rhs.fClsCharges[iCls];
    }
    fLabel[0] = rhs.fLabel[0];
    fLabel[1] = rhs.fLabel[1];
    fLabel[2] = rhs.fLabel[2];
}

AliTRDtrackletMCM::~AliTRDtrackletMCM()
{
  // destructor

  delete [] fResiduals;
  delete [] fClsCharges;
    delete fGeo;
}

Int_t AliTRDtrackletMCM::GetYbin() const {
  // returns (signed) value of Y
  if (fTrackletWord & 0x1000) {
    return -((~(fTrackletWord-1)) & 0x1fff);
  }
  else {
    return (fTrackletWord & 0x1fff);
  }
}

Int_t AliTRDtrackletMCM::GetdY() const
{
  // returns (signed) value of the deflection length
  if (fTrackletWord & (1 << 19)) {
    return -((~((fTrackletWord >> 13) - 1)) & 0x7f);
  }
  else {
    return ((fTrackletWord >> 13) & 0x7f);
  }
}

void AliTRDtrackletMCM::SetLabel(Int_t label[])
{
  // set the labels (up to 3)

  fLabel[0] = label[0];
  fLabel[1] = label[1];
  fLabel[2] = label[2];
}

void AliTRDtrackletMCM::SetClusters(Float_t *res, Float_t *q, Int_t n)
{
  fNClusters = n;
  delete [] fResiduals;
  delete [] fClsCharges;

  fResiduals = new Float_t[fNClusters];
  fClsCharges = new Float_t[fNClusters];

  for (Int_t iCls = 0; iCls < fNClusters; iCls++) {
    fResiduals[iCls] = res[iCls];
    fClsCharges[iCls] = q[iCls];
  }
}