---

mMuiRawUnpack.cxx

Go to the documentation of this file.

#include <iomanip>
#include <iostream>

#include "PHNode.h"
#include "PHCompositeNode.h"
#include "PHDataNode.h"
#include "PHIODataNode.h"

#include "TMuiAddressTable.hh"
#include "mMuiRawUnpack.h"

//MUIOO interfaces
#include <MUIOO.h>
#include <TMutNode.h>
#include <TMuiHitMapO.h>

using namespace std;

mMuiRawUnpack::mMuiRawUnpack():
  _timer( PHTimeServer::get()->insert_new( "mMuiRawUnpack" ) )
 {}

mMuiRawUnpack::~mMuiRawUnpack(){}


PHBoolean
mMuiRawUnpack::event(PHCompositeNode* baseNode)
{
  
  _timer.get()->restart();
  
  Event* evt = 0;
  
  // Find the prdf node.
  try{
    evt = TMutNode<Event>::find_node(baseNode,"PRDF");
  } catch(std::exception& e){
    MUIOO::TRACE(e.what());
    return False;
  }
  PHBoolean out = event(evt, baseNode);
  _timer.get()->stop();
  
  return out;
}

PHBoolean
mMuiRawUnpack::event(Event* evt, PHCompositeNode* baseNode)
{
  if(evt == 0)
  {
    cout<<PHWHERE<<" Null event pointer.\n";
    return False;
  }
  
   // module runtime parameters
  _mod_par = TMutNode<mMuiRawUnpackPar>::find_node(baseNode,"mMuiRawUnpackPar");

  //Find or create the new TMuiHitO objects
  TMuiHitMapO* _muihitmap = 0;
  try{
    _muihitmap = TMutNode<TMuiHitMapO>::find_node(baseNode,"TMuiHitMapO");
  } catch(std::exception& e){
    MUIOO::TRACE(e.what());
    return False;
  }
  
  TMuiAddressTable* pAddress = TMuiAddressTable::Table();
  if (!pAddress) return False;

  int packetdata[MUIOO::kWordsPerFEM]={0};
  
  unsigned int row = 0;
  //static const unsigned long dcm_words_max = 136;  // size of DCM array
  // Assign module ID to associated packet ID, then get Module ID
  // according to current packet ID instead extract it from
  // DCM word. Warning!!!! If assigned packet ID for MuID DCMs are
  // changed, we need to modifiy assignment in par file.

  for (int dcm=0; dcm<MUIOO::kFEMsTotal; dcm++) {

    unsigned int packetid = _mod_par->get_id_base() + dcm;
    Packet* muipacket = evt->getPacket(packetid);
    int numwords = 0;  //Number of words actually returned by packet accessor
    
    //All muid packets are not necessarily present
    if(muipacket<=0) continue;
    
    TMuiReadoutID idHard;
    TMuiChannelId idSoft;

    unsigned int moduleID = _mod_par->get_module_id(dcm);
    
    short Arm      = pAddress->Arm(moduleID);
    short orient   = pAddress->Orient(moduleID);

    muipacket->fillIntArray(packetdata,
                            MUIOO::kWordsPerFEM,
                            &numwords,
                            "SPARSE");
    delete muipacket;
    
    // Sanity check of the word count ...
    if (numwords > MUIOO::kWordsPerFEM) {
      cout << PHWHERE << packetid
           << " has " << numwords
           << " words, maximum allowed is " << MUIOO::kWordsPerFEM << endl;
      continue;
    }

    for (unsigned int j=0; j<(unsigned int)numwords; j++) {
      unsigned int pattern = packetdata[j];

        short roc  = getbits(pattern,30,7);
        short word = getbits(pattern,23,4);

        if (    (roc < 0) || (roc >= TMuiReadoutID::kROCsPerFEM)
             || (word < 0) || (word >= TMuiReadoutID::kWordsPerROC) ) {
          cout << "mMuiDCMUnpack-E1  ROC id " << roc
               << " or word-in-ROC id " << word << " out of range" << endl;
        }
        else {

          for (short bit=0; bit<16; bit++) {
            if (getbits(pattern,bit,1) == 1) {
              idHard.Set(moduleID,roc,word,bit);
              idSoft = pAddress->SoftwareAddress(idHard);

              if ( (idSoft.Arm() != Arm) || (idSoft.Orient() != orient) ) {
                cout << "mMuiDCMUnpack-E3  DCM arm, orient = "
                     << Arm << "  " << orient
                     << "  Lookup arm, orient = "
                     << idSoft.Arm() << "  " << idSoft.Orient() << endl;
              }
              else {
#ifdef MUIOO_DEBUG                
                TMuiHitMapO::iterator hitIter =
#endif
                _muihitmap->insert_new(
                                       idSoft.Arm(),
                                       idSoft.Plane(),
                                       idSoft.Panel(),
                                       idSoft.Orient(),
                                       idSoft.TwoPack()
                                       );
#ifdef MUIOO_DEBUG
                hitIter->get()->print();
#endif

                row++;
              }
            }

          } // for (bit=0 ...)
        } // if roc and word in range
    } // for (j=0 ...)
  } // for (dcm=0 ...)

  return True;
}

unsigned int mMuiRawUnpack::getbits(unsigned int word, int start, int length){
  return (word >> (start+1-length)) & ~(~0 << length);
}

---