Belle II Software development
Public Types | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes | List of all members
TransformDataGlobalAligned Class Reference

Transformation data (global, aligned): singleton version. More...

#include <TransformDataGlobalAligned.h>

Inheritance diagram for TransformDataGlobalAligned:
TransformData

Public Types

enum  Displacement {
  c_None ,
  c_Displacement ,
  c_Alignment
}
 Source of displacement (alignment) data. More...
 

Public Member Functions

const HepGeom::Transform3D * getSectionTransform (int section) const
 Get section transformation.
 
const HepGeom::Transform3D * getLayerTransform (int section, int layer) const
 Get layer transformation.
 
const HepGeom::Transform3D * getSectorTransform (int section, int layer, int sector) const
 Get sector transformation.
 
const HepGeom::Transform3D * getPlaneTransform (int section, int layer, int sector, int plane) const
 Get plane transformation.
 
const HepGeom::Transform3D * getPlaneDisplacement (int section, int layer, int sector, int plane) const
 Get additional displacement for plane internal volumes.
 
const HepGeom::Transform3D * getSegmentTransform (int section, int layer, int sector, int plane, int segment) const
 Get segment transformation.
 
const HepGeom::Transform3D * getStripLocalToGlobal (KLMDigit *hit) const
 Get strip local to global transformation by hit.
 
const HepGeom::Transform3D * getStripGlobalToLocal (KLMDigit *hit) const
 Get strip global to local transformation by hit.
 
const HepGeom::Transform3D * getStripGlobalToLocal (int section, int layer, int sector, int plane, int strip) const
 Get strip global to local transformation by hit.
 
const HepGeom::Transform3D * getStripTransform (int section, int layer, int sector, int plane, int strip) const
 Get strip transformation.
 
bool intersection (KLMDigit *hit1, KLMDigit *hit2, HepGeom::Point3D< double > *cross, double *d1, double *d2, double *sd, bool segments=true) const
 Check if strips intersect, and find intersection point if yes.
 
int getSectorByPosition (int section, const HepGeom::Point3D< double > &position) const
 Get sector by position.
 
int getStripsByIntersection (const HepGeom::Point3D< double > &intersection, int *strip1, int *strip2) const
 Find strips by intersection.
 

Static Public Member Functions

static const TransformDataGlobalAligned & Instance ()
 Instantiation.
 

Private Member Functions

 TransformDataGlobalAligned ()
 Constructor.
 
 ~TransformDataGlobalAligned ()
 Destructor.
 
void transformsToGlobal ()
 Make transformations global from local.
 

Private Attributes

const EKLMElementNumbers * m_ElementNumbers
 Element numbers.
 
const GeometryData * m_GeoDat
 Geometry data.
 
HepGeom::Transform3D * m_Section
 Section transformations.
 
HepGeom::Transform3D ** m_Layer
 Layer transformations.
 
HepGeom::Transform3D *** m_Sector
 Sector transformations.
 
HepGeom::Transform3D **** m_Plane
 Plane transformations.
 
HepGeom::Transform3D **** m_PlaneDisplacement
 Plane internal volumes displacements.
 
HepGeom::Transform3D ***** m_Segment
 Segment transformations.
 
HepGeom::Transform3D ***** m_Strip
 Strip transformations.
 
HepGeom::Transform3D ***** m_StripInverse
 Inverse strip transformations.
 

Detailed Description

Transformation data (global, aligned): singleton version.

Definition at line 25 of file TransformDataGlobalAligned.h.

Member Enumeration Documentation

◆ Displacement

enum Displacement
inherited

Source of displacement (alignment) data.

Enumerator
c_None 

Displacement is not used.

c_Displacement 

Use displacement data (for geometry).

c_Alignment 

Use alignment data (for everything else).

Definition at line 42 of file TransformData.h.

42 {
43 c_None,
44 c_Displacement,
45 c_Alignment,
46 };
Belle2::EKLM::TransformData::c_None
@ c_None
Displacement is not used.
Definition: TransformData.h:43
Belle2::EKLM::TransformData::c_Displacement
@ c_Displacement
Use displacement data (for geometry).
Definition: TransformData.h:44
Belle2::EKLM::TransformData::c_Alignment
@ c_Alignment
Use alignment data (for everything else).
Definition: TransformData.h:45

Constructor & Destructor Documentation

◆ TransformDataGlobalAligned()

