GGSPsdHitsReader.h

Go to the documentation of this file.

/*
 * GGSPsdHitsReader.h
 *
 *  Created on: 31 May 2021
 *      Author: Valerio Formato
 */

#ifndef HERD_GGSPSDHITSREADER_H_
#define HERD_GGSPSDHITSREADER_H_

#ifndef GGSSMARTLOG_H_
#define HERD_GGSPSDHITSREADER_H_UNDEF_LOGGING
#endif

// GGS headers
#include "montecarlo/readers/GGSTHitsReader.h"
#include "montecarlo/readers/GGSTRootReader.h"

#ifdef HERD_GGSPSDHITSREADER_H_UNDEF_LOGGING
// Remove the definitions of GGS output macros (would conflict with the same macros defined in EventAnalysis)
// We only do this if GGSSmartLog.h is being parsed in this TU, otherwise there would be nothing to undef
// (or we could even be undefining the EA macros)
#undef COUT
#undef CCOUT
#undef ENDL
#endif

// EA headers
#include "core/ObjectWrapper.h"
#include "core/RetrievalResult.h"
#include "utils/Memory.h"

// HerdSoftware headers
#include "common/DefaultValues.h"
#include "dataobjects/PsdGeoParamsColl.h"
#include "dataobjects/PsdHitsColl.h"
#include "dataobjects/PsdParticleHitsColl.h"

// c++ headers
#include <functional>
#include <string>
#include <utility>

using namespace EA;

namespace Herd {

namespace Readers {
// clang-format off
// clang-format on
class GGSPsdHitsReader {
public:
  enum class DetectorType : int { NONE = -1, TILES = 0, BARS = 1 };

  explicit GGSPsdHitsReader(std::shared_ptr<GGSTRootReader> ggsReader)
      : m_ggsReader{std::move(ggsReader)}, m_hitReader{m_ggsReader->GetReader<GGSTHitsReader>()} {};

  bool IsPsdHitsCollAvailable();

  RetrievalResult GetPsdHitsColl(ObjectWrapper &wrapper);

  bool IsPsdParticleHitsCollAvailable();

  RetrievalResult GetPsdParticleHitsColl(ObjectWrapper &wrapper);

  std::vector<std::string> FreeObjects(const std::vector<std::string> &objs, Memory::Status memStatus);

  void SetPsdGeoParamsColl(observer_ptr<PsdGeoParamsColl> psdGeoParamsColl);

  void SetFileName(const std::string &fileName) { m_fileName = fileName; };

  void SetPsdType(DetectorType type);

private:
  std::shared_ptr<GGSTRootReader> m_ggsReader;
  GGSTHitsReader *m_hitReader = nullptr;
  long long int m_lastHitsReadEvent = std::numeric_limits<long long int>::max();
  long long int m_lastPartHitsReadEvent = std::numeric_limits<long long int>::max();

  enum class SidePsdSegmentation : int { NONE = -1, H = 0, V = 1 };

  struct PsdIDVolInfo {
    unsigned int ID;
    unsigned int layer;
    RefFrame::Side side;
  };

  bool m_readPsdParticleHits = false;

  DetectorType m_type = DetectorType::NONE;
  observer_ptr<PsdGeoParamsColl> m_psdGeoParamsColl = nullptr;
  int m_topPsdDetIndex = std::numeric_limits<int>::lowest();
  std::vector<int> m_sidePsdDetIndexes;
  bool SetupPsdReadout(bool getParticleHits);
  PsdIDVolInfo SidePsdHitIDLayerSideFromGGSID(int iGGSPsdID, SidePsdSegmentation sidePsdSegmentation);

  PsdHitsColl m_psdHitsColl;
  PsdParticleHitsColl m_psdParticleHitsColl;

  std::string m_fileName;
};
} // namespace Readers
} // namespace Herd

#endif