mEmcGeaClusterEvalModule.C

Go to the documentation of this file.

#include <iostream>
#include "PHNode.h"
#include "PHCompositeNode.h"
#include "PHIODataNode.h"
#include "PHNodeIterator.h"
#include "PHTable.hh"
#include "mEmcGeaClusterEvalModule.h"

#include "dEmcEventWrapper.h"

#include "dEmcGeaTrackWrapper.h"

#include "dEmcGeaTowerTrackWrapper.h"

#include "dEmcClusterLocalExtWrapper.h"

#include "dEmcGeaTrackClusterWrapper.h"

#include "dEmcGeaClusterTrackWrapper.h"
typedef PHIODataNode<PHTable> TableNode_t;

#include "phool.h"

using namespace std;

// Default constructor and destructor to pacify CINT
mEmcGeaClusterEvalModule::mEmcGeaClusterEvalModule(){}

mEmcGeaClusterEvalModule::~mEmcGeaClusterEvalModule(){}

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

 emcNode = static_cast<PHCompositeNode*>(i.findFirst("PHCompositeNode", "EMC"));
 if (!emcNode) {
   cerr << PHWHERE << " Could not find EMC node" << endl;
   return False;
   //   emcNode = new PHCompositeNode("EMC");
   //root->addNode(emcNode);
 }

 dstNode = static_cast<PHCompositeNode*>(i.findFirst("PHCompositeNode", "DST"));
 if (!dstNode) {
   cerr << PHWHERE << " Could not find DST node " << endl;
   return False;
   //   dstNode = new PHCompositeNode("DST");
   //   root->addNode(dstNode);
 }

 evaNode = static_cast<PHCompositeNode*>(i.findFirst("PHCompositeNode", "EVA"));
 if (!evaNode) {
   cerr << PHWHERE << " Could not find EVA node " << endl;
   return False;
   //   evaNode = new PHCompositeNode("EVA");
   //   root->addNode(evaNode);
 }


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

  outNode = emcNode;
  n = i.findFirst("PHIODataNode", "dEmcEvent");
  if (!n) {
    cerr<< PHWHERE << "ERROR:  'in' parameter dEmcEvent not found" << endl;
//     w = new dEmcEventWrapper("dEmcEvent", 1);
//     if (!w) {
//       return 1;
//     }
//     n = new TableNode_t(w,"dEmcEvent");
//     outNode->addNode(n);
    return False;
  }
  nodes.append(n);

  outNode = evaNode;
  n = i.findFirst("PHIODataNode", "dEmcGeaTrack");
  if (!n) {
    cerr << PHWHERE << "ERROR:  'in' parameter dEmcGeaTrack not found" << endl;
//      w = new dEmcGeaTrackWrapper("dEmcGeaTrack", 10);
//      if (!w) {
//        return 1;
//      }
//      n = new TableNode_t(w,"dEmcGeaTrack");
//      outNode->addNode(n);
  }
  nodes.append(n);

  outNode = evaNode;
  n = i.findFirst("PHIODataNode", "dEmcGeaTowerTrack");
  if (!n) {
    cerr << PHWHERE << "ERROR:  'in' parameter dEmcGeaTowerTrack not found" << endl;
//      w = new dEmcGeaTowerTrackWrapper("dEmcGeaTowerTrack", 15000);
//      if (!w) {
//        return 1;
//      }
//      n = new TableNode_t(w,"dEmcGeaTowerTrack");
//      outNode->addNode(n);
    return False;
  }
  nodes.append(n);

  outNode = dstNode;
  n = i.findFirst("PHIODataNode", "dEmcClusterLocalExt");
  if (!n) {
    cerr << PHWHERE << "ERROR:  'in' parameter dEmcClusterLocalExt not found" << endl;
//      w = new dEmcClusterLocalExtWrapper("dEmcClusterLocalExt", 6000);
//      if (!w) {
//        return 1;
//      }
//      n = new TableNode_t(w,"dEmcClusterLocalExt");
//      outNode->addNode(n);
    return False;
  }
  nodes.append(n);

  outNode = evaNode;
  j = new PHNodeIterator(outNode);
  if (!(d = static_cast<TableNode_t*>(j->findFirst("PHIODataNode","dEmcGeaTrackCluster")))) {
    cerr << PHWHERE << " Could not find dEmcGeaTrackCluster node" << endl;
    return False;
 //     w = new dEmcGeaTrackClusterWrapper("dEmcGeaTrackCluster", 10);
//      if (!w) {
//        return 1;
//      }
//      d = new TableNode_t(w,"dEmcGeaTrackCluster");
//      outNode->addNode(d);
  }
  delete j;
  nodes.append(d);

  outNode = evaNode;
  j = new PHNodeIterator(outNode);
  if (!(d = static_cast<TableNode_t*>(j->findFirst("PHIODataNode","dEmcGeaClusterTrack")))) {
    cerr << PHWHERE << " Could not find dEmcGeaClusterTrack" << endl;
    return False;

//      w = new dEmcGeaClusterTrackWrapper("dEmcGeaClusterTrack", 10);
//      if (!w) {
//        return 1;
//      }
//      d = new TableNode_t(w,"dEmcGeaClusterTrack");
//      outNode->addNode(d);
  }
  delete j;
  nodes.append(d);

  return callPAM(nodes);
}