TOPGeoBarSegment.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
 
#pragma once
 
#include <top/dbobjects/TOPGeoBase.h>
#include <geometry/dbobjects/GeoOpticalSurface.h>
#include <string>
#include <vector>
#include <utility>
 
 
namespace Belle2 {
  class TOPGeoBarSegment: public TOPGeoBase {
  public:
 
    TOPGeoBarSegment()
    {}
 
    TOPGeoBarSegment(double width, double thickness, double length,
                     const std::string& material,
                     const std::string& name = "TOPBarSegment"): TOPGeoBase(name),
      m_width(width), m_thickness(thickness), m_length(length), m_material(material)
    {}
 
    virtual void setGlue(double thickness, const std::string& material)
    {
      m_glueThickness = thickness;
      m_glueMaterial = material;
    }
 
    virtual void setGlueDelamination(double fraction, double angle,
                                     const std::string& material)
    {
      m_brokenFraction = fraction;
      m_brokenAngle = angle;
      m_brokenGlueMaterial = material;
    }
 
    void setSurface(const GeoOpticalSurface& surface, double sigmaAlpha)
    {
      m_surface = surface;
      m_sigmaAlpha = sigmaAlpha;
    }
 
    void setVendorData(const std::string& vendor, const std::string& serialNumber)
    {
      m_vendor = vendor;
      m_serialNumber = serialNumber;
    }
 
    double getWidth() const {return m_width / s_unit;}
 
    double getThickness() const {return m_thickness / s_unit;}
 
    double getLength() const {return m_length / s_unit;}
 
    double getFullLength() const {return (m_length + m_glueThickness) / s_unit;}
 
    const std::string& getMaterial() const {return m_material;}
 
    virtual double getGlueThickness() const {return m_glueThickness / s_unit;}
 
    const std::string& getGlueMaterial() const {return m_glueMaterial;}
 
    double getBrokenGlueFraction() const {return m_brokenFraction;}
 
    double getBrokenGlueAngle() const {return m_brokenAngle;}
 
    const std::string& getBrokenGlueMaterial() const {return m_brokenGlueMaterial;}
 
    std::vector<std::pair<double, double> > getBrokenGlueContour() const;
 
    const GeoOpticalSurface& getSurface() const {return m_surface;}
 
    double getSurfaceReflectivity(double energy) const
    {
      return TOPGeoBase::getReflectivity(m_surface, energy);
    }
 
    double getSigmaAlpha() const {return m_sigmaAlpha;}
 
    const std::string& getVendor() const {return m_vendor;}
 
    const std::string& getSerialNumber() const {return m_serialNumber;}
 
    virtual bool isConsistent() const override;
 
 
    virtual void print(const std::string& title = "Bar segment geometry parameters") const override;
 
  protected:
 
    void constructContour(double A, double B, double fraction, double angle,
                          std::vector<std::pair<double, double> >& contour) const;
 
    float m_width = 0; 
    float m_thickness = 0; 
    float m_length = 0; 
    std::string m_material; 
    float m_glueThickness = 0; 
    std::string m_glueMaterial; 
    GeoOpticalSurface m_surface; 
    float m_sigmaAlpha = 0; 
    float m_brokenFraction = 0; 
    float m_brokenAngle = 0; 
    std::string m_brokenGlueMaterial; 
    std::string m_vendor; 
    std::string m_serialNumber; 
    ClassDefOverride(TOPGeoBarSegment, 3); 
  };
 
} // end namespace Belle2