WriteESAlignments Class Reference

#include <WriteESAlignments.h>

Public Types
typedef std::vector< AlignTransform >	AliVec
typedef std::vector< double >	DVec
typedef AlignTransform::Rotation	Rot
typedef AlignTransform::Translation	Trl

Public Member Functions
void	writeAlignments (const edm::EventSetup &eventSetup, const DVec &alphaVec, const DVec &betaVec, const DVec &gammaVec, const DVec &xtranslVec, const DVec &ytranslVec, const DVec &ztranslVec)
	WriteESAlignments (edm::ConsumesCollector &&cc)

Static Public Attributes
static const unsigned int	k_nA = EcalPreshowerGeometry::numberOfAlignments()

Private Member Functions
void	convert (const edm::EventSetup &eS, const DVec &a, const DVec &b, const DVec &g, const DVec &x, const DVec &y, const DVec &z, AliVec &va)
void	write (const Alignments &ali)

Private Attributes
edm::ESGetToken< Alignments, ESAlignmentRcd >	alignmentToken_
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	geometryToken_

Detailed Description

Definition at line 14 of file WriteESAlignments.h.

Member Typedef Documentation

AliVec

typedef std::vector<AlignTransform> WriteESAlignments::AliVec

Definition at line 16 of file WriteESAlignments.h.

DVec

typedef std::vector<double> WriteESAlignments::DVec

Definition at line 21 of file WriteESAlignments.h.

Rot

typedef AlignTransform::Rotation WriteESAlignments::Rot

Definition at line 19 of file WriteESAlignments.h.

Trl

typedef AlignTransform::Translation WriteESAlignments::Trl

Definition at line 18 of file WriteESAlignments.h.

Constructor & Destructor Documentation

WriteESAlignments()

WEA::WriteESAlignments

(

edm::ConsumesCollector &&

cc

)

Definition at line 17 of file WriteESAlignments.cc.

    : geometryToken_{cc.esConsumes<CaloGeometry, CaloGeometryRecord>(edm::ESInputTag{})},
      alignmentToken_{cc.esConsumes<Alignments, ESAlignmentRcd>(edm::ESInputTag{})} {}

Member Function Documentation

convert()

void WEA::convert ( const edm::EventSetup &  eS, const DVec &  a, const DVec &  b, const DVec &  g, const DVec &  x, const DVec &  y, const DVec &  z, WEA::AliVec &  va  )

private

Definition at line 54 of file WriteESAlignments.cc.

References a, alignmentToken_, alpha, b, HLT_2022v11_cff::beta, EcalPreshowerGeometry::detIdFromLocalAlignmentIndex(), DetId::Ecal, EcalPreshower, g, CustomPhysics_cfi::gamma, relativeConstraints::geom, geometryToken_, edm::EventSetup::getData(), mps_fire::i, triggerObjects_cff::id, k_nA, idealTransformation::rotation, x, y, and z.

