Line Code
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51
/** \file
 *
 *  $Date: 2008/03/26 21:59:18 $
 *  $Revision: 1.10 $
 *  \author Andre Sznajder - UERJ(Brazil)
 */

#include "Alignment/MuonAlignment/interface/AlignableCSCChamber.h"

AlignableCSCChamber::AlignableCSCChamber(const GeomDet* geomDet) : AlignableDet(geomDet) {
  theStructureType = align::AlignableCSCChamber;
  // DO NOT let the chamber position become an average of the layers
  // FIXME: is this redundant?
  theSurface = geomDet->surface();
  compConstraintType_ = Alignable::CompConstraintType::NONE;
}

void AlignableCSCChamber::update(const GeomDet* geomDet) {
  AlignableDet::update(geomDet);
  // DO NOT let the chamber position become an average of the layers
  // FIXME: is this redundant?
  theSurface = geomDet->surface();
}

/// Printout the DetUnits in the CSC chamber
std::ostream& operator<<(std::ostream& os, const AlignableCSCChamber& r) {
  const auto& theDets = r.components();

  os << "    This CSCChamber contains " << theDets.size() << " units" << std::endl;
  os << "    position = " << r.globalPosition() << std::endl;
  os << "    (phi, r, z)= (" << r.globalPosition().phi() << "," << r.globalPosition().perp() << ","
     << r.globalPosition().z();
  os << "), orientation:" << std::endl << r.globalRotation() << std::endl;

  os << "    total displacement and rotation: " << r.displacement() << std::endl;
  os << r.rotation() << std::endl;

  for (const auto& idet : theDets) {
    const auto& comp = idet->components();

    for (unsigned int i = 0; i < comp.size(); ++i) {
      os << "     Det position, phi, r: " << comp[i]->globalPosition() << " , " << comp[i]->globalPosition().phi()
         << " , " << comp[i]->globalPosition().perp() << std::endl;
      os << "     local  position, phi, r: " << r.surface().toLocal(comp[i]->globalPosition()) << " , "
         << r.surface().toLocal(comp[i]->globalPosition()).phi() << " , "
         << r.surface().toLocal(comp[i]->globalPosition()).perp() << std::endl;
    }
  }

  return os;
}