EventBuilder.h

Go to the documentation of this file.

/* Copyright (C) 2021 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Volker Friese [committer] */
 
#ifndef CBM_ALGO_EVENTBUILDER_H
#define CBM_ALGO_EVENTBUILDER_H 1
 
#include "CbmDefs.h"
#include "DigiData.h"
#include "DigiEventSelector.h"
#include "EventBuilderConfig.h"
 
#include <algorithm>
#include <gsl/span>
#include <map>
#include <optional>
#include <string>
#include <vector>
 
#include <xpu/host.h>
 
namespace cbm::algo::evbuild
{
 
  struct EventBuilderDetectorMonitorData {
    size_t nDigisInEvents = 0;  
    size_t nDigis         = 0;  
  };
 
 
  struct EventBuilderMonitorData {
    EventBuilderDetectorMonitorData sts;    
    EventBuilderDetectorMonitorData much;   
    EventBuilderDetectorMonitorData tof;    
    EventBuilderDetectorMonitorData bmon;   
    EventBuilderDetectorMonitorData trd;    
    EventBuilderDetectorMonitorData trd2d;  
    EventBuilderDetectorMonitorData rich;   
    EventBuilderDetectorMonitorData psd;    
    EventBuilderDetectorMonitorData fsd;    
    size_t numTriggers = 0;                 
    size_t numEvents   = 0;                 
    xpu::timings time;                      
  };
 
 
  class EventBuilder {
 
   public:
    typedef std::pair<std::vector<DigiEvent>, EventBuilderMonitorData> resultType;
 
    EventBuilder(const EventBuilderConfig& config) : fConfig(config) {}
 
 
    virtual ~EventBuilder(){};
 
 
    resultType operator()(const DigiData& ts, const std::vector<double> triggers,
                          std::optional<DigiEventSelector> selector) const;
 
 
    std::string ToString() const;
 
 
   private:  // methods
    DigiEvent BuildEvent(const DigiData& ts, EventBuilderMonitorData& monitor, double trigger) const;
 
 
    template<typename Vector>
    static Vector CopyRange(const Vector& source, double tMin, double tMax)
    {
      using Data = typename Vector::value_type;
      auto comp1 = [](const Data& obj, double value) { return obj.GetTime() < value; };
      auto comp2 = [](double value, const Data& obj) { return value < obj.GetTime(); };
      auto lower = std::lower_bound(source.begin(), source.end(), tMin, comp1);
      auto upper = std::upper_bound(lower, source.end(), tMax, comp2);
      return Vector(lower, upper);
    }
 
 
   private:                      // data members
    EventBuilderConfig fConfig;  
  };
 
}  // namespace cbm::algo::evbuild
 
#endif /* CBM_ALGO_EVENTBUILDER_H */