SafeDecorator.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef INDETTRACKPERFMON_SAFEDECORATOR_H
#define INDETTRACKPERFMON_SAFEDECORATOR_H
 
#include "AthContainers/AuxElement.h"
#include "StoreGate/WriteDecorHandleKey.h"
#include "StoreGate/WriteDecorHandle.h"
#include "StoreGate/ReadDecorHandleKey.h"
#include "StoreGate/ReadDecorHandle.h"
#include "GaudiKernel/EventContext.h"
#include "CxxUtils/fpcompare.h" // for equal
 
#include <iostream>
#include <utility>
#include <vector>
#include <cmath>
#include <cstdlib>
 
namespace IDTPM {
 
  template< class ContainerType, class VariableType >
  using WriteKeyAccessorPair =
      std::pair< SG::WriteDecorHandleKey<ContainerType>,
                 SG::AuxElement::ConstAccessor<VariableType> >;
  //
  template< class ContainerType, class VariableType >
  using WriteAccessorRefPair =
      std::pair< SG::WriteDecorHandle<ContainerType, VariableType>,
                 SG::AuxElement::ConstAccessor<VariableType>& >;
  //
  template< class ContainerType, class VariableType >
  using OptionalDecoration =
      std::pair< SG::WriteDecorHandle<ContainerType, VariableType>, bool >;
 
  template< class T_Parent, class T_Cont, class T >
  void createDecoratorKeysAndAccessor(
      T_Parent& parent,
      const SG::ReadHandleKey< T_Cont >& container_key,
      const std::string& prefix,
      const std::vector< std::string >& decor_names,
      std::vector< WriteKeyAccessorPair<T_Cont, T > >& decor_out )
  {
    decor_out.clear();
    decor_out.reserve( decor_names.size() );
    for( const std::string& a_decor_name: decor_names ) {
      decor_out.emplace_back(
          SG::WriteDecorHandleKey< T_Cont >( container_key.key()+"."+prefix+a_decor_name ),
          SG::AuxElement::ConstAccessor<T>( prefix+a_decor_name )
      );
      parent.declare( decor_out.back().first );
      decor_out.back().first.setOwner( &parent );
      decor_out.back().first.initialize().ignore();
    }
  }
 
  template< class T_Parent, class T_Cont >
  void createDecoratorKeys(
      T_Parent& parent,
      const SG::ReadHandleKey< T_Cont >& container_key,
      const std::string& prefix,
      const std::vector< std::string >& decor_names,
      std::vector< SG::WriteDecorHandleKey< T_Cont > >& decor_out)
  {
    decor_out.clear();
    decor_out.reserve( decor_names.size() );
    for (const std::string &a_decor_name : decor_names) {
      assert( !a_decor_name.empty() );
      decor_out.emplace_back( container_key.key()+"."+prefix+a_decor_name );
      // need to declare handles, otherwise the scheduler would not
      // pick up the data dependencies introduced by the decorations
      parent.declare( decor_out.back() );
      decor_out.back().setOwner(&parent);
      decor_out.back().initialize().ignore();
    }
  }
 
  /*template<class T_Parent, class T_Cont>
  void addReadDecoratorHandleKeys(T_Parent &parent,
                                  const SG::ReadHandleKey<T_Cont> &container_key,
                                  const std::string &prefix,
                                  const std::vector<std::string> &decor_names,
                                  std::vector<SG::ReadDecorHandleKey<T_Cont> > &decor_out) {
    decor_out.reserve(decor_out.size() + decor_names.size());
    for (const std::string &a_decor_name : decor_names) {
      decor_out.emplace_back( container_key.key()+"."+prefix+a_decor_name);
      parent.declare(decor_out.back());
      decor_out.back().setOwner(&parent);
      decor_out.back().initialize().ignore();
    }
  }*/
 
  template< class T_Cont, class T >
  std::vector< OptionalDecoration< T_Cont,T > >
  createDecoratorsIfNeeded(
      const T_Cont& container,
      const std::vector< WriteKeyAccessorPair<T_Cont, T > >& keys,
      const EventContext& ctx,
      bool verbose=false ) 
  {
    std::vector< OptionalDecoration< T_Cont, T > > out;
    bool all_available = true;
    if( !container.empty() ) {
      std::vector<bool> decorate;
      decorate.reserve( keys.size() );
      for( const WriteKeyAccessorPair< T_Cont, T >& a_key : keys ) {
        decorate.push_back( !a_key.second.isAvailable( *container[0] ) );
        all_available &= !decorate.back();
        if( verbose && !decorate.back() ) {
          std::cout << "WARNING IDTPM::createDecoratorsIfNeeded: Decoration "
                    << a_key.first.key() << " already exists; reject update."
                    << std::endl;
        } // close if( verbose && !decorate.back() )
      } // close WriteKeyAccessorPair loop
 
      if( !all_available ) {
        std::size_t idx = 0;
        out.reserve( keys.size() );
        for( const WriteKeyAccessorPair< T_Cont, T >& a_key : keys ) {
          assert( idx < decorate.size() );
          out.emplace_back(
            SG::WriteDecorHandle< T_Cont,T >( a_key.first, ctx ),
            decorate[idx++]
          );
          if( not out.back().first.isPresent() ) {
            std::stringstream msg;
            msg << "Container " << a_key.first.key()
                << " to be decorated does not exist.";
            throw std::runtime_error( msg.str() );
          } // close if( not out.back().first.isPresent() )
        } // close WriteKeyAccessorPair loop
      } // close if( !all_available )
    } // close if( !container.empty() )
    return out;
  }
 
  template< class T_Cont, class T >
  std::vector< SG::WriteDecorHandle< T_Cont, T > >
  createDecorators(
      const std::vector< SG::WriteDecorHandleKey< T_Cont > >& keys,
      const EventContext& ctx )
  {
    std::vector< SG::WriteDecorHandle< T_Cont, T > > out;
    out.reserve( keys.size() );
    for( const SG::WriteDecorHandleKey< T_Cont >& a_key : keys ) {
      out.emplace_back( a_key, ctx );
      if( not out.back().isValid() ) {
        std::stringstream msg;
        msg << "Failed to create decorator handle " << a_key.key();
        throw std::runtime_error( msg.str() );
      }
    } // close SG::WriteDecorHandleKey loop
    return out;
  }
 
  template< class T_Cont, class T_Cont_Elm, class T >
  void decorateOrWarnIfUnequal(
      const T_Cont_Elm& particle,
      WriteAccessorRefPair< T_Cont, T >& decorator,
      const T& value )
  {
    if( !decorator.second.isAvailable( particle ) ) {
      const T existing = decorator.second( particle );
      if( not CxxUtils::fpcompare::equal( existing, value ) ) {
        std::cout << "WARNING IDTPM::safeDecorator: " << decorator.first.decorKey()
                  << " Already exists on this object with a different value."
                  << std::endl;
      }
    } else {
      decorator.first( particle ) = value;
    }
  }
 
  template< class T_Cont, class T_Cont_Elm, class T >
  void decorateOrRejectQuietly(
      const T_Cont_Elm& particle,
      OptionalDecoration< T_Cont, T >& decorator,
      const T& value)
  {
    if( decorator.second ) {
      decorator.first( particle ) = value;
    }
  }
 
  template< class T_Cont, class T_Cont_Elm, class T >
  void decorate(
    const T_Cont_Elm& particle,
    OptionalDecoration< T_Cont, T >& decorator,
    const T& value)
  {
    decorator.first( particle ) = value;
  }
 
} // namespace IDTPM
 
#endif // > !INDETTRACKPERFMON_SAFEDECORATOR_H