HarpoDccEvent.cxx

Go to the documentation of this file.

// 
// File HarpoDccEvent.cxx
//
#include "HarpoDccEvent.h"
#include "HarpoDccPacket.h"
#include "HarpoDebug.h"
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <cassert>
// for swab 
#include <unistd.h> 

ClassImp(HarpoDccEvent)

HarpoDccEvent::HarpoDccEvent()
{
  fDet= XDCC;
  memset(&fEventHeader,0,sizeof(fEventHeader));
  fEventMap = (HarpoDccMap *) NULL;
  fEventData = (UShort_t *)NULL;
  decoder = new HarpoDccDecode();

  fZeroSuppr = false;
  fTimeStamp = 0;
  //  std::cout << "Dcc Event Created" << std::endl;
}

 HarpoDccEvent::~HarpoDccEvent()
{
  if(gHarpoDebug>1)
    std::cout << fEventMap << "   " << decoder << std::endl;
  delete [] fEventData;
  delete fEventMap;
  //  delete decoder;
  decoder->Delete();
  if(gHarpoDebug>1)
    std::cout << "Dcc Event Deleted" << std::endl;
}

void          HarpoDccEvent::SetHeader(UInt_t len, UInt_t evnum)
{
  fEventHeader.eventSize=len;
  fEventHeader.eventNumb=evnum;
  if(gHarpoDebug>1)
    std::cout << "Dcc Event SetHeader " 
              << fEventHeader.eventSize << " "
              << fEventHeader.eventNumb << std::endl;
  }

void          HarpoDccEvent::SetHeader(UShort_t *data)
{
  memcpy(&fEventHeader,data,sizeof(EventHeader_t));
}

void          HarpoDccEvent::SetData(UShort_t *data, size_t dlen)
{
  //if(fEventData!=NULL) delete [] fEventData;
  fEventData = new UShort_t[dlen/sizeof(UShort_t)];

  swab(data,fEventData,dlen);  // copy and  ntohs ( byte swapping)
  //  memcpy(&fEventData,data,dlen);

  //  std::cout << std::dec << fEventData[0] << " " << dlen <<  std::endl;
  //std::cout << std::hex << fEventData[1] << std::endl;
  //std::cout << std::hex << fEventData[2] << std::dec << std::endl;
  // Decode Frames
  HarpoDccPacket frm(fEventData,dlen);
  HarpoDccFrame tmpFrame;
  while (frm.Decode(tmpFrame) ) 
    {
      fFrames.push_back(tmpFrame);
    }
}

// void HarpoDccEvent::ResetData()
// {
//   fFrames.clear();
// }


void HarpoDccEvent::print() const
{
  std::cout << "---DccEvent Plane " << fDet  << std::endl;
  std::cout << "Header num : " << fEventHeader.eventNumb
            << " : size : " << fEventHeader.eventSize << std::endl;
  // std::cout << " Num Frames : " << fFrames.size() << std::endl;
  //  std::cout << " FRame : " << fFrames[0].compress << std::endl;
  //fFrames[0].print();
}

// Defined in DetEvent
// EventHeader_t *HarpoDccEvent::GetHeader()
// {
//  return &fEventHeader;
// }

UShort_t    *HarpoDccEvent::GetData()
{
 return fEventData;
}

inline Int_t isCtlWord(UShort_t w)
{
  w >>= 12;
  if (w == 3 ) w=2;
  return w;
}

inline UShort_t iAdcData(UShort_t w) { return (w & 0xfff); }

inline Int_t iBinIndex(UShort_t w) { return (w & 0x1ff); }

inline void iChannelWord(UShort_t w, UShort_t &a1, UShort_t &a2)
{
  w &= 0xDFFF;
  a2 = ( (w & 0x1e00) >> 9 );
  a1 = (w & 0x1FF);
  //  std::cout << "ArgRb " << std::hex << w << std::dec << " " << w << std::endl;
}

// and skip other valuse which coreespond to special cols 511
// all noise channels index ( 288 ..304) after decoding
inline void fillIf(HarpoDccMap *m, Int_t row, Int_t col, UShort_t word)
{
  //  if ( ( row >= 0 ) && (row < m->map.GetNrows() ) 
  //       && (col >= 0 ) && (col < m->map.GetNcols()))
  //    m->map(row,col) = iAdcData(word) ;
  m->SetRawData(row,col,iAdcData(word));
}

