|
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
|
/****************************************************************************
*
* This is a part of CTPPS offline software.
* Authors:
* Laurent Forthomme (laurent.forthomme@cern.ch)
* Nicola Minafra (nicola.minafra@cern.ch)
* Mateusz Szpyrka (mateusz.szpyrka@cern.ch)
*
****************************************************************************/
#include "DataFormats/CTPPSReco/interface/CTPPSTimingLocalTrack.h"
#include <cmath>
//----------------------------------------------------------------------------------------------------
//--- constructors
CTPPSTimingLocalTrack::CTPPSTimingLocalTrack() : num_hits_(0), num_planes_(0), valid_(true), t_(0.), t_sigma_(0.) {}
CTPPSTimingLocalTrack::CTPPSTimingLocalTrack(const math::XYZPoint& pos0,
const math::XYZPoint& pos0_sigma,
float t,
float t_sigma)
: pos0_(pos0), pos0_sigma_(pos0_sigma), num_hits_(0), num_planes_(0), valid_(false), t_(t), t_sigma_(t_sigma) {}
//--- interface member functions
bool CTPPSTimingLocalTrack::containsHit(const CTPPSTimingRecHit& recHit, float tolerance, CheckDimension check) const {
float xTolerance = pos0_sigma_.x() + (0.5 * recHit.xWidth()) + tolerance;
float yTolerance = pos0_sigma_.y() + (0.5 * recHit.yWidth()) + tolerance;
float zTolerance = pos0_sigma_.z() + (0.5 * recHit.zWidth()) + tolerance;
float xDiff = std::abs(pos0_.x() - recHit.x());
float yDiff = std::abs(pos0_.y() - recHit.y());
float zDiff = std::abs(pos0_.z() - recHit.z());
switch (check) {
case CheckDimension::x:
return xDiff < xTolerance;
case CheckDimension::y:
return yDiff < yTolerance;
case CheckDimension::all:
return xDiff < xTolerance && yDiff < yTolerance && zDiff < zTolerance;
}
return false;
}
//====================================================================================================
// Other methods implementation
//====================================================================================================
bool operator<(const CTPPSTimingLocalTrack& lhs, const CTPPSTimingLocalTrack& rhs) {
// start to sort by temporal coordinate
if (lhs.time() < rhs.time())
return true;
if (lhs.time() > rhs.time())
return false;
// then sort by x-position
if (lhs.x0() < rhs.x0())
return true;
if (lhs.x0() > rhs.x0())
return false;
// ...and y-position
if (lhs.y0() < rhs.y0())
return true;
if (lhs.y0() > rhs.y0())
return false;
// ...and z-position
return (lhs.z0() < rhs.z0());
}
|