Alignment/MuonAlignment/src/AlignableGEMSuperChamber.cc - master - githubmirror/cmssw

Alignment/MuonAlignment/src/AlignableGEMSuperChamber.cc

46 строк2 KB

hyunyong

add GEM to MuonAlignment

06 мар 2021, 01:53

06 мар 2021, 01:53375f214

/* AlignableGEMSuperChamber
 * \author Hyunyong Kim - TAMU
 */
#include "Alignment/MuonAlignment/interface/AlignableGEMSuperChamber.h"
#include "Alignment/MuonAlignment/interface/AlignableGEMChamber.h"

AlignableGEMSuperChamber::AlignableGEMSuperChamber(const GeomDet* geomDet) : AlignableDet(geomDet, false) {
  theStructureType = align::AlignableGEMSuperChamber;
  const std::vector<const GeomDet*>& geomDets = geomDet->components();
  for (std::vector<const GeomDet*>::const_iterator idet = geomDets.begin(); idet != geomDets.end(); ++idet) {
    addComponent(new AlignableGEMChamber(*idet));
  }
  this->theSurface = geomDet->surface();
}

void AlignableGEMSuperChamber::update(const GeomDet* geomDet) {
  AlignableDet::update(geomDet);
  theSurface = geomDet->surface();
}

std::ostream& operator<<(std::ostream& os, const AlignableGEMSuperChamber& r) {
  const auto& theDets = r.components();

  os << "    This GEMSuperChamber 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;
}