MdcSWire.cxx

Go to the documentation of this file.

#include "MdcGeom/MdcSWire.h"
#include "MdcGeom/BesAngle.h"
#include "MdcGeom/MdcLayer.h"
#include <iomanip>
#include <iostream>
using namespace std;
 
 
MdcSWire::MdcSWire(HepPoint3D &eastP, HepPoint3D &westP, double sag, int id, int cell)
: _traj(sag,eastP,westP), _east(eastP), _west(westP), _sag(sag), _id(id), _cell(cell)
{
  _rend = sqrt(xEast()*xEast() + yEast()*yEast());
  _phiend    = _east.phi();
  BesAngle westPphi(westP.phi());
  BesAngle eastPphi(eastP.phi());
  _twist = ( westPphi - eastPphi)* 0.5;
}
 
Hep3Vector
MdcSWire::yAxis( double z )
{
  Hep3Vector vec( xWireDC(z), yWireDC(z), z);
  return vec.unit();
}
 
void MdcSWire::print(std::ostream &o) const 
{
  o<< setw(6)<<Id()<<
    setw(10)<<xEast()<<
    setw(10)<<yEast()<<
    setw(10)<<zEast()<<
    setw(10)<<xWest()<<
    setw(10)<<yWest()<<
    setw(10)<<zWest()<<
    setw(10)<<xMid()<<
    setw(10)<<yMid()<<
    setw(10)<<rEnd()<<
    setw(10)<<getSag()<<
    setw(10)<<cell()<<
    setw(10)<<zLength()<<
    setw(10)<<dPhiz()<<
    setw(10)<<zEndDC()<<
    setw(10)<<rMid()<<
    setw(10)<<phiE()<<
    setw(10)<<phi()<<
    setw(10)<<stereo()<<
    setw(10)<<stDip()<<
    setw(10)<<radiusDC(0.)<<
    setw(10)<<phiDC(0.)<<
    setw(10)<<xWireDC(0.)<<
    setw(10)<<yWireDC(0.)<<
    setw(10)<<dPhizDC(0.)<<
    setw(10)<<getTraj()->sag()<<//
    setw(10)<<getTraj()->position(0.)<<
    setw(10)<<getTraj()->direction(0.)<<
    setw(10)<<getTraj()->curvature()<<//
    setw(10)<<getTraj()->rawDirection()<<
    std::endl;
}
 
void MdcSWire::printInfo(std::ostream &o) const 
{
  o << "east " << getEastPoint()->x() << " " << getEastPoint()->y() << " "
    << getEastPoint()->z() <<"\n"
    << "west " << getWestPoint()->x() <<" " << getWestPoint()->y()
    << " " << getWestPoint()->z() << "\n"
    << "sagitta: " << getSag()
    << " z length: " << zLength()
    << " twist: " << dPhiz()
    << " stereo: " << stereo() <<"\n"
    << "x-y mid: " << xMid() << " - " << yMid() << std::endl;
 
}
 
std::ostream&  operator << (std::ostream& o, const MdcSWire& w) 
{ 
  w.print(o);
  return o;
}