/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Event/RawDataCnv/RawDataCnv-00-04-35/src/Util/MdcConverter.cxx

Go to the documentation of this file.

#include "Identifier/MdcID.h"
#include "RawDataCnv/Util/MdcConverter.h"

MdcConverter* MdcConverter::s_instance = 0;

MdcConverter* MdcConverter::instance(int runMode)
{
   if ( s_instance == 0 ) {
      s_instance = new MdcConverter(runMode);
   }

   return s_instance;
}

void MdcConverter::destroy()
{
   if ( s_instance != 0 ) {
      delete s_instance;
      s_instance = 0;
   }
}

void MdcConverter::init(int runFrom, int runTo)
{
   m_runFrom = runFrom;
   m_runTo   = runTo;

   memset((void*)m_mdcTag, 0, 16384*4*4);
}

StatusCode MdcConverter::convert(const BufferHolder& src, MdcDigiCol* des)
{
   //int digiId = 0;
   uint32_t TEId, REId, TEData, TEOverflow, TETorQ;

   typedef pair<uint32_t, uint32_t> PairII;
   vector<uint32_t>  vHits;
   vector<PairII>    vmTDC;

   uint32_t* curTag;

   uint32_t nbuf = src.nBuf();

   for (uint32_t i = 0; i < nbuf; ++i)
   {
      uint32_t* buf = src(i);
      uint32_t* bufend= buf + src.bufSize(i);
      for ( ; buf < bufend; ++buf /*, ++digiId*/)
      {
         uint32_t digi = *buf;
         if (((digi>>18) & 0x3FFF) == 0) {
            //digiId--;
            continue;
         }

         m_builder.unPack(digi, REId, TEData, TEOverflow, TETorQ);

         TEId = m_builder.getTEID( REId );
         if ( TEId == 0xFFFFFFFF ) continue;

         curTag = m_mdcTag[REId];

         if ( curTag[0] == 0 ) {  //1st signal from detector-unit-REId
            curTag[1] = 0x7FFFFFFF;
            curTag[2] = 0x7FFFFFFF;
            curTag[3] = 0;

            int layer = MdcID::layer(Identifier(TEId));
            int wire = MdcID::wire(Identifier(TEId));
            //yzhang exchange pre-amplifier
            if ( m_runId >= m_runFrom && m_runId <= m_runTo ) {
               if (layer == 20 && wire >= 0 && wire <= 7){ 
                  TEId = 268456960 + wire;
                  curTag[3] |= 0x10;  //16
               }
               else if (layer == 42 && wire >= 0 && wire <= 7){ 
                  TEId = 268478464 + wire;
                  curTag[3] |= 0x10;  //16
               }
               else if (layer == 40 && wire >= 200 && wire <= 207){ 
                  TEId += 8;
                  curTag[3] |= 0x10;  //16
               }
               else if (layer == 40 && wire >= 208 && wire <= 215){ 
                  TEId -= 8;
                  curTag[3] |= 0x10;  //16
               }
            }

            curTag[0] = (TEId << 2);

            vHits.push_back( REId );
         }

         if ( TETorQ == 0 ) {  //T
            if ( (curTag[0]&1) == 0 ) {  //1st-found-T
               curTag[0] |= 1;
               curTag[1] = TEData;
               curTag[3] |= TEOverflow;  //0 or 1
            }
            else {  //multi-T
               curTag[3] |= 0xC;  //12, set the multi-T and 1st-T bits
               if ( TEData >= curTag[1] ) {
                  if ( TEOverflow ) TEData |= (1<<31);  //LastBit
                  vmTDC.push_back( make_pair(REId, TEData) );
               }
               else {
                  if ( curTag[3] & 1 ) curTag[1] |= (1<<31);
                  vmTDC.push_back( make_pair(REId, curTag[1]) );
                  curTag[1] = TEData;
                  curTag[3] &= (0xFFFFFFFF-1);
                  curTag[3] |= TEOverflow;
               }
            }
         }
         else {  //Q
            curTag[0] |= 2;
            curTag[2] = TEData;
            if ( TEOverflow ) curTag[3] |= 2;
         }
      }
   }

   MdcDigi* mdcDigi;
   if ( m_runMode == 1 ) {  //OnlineMode, discard the hits that only have T or Q
      for ( vector<PairII>::iterator it = vmTDC.begin(); it != vmTDC.end(); ++it ) {
         curTag   = m_mdcTag[it->first];
         if ( (curTag[0]&3) == 3 ) {
            uint32_t  data  = it->second;
            if(data>>31) continue;
            mdcDigi = new MdcDigi(Identifier(curTag[0]>>2), (data&0x7FFFFFFF), curTag[2]);
            mdcDigi->setOverflow( (curTag[3] & 0x16) | (data>>31) );
            des->push_back(mdcDigi);
         }
      }

      for ( vector<uint32_t>::iterator it = vHits.begin(); it != vHits.end(); ++it) {
         curTag   = m_mdcTag[*it];
         if ( (curTag[0]&3) == 3 ) {
            if( ((curTag[3]&3)>0) || (((curTag[3]&12)!=12) && ((curTag[3]&12)!=0)) ) {
               curTag[0] = 0;
               continue; 
            }
            mdcDigi = new MdcDigi(Identifier(curTag[0]>>2), curTag[1], curTag[2]);
            mdcDigi->setOverflow(curTag[3]);
            des->push_back(mdcDigi);
         }
         curTag[0] = 0;
      }
   }
   else {  //Not OnlineMode, keep all the hits found
      for ( vector<PairII>::iterator it = vmTDC.begin(); it != vmTDC.end(); ++it ) {
         curTag   = m_mdcTag[it->first];
         uint32_t  data     = it->second;
         mdcDigi = new MdcDigi(Identifier(curTag[0]>>2), (data&0x7FFFFFFF), curTag[2]);
         mdcDigi->setOverflow( (curTag[3] & 0x16) | (data>>31) );
         des->push_back(mdcDigi);
      }

      for ( vector<uint32_t>::iterator it = vHits.begin(); it != vHits.end(); ++it) {
         curTag   = m_mdcTag[*it];
         mdcDigi = new MdcDigi(Identifier(curTag[0]>>2), curTag[1], curTag[2]);
         mdcDigi->setOverflow(curTag[3]);
         des->push_back(mdcDigi);
         curTag[0] = 0;
      }
   }

   return StatusCode::SUCCESS;
}

StatusCode MdcConverter::convert(MdcDigiCol* src, WriteRawEvent*& des)
{
   return m_builder.pack(src, des);
}

MdcConverter::MdcConverter(int runMode)
   : m_runMode( runMode )
{
}

MdcConverter::~MdcConverter()
{
}

---