PHG4OHCalDetector.h

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

// Tell emacs that this is a C++ source
//  -*- C++ -*-.
#ifndef G4OHCAL_PHG4OHCALDETECTOR_H
#define G4OHCAL_PHG4OHCALDETECTOR_H

#include <g4main/PHG4Detector.h>

#include <cmath>  // for NAN
#include <map>
#include <set>
#include <string>  // for string
#include <tuple>

class G4AssemblyVolume;
class G4LogicalVolume;
class G4VPhysicalVolume;
class PHCompositeNode;
class PHG4OHCalDisplayAction;
class PHG4OHCalFieldSetup;
class PHParameters;
class PHG4Subsystem;
class PHG4GDMLConfig;
class RawTowerGeomContainer;

class PHG4OHCalDetector : public PHG4Detector
{
 public:
  PHG4OHCalDetector(PHG4Subsystem *subsys, PHCompositeNode *Node, PHParameters *params, const std::string &dnam);

  ~PHG4OHCalDetector() override;

  void ConstructMe(G4LogicalVolume *world) override;

  void Print(const std::string &what = "ALL") const override;


  int IsInOHCal(G4VPhysicalVolume *) const;

  void SuperDetector(const std::string &name) { m_SuperDetector = name; }
  const std::string SuperDetector() const { return m_SuperDetector; }
  int get_Layer() const { return m_Layer; }
  int ConsistencyCheck() const;
  std::tuple<int, int, int> GetRowColumnId(G4VPhysicalVolume *volume) const;

 private:
  void AddGeometryNode();
  int ConstructOHCal(G4LogicalVolume *hcalenvelope);
  int map_towerid(const int tower_id);
  int map_layerid(const unsigned int isector, const int layer_id);
  std::tuple<int, int, int> ExtractLayerTowerId(const unsigned int isector, G4VPhysicalVolume *volume);
  PHG4OHCalDisplayAction *m_DisplayAction = nullptr;
  PHG4OHCalFieldSetup *m_FieldSetup = nullptr;
  PHParameters *m_Params = nullptr;
  G4AssemblyVolume *m_ScintiMotherAssembly = nullptr;
  G4AssemblyVolume *m_ChimScintiMotherAssembly = nullptr;
  double m_InnerRadius = NAN;
  double m_OuterRadius = NAN;
  double m_SizeZ = NAN;
  double m_VolumeEnvelope = NAN;
  double m_VolumeSteel = 0;
  double m_VolumeScintillator = 0;

  int m_NumScintiPlates = -9999;

  int m_ActiveFlag = 0;
  int m_AbsorberActiveFlag = 0;

  int m_Layer = 0;
  std::string m_SuperDetector;
  std::set<G4LogicalVolume *> m_SteelAbsorberLogVolSet;
  std::set<G4LogicalVolume *> m_ScintiTileLogVolSet;
  std::map<G4VPhysicalVolume *, std::tuple<int, int, int>> m_ScintiTilePhysVolMap;

  PHG4GDMLConfig *gdml_config = nullptr;

  std::string m_GDMPath;
  RawTowerGeomContainer *m_RawTowerGeom = nullptr;
  std::string m_TowerGeomNodeName;
};

#endif  // G4OHCAL_PHG4OHCALDETECTOR_H