mEmcPRDFToRaw.C

Go to the documentation of this file.

#include "mEmcPRDFToRaw.h"
#include "dEmcRawDataWrapper.h"

#include "Event.h"

#include "phool.h"
#include "PHIODataNode.h"
#include "PHCompositeNode.h"
#include "PHNodeIterator.h"
#include <iostream>
using namespace std;

typedef PHDataNode<Event> EventNode_t;
typedef PHIODataNode<PHTable> TableNode_t;
typedef PHDataNode<dEmcRawDataWrapper> dEmcRawDataNode_t;

static const int FeeToSoftwareMap[144] = 
{ 121, 120, 133, 132,  97,  96, 109, 108, 73, 72, 85, 84,
   49,  48,  61,  60,  25,  24,  37,  36,  1,  0, 13, 12,
  123, 122, 135, 134,  99,  98, 111, 110, 75, 74, 87, 86,
   51,  50,  63,  62,  27,  26,  39,  38,  3,  2, 15, 14,
  125, 124, 137, 136, 101, 100, 113, 112, 77, 76, 89, 88,
   53,  52,  65,  64,  29,  28,  41,  40,  5,  4, 17, 16,
  127, 126, 139, 138, 103, 102, 115, 114, 79, 78, 91, 90,
   55,  54,  67,  66,  31,  30,  43,  42,  7,  6, 19, 18,
  129, 128, 141, 140, 105, 104, 117, 116, 81, 80, 93, 92,
   57,  56,  69,  68,  33,  32,  45,  44,  9,  8, 21, 20,
  131, 130, 143, 142, 107, 106, 119, 118, 83, 82, 95, 94,
   59,  58,  71,  70,  35,  34,  47,  46, 11, 10, 23, 22
};

int mEmcPRDFToRaw(PHCompositeNode* topNode)
{
  PHNodeIterator iter(topNode), *j;
  PHNode *n1,*n2;

  // Find the node with the RawData object
  n2 = iter.findFirst("PHCompositeNode", "EMC");
  if (!n2) {
    cout << "mEmcPRDFToRAW: error, couldn't find EMC Node" << endl;
    return 0;
  }
  PHCompositeNode *emcNode = static_cast<PHCompositeNode*>(n2);
  dEmcRawDataWrapper *demcraw;
  j = new PHNodeIterator(emcNode);

  dEmcRawDataNode_t *demcrawNode = static_cast<dEmcRawDataNode_t*>(j->findFirst("PHIODataNode","dEmcRawData"));
  if (!demcrawNode) {
    cout << "mEmcPRDFToRAW: error, couldn't find dEmcRawData" << endl;
    return 0;
  }
  demcraw = static_cast<dEmcRawDataWrapper*>(demcrawNode->getData());
  delete j;

  // Find the node with the Event object; it is called "PRDF".
  n1 = iter.findFirst("PHDataNode", "PRDF");
  if (!n1) {
    return 1;
  }

  // Extract the Event object from the node.
  EventNode_t* node = static_cast<EventNode_t*>(n1);
  Event* event = node->getData();
  if (!event) {
    return 1;
  }

  // Get the relevant packets from the Event object and transfer the
  // data to the subsystem-specific table.
  Packet* p;
  static const int id_base     = 8000;
  static const int max_packets = 172;
  short id = 0; // unique identifier of the RAW entry

  for (int ipkt = 1; ipkt <= max_packets; ipkt++)
    {
      int pktid = id_base + ipkt;

      // get packet
      p = event->getPacket(pktid);
      if ( p != 0 )
        {
          // go through all channels, save those that have a hit
          for (int ifemch=0; ifemch<144; ifemch++)
            {
              // check if there is a hit (ie, if all data values are 0)
              if ( (p->iValue(ifemch, 0) + p->iValue(ifemch, 0) + p->iValue(ifemch, 0) +
                    p->iValue(ifemch, 0) + p->iValue(ifemch, 0)) == 0 )
                {
                  continue;
                }

              // now get the swkey and detector type (hwkey is not used?)
              short detector_type = -1;
              int arm, sector, iy, iz;
              int swkey = -1;
              int module = p->iValue(0, "MODULE");
              if (module != ipkt)
                {
                  cout << "mEmcPRDFToRaw: error, FEM and packetid mismatch" << endl;
                }

              // module number from packet should always match packetid-8000, but if not
              // we use the packetid as the correct identity
              if (ipkt <= 72)
                {
                  arm = 0;
                  sector = (ipkt-1)/18;
                }
              else
                {
                  arm = 1;
                  if (ipkt <= 108)
                    {
                      sector = (ipkt-73)/18 + 2;
                      iy = ((ipkt-(sector*18)-1)/6)*12 + FeeToSoftwareMap[ifemch]/12;
                      iz = ((ipkt-(sector*18)-1)%6)*12 + FeeToSoftwareMap[ifemch]%12;
                    }
                  else             sector = (ipkt-109)/32;
                }

              if (ipkt <= 108)
                {
                  iy = ( ((ipkt-1)%18)/6 )*12 + FeeToSoftwareMap[ifemch]/12;
                  iz = ((ipkt - 1)%6)*12 + FeeToSoftwareMap[ifemch]%12;
                }
              else
                {
                  iy = ( ((ipkt-109)%32)/8 )*12 + FeeToSoftwareMap[ifemch]/12;
                  iz = ((ipkt - 109)%8)*12 + FeeToSoftwareMap[ifemch]%12;
                }

              swkey = arm*100000 + sector*10000 + iy*100 + iz;

              if (ipkt < 109) detector_type = 1;        // PbSc
              else            detector_type = 2;        // PbGl

              demcraw->set_id(id, id+1);        // starts counting at 1
              demcraw->set_evno(id, p->iValue(0,"EVTNR") );
              demcraw->set_hwkey(id, 0);
              demcraw->set_swkey(id, swkey);
              demcraw->set_type(id, detector_type);
              demcraw->set_adclopre(id,  p->iValue(ifemch, 3) ^ 0xfff );
              demcraw->set_adclopost(id, p->iValue(ifemch, 1) ^ 0xfff );
              demcraw->set_adchipre(id,  p->iValue(ifemch, 4) ^ 0xfff );
              demcraw->set_adchipost(id, p->iValue(ifemch, 2) ^ 0xfff );
              demcraw->set_tdc(id,       p->iValue(ifemch, 0) ^ 0xfff );

              id++;
            }

          delete p;     // free the sucker
        }

    }

  demcraw->SetRowCount( id );

  return 0;
}