TransformDataGlobalAligned ( )
private

Constructor.

Definition at line 21 of file TransformDataGlobalAligned.cc.

21 :
22 EKLM::TransformData(true, EKLM::TransformData::c_Alignment)
23{
24}
Belle2::EKLM::TransformData
Transformation data.
Definition: TransformData.h:35

◆ ~TransformDataGlobalAligned()

~TransformDataGlobalAligned ( )
private

Destructor.

Definition at line 26 of file TransformDataGlobalAligned.cc.

27{
28}

Member Function Documentation

◆ getLayerTransform()

const HepGeom::Transform3D * getLayerTransform ( int  section,
int  layer 
) const
inherited

Get layer transformation.

Parameters
[in]sectionSection number.
[in]layerLayer number.

Definition at line 290 of file TransformData.cc.

291{
292 return &m_Layer[section - 1][layer - 1];
293}
Belle2::EKLM::TransformData::m_Layer
HepGeom::Transform3D ** m_Layer
Layer transformations.
Definition: TransformData.h:220

◆ getPlaneDisplacement()

const HepGeom::Transform3D * getPlaneDisplacement ( int  section,
int  layer,
int  sector,
int  plane 
) const
inherited

Get additional displacement for plane internal volumes.

Parameters
[in]sectionSection number.
[in]layerLayer number.
[in]sectorSector number.
[in]planePlane number.

Definition at line 307 of file TransformData.cc.

309{
310 return &m_PlaneDisplacement[section - 1][layer - 1][sector - 1][plane - 1];
311}
Belle2::EKLM::TransformData::m_PlaneDisplacement
HepGeom::Transform3D **** m_PlaneDisplacement
Plane internal volumes displacements.
Definition: TransformData.h:229

◆ getPlaneTransform()

const HepGeom::Transform3D * getPlaneTransform ( int  section,
int  layer,
int  sector,
int  plane 
) const
inherited

Get plane transformation.

Parameters
[in]sectionSection number.
[in]layerLayer number.
[in]sectorSector number.
[in]planePlane number.

Definition at line 301 of file TransformData.cc.

303{
304 return &m_Plane[section - 1][layer - 1][sector - 1][plane - 1];
305}
Belle2::EKLM::TransformData::m_Plane
HepGeom::Transform3D **** m_Plane
Plane transformations.
Definition: TransformData.h:226

◆ getSectionTransform()

const HepGeom::Transform3D * getSectionTransform ( int  section) const
inherited

Get section transformation.

Parameters
[in]sectionSection number.

Definition at line 284 of file TransformData.cc.

285{
286 return &m_Section[section - 1];
287}
Belle2::EKLM::TransformData::m_Section
HepGeom::Transform3D * m_Section
Section transformations.
Definition: TransformData.h:217

◆ getSectorByPosition()

int getSectorByPosition ( int  section,
const HepGeom::Point3D< double > &  position 
) const
inherited

Get sector by position.

Parameters
[in]sectionSection number.
[in]positionPosition.

Definition at line 416 of file TransformData.cc.

418{
419 int sector;
420 double x;
421 if (section == 1) {
422 x = position.x();
423 } else {
424 x = -position.x();
425 }
426 if (position.y() > 0) {
427 if (x > 0)
428 sector = 1;
429 else
430 sector = 2;
431 } else {
432 if (x > 0)
433 sector = 4;
434 else
435 sector = 3;
436 }
437 return sector;
438}

◆ getSectorTransform()

const HepGeom::Transform3D * getSectorTransform ( int  section,
int  layer,
int  sector 
) const
inherited

Get sector transformation.

Parameters
[in]sectionSection number.
[in]layerLayer number.
[in]sectorSector number.

Definition at line 295 of file TransformData.cc.

297{
298 return &m_Sector[section - 1][layer - 1][sector - 1];
299}
Belle2::EKLM::TransformData::m_Sector
HepGeom::Transform3D *** m_Sector
Sector transformations.
Definition: TransformData.h:223

◆ getSegmentTransform()

const HepGeom::Transform3D * getSegmentTransform ( int  section,
int  layer,
int  sector,
int  plane,
int  segment 
) const
inherited

Get segment transformation.

Parameters
[in]sectionSection number.
[in]layerLayer number.
[in]sectorSector number.
[in]planePlane number.
[in]segmentSegment number.

Definition at line 313 of file TransformData.cc.

