/**************************************************************************
 * 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$ */

////////////////////////////////////////////////////////////////////////////
//                                                                        //
//  Pre-Trigger Control-Box A or C for simulation                         //
//                                                                        //
//  Authors: F. Reidt (Felix.Reidt@cern.ch)                               //
//                                                                        //
////////////////////////////////////////////////////////////////////////////

#include <stdio.h>

#include "AliRun.h"
#include "AliRunLoader.h"

#include "AliLog.h"

#include "AliTRDptrgParam.h"
#include "AliTRDptrgFEB.h"
#include "AliTRDptrgLUT.h"

#include "AliTRDptrgCBAC.h"

ClassImp(AliTRDptrgCBAC)

//______________________________________________________________________________
AliTRDptrgCBAC::AliTRDptrgCBAC(AliRunLoader *rl) 
  : TObject(),
  fRunLoader(rl),
  fLUTArray(),
  fFEBArray(),
  fPosition(AliTRDptrgParam::kUnknown),
  fOperatingMode(AliTRDptrgParam::kDigits),
  fParam(0x0)
{
  //
  // ctor
  //

  AliError("default ctor - usage not recommended");

}

//______________________________________________________________________________
AliTRDptrgCBAC::AliTRDptrgCBAC(AliRunLoader *rl, 
                               AliTRDptrgParam::AliTRDptrgFEBPosition_t position,
                               AliTRDptrgParam::AliTRDptrgOperatingMode_t operatingMode,
                               AliTRDptrgParam *param) 
  : TObject(),
  fRunLoader(rl),
  fLUTArray(),
  fFEBArray(),
  fPosition(position),
  fOperatingMode(operatingMode),
  fParam(param)
{
  //
  // ctor  
  //

  this->LoadParams(); // load parameters
 
  // T0
  AliTRDptrgFEB *FEB = new AliTRDptrgFEB(this->fRunLoader, AliTRDptrgParam::kTZERO, 
                                         this->fOperatingMode, this->fPosition,
                                         0, this->fParam);
  this->fFEBArray.AddLast(FEB);

  // V0-1
  FEB = new AliTRDptrgFEB(this->fRunLoader, AliTRDptrgParam::kVZERO, this->fOperatingMode, 
                          this->fPosition, 1, this->fParam);
  this->fFEBArray.AddLast(FEB);

  // V0-2
  FEB = new AliTRDptrgFEB(this->fRunLoader, AliTRDptrgParam::kVZERO, this->fOperatingMode, 
                          this->fPosition, 2, this->fParam);
  this->fFEBArray.AddLast(FEB);

  // V0-3
  FEB = new AliTRDptrgFEB(this->fRunLoader, AliTRDptrgParam::kVZERO, this->fOperatingMode, 
                          this->fPosition, 3, this->fParam);
  this->fFEBArray.AddLast(FEB);

  // V0-4
  FEB = new AliTRDptrgFEB(this->fRunLoader, AliTRDptrgParam::kVZERO, this->fOperatingMode, 
                          this->fPosition, 4, this->fParam);
  this->fFEBArray.AddLast(FEB);

}

//______________________________________________________________________________
Bool_t AliTRDptrgCBAC::LoadParams() 
{
  //
  // Load configuration parameters
  //

  if (this->fParam != 0x0) {
    // read AliTRDptrgParam

    // get LUTs
    AliTRDptrgLUT* lut = 0x0;
    // 0
    lut = new AliTRDptrgLUT();
    lut->InitTable(10, 10, this->fParam->GetCBLUT(this->fPosition, 0), kFALSE); 
    // do not copy table data 
    this->fLUTArray.AddLast(lut);
    // 1
    lut = new AliTRDptrgLUT();
    lut->InitTable(10, 10, this->fParam->GetCBLUT(this->fPosition, 1), kFALSE); 
    // do not copy table data 
    this->fLUTArray.AddLast(lut);
  }
  else {
    // load default parameters 
    AliTRDptrgLUT* lut = new AliTRDptrgLUT();
    this->fLUTArray.AddLast(lut);
    lut = new AliTRDptrgLUT();
    this->fLUTArray.AddLast(lut);
    // the following lines are only needed for test reasons
    lut = dynamic_cast<AliTRDptrgLUT*>(this->fLUTArray.At(0));
    if (lut) {
      Int_t* initData = new Int_t[1024]; // 2^10
      for (Int_t i = 0; i < 1024; i++ ) {
        initData[i] = i;
      }
      lut->InitTable(10, 10, initData, kTRUE); // copy initData
      lut = dynamic_cast<AliTRDptrgLUT*>(this->fLUTArray.At(1));
      if (lut) {
        for (Int_t i = 1023; i >= 0; i--) {
          initData[31 - i] = i;  // inverse ramp
        }
        lut->InitTable(10, 10, initData, kTRUE); // copy initData
      }
    }
  }  

  return false;

}

//______________________________________________________________________________
AliTRDptrgCBAC::~AliTRDptrgCBAC() 
{
  //
  // Destructor
  //

  this->fLUTArray.Delete();
  this->fFEBArray.Delete();

}

//______________________________________________________________________________
Int_t* AliTRDptrgCBAC::Simulate()
{ 
  //
  // Simulate the CBAC behavior of event
  //

  Int_t nFEBs = this->fFEBArray.GetEntries();
  Int_t nLUTs = this->fLUTArray.GetEntries();

  Int_t inputVector = 0x0;

  Int_t** partResults = 0x0;  
  partResults = new Int_t* [nFEBs];

  for (Int_t iFEB = 0; iFEB < nFEBs; iFEB++) {
    AliTRDptrgFEB *feb = 0x0;
    if ((feb = dynamic_cast<AliTRDptrgFEB*>(this->fFEBArray.At(iFEB)))) {
      partResults[iFEB] = feb->Simulate();
    }
  }
  
  // combine partResults and create inputVector  
  Int_t iBit = 0;
  Int_t mask = 0x1;
  for (Int_t iFEB = 0; iFEB < nFEBs ; iFEB++) {
    for (Int_t j = 1; j <= partResults[iFEB][0]; j++) {
      AliTRDptrgLUT *lut0 = dynamic_cast<AliTRDptrgLUT*>(this->fLUTArray[0]);
      AliTRDptrgLUT *lut1 = dynamic_cast<AliTRDptrgLUT*>(this->fLUTArray[1]);
      if (lut0 && lut1) {
        if ((iBit > lut0->GetInputWidth()) 
         || (iBit > lut1->GetInputWidth())) {
          AliError("FEB result width does not match CB-A/C input with!");
	}
      }
      iBit++;
      if (partResults[iFEB][j] > 0) {
        inputVector |= mask; // Add bit to the corresponding inputVector
        mask <<= 1;
      } 
    }
  }

  AliDebug(5, Form("Inputvector: 0x%x", inputVector));
    
  // perform look up
  Int_t* result = new Int_t[nLUTs + 1]; // generate new return array
  result[0] = nLUTs; // storage array length in the first array value
  for (Int_t iLUT = 0; iLUT < nLUTs; iLUT++) { 
    // process the return value for each LUT and store the result in the array
    AliTRDptrgLUT *lutTmp = 0x0;
    if ((lutTmp = dynamic_cast<AliTRDptrgLUT*>(this->fLUTArray[iLUT]))) {
      result[iLUT + 1] = lutTmp->LookUp(inputVector);
    }
    AliDebug(4, Form("CBAC result[%d] = 0x%x",(iLUT + 1),result[iLUT + 1])); 
  }

  // Clean up
  delete [] partResults;

  return result;

}