CgemGeomSvc.h

Go to the documentation of this file.

// length unit is mm
 
#ifndef CGEMGEOMSVC_H
#define CGEMGEOMSVC_H 
 
#include <CLHEP/Units/PhysicalConstants.h>
#include <CLHEP/Geometry/Point3D.h>
#include <vector>
#include <iostream>
#include <string>
 
#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/Property.h"
#include "GaudiKernel/MsgStream.h"
 
#include "GaudiKernel/IIncidentListener.h"
#include "GaudiKernel/Service.h"
#include "GaudiKernel/IInterface.h"
#include "CgemGeomSvc/ICgemGeomSvc.h"
#include "GaudiKernel/IDataProviderSvc.h"
#include "GaudiKernel/IIncidentListener.h"
 
#include "CLHEP/Vector/ThreeVector.h"
 
class CgemGeomSvc : public Service, virtual public ICgemGeomSvc, virtual public IIncidentListener{
 
    public:
        CgemGeomSvc(const std::string& name, ISvcLocator* pSvcLocator);
        ~CgemGeomSvc(); 
 
        virtual StatusCode queryInterface(const InterfaceID& riid, void** ppvUnknown);
        virtual StatusCode initialize();
        virtual StatusCode finalize();
        void handle(const Incident& inc);
        bool initGeom();
 
        void readAlignPar();
 
        /* Information of Cgem */
        double getInnerROfCgem() const {return m_R_i_Cgem;}
        double getOuterROfCgem() const {return m_R_o_Cgem;}
        double getLengthOfCgem() const {return m_L_Cgem;}
 
        /* Information of CgemLayer */
        int getNumberOfCgemLayer()     const {return m_N_CgemLayer ;}
        int getNumberOfCgemFoil()       const {return m_N_GemFoil   ;}
        double getThicknessOfCgemLayer(int i) const {return m_T_CgemLayer[i] ;}
        CgemGeoLayer* getCgemLayer(int i) const {return m_CgemLayer[i];}
 
        double getThicknessOfGapD(int i)        const {return m_T_Gap_D[i]  ;}
        double getThicknessOfCathode(int i)          const {return m_T_Cathode[i]   ;}
        double getThicknessOfAnode(int i)          const {return m_T_Anode[i]   ;}
 
        CgemGeoReadoutPlane* getReadoutPlane(int iLayer, int iSheet) const;
        /* information about strip description */
        bool isXStripDescriptionOn()        const {return m_B_strip_x_on; } 
        bool isVStripDescriptionOn()        const {return m_B_strip_v_on; } 
 
        void setDensityOfCable(double rho) {m_Rho_Cable = rho;}
        double getDensityOfCable() const {return m_Rho_Cable;}
        string getMaterialOfCgemFoil(int i) const {return nmg_Materials[i+1];}
        string getMaterialOfCathode(int i) const {return nmc_Materials[i+1];}
        string getMaterialOfAnode(int i) const {return nma_Materials[i+1];}
 
        /* Information of CGEM/MDC Separator */
        CgemGeoSeparator* getCgemSeparator() const {return m_CgemSeparator;}
 
        /* get Alignment pointer */
        CgemGeoAlign* getAlignPtr() const{ return m_align; }
 
        /* get the pointer of mid-drift plane */
        CgemMidDriftPlane* getMidDriftPtr() const{ return m_midDrift; }
 
        bool isPassive() const { return m_passive; }
        bool isEffDensityHoles() const {return m_effDensityHoles;}
        bool isEffDensityStrips() const {return m_effDensityStrips;}
 
        static const int NLAYER_MAX = 3;
        static const int NFOIL_MAX = 3;
        static const int NSHEET_MAX = 2;
        static const int N_MATERIAL_MAX = 30;
        static const int NMATERIALS_CATHODE = 14;
        static const int NMATERIALS_ANODE = 20;
        static const int NMATERIALS_FOIL = 3;
        int getNMaterialsCathode() const {return NMATERIALS_CATHODE;}
        int getNMaterialsAnode() const {return NMATERIALS_ANODE;}
        int getNMaterialsFoil() const {return NMATERIALS_FOIL;}
 
    private:
        int m_iFile;
        bool m_debugCout;
        bool m_testReadoutPlane;
        /*      std::string m_geomFile; */
        std::string m_alignFile;
        CgemGeoAlign* m_align;
        CgemMidDriftPlane* m_midDrift;
 
        /* passive elements on/off */
        bool m_passive;
 
        /* Information of Cgem */
        double m_R_i_Cgem, m_R_o_Cgem, m_L_Cgem;
        string nmc_Materials[NMATERIALS_CATHODE+1];
        string nma_Materials[NMATERIALS_ANODE+1];
        string nmg_Materials[NMATERIALS_FOIL+1];
 
        double  m_T_Gap_D[NLAYER_MAX]  ;
        double m_T_Cathode[NLAYER_MAX]  ;
        double m_T_Anode[NLAYER_MAX] ;
 
        /* Information of CgemLayer */
        int       m_N_CgemLayer ;
        int       m_N_GemFoil   ;
        double    m_T_CgemLayer[NLAYER_MAX] ;
        CgemGeoLayer* m_CgemLayer[NLAYER_MAX];
        CgemGeoFoil* m_foil[NLAYER_MAX][NFOIL_MAX];
        CgemGeoSeparator *m_CgemSeparator;
 
 
        /* Orientation 
           0 = straight
           1 = reversed */
        bool m_Orientation;
        bool m_effDensityHoles;
        bool m_effDensityStrips;
 
 
        /* Do we want to use the strip description in the anode? */
        bool m_B_strip_x_on;
        bool m_B_strip_v_on;
        CgemGeoReadoutPlane* m_ReadoutPlane[NLAYER_MAX][NSHEET_MAX];
        void TestReadoutPlane();
        /** density of cables
          taking into account fill factor and weight/volume of HV, LV, signal cables
         **/
        double m_Rho_Cable;
 
};
 
inline CgemGeoReadoutPlane* CgemGeomSvc::getReadoutPlane(int iLayer, int iSheet) const {
    if( (iLayer<NLAYER_MAX) && (iSheet<NSHEET_MAX) ) return m_ReadoutPlane[iLayer][iSheet];
    else return 0;
}
 
 
#endif  /* CGEMGEOMSVC_H */