316{
317 return &m_Segment[section - 1][layer - 1][sector - 1][plane - 1][segment - 1];
318}
Belle2::EKLM::TransformData::m_Segment
HepGeom::Transform3D ***** m_Segment
Segment transformations.
Definition: TransformData.h:232

◆ getStripGlobalToLocal() [1/2]

const HepGeom::Transform3D * getStripGlobalToLocal ( int  section,
int  layer,
int  sector,
int  plane,
int  strip 
) const
inherited

Get strip global to local transformation by hit.

Parameters
[in]sectionSection number.
[in]layerLayer number.
[in]sectorSector number.
[in]planePlane number.
[in]stripStrip number.
Returns
Transformation.

Definition at line 341 of file TransformData.cc.

343{
344 return &(m_StripInverse[section - 1][layer - 1][sector - 1][plane - 1]
345 [strip - 1]);
346}
Belle2::EKLM::TransformData::m_StripInverse
HepGeom::Transform3D ***** m_StripInverse
Inverse strip transformations.
Definition: TransformData.h:238

◆ getStripGlobalToLocal() [2/2]

const HepGeom::Transform3D * getStripGlobalToLocal ( KLMDigit *  hit) const
inherited

Get strip global to local transformation by hit.

Parameters
[in]hitHit.
Returns
Transformation.

Definition at line 334 of file TransformData.cc.

335{
336 return &(m_StripInverse[hit->getSection() - 1][hit->getLayer() - 1]
337 [hit->getSector() - 1][hit->getPlane() - 1][hit->getStrip() - 1]);
338}

◆ getStripLocalToGlobal()

const HepGeom::Transform3D * getStripLocalToGlobal ( KLMDigit *  hit) const
inherited

Get strip local to global transformation by hit.

Parameters
[in]hitHit.
Returns
Transformation.

Definition at line 327 of file TransformData.cc.

328{
329 return &(m_Strip[hit->getSection() - 1][hit->getLayer() - 1]
330 [hit->getSector() - 1][hit->getPlane() - 1][hit->getStrip() - 1]);
331}
Belle2::EKLM::TransformData::m_Strip
HepGeom::Transform3D ***** m_Strip
Strip transformations.
Definition: TransformData.h:235

◆ getStripsByIntersection()

int getStripsByIntersection ( const HepGeom::Point3D< double > &  intersection,
int *  strip1,
int *  strip2 
) const
inherited

Find strips by intersection.

Parameters
[in]intersectionIntersection point.
[out]strip1Strip 1 global number.
[out]strip2Strip 2 global number.
Returns
0 on success, -1 on error.

Definition at line 440 of file TransformData.cc.

