-
Notifications
You must be signed in to change notification settings - Fork 496
Expand file tree
/
Copy pathAlignableSensor.h
More file actions
114 lines (105 loc) · 4.1 KB
/
AlignableSensor.h
File metadata and controls
114 lines (105 loc) · 4.1 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
//
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.
/// @file AlignableSensor.h
/// @author ruben.shahoyan@cern.ch, michael.lettrich@cern.ch
/// @since 2021-02-01
/// @brief End-chain alignment volume in detector branch, where the actual measurement is done.
#ifndef ALIGNABLESENSOR_H
#define ALIGNABLESENSOR_H
#include <TMath.h>
#include <TObjArray.h>
#include "ReconstructionDataFormats/BaseCluster.h"
#include "Align/AlignableVolume.h"
#include "Align/utils.h"
#include <vector>
namespace o2
{
namespace align
{
class AlignableDetector;
class AlignmentPoint;
class AlignableSensor : public AlignableVolume
{
public:
//
using ClusterF = o2::BaseCluster<float>;
AlignableSensor() = default;
AlignableSensor(const char* name, int vid, int iid, Controller* ctr);
~AlignableSensor() override = default;
//
void addChild(AlignableVolume*) override;
//
void setDetector(AlignableDetector* det) { mDet = det; }
AlignableDetector* getDetector() const { return mDet; }
//
int getSID() const { return mSID; }
void setSID(int s) { mSID = s; }
//
void incrementStat() { mNProcPoints++; }
//
// derivatives calculation
virtual void dPosTraDParCalib(const AlignmentPoint* pnt, double* deriv, int calibID, const AlignableVolume* parent = nullptr) const;
virtual void dPosTraDParGeom(const AlignmentPoint* pnt, double* deriv, const AlignableVolume* parent = nullptr) const;
//
virtual void dPosTraDParGeomLOC(const AlignmentPoint* pnt, double* deriv) const;
virtual void dPosTraDParGeomTRA(const AlignmentPoint* pnt, double* deriv) const;
virtual void dPosTraDParGeomLOC(const AlignmentPoint* pnt, double* deriv, const AlignableVolume* parent) const;
virtual void dPosTraDParGeomTRA(const AlignmentPoint* pnt, double* deriv, const AlignableVolume* parent) const;
//
void getModifiedMatrixT2LmodLOC(TGeoHMatrix& matMod, const double* delta) const;
void getModifiedMatrixT2LmodTRA(TGeoHMatrix& matMod, const double* delta) const;
//
virtual void applyAlignmentFromMPSol();
//
void setAddError(double y, double z)
{
mAddError[0] = y;
mAddError[1] = z;
}
const double* getAddError() const { return mAddError; }
//
bool isSensor() const override { return true; }
void Print(const Option_t* opt = "") const override;
//
virtual void updatePointByTrackInfo(AlignmentPoint* pnt, const trackParam_t* t) const;
void updateL2GRecoMatrices(const std::vector<o2::detectors::AlignParam>& algArr, const TGeoHMatrix* cumulDelta) override;
//
int finalizeStat() override;
//
virtual void prepareMatrixClAlg();
virtual void prepareMatrixClAlgReco();
const TGeoHMatrix& getMatrixClAlg() const { return mMatClAlg; }
const TGeoHMatrix& getMatrixClAlgReco() const { return mMatClAlgReco; }
void setMatrixClAlg(const TGeoHMatrix& m) { mMatClAlg = m; }
void setMatrixClAlgReco(const TGeoHMatrix& m) { mMatClAlgReco = m; }
//
protected:
//
bool IsSortable() const override { return true; }
int Compare(const TObject* a) const override;
//
// --------- dummies -----------
AlignableSensor(const AlignableSensor&);
AlignableSensor& operator=(const AlignableSensor&);
//
protected:
//
int mSID = -1; // sensor id in detector
double mAddError[2] = {0, 0}; // additional error increment for measurement
AlignableDetector* mDet = nullptr; // pointer on detector
TGeoHMatrix mMatClAlg; // reference cluster alignment matrix in tracking frame, i.e. the matrix we want to correct
TGeoHMatrix mMatClAlgReco; // reco-time cluster alignment matrix in tracking frame
//
ClassDefOverride(AlignableSensor, 1)
};
} // namespace align
} // namespace o2
#endif