External instead of internal covariance matrix used

[u/mrichter/AliRoot.git] / HLT / comp / AliL3DataCompressorHelper.cxx

// @(#) $Id$

// Author: Anders Vestbo <mailto:vestbo@fi.uib.no>
//*-- Copyright &copy ALICE HLT Group

#include "AliL3StandardIncludes.h"

#include "AliL3RootTypes.h"
#include "AliL3Transform.h"

#include "AliL3DataCompressorHelper.h"

#if GCCVERSION == 3
using namespace std;
#endif

//_____________________________________________________________
//
//  AliL3DataCompression
//
// Interface class; binary <-> AliROOT handling of TPC data compression classes.
//


ClassImp(AliL3DataCompressorHelper)


Int_t AliL3DataCompressorHelper::fNumTimeBits = 12;
Int_t AliL3DataCompressorHelper::fNumPadBits = 12;
Int_t AliL3DataCompressorHelper::fNumChargeBits = 14;
Int_t AliL3DataCompressorHelper::fNumShapeBits = 14;
Float_t AliL3DataCompressorHelper::fXYResidualStep1 = 0.03;
Float_t AliL3DataCompressorHelper::fXYResidualStep2 = 0.03;
Float_t AliL3DataCompressorHelper::fXYResidualStep3 = 0.03;
Float_t AliL3DataCompressorHelper::fZResidualStep1 = 0.05;
Float_t AliL3DataCompressorHelper::fZResidualStep2 = 0.05;
Float_t AliL3DataCompressorHelper::fZResidualStep3 = 0.05;
Float_t AliL3DataCompressorHelper::fXYWidthStep = 0.005;
Float_t AliL3DataCompressorHelper::fZWidthStep = 0.005;
Int_t AliL3DataCompressorHelper::fClusterCharge = 100;


void AliL3DataCompressorHelper::SetBitNumbers(Int_t pad,Int_t time,Int_t charge,Int_t shape)
{
  fNumPadBits = pad;
  fNumTimeBits = time;
  fNumChargeBits = charge;
  fNumShapeBits = shape;
}

void AliL3DataCompressorHelper::SetTransverseResolutions(Float_t res1,Float_t res2,Float_t res3,Float_t width)
{
  fXYResidualStep1 = res1;
  fXYResidualStep2 = res2;
  fXYResidualStep3 = res3;
  fXYWidthStep = width;
}

void AliL3DataCompressorHelper::SetLongitudinalResolutions(Float_t res1,Float_t res2,Float_t res3,Float_t width)
{
  fZResidualStep1 = res1;
  fZResidualStep2 = res2;
  fZResidualStep3 = res3;
  fZWidthStep = width;
}

const Float_t AliL3DataCompressorHelper::GetXYResidualStep(Int_t row) 
{
  if(row < AliL3Transform::GetNRowLow())
    return fXYResidualStep1;
  else if(row < AliL3Transform::GetNRowLow() + AliL3Transform::GetNRowUp1())
    return fXYResidualStep2;
  else if(row < AliL3Transform::GetNRowLow() + AliL3Transform::GetNRowUp1() + AliL3Transform::GetNRowUp2())
    return fXYResidualStep3;
  else
    {
      cerr<<"AliL3DataCompressorHelper::GetXYResidualStep : Wrong row number "<<row<<endl;
      return -1;
    }
}

const Float_t AliL3DataCompressorHelper::GetZResidualStep(Int_t row) 
{
  if(row < AliL3Transform::GetNRowLow())
    return fZResidualStep1;
  else if(row < AliL3Transform::GetNRowLow() + AliL3Transform::GetNRowUp1())
    return fZResidualStep2;
  else if(row < AliL3Transform::GetNRowLow() + AliL3Transform::GetNRowUp1() + AliL3Transform::GetNRowUp2())
    return fZResidualStep3;
  else
    {
      cerr<<"AliL3DataCompressorHelper::GetXYResidualStep : Wrong row number "<<row<<endl;
      return -1;
    }
}