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 / 2;  // store half-value
    }

    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);
    }

    // cppcheck-suppress truncLongCastReturn
    double getSigmaAlpha() const {return m_sigmaAlpha * 2;}  // multiply by 2 since half-value stored

    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