EventInfo.cxx

Go to the documentation of this file.

/*
   Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
 */
 
#include <memory>
 
#include "EventInfo/EventInfo.h"
#include "EventInfo/EventID.h"
#include "EventInfo/EventType.h"
#include "EventInfo/TriggerInfo.h"
 
EventInfo::EventInfo()
    : m_event_ID(std::make_unique<EventID>(0, 0)),
      m_event_type(std::make_unique<EventType>()),
      m_trigger_info(std::make_unique<TriggerInfo>()) {}
 
EventInfo::EventInfo(owner<EventID*> id, owner<EventType*> type)
    : EventInfo(std::unique_ptr<EventID>(id),
                std::unique_ptr<EventType>(type)) {}
 
EventInfo::EventInfo(std::unique_ptr<EventID> id,
                     std::unique_ptr<EventType> type)
    : m_event_ID(std::move(id)),
      m_event_type(std::move(type)),
      m_trigger_info(std::make_unique<TriggerInfo>()) {}
 
EventInfo::EventInfo(owner<EventID*> id, owner<EventType*> type,
                     owner<TriggerInfo*> trig_info)
    : EventInfo(std::unique_ptr<EventID>(id), std::unique_ptr<EventType>(type),
                std::unique_ptr<TriggerInfo>(trig_info)) {}
 
EventInfo::EventInfo(std::unique_ptr<EventID> id,
                     std::unique_ptr<EventType> type,
                     std::unique_ptr<TriggerInfo> trig_info)
    : m_event_ID(std::move(id)),
      m_event_type(std::move(type)),
      m_trigger_info(std::move(trig_info)) {}
 
EventInfo::EventInfo(const EventInfo& e)
    : m_event_ID(e.m_event_ID ? std::make_unique<EventID>(*(e.m_event_ID))
                              : std::make_unique<EventID>(0, 0)),
      m_event_type(e.m_event_type
                       ? std::make_unique<EventType>(*(e.m_event_type))
                       : std::make_unique<EventType>()),
      m_trigger_info(e.m_trigger_info
                         ? std::make_unique<TriggerInfo>(*(e.m_trigger_info))
                         : std::make_unique<TriggerInfo>()),
      m_event_flags(e.m_event_flags) {}
 
EventInfo& EventInfo::operator = (const EventInfo& e) {
  if (this != &e) {
    setEventID(e.m_event_ID ? std::make_unique<EventID>(*(e.m_event_ID))
                            : std::make_unique<EventID>(0, 0));
 
    setEventType(e.m_event_type ? std::make_unique<EventType>(*(e.m_event_type))
                                : std::make_unique<EventType>());
 
    setTriggerInfo(e.m_trigger_info
                       ? std::make_unique<TriggerInfo>(*(e.m_trigger_info))
                       : std::make_unique<TriggerInfo>());
 
    m_event_flags = e.m_event_flags;
  }
  return *this;
}
 
// Needs to be defined here as unique_ptr's deleter needs the pointee's type to
// be a complete type when called
EventInfo::~EventInfo() = default;
EventInfo::EventInfo(EventInfo&& e) noexcept = default;
EventInfo& EventInfo::operator=(EventInfo&& e) noexcept = default;
 
void EventInfo::setEventID(owner<EventID*> pid) {
  m_event_ID.reset(pid);
}
 
void EventInfo::setEventID(std::unique_ptr<EventID> pid) {
  m_event_ID = std::move(pid);
}
 
void EventInfo::setEventType(owner<EventType*> pType) {
  m_event_type.reset(pType);
}
 
void EventInfo::setEventType(std::unique_ptr<EventType> pType) {
  m_event_type = std::move(pType);
}
 
void EventInfo::setTriggerInfo(owner<TriggerInfo*> pTrig) {
  m_trigger_info.reset(pTrig);
}
 
void EventInfo::setTriggerInfo(std::unique_ptr<TriggerInfo> pTrig) {
  m_trigger_info = std::move(pTrig);
}
 
