MuonGeometryWriter.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef JIVEXML_MUONGEOMETRYWRITER
#define JIVEXML_MUONGEOMETRYWRITER
 
#include "JiveXML/IGeometryWriter.h"
#include "AthenaBaseComps/AthAlgTool.h"
#include "CLHEP/Geometry/Point3D.h"
#include "GeoPrimitives/GeoPrimitives.h"
#include "GeoModelKernel/GeoVolumeCursor.h"
#include "GeoModelKernel/GeoSimplePolygonBrep.h"
 
#include <string>
#include <fstream>
 
namespace MuonGM {
  class MuonDetectorManager;
  class MuonStation;
}
 
namespace JiveXML
{
  class MuonGeometryWriter : virtual public IGeometryWriter, public AthAlgTool {
 
  public:

    MuonGeometryWriter(const std::string& t,const std::string& n,const IInterface* p):
      AthAlgTool(t,n,p){declareInterface<IGeometryWriter>(this);};

    virtual StatusCode writeGeometry();
 
  private:

    void writeHeader(std::ofstream &out) const ;

    void writeStations(std::ofstream &out) const ;

    void processNSW(std::ofstream &out);

    void writeFooter(std::ofstream &out) const ;

    void writeATrd(std::ofstream &out,
                   const std::string& stationTech,
                   const std::string& stationName,
                   double zi, double zo,
                   double ri, double ro,
                   double wi, double wo, int eta,
                   const std::string& phiString,
                   double dphi, double shift, double alpha) const ;
 

    HepGeom::Point3D<double> getPosition(const MuonGM::MuonStation *station, int maxPhi) const ;

    double getDeltaPhi(const HepGeom::Point3D<double> &pos, int maxPhi) const ;

    double getShift(const HepGeom::Point3D<double> &pos, double dphi) const ;

    double getAlpha(const HepGeom::Transform3D &trans) const ;

    HepGeom::Point3D<double> getPositionNSW(Amg::Vector3D pos, int maxPhi) const ;

    int getPhiIndex(double phi, int maxPhi) const ;

    void readNSWMMPars(const GeoVolumeCursor *pv, int maxPhi, std::string& chamberName, HepGeom::Point3D<double>& pos_rot,
          double& zi, double& zo, double& ri, double& ro, double& wi, double& wo, double& dphi, double& shift, int& phiIndex) const ;

    void readNSWSTGCPars(const GeoVolumeCursor *pv, int maxPhi,
          std::string& chamberName, HepGeom::Point3D<double>& pos_rot, const GeoSimplePolygonBrep*& theBrep,
          int& nvtx, double& dz, double& dphi, double& shift, int& phiIndex) const ;

    void readBrepAsATrd(const GeoSimplePolygonBrep* theBrep, double rho, const int* vtx, double& ri, double& ro, double& wi, double& wo) const ;

    bool equalAngle(double a, double b) const ;

    bool equalLength(double a, double b) const ;

    const MuonGM::MuonDetectorManager* m_muon_manager{nullptr} ;

    static const double m_smallAngle, m_smallDistance ;
  };
}
#endif