GGSScdHitsReader.h

Go to the documentation of this file.

/*
 * GGSScdHitsReader.h
 *
 *  Created on: 27 May 2021
 *      Author: Valerio Formato
 */
 
#ifndef HERD_GGSSCDHITSREADER_H_
#define HERD_GGSSCDHITSREADER_H_
 
#ifndef GGSSMARTLOG_H_
#define HERD_GGSSCDHITSREADER_H_UNDEF_LOGGING
#endif
 
// GGS headers
#include "montecarlo/readers/GGSTHitsReader.h"
#include "montecarlo/readers/GGSTRootReader.h"
 
#ifdef HERD_GGSSCDHITSREADER_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 "herd/analysis/common/DefaultValues.h"
#include "herd/analysis/dataobjects/ScdGeoParamsColl.h"
#include "herd/analysis/dataobjects/ScdHitsColl.h"
#include "herd/analysis/dataobjects/ScdParticleHitsColl.h"
 
// c++ headers
#include <string>
#include <utility>
 
namespace Herd::Readers {
// clang-format off
// clang-format on
class GGSScdHitsReader {
public:
  explicit GGSScdHitsReader(std::shared_ptr<GGSTRootReader> ggsReader)
      : m_ggsReader{std::move(ggsReader)}, m_hitReader{m_ggsReader->GetReader<GGSTHitsReader>()} {};
 
  bool IsScdHitsCollAvailable();
  EA::RetrievalResult GetScdHitsColl(EA::ObjectWrapper &wrapper);
 
  bool IsScdWaferParticleHitsCollAvailable() { return IsScdHitsCollAvailable(); };
 
  EA::RetrievalResult GetScdWaferParticleHitsColl(EA::ObjectWrapper &wrapper);
 
  std::vector<std::string> FreeObjects(const std::vector<std::string> &objs, EA::Memory::Status memStatus);
 
  void SetScdGeoParamsColl(EA::observer_ptr<ScdGeoParamsColl> sdGeoParamsColl);
 
  void SetSideDetectorFirstWaferID(std::array<int, 4> sideScdFirstWaferID) {
    m_sideScdFirstWaferID = sideScdFirstWaferID;
  }
 
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();
 
  EA::observer_ptr<ScdGeoParamsColl> m_sdGeoParamsColl = nullptr;
  int m_topScdDetIndex = std::numeric_limits<int>::lowest();
  int m_sideScdHDetIndex = std::numeric_limits<int>::lowest();
  int m_sideScdVDetIndex = std::numeric_limits<int>::lowest();
  std::array<int, 4> m_sideScdFirstWaferID{}; // ID of the first Wafer on each lateral Scd, indexed as in GGS files
 
  ScdHitsColl m_sdHitsColl;
  ScdParticleHitsColl m_sdWaferParticleHitsColl;
 
  bool m_firstRun = true;
  bool SetupScdReadout();
};
} // namespace Herd::Readers
 
#endif