KVTrapezoidalStripDetector.h

#ifndef __KVTRAPEZOIDALSTRIPDETECTOR_H
#define __KVTRAPEZOIDALSTRIPDETECTOR_H
 
#include "KVMaterial.h"
#include "KVUnits.h"
#include "TGeoArb8.h"
#include "TGeoCompositeShape.h"
#include "TGeoMatrix.h"
#include "TGeoVolume.h"
class KVTrapezoidalStripDetector : public TGeoVolumeAssembly {
   double frame_thickness = 3 * KVUnits::mm;
   double dX_IF;
 
   TGeoArb8* make_arb8(const TString& name, double A, double B, double C, double dZ, double dX = 0);
   TGeoArb8* make_vert_strip_arb8(const TString& name, double y_A, double y_B, double x_AB, double x_C, double x_D, double dZ, double dX = 0);
   TGeoArb8* make_hori_strip_arb8(const TString& name, double x_AB, double y_A, double x_CD, double y_D, double dZ, double dX = 0);
   double Xcoord_of_frame_bottom(double Ycoord, double A, double B, double C);
   double Ycoord_of_frame_side(double Xcoord, double A, double B, double C);
public:
   enum class StripDirection {
      horizontal,
      vertical
   };
   double get_frame_thickness() const
   {
      return frame_thickness;
   }
 
   KVTrapezoidalStripDetector() : TGeoVolumeAssembly() {}
   KVTrapezoidalStripDetector(const TString& name, double OF_A, double OF_B, double OF_C,
                              double IF_A, double IF_B, double IF_C, double a,
                              double dZ,
                              StripDirection strip_direction, int N, double inter_strip = 0.01);
   double get_X_offset() const
   {
      return dX_IF;
   }
 
   ClassDef(KVTrapezoidalStripDetector, 1) //Single-sided stripped trapezoidal silicon detector
};
 
#endif