MdcGeomSvc-00-01-42/MdcGeomSvc/MdcGeomSvc.h

Go to the documentation of this file.

// header file for a class called "MdcGeomSvc"
#ifndef MDCGEOMSVC_H
#define MDCGEOMSVC_H
 
#include <CLHEP/Units/PhysicalConstants.h>  
#include <CLHEP/Geometry/Point3D.h> 
#include <vector>
#include <iostream>
#include "GaudiKernel/Algorithm.h"
 
#include "MdcGeomSvc/MdcGeoWire.h"
#include "MdcGeomSvc/MdcGeoLayer.h"
#include "MdcGeomSvc/MdcGeoGeneral.h"
#include "MdcGeomSvc/MdcGeoSuper.h"
#include "MdcGeomSvc/MdcGeoMisc.h"
#include "MdcGeomSvc/MdcGeoEnd.h"  
#include "GaudiKernel/IIncidentListener.h"
#include "GaudiKernel/Service.h"
#include "GaudiKernel/IInterface.h"
#include "MdcGeomSvc/IMdcGeomSvc.h"
#include "GaudiKernel/IDataProviderSvc.h"
#include "GaudiKernel/IIncidentListener.h"
#include "CLHEP/Vector/TwoVector.h"
 
//class MdcGeomSvc : public Service, virtual public IMdcGeomSvc,virtual public IIncidentListener
class MdcGeomSvc : public Service, virtual public IMdcGeomSvc, virtual public IIncidentListener
{
 public:
   // Constructors and destructors
   MdcGeomSvc( const std::string& name, ISvcLocator* svcloc );
   ~MdcGeomSvc();
 
   virtual StatusCode queryInterface(const InterfaceID& riid, void** ppvUnknown);
   virtual StatusCode initialize ( );
   virtual StatusCode finalize ( );
     void handle(const Incident& inc);
 
   // Access by others
   const MdcGeoWire * const Wire(unsigned id);
   const MdcGeoWire * const Wire(unsigned lyrid,unsigned wirid);  
   const MdcGeoLayer * const Layer(unsigned id);
   const MdcGeoSuper * const SuperLayer(unsigned id);
   const MdcGeoGeneral * const GeneralLayer(unsigned id);
   const MdcGeoMisc * const Misc(void);    
   const MdcGeoEnd * const End(unsigned id);    
   const int getWireSize();
   const int getLayerSize();
   const int getSuperLayerSize();
   const int getGeneralLayerSize(); 
   const int getSegmentNo();    
   void Dump();
   static bool getSagFlag(void); 
 
   // CrossPoint relevant
   /**
    * @brief forms a line in x-y space, ax + by + c = 0 
    */
   struct AdjCandiTgtLine2D {
       const static double INFTY;
       double a, b, c,  // ax + by + c = 0
          x0,y0,z0, //a point coordinate on the line, center of two tangant points
          angle2D,  //angle of the line and the 'plane' from two wires at z=z0 plane given
          angle3D;  //angle of the line and the 'plane' from two wires
       AdjCandiTgtLine2D ():a(INFTY),b(INFTY),c(INFTY),x0(INFTY),y0(INFTY),z0(INFTY),angle2D(INFTY),angle3D(INFTY){};// float point NaN or Infinity
       /**
        * @brief get direction vector of a line.  garanteed to be outwards
        */
       CLHEP::Hep2Vector direction(){ 
           CLHEP::Hep2Vector ret(b,-a);
           if (ret.dot(location())<0)
               ret=-ret;
           return ret;
       }
       /**
        * @brief returns the cosine of angle of the line with R at its location (x0,y0). note that location is the mid point of two wires at selected plane.
        * 
        * @return double 
        */
       double angleWithR(){
           CLHEP::Hep2Vector directionR(x0,y0);
           CLHEP::Hep2Vector dirRRot90(-y0,x0);
           CLHEP::Hep2Vector dirL=direction();
           int sgn=dirRRot90.dot(dirL)>0; // bad-optimized code for telling rotate L-R.
           return sgn*directionR.angle(dirL);
       }
 
       CLHEP::Hep2Vector location(){
           CLHEP::Hep2Vector ret(x0,y0);
           return ret;
       }
       HepPoint3D location3(){
           HepPoint3D ret(x0,y0,z0);
           return ret;
       }
       bool isLocationValid(){
           return x0<99999 && x0>-99999; 
       }
   };
 
   vector<HepPoint3D> getAdjacentIntersectionShape(const MdcGeoWire* wireA,const MdcGeoWire* wireB) ;
   vector<double> getAdjacentIntersectionBounds(const MdcGeoWire* wireA,const MdcGeoWire* wireB);
   HepPoint3D getAdjacentIntersectionPointFast(const MdcGeoWire* wireA,const MdcGeoWire* wireB);//dont use! 
   HepPoint3D getAdjacentIntersectionPointSlower(const MdcGeoWire* wireA,const MdcGeoWire* wireB);
   vector<AdjCandiTgtLine2D> getAdjacentCandidateWithDriftCircle(const MdcGeoWire* wireA,const MdcGeoWire* wireB, double driftA, double driftB);
   vector<AdjCandiTgtLine2D> getAdjacentCandidateWithDriftCircle(const MdcGeoWire* wireA,const MdcGeoWire* wireB, double driftA, double driftB, double expectedZ);
 
   void testAdjacentIntersection();
 
 
 
 
 public:
   static bool m_doSag;
   static bool m_readAlignParDataBase;
   static bool m_nomcalignment;
 
   // private method
 private:
   void clean();  // clean the fWires,fGenerals,fLayers,fSupers,,fEnd,etc... 
   void ReadFilePar();        //get geometry data from file SimUtil/dat/Mdc.txt
   void ReadTensionDataBase(std::vector<double> & wireTensionVec);
   void ReadWirePosDataBase(std::vector<vector<double> > & wirePosVec);
   void ReadAliParDataBase(vector<double>& Sx, vector<double>& Sy, vector<double>& Sz,
           vector<double>& Rx, vector<double>& Ry, vector<double>& Rz);
   void Fill();  //get geometry data from Database
   void AddWireNeighbors();
   void RenewWireNeighbors(double range);
   const int getAlignParIndexEast(int lyr) const;
   const int getAlignParIndexWest(int lyr) const;
 
   const double CalculatePhi(double phi) const;
 
 
 
   // private data members
 private:
   double m_wholeShiftX; 
   double m_wholeShiftY;
   double m_wholeShiftZ;
   double m_wholeRotatX;
   double m_wholeRotatY;
   double m_wholeRotatZ; 
   int    m_mindex;
   std::string m_alignFilePath;
   std::string m_wirePosFilePath;
   std::string m_wireTensionFilePath;
   bool m_updataalign;
   bool m_useCgem;
 
   double m_NeighborRange;//
 
   bool m_conf_testAdjacentIntersection;
 
   vector <MdcGeoWire*> fWires;  //wire list (6860)
   vector <MdcGeoGeneral> fGenerals; //general layer list (90)
   vector <MdcGeoLayer*> fLayers;    //sense wire layer list (43)
   vector <MdcGeoSuper*> fSupers;    //super layer list (11)
   MdcGeoMisc fMisc;             // parms. of MDC misc
   vector <MdcGeoEnd*> fEnd;       // tube and endplate list (67)
   IDataProviderSvc* m_pCalibDataSvc;
 
};
 
 
#endif