VXDGlobalParamInterface Class Reference

Very tentative interface for VXD. More...

#include <GlobalParam.h>

Inheritance diagram for VXDGlobalParamInterface:

Public Types
enum	E_VXDHierarchyType {   c_None = 0 ,   c_Flat = 1 ,   c_HalfShells = 2 ,   c_Full = 3 }
	Type of VXD hierarchy. More...

Public Member Functions
virtual void	setupAlignmentHierarchy (GlobalDerivativesHierarchy &hierarchy) override final
	Very tentative function: not yet used.
virtual void	readFromResult (std::vector< std::tuple< unsigned short, unsigned short, unsigned short, double > > &, GlobalParamVector &)
	Implement this to be called after Millepede finishes and all global parameters are filled in DB objects (which allow this) in GlobalParam.

Static Public Attributes
static E_VXDHierarchyType	s_hierarchyType
	What type of hierarchy to use for VXD?
static bool	s_enablePXD = true
	Enable PXD in hierarchy?
static bool	s_enableSVD = true
	Enable SVD in hierarchy?

Detailed Description

Very tentative interface for VXD.

Definition at line 73 of file GlobalParam.h.

Member Enumeration Documentation

E_VXDHierarchyType

enum E_VXDHierarchyType

Type of VXD hierarchy.

Definition at line 76 of file GlobalParam.h.

                              {
        c_None = 0,
        c_Flat = 1,
        c_HalfShells = 2,
        c_Full = 3
      };

Member Function Documentation

readFromResult()

virtual void readFromResult ( std::vector< std::tuple< unsigned short, unsigned short, unsigned short, double > > & , GlobalParamVector &  )

inlinevirtualinherited

Implement this to be called after Millepede finishes and all global parameters are filled in DB objects (which allow this) in GlobalParam.

The function takes two parameters, both of which are unused thus no name is given First parameter: Vector with a tuple representing <unique id, element, param, value> per result Second parameter: The GlobalParamVector which you have to update, extend or just read jsu after it is filled automatically (for DB objects which support it) and before its DBObjects are stored back into DB.

Definition at line 60 of file GlobalParam.h.

                                                      {};

setupAlignmentHierarchy()

void setupAlignmentHierarchy ( GlobalDerivativesHierarchy & hierarchy )

finaloverridevirtual

Very tentative function: not yet used.

This simpler hierarchy is as follows: Belle 2 / | | \ Ying Yang Pat Mat ... other sub-detectors / / | / | | \ | \ \ ......... sensors ........

So the hierarchy is as follows: Belle 2 / | | \ Ying Yang Pat Mat ... other sub-detectors / | / | | \ | \ ...... ladders ...... / / | / | | \ | \ \ ......... sensors ........

Reimplemented from IGlobalParamInterface.

