-
Notifications
You must be signed in to change notification settings - Fork 499
Expand file tree
/
Copy pathAlignableSensorHMPID.cxx
More file actions
73 lines (64 loc) · 2.07 KB
/
AlignableSensorHMPID.cxx
File metadata and controls
73 lines (64 loc) · 2.07 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
// 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 AlignableSensorHMPID.h
/// @author ruben.shahoyan@cern.ch, michael.lettrich@cern.ch
/// @since 2021-02-01
/// @brief HMPID sensor (chamber)
#include "Align/AlignableSensorHMPID.h"
#include "Align/utils.h"
#include "Framework/Logger.h"
//#include "AliTrackPointArray.h"
//#include "AliESDtrack.h"
#include "Align/AlignmentPoint.h"
#include "Align/AlignableDetector.h"
ClassImp(o2::align::AlignableSensorHMPID);
using namespace o2::align::utils;
using namespace TMath;
namespace o2
{
namespace align
{
//_________________________________________________________
AlignableSensorHMPID::AlignableSensorHMPID(const char* name, int vid, int iid, int isec)
: AlignableSensor(name, vid, iid)
{
// def c-tor
}
//_________________________________________________________
AlignableSensorHMPID::~AlignableSensorHMPID()
{
// d-tor
}
//____________________________________________
void AlignableSensorHMPID::prepareMatrixT2L()
{
// creat T2L matrix
double loc[3] = {0, 0, 0}, glo[3];
getMatrixL2GIdeal().LocalToMaster(loc, glo);
double alp = ATan2(glo[1], glo[0]);
double x = Sqrt(glo[0] * glo[0] + glo[1] * glo[1]);
TGeoHMatrix t2l;
t2l.SetDx(x);
t2l.RotateZ(alp * RadToDeg());
const TGeoHMatrix& l2gi = getMatrixL2GIdeal().Inverse();
t2l.MultiplyLeft(&l2gi);
/*
const TGeoHMatrix* t2l = AliGeomManager::GetTracking2LocalMatrix(getVolID());
if (!t2l) {
Print("long");
AliFatalF("Failed to find T2L matrix for VID:%d %s",getVolID(),getSymName());
}
*/
setMatrixT2L(t2l);
//
}
} // namespace align
} // namespace o2