Public Member Functions
	MeasuredGeometryProducer (const edm::ESHandle< DDCompactView > &idealCV, const edm::ESHandle< RPAlignmentCorrectionsData > &alignments)

DDCompactView *&	produce ()

Private Member Functions
void	applyCorrection (const DDLogicalPart &parent, const DDLogicalPart &child, const RPAlignmentCorrectionData &correction, DDTranslation &translation, DDRotationMatrix &rotation, const bool useMeasuredParent=true)

void	positionDetectors (void)

void	positionEverythingButDetectors (void)

Static Private Member Functions
static void	applyCorrectionToTransform (const RPAlignmentCorrectionData &correction, TMatrixD &transform)

static DDExpandedView *	delExpandedView (DDExpandedView *expandedView)

static int	getDetectorId (const int rpId, const int detNo)

static void	getGlobalTransform (const DDLogicalPart &part, const DDCompactView &compactView, TMatrixD &transform)

static int	getRPIdFromNamespace (const DDLogicalPart &part)

static bool	isDetector (const DDLogicalPart &part)

static bool	isRPBox (const DDLogicalPart &part)

static DDExpandedView *	newExpandedView (const DDCompactView &compactView, const DDLogicalPart &part)

static void	translRotFromTransform (DDTranslation &translation, DDRotationMatrix &rotation, const TMatrixD &transform)

static void	translRotToTransform (const DDTranslation &translation, const DDRotationMatrix &rotation, TMatrixD &transform)

Private Attributes
const RPAlignmentCorrectionsData *const	alignments

const DDCompactView &	idealCV

DDCompactView *	measuredCV

DDLogicalPart	root

Static Private Attributes
static bool	evRotationStoreState

Detailed Description

Definition at line 167 of file TotemRPGeometryESModule.cc.

Constructor & Destructor Documentation

MeasuredGeometryProducer::MeasuredGeometryProducer	(	const edm::ESHandle< DDCompactView > &	idealCV,
		const edm::ESHandle< RPAlignmentCorrectionsData > &	alignments
	)

inline

Definition at line 388 of file TotemRPGeometryESModule.cc.

References pyrootRender::root.

                                                                            : idealCV(*idealCV), alignments(alignments.isValid() ? &(*alignments) : NULL) {
             root = this->idealCV.root();
         }

Member Function Documentation

void MeasuredGeometryProducer::applyCorrection	(	const DDLogicalPart &	parent,
		const DDLogicalPart &	child,
		const RPAlignmentCorrectionData &	correction,
		DDTranslation &	translation,
		DDRotationMatrix &	rotation,
		const bool	useMeasuredParent = `true`
	)

inlineprivate

Definition at line 309 of file TotemRPGeometryESModule.cc.

References funct::C, and F().

                                                                                                              {
             TMatrixD C(4,4);    // child relative transform
             TMatrixD iP(4,4);   // ideal parent global transform
             TMatrixD mP(4,4);   // measured parent global transform
             TMatrixD F(4,4);    // final child transform
 
             translRotToTransform(translation, rotation, C);
 
             if (useMeasuredParent) 
                 getGlobalTransform(parent, *measuredCV, mP);
             else 
                 getGlobalTransform(parent, idealCV, mP);
             getGlobalTransform(parent, idealCV, iP);
 
             // global final transform
             F = iP * C;
             applyCorrectionToTransform(correction, F);
             // relative final transform
             mP.Invert();
             F = mP * F;
 
             translRotFromTransform(translation, rotation, F);
         }

static void MeasuredGeometryProducer::applyCorrectionToTransform	(	const RPAlignmentCorrectionData &	correction,
		TMatrixD &	transform
	)

inlinestaticprivate

Definition at line 296 of file TotemRPGeometryESModule.cc.

References RPAlignmentCorrectionData::getRotationMatrix(), RPAlignmentCorrectionData::getTranslation(), and idealTransformation::rotation.

                                                                                                                  {
             DDTranslation translation;
             DDRotationMatrix rotation;
 
             translRotFromTransform(translation, rotation, transform);
 
             translation = correction.getTranslation() + translation;
             rotation    = correction.getRotationMatrix() * rotation;
 
             translRotToTransform(translation, rotation, transform);
         }

