GGSFitHitsReader.h

Go to the documentation of this file.

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

#ifndef HERD_GGSFITHITSREADER_H_
#define HERD_GGSFITHITSREADER_H_

#ifndef GGSSMARTLOG_H_
#define HERD_GGSFITHITSREADER_H_UNDEF_LOGGING
#endif

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

#ifdef HERD_GGSFITHITSREADER_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/FitGeoParamsColl.h"
#include "dataobjects/FitHitsColl.h"
#include "dataobjects/FitParticleHitsColl.h"

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

using namespace EA;

namespace Herd {

namespace Readers {
// clang-format off
// clang-format on
class GGSFitHitsReader {
public:
  explicit GGSFitHitsReader(std::shared_ptr<GGSTRootReader> ggsReader)
      : m_ggsReader{std::move(ggsReader)}, m_hitReader{m_ggsReader->GetReader<GGSTHitsReader>()} {};

  bool IsFitHitsCollAvailable();

  RetrievalResult GetFitHitsColl(ObjectWrapper &wrapper);

  bool IsFitParticleHitsCollAvailable();

  RetrievalResult GetFitParticleHitsColl(ObjectWrapper &wrapper);

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

  void SetFitGeoParamsColl(observer_ptr<FitGeoParamsColl> fitGeoParamsColl);

  void SetSideDetectorFirstVolumeID(std::array<int, 4> sideFitFirstVolumeID) {
    m_sideFitFirstVolumeID = sideFitFirstVolumeID;
  }

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();

  observer_ptr<FitGeoParamsColl> m_fitGeoParamsColl = nullptr;
  int m_topFitDetIndex = std::numeric_limits<int>::lowest();
  int m_sideFitDetIndex = std::numeric_limits<int>::lowest();
  std::array<int, 4> m_sideFitFirstVolumeID{}; // ID of the first MAT on each lateral FIT, indexed as in GGS files

  FitHitsColl m_fitHitsColl;
  FitParticleHitsColl m_fitParticleHitsColl;

  bool SetupFitReadout();
};
} // namespace Readers
} // namespace Herd

#endif