Back to home page

sPhenix code displayed by LXR

 
 

    


File indexing completed on 2025-08-06 08:17:29

0001 #include "RawTowerGeomv4.h"
0002 
0003 #include <cmath>
0004 #include <iostream>
0005 
0006 RawTowerGeomv4::RawTowerGeomv4(RawTowerDefs::keytype id)
0007   : _towerid(id)
0008 {
0009 }
0010 
0011 double RawTowerGeomv4::get_center_radius() const
0012 {
0013   return std::sqrt(_center_x * _center_x +
0014                    _center_y * _center_y);
0015 }
0016 
0017 double RawTowerGeomv4::get_theta() const
0018 {
0019   double radius = sqrt(_center_x * _center_x + _center_y * _center_y);
0020   double theta = atan2(radius, _center_z);
0021   return theta;
0022 }
0023 
0024 double RawTowerGeomv4::get_eta() const
0025 {
0026   double theta = get_theta();
0027   double eta = -log(tan(theta / 2.));
0028   return eta;
0029 }
0030 
0031 double RawTowerGeomv4::get_phi() const
0032 {
0033   return atan2(_center_y, _center_x);
0034 }
0035 
0036 void RawTowerGeomv4::identify(std::ostream& os) const
0037 {
0038   os << "RawTowerGeomv4:  x: " << get_center_x() << "  y: " << get_center_y() << "  z: " << get_center_z()
0039      << "\n           dx: " << get_size_x() << " dy: " << get_size_y() << " dz: " << get_size_z()
0040      << "\n           tower_type = " << _tower_type << std::endl;
0041 }