CylinderGeomMicromegas.h

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file CylinderGeomMicromegas.h

// Tell emacs that this is a C++ source
//  -*- C++ -*-.

#ifndef MICROMEGAS_CYLINDERGEOMMICROMEGAS_H
#define MICROMEGAS_CYLINDERGEOMMICROMEGAS_H

#include "MicromegasDefs.h"
#include "MicromegasTile.h"

#include <g4detectors/PHG4CylinderGeom.h>

#include <TGeoMatrix.h>

#include <cmath>
#include <iostream>

class TVector3;
class PHG4Hit;

class CylinderGeomMicromegas : public PHG4CylinderGeom
{
  public:

  //* empty constructor
  explicit CylinderGeomMicromegas()
  {}

  //* constructor
  CylinderGeomMicromegas(int layer )
    : m_layer( layer )
  {}

  void identify(std::ostream&) const override;


  int get_layer() const override {return m_layer;}
  double get_radius() const override {return m_radius;}
  double get_thickness() const override { return m_thickness;}
  double get_zmin() const override {return m_zmin;}
  double get_zmax() const override {return m_zmax;}
  double get_pitch() const { return m_pitch; }

  MicromegasDefs::SegmentationType get_segmentation_type() const {return m_segmentation_type;}

  MicromegasDefs::DriftDirection get_drift_direction() const {return m_drift_direction;}

  // check if hit radius matches this cylinder
  bool check_radius( const TVector3& ) const;


  TVector3 get_local_from_world_coords( uint tileid, const TVector3& ) const;

  TVector3 get_local_from_world_vect( uint tileid, const TVector3& ) const;


  TVector3 get_world_from_local_coords( uint tileid, const TVector3& ) const;

  TVector3 get_world_from_local_vect( uint tileid, const TVector3& ) const;

  int find_tile_cylindrical( const TVector3& ) const;

  int find_tile_planar( const TVector3& ) const;

  size_t get_tiles_count() const { return m_tiles.size(); }

  const MicromegasTile& get_tile( uint tileid ) const
  {
    assert( tileid < m_tiles.size() );
    return m_tiles[tileid];
  }


  void convert_to_planar( uint tileid, PHG4Hit* ) const;

  int find_strip_from_world_coords( uint tileid, const TVector3& ) const;

  int find_strip_from_local_coords( uint tileid, const TVector3& ) const;

  double get_strip_length( uint tileid ) const;

  uint get_strip_count( uint tileid ) const;

  TVector3 get_local_coordinates( uint tileid, uint stripnum ) const;

  TVector3 get_world_coordinates( uint tileid, uint stripnum ) const;

  double get_center_phi( uint tileid ) const
  {
    assert( tileid < m_tiles.size() );
    return m_tiles[tileid].m_centerPhi;
  }

  static constexpr double reference_radius = 82;



  void set_layer(const int i) override {m_layer = i;}
  void set_radius(const double value) override {m_radius = value;}
  void set_thickness(const double value) override {m_thickness = value;}
  void set_zmin(const double value) override {m_zmin = value;}
  void set_zmax(const double value) override {m_zmax = value;}
  void set_pitch( double value ) { m_pitch = value; }

  void set_tiles( const MicromegasTile::List& tiles ) {m_tiles = tiles;}

  void set_segmentation_type( MicromegasDefs::SegmentationType value ) {m_segmentation_type = value;}

  void set_drift_direction( MicromegasDefs::DriftDirection value ) {m_drift_direction = value;}

  private:

  // get local to master transformation matrix for a given tile
  TGeoHMatrix transformation_matrix( uint tileid ) const;

  int m_layer = 0;

  MicromegasDefs::SegmentationType m_segmentation_type = MicromegasDefs::SegmentationType::SEGMENTATION_PHI;

  MicromegasDefs::DriftDirection m_drift_direction = MicromegasDefs::DriftDirection::OUTWARD;

  double m_radius = 0;

  double m_thickness = 0;

  double m_zmin = 0;

  double m_zmax = 0;

  double m_pitch = 0.1;

  MicromegasTile::List m_tiles;

  ClassDefOverride(CylinderGeomMicromegas, 1)
};

#endif