Definition at line 174 of file GlobalParam.cc.

    {
      try {
        auto& rigidBodyHierarchy = dynamic_cast<RigidBodyHierarchy&>(hierarchy);
        auto& geo = VXD::GeoCache::getInstance();
        // Set-up hierarchy
        DBObjPtr<VXDAlignment> vxdAlignments;
        if (s_hierarchyType == c_None) {
          // Nothing to be done
          return;
 
        } else if (s_hierarchyType == c_Flat) {
          for (auto& sensor : geo.getListOfSensors()) {
            //TODO: Don't we have better way to distinguish PXD and SVD?
            // For PXD
            if (sensor.getLayerNumber() <= 2 and !s_enablePXD)
              continue;
            // SVD layer > 2
            if (sensor.getLayerNumber() > 2 and !s_enableSVD)
              continue;
 
            rigidBodyHierarchy.insertTGeoTransform<VXDAlignment, VXDAlignment>(sensor, VxdID(0, 0, 0, 1),
                geo.getSensorInfo(sensor).getTransformation(true));
          }
        } else if (s_hierarchyType == c_HalfShells) {
 
          for (auto& halfShellPlacement : geo.getHalfShellPlacements()) {
            // For PXD
            if (halfShellPlacement.first.getLayerNumber() <= 2 and !s_enablePXD)
              continue;
            // for SVD - layer > 2
            if (halfShellPlacement.first.getLayerNumber() > 2 and !s_enableSVD)
              continue;
 
            TGeoHMatrix trafoHalfShell = halfShellPlacement.second;
            trafoHalfShell *= geo.getTGeoFromRigidBodyParams(
                                vxdAlignments->get(halfShellPlacement.first, VXDAlignment::dU),
                                vxdAlignments->get(halfShellPlacement.first, VXDAlignment::dV),
                                vxdAlignments->get(halfShellPlacement.first, VXDAlignment::dW),
                                vxdAlignments->get(halfShellPlacement.first, VXDAlignment::dAlpha),
                                vxdAlignments->get(halfShellPlacement.first, VXDAlignment::dBeta),
                                vxdAlignments->get(halfShellPlacement.first, VXDAlignment::dGamma)
                              );
            rigidBodyHierarchy.insertTGeoTransform<VXDAlignment, alignment::EmptyGlobalParamSet>(halfShellPlacement.first, 0, trafoHalfShell);
 
            for (auto& ladderPlacement : geo.getLadderPlacements(halfShellPlacement.first)) {
              // Updated trafo
              TGeoHMatrix trafoLadder = ladderPlacement.second;
              trafoLadder *= geo.getTGeoFromRigidBodyParams(
                               vxdAlignments->get(ladderPlacement.first, VXDAlignment::dU),
                               vxdAlignments->get(ladderPlacement.first, VXDAlignment::dV),
                               vxdAlignments->get(ladderPlacement.first, VXDAlignment::dW),
                               vxdAlignments->get(ladderPlacement.first, VXDAlignment::dAlpha),
                               vxdAlignments->get(ladderPlacement.first, VXDAlignment::dBeta),
                               vxdAlignments->get(ladderPlacement.first, VXDAlignment::dGamma)
                             );
              // Do not insert ladder, but insert sensor directly into half-shells below, with the additional transformation from ladder
              //rigidBodyHierarchy.insertTGeoTransform<VXDAlignment, VXDAlignment>(ladderPlacement.first, halfShellPlacement.first, trafoLadder);
 
              for (auto& sensorPlacement : geo.getSensorPlacements(ladderPlacement.first)) {
                // Updated trafo
                TGeoHMatrix trafoSensor = sensorPlacement.second;
                trafoSensor *= geo.getTGeoFromRigidBodyParams(
                                 vxdAlignments->get(sensorPlacement.first, VXDAlignment::dU),
                                 vxdAlignments->get(sensorPlacement.first, VXDAlignment::dV),
                                 vxdAlignments->get(sensorPlacement.first, VXDAlignment::dW),
                                 vxdAlignments->get(sensorPlacement.first, VXDAlignment::dAlpha),
                                 vxdAlignments->get(sensorPlacement.first, VXDAlignment::dBeta),
                                 vxdAlignments->get(sensorPlacement.first, VXDAlignment::dGamma)
                               );
                rigidBodyHierarchy.insertTGeoTransform<VXDAlignment, VXDAlignment>(sensorPlacement.first, halfShellPlacement.first,
                    trafoLadder * trafoSensor);
 
 
              }
            }
          }
        } else if (s_hierarchyType == c_Full) {
 
          for (auto& halfShellPlacement : geo.getHalfShellPlacements()) {
            // For PXD
            if (halfShellPlacement.first.getLayerNumber() <= 2 and !s_enablePXD)
              continue;
            // for SVD - layer > 2
            if (halfShellPlacement.first.getLayerNumber() > 2 and !s_enableSVD)
              continue;
 
            TGeoHMatrix trafoHalfShell = halfShellPlacement.second;
            trafoHalfShell *= geo.getTGeoFromRigidBodyParams(
                                vxdAlignments->get(halfShellPlacement.first, VXDAlignment::dU),
                                vxdAlignments->get(halfShellPlacement.first, VXDAlignment::dV),
                                vxdAlignments->get(halfShellPlacement.first, VXDAlignment::dW),
                                vxdAlignments->get(halfShellPlacement.first, VXDAlignment::dAlpha),
                                vxdAlignments->get(halfShellPlacement.first, VXDAlignment::dBeta),
                                vxdAlignments->get(halfShellPlacement.first, VXDAlignment::dGamma)
                              );
            rigidBodyHierarchy.insertTGeoTransform<VXDAlignment, alignment::EmptyGlobalParamSet>(halfShellPlacement.first, 0, trafoHalfShell);
 
            for (auto& ladderPlacement : geo.getLadderPlacements(halfShellPlacement.first)) {
              // Updated trafo
              TGeoHMatrix trafoLadder = ladderPlacement.second;
              trafoLadder *= geo.getTGeoFromRigidBodyParams(
                               vxdAlignments->get(ladderPlacement.first, VXDAlignment::dU),
                               vxdAlignments->get(ladderPlacement.first, VXDAlignment::dV),
                               vxdAlignments->get(ladderPlacement.first, VXDAlignment::dW),
                               vxdAlignments->get(ladderPlacement.first, VXDAlignment::dAlpha),
                               vxdAlignments->get(ladderPlacement.first, VXDAlignment::dBeta),
                               vxdAlignments->get(ladderPlacement.first, VXDAlignment::dGamma)
                             );
              rigidBodyHierarchy.insertTGeoTransform<VXDAlignment, VXDAlignment>(ladderPlacement.first, halfShellPlacement.first, trafoLadder);
 
              for (auto& sensorPlacement : geo.getSensorPlacements(ladderPlacement.first)) {
                // Updated trafo
                TGeoHMatrix trafoSensor = sensorPlacement.second;
                trafoSensor *= geo.getTGeoFromRigidBodyParams(
                                 vxdAlignments->get(sensorPlacement.first, VXDAlignment::dU),
                                 vxdAlignments->get(sensorPlacement.first, VXDAlignment::dV),
                                 vxdAlignments->get(sensorPlacement.first, VXDAlignment::dW),
                                 vxdAlignments->get(sensorPlacement.first, VXDAlignment::dAlpha),
                                 vxdAlignments->get(sensorPlacement.first, VXDAlignment::dBeta),
                                 vxdAlignments->get(sensorPlacement.first, VXDAlignment::dGamma)
                               );
                rigidBodyHierarchy.insertTGeoTransform<VXDAlignment, VXDAlignment>(sensorPlacement.first, ladderPlacement.first, trafoSensor);
 
 
              }
            }
          }
        }
      } catch (const std::bad_cast& e) {
        // Bad cast
      }
    }

Member Data Documentation

s_enablePXD

bool s_enablePXD = true

static

Enable PXD in hierarchy?

Definition at line 85 of file GlobalParam.h.

s_enableSVD

bool s_enableSVD = true

static

Enable SVD in hierarchy?

Definition at line 87 of file GlobalParam.h.

s_hierarchyType

VXDGlobalParamInterface::E_VXDHierarchyType s_hierarchyType

static

Initial value:

=
      VXDGlobalParamInterface::E_VXDHierarchyType::c_Full

What type of hierarchy to use for VXD?

Definition at line 83 of file GlobalParam.h.

---

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

• alignment/include/GlobalParam.h
• alignment/src/GlobalParam.cc