442{
443 /* cppcheck-suppress variableScope */
444 int section, layer, sector, plane, segment, strip, stripSegment, stripGlobal;
445 /* cppcheck-suppress variableScope */
446 int nLayers, nPlanes, nSegments, nStripsSegment, minDistanceSegment;
447 double solenoidCenter, firstLayerCenter, layerShift;
448 /* cppcheck-suppress variableScope */
449 double x, y, z, l, minY, maxY;
450 double minDistance = 0, minDistanceNew, stripWidth;
451 HepGeom::Point3D<double> intersectionClhep, intersectionLocal;
452 intersectionClhep = intersection * CLHEP::cm / Unit::cm;
453 solenoidCenter = m_GeoDat->getSolenoidZ() / CLHEP::cm * Unit::cm;
454 if (intersection.z() < solenoidCenter)
455 section = 1;
456 else
457 section = 2;
458 firstLayerCenter =
459 (m_GeoDat->getSectionPosition()->getZ()
460 - 0.5 * m_GeoDat->getSectionPosition()->getLength()
461 + m_GeoDat->getLayerShiftZ()
462 - 0.5 * m_GeoDat->getLayerPosition()->getLength()) /
463 CLHEP::cm * Unit::cm;
464 layerShift = m_GeoDat->getLayerShiftZ() / CLHEP::cm * Unit::cm;
465 z = fabs(intersection.z() - solenoidCenter);
466 layer = round((z - firstLayerCenter) / layerShift) + 1;
467 if (layer <= 0)
468 layer = 1;
469 nLayers = m_GeoDat->getNDetectorLayers(section);
470 if (layer > nLayers)
471 layer = nLayers;
472 sector = getSectorByPosition(section, intersection);
473 nPlanes = m_GeoDat->getNPlanes();
474 nSegments = m_GeoDat->getNSegments();
475 nStripsSegment = m_ElementNumbers->getNStripsSegment();
476 stripWidth = m_GeoDat->getStripGeometry()->getWidth() / CLHEP::cm * Unit::cm;
477 minY = -stripWidth / 2;
478 maxY = (double(nStripsSegment) - 0.5) * stripWidth;
479 for (plane = 1; plane <= nPlanes; plane++) {
480 minDistanceSegment = 1;
481 for (segment = 1; segment <= nSegments; segment++) {
482 strip = (segment - 1) * nStripsSegment;
483 intersectionLocal = m_StripInverse[section - 1][layer - 1]
484 [sector - 1][plane - 1][strip] * intersectionClhep;
485 y = intersectionLocal.y() / CLHEP::cm * Unit::cm;
486 if (y < minY) {
487 minDistanceNew = minY - y;
488 } else if (y > maxY) {
489 minDistanceNew = y - maxY;
490 } else {
491 minDistance = 0;
492 minDistanceSegment = segment;
493 break;
494 }
495 if (segment == 1) {
496 minDistance = minDistanceNew;
497 } else if (minDistanceNew < minDistance) {
498 minDistance = minDistanceNew;
499 minDistanceSegment = segment;
500 }
501 }
502 /*
503 * The intersection is required to be strictly within a segment,
504 * this condition might be adjusted later.
505 */
506 if (minDistance > 0)
507 return -1;
508 strip = (minDistanceSegment - 1) * nStripsSegment;
509 intersectionLocal = m_StripInverse[section - 1][layer - 1]
510 [sector - 1][plane - 1][strip] * intersectionClhep;
511 y = intersectionLocal.y() / CLHEP::cm * Unit::cm;
512 stripSegment = ceil((y - 0.5 * stripWidth) / stripWidth);
513 if (stripSegment <= 0)
514 stripSegment = 1;
515 else if (stripSegment > nStripsSegment)
516 stripSegment = nStripsSegment;
517 strip = stripSegment + (minDistanceSegment - 1) * nStripsSegment;
518 intersectionLocal = m_StripInverse[section - 1][layer - 1]
519 [sector - 1][plane - 1][strip - 1] * intersectionClhep;
520 x = intersectionLocal.x();
521 l = m_GeoDat->getStripLength(strip);
522 /*
523 * The intersection is required to be strictly within the strip length,
524 * this condition might be adjusted later.
525 */
526 if (fabs(x) > 0.5 * l)
527 return -1;
528 stripGlobal = m_ElementNumbers->stripNumber(
529 section, layer, sector, plane, strip);
530 if (plane == 1)
531 *strip1 = stripGlobal;
532 else
533 *strip2 = stripGlobal;
534 }
535 return 0;
536}
Belle2::EKLMElementNumbers::stripNumber
int stripNumber(int section, int layer, int sector, int plane, int strip) const
Get strip number.
Definition: EKLMElementNumbers.cc:210
Belle2::EKLMElementNumbers::getNStripsSegment
static constexpr int getNStripsSegment()
Get number of strips in a segment.
Definition: EKLMElementNumbers.h:391
Belle2::EKLMGeometry::ElementPosition::getZ
double getZ() const
Get Z coordinate.
Definition: EKLMGeometry.h:197
Belle2::EKLMGeometry::ElementPosition::getLength
double getLength() const
Get length.
Definition: EKLMGeometry.h:146
Belle2::EKLMGeometry::StripGeometry::getWidth
double getWidth() const
Get width.
Definition: EKLMGeometry.h:1189
Belle2::EKLMGeometry::getLayerPosition
const ElementPosition * getLayerPosition() const
Get position data for layers.
Definition: EKLMGeometry.h:1799
Belle2::EKLMGeometry::getNPlanes
int getNPlanes() const
Get number of planes.
Definition: EKLMGeometry.h:1715
Belle2::EKLMGeometry::getLayerShiftZ
double getLayerShiftZ() const
Get Z distance between two layers.
Definition: EKLMGeometry.h:1807
Belle2::EKLMGeometry::getNDetectorLayers
int getNDetectorLayers(int section) const
Get number of detector layers.
Definition: EKLMGeometry.cc:293
Belle2::EKLMGeometry::getNSegments
int getNSegments() const
Get number of segments.
Definition: EKLMGeometry.h:1723
Belle2::EKLMGeometry::getSolenoidZ
double getSolenoidZ() const
Get solenoid center Z coordinate.
Definition: EKLMGeometry.h:1775
Belle2::EKLMGeometry::getStripGeometry
const StripGeometry * getStripGeometry() const
Get strip geometry data.
Definition: EKLMGeometry.h:1871
Belle2::EKLMGeometry::getSectionPosition
const ElementPosition * getSectionPosition() const
Get position data for sections.
Definition: EKLMGeometry.h:1791
Belle2::EKLM::GeometryData::getStripLength
double getStripLength(int strip) const
Get strip length.
Definition: GeometryData.h:71
Belle2::EKLM::TransformData::getSectorByPosition
int getSectorByPosition(int section, const HepGeom::Point3D< double > &position) const
Get sector by position.
Definition: TransformData.cc:416
Belle2::EKLM::TransformData::m_ElementNumbers
const EKLMElementNumbers * m_ElementNumbers
Element numbers.
Definition: TransformData.h:211
Belle2::EKLM::TransformData::m_GeoDat
const GeometryData * m_GeoDat
Geometry data.
Definition: TransformData.h:214
Belle2::EKLM::TransformData::intersection
bool intersection(KLMDigit *hit1, KLMDigit *hit2, HepGeom::Point3D< double > *cross, double *d1, double *d2, double *sd, bool segments=true) const
Check if strips intersect, and find intersection point if yes.
Definition: TransformData.cc:348
Belle2::Unit::cm
static const double cm
Standard units with the value = 1.
Definition: Unit.h:47

