V11Geometry.h

Go to the documentation of this file.

// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
//
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.
 
 
#ifndef ALICEO2_ITS_V11GEOMETRY_H_
#define ALICEO2_ITS_V11GEOMETRY_H_
 
#include <TMath.h>   // for DegToRad, Cos, Sin, Tan
#include <TObject.h> // for TObject
#include "Rtypes.h"  // for Double_t, Int_t, Bool_t, V11Geometry::Class, etc
 
class TGeoArb8;    // lines 11-11
class TGeoBBox;    // lines 16-16
class TGeoConeSeg; // lines 15-15
class TGeoPcon;    // lines 12-12
class TGeoTube;    // lines 13-13
class TGeoTubeSeg; // lines 14-14
 
namespace o2
{
namespace its
{
 
class V11Geometry : public TObject
{
 
 public:
  V11Geometry() : mDebug(), mDetName("ITS"){};
 
  V11Geometry(Int_t debug, const char* name = "ITS") : mDebug(debug), mDetName(name){};
 
  ~V11Geometry()
    override = default;
 
  void SetDetName(const char* name)
  {
    mDetName = name;
  }
 
  const char* GetDetName() const
  {
    return mDetName;
  }
 
  void setDebug(Int_t level = 5)
  {
    mDebug = level;
  }
 
  void setNoDebug()
  {
    mDebug = 0;
  }
 
  Bool_t getDebug(Int_t level = 1) const
  {
    return mDebug >= level;
  }
 
  // Static functions
 
  Double_t sinD(Double_t deg) const
  {
    return TMath::Sin(deg * TMath::DegToRad());
  }
 
  Double_t cosD(Double_t deg) const
  {
    return TMath::Cos(deg * TMath::DegToRad());
  }
 
  Double_t tanD(Double_t deg) const
  {
    return TMath::Tan(deg * TMath::DegToRad());
  }
 
  void intersectLines(Double_t m, Double_t x0, Double_t y0, Double_t n, Double_t x1, Double_t y1, Double_t& xi,
                      Double_t& yi) const;
 
  static Bool_t intersectCircle(Double_t m, Double_t x0, Double_t y0, Double_t rr, Double_t xc, Double_t yc,
                                Double_t& xi1, Double_t& yi1, Double_t& xi2, Double_t& yi2);
 
  Double_t yFrom2Points(Double_t x0, Double_t y0, Double_t x1, Double_t y1, Double_t x) const;
 
  Double_t xFrom2Points(Double_t x0, Double_t y0, Double_t x1, Double_t y1, Double_t y) const;
 
  Double_t rMaxFrom2Points(const TGeoPcon* p, Int_t i1, Int_t i2, Double_t z) const;
 
  Double_t rMinFrom2Points(const TGeoPcon* p, Int_t i1, Int_t i2, Double_t z) const;
 
  Double_t rFrom2Points(const Double_t* ar, const Double_t* az, Int_t i1, Int_t i2, Double_t z) const;
 
  Double_t zFrom2MinPoints(const TGeoPcon* p, Int_t i1, Int_t i2, Double_t r) const;
 
  Double_t zFrom2MaxPoints(const TGeoPcon* p, Int_t i1, Int_t i2, Double_t r) const;
 
  Double_t zFrom2Points(const Double_t* az, const Double_t* ar, Int_t i1, Int_t i2, Double_t r) const;
 
  Double_t rMaxFromZpCone(const TGeoPcon* p, int ip, Double_t tc, Double_t z, Double_t th = 0.0) const;
 
  // General Cone surface equation R(z). Returns the value R correstponding to the line at
  // angle th, offset by th, and the point p->GetZ/Rmax[ip] at the location z.
  // \param Double_t ar The array of R values
  // \param Double_t az The array of Z values
  // \param Int_t    ip The index in p to get the point location
  // \param Double_t tc The angle of that part of the cone is at
  // \param Double_t  z The value of z to compute R from
  // \param Double_t th The perpendicular distance the parralell line is from the point ip
  Double_t rFromZpCone(const Double_t* ar, const Double_t* az, int ip, Double_t tc, Double_t z,
                       Double_t th = 0.0) const;
 
