PHG4EPDModuleReco.h

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

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

#include <phparameter/PHParameterInterface.h>

#include <fun4all/SubsysReco.h>

#include <array>
#include <cmath>
#include <string>

class PHCompositeNode;

class PHG4EPDModuleReco : public SubsysReco, public PHParameterInterface
{
 public:
  PHG4EPDModuleReco(const std::string &name = "EpdModuleReco");

  ~PHG4EPDModuleReco() override {}

  int InitRun(PHCompositeNode *topNode) override;

  int process_event(PHCompositeNode *topNode) override;

  int ResetEvent(PHCompositeNode * /*topNode*/) override;

  void SetDefaultParameters() override;

  void Detector(const std::string &detector);

  std::string
  get_epd_sim_tower_node_prefix() const
  {
    return m_EPDSimTowerNodePrefix;
  }

  void
  set_epd_sim_tower_node_prefix(const std::string &epdsimTowerNodePrefix)
  {
    m_EPDSimTowerNodePrefix = epdsimTowerNodePrefix;
  }

  std::string
  get_epd_calib_tower_node_prefix() const
  {
    return m_EPDCalibTowerNodePrefix;
  }

  void
  set_epd_calib_tower_node_prefix(const std::string &epdcalibTowerNodePrefix)
  {
    m_EPDCalibTowerNodePrefix = epdcalibTowerNodePrefix;
  }

  void set_timing_window(const double tmi, const double tma);

 private:
  int Getrmap(int rindex);
  int Getphimap(int phiindex);
  float GetTilePhi(int thisphi);
  float GetTilePhi0(int thisphi0);
  float GetTileR(int thisr);
  float GetTileZ(int thisz);
  void CreateNodes(PHCompositeNode *topNode);

  std::string m_Detector;
  std::string m_Hitnodename;
  std::string m_EPDSimTowerNodePrefix = "SIM";
  std::string m_EPDCalibTowerNodePrefix = "CALIB";

  std::string m_TowerInfoNodeName;
  std::string m_TowerInfoNodeName_calib;
  std::array<std::array<std::array<double, 31>, 12>, 2> m_EpdTile_e = {};
  std::array<std::array<std::array<double, 31>, 12>, 2> m_EpdTile_Calib_e = {};

  double m_EpdMpv = NAN;

  double tmin = NAN;
  double tmax = NAN;
  double m_DeltaT = NAN;
};

#endif