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
typedef ROOT::Math::Rotation3D	RotationMatrix

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...

DDRotationMatrix	rotation_

DDTranslation	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 28 of file RPDisplacementGenerator.h.

Member Typedef Documentation

◆ RotationMatrix

typedef ROOT::Math::Rotation3D RPDisplacementGenerator::RotationMatrix

Definition at line 30 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.

     : 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 DDRotationMatrix &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";
     }
   });
 }

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

Member Function Documentation

◆ displace()

PSimHit RPDisplacementGenerator::displace ( const PSimHit & input )

returns displaced PSimHit

Definition at line 69 of file RPDisplacementGenerator.cc.

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

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

◆ 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.

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

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

Referenced by displace().

◆ rawToDecId()

uint32_t RPDisplacementGenerator::rawToDecId ( uint32_t raw )

static

Definition at line 92 of file RPDisplacementGenerator.cc.

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