RPDisplacementGenerator Class Reference

This class introduces displacements of RP. It actually shifts and rotates PSimHit positions. It doesn't test whether the displaced hit is still on the detector's surface. This check takes place later in the process. It is done via edge effectivity. More...

#include <RPDisplacementGenerator.h>

Public Types
using	RotationMatrix = ROOT::Math::Rotation3D
using	Translation = ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D< double > >

Public Member Functions
PSimHit	displace (const PSimHit &)
	returns displaced PSimHit More...
	RPDisplacementGenerator (const edm::ParameterSet &, RPDetId, const edm::EventSetup &)

Static Public Member Functions
static uint32_t	rawToDecId (uint32_t raw)

Private Member Functions
Local3DPoint	displacePoint (const Local3DPoint &)
	displaces a point More...

Private Attributes
RPDetId	detId_
	ID of the detector. More...
bool	isOn_
	set to false to bypass displacements More...
RotationMatrix	rotation_
Translation	shift_
	displacement More...

Detailed Description

This class introduces displacements of RP. It actually shifts and rotates PSimHit positions. It doesn't test whether the displaced hit is still on the detector's surface. This check takes place later in the process. It is done via edge effectivity.

PSimHit points are given in the "local Det frame" (PSimHit.h)

Definition at line 26 of file RPDisplacementGenerator.h.

Member Typedef Documentation

RotationMatrix

using RPDisplacementGenerator::RotationMatrix = ROOT::Math::Rotation3D

Definition at line 28 of file RPDisplacementGenerator.h.

Translation

using RPDisplacementGenerator::Translation = ROOT::Math::DisplacementVector3D<ROOT::Math::Cartesian3D<double> >

Definition at line 29 of file RPDisplacementGenerator.h.

Constructor & Destructor Documentation

RPDisplacementGenerator()

RPDisplacementGenerator::RPDisplacementGenerator	(	const edm::ParameterSet &	ps,
		RPDetId	_detId,
		const edm::EventSetup &	iSetup
	)

Definition at line 19 of file RPDisplacementGenerator.cc.

References detId_, g, relativeConstraints::geom, edm::EventSetup::get(), edm::ParameterSet::getParameter(), CTPPSRPAlignmentCorrectionData::getRotationMatrix(), CTPPSRPAlignmentCorrectionData::getTranslation(), isOn_, dqm-mbProfile::log, LogDebug, rawToDecId(), rotation_, shift_, and edm::EventSetup::tryToGet().

    : detId_(_detId) {
  isOn_ = ps.getParameter<bool>("RPDisplacementOn");

  // read the alignment correction
  ESHandle<CTPPSRPAlignmentCorrectionsData> alignments;
  if (auto rec = iSetup.tryToGet<VeryForwardMisalignedGeometryRecord>()) {
    iSetup.get<VeryForwardMisalignedGeometryRecord>().get(alignments);
  }

  unsigned int decId = rawToDecId(detId_);

  math::XYZVectorD S_m;
  RotationMatrix R_m;

  if (alignments.isValid()) {
    const CTPPSRPAlignmentCorrectionData &ac = alignments->getFullSensorCorrection(decId);
    S_m = ac.getTranslation();
    R_m = ac.getRotationMatrix();
  } else
    isOn_ = false;

  // transform shift and rotation to the local coordinate frame
  ESHandle<CTPPSGeometry> geom;
  iSetup.get<VeryForwardRealGeometryRecord>().get(geom);
  const DetGeomDesc *g = geom->sensor(detId_);
  const RotationMatrix &R_l = g->rotation();
  rotation_ = R_l.Inverse() * R_m.Inverse() * R_l;
  shift_ = R_l.Inverse() * R_m.Inverse() * S_m;

  LogDebug("RPDisplacementGenerator").log([&](auto &log) {
    log << " det id = " << decId << ", isOn = " << isOn_ << "\n";
    if (isOn_) {
      log << " shift = " << shift_ << "\n";
      log << " rotation = " << rotation_ << "\n";
    }
  });
}

Member Function Documentation

displace()

PSimHit RPDisplacementGenerator::displace

(

const PSimHit &

input

)

returns displaced PSimHit

Definition at line 69 of file RPDisplacementGenerator.cc.

References displacePoint(), SiStripBadComponentsDQMServiceTemplate_cfg::ep, input, isOn_, and LogDebug.

                                                              {
  if (!isOn_)
    return input;

  const Local3DPoint &ep = input.entryPoint(), &xp = input.exitPoint();
  const Local3DPoint &dep = displacePoint(ep), &dxp = displacePoint(xp);

  LogDebug("RPDisplacementGenerator::displace\n") << " entry point: " << ep << " -> " << dep << "\n"
                                                  << " exit point : " << xp << " -> " << dxp << "\n";

  return PSimHit(dep,
                 dxp,
                 input.pabs(),
                 input.tof(),
                 input.energyLoss(),
                 input.particleType(),
                 input.detUnitId(),
                 input.trackId(),
                 input.thetaAtEntry(),
                 input.phiAtEntry(),
                 input.processType());
}

displacePoint()

Local3DPoint RPDisplacementGenerator::displacePoint ( const Local3DPoint &  p )

private

displaces a point

input is in mm, shifts are in mm too

Definition at line 60 of file RPDisplacementGenerator.cc.

References AlCaHLTBitMon_ParallelJobs::p, rotation_, shift_, and findQualityFiles::v.

Referenced by displace().

                                                                         {

  Translation v(p.x(), p.y(), p.z());
  v = rotation_ * v - shift_;

  return Local3DPoint(v.x(), v.y(), v.z());
}

rawToDecId()

uint32_t RPDisplacementGenerator::rawToDecId ( uint32_t  raw )

static

Definition at line 92 of file RPDisplacementGenerator.cc.

References CTPPSDetId::maskArm, TotemRPDetId::maskPlane, CTPPSDetId::maskRP, CTPPSDetId::maskStation, CTPPSDetId::startArmBit, TotemRPDetId::startPlaneBit, CTPPSDetId::startRPBit, and CTPPSDetId::startStationBit.

Referenced by DeadChannelsManager::isChannelDead(), and RPDisplacementGenerator().

                                                         {
  return ((raw >> CTPPSDetId::startArmBit) & CTPPSDetId::maskArm) * 1000 +
         ((raw >> CTPPSDetId::startStationBit) & CTPPSDetId::maskStation) * 100 +
         ((raw >> CTPPSDetId::startRPBit) & CTPPSDetId::maskRP) * 10 +
         ((raw >> TotemRPDetId::startPlaneBit) & TotemRPDetId::maskPlane);
}

Member Data Documentation

detId_

RPDetId RPDisplacementGenerator::detId_

private

ID of the detector.

Definition at line 40 of file RPDisplacementGenerator.h.

Referenced by RPDisplacementGenerator().

isOn_

bool RPDisplacementGenerator::isOn_

private

set to false to bypass displacements

Definition at line 47 of file RPDisplacementGenerator.h.

Referenced by displace(), and RPDisplacementGenerator().

rotation_

RotationMatrix RPDisplacementGenerator::rotation_

private

Definition at line 44 of file RPDisplacementGenerator.h.

Referenced by displacePoint(), and RPDisplacementGenerator().

shift_

Translation RPDisplacementGenerator::shift_

private

displacement

Definition at line 43 of file RPDisplacementGenerator.h.

Referenced by displacePoint(), and RPDisplacementGenerator().