  Double_t rMinFromZpCone(const TGeoPcon* p, Int_t ip, Double_t tc, Double_t z, Double_t th = 0.0) const;
 
  Double_t zFromRMaxpCone(const TGeoPcon* p, int ip, Double_t tc, Double_t r, Double_t th = 0.0) const;
 
  Double_t zFromRMaxpCone(const Double_t* ar, const Double_t* az, Int_t ip, Double_t tc, Double_t r,
                          Double_t th = 0.0) const;
 
  Double_t zFromRMinpCone(const TGeoPcon* p, int ip, Double_t tc, Double_t r, Double_t th = 0.0) const;
 
  void insidePoint(const TGeoPcon* p, Int_t i1, Int_t i2, Int_t i3, Double_t Cthick, TGeoPcon* q, Int_t j1,
                   Bool_t max) const;
 
  void insidePoint(Double_t x0, Double_t y0, Double_t x1, Double_t y1, Double_t x2, Double_t y2, Double_t c,
                   Double_t& x, Double_t& y) const;
 
  void radiusOfCurvature(Double_t rc, Double_t theta0, Double_t z0, Double_t r0, Double_t theta1, Double_t& z1,
                         Double_t& r1) const;
 
  // Output functions for debugging
 
  void printArb8(const TGeoArb8* a) const;
 
  void printPcon(const TGeoPcon* a) const;
 
  void printTube(const TGeoTube* a) const;
 
  void printTubeSeg(const TGeoTubeSeg* a) const;
 
  void printConeSeg(const TGeoConeSeg* a) const;
 
  void printBBox(const TGeoBBox* a) const;
 
  void drawCrossSection(const TGeoPcon* p, Int_t fillc = 7, Int_t fills = 4050, Int_t linec = 3, Int_t lines = 1,
                        Int_t linew = 4, Int_t markc = 2, Int_t marks = 4, Float_t marksize = 1.0) const;
 
  Bool_t angleOfIntersectionWithLine(Double_t x0, Double_t y0, Double_t x1, Double_t y1, Double_t xc, Double_t yc,
                                     Double_t rc, Double_t& t0, Double_t& t1) const;
 
  /*
      <img src="picts/ITS/V11Geometry_AnglesForRoundedCorners.gif">
     */
  // End_Html
  void anglesForRoundedCorners(Double_t x0, Double_t y0, Double_t r0, Double_t x1, Double_t y1, Double_t r1,
                               Double_t& t0, Double_t& t1) const;
 
  void createDefaultMaterials();
 
  virtual void createMaterials(){};
 
  void makeFigure1(Double_t x0 = 0.0, Double_t y0 = 0.0, Double_t r0 = 2.0, Double_t x1 = -4.0, Double_t y1 = -2.0,
                   Double_t r1 = 1.0);
 
 protected:
  // Units, Convert from k?? to cm,degree,GeV,seconds,
  static const Double_t sMicron;  
  static const Double_t sMm;      
  static const Double_t sCm;      
  static const Double_t sDegree;  
  static const Double_t sRadian;  
  static const Double_t sGCm3;    
  static const Double_t sKgm3;    
  static const Double_t sKgdm3;   
  static const Double_t sCelsius; 
  static const Double_t sPascal;  
  static const Double_t sKPascal; 
  static const Double_t sEV;      
  static const Double_t sKEV;     
  static const Double_t sMEV;     
  static const Double_t sGEV;     
 
 private:
  Double_t angleForRoundedCorners0(Double_t dx, Double_t dy, Double_t sdr) const;
 
  Double_t angleForRoundedCorners1(Double_t dx, Double_t dy, Double_t sdr) const;
 
  Int_t mDebug;                     
  const char* mDetName;             
  ClassDefOverride(V11Geometry, 1); // Base class for ITS v11 geometry
};
} // namespace its
} // namespace o2
 
#endif