PhotoMultiplierHitDigi_factory.h

Go to the documentation of this file.

// SPDX-License-Identifier: LGPL-3.0-or-later
// Copyright (C) 2022, 2023, Christopher Dilks
 
#pragma once
 
#include <JANA/JEvent.h>
#include <edm4eic/MCRecoTrackerHitAssociationCollection.h>
#include <edm4eic/EDM4eicVersion.h>
#include <edm4eic/RawTrackerHitCollection.h>
#include <memory>
#include <string>
#include <utility>
#include <vector>
 
// algorithms
#include "algorithms/digi/PhotoMultiplierHitDigi.h"
#include "algorithms/digi/PhotoMultiplierHitDigiConfig.h"
// JANA
#include "extensions/jana/JOmniFactory.h"
// services
#include "services/algorithms_init/AlgorithmsInit_service.h"
#include "services/geometry/dd4hep/DD4hep_service.h"
#include "services/geometry/richgeo/RichGeo_service.h"
#include "services/geometry/richgeo/ReadoutGeo.h"
 
namespace eicrecon {
 
class PhotoMultiplierHitDigi_factory
    : public JOmniFactory<PhotoMultiplierHitDigi_factory, PhotoMultiplierHitDigiConfig> {
 
public:
  using AlgoT = eicrecon::PhotoMultiplierHitDigi;
 
private:
  std::unique_ptr<AlgoT> m_algo;
 
  PodioInput<edm4hep::EventHeader> m_event_headers_input{this};
  PodioInput<edm4hep::SimTrackerHit> m_sim_hits_input{this};
  PodioOutput<edm4eic::RawTrackerHit> m_raw_hits_output{this};
#if EDM4EIC_BUILD_VERSION >= EDM4EIC_VERSION(8, 7, 0)
  PodioOutput<edm4eic::MCRecoTrackerHitLink> m_links_output{this};
#endif
  PodioOutput<edm4eic::MCRecoTrackerHitAssociation> m_raw_assocs_output{this};
 
  ParameterRef<std::string> m_detectorName{this, "detectorName", config().detectorName, ""};
  ParameterRef<std::string> m_readoutClass{this, "readoutClass", config().readoutClass, ""};
  ParameterRef<double> m_hitTimeWindow{this, "hitTimeWindow", config().hitTimeWindow, ""};
  ParameterRef<double> m_timeResolution{this, "timeResolution", config().timeResolution, ""};
  ParameterRef<double> m_speMean{this, "speMean", config().speMean, ""};
  ParameterRef<double> m_speError{this, "speError", config().speError, ""};
  ParameterRef<double> m_pedMean{this, "pedMean", config().pedMean, ""};
  ParameterRef<double> m_pedError{this, "pedError", config().pedError, ""};
  ParameterRef<bool> m_enablePixelGaps{this, "enablePixelGaps", config().enablePixelGaps,
                                       "enable/disable removal of hits in gaps between pixels"};
  ParameterRef<double> m_safetyFactor{this, "safetyFactor", config().safetyFactor,
                                      "overall safety factor"};
  ParameterRef<bool> m_enableNoise{this, "enableNoise", config().enableNoise, ""};
  ParameterRef<double> m_noiseRate{this, "noiseRate", config().noiseRate, ""};
  ParameterRef<double> m_noiseTimeWindow{this, "noiseTimeWindow", config().noiseTimeWindow, ""};
  //ParameterRef<std::vector<std::pair<double, double>>> m_quantumEfficiency {this, "quantumEfficiency", config().quantumEfficiency, ""};
 
  Service<AlgorithmsInit_service> m_algorithmsInit{this};
  Service<RichGeo_service> m_RichGeoSvc{this};
 
public:
  void Configure() {
    m_algo = std::make_unique<AlgoT>(GetPrefix());
    m_algo->level(static_cast<algorithms::LogLevel>(logger()->level()));
 
    // Initialize richgeo ReadoutGeo and set random CellID visitor lambda (if a RICH)
    if (GetPluginName() == "DRICH" || GetPluginName() == "PFRICH") {
      m_RichGeoSvc().GetReadoutGeo(config().detectorName, config().readoutClass);
      m_algo->SetVisitRngCellIDs(
          [this](std::function<void(PhotoMultiplierHitDigi::CellIDType)> lambda, float p) {
            m_RichGeoSvc()
                .GetReadoutGeo(config().detectorName, config().readoutClass)
                ->VisitAllRngPixels(lambda, p);
          });
      m_algo->SetPixelGapMask(
          [this](PhotoMultiplierHitDigi::CellIDType cellID, dd4hep::Position pos) {
            return m_RichGeoSvc()
                .GetReadoutGeo(config().detectorName, config().readoutClass)
                ->PixelGapMask(cellID, pos);
          });
    }
 
    m_algo->applyConfig(config());
    m_algo->init();
  }
 
  void Process(int32_t /* run_number */, uint64_t /* event_number */) {
    m_algo->process({m_event_headers_input(), m_sim_hits_input()}, {m_raw_hits_output().get(),
#if EDM4EIC_BUILD_VERSION >= EDM4EIC_VERSION(8, 7, 0)
                                                                    m_links_output().get(),
#endif
                                                                    m_raw_assocs_output().get()});
  }
};
 
} // namespace eicrecon