Referenced by writeAlignments().

                                 {
  const auto& pG = eS.getData(geometryToken_);

  const CaloSubdetectorGeometry* geom(pG.getSubdetectorGeometry(DetId::Ecal, EcalPreshower));

  const auto& pA = eS.getData(alignmentToken_);
  const AliVec& vaPrev(pA.m_align);

  va.reserve(k_nA);
  for (unsigned int i(0); i != k_nA; ++i) {
    // ordering of i is: left, right, left, right,...
    // starting at ES- rear, then ES- front,
    // then ES+ front, then ES+ rear

    const ESDetId id(EcalPreshowerGeometry::detIdFromLocalAlignmentIndex(i));

    double zPlanePrev(geom->getGeometry(id)->getPosition().z());
    const double zAlignPrev(vaPrev[i].translation().z());
    const Trl q_I(0, 0, zPlanePrev - zAlignPrev);
    const Trl& s_p(vaPrev[i].translation());
    const Trl t_n(x[i], y[i], z[i]);
    const Rot G_p(vaPrev[i].rotation());
    const double gamma(g[i]);
    const double alpha(a[i]);
    const double beta(b[i]);

    const Rot L_n(  // New rotation in local frame!
        Rot(Rot(Rot().rotateZ(-gamma)).rotateX(-alpha)).rotateY(-beta));

    const Rot InvL_n(L_n.inverse());

    const Rot G_n(InvL_n * G_p);

    const Trl s_n(t_n + s_p + q_I - InvL_n * q_I);

    edm::LogVerbatim("WriteESAlignments") << "For i = " << i << ", q_I=" << q_I;
    edm::LogVerbatim("WriteESAlignments") << "For i = " << i << ", s_p=" << s_p;
    edm::LogVerbatim("WriteESAlignments") << "For i = " << i << ", alpha = " << 1000. * alpha << " mr";
    edm::LogVerbatim("WriteESAlignments") << "For i = " << i << ", beta  = " << 1000. * beta << " mr";
    edm::LogVerbatim("WriteESAlignments") << "For i = " << i << ", gamma = " << 1000. * gamma << " mr";
    edm::LogVerbatim("WriteESAlignments")
        << " For i = " << i << ", L_n = " << L_n << "   Euler angles=" << InvL_n.eulerAngles() << "\n";
    edm::LogVerbatim("WriteESAlignments") << "For i = " << i << ", t_n=" << t_n;
    edm::LogVerbatim("WriteESAlignments")
        << "For i = " << i << ", G_p=" << G_p << "   Euler angles=" << G_p.eulerAngles() << "\n";
    edm::LogVerbatim("WriteESAlignments")
        << " For i = " << i << ", InvL_n = " << InvL_n << "   Euler angles=" << InvL_n.eulerAngles() << "\n";
    edm::LogVerbatim("WriteESAlignments")
        << " For i =" << i << ", G_n = " << G_n << "    Euler angles=" << G_n.eulerAngles() << "\n";
    edm::LogVerbatim("WriteESAlignments") << " For i =" << i << ", s_n = " << s_n;
    edm::LogVerbatim("WriteESAlignments") << "++++++++++++++++++++++++++\n\n";

    va.emplace_back(AlignTransform(s_n, G_n, id));
  }
}

write()

void WEA::write ( const Alignments &  ali )

private

Definition at line 43 of file WriteESAlignments.cc.

References cond::service::PoolDBOutputService::currentTime(), edm::Service< T >::isAvailable(), and cond::service::PoolDBOutputService::writeOneIOV().

Referenced by writeAlignments().

                                     {
  edm::LogVerbatim("WriteESAlignments") << "Uploading ES alignments to the database";

  edm::Service<cond::service::PoolDBOutputService> poolDbService;

  if (!poolDbService.isAvailable())
    throw cms::Exception("NotAvailable") << "PoolDBOutputService not available";

  poolDbService->writeOneIOV<Alignments>(ali, poolDbService->currentTime(), "ESAlignmentRcd");
}

writeAlignments()

void WEA::writeAlignments	(	const edm::EventSetup &	eventSetup,
		const DVec &	alphaVec,
		const DVec &	betaVec,
		const DVec &	gammaVec,
		const DVec &	xtranslVec,
		const DVec &	ytranslVec,
		const DVec &	ztranslVec
	)

Definition at line 21 of file WriteESAlignments.cc.

References cms::cuda::assert(), convert(), options_cfi::eventSetup, k_nA, and write().

                                                     {
  assert(alphaVec.size() == k_nA);
  assert(betaVec.size() == k_nA);
  assert(gammaVec.size() == k_nA);
  assert(xtranslVec.size() == k_nA);
  assert(ytranslVec.size() == k_nA);
  assert(ztranslVec.size() == k_nA);

  Alignments ali;
  AliVec& vali(ali.m_align);

  convert(eventSetup, alphaVec, betaVec, gammaVec, xtranslVec, ytranslVec, ztranslVec, vali);

  write(ali);
}

Member Data Documentation

alignmentToken_

edm::ESGetToken<Alignments, ESAlignmentRcd> WriteESAlignments::alignmentToken_

private

Definition at line 48 of file WriteESAlignments.h.

Referenced by convert().

geometryToken_

edm::ESGetToken<CaloGeometry, CaloGeometryRecord> WriteESAlignments::geometryToken_

private

Definition at line 47 of file WriteESAlignments.h.

Referenced by convert().

k_nA

const unsigned int WEA::k_nA = EcalPreshowerGeometry::numberOfAlignments()

static

Definition at line 23 of file WriteESAlignments.h.

Referenced by convert(), and writeAlignments().