◆ getStripTransform()

const HepGeom::Transform3D * getStripTransform ( int  section,
int  layer,
int  sector,
int  plane,
int  strip 
) const
inherited

Get strip transformation.

Parameters
[in]sectionSection number.
[in]layerLayer number.
[in]sectorSector number.
[in]planePlane number.
[in]stripStrip number.

Definition at line 320 of file TransformData.cc.

322{
323 return &m_Strip[section - 1][layer - 1][sector - 1][plane - 1][strip - 1];
324}

◆ Instance()

const EKLM::TransformDataGlobalAligned & Instance ( )
static

Instantiation.

Definition at line 15 of file TransformDataGlobalAligned.cc.

16{
17 static EKLM::TransformDataGlobalAligned t;
18 return t;
19}
Belle2::EKLM::TransformDataGlobalAligned
Transformation data (global, aligned): singleton version.
Definition: TransformDataGlobalAligned.h:25

◆ intersection()

bool intersection ( KLMDigit *  hit1,
KLMDigit *  hit2,
HepGeom::Point3D< double > *  cross,
double *  d1,
double *  d2,
double *  sd,
bool  segments = true 
) const
inherited

Check if strips intersect, and find intersection point if yes.

Parameters
[in]hit1First hit.
[in]hit2Second hit.
[out]crossCrossing point (coordinate unit is cm).
[out]d1Distance from hit to SiPM of strip 1, cm.
[out]d2Distance from hit to SiPM of strip 2, cm.
[out]sdShortest distance between strips, cm. Or if second strip is closer to interaction point, then (- shortest distance).
[in]segmentsCheck if segments intersect (may need to turn this check off for debugging).
Returns
True if strips intersect.

Line parametric equations: (s1_1g) + (s1_2g - s1_1g) * t1 = a1 + v1 * t1, (s2_1g) + (s2_2g - s2_1g) * t2 = a2 + v2 * t2. Points of closest approach: t1 = ((v1,v2)*(d,v2) - v2^2*(d,v1)) / (v1^2*v2^2 - (v1,v2)^2) t2 = - ((v1,v2)*(d,v1) - v1^2*(d,v2)) / (v1^2*v2^2 - (v1,v2)^2) where d = a1 - a2.

Definition at line 348 of file TransformData.cc.

