EventT0.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
#pragma once
 
#include <framework/gearbox/Const.h>
#include <TObject.h>
 
#include <cmath>
#include <optional>
#include <vector>
 
namespace Belle2 {
  class EventT0 : public TObject {
 
  public:
    struct EventT0Component {
      EventT0Component() = default;
 
      EventT0Component(double eventT0_, double eventT0Uncertainty_, const Const::DetectorSet& detectorSet_,
                       const std::string& algorithm_ = "", double quality_ = NAN) :
        eventT0(eventT0_), eventT0Uncertainty(eventT0Uncertainty_), detectorSet(detectorSet_),
        algorithm(algorithm_), quality(quality_) {}
 
      double eventT0 = NAN;
      double eventT0Uncertainty = NAN;
      Const::DetectorSet detectorSet;
      std::string algorithm = "";
      double quality = NAN;
 
      ClassDefNV(EventT0Component, 3) 
    };
 
    bool hasEventT0() const;
 
    double getEventT0() const;
 
    std::optional<EventT0Component> getEventT0Component() const;
 
    double getEventT0Uncertainty() const;
 
    void setEventT0(double eventT0, double eventT0Uncertainty, const Const::DetectorSet& detector, const std::string& algorithm = "");
 
    void setEventT0(const EventT0Component& eventT0);
 
    void addTemporaryEventT0(const EventT0Component& eventT0);
 
    const std::vector<EventT0Component>& getTemporaryEventT0s() const;
 
    const std::vector<EventT0Component> getTemporaryEventT0s(Const::EDetector detector) const;
 
    bool hasTemporaryEventT0(const Const::DetectorSet& detectorSet = Const::allDetectors) const;
 
    Const::DetectorSet getTemporaryDetectors() const;
 
    unsigned long getNumberOfTemporaryEventT0s() const;
 
    std::optional<EventT0Component> getBestSVDTemporaryEventT0() const;
    std::optional<EventT0Component> getBestCDCTemporaryEventT0() const;
    std::optional<EventT0Component> getBestTOPTemporaryEventT0() const;
    std::optional<EventT0Component> getBestECLTemporaryEventT0() const;
 
    bool isEventT0Of(Const::EDetector detector) const
    {
      if (hasTemporaryEventT0(detector)) {
        return m_eventT0.detectorSet.contains(detector);
      }
      return false;
    }
 
    bool isSVDEventT0() const
    {
      return isEventT0Of(Const::EDetector::SVD);
    }
    bool isCDCEventT0() const
    {
      return isEventT0Of(Const::EDetector::CDC);
    }
    bool isTOPEventT0() const
    {
      return isEventT0Of(Const::EDetector::TOP);
    }
    bool isECLEventT0() const
    {
      return isEventT0Of(Const::EDetector::ECL);
    }
 
 
    void clearTemporaries();
 
    void clearEventT0();
 
  private:
    std::vector<EventT0Component> m_temporaryEventT0List;
    EventT0Component m_eventT0;
    bool m_hasEventT0 = false;
 
    ClassDef(EventT0, 5) 
  };
}