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

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

Definition at line 28 of file RPDisplacementGenerator.h.

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

Definition at line 29 of file RPDisplacementGenerator.h.

Constructor & Destructor Documentation

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_, log, LogDebug, rawToDecId(), DetGeomDesc::rotation(), 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

PSimHit RPDisplacementGenerator::displace ( const PSimHit & input )

returns displaced PSimHit

Definition at line 69 of file RPDisplacementGenerator.cc.

References PSimHit::detUnitId(), displacePoint(), PSimHit::energyLoss(), PSimHit::entryPoint(), PSimHit::exitPoint(), input, isOn_, LogDebug, PSimHit::pabs(), PSimHit::particleType(), PSimHit::phiAtEntry(), PSimHit::processType(), PSimHit::thetaAtEntry(), PSimHit::tof(), and PSimHit::trackId().

                                                               {
   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());
 }

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 rotation_, shift_, findQualityFiles::v, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by displace().

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

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);
 }