352{
353 /* Hits must be from the same sector, */
354 if (hit1->getSection() != hit2->getSection())
355 return false;
356 if (hit1->getLayer() != hit2->getLayer())
357 return false;
358 if (hit1->getSector() != hit2->getSector())
359 return false;
360 /* but different planes. */
361 if (hit1->getPlane() == hit2->getPlane())
362 return false;
363 /* Coordinates of strip 1 ends. */
364 double l1 = m_GeoDat->getStripLength(hit1->getStrip());
365 HepGeom::Point3D<double> s1_1(-0.5 * l1, 0.0, 0.0);
366 HepGeom::Point3D<double> s1_2(0.5 * l1, 0.0, 0.0);
367 const HepGeom::Transform3D* tr1 = getStripLocalToGlobal(hit1);
368 HepGeom::Point3D<double> s1_1g = (*tr1) * s1_1;
369 HepGeom::Point3D<double> s1_2g = (*tr1) * s1_2;
370 /* Coordinates of strip 2 ends. */
371 double l2 = m_GeoDat->getStripLength(hit2->getStrip());
372 HepGeom::Point3D<double> s2_1(-0.5 * l2, 0.0, 0.0);
373 HepGeom::Point3D<double> s2_2(0.5 * l2, 0.0, 0.0);
374 const HepGeom::Transform3D* tr2 = getStripLocalToGlobal(hit2);
375 HepGeom::Point3D<double> s2_1g = (*tr2) * s2_1;
376 HepGeom::Point3D<double> s2_2g = (*tr2) * s2_2;
386 HepGeom::Vector3D<double> v1 = s1_2g - s1_1g;
387 HepGeom::Vector3D<double> v2 = s2_2g - s2_1g;
388 HepGeom::Vector3D<double> d = s1_1g - s2_1g;
389 double v1sq = v1.mag2();
390 double v2sq = v2.mag2();
391 double v1dv2 = v1.dot(v2);
392 double ddv1 = d.dot(v1);
393 double ddv2 = d.dot(v2);
394 double den = v1sq * v2sq - v1dv2 * v1dv2;
395 double t1 = (v1dv2 * ddv2 - v2sq * ddv1) / den;
396 double t2 = (- v1dv2 * ddv1 + v1sq * ddv2) / den;
397 /* Segments do not intersect. */
398 if (segments) {
399 if (t1 < 0.0 || t1 > 1.0)
400 return false;
401 if (t2 < 0.0 || t2 > 1.0)
402 return false;
403 }
404 /* Segments intersect, set return values. */
405 HepGeom::Point3D<double> s1_cg = s1_1g + v1 * t1;
406 HepGeom::Point3D<double> s2_cg = s2_1g + v2 * t2;
407 *d1 = s1_2g.distance(s1_cg) / CLHEP::mm * Unit::mm;
408 *d2 = s2_2g.distance(s2_cg) / CLHEP::mm * Unit::mm;
409 *cross = 0.5 * (s1_cg + s2_cg) / CLHEP::mm * Unit::mm;
410 *sd = s1_cg.distance(s2_cg) / CLHEP::mm * Unit::mm;
411 if (s2_cg.mag2() < s1_cg.mag2())
412 *sd = - *sd;
413 return true;
414}
Belle2::EKLM::TransformData::getStripLocalToGlobal
const HepGeom::Transform3D * getStripLocalToGlobal(KLMDigit *hit) const
Get strip local to global transformation by hit.
Definition: TransformData.cc:327
Belle2::KLMDigit::getLayer
int getLayer() const
Get layer number.
Definition: KLMDigit.h:126
Belle2::KLMDigit::getSection
int getSection() const
Get section number.
Definition: KLMDigit.h:90
Belle2::KLMDigit::getPlane
int getPlane() const
Get plane number.
Definition: KLMDigit.h:144
Belle2::KLMDigit::getStrip
int getStrip() const
Get strip number.
Definition: KLMDigit.h:162
Belle2::KLMDigit::getSector
int getSector() const
Get sector number.
Definition: KLMDigit.h:108
Belle2::Unit::mm
static const double mm
[millimeters]
Definition: Unit.h:70
VXDTFFilterTest::v2
const std::vector< double > v2
MATLAB generated random vector.
Definition: decorrelationMatrix.cc:42
VXDTFFilterTest::v1
const std::vector< double > v1
MATLAB generated random vector.
Definition: decorrelationMatrix.cc:28

◆ transformsToGlobal()

void transformsToGlobal ( )
privateinherited

Make transformations global from local.

Definition at line 239 of file TransformData.cc.

