mEmcFEMToDCMModule.C

Go to the documentation of this file.

#include "mEmcFEMToDCMModule.h"

#include "dEmcFEMDataWrapper.h"

#include "dEmcDCMDataWrapper.h"

#include "PHNode.h"
#include "PHCompositeNode.h"
#include "PHIODataNode.h"
#include "PHNodeIterator.h"
#include "PHTable.hh"

#include <iostream>
using namespace std;

typedef PHIODataNode<PHTable> TableNode_t;

PHBoolean
mEmcFEMToDCMModule::event(PHCompositeNode *root) {
 PHPointerList<PHNode> nodes;
 PHNodeIterator i(root), *j;
 PHNode *n;
 TableNode_t *d;
 PHTable *w;
 PHCompositeNode *dcmNode, *emcNode, *outNode;

 dcmNode = static_cast<PHCompositeNode*>(i.findFirst("PHCompositeNode", "DCM"));
 if (!dcmNode) {
   dcmNode = new PHCompositeNode("DCM");
   root->addNode(dcmNode);
 }

 emcNode = static_cast<PHCompositeNode*>(i.findFirst("PHCompositeNode", "EMC"));
 if (!emcNode) {
   emcNode = new PHCompositeNode("EMC");
   root->addNode(emcNode);
 }

// Insert code here to navigate node hierarchy and find
// or create specific nodes to pass to physics module...

  outNode = emcNode;
  n = i.findFirst("PHIODataNode", "dEmcFEMData");
  if (!n) {
    cout << "ERROR:  'in' parameter dEmcFEMData not found" << endl;
     w = new dEmcFEMDataWrapper("dEmcFEMData", 500);
     if (!w) {
       return 1;
     }
     n = new TableNode_t(w,"dEmcFEMData");
     outNode->addNode(n);
  }
  nodes.append(n);

  outNode = dcmNode;
  j = new PHNodeIterator(outNode);
  if (!(d = static_cast<TableNode_t*>(j->findFirst("PHIODataNode","dEmcDCMData")))) {
     w = new dEmcDCMDataWrapper("dEmcDCMData", 500);
     if (!w) {
       return 1;
     }
     d = new TableNode_t(w,"dEmcDCMData");
     outNode->addNode(d);
  }
  delete j;
  nodes.append(d);

  return callPAM(nodes);
}