G4_EPD.C

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

#ifndef COMMON_G4EPD_C
#define COMMON_G4EPD_C

#include <GlobalVariables.C>

#include <g4epd/PHG4EPDSubsystem.h>
#include <g4epd/PHG4EPDModuleReco.h>

#include <g4main/PHG4Reco.h>

#include <fun4all/Fun4AllServer.h>

R__LOAD_LIBRARY(libg4epd.so)

namespace Enable
{
  bool EPD = false;
  bool EPD_TILE = false;
  bool EPD_OVERLAPCHECK = false;
}  // namespace Enable

namespace G4EPD
{
  double dz = 6.;
  double place_z = 316.;
}

void EPDInit()
{
  BlackHoleGeometry::max_radius = std::max(BlackHoleGeometry::max_radius, 91.);
  // using default z-position and add 10 cm for tile thickness
  BlackHoleGeometry::min_z = std::min(BlackHoleGeometry::min_z, -G4EPD::place_z - G4EPD::dz/2. - no_overlapp);
  BlackHoleGeometry::max_z = std::max(BlackHoleGeometry::max_z, G4EPD::place_z + G4EPD::dz/2. + no_overlapp);
}

void EPD(PHG4Reco* g4Reco)
{
  bool overlap_check = Enable::OVERLAPCHECK || Enable::EPD_OVERLAPCHECK;

  PHG4EPDSubsystem* epd = new PHG4EPDSubsystem("EPD");

  epd->SuperDetector("EPD");
  epd->OverlapCheck(overlap_check);
  epd->SetActive();

  g4Reco->registerSubsystem(epd);
}

void EPD_Tiles()
{
    Fun4AllServer *se = Fun4AllServer::instance();
    PHG4EPDModuleReco *TileBuilder = new PHG4EPDModuleReco("EPDTileBuilder");
    TileBuilder->Detector("EPD");
    TileBuilder->set_double_param("epdmpv",2.064e-6);
    se->registerSubsystem(TileBuilder);
}

#endif /* COMMON_G4EPD_C */