PHG4FPbScSubsystem.cc

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file PHG4FPbScSubsystem.cc

#include "PHG4FPbScSubsystem.h"

#include <g4detectors/PHG4EventActionClearZeroEdep.h>
#include "PHG4FPbScDetector.h"
#include "PHG4FPbScRegionSteppingAction.h"
#include "PHG4FPbScSteppingAction.h"

#include <g4main/PHG4HitContainer.h>
#include <g4main/PHG4SteppingAction.h>  // for PHG4SteppingAction
#include <g4main/PHG4Subsystem.h>       // for PHG4Subsystem

#include <phool/PHCompositeNode.h>
#include <phool/PHIODataNode.h>    // for PHIODataNode
#include <phool/PHNode.h>          // for PHNode
#include <phool/PHNodeIterator.h>  // for PHNodeIterator
#include <phool/PHObject.h>        // for PHObject
#include <phool/getClass.h>

#include <Geant4/G4UserSteppingAction.hh>  // for G4UserSteppingAction

#include <sstream>

class PHG4Detector;

using namespace std;

//_______________________________________________________________________
PHG4FPbScSubsystem::PHG4FPbScSubsystem(const string& name)
  : PHG4Subsystem(name)
  , detector_(nullptr)
  , steppingAction_(nullptr)
  , eventAction_(nullptr)
  , x_position(0)
  , y_position(0)
  , z_position(0)
{
}

//_______________________________________________________________________
int PHG4FPbScSubsystem::Init(PHCompositeNode* topNode)
{
  std::ostringstream hitnodename;
  hitnodename << "G4HIT_" << Name();

  // create hit list
  PHG4HitContainer* fcal_hits = findNode::getClass<PHG4HitContainer>(topNode, ("G4HIT_" + Name()).c_str());
  if (!fcal_hits)
  {
    PHNodeIterator iter(topNode);
    PHCompositeNode* dstNode = dynamic_cast<PHCompositeNode*>(iter.findFirst("PHCompositeNode", "DST"));
    dstNode->addNode(new PHIODataNode<PHObject>(fcal_hits = new PHG4HitContainer(("G4HIT_" + Name())), ("G4HIT_" + Name()), "PHObject"));
  }

  // create detector
  detector_ = new PHG4FPbScDetector(this, topNode, Name());
  detector_->set_Place(x_position, y_position, z_position);
  detector_->Verbosity(Verbosity());

  // create stepping action
  steppingAction_ = new PHG4FPbScSteppingAction(detector_);

  eventAction_ = new PHG4EventActionClearZeroEdep(topNode, hitnodename.str());

  return 0;
}

//_______________________________________________________________________
int PHG4FPbScSubsystem::process_event(PHCompositeNode* topNode)
{
  if (PHG4FPbScRegionSteppingAction* p = dynamic_cast<PHG4FPbScRegionSteppingAction*>(detector_->GetSteppingAction()))
    p->SetInterfacePointers(topNode);
  else
    steppingAction_->SetInterfacePointers(topNode);

  return 0;
}

//_______________________________________________________________________
PHG4Detector* PHG4FPbScSubsystem::GetDetector(void) const
{
  return detector_;
}

//_______________________________________________________________________
PHG4SteppingAction* PHG4FPbScSubsystem::GetSteppingAction(void) const
{
  return steppingAction_;
}