ClusterMakerTool.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
//ClusterMaker.h
// (c) ATLAS Detector software
// Fill the global position fields of the PrepRawData
//
// First version 04/08/2003 Tommaso Lari
//
 
#ifndef INDETRIOMAKER_CLUSTERMAKERTOOL_H
#define INDETRIOMAKER_CLUSTERMAKERTOOL_H
 
#include "AthenaBaseComps/AthAlgTool.h"
 
#include "GeoPrimitives/GeoPrimitives.h"
#include "InDetCondTools/ISiLorentzAngleTool.h"
#include "PixelConditionsData/PixelChargeCalibCondData.h"
#include "PixelConditionsData/PixelOfflineCalibData.h"
#include "PixelReadoutGeometry/IPixelReadoutManager.h"
#include "StoreGate/ReadCondHandleKey.h"
#include "xAODInDetMeasurement/PixelCluster.h"
#include "xAODInDetMeasurement/PixelClusterContainer.h"
 
#include "GaudiKernel/ToolHandle.h"
 
#include <atomic>
#include <vector>
 
template <class T> class ServiceHandle;
class Identifier;
class StatusCode;
 
namespace InDetDD {
  class SiDetectorElement;
}
 
 
class PixelID;
 
namespace InDet {
 
static const InterfaceID IID_ClusterMakerTool("InDet::ClusterMakerTool", 1, 0);
 
class PixelCluster;
class SCT_Cluster;
class SiWidth;
 
class ClusterMakerTool : public AthAlgTool {
 
  // Public methods:
public:
 
  ClusterMakerTool(const std::string &type,
                   const std::string &name,
                   const IInterface *parent);
  ~ClusterMakerTool() = default;
 
  static const InterfaceID& interfaceID() { return IID_ClusterMakerTool; };
 
  StatusCode initialize();
 
 
  // Compute the pixel cluster global position, and the error associated
  // to the position.
  // Called by the pixel clustering tools
  //
  // Input parameters
  // - the cluster Identifier
  // - the position in local reference frame
  // - the list of identifiers of the Raw Data Objects belonging to the cluster
  // - the width of the cluster
  // - the module the cluster belongs to
  // - wheter the cluster contains ganged pixels
  // - the error strategy, currently
  //    0: cluster width/sqrt(12.)
  //    1: pixel pitch/sqrt(12.)
  //    2: parametrized as a function ofpseudorapidity and cluster size
  //       (default)
  //   10: CTB parametrization (as a function of module and cluster size)
  //       no magnetic field
  // - const reference to a PixelID helper class
  PixelCluster  pixelCluster(const Identifier& clusterID,
                             const Amg::Vector2D& localPos,
                             std::vector<Identifier>&& rdoList,
                             const int lvl1a,
                             std::vector<int>&& totList,
                             const SiWidth& width,
                             const InDetDD::SiDetectorElement* element,
                             bool ganged,
                             int errorStrategy,
                             const PixelID& pixelID,
                             bool split,
                             double splitProb1,
                             double splitProb2,
                             const PixelChargeCalibCondData *calibData,
                             const PixelCalib::PixelOfflineCalibData *offlineCalibData,
                             const EventContext& ctx) const;
 
  xAOD::PixelCluster* xAODpixelCluster(xAOD::PixelCluster& cluster,
                       const Amg::Vector2D& localPos,
                       const std::vector<Identifier>& rdoList,
                       const int lvl1a,
                       const std::vector<int>& totList,
                       const SiWidth& width,
                       const InDetDD::SiDetectorElement* element,
                       bool ganged,
                       int errorStrategy,
                       const PixelID& pixelID,
               bool split,
               double splitProb1,
               double splitProb2,
               const PixelChargeCalibCondData *calibData,
               const PixelCalib::PixelOfflineCalibData *offlineCalibData,
               const EventContext& ctx) const;
 
  // Computes global position and errors for SCT cluster.
  // Called by SCT Clustering tools
  //
  // Input parameters
  // - the cluster Identifier
  // - the position in local reference frame
  // - the list of identifiers of the Raw Data Objects belonging to the cluster
  // - the width of the cluster
  // - the module the cluster belongs to
  // - the error strategy, currently
  //    0: Cluster Width/sqrt(12.)
  //    1: Set to a different values for one and two-strip clusters (def.)
 
  SCT_Cluster sctCluster(const Identifier& clusterID,
                         const Amg::Vector2D& localPos,
                         std::vector<Identifier>&& rdoList,
                         const SiWidth& width,
                         const InDetDD::SiDetectorElement* element,
                         int errorStrategy) const;
 
private:
 
  template <typename ClusterType, typename  IdentifierVec, typename ToTList>
  ClusterType makePixelCluster(const Identifier& clusterID,
             const Amg::Vector2D& localPos,
             IdentifierVec&& rdoList,
             const int lvl1a,
             ToTList&& totList,
             const SiWidth& width,
             const InDetDD::SiDetectorElement* element,
             bool ganged,
             int errorStrategy,
             const PixelID& pixelID,
             bool split,
             double splitProb1,
             double splitProb2,
             const PixelChargeCalibCondData *calibData,
             const PixelCalib::PixelOfflineCalibData *offlineCalibData,
             const EventContext& ctx,
             xAOD::PixelCluster* cluster = nullptr) const;
 
  ServiceHandle<InDetDD::IPixelReadoutManager> m_pixelReadout
  {this, "PixelReadoutManager", "PixelReadoutManager", "Pixel readout manager" };
 
  ToolHandle<ISiLorentzAngleTool> m_pixelLorentzAngleTool
  {this, "PixelLorentzAngleTool", "SiLorentzAngleTool/PixelLorentzAngleTool", "Tool to retreive Lorentz angle of Pixel"};
 
  ToolHandle<ISiLorentzAngleTool> m_sctLorentzAngleTool
  {this, "SCTLorentzAngleTool", "SiLorentzAngleTool/SCTLorentzAngleTool", "Tool to retreive Lorentz angle of SCT"};
 
  bool m_forceErrorStrategy1B{false};
 
  // Parametrization of the Pixel errors
  // now moved in PixelConditionsData, except for CTB parametrization
 
  double getPixelCTBPhiError(int layer, int phi, int PhiClusterSize) const;
 
};
 
}
 
#endif // INDETRIOMAKER_CLUSTERMAKERTOOL_H