CgemMidDriftPlane.cxx

Go to the documentation of this file.

#include "CgemGeomSvc/CgemMidDriftPlane.h"
 
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/IMessageSvc.h"
#include "GaudiKernel/StatusCode.h"
#include "GaudiKernel/ISvcLocator.h"
#include "GaudiKernel/Bootstrap.h"
#include "GaudiKernel/SmartDataPtr.h"
#include "GaudiKernel/IDataProviderSvc.h"
#include "GaudiKernel/PropertyMgr.h"
 
#include <iostream>
#include <fstream>
#include <iomanip>
 
using namespace std;
 
CgemMidDriftPlane::CgemMidDriftPlane(){
}
 
CgemMidDriftPlane::~CgemMidDriftPlane(){
}
 
bool CgemMidDriftPlane::getPointIdealGeom(int layer, StraightLine pLine,
                      HepPoint3D& posUp, HepPoint3D& posDown,
                      double phiVUp[], double phiVDown[]){
 
     posDown = pLine.xAtR(m_rad[layer], 1);
     xyzToPhiVIdealGeom(layer, posDown, phiVDown);
 
     posUp = pLine.xAtR(m_rad[layer], -1);
     xyzToPhiVIdealGeom(layer, posUp, phiVUp);
 
     return true;
}
 
bool CgemMidDriftPlane::getPointIdealGeom(int layer, StraightLine* pLine,
                      HepPoint3D& posUp, HepPoint3D& posDown,
                      double phiVUp[], double phiVDown[]){
     return getPointIdealGeom(layer, *pLine, posUp, posDown, phiVUp, phiVDown);
}
 
bool CgemMidDriftPlane::getPointAligned(int layer, StraightLine pLine,
                    HepPoint3D& posUp, HepPoint3D& posDown,
                    double phiVUp[], double phiVDown[]){
     double trk[4];
     trk[0] = pLine.dr();
     trk[1] = pLine.phi0();
     trk[2] = pLine.dz();
     trk[3] = pLine.tanl();
 
     double trkCnv[4];
     StraightLine lineCnv = m_align->StraightLineConversion(layer, pLine);
 
     HepPoint3D posDownLocal = lineCnv.xAtR(m_rad[layer], 1);
     xyzToPhiVIdealGeom(layer, posDownLocal, phiVDown);
     posDown = m_align->point_invTransform(layer, posDownLocal);
 
     HepPoint3D posUpLocal = lineCnv.xAtR(m_rad[layer], -1);
     xyzToPhiVIdealGeom(layer, posUpLocal, phiVUp);
     posUp = m_align->point_invTransform(layer, posUpLocal);
 
     return true;
}
 
bool CgemMidDriftPlane::getPointAligned(int layer, StraightLine* pLine,
                    HepPoint3D& posUp, HepPoint3D& posDown,
                    double phiVUp[], double phiVDown[]){
     return getPointAligned(layer, *pLine, posUp, posDown, phiVUp, phiVDown);
}
 
bool CgemMidDriftPlane::xyzToPhiVIdealGeom(int layer, HepPoint3D pos, double phiV[]){
     double phi=pos.phi();
     double z  =pos.z();
     phiV[0]=phi;
     phiV[1]=99999.;
     int nSheet=2; if(layer==0) nSheet=1;
     for(int iSheet=0; iSheet<nSheet; iSheet++)
     {
      if(m_readoutPlane[layer][iSheet]->OnThePlane(phi,z))
      {
           phiV[1]=m_readoutPlane[layer][iSheet]->getVFromPhiZ(phi,z);
           break;
      }
     }
 
     return true;
}