AFPSiRowColToLocalCSTool.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef AFP_SICLUSTERTOOLS_AFPSIROWCOLTOLOCALCSTOOL_H
#define AFP_SICLUSTERTOOLS_AFPSIROWCOLTOLOCALCSTOOL_H 1

 
 
 
// FrameWork includes
#include "AthenaBaseComps/AthAlgTool.h"
#include "GaudiKernel/ServiceHandle.h"
 
#include "AFP_SiClusterTools/IAFPSiRowColToLocalCSTool.h"
#include "AFP_DBTools/SiLocAlignData.h"
#include "AFP_DBTools/SiGlobAlignData.h"
 
// using root transformations because CLHEP does not allow to directly set transformation matrix 
#include <Math/Point3D.h>
#include <Math/RotationY.h>
#include <Math/RotationZYX.h>
#include <Math/Translation3D.h>
#include <Math/Transform3D.h>
 
 
#include "CLHEP/Units/SystemOfUnits.h"

class AFPSiRowColToLocalCSTool
  : public extends<AthAlgTool, IAFPSiRowColToLocalCSTool>
{
public:
  AFPSiRowColToLocalCSTool (const std::string& type,
                        const std::string& name,
                        const IInterface* parent);

  virtual StatusCode initialize() override;

  virtual StatusCode finalize() override;
 
  
  virtual xAOD::AFPSiHitsCluster* newXAODLocal (const int stationID, const int layerID, const AFP::SiLocAlignData& LA, const AFP::SiGlobAlignData& GA, const AFPSiClusterBasicObj& cluster, std::unique_ptr<xAOD::AFPSiHitsClusterContainer>& xAODContainer) const override;
 

  void fillXAOD (const int stationID, const int layerID, const ROOT::Math::XYZPoint& position, const ROOT::Math::XYZPoint& positionError, xAOD::AFPSiHitsCluster* xAODCluster) const;
  
private:
  Gaudi::Property<std::vector<int>> m_numberOfLayersInStations {this, "numberOfLayersInStations",{4,4,4,4}, "Array defining number of stations and layers in each station. The size of the array in argument defines number of stations, the value in each entry defines number of layers in the station."};
  
  std::vector< std::vector<ROOT::Math::Transform3D> > m_transformations; 

  Gaudi::Property<std::vector<double>> m_transformationsInit{this, "transformations", {}, "Initialisation of the transformation matrices (4x3). It is not full transformation only part. Must be defined simultaneously with transformationsErr."};

  std::vector< std::vector<ROOT::Math::Transform3D> > m_transformationsErr;

  Gaudi::Property<std::vector<double>> m_transformationsErrInit{this, "transformationsErr", {}, "Initialisation of the error transformation matrices (4x3). It is not full transformation only part. Must be defined simultaneously with transformations."};
 

  Gaudi::Property<double> m_pixelHorizSize {this, "pixelHorizSize", 0.05, "Size of the pixel in horizontal direction when mounted (default = 0.05 mm)"};

  Gaudi::Property<double>  m_pixelVertSize {this, "pixelVertSize", 0.25, "Size of the pixel in vertical direction when mounted (default = 0.25 mm)"};

  void initTransformationMatricesSize (std::list<std::vector< std::vector<ROOT::Math::Transform3D> >* >& matrices, const std::vector<int>& layersInStations);
};
 
#endif