240{
241 int iSection, iLayer, iSector, iPlane, iSegment, iStrip;
242 /* cppcheck-suppress variableScope */
243 int nSections, nLayers, nDetectorLayers, nSectors, nPlanes, nSegments, nStrips;
244 nSections = m_GeoDat->getNSections();
245 nLayers = m_GeoDat->getNLayers();
246 nSectors = m_GeoDat->getNSectors();
247 nPlanes = m_GeoDat->getNPlanes();
248 nSegments = m_GeoDat->getNSegments();
249 nStrips = m_GeoDat->getNStrips();
250 for (iSection = 0; iSection < nSections; iSection++) {
251 nDetectorLayers = m_GeoDat->getNDetectorLayers(iSection + 1);
252 for (iLayer = 0; iLayer < nLayers; iLayer++) {
253 m_Layer[iSection][iLayer] = m_Section[iSection] * m_Layer[iSection][iLayer];
254 for (iSector = 0; iSector < nSectors; iSector++) {
255 m_Sector[iSection][iLayer][iSector] =
256 m_Layer[iSection][iLayer] * m_Sector[iSection][iLayer][iSector];
257 if (iLayer >= nDetectorLayers)
258 continue;
259 for (iPlane = 0; iPlane < nPlanes; iPlane++) {
260 m_Plane[iSection][iLayer][iSector][iPlane] =
261 m_Sector[iSection][iLayer][iSector] *
262 m_Plane[iSection][iLayer][iSector][iPlane];
263 for (iSegment = 0; iSegment < nSegments; iSegment++) {
264 m_Segment[iSection][iLayer][iSector][iPlane][iSegment] =
265 m_Plane[iSection][iLayer][iSector][iPlane] *
266 m_PlaneDisplacement[iSection][iLayer][iSector][iPlane] *
267 m_Segment[iSection][iLayer][iSector][iPlane][iSegment];
268 }
269 for (iStrip = 0; iStrip < nStrips; iStrip++) {
270 m_Strip[iSection][iLayer][iSector][iPlane][iStrip] =
271 m_Plane[iSection][iLayer][iSector][iPlane] *
272 m_PlaneDisplacement[iSection][iLayer][iSector][iPlane] *
273 m_Strip[iSection][iLayer][iSector][iPlane][iStrip];
274 m_StripInverse[iSection][iLayer][iSector][iPlane][iStrip] =
275 m_Strip[iSection][iLayer][iSector][iPlane][iStrip].inverse();
276 }
277 }
278 }
279 }
280 }
281}
Belle2::EKLMGeometry::getNSections
int getNSections() const
Get number of sections.
Definition: EKLMGeometry.h:1685
Belle2::EKLMGeometry::getNLayers
int getNLayers() const
Get number of layers.
Definition: EKLMGeometry.h:1693
Belle2::EKLMGeometry::getNStrips
int getNStrips() const
Get number of strips.
Definition: EKLMGeometry.h:1739
Belle2::EKLMGeometry::getNSectors
int getNSectors() const
Get number of sectors.
Definition: EKLMGeometry.h:1707

Member Data Documentation

◆ m_ElementNumbers

const EKLMElementNumbers* m_ElementNumbers
privateinherited

Element numbers.

Definition at line 211 of file TransformData.h.

◆ m_GeoDat

const GeometryData* m_GeoDat
privateinherited

Geometry data.

Definition at line 214 of file TransformData.h.

◆ m_Layer

HepGeom::Transform3D** m_Layer
privateinherited

Layer transformations.

Definition at line 220 of file TransformData.h.

◆ m_Plane

HepGeom::Transform3D**** m_Plane
privateinherited

Plane transformations.

Definition at line 226 of file TransformData.h.

◆ m_PlaneDisplacement

HepGeom::Transform3D**** m_PlaneDisplacement
privateinherited

Plane internal volumes displacements.

Definition at line 229 of file TransformData.h.

◆ m_Section

HepGeom::Transform3D* m_Section
privateinherited

Section transformations.

Definition at line 217 of file TransformData.h.

◆ m_Sector

HepGeom::Transform3D*** m_Sector
privateinherited

Sector transformations.

Definition at line 223 of file TransformData.h.

◆ m_Segment

HepGeom::Transform3D***** m_Segment
privateinherited

Segment transformations.

Definition at line 232 of file TransformData.h.

◆ m_Strip

HepGeom::Transform3D***** m_Strip
privateinherited

Strip transformations.

Definition at line 235 of file TransformData.h.

◆ m_StripInverse

HepGeom::Transform3D***** m_StripInverse
privateinherited

Inverse strip transformations.

Definition at line 238 of file TransformData.h.

The documentation for this class was generated from the following files: