TBCaloCoordinate.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
// ***************************************************************************
// Liquid Argon H8 Combined TestBeam Geometry 
// -----------------------------------------
//
// 30.9.2004 Creation of the class by claire.bourdarios@cern.ch
//
// ***************************************************************************
 
#ifndef TBCALOGEOMETRY_TBCALOCOORDINATE_H
#define TBCALOGEOMETRY_TBCALOCOORDINATE_H
 
#include "CaloDetDescr/ICaloCoordinateTool.h"
#include "AthenaBaseComps/AthAlgTool.h"
 
class CaloPhiRange;
class ITBCaloPosTool;
class IMessageSvc;
class TBDetDescrManager;
 
//done in the interface : #include "CLHEP/Geometry/Transform3D.h"
 
//<<<<<< INCLUDES >>>>>>

 
class TBCaloCoordinate : public extends<AthAlgTool, ICaloCoordinateTool>
{
public:
  
  // Constructors:
 
  TBCaloCoordinate(const std::string& type, 
          const std::string& name, 
          const IInterface* parent); 
 
  ~TBCaloCoordinate();
  StatusCode initialize();
  StatusCode finalize();
 
  // Methods for general use: 
 
  // 1) basic translation between the 2 cartesian coordinate systems:
  
  void ctb_to_local( Amg::Vector3D& pt_ctb,
                     Amg::Vector3D& pt_local);
  void local_to_ctb( Amg::Vector3D& pt_local,
                     Amg::Vector3D& pt_ctb);
  
  void ctb_to_local( double& x_ctb, double& y_ctb, double& z_ctb, 
             double& x_local, double& y_local, double& z_local);
  void local_to_ctb( double& x_local, double& y_local, double& z_local, 
             double& x_ctb, double& y_ctb, double& z_ctb);
 
  // 2) the most common use-case : 
  //    the H8 beam goes along the x axis -> that are the local
  //    eta and phi seen by the calorimeter ?
 
  double beam_local_eta();
  double beam_local_phi();
 
 
  // Methods of less general interest : ( see ICaloCoordinate.h for comments )
 
  // real table reading :
  virtual void read_table_position();
 
  // hand-coded value :
  virtual void read_fake_table_position();
 
  Amg::Transform3D* transform_calo_to_ctb();
  Amg::Transform3D* transform_ctb_to_calo();
 
  void print_transform(Amg::Transform3D& htrans);
 
private:
  // Avoid coverity warnings
  TBCaloCoordinate (const TBCaloCoordinate&);
  TBCaloCoordinate& operator= (const TBCaloCoordinate&);
 
  // these methods are used by read_table_position to get numbers :
 
  bool read_data_position();   // for H8 data
  bool read_MC_position();     // for H8 MC where numbers are in POOL
  void read_user_position();   // for H8 MC where numbers are not avail
  void read_neutral_position(); // for Atlas
 
  // these are the famous key numbers :
 
  // 1) position of the table axis (oral tradition)
  double m_table_axis_data;
  double m_table_axis_MC;
  // 2) distance between the table axis and projective point (oral tradition)
  double m_table_proj_data;
  double m_table_proj_MC;
  // 3) angle of rotation : default is set in jobOpt (property)
  //                        + will be overwritten if DB/Pool is used
  //                        + will be ignored "if atlas" 
  double m_table_eta;
 
  // Another layer on top, "in case" the user wants to add an
  // extra HepTransform on top of these known moovements :
  double m_calo_phi_shift;
  double m_calo_theta_shift;
  double m_calo_psi_shift;
  double m_calo_x_shift;
  double m_calo_y_shift;
  double m_calo_z_shift;
 
  // these will be read in the DB :
  double m_table_theta;
  double m_table_z;
  double m_table_delta;
  // this is calculated from the DB info :
  double m_table_calc_theta;
  double m_table_calc_x;
 
  // these properties define what we are reading (see jobOpt)
  int m_DBRead;
  int m_PoolRead;
 
  // keep the 2 steps separated : table rotation, then translation 
  //HepGeom::Rotate3D* m_table_rotate;
  Amg::RotationMatrix3D* m_table_rotate;
  //HepGeom::Translate3D* m_table_shift;
  Amg::Translation3D* m_table_shift;
 
  // the user defined one on top :
 
  Amg::Transform3D m_rotz_extra_calo_to_ctb;
  Amg::Transform3D m_roty_extra_calo_to_ctb;
  Amg::Transform3D m_rotx_extra_calo_to_ctb;
  Amg::Transform3D m_translxyz_extra_calo_to_ctb;
 
  // the final transform from local to ctb :
  Amg::Transform3D* m_transform_calo_to_ctb;
  Amg::Transform3D* m_transform_ctb_to_calo;
 
  // Run number
  long m_runNumber;
 
  // to limit the printout
  int m_firstevt;
 
  ITBCaloPosTool* m_postool;
  const TBDetDescrManager* m_MCmgr;
  CaloPhiRange* m_range;
};
 
 
#endif // TBCALOGEOMETRY_TBCALOCOORDINATE_H