#include "evt/Event.hh"
#include <Bbos/bank_interface.hh>

#include <Stntuple/obj/TClcDataBlock.hh>

#include <Stntuple/mod/StntupleUtilities.hh>
#include <Stntuple/mod/InitStntupleDataBlocks.hh>

//_____________________________________________________________________________
Int_t StntupleInitClcDataBlock(TStnDataBlock* block, AbsEvent* event, int mode) 
{
  // initialize CLC data block with the `event' data

  int ev_number, rn_number;
					// check if block has already been 
					// initialized

  ev_number = AbsEnv::instance()->trigNumber();
  rn_number = AbsEnv::instance()->runNumber();

  if (block->Initialized(ev_number,rn_number)) return 0;
  
  TClcDataBlock* clc_block = (TClcDataBlock*) block;
  clc_block->Clear();

  TRYBOS_BANK *clad, *cltd;
  int  _nIsolated = 0;
  int  tmp[96];

  clad = TRYBOS_BANK::Find(event,"CLAD");

  if (clad) {
    for (int i=0; i<TClcDataBlock::kNChannels; i++) {
      clc_block->fAdcData[i] = clad->sw(5+2*i);
    }
//-----------------------------------------------------------------------------
// check if channel is isolated:  loop over channels and identify neighborhs
//-----------------------------------------------------------------------------
    for (int i=0; i<TClcDataBlock::kNChannels; i++) {
      int nin, ind[6];
      int key = i/16 + 1;
      switch (key) {
             case 1:
               nin = 4;
               ind[0] = (i+14)%16 + 1;
               if(i==15) ind[1]=0; else ind[1] = i%16 + 1;
               if(i==1) ind[2]=0; else ind[2] = (i+14)%16 + 17;
               ind[3] = i+16;
               break;

             case 2:
               nin = 6;
               ind[0] = i - 16;
               ind[1] = i%16 + 1;
               if(i==17) ind[2]=16; else ind[2] = (i-2)%16 + 17;
               ind[3] = i%16 + 17;
               ind[4] = i + 16;
               if(i==31) ind[5]=0; else ind[5] = i%16 + 33;
               break;

             case 3:
               nin = 4;
               ind[0] = (i-2)%16 + 17;
               ind[1] = i - 16;
               if(i==33) ind[2]=16; else ind[2] = (i-2)%16 + 33;
               if(i==47) ind[3]=32; else ind[3] = i%16 + 33;
               break;

             case 4:
               nin = 4;
               ind[0] = (i-2)%16 + 49;
               if(i==63) ind[1]=48; else ind[1] = i%16 + 49;
               if(i==49) ind[2]=48; else ind[2] = (i-2)%16 + 65;
               ind[3] = i+16;
               break;

             case 5:
               nin = 6;
               ind[0] = i - 16;
               ind[1] = i%16 + 49;
               if(i==65) ind[2]=64; else ind[2] = (i-2)%16 + 65;
               ind[3] = i%16 + 65;
               ind[4] = i + 16;
               if(i==79) ind[5]=48; else ind[5] = i%16 + 81;
               break;

             case 6:
               nin = 4;
               ind[0] = (i-2)%16 + 65;
               ind[1] = i - 16;
               if(i==81) ind[2]=64; else ind[2] = (i-2)%16 + 81;
               if(i==95) ind[3]=80; else ind[3] = i%16+ 81;
               break;
      }
	
      TClcChannel *ch = clc_block->Channel(i);	
      for ( int ich = 0; ich < nin; ich++) {
				// thersholds for isolation are hardcoded for now
				// probably could change to set it trhough talk to
				// in the future 
	if (! clc_block->fAdcData[i] > 0 || clc_block->fAdcData[ind[ich]] > 100) {
				// channel is not isolated 
				// break out from loop
	  ch->SetIsolated(0);            
	  break;                         
	}
      }
				// if it is isolated put it in temporary structure
      if ( ch->IsIsolated()) {
	tmp[_nIsolated] = i;
	_nIsolated++;
      }
    }
  }
//-----------------------------------------------------------------------------
// put it in the resizable container later will be stored in stntuple
//-----------------------------------------------------------------------------
  clc_block->fListOfIsolatedChannels.Reserve(_nIsolated);
  clc_block->fListOfIsolatedChannels.Append(tmp,_nIsolated); 

  cltd = TRYBOS_BANK::Find(event,"CLTD");
  StntupleUnpackTdcBank(cltd,clc_block->fTdcData);
  clc_block->InitTdcChannels();

  clc_block->f_RunNumber   = rn_number;
  clc_block->f_EventNumber = ev_number;

  return 0;
}