Bool_t        HarpoDccEvent::UnPackData()
{
  Bool_t res = true;
  Int_t rchannel, channel, i;
  
  if(fEventMap!=NULL) return kTRUE;

  if(gHarpoDebug>1) 
    Info("UnPackData","fEventMap does not exist. Unpacking RAW data");

  /* if(fEventMap==NULL) */ fEventMap = new HarpoDccMap();
  //  fEventMap->ResetMaps();

   std::list<HarpoDccFrame>::const_iterator frm;
  
   // std::cout << "<Dcc Map> raws "<< fEventMap->map.GetNrows()
   //        << " cols " << fEventMap->map.GetNcols() << std::endl;

    for(frm = fFrames.begin();frm != fFrames.end(); frm++)
     {
       if (frm->pkgtype == 0 ) {
         if ( frm->compress == 0 )
           {
             //dbg
             // frm->print();
             // for (i=0;i<12;i++) {
             //   std::cout<< "dccData " << i << " " << std::hex 
             //              << frm->dccData[i] << std::dec << std::endl;
             // }
             //dbg
             if ( frm->size == 1044 ) { // All 512 words
             rchannel = frm->ChannelNum() ; // channel from header
             // std::cout << "0 rChannel " << rchannel
             //                << " size " << frm->size
             //                << " type : " << frm->msgtype << std::endl;
             //
             assert(rchannel >= 0 && rchannel < NALL); 
             channel = decoder->decode(rchannel,fDet);
             //
             for (i=0 ; i<NADC;i++)// Data[512] are special
               {
                 fillIf(fEventMap,channel,i,frm->dccData[i]);
                 //              fEventMap->map(channel,i) = iAdcData(frm->dccData[i]);
               }
             } else {
               std::cout << " No Zero Supperssion frame size " << 
                 frm->size << " != 1044" << std::endl;
               res = false;
             }
           }
         else // Compressed data
           {
             rchannel = frm->ChannelNum();
             assert(rchannel >= 0 && rchannel < NALL); 
             // std::cout << "z rChannel " << rchannel 
             //                << " size " << frm->size
             //                << " type : " << frm->msgtype << std::endl;
             //
             channel = decoder->decode(rchannel,fDet);
             //channel = rchannel;
             //
             //dbg
             //frm->print();
             // for (i=0;i<12;i++) {
             //   std::cout<< "dccData " << i << " " << std::hex 
             //              << frm->dccData[i] << std::dec << std::endl;
             // }
             //dbg
             UShort_t *samp;
             UShort_t iival;
             Int_t iibin = -1;
             Int_t iitype = -1;
             Bool_t flag512 = false;
             int jj;
             for (samp = frm->dccData, jj= 0; samp != frm->dccTrailer; samp++,jj++)
             {
               iival = *samp;
               // std::cout << "Sampl(" << jj << ")="
               //                << std::hex <<  iival << std::dec << std::endl;
               iitype = isCtlWord(iival);
               if ( iitype == 0 ) { // Data Word
                 if( flag512 ) {
                   // std::cout << "Flag 512 data " << iibin
                   //        << " = " << iival << std::endl;
                 } else {
                   if (iibin >= 0 && iibin < (511 + 10)) {
                     if ( iAdcData(iival) > 0) { // skeep padding and firs zeros
                       fillIf(fEventMap,channel,iibin-10,iival);
                       // fEventMap->map(channel,iibin-10) = iAdcData(iival);
                     }
                   } else {
                     // std::cout << "iibin index range " << iibin << std::endl;
                   }
                 }
                 iibin++;
               } else if ( iitype == 1 ) { // New cluster index
                 iibin = iBinIndex(iival);
                 //std::cout << "iibin  " << iibin << std::endl;
                 if ( iibin == 511 ) {
                   flag512 = true;
                 } else {
                   flag512 = false;
                 }
               } else if ( iitype == 2 ) { // Found arg word
                 UShort_t a1;
                 UShort_t a2;
                 iChannelWord(iival,a1,a2);
                 // std::cout << std::hex << "Addr "
                 //        << iival << " " << a1 << " " << a2 
                 //        << std::dec << std::endl;
                 //
                 rchannel = frm->ChannelNum(a1,a2);
                 channel = decoder->decode(rchannel,fDet);
                 //channel = rchannel;
                 flag512 = false;
               } else if ( iitype == 4 ) { // channel word counter
                 // std::cout << "wcnt  " << iival << std::endl;
                 // just skeep it
               } else {  // format error
               std::cout << "Zero Supperssion frame bad data "
                         << std::hex << iival << std::dec << std::endl;
               res = false;
               break ;
               }
             }
           }  // zero sup
       }
       else
         {
           // std::cout << "Warning unsupported frame" << std::endl;
           // frm->print();
         }
     }
    //    fEventMap->rmap.Draw();
    //  std::cout << "Unpack Data" << std::endl;

    fZeroSuppr = (fFrames.front().compress != 0);
    fTimeStamp = fFrames.front().timestamp;

  return res;
}

HarpoDccMap   *HarpoDccEvent::GetMap()
{
  if(gHarpoDebug>1) std::cout << "HarpoDccEvent::GetMap" << std::endl;

  if (UnPackData()) {
    if(gHarpoDebug>1) 
      Info("GetMap","Returning %p",(void *) fEventMap);
    return fEventMap;
  }
  else {
    return NULL;
  }
}

Bool_t      HarpoDccEvent::IsZeroSuppressed(Int_t /* after */ )
{
  return fZeroSuppr;
  //  return (fFrames.front().compress != 0);
  //  return (fFrames[after].compress != 0);
}

ULong_t      HarpoDccEvent::GetTimeStamp(Int_t /* after */)
{
  if(gHarpoDebug>2) Info("GetTimeStamp","%d",fFrames.front().timestamp);
  return fTimeStamp;
  // return (ULong_t) fFrames.front().timestamp;
}