bool
EventInfo::setEventFlags(EventFlagSubDet subDet, unsigned int flags) {
  // If event_flags has been cleared on input, reset size to just
  // before lumi
  if (m_event_flags.size() < (unsigned int) Lumi) m_event_flags.resize(Lumi, 0);
  if ((unsigned) subDet < m_event_flags.size()) {
    // Combine potentially-existing error states with new flags
    m_event_flags[subDet] = (EF_ERROR_BITS & m_event_flags[subDet]) | (EF_BITS & flags);
    return(true);
  }
  return(false);
}
 
bool
EventInfo::setEventFlagBit(EventFlagSubDet subDet, unsigned char bit, bool set) {
  // If event_flags has been cleared on input, reset size to just
  // before lumi
  if (m_event_flags.size() < (unsigned int) Lumi) m_event_flags.resize(Lumi, 0);
  if ((unsigned) subDet < m_event_flags.size() && bit < EF_ERROR_SHIFT) {
    if (set) {
      m_event_flags[subDet] = (EF_ERROR_BITS & m_event_flags[subDet]) |
                              ((1u << bit) | m_event_flags[subDet]);
    } else {
      unsigned int mask = EF_BITS ^ (1u << bit);
      m_event_flags[subDet] = (EF_ERROR_BITS & m_event_flags[subDet]) |
                              (mask & m_event_flags[subDet]);
    }
    return(true);
  }
  return(false);
}
 
bool
EventInfo::setErrorState(EventFlagSubDet subDet, EventFlagErrorState errorState) {
  // If event_flags has been cleared on input, reset size to just
  // before lumi
  if (m_event_flags.size() < (unsigned int) Lumi) m_event_flags.resize(Lumi, 0);
  // check that subdet index is in range and that we have a valid error state
  if ((unsigned) subDet < m_event_flags.size() &&
      (errorState == NotSet || errorState == Warning || errorState == Error)) {
    // Combine potentially-existing flags with new error states
    m_event_flags[subDet] =
      (unsigned int) ((EF_BITS & m_event_flags[subDet]) | ((unsigned int) errorState << EF_ERROR_SHIFT));
    return(true);
  }
  return(false);
}
 
void
EventInfo::setActualInteractionsPerCrossing(float interactions) {
  // pack interaction into 16 bits saving the precision to the 0.01
  // level
  unsigned int ivalue = (unsigned int) (interactions * 100.f);
 
  if (ivalue > ACTUAL_INTERACTS_BITS) {
    std::cout <<
    "EventInfo::setActualInteractionsPerCrossing ***  ERROR ***: trying to store too large of value: interactions/ivalue/max value: "
              << interactions << " " << ivalue << " " << ACTUAL_INTERACTS_BITS << std::endl;
  }
  ivalue &= ACTUAL_INTERACTS_BITS;
  // If event_flags has been cleared on input, reset size AND INCLUDE lumi
  if (m_event_flags.size() < nDets) m_event_flags.resize(nDets, 0);
  m_event_flags[Lumi] |= ivalue; // store lower 16 bits
}
 
void
EventInfo::setAverageInteractionsPerCrossing(float interactions) {
  // pack interaction into 16 bits saving the precision to the 0.01
  // level
  unsigned int ivalue = (unsigned int) (interactions * 100.f);
 
  if (ivalue > ACTUAL_INTERACTS_BITS) {
    std::cout <<
    "EventInfo::setAverageInteractionsPerCrossing ***  ERROR ***: trying to store too large of value: interactions/ivalue/max value: "
              << interactions << " " << ivalue << " " << ACTUAL_INTERACTS_BITS << std::endl;
  }
  ivalue = (ivalue << AVE_INTERACTS_SHIFT) & AVE_INTERACTS_BITS;
  // If event_flags has been cleared on input, reset size AND INCLUDE lumi
  if (m_event_flags.size() < nDets) m_event_flags.resize(nDets, 0);
  m_event_flags[Lumi] |= ivalue; // store upper 16 bits
}