static DDExpandedView* MeasuredGeometryProducer::delExpandedView ( DDExpandedView * expandedView )

inlinestaticprivate

Definition at line 209 of file TotemRPGeometryESModule.cc.

References DDI::Singleton< I >::instance(), and NULL.

                                                                              {
             delete expandedView;
             DDRotation::StoreT::instance().setReadOnly(evRotationStoreState);
             return NULL;
         }

static int MeasuredGeometryProducer::getDetectorId	(	const int	rpId,
		const int	detNo
	)

inlinestaticprivate

Definition at line 286 of file TotemRPGeometryESModule.cc.

                                                                          {
             return rpId * 10 + detNo;
         }

static void MeasuredGeometryProducer::getGlobalTransform	(	const DDLogicalPart &	part,
		const DDCompactView &	compactView,
		TMatrixD &	transform
	)

inlinestaticprivate

Definition at line 261 of file TotemRPGeometryESModule.cc.

References DDExpandedView::rotation(), create_public_lumi_plots::transform, and DDExpandedView::translation().

                                                                                                                          {
             DDExpandedView *expandedView = newExpandedView(compactView, part);
             translRotToTransform(expandedView->translation(), expandedView->rotation(), transform);
             expandedView = delExpandedView(expandedView);
         }

static int MeasuredGeometryProducer::getRPIdFromNamespace ( const DDLogicalPart & part )

inlinestaticprivate

Definition at line 280 of file TotemRPGeometryESModule.cc.

References DDBase< N, C >::name(), and DDName::ns().

                                                                           {
             int nsLength = part.name().ns().length();
             return atoi(part.name().ns().substr(nsLength - 3, nsLength).c_str());
         }

static bool MeasuredGeometryProducer::isDetector ( const DDLogicalPart & part )

inlinestaticprivate

Definition at line 275 of file TotemRPGeometryESModule.cc.

References DDD_TOTEM_RP_DETECTOR_NAME, DDName::name(), and DDBase< N, C >::name().

                                                                  {
             return (! part.name().name().compare(DDD_TOTEM_RP_DETECTOR_NAME));
         }

static bool MeasuredGeometryProducer::isRPBox ( const DDLogicalPart & part )

inlinestaticprivate

Definition at line 270 of file TotemRPGeometryESModule.cc.

References DDD_TOTEM_RP_PRIMARY_VACUUM_NAME, DDName::name(), and DDBase< N, C >::name().

                                                               {
             return (! part.name().name().compare(DDD_TOTEM_RP_PRIMARY_VACUUM_NAME));
         }

static DDExpandedView* MeasuredGeometryProducer::newExpandedView	(	const DDCompactView &	compactView,
		const DDLogicalPart &	part
	)

inlinestaticprivate

Definition at line 186 of file TotemRPGeometryESModule.cc.

References DDI::Singleton< I >::instance(), DDExpandedView::logicalPart(), DDName::name(), mergeVDriftHistosByStation::name, DDBase< N, C >::name(), DDExpandedView::next(), and DDName::ns().

                                                                                                             {
             evRotationStoreState = DDRotation::StoreT::instance().readOnly();
             DDExpandedView *expandedView = new DDExpandedView(compactView);
             // traverse the tree until name and ns are mached
             const string &name = part.name().name();
             const string &ns = part.name().ns();
             bool noMatch = true;
 
             noMatch = false;
             noMatch |= expandedView->logicalPart().name().name().compare(name);
             noMatch |= expandedView->logicalPart().name().ns().compare(ns);
             while (noMatch) {
                 expandedView->next();
                 noMatch = false;
                 noMatch |= expandedView->logicalPart().name().name().compare(name);
                 noMatch |= expandedView->logicalPart().name().ns().compare(ns);
             }
             return expandedView;
         }

void MeasuredGeometryProducer::positionDetectors ( void )

inlineprivate

Definition at line 361 of file TotemRPGeometryESModule.cc.

