PHG4BlockCellReco.h

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

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

#include <phparameter/PHParameterContainerInterface.h>

#include <fun4all/SubsysReco.h>

#include <map>
#include <set>
#include <string>
#include <utility>  // for pair

class PHCompositeNode;

class PHG4BlockCellReco : public SubsysReco, public PHParameterContainerInterface
{
 public:
  PHG4BlockCellReco(const std::string &name = "BLOCKRECO");

  ~PHG4BlockCellReco() 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 &d) { detector = d; }
  void etaxsize(const int i, const double deltaeta, const double deltax);
  void checkenergy(const int i = 1) { chkenergyconservation = i; }

  void set_timing_window(const int detid, const double tmin, const double tmax);

 protected:
  void set_size(const int i, const double sizeA, const double sizeB, const int what);
  int CheckEnergy(PHCompositeNode *topNode);

  double sum_energy_g4hit;
  std::map<int, int> binning;
  std::map<int, std::pair<double, double> > cell_size;  // cell size in x/z
  std::map<int, std::pair<double, double> > zmin_max;   // zmin/zmax for each layer for faster lookup
  std::map<int, double> xstep;
  std::map<int, double> etastep;
  std::map<int, std::pair<double, double> > tmin_max;
  std::set<int> implemented_detid;
  std::string detector;
  std::string hitnodename;
  std::string cellnodename;
  std::string geonodename;
  std::string seggeonodename;
  std::map<int, std::pair<int, int> > n_x_z_bins;
  int chkenergyconservation;
};

#endif