SRBlockFormatter.cc

Go to the documentation of this file.

#include <memory>

#include "EventFilter/EcalDigiToRaw/interface/SRBlockFormatter.h"


using namespace std;


SRBlockFormatter::SRBlockFormatter() {

}

SRBlockFormatter::~SRBlockFormatter() {

}


void SRBlockFormatter::StartEvent() {
  if (debug_) cout << "enter in StartEvent. header_ size is  " << header_.size() << endl;
  header_ .clear() ;
  if (debug_) cout << "after empty : header_ size is  " << header_.size() << endl;
  return;
}

void SRBlockFormatter::DigiToRaw(int dccid, int dcc_channel, int flag, FEDRawData& rawdata)
{

  if (debug_) cout << "enter in SRBlockFormatter::DigiToRaw " << endl;
  if (debug_) print(rawdata);

  int bx = *pbx_;
  int lv1 = *plv1_;

  int Nrows_SRP = 5;   // Both for Barrel and EndCap (without the header row)
  int SRid = (dccid -1) / 3 +1;


    // int Number_SRP_Flags = 68;

    int SRP_index;
    int icode = 1000 * dccid +  SRid;
    if (debug_) cout << "size of header_ map is " << header_.size() << endl;

    std::map<int, int>::const_iterator it_header = header_.find(icode);

    if ( it_header != header_.end() ) {
        SRP_index = rawdata.size() / 8 - Nrows_SRP;
        if (debug_) cout << "This SRid is already there." << endl; 
    }
    else {
        if (debug_) cout << "New SR Block added on Raw data " << endl;
        header_[icode] = 1;
        SRP_index = rawdata.size() / 8;
        rawdata.resize (rawdata.size() + 8 + 8*Nrows_SRP);  // 1 line for SRP header, 5 lines of data
        unsigned char* ppData = rawdata.data();
        ppData[8*SRP_index] = SRid & 0xFF;
        ppData[8*SRP_index+2] = bx & 0xFF; 
                ppData[8*SRP_index+3] = (bx & 0xF00)>>8;
        ppData[8*SRP_index+3] |= 0x80;
        ppData[8*SRP_index+4] = lv1 & 0xFF;
        ppData[8*SRP_index+5] = (lv1 & 0xF00)>>8;
                // ppData[8*SRP_index+6] = Number_SRP_Flags;  
        ppData[8*SRP_index+6] = 0;
                ppData[8*SRP_index+7] = 0x80;
                SRP_index ++;
        if (debug_) cout << "Added headers and empty lines : " << endl;
        if (debug_) print(rawdata);

        // -- put the B011 and B100 already, since for Endcap there can be empty
        // -- lines in the TCC and the SRP blocks
        unsigned char* Data = rawdata.data();
        for (int iline=SRP_index; iline < SRP_index+Nrows_SRP; iline++) {
         Data[8*iline + 7] |= 0x80;
         Data[8*iline + 3] |= 0x80;
        }
    }

    unsigned char* pData = rawdata.data();

    // -- Now the TCC Block :

    int nflags = pData[8*(SRP_index-1) +6] & 0x7F;
    nflags ++;
    pData[8*(SRP_index-1) + 6] = nflags & 0x7F;

    int jTT = (dcc_channel-1);
    int irow = jTT/16 ;
    int ival = jTT % 4;
        int kval = (jTT % 16) / 4;
    SRP_index += irow;

        if (debug_) cout << "Now add SC to SRBlock " << dec << dcc_channel << " irow ival " << dec << irow << " " << dec << ival << endl;
        if (debug_) cout << "new data will be added at line " << dec << SRP_index << endl;

 
    unsigned char* buff = &pData[8*SRP_index];
    Word64* pw = reinterpret_cast<Word64*>(const_cast<unsigned char*>(buff));
        int nbits = kval*16 + 3*ival;
     Word64 wflag = (Word64)((short int)flag & 0x7) << nbits;
     *pw |= wflag; 
         Word64 b1 = (Word64)((Word64)0x80 <<56);
     *pw |= b1 ;
     Word64 b2 = (Word64)((Word64)0x80 <<24);
     *pw |= b2 ;

        if (debug_) print(rawdata);

    
}