GRSISort "v4.0.0.5"
An extension of the ROOT analysis Framework
Loading...
Searching...
No Matches
TTip.cxx
Go to the documentation of this file.
1#include "TTip.h"
2
3#include <iostream>
4
5#include "TRandom.h"
6#include "TMath.h"
7#include "TClass.h"
8
9// Coordinates for the TIP detectors (in mm)
10// Assumes the target position is at (0,0,0)
11std::array<TVector3, 128> TTip::fPositionVectors = {
12 TVector3(12.8039, 0.0000, 91.1047),
13 TVector3(0.0000, 12.8039, 91.1047),
14 TVector3(-12.8039, 0.0000, 91.1047),
15 TVector3(0.0000, -12.8039, 91.1047),
16 TVector3(17.2004, 0.0000, 54.5526),
17 TVector3(8.60019, 14.896, 54.5526),
18 TVector3(-8.60019, 14.896, 54.5526),
19 TVector3(-17.2004, 0.0000, 54.5526),
20 TVector3(-8.60019, -14.896, 54.5526),
21 TVector3(8.60019, -14.896, 54.5526),
22 TVector3(31.1534, 0.0000, 47.972),
23 TVector3(26.9796, 15.5767, 47.972),
24 TVector3(15.5767, 26.9796, 47.972),
25 TVector3(0.0000, 31.1534, 47.972),
26 TVector3(-15.5767, 26.9796, 47.972),
27 TVector3(-26.9796, 15.5767, 47.972),
28 TVector3(-31.1534, 0.0000, 47.972),
29 TVector3(-26.9796, -15.5767, 47.972),
30 TVector3(-15.5767, -26.9796, 47.972),
31 TVector3(0.0000, -31.1534, 47.972),
32 TVector3(15.5767, -26.9796, 47.972),
33 TVector3(26.9796, -15.5767, 47.972),
34 TVector3(42.8403, 0.0000, 37.9019),
35 TVector3(39.5792, 16.3943, 37.9019),
36 TVector3(30.2926, 30.2926, 37.9019),
37 TVector3(16.3943, 39.5792, 37.9019),
38 TVector3(0.0000, 42.8403, 37.9019),
39 TVector3(-16.3943, 39.5792, 37.9019),
40 TVector3(-30.2926, 30.2926, 37.9019),
41 TVector3(-39.5792, 16.3943, 37.9019),
42 TVector3(-42.8403, 0.0000, 37.9019),
43 TVector3(-39.5792, -16.3943, 37.9019),
44 TVector3(-30.2926, -30.2926, 37.9019),
45 TVector3(-16.3943, -39.5792, 37.9019),
46 TVector3(0.0000, -42.8403, 37.9019),
47 TVector3(16.3943, -39.5792, 37.9019),
48 TVector3(30.2926, -30.2926, 37.9019),
49 TVector3(39.5792, -16.3943, 37.9019),
50 TVector3(51.411, 0.0000, 25.0748),
51 TVector3(48.8948, 15.8869, 25.0748),
52 TVector3(41.5924, 30.2186, 25.0748),
53 TVector3(30.2186, 41.5924, 25.0748),
54 TVector3(15.8869, 48.8948, 25.0748),
55 TVector3(0.0000, 51.411, 25.0748),
56 TVector3(-15.8869, 48.8948, 25.0748),
57 TVector3(-30.2186, 41.5924, 25.0748),
58 TVector3(-41.5924, 30.2186, 25.0748),
59 TVector3(-48.8948, 15.8869, 25.0748),
60 TVector3(-51.411, 0.0000, 25.0748),
61 TVector3(-48.8948, -15.8869, 25.0748),
62 TVector3(-41.5924, -30.2186, 25.0748),
63 TVector3(-30.2186, -41.5924, 25.0748),
64 TVector3(-15.8869, -48.8948, 25.0748),
65 TVector3(0.0000, -51.411, 25.0748),
66 TVector3(15.8869, -48.8948, 25.0748),
67 TVector3(30.2186, -41.5924, 25.0748),
68 TVector3(41.5924, -30.2186, 25.0748),
69 TVector3(48.8948, -15.8869, 25.0748),
70 TVector3(56.2422, 0.0000, 10.4239),
71 TVector3(52.8504, 19.236, 10.4239),
72 TVector3(43.084, 36.1518, 10.4239),
73 TVector3(28.1211, 48.7072, 10.4239),
74 TVector3(9.76635, 55.3877, 10.4239),
75 TVector3(-9.76635, 55.3877, 10.4239),
76 TVector3(-28.1211, 48.7072, 10.4239),
77 TVector3(-43.084, 36.1518, 10.4239),
78 TVector3(-52.8504, 19.236, 10.4239),
79 TVector3(-56.2422, 0.0000, 10.4239),
80 TVector3(-52.8504, -19.236, 10.4239),
81 TVector3(-43.084, -36.1518, 10.4239),
82 TVector3(-28.1211, -48.7072, 10.4239),
83 TVector3(-9.76635, -55.3877, 10.4239),
84 TVector3(9.76635, -55.3877, 10.4239),
85 TVector3(28.1211, -48.7072, 10.4239),
86 TVector3(43.084, -36.1518, 10.4239),
87 TVector3(52.8504, -19.236, 10.4239),
88 TVector3(56.9823, 0.0000, -4.98531),
89 TVector3(53.5459, 19.4891, -4.98531),
90 TVector3(43.651, 36.6275, -4.98531),
91 TVector3(28.4912, 49.3482, -4.98531),
92 TVector3(9.89488, 56.1166, -4.98531),
93 TVector3(-9.89488, 56.1166, -4.98531),
94 TVector3(-28.4912, 49.3482, -4.98531),
95 TVector3(-43.651, 36.6275, -4.98531),
96 TVector3(-53.5459, 19.4891, -4.98531),
97 TVector3(-56.9823, 0.0000, -4.98531),
98 TVector3(-53.5459, -19.4891, -4.98531),
99 TVector3(-43.651, -36.6275, -4.98531),
100 TVector3(-28.4912, -49.3482, -4.98531),
101 TVector3(-9.89488, -56.1166, -4.98531),
102 TVector3(9.89488, -56.1166, -4.98531),
103 TVector3(28.4912, -49.3482, -4.98531),
104 TVector3(43.651, -36.6275, -4.98531),
105 TVector3(53.5459, -19.4891, -4.98531),
106 TVector3(53.0722, 0.0000, -21.3349),
107 TVector3(47.8164, 23.0272, -21.3349),
108 TVector3(33.09, 41.4935, -21.3349),
109 TVector3(11.8097, 51.7416, -21.3349),
110 TVector3(-11.8097, 51.7416, -21.3349),
111 TVector3(-33.09, 41.4935, -21.3349),
112 TVector3(-47.8164, 23.0272, -21.3349),
113 TVector3(-53.0722, 0.0000, -21.3349),
114 TVector3(-47.8164, -23.0272, -21.3349),
115 TVector3(-33.09, -41.4935, -21.3349),
116 TVector3(-11.8097, -51.7416, -21.3349),
117 TVector3(11.8097, -51.7416, -21.3349),
118 TVector3(33.09, -41.4935, -21.3349),
119 TVector3(47.8164, -23.0272, -21.3349),
120 TVector3(43.6891, 0.0000, -36.9202),
121 TVector3(37.8359, 21.8446, -36.9202),
122 TVector3(21.8446, 37.8359, -36.9202),
123 TVector3(0.0000, 43.6891, -36.9202),
124 TVector3(-21.8446, 37.8359, -36.9202),
125 TVector3(-37.8359, 21.8446, -36.9202),
126 TVector3(-43.6891, 0.0000, -36.9202),
127 TVector3(-37.8359, -21.8446, -36.9202),
128 TVector3(-21.8446, -37.8359, -36.9202),
129 TVector3(0.0000, -43.6891, -36.9202),
130 TVector3(21.8446, -37.8359, -36.9202),
131 TVector3(37.8359, -21.8446, -36.9202),
132 TVector3(29.8869, 0.0000, -48.771),
133 TVector3(21.1332, 21.1332, -48.771),
134 TVector3(0.0000, 29.8869, -48.771),
135 TVector3(-21.1332, 21.1332, -48.771),
136 TVector3(-29.8869, 0.0000, -48.771),
137 TVector3(-21.1332, -21.1332, -48.771),
138 TVector3(0.0000, -29.8869, -48.771),
139 TVector3(21.1332, -21.1332, -48.771)};
140
141TTip::TTip(const TTip& rhs) : TDetector(rhs)
142{
143 rhs.Copy(*this);
144}
145
146TTip& TTip::operator=(const TTip& rhs)
147{
148 rhs.Copy(*this);
149 return *this;
150}
151
152void TTip::AddFragment(const std::shared_ptr<const TFragment>& frag, TChannel* chan)
153{
154 if(frag == nullptr || chan == nullptr) {
155 return;
156 }
157
158 auto* hit = new TTipHit(*frag);
159 hit->SetUpNumbering(chan); // Think about moving this to ctor
160 AddHit(hit);
161}
162
163void TTip::Print(Option_t*) const
164{
165 /// Prints out TTip members, currently only prints the multiplicity.
166 Print(std::cout);
167}
168
169void TTip::Print(std::ostream& out) const
170{
171 std::ostringstream str;
172 str << GetMultiplicity() << " hits" << std::endl;
173 out << str.str();
174}
175
176TVector3 TTip::GetPosition(const TTipHit& hit)
177{
178 return TTip::GetPosition(hit.GetTipChannel() - 1);
179}
180
181TVector3 TTip::GetPosition(int DetNbr)
182{
183 if(DetNbr > 127) {
184 return {0, 0, 1}; //handle invalid detector index as in TGriffin
185 }
186 return fPositionVectors[DetNbr];
187}
TDetector::Copy
void Copy(TObject &) const override
!
Definition TDetector.cxx:24
TDetector::GetMultiplicity
virtual Short_t GetMultiplicity() const
Definition TDetector.h:73
TDetector::AddHit
virtual void AddHit(TDetectorHit *hit)
Definition TDetector.h:63
TTipHit::GetTipChannel
Int_t GetTipChannel() const
!
Definition TTipHit.h:58
TTip
Definition TTip.h:22
TTip::operator=
TTip & operator=(const TTip &)
!
Definition TTip.cxx:146
TTip::fPositionVectors
static std::array< TVector3, 128 > fPositionVectors
!
Definition TTip.h:11
TTip::Print
void Print(Option_t *opt="") const override
Definition TTip.cxx:163
TTip::AddFragment
void AddFragment(const std::shared_ptr< const TFragment > &, TChannel *) override
!
Definition TTip.cxx:152
TTip::TTip
TTip()=default
TTip::GetPosition
static TVector3 GetPosition(int DetNbr)
!
Definition TTip.cxx:181