HcalTBEventPosition.cc

Go to the documentation of this file.

#include "TBDataFormats/HcalTBObjects/interface/HcalTBEventPosition.h"
#include <cstdio>

using namespace std;

  HcalTBEventPosition::HcalTBEventPosition() :
    hfTableX_(-1e24),
    hfTableY_(-1e24),
    hfTableV_(-1e24),
    hbheTableEta_(-1000),
    hbheTablePhi_(-1000),
    ax_(),
    ay_(),
    bx_(),
    by_(),
    cx_(),
    cy_(),
    dx_(),
    dy_(),
    ex_(),
    ey_(),
    fx_(),
    fy_(),
    gx_(),
    gy_(),
    hx_(),
    hy_() {
  }

  void HcalTBEventPosition::getChamberHits ( char chamberch,
                         vector<double>& xvec,
                         vector<double>& yvec ) const {
    switch (chamberch) {
    case 'A': xvec = ax_; yvec = ay_; break;
    case 'B': xvec = bx_; yvec = by_; break;
    case 'C': xvec = cx_; yvec = cy_; break;
    case 'D': xvec = dx_; yvec = dy_; break;
    case 'E': xvec = ex_; yvec = ey_; break;
    case 'F': xvec = fx_; yvec = fy_; break;
    case 'G': xvec = gx_; yvec = gy_; break;
    case 'H': xvec = hx_; yvec = hy_; break;
    default:
      cerr << "Unrecognized chamber character " << chamberch << endl;
      return;
    }
  }

  void HcalTBEventPosition::setChamberHits ( char chamberch,
                         const vector<double>& xvec,
                         const vector<double>& yvec ) {
    switch (chamberch) {
    case 'A': ax_ = xvec; ay_ = yvec; break;
    case 'B': bx_ = xvec; by_ = yvec; break;
    case 'C': cx_ = xvec; cy_ = yvec; break;
    case 'D': dx_ = xvec; dy_ = yvec; break;
    case 'E': ex_ = xvec; ey_ = yvec; break;
    case 'F': fx_ = xvec; fy_ = yvec; break;
    case 'G': gx_ = xvec; gy_ = yvec; break;
    case 'H': hx_ = xvec; hy_ = yvec; break;
    default:
      cerr << "Unrecognized chamber character " << chamberch << endl;
      return;
    }
  }

  void HcalTBEventPosition::setHFtableCoords( double x, double y, double v ) {
    hfTableX_ = x;
    hfTableY_ = y;
    hfTableV_ = v;
  }

  void HcalTBEventPosition::setHBHEtableCoords( double eta, double phi ) {
    hbheTableEta_ = eta;
    hbheTablePhi_ = phi;
  }

  ostream& operator<<(ostream& s, const HcalTBEventPosition& htbep) {
    s << "HF Table (X,Y,V) = (" << htbep.hfTableX() << "," << htbep.hfTableY() << "," << htbep.hfTableV() << ")" << std::endl;
    s << "HB/HE Table (eta,phi) = (" << htbep.hbheTableEta() << "," << htbep.hbheTablePhi() << ")" << std::endl;

    vector<double> xvec, yvec;
    vector<double>::const_iterator j;

    htbep.getChamberHits('A', xvec, yvec);

    s << "WC Ax: ";
    for (j=xvec.begin(); j!=xvec.end(); j++) {
      if (j!=xvec.begin()) s << ", ";
      s << *j;
    }
    s << endl;

    s << "WC Ay: ";
    for (j=yvec.begin(); j!=yvec.end(); j++) {
      if (j!=yvec.begin()) s << ", ";
      s << *j;
    }
    s << endl;

    htbep.getChamberHits('B', xvec, yvec);

    s << "WC Bx: ";
    for (j=xvec.begin(); j!=xvec.end(); j++) {
      if (j!=xvec.begin()) s << ", ";
      s << *j;
    }
    s << endl;

    s << "WC By: ";
    for (j=yvec.begin(); j!=yvec.end(); j++) {
      if (j!=yvec.begin()) s << ", ";
      s << *j;
    }
    s << endl;

    htbep.getChamberHits('C', xvec, yvec);

    s << "WC Cx: ";
    for (j=xvec.begin(); j!=xvec.end(); j++) {
      if (j!=xvec.begin()) s << ", ";
      s << *j;
    }
    s << endl;

    s << "WC Cy: ";
    for (j=yvec.begin(); j!=yvec.end(); j++) {
      if (j!=yvec.begin()) s << ", ";
      s << *j;
    }
    s << endl;

    htbep.getChamberHits('D', xvec, yvec);

    s << "WC Dx: ";
    for (j=xvec.begin(); j!=xvec.end(); j++) {
      if (j!=xvec.begin()) s << ", ";
      s << *j;
    }
    s << endl;

    s << "WC Dy: ";
    for (j=yvec.begin(); j!=yvec.end(); j++) {
      if (j!=yvec.begin()) s << ", ";
      s << *j;
    }
    s << endl;

    htbep.getChamberHits('E', xvec, yvec);

    s << "WC Ex: ";
    for (j=xvec.begin(); j!=xvec.end(); j++) {
      if (j!=xvec.begin()) s << ", ";
      s << *j;
    }
    s << endl;

    s << "WC Ey: ";
    for (j=yvec.begin(); j!=yvec.end(); j++) {
      if (j!=yvec.begin()) s << ", ";
      s << *j;
    }
    s << endl;

    htbep.getChamberHits('F', xvec, yvec);

    s << "WC Fx: ";
    for (j=xvec.begin(); j!=xvec.end(); j++) {
      if (j!=xvec.begin()) s << ", ";
      s << *j;
    }
    s << endl;

    s << "WC Fy: ";
    for (j=yvec.begin(); j!=yvec.end(); j++) {
      if (j!=yvec.begin()) s << ", ";
      s << *j;
    }
    s << endl;

    htbep.getChamberHits('G', xvec, yvec);

    s << "WC Gx: ";
    for (j=xvec.begin(); j!=xvec.end(); j++) {
      if (j!=xvec.begin()) s << ", ";
      s << *j;
    }
    s << endl;

    s << "WC Gy: ";
    for (j=yvec.begin(); j!=yvec.end(); j++) {
      if (j!=yvec.begin()) s << ", ";
      s << *j;
    }
    s << endl;

    htbep.getChamberHits('H', xvec, yvec);

    s << "WC Hx: ";
    for (j=xvec.begin(); j!=xvec.end(); j++) {
      if (j!=xvec.begin()) s << ", ";
      s << *j;
    }
    s << endl;

    s << "WC Hy: ";
    for (j=yvec.begin(); j!=yvec.end(); j++) {
      if (j!=yvec.begin()) s << ", ";
      s << *j;
    }
    s << endl;

    return s;
  }