CaloFillRectangularCluster.h

Go to the documentation of this file.

// This file's extension implies that it's C, but it's really -*- C++ -*-.
 
/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/
 
// $Id: CaloFillRectangularCluster.h,v 1.7 2009-04-25 17:57:00 ssnyder Exp $
 
 
#ifndef CALOCLUSTERCORRECTION_CALOFILLRECTANGULARCLUSTER_H
#define CALOCLUSTERCORRECTION_CALOFILLRECTANGULARCLUSTER_H
 
 
#include "CaloClusterCorrection/CaloClusterCorrection.h"
#include "CaloInterface/ISetCaloCellContainerName.h"
#include <string>
#include "AthenaKernel/IOVSvcDefs.h"
#include "StoreGate/ReadHandleKey.h"
#include "GaudiKernel/ThreadLocalContext.h"
#include "StoreGate/ReadCondHandleKey.h"
#include "CaloDetDescr/CaloDetDescrManager.h"
 
 
 
// Helper object used for the per-sampling calculations.
namespace CaloClusterCorr {
class SamplingHelper;
}
//class CaloCluster;
class CaloCell;
class CaloCellContainer;
 
 
class CaloFillRectangularCluster
  : public CaloClusterCorrection,
    public ISetCaloCellContainerName
{
public:
  CaloFillRectangularCluster(const std::string& type,
                             const std::string& name,
                             const IInterface* parent);

  virtual StatusCode initialize() override;

  virtual void makeCorrection (const Context& myctx,
                               xAOD::CaloCluster* cluster) const override;
 
 
  // Alternate version that takes an EventContext.
  void makeCorrection (const EventContext& ctx,
                       xAOD::CaloCluster* cluster) const
  {
    return makeCorrection (context(ctx), cluster);
  }
 
 
  /*
   * @brief Return the seed position of a cluster.
   * @param helper Sampling calculation helper object.
   * @param cluster The cluster on which to operate.
   * @param[out] eta The @f$\eta@f$ location of the cluster seed.
   * @param[out] phi The @f$\phi@f$ location of the cluster seed.
   *
   * The cluster seed is the center of rectangular cluster windows.
   * This may be overridden by derived classes to change the seed definition.
   */
  virtual void get_seed (CaloClusterCorr::SamplingHelper& helper,
                         const xAOD::CaloCluster* cluster,
                         double& eta, double& phi) const;
 

  virtual StatusCode setCaloCellContainerName (const std::string& name) override;
 

  typedef std::array<std::pair<double, double>, 4> WindowArray_t;
 

  virtual
  WindowArray_t initWindows (const int neta,
                             const int nphi,
                             const double detas2,
                             const double dphis2) const;
 
 
 
private:

  CaloFillRectangularCluster() = delete;
 
 
  /*
   * @brief Actually make the correction for one region (barrel or endcap).
   * @param ctx     The event context.
   * @param helper Sampling calculation helper object.
   * @param eta The @f$\eta@f$ seed of the cluster.
   * @param phi The @f$\phi@f$ seed of the cluster.
   * @param samplings List of samplings for this region.
   */
  void makeCorrection1 (const EventContext& ctx,
                        const CaloDetDescrManager& dd_man,
                        CaloClusterCorr::SamplingHelper& helper,
                        double eta,
                        double phi,
                        const CaloSampling::CaloSample samplings[4]) const;
 
 
  /*
   * @brief Execute the correction, given a helper object.
   * @param ctx     The event context.
   * @param helper Sampling calculation helper object.
   */
  void makeCorrection2 (const EventContext& ctx,
                        const CaloDetDescrManager& dd_man,
                        CaloClusterCorr::SamplingHelper& helper) const;
 
 
protected:
  double m_deta0 = 0.0, m_deta1 = 0.0, m_deta2 = 0.0, m_deta3 = 0.0;
  double m_dphi0 = 0.0, m_dphi1 = 0.0, m_dphi2 = 0.0, m_dphi3 = 0.0;

  int m_neta, m_nphi;

  bool m_fill_cluster;
 

  bool m_setRawState;
 
private:
  SG::ReadHandleKey<CaloCellContainer> m_cellsName;
  
  
};
 
#endif // not CALOCLUSTERCORRECTION_CALOFILLRECTANGULARCLUSTER_H