References DDanonymousRot(), NULL, and idealTransformation::rotation.

                                      {
             DDCompactView::graph_type::const_iterator it = idealCV.graph().begin_iter();
             DDCompactView::graph_type::const_iterator itEnd = idealCV.graph().end_iter();
             for (; it != itEnd; ++it) {
                 if (isDetector(it->to())) {
                     const DDLogicalPart from    = it->from();
                     const DDLogicalPart to      = it->to();
                     const int           copyNo      = it->edge()->copyno_;
                     const DDDivision    &division   = it->edge()->division();
                     DDTranslation       translation(it->edge()->trans());
                     DDRotationMatrix    &rotationMatrix = *(new DDRotationMatrix(it->edge()->rot()));
 
                     const int rpId  = getRPIdFromNamespace(from);
                     const int detId = getDetectorId(rpId, copyNo);
                     if (alignments != NULL) {
                         const RPAlignmentCorrectionData correction = alignments->GetFullSensorCorrection(detId);
                         applyCorrection(from, to, correction, 
                                 translation, rotationMatrix);
                     }
 
                     const DDRotation rotation = DDanonymousRot(&rotationMatrix);
                     measuredCV->position(to, from, copyNo, translation, rotation, &division);
                 }
             }
         }

void MeasuredGeometryProducer::positionEverythingButDetectors ( void )

inlineprivate

Definition at line 334 of file TotemRPGeometryESModule.cc.

References DDanonymousRot(), NULL, and idealTransformation::rotation.

                                                   {
             DDCompactView::graph_type::const_iterator it = idealCV.graph().begin_iter();
             DDCompactView::graph_type::const_iterator itEnd = idealCV.graph().end_iter();
             for (; it != itEnd; ++it) {
                 if (!isDetector(it->to())) {
                     const DDLogicalPart from    = it->from();
                     const DDLogicalPart to      = it->to();
                     const int           copyNo      = it->edge()->copyno_;
                     const DDDivision    &division   = it->edge()->division();
                     DDTranslation       translation(it->edge()->trans());
                     DDRotationMatrix    &rotationMatrix = *(new DDRotationMatrix(it->edge()->rot()));
 
                     if (isRPBox(to)) {
                         const int rpId = getRPIdFromNamespace(to);
                         if (alignments != NULL) {
                             const RPAlignmentCorrectionData correction = alignments->GetRPCorrection(rpId);
                             applyCorrection(from, to, correction,
                                     translation, rotationMatrix, false);
                         }
                     }
 
                     const DDRotation rotation = DDanonymousRot(&rotationMatrix);
                     measuredCV->position(to, from, copyNo, translation, rotation, &division);
                 }
             }
         }

DDCompactView*& MeasuredGeometryProducer::produce ( )

inline

Definition at line 393 of file TotemRPGeometryESModule.cc.

References pyrootRender::root.

Referenced by TotemRPGeometryESModule::ApplyAlignments().

                                   {
             // create DDCompactView for measured geometry
             // notice that this class is not responsible for deleting this object
             measuredCV = new DDCompactView(root);
             // CMSSW/DetectorDescription graph interface sucks, so instead of doing a one bfs
             // we go over the tree twice (this is needed, as final detector postions are
             // dependent on new positions of RP units).
             positionEverythingButDetectors();
             positionDetectors();
             return measuredCV;
         }

static void MeasuredGeometryProducer::translRotFromTransform	(	DDTranslation &	translation,
		DDRotationMatrix &	rotation,
		const TMatrixD &	transform
	)

inlinestaticprivate

Definition at line 248 of file TotemRPGeometryESModule.cc.

References i, and makeHLTPrescaleTable::values.

                                            {
             // set rotation
             double values[9];
             for (int i = 0; i < 9; ++i) {
                 values[i] = transform[i / 3][i % 3];
             }
             rotation.SetComponents(values, values + 9);
             // set translation
             translation.SetXYZ(transform[0][3], transform[1][3], transform[2][3]);
         }

static void MeasuredGeometryProducer::translRotToTransform	(	const DDTranslation &	translation,
		const DDRotationMatrix &	rotation,
		TMatrixD &	transform
	)

inlinestaticprivate

Definition at line 232 of file TotemRPGeometryESModule.cc.

References i, and makeHLTPrescaleTable::values.

                                      {
             // set rotation
             double values[9];
             rotation.GetComponents(values);
             for (int i = 0; i < 9; ++i) {
                 transform[i / 3][i % 3] = values[i];
             }
             // set translation
             transform[0][3] = translation.X();
             transform[1][3] = translation.Y();
             transform[2][3] = translation.Z();
             transform[3][3] = 1.;
         }