MdcLayer.cxx

Go to the documentation of this file.

#include "MdcGeom/MdcLayer.h"
#include "MdcGeom/MdcDetector.h"
#include <iomanip>
#include <iostream>
using namespace std;
 
    // Constructor.
 
MdcLayer::MdcLayer(int lay, int nwir, MdcSWire** wires, const MdcDetector &theDet)
  : _mdcDet(theDet), _exist(1),  _layer(lay), _nwires(nwir),
    _wires(wires)
{
  if (_nwires == 0) _exist = 0;
  // set pointer to layer for wires
  for (int wire=0; wire<_nwires; wire++) {
    getWire(wire)->setLayerPtr(this);
//    _cells.push_back(buildCell(wire));
  }
 
  // get nominal values from wire # 0
  _rend = wires[0]->rEnd();
  _rmid = wires[0]->rMid();
  _stdip = _rend - _rmid;
  // in the local (DC) system frame
  _zend = wires[0]->zEndDC();
  _stereo = wires[0]->stereo();
  // half twist angle
  _delphi = wires[0]->dPhiz();
  // stereo view
  if((_layer>=0 && _layer<=7)||(_layer>=20&&_layer<=35)){
    _view = (int)(_stereo/(fabs(_stereo)));
  } else {
    _view = 0;
  }
  _phiOffset = wires[0]->phiE();
}
 
 
 
void 
MdcLayer::print( std::ostream& o ) const
{
  o<< setw(3)<<layNum()<<
    setw(3)<<view()<<
    setw(4)<<nWires()<<
    setw(3)<<superLayer()<<
    setw(3)<<subLayer()<<
    setw(10)<<stDip()<<
    setw(10)<<zEnd()<<
    setw(10)<<rMid()<<
    setw(10)<<rEnd()<<
    setw(10)<<rIn()<<
    setw(10)<<rOut()<<
    setw(10)<<stereo()<<
    setw(10)<<dPhiz()<<
    setw(7)<<zLength()<<
    setw(10)<<sag()<<
    setw(10)<<dPhi()<<
    setw(10)<<phiOffset()<<
    setw(10)<<phiEPOffset()<<
    setw(10)<<dPhizDC(0.)<<
    setw(10)<<radiusDC(0.)<<
    setw(10)<<phiWireDC(0,0.)<<
    setw(10)<<phiWire(0)<<
    setw(10)<<xWire(0)<<
    setw(10)<<yWire(0)<<
    setw(10)<<cellWidth()<<
    setw(10)<<cellWidth(0.)<<
    std::endl;
}
 
 
double 
MdcLayer::phiWireDC(int cell, double z) const 
{
  if ( cell >= 0 && cell < nWires() ) {
    return getWire(cell)->phiDC(z);
  } else {
    std::cout << "phiWireDC: wrong cell # " << cell << "\n"
          << "  number of cells in this layer is " << nWires() 
          << std::endl;
    return 0.;
  }
}
 
double 
MdcLayer::phiWire(int cell) const 
{
  // in the middle of the chamber
  if ( cell >= 0 && cell < nWires() ) {
    return getWire(cell)->phi();
  } else {
    std::cout  << "phiWire: wrong cell # " << cell << "\n"
          << "  number of cells in this layer is " << nWires() 
          << std::endl;
    return 0.;
  }
}
 
double 
MdcLayer::xWire(int cell) const 
{
  // in the middle of the chamber
  if ( cell >= 0 && cell < nWires() ) {
    return getWire(cell)->xMid();
  } else {
    std::cout << "xWire: wrong cell # " << cell << "\n"
          << "  number of cells in this layer is " << nWires() 
          << std::endl;
    return 0.;
  }
}
 
double 
MdcLayer::yWire(int cell) const 
{
  // in the middle of the chamber 
  if ( cell >= 0 && cell < nWires() ) {
    return getWire(cell)->yMid(); 
  } else {
    std::cout << "YWire: wrong cell # " << cell << "\n"
          << "  number of cells in this layer is " << nWires() 
          << std::endl;
    return 0.;
  }
}
 
 
 
 
const Trajectory*
MdcLayer::makeHitTrajInGlobalCoords(int wire, double /*z*/) const
{
  if ( wire >= 0 && wire < nWires() ) {
    return getWire(wire)->getTraj();
  } else {
    std::cout << "Error In makeHitTrajInGlobalCoords: wrong cell # "<< wire << "\n"
                 <<"  number of cells in this layer is " << nWires() << std::endl;
    return getWire(0)->getTraj();
  }
}