Build your own module - E1039-Collaboration/e1039-wiki GitHub Wiki

#include <g4hit/PHG4HitContainer.h>
#include <fun4all/getClass.h>

Myanalysis::InitRun(PHCompositeNode *topNode) {
    PHNodeIterator iter(topNode);
    PHCompositeNode* eventNode = static_cast<PHCompositeNode*>(iter.findFirst("PHCompositeNode", "DST"));
    if(!eventNode) {/*add some protection*/}

    PHG4HitContainer *g4hits = new PHG4HitContainer();
    PHIODataNode<PHObject>* g4hits_node = new PHIODataNode<PHObject>(g4hits, " G4HIT_HCALIN ", "PHObject");

    eventNode->addNode(g4hits_node);
}

Myanalysis::process_event(PHCompositeNode *topNode){
    PHG4HitContainer *g4hits = findNode::getClass<PHG4HitContainer>(topNode,”G4HIT_HCALIN”);
    if(g4hits) {/*analysis code*/}
}