PHG4RICHDetector.cc

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

// $$Id: PHG4RICHDetector.cc,v 1.1 2013/10/01 00:33:00 jinhuang Exp $$

#include "PHG4RICHDetector.h"
#include "PHG4RICHSteppingAction.h"

#include <g4main/PHG4Detector.h>  // for PHG4Detector

#include <Geant4/G4ProductionCuts.hh>  // for G4ProductionCuts
#include <Geant4/G4Region.hh>  // for G4Region
#include <Geant4/G4RegionStore.hh>

#include <boost/foreach.hpp>

#include <map>
#include <sstream>

class PHCompositeNode;

using namespace std;
using namespace ePHENIXRICH;

PHG4RICHDetector::PHG4RICHDetector(PHG4Subsystem *subsys, PHCompositeNode *Node, const string &dnam, const RICH_Geometry &g)
  : PHG4Detector(subsys, Node, dnam)
  , ePHENIXRICHConstruction(subsys, g)
  , _region(nullptr)
{
}

PHG4RICHDetector::PHG4RICHDetector(PHG4Subsystem *subsys, PHCompositeNode *Node, const string &dnam)
  : PHG4Detector(subsys, Node, dnam)
  , ePHENIXRICHConstruction(subsys)
  , _region(nullptr)
{
}

void PHG4RICHDetector::ConstructMe(G4LogicalVolume *logicWorld)
{
  const G4RegionStore *theRegionStore = G4RegionStore::GetInstance();
  G4ProductionCuts *gcuts = new G4ProductionCuts(*(theRegionStore->GetRegion("DefaultRegionForTheWorld")->GetProductionCuts()));
  _region = new G4Region("REGION_RICH");
  _region->SetProductionCuts(gcuts);
  _region->SetRegionalSteppingAction(new PHG4RICHSteppingAction(this));

  ePHENIXRICHConstruction::Construct_RICH(logicWorld);

  BOOST_FOREACH (map_log_vol_t::value_type &vol_pair, map_log_vol)
  {
    _region->AddRootLogicalVolume(vol_pair.second);
  }
}