LArClusterCollisionTimeAlg.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef LARCLUSTERCOLLISIONTIMEALG_H
#define LARCLUSTERCOLLISIONTIMEALG_H
 
#include "Gaudi/Property.h"
 
#include "AthenaBaseComps/AthReentrantAlgorithm.h"
#include "xAODCaloEvent/CaloClusterContainer.h"
 
#include "LArRecEvent/LArCollisionTime.h"
 
#include "StoreGate/ReadHandleKey.h"
#include "StoreGate/WriteHandleKey.h"
 
 
class LArClusterCollisionTimeAlg : public AthReentrantAlgorithm {
 public:
  //LArClusterCollisionTimeAlg(const std::string& name, ISvcLocator* pSvcLocator);
 
  using AthReentrantAlgorithm::AthReentrantAlgorithm;
 
  ~LArClusterCollisionTimeAlg() = default;
    
  StatusCode          initialize() override final;
  StatusCode          execute(const EventContext& ctx) const override final;
  StatusCode          finalize() override final;
 
 private:
 
  struct perSide_t {
    float time;
    float energy;
    unsigned nClusters;
  public:
  perSide_t() : time(0.), energy(0.),nClusters(0) {} ;
  perSide_t(const float t, const float e, const unsigned nC) : time(t),energy(e),nClusters(nC) {} ;
  };
 
  perSide_t analyseClustersPerSide(std::vector<const xAOD::CaloCluster*>& clusters) const;
 
 
    
  //---------------------------------------------------
  // Member variables
  //---------------------------------------------------
  mutable std::atomic<unsigned> m_nEvt{0};
  mutable std::atomic<unsigned> m_nCollEvt{0};
 
  //---------------------------------------------------
  // Properties
  //--------------------------------------------------- 
  Gaudi::Property<float> m_timeCut { this, "timeDiffCut", 2., "max |A-C| time difference tu pass the filter" };
  Gaudi::Property<size_t> m_maxClusters { this, "maxNClusters", 3, "how many clusters taken into sum" };
 
  SG::ReadHandleKey<xAOD::CaloClusterContainer> m_clusterContainerName{this, "InputName", "LArClusterEM"};
  SG::WriteHandleKey<LArCollisionTime> m_outputName{this, "OutputName", "LArClusterCollTime"};
  
};
#endif