HighLAND
ND280GeomId.cxx
1 //#include "TND280Log.hxx"
2 //#include "TGeometryId.hxx"
3 #include "ND280GeomId.hxx"
4 #include "ND280GeomIdDef.hxx"
5 #include <iostream>
6 
7 #define ND280Warn(v) do {std::cout << "WARNING:" << v << std::endl;} while(0)
8 #define ND280Error(v) do {std::cout << "ERROR:" << v << std::endl;} while(0)
9 
10 /*
11  This file is an exact copy of the file in oaEvent except:
12  - The two commented includes
13  - The definition of ND280Warm and ND280Error
14  - a simplification of the class TSettableGeometryId to avoid using TGeometryId
15 */
16 
17 /*
18 namespace {
19  // Provide access to the protected setter and getter methods.
20  class TSettableGeometryId: public ND::TGeometryId {
21  public:
22  TSettableGeometryId() {}
23  explicit TSettableGeometryId(const ND::TGeometryId& id)
24  : ND::TGeometryId(id) {}
25  void SetField(int v, int m, int b) {
26  ND::TGeometryId::SetFieldSafe(v,m,b);
27  }
28  int GetField(int m, int b) {
29  return ND::TGeometryId::GetFieldSafe(m,b);
30  }
31  };
32 }
33 */
34 
35 namespace {
36  // Provide access to the protected setter and getter methods.
37  class TSettableGeometryId {
38  ND::TGeometryId fGeometryId;
39  public:
40  TSettableGeometryId() : fGeometryId(0) {}
41  explicit TSettableGeometryId(ND::TGeometryId id) : fGeometryId(id) {}
42  void SetField(int val, int msb, int lsb) {
43  // Build a mask for this field
44  long mask = ((1<<(msb-lsb+1))-1) << lsb;
45  // Clear the bit field.
46  fGeometryId &= ~mask;
47  // Shift the value and set the field.
48  fGeometryId |= ((val << lsb)&mask);
49  }
50 
51  int GetField(int msb, int lsb) const {
52  int field = fGeometryId >> lsb;
53  int mask = (1<<(msb-lsb+1))-1;
54  return (field & mask);
55  }
56 
57  operator ND::TGeometryId() {
58  return fGeometryId;
59  }
60  };
61 }
62 
63 ND::TGeometryId ND::GeomId::Empty() {
64  return TGeometryId();
65 }
66 
67 int ND::GeomId::GetSubSystem(TGeometryId i) {
68  TSettableGeometryId id(i);
69  return id.GetField(ND::GeomId::Def::kDetectorIdMSB,
70  ND::GeomId::Def::kDetectorIdLSB);
71 }
72 
73 ////////////////////////////////////////////////////////////
74 // P0D Geometry Identifiers.
75 ////////////////////////////////////////////////////////////
76 
77 bool ND::GeomId::P0D::IsP0D(TGeometryId i) {
78  TSettableGeometryId id(i);
79  return (ND::GeomId::Def::kP0D
80  == id.GetField(ND::GeomId::Def::kDetectorIdMSB,
81  ND::GeomId::Def::kDetectorIdLSB));
82 }
83 
84 ND::TGeometryId ND::GeomId::P0D::Detector() {
85  TSettableGeometryId id;
86  id.SetField(ND::GeomId::Def::kP0D,
87  ND::GeomId::Def::kDetectorIdMSB,
88  ND::GeomId::Def::kDetectorIdLSB);
89  id.SetField(ND::GeomId::Def::P0D::kGlobal,
90  ND::GeomId::Def::P0D::kSeqIdMSB,
91  ND::GeomId::Def::P0D::kSeqIdLSB);
92  id.SetField(ND::GeomId::Def::P0D::Global::kP0D,
93  ND::GeomId::Def::P0D::Global::kSeqIdMSB,
94  ND::GeomId::Def::P0D::Global::kSeqIdLSB);
95  id.SetField(0,
96  ND::GeomId::Def::P0D::Global::kFieldMSB,
97  ND::GeomId::Def::P0D::Global::kFieldLSB);
98  return id;
99 }
100 
101 ND::TGeometryId ND::GeomId::P0D::SuperP0Dule(int sP0Dule) {
102  TSettableGeometryId id;
103  id.SetField(ND::GeomId::Def::kP0D,
104  ND::GeomId::Def::kDetectorIdMSB,
105  ND::GeomId::Def::kDetectorIdLSB);
106  id.SetField(ND::GeomId::Def::P0D::kGlobal,
107  ND::GeomId::Def::P0D::kSeqIdMSB,
108  ND::GeomId::Def::P0D::kSeqIdLSB);
109  id.SetField(ND::GeomId::Def::P0D::Global::kSP0Dule,
110  ND::GeomId::Def::P0D::Global::kSeqIdMSB,
111  ND::GeomId::Def::P0D::Global::kSeqIdLSB);
112  id.SetField(sP0Dule,
113  ND::GeomId::Def::P0D::Global::kFieldMSB,
114  ND::GeomId::Def::P0D::Global::kFieldLSB);
115  return id;
116 }
117 
118 int ND::GeomId::P0D::GetSuperP0Dule(TGeometryId i) {
119  TSettableGeometryId id(i);
120  if (ND::GeomId::Def::kP0D
121  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
122  ND::GeomId::Def::kDetectorIdLSB)) return -1;
123  if (ND::GeomId::Def::P0D::kGlobal
124  != id.GetField(ND::GeomId::Def::P0D::kSeqIdMSB,
125  ND::GeomId::Def::P0D::kSeqIdLSB)) return -1;
126  if (ND::GeomId::Def::P0D::Global::kSP0Dule
127  != id.GetField(ND::GeomId::Def::P0D::Global::kSeqIdMSB,
128  ND::GeomId::Def::P0D::Global::kSeqIdLSB)) return -1;
129  return id.GetField(ND::GeomId::Def::P0D::Global::kFieldMSB,
130  ND::GeomId::Def::P0D::Global::kFieldLSB);
131 }
132 
133 
134 ND::TGeometryId ND::GeomId::P0D::P0Dule(int p0dule) {
135  TSettableGeometryId id;
136  id.SetField(ND::GeomId::Def::kP0D,
137  ND::GeomId::Def::kDetectorIdMSB,
138  ND::GeomId::Def::kDetectorIdLSB);
139  id.SetField(ND::GeomId::Def::P0D::kGlobal,
140  ND::GeomId::Def::P0D::kSeqIdMSB,
141  ND::GeomId::Def::P0D::kSeqIdLSB);
142  id.SetField(ND::GeomId::Def::P0D::Global::kP0Dule,
143  ND::GeomId::Def::P0D::Global::kSeqIdMSB,
144  ND::GeomId::Def::P0D::Global::kSeqIdLSB);
145  id.SetField(p0dule,
146  ND::GeomId::Def::P0D::Global::kFieldMSB,
147  ND::GeomId::Def::P0D::Global::kFieldLSB);
148  return id;
149 }
150 
151 int ND::GeomId::P0D::GetP0Dule(TGeometryId i) {
152  TSettableGeometryId id(i);
153  if (ND::GeomId::Def::kP0D
154  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
155  ND::GeomId::Def::kDetectorIdLSB)) return -1;
156  if (ND::GeomId::Def::P0D::kGlobal
157  != id.GetField(ND::GeomId::Def::P0D::kSeqIdMSB,
158  ND::GeomId::Def::P0D::kSeqIdLSB)) return -1;
159  if (ND::GeomId::Def::P0D::Global::kP0Dule
160  != id.GetField(ND::GeomId::Def::P0D::Global::kSeqIdMSB,
161  ND::GeomId::Def::P0D::Global::kSeqIdLSB)) return -1;
162  return id.GetField(ND::GeomId::Def::P0D::Global::kFieldMSB,
163  ND::GeomId::Def::P0D::Global::kFieldLSB);
164 }
165 
166 ND::TGeometryId ND::GeomId::P0D::Target(int target) {
167  TSettableGeometryId id;
168  id.SetField(ND::GeomId::Def::kP0D,
169  ND::GeomId::Def::kDetectorIdMSB,
170  ND::GeomId::Def::kDetectorIdLSB);
171  id.SetField(ND::GeomId::Def::P0D::kGlobal,
172  ND::GeomId::Def::P0D::kSeqIdMSB,
173  ND::GeomId::Def::P0D::kSeqIdLSB);
174  id.SetField(ND::GeomId::Def::P0D::Global::kTarget,
175  ND::GeomId::Def::P0D::Global::kSeqIdMSB,
176  ND::GeomId::Def::P0D::Global::kSeqIdLSB);
177  id.SetField(target,
178  ND::GeomId::Def::P0D::Global::kFieldMSB,
179  ND::GeomId::Def::P0D::Global::kFieldLSB);
180  return id;
181 }
182 
183 int ND::GeomId::P0D::GetTarget(TGeometryId i) {
184  TSettableGeometryId id(i);
185  if (ND::GeomId::Def::kP0D
186  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
187  ND::GeomId::Def::kDetectorIdLSB)) return -1;
188  if (ND::GeomId::Def::P0D::kGlobal
189  != id.GetField(ND::GeomId::Def::P0D::kSeqIdMSB,
190  ND::GeomId::Def::P0D::kSeqIdLSB)) return -1;
191  if (ND::GeomId::Def::P0D::Global::kTarget
192  != id.GetField(ND::GeomId::Def::P0D::Global::kSeqIdMSB,
193  ND::GeomId::Def::P0D::Global::kSeqIdLSB)) return -1;
194  return id.GetField(ND::GeomId::Def::P0D::Global::kFieldMSB,
195  ND::GeomId::Def::P0D::Global::kFieldLSB);
196 }
197 
198 ND::TGeometryId ND::GeomId::P0D::ECalRadiator(int radiator) {
199  TSettableGeometryId id;
200  id.SetField(ND::GeomId::Def::kP0D,
201  ND::GeomId::Def::kDetectorIdMSB,
202  ND::GeomId::Def::kDetectorIdLSB);
203  id.SetField(ND::GeomId::Def::P0D::kGlobal,
204  ND::GeomId::Def::P0D::kSeqIdMSB,
205  ND::GeomId::Def::P0D::kSeqIdLSB);
206  id.SetField(ND::GeomId::Def::P0D::Global::kECalRadiator,
207  ND::GeomId::Def::P0D::Global::kSeqIdMSB,
208  ND::GeomId::Def::P0D::Global::kSeqIdLSB);
209  id.SetField(radiator,
210  ND::GeomId::Def::P0D::Global::kFieldMSB,
211  ND::GeomId::Def::P0D::Global::kFieldLSB);
212  return id;
213 }
214 
215 int ND::GeomId::P0D::GetECalRadiator(TGeometryId i) {
216  TSettableGeometryId id(i);
217  if (ND::GeomId::Def::kP0D
218  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
219  ND::GeomId::Def::kDetectorIdLSB)) return -1;
220  if (ND::GeomId::Def::P0D::kGlobal
221  != id.GetField(ND::GeomId::Def::P0D::kSeqIdMSB,
222  ND::GeomId::Def::P0D::kSeqIdLSB)) return -1;
223  if (ND::GeomId::Def::P0D::Global::kECalRadiator
224  != id.GetField(ND::GeomId::Def::P0D::Global::kSeqIdMSB,
225  ND::GeomId::Def::P0D::Global::kSeqIdLSB)) return -1;
226  return id.GetField(ND::GeomId::Def::P0D::Global::kFieldMSB,
227  ND::GeomId::Def::P0D::Global::kFieldLSB);
228 }
229 
230 ND::TGeometryId ND::GeomId::P0D::TargetRadiator(int targetRadiator) {
231  TSettableGeometryId id;
232  id.SetField(ND::GeomId::Def::kP0D,
233  ND::GeomId::Def::kDetectorIdMSB,
234  ND::GeomId::Def::kDetectorIdLSB);
235  id.SetField(ND::GeomId::Def::P0D::kGlobal,
236  ND::GeomId::Def::P0D::kSeqIdMSB,
237  ND::GeomId::Def::P0D::kSeqIdLSB);
238  id.SetField(ND::GeomId::Def::P0D::Global::kTargetRadiator,
239  ND::GeomId::Def::P0D::Global::kSeqIdMSB,
240  ND::GeomId::Def::P0D::Global::kSeqIdLSB);
241  id.SetField(targetRadiator,
242  ND::GeomId::Def::P0D::Global::kFieldMSB,
243  ND::GeomId::Def::P0D::Global::kFieldLSB);
244  return id;
245 }
246 
247 int ND::GeomId::P0D::GetTargetRadiator(TGeometryId i) {
248  TSettableGeometryId id(i);
249  if (ND::GeomId::Def::kP0D
250  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
251  ND::GeomId::Def::kDetectorIdLSB)) return -1;
252  if (ND::GeomId::Def::P0D::kGlobal
253  != id.GetField(ND::GeomId::Def::P0D::kSeqIdMSB,
254  ND::GeomId::Def::P0D::kSeqIdLSB)) return -1;
255  if (ND::GeomId::Def::P0D::Global::kTargetRadiator
256  != id.GetField(ND::GeomId::Def::P0D::Global::kSeqIdMSB,
257  ND::GeomId::Def::P0D::Global::kSeqIdLSB)) return -1;
258  return id.GetField(ND::GeomId::Def::P0D::Global::kFieldMSB,
259  ND::GeomId::Def::P0D::Global::kFieldLSB);
260 }
261 
262 ND::TGeometryId ND::GeomId::P0D::Bar(int P0Dule,
263  int layer,
264  int bar) {
265  TSettableGeometryId id;
266  id.SetField(ND::GeomId::Def::kP0D,
267  ND::GeomId::Def::kDetectorIdMSB,
268  ND::GeomId::Def::kDetectorIdLSB);
269  id.SetField(ND::GeomId::Def::P0D::kScintillator,
270  ND::GeomId::Def::P0D::kSeqIdMSB,
271  ND::GeomId::Def::P0D::kSeqIdLSB);
272  id.SetField(0,
273  ND::GeomId::Def::P0D::Bar::kSP0DuleMSB,
274  ND::GeomId::Def::P0D::Bar::kSP0DuleLSB);
275  id.SetField(P0Dule,
276  ND::GeomId::Def::P0D::Bar::kP0DuleMSB,
277  ND::GeomId::Def::P0D::Bar::kP0DuleLSB);
278  id.SetField(layer,
279  ND::GeomId::Def::P0D::Bar::kLayerMSB,
280  ND::GeomId::Def::P0D::Bar::kLayerLSB);
281  id.SetField(bar,
282  ND::GeomId::Def::P0D::Bar::kBarMSB,
283  ND::GeomId::Def::P0D::Bar::kBarLSB);
284  return id;
285 }
286 
287 int ND::GeomId::P0D::GetBarP0Dule(TGeometryId i) {
288  TSettableGeometryId id(i);
289  if (ND::GeomId::Def::kP0D
290  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
291  ND::GeomId::Def::kDetectorIdLSB)) return -1;
292  if (ND::GeomId::Def::P0D::kScintillator
293  != id.GetField(ND::GeomId::Def::P0D::kSeqIdMSB,
294  ND::GeomId::Def::P0D::kSeqIdLSB)) return -1;
295  return id.GetField(ND::GeomId::Def::P0D::Bar::kP0DuleMSB,
296  ND::GeomId::Def::P0D::Bar::kP0DuleLSB);
297 }
298 
299 int ND::GeomId::P0D::GetBarLayer(TGeometryId i) {
300  TSettableGeometryId id(i);
301  if (ND::GeomId::Def::kP0D
302  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
303  ND::GeomId::Def::kDetectorIdLSB)) return -1;
304  if (ND::GeomId::Def::P0D::kScintillator
305  != id.GetField(ND::GeomId::Def::P0D::kSeqIdMSB,
306  ND::GeomId::Def::P0D::kSeqIdLSB)) return -1;
307  return id.GetField(ND::GeomId::Def::P0D::Bar::kLayerMSB,
308  ND::GeomId::Def::P0D::Bar::kLayerLSB);
309 }
310 
311 int ND::GeomId::P0D::GetBarNumber(TGeometryId i) {
312  TSettableGeometryId id(i);
313  if (ND::GeomId::Def::kP0D
314  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
315  ND::GeomId::Def::kDetectorIdLSB)) return -1;
316  if (ND::GeomId::Def::P0D::kScintillator
317  != id.GetField(ND::GeomId::Def::P0D::kSeqIdMSB,
318  ND::GeomId::Def::P0D::kSeqIdLSB)) return -1;
319  return id.GetField(ND::GeomId::Def::P0D::Bar::kBarMSB,
320  ND::GeomId::Def::P0D::Bar::kBarLSB);
321 }
322 
323 ////////////////////////////////////////////////////////////
324 // TPC Geometry Identifiers
325 ////////////////////////////////////////////////////////////
326 
327 bool ND::GeomId::TPC::IsTPC(TGeometryId i) {
328  TSettableGeometryId id(i);
329  return (ND::GeomId::Def::kTPC
330  == id.GetField(ND::GeomId::Def::kDetectorIdMSB,
331  ND::GeomId::Def::kDetectorIdLSB));
332 }
333 
334 ND::TGeometryId ND::GeomId::TPC::Module(int tpc) {
335  TSettableGeometryId id;
336  if (tpc<0 || tpc>2) {
337  ND280Warn("TPC module out of range [0,2]: " << tpc);
338  return id;
339  }
340  id.SetField(ND::GeomId::Def::kTPC,
341  ND::GeomId::Def::kDetectorIdMSB,
342  ND::GeomId::Def::kDetectorIdLSB);
343  id.SetField(ND::GeomId::Def::TPC::kGlobal,
344  ND::GeomId::Def::TPC::kSeqIdMSB,
345  ND::GeomId::Def::TPC::kSeqIdLSB);
346  id.SetField(ND::GeomId::Def::TPC::Global::kTPC,
347  ND::GeomId::Def::TPC::Global::kSeqIdMSB,
348  ND::GeomId::Def::TPC::Global::kSeqIdLSB);
349  id.SetField(tpc,
350  ND::GeomId::Def::TPC::Global::kFieldMSB,
351  ND::GeomId::Def::TPC::Global::kFieldLSB);
352  return id;
353 }
354 
355 int ND::GeomId::TPC::GetModule(TGeometryId i) {
356  TSettableGeometryId id(i);
357  if (ND::GeomId::Def::kTPC
358  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
359  ND::GeomId::Def::kDetectorIdLSB)) return -1;
360  if (ND::GeomId::Def::TPC::kGlobal
361  != id.GetField(ND::GeomId::Def::TPC::kSeqIdMSB,
362  ND::GeomId::Def::TPC::kSeqIdLSB)) return -1;
363  if (ND::GeomId::Def::TPC::Global::kTPC
364  != id.GetField(ND::GeomId::Def::TPC::Global::kSeqIdMSB,
365  ND::GeomId::Def::TPC::Global::kSeqIdLSB)) return -1;
366  return id.GetField(ND::GeomId::Def::TPC::Global::kFieldMSB,
367  ND::GeomId::Def::TPC::Global::kFieldLSB);
368 }
369 
370 ND::TGeometryId ND::GeomId::TPC::TPC1() {return Module(0);}
371 ND::TGeometryId ND::GeomId::TPC::TPC2() {return Module(1);}
372 ND::TGeometryId ND::GeomId::TPC::TPC3() {return Module(2);}
373 
374 bool ND::GeomId::TPC::IsTPC1(TGeometryId i) {
375  return (ND::GeomId::TPC::GetModule(i) == 0);
376 }
377 
378 bool ND::GeomId::TPC::IsTPC2(TGeometryId i) {
379  return (ND::GeomId::TPC::GetModule(i) == 1);
380 }
381 
382 bool ND::GeomId::TPC::IsTPC3(TGeometryId i) {
383  return (ND::GeomId::TPC::GetModule(i) == 2);
384 }
385 
386 ND::TGeometryId ND::GeomId::TPC::MicroMega(int tpc, int half, int mm) {
387  TSettableGeometryId id;
388  if (tpc<0 || tpc>2) {
389  ND280Warn("TPC module out of range [0,2]: " << tpc);
390  return id;
391  }
392  if (half<0 || half>1) {
393  ND280Warn("TPC half out of range [0,1]: " << half);
394  return id;
395  }
396  if (mm<0 || mm>11) {
397  ND280Warn("TPC micromega out of range [0,11]: " << mm);
398  return id;
399  }
400  id.SetField(ND::GeomId::Def::kTPC,
401  ND::GeomId::Def::kDetectorIdMSB,
402  ND::GeomId::Def::kDetectorIdLSB);
403  id.SetField(ND::GeomId::Def::TPC::kPad,
404  ND::GeomId::Def::TPC::kSeqIdMSB,
405  ND::GeomId::Def::TPC::kSeqIdLSB);
406  id.SetField(tpc,
407  ND::GeomId::Def::TPC::Pad::kTPCMSB,
408  ND::GeomId::Def::TPC::Pad::kTPCLSB);
409  id.SetField(half,
410  ND::GeomId::Def::TPC::Pad::kHalfMSB,
411  ND::GeomId::Def::TPC::Pad::kHalfLSB);
412  id.SetField(mm,
413  ND::GeomId::Def::TPC::Pad::kMMegaMSB,
414  ND::GeomId::Def::TPC::Pad::kMMegaLSB);
415  id.SetField(0,
416  ND::GeomId::Def::TPC::Pad::kPadFlagMSB,
417  ND::GeomId::Def::TPC::Pad::kPadFlagLSB);
418  id.SetField(0,
419  ND::GeomId::Def::TPC::Pad::kPadMSB,
420  ND::GeomId::Def::TPC::Pad::kPadLSB);
421  return id;
422 }
423 
424 bool ND::GeomId::TPC::IsMicroMega(ND::TGeometryId i) {
425  TSettableGeometryId id(i);
426  int detector = id.GetField(ND::GeomId::Def::kDetectorIdMSB,
427  ND::GeomId::Def::kDetectorIdLSB);
428  if (detector != ND::GeomId::Def::kTPC) return false;
429  /// This may be a pad.
430  int seqId = id.GetField(ND::GeomId::Def::TPC::kSeqIdMSB,
431  ND::GeomId::Def::TPC::kSeqIdLSB);
432  if (seqId != ND::GeomId::Def::TPC::kPad) return false;
433  int padFlag = id.GetField(ND::GeomId::Def::TPC::Pad::kPadFlagMSB,
434  ND::GeomId::Def::TPC::Pad::kPadFlagLSB);
435  if (padFlag) return false;
436  return true;
437 }
438 
439 ND::TGeometryId ND::GeomId::TPC::Pad(int tpc, int half, int mm, int pad) {
440  TSettableGeometryId id;
441  if (tpc<0 || tpc>2) {
442  ND280Warn("TPC module out of range [0,2]: " << tpc);
443  return id;
444  }
445  if (half<0 || half>1) {
446  ND280Warn("TPC half out of range [0,1]: " << half);
447  return id;
448  }
449  if (mm<0 || mm>11) {
450  ND280Warn("TPC micromega out of range [0,11]: " << mm);
451  return id;
452  }
453  if (pad<0 || pad>1727) {
454  ND280Warn("TPC pad out of range [0,1727]: " << pad);
455  return id;
456  }
457  id.SetField(ND::GeomId::Def::kTPC,
458  ND::GeomId::Def::kDetectorIdMSB,
459  ND::GeomId::Def::kDetectorIdLSB);
460  id.SetField(ND::GeomId::Def::TPC::kPad,
461  ND::GeomId::Def::TPC::kSeqIdMSB,
462  ND::GeomId::Def::TPC::kSeqIdLSB);
463  id.SetField(tpc,
464  ND::GeomId::Def::TPC::Pad::kTPCMSB,
465  ND::GeomId::Def::TPC::Pad::kTPCLSB);
466  id.SetField(half,
467  ND::GeomId::Def::TPC::Pad::kHalfMSB,
468  ND::GeomId::Def::TPC::Pad::kHalfLSB);
469  id.SetField(mm,
470  ND::GeomId::Def::TPC::Pad::kMMegaMSB,
471  ND::GeomId::Def::TPC::Pad::kMMegaLSB);
472  id.SetField(1,
473  ND::GeomId::Def::TPC::Pad::kPadFlagMSB,
474  ND::GeomId::Def::TPC::Pad::kPadFlagLSB);
475  id.SetField(pad,
476  ND::GeomId::Def::TPC::Pad::kPadMSB,
477  ND::GeomId::Def::TPC::Pad::kPadLSB);
478  return id;
479 }
480 
481 int ND::GeomId::TPC::GetMicroMegaTPC(TGeometryId i) {
482  return ND::GeomId::TPC::GetPadTPC(i);
483 }
484 
485 int ND::GeomId::TPC::GetPadTPC(TGeometryId i) {
486  TSettableGeometryId id(i);
487  if (ND::GeomId::Def::kTPC
488  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
489  ND::GeomId::Def::kDetectorIdLSB)) return -1;
490  if (ND::GeomId::Def::TPC::kPad
491  != id.GetField(ND::GeomId::Def::TPC::kSeqIdMSB,
492  ND::GeomId::Def::TPC::kSeqIdLSB)) return -1;
493  return id.GetField(ND::GeomId::Def::TPC::Pad::kTPCMSB,
494  ND::GeomId::Def::TPC::Pad::kTPCLSB);
495 }
496 
497 int ND::GeomId::TPC::GetMicroMegaHalf(TGeometryId i) {
498  return ND::GeomId::TPC::GetPadHalf(i);
499 }
500 
501 int ND::GeomId::TPC::GetPadHalf(TGeometryId i) {
502  TSettableGeometryId id(i);
503  if (ND::GeomId::Def::kTPC
504  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
505  ND::GeomId::Def::kDetectorIdLSB)) return -1;
506  if (ND::GeomId::Def::TPC::kPad
507  != id.GetField(ND::GeomId::Def::TPC::kSeqIdMSB,
508  ND::GeomId::Def::TPC::kSeqIdLSB)) return -1;
509  return id.GetField(ND::GeomId::Def::TPC::Pad::kHalfMSB,
510  ND::GeomId::Def::TPC::Pad::kHalfLSB);
511 }
512 
513 int ND::GeomId::TPC::GetMicroMegaNumber(TGeometryId i) {
514  return ND::GeomId::TPC::GetPadMicroMega(i);
515 }
516 
517 int ND::GeomId::TPC::GetPadMicroMega(TGeometryId i) {
518  TSettableGeometryId id(i);
519  if (ND::GeomId::Def::kTPC
520  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
521  ND::GeomId::Def::kDetectorIdLSB)) return -1;
522  if (ND::GeomId::Def::TPC::kPad
523  != id.GetField(ND::GeomId::Def::TPC::kSeqIdMSB,
524  ND::GeomId::Def::TPC::kSeqIdLSB)) return -1;
525  return id.GetField(ND::GeomId::Def::TPC::Pad::kMMegaMSB,
526  ND::GeomId::Def::TPC::Pad::kMMegaLSB);
527 }
528 
529 int ND::GeomId::TPC::GetPadNumber(TGeometryId i) {
530  TSettableGeometryId id(i);
531  if (ND::GeomId::Def::kTPC
532  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
533  ND::GeomId::Def::kDetectorIdLSB)) return -1;
534  if (ND::GeomId::Def::TPC::kPad
535  != id.GetField(ND::GeomId::Def::TPC::kSeqIdMSB,
536  ND::GeomId::Def::TPC::kSeqIdLSB)) return -1;
537  if (1 != id.GetField(ND::GeomId::Def::TPC::Pad::kPadFlagMSB,
538  ND::GeomId::Def::TPC::Pad::kPadFlagLSB)) return -1;
539  return id.GetField(ND::GeomId::Def::TPC::Pad::kPadMSB,
540  ND::GeomId::Def::TPC::Pad::kPadLSB);
541 }
542 
543 bool ND::GeomId::TPC::IsPad(ND::TGeometryId i) {
544  TSettableGeometryId id(i);
545  int detector = id.GetField(ND::GeomId::Def::kDetectorIdMSB,
546  ND::GeomId::Def::kDetectorIdLSB);
547  if (detector != ND::GeomId::Def::kTPC) return false;
548  /// This may be a pad.
549  int seqId = id.GetField(ND::GeomId::Def::TPC::kSeqIdMSB,
550  ND::GeomId::Def::TPC::kSeqIdLSB);
551  if (seqId != ND::GeomId::Def::TPC::kPad) return false;
552  int padFlag = id.GetField(ND::GeomId::Def::TPC::Pad::kPadFlagMSB,
553  ND::GeomId::Def::TPC::Pad::kPadFlagLSB);
554  if (!padFlag) return false;
555  return true;
556 }
557 
558 ////////////////////////////////////////////////////////////
559 // FGD Geometry Identifiers.
560 ////////////////////////////////////////////////////////////
561 bool ND::GeomId::FGD::IsFGD(TGeometryId i) {
562  TSettableGeometryId id(i);
563  return (ND::GeomId::Def::kFGD
564  == id.GetField(ND::GeomId::Def::kDetectorIdMSB,
565  ND::GeomId::Def::kDetectorIdLSB));
566 }
567 
568 ND::TGeometryId ND::GeomId::FGD::FGD(int fgd) {
569  TSettableGeometryId id;
570  if (fgd<0 || fgd>1) {
571  ND280Warn("FGD out of range [0,1]: " << fgd);
572  return id;
573  }
574  id.SetField(ND::GeomId::Def::kFGD,
575  ND::GeomId::Def::kDetectorIdMSB,
576  ND::GeomId::Def::kDetectorIdLSB);
577  id.SetField(ND::GeomId::Def::FGD::kGlobal,
578  ND::GeomId::Def::FGD::kSeqIdMSB,
579  ND::GeomId::Def::FGD::kSeqIdLSB);
580  id.SetField(ND::GeomId::Def::FGD::Global::kFGD,
581  ND::GeomId::Def::FGD::Global::kSeqIdMSB,
582  ND::GeomId::Def::FGD::Global::kSeqIdLSB);
583  id.SetField(fgd,
584  ND::GeomId::Def::FGD::Global::kFieldMSB,
585  ND::GeomId::Def::FGD::Global::kFieldLSB);
586  return id;
587 }
588 
589 ND::TGeometryId ND::GeomId::FGD::FGD1() {return FGD(0);}
590 ND::TGeometryId ND::GeomId::FGD::FGD2() {return FGD(1);}
591 
592 int ND::GeomId::FGD::GetFGD(TGeometryId i) {
593  TSettableGeometryId id(i);
594  if (ND::GeomId::Def::kFGD
595  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
596  ND::GeomId::Def::kDetectorIdLSB)) return -1;
597  if (ND::GeomId::Def::FGD::kGlobal
598  != id.GetField(ND::GeomId::Def::FGD::kSeqIdMSB,
599  ND::GeomId::Def::FGD::kSeqIdLSB)) return -1;
600  if (ND::GeomId::Def::FGD::Global::kFGD
601  != id.GetField(ND::GeomId::Def::FGD::Global::kSeqIdMSB,
602  ND::GeomId::Def::FGD::Global::kSeqIdLSB)) return -1;
603  return id.GetField(ND::GeomId::Def::FGD::Global::kFieldMSB,
604  ND::GeomId::Def::FGD::Global::kFieldLSB);
605 }
606 
607 bool ND::GeomId::FGD::IsFGD1(TGeometryId id) {
608  return (ND::GeomId::FGD::GetFGD(id) == 0);
609 }
610 
611 bool ND::GeomId::FGD::IsFGD2(TGeometryId id) {
612  return (ND::GeomId::FGD::GetFGD(id) == 1);
613 }
614 
615 ND::TGeometryId ND::GeomId::FGD::Bar(int fgd,
616  int module,
617  int layer,
618  int bar) {
619  TSettableGeometryId id;
620  if (fgd<0 || fgd>1) {
621  ND280Warn("FGD out of range [0,1]: " << fgd);
622  return id;
623  }
624  id.SetField(ND::GeomId::Def::kFGD,
625  ND::GeomId::Def::kDetectorIdMSB,
626  ND::GeomId::Def::kDetectorIdLSB);
627  id.SetField(ND::GeomId::Def::FGD::kScintillator,
628  ND::GeomId::Def::FGD::kSeqIdMSB,
629  ND::GeomId::Def::FGD::kSeqIdLSB);
630  id.SetField(fgd,
631  ND::GeomId::Def::FGD::Bar::kModuleMSB,
632  ND::GeomId::Def::FGD::Bar::kModuleLSB);
633  id.SetField(module,
634  ND::GeomId::Def::FGD::Bar::kPlaneMSB,
635  ND::GeomId::Def::FGD::Bar::kPlaneLSB);
636  id.SetField(layer,
637  ND::GeomId::Def::FGD::Bar::kLayerMSB,
638  ND::GeomId::Def::FGD::Bar::kLayerLSB);
639  id.SetField(bar,
640  ND::GeomId::Def::FGD::Bar::kBarMSB,
641  ND::GeomId::Def::FGD::Bar::kBarLSB);
642  return id;
643 }
644 
645 int ND::GeomId::FGD::GetBarFGD(TGeometryId i) {
646  TSettableGeometryId id(i);
647  if (ND::GeomId::Def::kFGD
648  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
649  ND::GeomId::Def::kDetectorIdLSB)) return -1;
650  if (ND::GeomId::Def::FGD::kScintillator
651  != id.GetField(ND::GeomId::Def::FGD::kSeqIdMSB,
652  ND::GeomId::Def::FGD::kSeqIdLSB)) return -1;
653  return id.GetField(ND::GeomId::Def::FGD::Bar::kModuleMSB,
654  ND::GeomId::Def::FGD::Bar::kModuleLSB);
655 }
656 
657 int ND::GeomId::FGD::GetBarModule(TGeometryId i) {
658  TSettableGeometryId id(i);
659  if (ND::GeomId::Def::kFGD
660  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
661  ND::GeomId::Def::kDetectorIdLSB)) return -1;
662  if (ND::GeomId::Def::FGD::kScintillator
663  != id.GetField(ND::GeomId::Def::FGD::kSeqIdMSB,
664  ND::GeomId::Def::FGD::kSeqIdLSB)) return -1;
665  return id.GetField(ND::GeomId::Def::FGD::Bar::kPlaneMSB,
666  ND::GeomId::Def::FGD::Bar::kPlaneLSB);
667 }
668 
669 int ND::GeomId::FGD::GetBarLayer(TGeometryId i) {
670  TSettableGeometryId id(i);
671  if (ND::GeomId::Def::kFGD
672  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
673  ND::GeomId::Def::kDetectorIdLSB)) return -1;
674  if (ND::GeomId::Def::FGD::kScintillator
675  != id.GetField(ND::GeomId::Def::FGD::kSeqIdMSB,
676  ND::GeomId::Def::FGD::kSeqIdLSB)) return -1;
677  return id.GetField(ND::GeomId::Def::FGD::Bar::kLayerMSB,
678  ND::GeomId::Def::FGD::Bar::kLayerLSB);
679 }
680 
681 int ND::GeomId::FGD::GetBarNumber(TGeometryId i) {
682  TSettableGeometryId id(i);
683  if (ND::GeomId::Def::kFGD
684  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
685  ND::GeomId::Def::kDetectorIdLSB)) return -1;
686  if (ND::GeomId::Def::FGD::kScintillator
687  != id.GetField(ND::GeomId::Def::FGD::kSeqIdMSB,
688  ND::GeomId::Def::FGD::kSeqIdLSB)) return -1;
689  return id.GetField(ND::GeomId::Def::FGD::Bar::kBarMSB,
690  ND::GeomId::Def::FGD::Bar::kBarLSB);
691 }
692 
693 ND::TGeometryId ND::GeomId::FGD::Target(int target) {
694  TSettableGeometryId id;
695  id.SetField(ND::GeomId::Def::kFGD,
696  ND::GeomId::Def::kDetectorIdMSB,
697  ND::GeomId::Def::kDetectorIdLSB);
698  id.SetField(ND::GeomId::Def::FGD::kGlobal,
699  ND::GeomId::Def::FGD::kSeqIdMSB,
700  ND::GeomId::Def::FGD::kSeqIdLSB);
701  id.SetField(ND::GeomId::Def::FGD::Global::kTarget,
702  ND::GeomId::Def::FGD::Global::kSeqIdMSB,
703  ND::GeomId::Def::FGD::Global::kSeqIdLSB);
704  id.SetField(target,
705  ND::GeomId::Def::FGD::Global::kFieldMSB,
706  ND::GeomId::Def::FGD::Global::kFieldLSB);
707  return id;
708 }
709 
710 int ND::GeomId::FGD::GetTarget(TGeometryId i) {
711  TSettableGeometryId id(i);
712  if (ND::GeomId::Def::kFGD
713  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
714  ND::GeomId::Def::kDetectorIdLSB)) return -1;
715  if (ND::GeomId::Def::FGD::kGlobal
716  != id.GetField(ND::GeomId::Def::FGD::kSeqIdMSB,
717  ND::GeomId::Def::FGD::kSeqIdLSB)) return -1;
718  if (ND::GeomId::Def::FGD::Global::kTarget
719  != id.GetField(ND::GeomId::Def::FGD::Global::kSeqIdMSB,
720  ND::GeomId::Def::FGD::Global::kSeqIdLSB)) return -1;
721  return id.GetField(ND::GeomId::Def::FGD::Global::kFieldMSB,
722  ND::GeomId::Def::FGD::Global::kFieldLSB);
723 }
724 
725 ND::TGeometryId ND::GeomId::FGD::Layer(int fgd,
726  int module,
727  int layer) {
728  TSettableGeometryId id;
729  if (fgd<0 || fgd>1) {
730  ND280Warn("FGD out of range [0,1]: " << fgd);
731  return id;
732  }
733  id.SetField(ND::GeomId::Def::kFGD,
734  ND::GeomId::Def::kDetectorIdMSB,
735  ND::GeomId::Def::kDetectorIdLSB);
736  id.SetField(ND::GeomId::Def::FGD::kGlobal,
737  ND::GeomId::Def::FGD::kSeqIdMSB,
738  ND::GeomId::Def::FGD::kSeqIdLSB);
739  id.SetField(ND::GeomId::Def::FGD::Global::kLayer,
740  ND::GeomId::Def::FGD::Global::kSeqIdMSB,
741  ND::GeomId::Def::FGD::Global::kSeqIdLSB);
742  id.SetField(fgd,
743  ND::GeomId::Def::FGD::Global::Layer::kFGDMSB,
744  ND::GeomId::Def::FGD::Global::Layer::kFGDLSB);
745  id.SetField(layer,
746  ND::GeomId::Def::FGD::Global::Layer::kLayerMSB,
747  ND::GeomId::Def::FGD::Global::Layer::kLayerLSB);
748  id.SetField(module,
749  ND::GeomId::Def::FGD::Global::Layer::kModuleMSB,
750  ND::GeomId::Def::FGD::Global::Layer::kModuleLSB);
751  return id;
752 }
753 
754 int ND::GeomId::FGD::GetLayerFGD(TGeometryId i) {
755  TSettableGeometryId id(i);
756  if (ND::GeomId::Def::kFGD
757  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
758  ND::GeomId::Def::kDetectorIdLSB)) return -1;
759  if (ND::GeomId::Def::FGD::kGlobal
760  != id.GetField(ND::GeomId::Def::FGD::kSeqIdMSB,
761  ND::GeomId::Def::FGD::kSeqIdLSB)) return -1;
762  if (ND::GeomId::Def::FGD::Global::kLayer
763  != id.GetField(ND::GeomId::Def::FGD::Global::kSeqIdMSB,
764  ND::GeomId::Def::FGD::Global::kSeqIdLSB)) return -1;
765  return id.GetField(ND::GeomId::Def::FGD::Global::Layer::kFGDMSB,
766  ND::GeomId::Def::FGD::Global::Layer::kFGDLSB);
767 }
768 
769 int ND::GeomId::FGD::GetLayerModule(TGeometryId i) {
770  TSettableGeometryId id(i);
771  if (ND::GeomId::Def::kFGD
772  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
773  ND::GeomId::Def::kDetectorIdLSB)) return -1;
774  if (ND::GeomId::Def::FGD::kGlobal
775  != id.GetField(ND::GeomId::Def::FGD::kSeqIdMSB,
776  ND::GeomId::Def::FGD::kSeqIdLSB)) return -1;
777  if (ND::GeomId::Def::FGD::Global::kLayer
778  != id.GetField(ND::GeomId::Def::FGD::Global::kSeqIdMSB,
779  ND::GeomId::Def::FGD::Global::kSeqIdLSB)) return -1;
780  return id.GetField(ND::GeomId::Def::FGD::Global::Layer::kModuleMSB,
781  ND::GeomId::Def::FGD::Global::Layer::kModuleLSB);
782 }
783 
784 int ND::GeomId::FGD::GetLayerNumber(TGeometryId i) {
785  TSettableGeometryId id(i);
786  if (ND::GeomId::Def::kFGD
787  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
788  ND::GeomId::Def::kDetectorIdLSB)) return -1;
789  if (ND::GeomId::Def::FGD::kGlobal
790  != id.GetField(ND::GeomId::Def::FGD::kSeqIdMSB,
791  ND::GeomId::Def::FGD::kSeqIdLSB)) return -1;
792  if (ND::GeomId::Def::FGD::Global::kLayer
793  != id.GetField(ND::GeomId::Def::FGD::Global::kSeqIdMSB,
794  ND::GeomId::Def::FGD::Global::kSeqIdLSB)) return -1;
795  return id.GetField(ND::GeomId::Def::FGD::Global::Layer::kLayerMSB,
796  ND::GeomId::Def::FGD::Global::Layer::kLayerLSB);
797 }
798 
799 ////////////////////////////////////////////////////////////
800 // ECal Geometry Identifiers.
801 ////////////////////////////////////////////////////////////
802 
803 ND::TGeometryId ND::GeomId::ECal::Module(int ecal, int clam, int module) {
804  TSettableGeometryId id;
805  if (ecal == ND::GeomId::Def::kDSECal) {
806  if (clam != 0) {
807  ND280Error("Downstream ECal with invalid clam: " << clam);
808  return id;
809  }
810  if (module != 0) {
811  ND280Error("Downstream ECal with invalid module: " << module);
812  return id;
813  }
814  }
815  else if (ecal == ND::GeomId::Def::kTECal) {
816  if (clam != 0 && clam != 1) {
817  ND280Error("Tracker ECal with invalid clam: " << clam);
818  return id;
819  }
820  if (module<0 || 2<module) {
821  ND280Error("Tracker ECal with invalid module: " << module);
822  return id;
823  }
824  }
825  else if (ecal == ND::GeomId::Def::kPECal) {
826  if (clam != 0 && clam != 1) {
827  ND280Error("P0D ECal with invalid clam: " << clam);
828  return id;
829  }
830  if (module<0 || 2<module) {
831  ND280Error("P0D ECal with invalid module: " << module);
832  return id;
833  }
834  }
835  id.SetField(ecal,
836  ND::GeomId::Def::kDetectorIdMSB,
837  ND::GeomId::Def::kDetectorIdLSB);
838  id.SetField(ND::GeomId::Def::ECal::kGlobal,
839  ND::GeomId::Def::ECal::kSeqIdMSB,
840  ND::GeomId::Def::ECal::kSeqIdLSB);
841  id.SetField(module*2+clam,
842  ND::GeomId::Def::ECal::kModuleMSB,
843  ND::GeomId::Def::ECal::kModuleLSB);
844 
845  id.SetField(ND::GeomId::Def::ECal::Global::kECal,
846  ND::GeomId::Def::ECal::Global::kSeqIdMSB,
847  ND::GeomId::Def::ECal::Global::kSeqIdLSB);
848  id.SetField(0,
849  ND::GeomId::Def::ECal::Global::kFieldMSB,
850  ND::GeomId::Def::ECal::Global::kFieldLSB);
851 
852  return id;
853 }
854 
855 ND::TGeometryId ND::GeomId::ECal::Container(int ecal, int clam, int module) {
856  TSettableGeometryId id;
857  if (ecal == ND::GeomId::Def::kTECal) {
858  if (clam != 0 && clam != 1) {
859  ND280Error("Tracker ECal container with invalid clam: " << clam);
860  return id;
861  }
862  if (module<0 || 2<module) {
863  ND280Error("Tracker ECal container with invalid module: " << module);
864  return id;
865  }
866  }
867  else if (ecal == ND::GeomId::Def::kPECal) {
868  if (clam != 0 && clam != 1) {
869  ND280Error("P0D ECal container with invalid clam: " << clam);
870  return id;
871  }
872  if (module<0 || 2<module) {
873  ND280Error("P0D ECal container with invalid module: " << module);
874  return id;
875  }
876  }
877  else {
878  ND280Error("Container for non-contained ECal: " << ecal);
879  return id;
880  }
881 
882  id.SetField(ecal,
883  ND::GeomId::Def::kDetectorIdMSB,
884  ND::GeomId::Def::kDetectorIdLSB);
885  id.SetField(ND::GeomId::Def::ECal::kGlobal,
886  ND::GeomId::Def::ECal::kSeqIdMSB,
887  ND::GeomId::Def::ECal::kSeqIdLSB);
888  id.SetField(module*2+clam,
889  ND::GeomId::Def::ECal::kModuleMSB,
890  ND::GeomId::Def::ECal::kModuleLSB);
891 
892  id.SetField(ND::GeomId::Def::ECal::Global::kContainer,
893  ND::GeomId::Def::ECal::Global::kSeqIdMSB,
894  ND::GeomId::Def::ECal::Global::kSeqIdLSB);
895  id.SetField(0,
896  ND::GeomId::Def::ECal::Global::kFieldMSB,
897  ND::GeomId::Def::ECal::Global::kFieldLSB);
898 
899  return id;
900 }
901 
902 
903 int ND::GeomId::ECal::GetModuleECal(TGeometryId i) {
904  TSettableGeometryId id(i);
905  int det = ND::GeomId::GetSubSystem(id);
906  if (det == ND::GeomId::Def::kDSECal) return det;
907  else if (det == ND::GeomId::Def::kTECal) return det;
908  else if (det == ND::GeomId::Def::kPECal) return det;
909  else return -1;
910 }
911 
912 int ND::GeomId::ECal::GetModuleNumber(TGeometryId i) {
913  TSettableGeometryId id(i);
914  if (ND::GeomId::ECal::GetModuleECal(i)<0) return -1;
915  if (ND::GeomId::Def::ECal::kGlobal
916  != id.GetField(ND::GeomId::Def::ECal::kSeqIdMSB,
917  ND::GeomId::Def::ECal::kSeqIdLSB)) return -1;
918  if (ND::GeomId::Def::ECal::Global::kECal
919  != id.GetField(ND::GeomId::Def::ECal::Global::kSeqIdMSB,
920  ND::GeomId::Def::ECal::Global::kSeqIdLSB)) return -1;
921 
922  return id.GetField(ND::GeomId::Def::ECal::kModuleMSB,
923  ND::GeomId::Def::ECal::kModuleLSB) / 2;
924 }
925 
926 int ND::GeomId::ECal::GetModuleClam(TGeometryId i) {
927  TSettableGeometryId id(i);
928  if (ND::GeomId::ECal::GetModuleECal(i)<0) return -1;
929  if (ND::GeomId::Def::ECal::kGlobal
930  != id.GetField(ND::GeomId::Def::ECal::kSeqIdMSB,
931  ND::GeomId::Def::ECal::kSeqIdLSB)) return -1;
932  if (ND::GeomId::Def::ECal::Global::kECal
933  != id.GetField(ND::GeomId::Def::ECal::Global::kSeqIdMSB,
934  ND::GeomId::Def::ECal::Global::kSeqIdLSB)) return -1;
935 
936  return id.GetField(ND::GeomId::Def::ECal::kModuleMSB,
937  ND::GeomId::Def::ECal::kModuleLSB) % 2;
938 }
939 
940 ND::TGeometryId ND::GeomId::ECal::Layer(int ecal, int clam,
941  int module, int layer) {
942  TSettableGeometryId id;
943  if (ecal == ND::GeomId::Def::kDSECal) {
944  if (clam != 0) {
945  ND280Error("Downstream ECal with invalid clam: " << clam);
946  return id;
947  }
948  if (module != 0) {
949  ND280Error("Downstream ECal with invalid module: " << module);
950  return id;
951  }
952  }
953  else if (ecal == ND::GeomId::Def::kTECal) {
954  if (clam != 0 && clam != 1) {
955  ND280Error("Tracker ECal with invalid clam: " << clam);
956  return id;
957  }
958  if (module<0 || 2<module) {
959  ND280Error("Tracker ECal with invalid module: " << module);
960  return id;
961  }
962  }
963  else if (ecal == ND::GeomId::Def::kPECal) {
964  if (clam != 0 && clam != 1) {
965  ND280Error("P0D ECal with invalid clam: " << clam);
966  return id;
967  }
968  if (module<0 || 2<module) {
969  ND280Error("P0D ECal with invalid module: " << module);
970  return id;
971  }
972  }
973  id.SetField(ecal,
974  ND::GeomId::Def::kDetectorIdMSB,
975  ND::GeomId::Def::kDetectorIdLSB);
976  id.SetField(ND::GeomId::Def::ECal::kGlobal,
977  ND::GeomId::Def::ECal::kSeqIdMSB,
978  ND::GeomId::Def::ECal::kSeqIdLSB);
979  id.SetField(module*2+clam,
980  ND::GeomId::Def::ECal::kModuleMSB,
981  ND::GeomId::Def::ECal::kModuleLSB);
982  id.SetField(ND::GeomId::Def::ECal::Global::kLayer,
983  ND::GeomId::Def::ECal::Global::kSeqIdMSB,
984  ND::GeomId::Def::ECal::Global::kSeqIdLSB);
985  id.SetField(layer,
986  ND::GeomId::Def::ECal::Global::kFieldMSB,
987  ND::GeomId::Def::ECal::Global::kFieldLSB);
988  return id;
989 }
990 
991 int ND::GeomId::ECal::GetLayerECal(TGeometryId i) {
992  TSettableGeometryId id(i);
993  int det = ND::GeomId::GetSubSystem(id);
994  if (det == ND::GeomId::Def::kDSECal) return det;
995  else if (det == ND::GeomId::Def::kTECal) return det;
996  else if (det == ND::GeomId::Def::kPECal) return det;
997  else return -1;
998 }
999 
1000 int ND::GeomId::ECal::GetLayerClam(TGeometryId i) {
1001  TSettableGeometryId id(i);
1002  if (ND::GeomId::ECal::GetLayerECal(i)<0) return -1;
1003  if (ND::GeomId::Def::ECal::kGlobal
1004  != id.GetField(ND::GeomId::Def::ECal::kSeqIdMSB,
1005  ND::GeomId::Def::ECal::kSeqIdLSB)) return -1;
1006  if (ND::GeomId::Def::ECal::Global::kLayer
1007  != id.GetField(ND::GeomId::Def::ECal::Global::kSeqIdMSB,
1008  ND::GeomId::Def::ECal::Global::kSeqIdLSB)) return -1;
1009 
1010  return id.GetField(ND::GeomId::Def::ECal::kModuleMSB,
1011  ND::GeomId::Def::ECal::kModuleLSB) % 2;
1012 }
1013 
1014 int ND::GeomId::ECal::GetLayerModule(TGeometryId i) {
1015  TSettableGeometryId id(i);
1016  if (ND::GeomId::ECal::GetLayerECal(i)<0) return -1;
1017  if (ND::GeomId::Def::ECal::kGlobal
1018  != id.GetField(ND::GeomId::Def::ECal::kSeqIdMSB,
1019  ND::GeomId::Def::ECal::kSeqIdLSB)) return -1;
1020  if (ND::GeomId::Def::ECal::Global::kLayer
1021  != id.GetField(ND::GeomId::Def::ECal::Global::kSeqIdMSB,
1022  ND::GeomId::Def::ECal::Global::kSeqIdLSB)) return -1;
1023 
1024  return id.GetField(ND::GeomId::Def::ECal::kModuleMSB,
1025  ND::GeomId::Def::ECal::kModuleLSB) / 2;
1026 }
1027 
1028 int ND::GeomId::ECal::GetLayerNumber(TGeometryId i) {
1029  TSettableGeometryId id(i);
1030  if (ND::GeomId::ECal::GetLayerECal(i)<0) return -1;
1031  if (ND::GeomId::Def::ECal::kGlobal
1032  != id.GetField(ND::GeomId::Def::ECal::kSeqIdMSB,
1033  ND::GeomId::Def::ECal::kSeqIdLSB)) return -1;
1034  if (ND::GeomId::Def::ECal::Global::kLayer
1035  != id.GetField(ND::GeomId::Def::ECal::Global::kSeqIdMSB,
1036  ND::GeomId::Def::ECal::Global::kSeqIdLSB)) return -1;
1037 
1038  return id.GetField(ND::GeomId::Def::ECal::Global::kFieldMSB,
1039  ND::GeomId::Def::ECal::Global::kFieldLSB);
1040 }
1041 
1042 ND::TGeometryId ND::GeomId::ECal::Radiator(int ecal, int clam,
1043  int module, int rad) {
1044  TSettableGeometryId id;
1045  if (ecal == ND::GeomId::Def::kDSECal) {
1046  if (clam != 0) {
1047  ND280Error("Downstream ECal with invalid clam: " << clam);
1048  return id;
1049  }
1050  if (module != 0) {
1051  ND280Error("Downstream ECal with invalid module: " << module);
1052  return id;
1053  }
1054  }
1055  else if (ecal == ND::GeomId::Def::kTECal) {
1056  if (clam != 0 && clam != 1) {
1057  ND280Error("Tracker ECal with invalid clam: " << clam);
1058  return id;
1059  }
1060  if (module<0 || 2<module) {
1061  ND280Error("Tracker ECal with invalid module: " << module);
1062  return id;
1063  }
1064  }
1065  else if (ecal == ND::GeomId::Def::kPECal) {
1066  if (clam != 0 && clam != 1) {
1067  ND280Error("P0D ECal with invalid clam: " << clam);
1068  return id;
1069  }
1070  if (module<0 || 2<module) {
1071  ND280Error("P0D ECal with invalid module: " << module);
1072  return id;
1073  }
1074  }
1075  id.SetField(ecal,
1076  ND::GeomId::Def::kDetectorIdMSB,
1077  ND::GeomId::Def::kDetectorIdLSB);
1078  id.SetField(ND::GeomId::Def::ECal::kGlobal,
1079  ND::GeomId::Def::ECal::kSeqIdMSB,
1080  ND::GeomId::Def::ECal::kSeqIdLSB);
1081  id.SetField(module*2+clam,
1082  ND::GeomId::Def::ECal::kModuleMSB,
1083  ND::GeomId::Def::ECal::kModuleLSB);
1084  id.SetField(ND::GeomId::Def::ECal::Global::kRadiator,
1085  ND::GeomId::Def::ECal::Global::kSeqIdMSB,
1086  ND::GeomId::Def::ECal::Global::kSeqIdLSB);
1087  id.SetField(rad,
1088  ND::GeomId::Def::ECal::Global::kFieldMSB,
1089  ND::GeomId::Def::ECal::Global::kFieldLSB);
1090  return id;
1091 }
1092 
1093 int ND::GeomId::ECal::GetRadiatorECal(TGeometryId i) {
1094  TSettableGeometryId id(i);
1095  int det = ND::GeomId::GetSubSystem(id);
1096  if (det == ND::GeomId::Def::kDSECal) return det;
1097  else if (det == ND::GeomId::Def::kTECal) return det;
1098  else if (det == ND::GeomId::Def::kPECal) return det;
1099  else return -1;
1100 }
1101 
1102 int ND::GeomId::ECal::GetRadiatorClam(TGeometryId i) {
1103  TSettableGeometryId id(i);
1104  if (ND::GeomId::ECal::GetRadiatorECal(i)<0) return -1;
1105  if (ND::GeomId::Def::ECal::kGlobal
1106  != id.GetField(ND::GeomId::Def::ECal::kSeqIdMSB,
1107  ND::GeomId::Def::ECal::kSeqIdLSB)) return -1;
1108  if (ND::GeomId::Def::ECal::Global::kRadiator
1109  != id.GetField(ND::GeomId::Def::ECal::Global::kSeqIdMSB,
1110  ND::GeomId::Def::ECal::Global::kSeqIdLSB)) return -1;
1111 
1112  return id.GetField(ND::GeomId::Def::ECal::kModuleMSB,
1113  ND::GeomId::Def::ECal::kModuleLSB) % 2;
1114 }
1115 
1116 int ND::GeomId::ECal::GetRadiatorModule(TGeometryId i) {
1117  TSettableGeometryId id(i);
1118  if (ND::GeomId::ECal::GetRadiatorECal(i)<0) return -1;
1119  if (ND::GeomId::Def::ECal::kGlobal
1120  != id.GetField(ND::GeomId::Def::ECal::kSeqIdMSB,
1121  ND::GeomId::Def::ECal::kSeqIdLSB)) return -1;
1122  if (ND::GeomId::Def::ECal::Global::kRadiator
1123  != id.GetField(ND::GeomId::Def::ECal::Global::kSeqIdMSB,
1124  ND::GeomId::Def::ECal::Global::kSeqIdLSB)) return -1;
1125 
1126  return id.GetField(ND::GeomId::Def::ECal::kModuleMSB,
1127  ND::GeomId::Def::ECal::kModuleLSB) / 2;
1128 }
1129 
1130 int ND::GeomId::ECal::GetRadiatorNumber(TGeometryId i) {
1131  TSettableGeometryId id(i);
1132  if (ND::GeomId::ECal::GetRadiatorECal(i)<0) return -1;
1133  if (ND::GeomId::Def::ECal::kGlobal
1134  != id.GetField(ND::GeomId::Def::ECal::kSeqIdMSB,
1135  ND::GeomId::Def::ECal::kSeqIdLSB)) return -1;
1136  if (ND::GeomId::Def::ECal::Global::kRadiator
1137  != id.GetField(ND::GeomId::Def::ECal::Global::kSeqIdMSB,
1138  ND::GeomId::Def::ECal::Global::kSeqIdLSB)) return -1;
1139 
1140  return id.GetField(ND::GeomId::Def::ECal::Global::kFieldMSB,
1141  ND::GeomId::Def::ECal::Global::kFieldLSB);
1142 }
1143 
1144 ND::TGeometryId ND::GeomId::ECal::Bar(int ecal, int clam,
1145  int module, int layer, int bar) {
1146  TSettableGeometryId id;
1147  if (ecal == ND::GeomId::Def::kDSECal) {
1148  if (clam != 0) {
1149  ND280Error("Downstream ECal with invalid clam: " << clam);
1150  return id;
1151  }
1152  if (module != 0) {
1153  ND280Error("Downstream ECal with invalid module: " << module);
1154  return id;
1155  }
1156  }
1157  else if (ecal == ND::GeomId::Def::kTECal) {
1158  if (clam != 0 && clam != 1) {
1159  ND280Error("Tracker ECal with invalid clam: " << clam);
1160  return id;
1161  }
1162  if (module<0 || 2<module) {
1163  ND280Error("Tracker ECal with invalid module: " << module);
1164  return id;
1165  }
1166  }
1167  else if (ecal == ND::GeomId::Def::kPECal) {
1168  if (clam != 0 && clam != 1) {
1169  ND280Error("P0D ECal with invalid clam: " << clam);
1170  return id;
1171  }
1172  if (module<0 || 2<module) {
1173  ND280Error("P0D ECal with invalid module: " << module);
1174  return id;
1175  }
1176  }
1177  id.SetField(ecal,
1178  ND::GeomId::Def::kDetectorIdMSB,
1179  ND::GeomId::Def::kDetectorIdLSB);
1180  id.SetField(ND::GeomId::Def::ECal::kScintillator,
1181  ND::GeomId::Def::ECal::kSeqIdMSB,
1182  ND::GeomId::Def::ECal::kSeqIdLSB);
1183  id.SetField(module*2+clam,
1184  ND::GeomId::Def::ECal::kModuleMSB,
1185  ND::GeomId::Def::ECal::kModuleLSB);
1186  id.SetField(layer,
1187  ND::GeomId::Def::ECal::Bar::kLayerMSB,
1188  ND::GeomId::Def::ECal::Bar::kLayerLSB);
1189  id.SetField(bar,
1190  ND::GeomId::Def::ECal::Bar::kBarMSB,
1191  ND::GeomId::Def::ECal::Bar::kBarLSB);
1192  return id;
1193 }
1194 
1195 int ND::GeomId::ECal::GetBarECal(TGeometryId i) {
1196  TSettableGeometryId id(i);
1197  int det = ND::GeomId::GetSubSystem(id);
1198  if (det == ND::GeomId::Def::kDSECal) return det;
1199  else if (det == ND::GeomId::Def::kTECal) return det;
1200  else if (det == ND::GeomId::Def::kPECal) return det;
1201  else return -1;
1202 }
1203 
1204 int ND::GeomId::ECal::GetBarClam(TGeometryId i) {
1205  TSettableGeometryId id(i);
1206  if (ND::GeomId::ECal::GetBarECal(i)<0) return -1;
1207  if (ND::GeomId::Def::ECal::kScintillator
1208  != id.GetField(ND::GeomId::Def::ECal::kSeqIdMSB,
1209  ND::GeomId::Def::ECal::kSeqIdLSB)) return -1;
1210 
1211  return id.GetField(ND::GeomId::Def::ECal::kModuleMSB,
1212  ND::GeomId::Def::ECal::kModuleLSB) % 2;
1213 }
1214 
1215 int ND::GeomId::ECal::GetBarModule(TGeometryId i) {
1216  TSettableGeometryId id(i);
1217  if (ND::GeomId::ECal::GetBarECal(i)<0) return -1;
1218  if (ND::GeomId::Def::ECal::kScintillator
1219  != id.GetField(ND::GeomId::Def::ECal::kSeqIdMSB,
1220  ND::GeomId::Def::ECal::kSeqIdLSB)) return -1;
1221 
1222  return id.GetField(ND::GeomId::Def::ECal::kModuleMSB,
1223  ND::GeomId::Def::ECal::kModuleLSB) / 2;
1224 }
1225 
1226 int ND::GeomId::ECal::GetBarLayer(TGeometryId i) {
1227  TSettableGeometryId id(i);
1228  if (ND::GeomId::ECal::GetBarECal(i)<0) return -1;
1229  if (ND::GeomId::Def::ECal::kScintillator
1230  != id.GetField(ND::GeomId::Def::ECal::kSeqIdMSB,
1231  ND::GeomId::Def::ECal::kSeqIdLSB)) return -1;
1232 
1233  return id.GetField(ND::GeomId::Def::ECal::Bar::kLayerMSB,
1234  ND::GeomId::Def::ECal::Bar::kLayerLSB);
1235 }
1236 
1237 int ND::GeomId::ECal::GetBarNumber(TGeometryId i) {
1238  TSettableGeometryId id(i);
1239  if (ND::GeomId::ECal::GetBarECal(i)<0) return -1;
1240  if (ND::GeomId::Def::ECal::kScintillator
1241  != id.GetField(ND::GeomId::Def::ECal::kSeqIdMSB,
1242  ND::GeomId::Def::ECal::kSeqIdLSB)) return -1;
1243 
1244  return id.GetField(ND::GeomId::Def::ECal::Bar::kBarMSB,
1245  ND::GeomId::Def::ECal::Bar::kBarLSB);
1246 }
1247 
1248 ////////////////////////////////////////////////////////////
1249 // Downstream ECal Geometry Identifiers.
1250 ////////////////////////////////////////////////////////////
1251 bool ND::GeomId::DSECal::IsDSECal(TGeometryId i) {
1252  TSettableGeometryId id(i);
1253  return (ND::GeomId::Def::kDSECal
1254  == id.GetField(ND::GeomId::Def::kDetectorIdMSB,
1255  ND::GeomId::Def::kDetectorIdLSB));
1256 }
1257 
1258 ND::TGeometryId ND::GeomId::DSECal::Detector() {
1259  return ND::GeomId::ECal::Module(ND::GeomId::Def::kDSECal,0,0);
1260 }
1261 
1262 ND::TGeometryId ND::GeomId::DSECal::Layer(int layer) {
1263  return ND::GeomId::ECal::Layer(ND::GeomId::Def::kDSECal,0,0,layer);
1264 }
1265 
1266 int ND::GeomId::DSECal::GetLayer(TGeometryId i) {
1267  if (!ND::GeomId::DSECal::IsDSECal(i)) return -1;
1268  return ND::GeomId::ECal::GetLayerNumber(i);
1269 }
1270 
1271 ND::TGeometryId ND::GeomId::DSECal::Radiator(int rad) {
1272  return ND::GeomId::ECal::Radiator(ND::GeomId::Def::kDSECal,0,0,rad);
1273 }
1274 
1275 int ND::GeomId::DSECal::GetRadiator(TGeometryId i) {
1276  if (!ND::GeomId::DSECal::IsDSECal(i)) return -1;
1277  return ND::GeomId::ECal::GetRadiatorNumber(i);
1278 }
1279 
1280 ND::TGeometryId ND::GeomId::DSECal::Bar(int layer, int bar) {
1281  return ND::GeomId::ECal::Bar(ND::GeomId::Def::kDSECal,
1282  0,0,layer,bar);
1283 }
1284 
1285 int ND::GeomId::DSECal::GetBarLayer(TGeometryId i) {
1286  if (!ND::GeomId::DSECal::IsDSECal(i)) return -1;
1287  return ND::GeomId::ECal::GetBarLayer(i);
1288 }
1289 
1290 int ND::GeomId::DSECal::GetBarNumber(TGeometryId i) {
1291  if (!ND::GeomId::DSECal::IsDSECal(i)) return -1;
1292  return ND::GeomId::ECal::GetBarNumber(i);
1293 }
1294 
1295 ////////////////////////////////////////////////////////////
1296 // Tracker ECal Geometry Identifiers.
1297 ////////////////////////////////////////////////////////////
1298 bool ND::GeomId::TECal::IsTECal(TGeometryId i) {
1299  TSettableGeometryId id(i);
1300  return (ND::GeomId::Def::kTECal
1301  == id.GetField(ND::GeomId::Def::kDetectorIdMSB,
1302  ND::GeomId::Def::kDetectorIdLSB));
1303 }
1304 
1305 ND::TGeometryId ND::GeomId::TECal::Module(int clam, int module) {
1306  return ND::GeomId::ECal::Module(ND::GeomId::Def::kTECal,
1307  clam,module);
1308 }
1309 
1310 ND::TGeometryId ND::GeomId::TECal::Layer(int clam, int module, int layer) {
1311  return ND::GeomId::ECal::Layer(ND::GeomId::Def::kTECal,
1312  clam,module,layer);
1313 }
1314 
1315 ND::TGeometryId ND::GeomId::TECal::Radiator(int clam, int module, int rad) {
1316  return ND::GeomId::ECal::Radiator(ND::GeomId::Def::kTECal,
1317  clam,module,rad);
1318 }
1319 
1320 ND::TGeometryId ND::GeomId::TECal::Bar(int clam, int mod, int layer, int bar) {
1321  return ND::GeomId::ECal::Bar(ND::GeomId::Def::kTECal,
1322  clam,mod,layer,bar);
1323 }
1324 
1325 int ND::GeomId::TECal::GetModuleNumber(TGeometryId i) {
1326  if (!ND::GeomId::TECal::IsTECal(i)) return -1;
1327  return ND::GeomId::ECal::GetModuleNumber(i);
1328 }
1329 
1330 int ND::GeomId::TECal::GetModuleClam(TGeometryId i) {
1331  if (!ND::GeomId::TECal::IsTECal(i)) return -1;
1332  return ND::GeomId::ECal::GetModuleClam(i);
1333 }
1334 
1335 int ND::GeomId::TECal::GetLayerNumber(TGeometryId i) {
1336  if (!ND::GeomId::TECal::IsTECal(i)) return -1;
1337  return ND::GeomId::ECal::GetLayerNumber(i);
1338 }
1339 
1340 int ND::GeomId::TECal::GetLayerClam(TGeometryId i) {
1341  if (!ND::GeomId::TECal::IsTECal(i)) return -1;
1342  return ND::GeomId::ECal::GetLayerClam(i);
1343 }
1344 
1345 int ND::GeomId::TECal::GetLayerModule(TGeometryId i) {
1346  if (!ND::GeomId::TECal::IsTECal(i)) return -1;
1347  return ND::GeomId::ECal::GetLayerModule(i);
1348 }
1349 
1350 int ND::GeomId::TECal::GetRadiatorNumber(TGeometryId i) {
1351  if (!ND::GeomId::TECal::IsTECal(i)) return -1;
1352  return ND::GeomId::ECal::GetRadiatorNumber(i);
1353 }
1354 
1355 int ND::GeomId::TECal::GetRadiatorClam(TGeometryId i) {
1356  if (!ND::GeomId::TECal::IsTECal(i)) return -1;
1357  return ND::GeomId::ECal::GetRadiatorClam(i);
1358 }
1359 
1360 int ND::GeomId::TECal::GetRadiatorModule(TGeometryId i) {
1361  if (!ND::GeomId::TECal::IsTECal(i)) return -1;
1362  return ND::GeomId::ECal::GetRadiatorModule(i);
1363 }
1364 
1365 int ND::GeomId::TECal::GetBarNumber(TGeometryId i) {
1366  if (!ND::GeomId::TECal::IsTECal(i)) return -1;
1367  return ND::GeomId::ECal::GetBarNumber(i);
1368 }
1369 
1370 int ND::GeomId::TECal::GetBarLayer(TGeometryId i) {
1371  if (!ND::GeomId::TECal::IsTECal(i)) return -1;
1372  return ND::GeomId::ECal::GetBarLayer(i);
1373 }
1374 
1375 int ND::GeomId::TECal::GetBarClam(TGeometryId i) {
1376  if (!ND::GeomId::TECal::IsTECal(i)) return -1;
1377  return ND::GeomId::ECal::GetBarClam(i);
1378 }
1379 
1380 int ND::GeomId::TECal::GetBarModule(TGeometryId i) {
1381  if (!ND::GeomId::TECal::IsTECal(i)) return -1;
1382  return ND::GeomId::ECal::GetBarModule(i);
1383 }
1384 
1385 ////////////////////////////////////////////////////////////
1386 // P0D ECal Geometry Identifiers.
1387 ////////////////////////////////////////////////////////////
1388 bool ND::GeomId::PECal::IsPECal(TGeometryId i) {
1389  TSettableGeometryId id(i);
1390  return (ND::GeomId::Def::kPECal
1391  == id.GetField(ND::GeomId::Def::kDetectorIdMSB,
1392  ND::GeomId::Def::kDetectorIdLSB));
1393 }
1394 
1395 ND::TGeometryId ND::GeomId::PECal::Module(int clam, int module) {
1396  return ND::GeomId::ECal::Module(ND::GeomId::Def::kPECal,
1397  clam,module);
1398 }
1399 
1400 ND::TGeometryId ND::GeomId::PECal::Layer(int clam, int module, int layer) {
1401  return ND::GeomId::ECal::Layer(ND::GeomId::Def::kPECal,
1402  clam,module,layer);
1403 }
1404 
1405 ND::TGeometryId ND::GeomId::PECal::Radiator(int clam, int module, int rad) {
1406  return ND::GeomId::ECal::Radiator(ND::GeomId::Def::kPECal,
1407  clam,module,rad);
1408 }
1409 
1410 ND::TGeometryId ND::GeomId::PECal::Bar(int clam, int mod, int layer, int bar) {
1411  return ND::GeomId::ECal::Bar(ND::GeomId::Def::kPECal,
1412  clam,mod,layer,bar);
1413 }
1414 
1415 int ND::GeomId::PECal::GetModuleNumber(TGeometryId i) {
1416  if (!ND::GeomId::PECal::IsPECal(i)) return -1;
1417  return ND::GeomId::ECal::GetModuleNumber(i);
1418 }
1419 
1420 int ND::GeomId::PECal::GetModuleClam(TGeometryId i) {
1421  if (!ND::GeomId::PECal::IsPECal(i)) return -1;
1422  return ND::GeomId::ECal::GetModuleClam(i);
1423 }
1424 
1425 int ND::GeomId::PECal::GetLayerNumber(TGeometryId i) {
1426  if (!ND::GeomId::PECal::IsPECal(i)) return -1;
1427  return ND::GeomId::ECal::GetLayerNumber(i);
1428 }
1429 
1430 int ND::GeomId::PECal::GetLayerClam(TGeometryId i) {
1431  if (!ND::GeomId::PECal::IsPECal(i)) return -1;
1432  return ND::GeomId::ECal::GetLayerClam(i);
1433 }
1434 
1435 int ND::GeomId::PECal::GetLayerModule(TGeometryId i) {
1436  if (!ND::GeomId::PECal::IsPECal(i)) return -1;
1437  return ND::GeomId::ECal::GetLayerModule(i);
1438 }
1439 
1440 int ND::GeomId::PECal::GetRadiatorNumber(TGeometryId i) {
1441  if (!ND::GeomId::PECal::IsPECal(i)) return -1;
1442  return ND::GeomId::ECal::GetRadiatorNumber(i);
1443 }
1444 
1445 int ND::GeomId::PECal::GetRadiatorClam(TGeometryId i) {
1446  if (!ND::GeomId::PECal::IsPECal(i)) return -1;
1447  return ND::GeomId::ECal::GetRadiatorClam(i);
1448 }
1449 
1450 int ND::GeomId::PECal::GetRadiatorModule(TGeometryId i) {
1451  if (!ND::GeomId::PECal::IsPECal(i)) return -1;
1452  return ND::GeomId::ECal::GetRadiatorModule(i);
1453 }
1454 
1455 int ND::GeomId::PECal::GetBarNumber(TGeometryId i) {
1456  if (!ND::GeomId::PECal::IsPECal(i)) return -1;
1457  return ND::GeomId::ECal::GetBarNumber(i);
1458 }
1459 
1460 int ND::GeomId::PECal::GetBarLayer(TGeometryId i) {
1461  if (!ND::GeomId::PECal::IsPECal(i)) return -1;
1462  return ND::GeomId::ECal::GetBarLayer(i);
1463 }
1464 
1465 int ND::GeomId::PECal::GetBarClam(TGeometryId i) {
1466  if (!ND::GeomId::PECal::IsPECal(i)) return -1;
1467  return ND::GeomId::ECal::GetBarClam(i);
1468 }
1469 
1470 int ND::GeomId::PECal::GetBarModule(TGeometryId i) {
1471  if (!ND::GeomId::PECal::IsPECal(i)) return -1;
1472  return ND::GeomId::ECal::GetBarModule(i);
1473 }
1474 
1475 ////////////////////////////////////////////////////////////
1476 // SMRD Geometry Identifiers.
1477 ////////////////////////////////////////////////////////////
1478 bool ND::GeomId::SMRD::IsSMRD(TGeometryId i) {
1479  TSettableGeometryId id(i);
1480  return (ND::GeomId::Def::kSMRD
1481  == id.GetField(ND::GeomId::Def::kDetectorIdMSB,
1482  ND::GeomId::Def::kDetectorIdLSB));
1483 }
1484 
1485 ND::TGeometryId ND::GeomId::SMRD::Module(int clam, int yoke,
1486  int layer, int slot) {
1487  TSettableGeometryId id;
1488  id.SetField(ND::GeomId::Def::kSMRD,
1489  ND::GeomId::Def::kDetectorIdMSB,
1490  ND::GeomId::Def::kDetectorIdLSB);
1491  id.SetField(ND::GeomId::Def::SMRD::kGlobal,
1492  ND::GeomId::Def::SMRD::kSeqIdMSB,
1493  ND::GeomId::Def::SMRD::kSeqIdLSB);
1494  id.SetField(ND::GeomId::Def::SMRD::Global::kModule,
1495  ND::GeomId::Def::SMRD::Global::kSeqIdMSB,
1496  ND::GeomId::Def::SMRD::Global::kSeqIdLSB);
1497  id.SetField(clam,
1498  ND::GeomId::Def::SMRD::Global::Module::kClamMSB,
1499  ND::GeomId::Def::SMRD::Global::Module::kClamLSB);
1500  id.SetField(yoke,
1501  ND::GeomId::Def::SMRD::Global::Module::kYokeMSB,
1502  ND::GeomId::Def::SMRD::Global::Module::kYokeLSB);
1503  id.SetField(layer,
1504  ND::GeomId::Def::SMRD::Global::Module::kLayerMSB,
1505  ND::GeomId::Def::SMRD::Global::Module::kLayerLSB);
1506  id.SetField(slot,
1507  ND::GeomId::Def::SMRD::Global::Module::kSlotMSB,
1508  ND::GeomId::Def::SMRD::Global::Module::kSlotLSB);
1509  return id;
1510 }
1511 
1512 int ND::GeomId::SMRD::GetModuleClam(TGeometryId i) {
1513  TSettableGeometryId id(i);
1514  if (ND::GeomId::Def::kSMRD
1515  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
1516  ND::GeomId::Def::kDetectorIdLSB)) return -1;
1517  if (ND::GeomId::Def::SMRD::kGlobal
1518  != id.GetField(ND::GeomId::Def::SMRD::kSeqIdMSB,
1519  ND::GeomId::Def::SMRD::kSeqIdLSB)) return -1;
1520  if (ND::GeomId::Def::SMRD::Global::kModule
1521  != id.GetField(ND::GeomId::Def::SMRD::Global::kSeqIdMSB,
1522  ND::GeomId::Def::SMRD::Global::kSeqIdLSB)) return -1;
1523 
1524  return id.GetField(ND::GeomId::Def::SMRD::Global::Module::kClamMSB,
1525  ND::GeomId::Def::SMRD::Global::Module::kClamLSB);
1526 }
1527 
1528 int ND::GeomId::SMRD::GetModuleYoke(TGeometryId i) {
1529  TSettableGeometryId id(i);
1530  if (ND::GeomId::Def::kSMRD
1531  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
1532  ND::GeomId::Def::kDetectorIdLSB)) return -1;
1533  if (ND::GeomId::Def::SMRD::kGlobal
1534  != id.GetField(ND::GeomId::Def::SMRD::kSeqIdMSB,
1535  ND::GeomId::Def::SMRD::kSeqIdLSB)) return -1;
1536  if (ND::GeomId::Def::SMRD::Global::kModule
1537  != id.GetField(ND::GeomId::Def::SMRD::Global::kSeqIdMSB,
1538  ND::GeomId::Def::SMRD::Global::kSeqIdLSB)) return -1;
1539 
1540 
1541  return id.GetField(ND::GeomId::Def::SMRD::Global::Module::kYokeMSB,
1542  ND::GeomId::Def::SMRD::Global::Module::kYokeLSB);
1543 }
1544 
1545 int ND::GeomId::SMRD::GetModuleLayer(TGeometryId i) {
1546  TSettableGeometryId id(i);
1547  if (ND::GeomId::Def::kSMRD
1548  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
1549  ND::GeomId::Def::kDetectorIdLSB)) return -1;
1550  if (ND::GeomId::Def::SMRD::kGlobal
1551  != id.GetField(ND::GeomId::Def::SMRD::kSeqIdMSB,
1552  ND::GeomId::Def::SMRD::kSeqIdLSB)) return -1;
1553  if (ND::GeomId::Def::SMRD::Global::kModule
1554  != id.GetField(ND::GeomId::Def::SMRD::Global::kSeqIdMSB,
1555  ND::GeomId::Def::SMRD::Global::kSeqIdLSB)) return -1;
1556 
1557 
1558  return id.GetField(ND::GeomId::Def::SMRD::Global::Module::kLayerMSB,
1559  ND::GeomId::Def::SMRD::Global::Module::kLayerLSB);
1560 }
1561 
1562 int ND::GeomId::SMRD::GetModuleSlot(TGeometryId i) {
1563  TSettableGeometryId id(i);
1564  if (ND::GeomId::Def::kSMRD
1565  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
1566  ND::GeomId::Def::kDetectorIdLSB)) return -1;
1567  if (ND::GeomId::Def::SMRD::kGlobal
1568  != id.GetField(ND::GeomId::Def::SMRD::kSeqIdMSB,
1569  ND::GeomId::Def::SMRD::kSeqIdLSB)) return -1;
1570  if (ND::GeomId::Def::SMRD::Global::kModule
1571  != id.GetField(ND::GeomId::Def::SMRD::Global::kSeqIdMSB,
1572  ND::GeomId::Def::SMRD::Global::kSeqIdLSB)) return -1;
1573 
1574 
1575  return id.GetField(ND::GeomId::Def::SMRD::Global::Module::kSlotMSB,
1576  ND::GeomId::Def::SMRD::Global::Module::kSlotLSB);
1577 }
1578 
1579 ND::TGeometryId ND::GeomId::SMRD::Bar(int clam, int yoke,
1580  int layer, int slot, int bar) {
1581  TSettableGeometryId id;
1582  id.SetField(ND::GeomId::Def::kSMRD,
1583  ND::GeomId::Def::kDetectorIdMSB,
1584  ND::GeomId::Def::kDetectorIdLSB);
1585  id.SetField(ND::GeomId::Def::SMRD::kScintillator,
1586  ND::GeomId::Def::SMRD::kSeqIdMSB,
1587  ND::GeomId::Def::SMRD::kSeqIdLSB);
1588  id.SetField(clam,
1589  ND::GeomId::Def::SMRD::Bar::kClamMSB,
1590  ND::GeomId::Def::SMRD::Bar::kClamLSB);
1591  id.SetField(yoke,
1592  ND::GeomId::Def::SMRD::Bar::kYokeMSB,
1593  ND::GeomId::Def::SMRD::Bar::kYokeLSB);
1594  id.SetField(layer,
1595  ND::GeomId::Def::SMRD::Bar::kLayerMSB,
1596  ND::GeomId::Def::SMRD::Bar::kLayerLSB);
1597  id.SetField(slot,
1598  ND::GeomId::Def::SMRD::Bar::kSlotMSB,
1599  ND::GeomId::Def::SMRD::Bar::kSlotLSB);
1600  id.SetField(bar,
1601  ND::GeomId::Def::SMRD::Bar::kBarMSB,
1602  ND::GeomId::Def::SMRD::Bar::kBarLSB);
1603  return id;
1604 }
1605 
1606 int ND::GeomId::SMRD::GetBarClam(TGeometryId i) {
1607  TSettableGeometryId id(i);
1608  if (ND::GeomId::Def::kSMRD
1609  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
1610  ND::GeomId::Def::kDetectorIdLSB)) return -1;
1611  if (ND::GeomId::Def::SMRD::kScintillator
1612  != id.GetField(ND::GeomId::Def::SMRD::kSeqIdMSB,
1613  ND::GeomId::Def::SMRD::kSeqIdLSB)) return -1;
1614 
1615  return id.GetField(ND::GeomId::Def::SMRD::Bar::kClamMSB,
1616  ND::GeomId::Def::SMRD::Bar::kClamLSB);
1617 }
1618 
1619 int ND::GeomId::SMRD::GetBarYoke(TGeometryId i) {
1620  TSettableGeometryId id(i);
1621  if (ND::GeomId::Def::kSMRD
1622  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
1623  ND::GeomId::Def::kDetectorIdLSB)) return -1;
1624  if (ND::GeomId::Def::SMRD::kScintillator
1625  != id.GetField(ND::GeomId::Def::SMRD::kSeqIdMSB,
1626  ND::GeomId::Def::SMRD::kSeqIdLSB)) return -1;
1627 
1628  return id.GetField(ND::GeomId::Def::SMRD::Bar::kYokeMSB,
1629  ND::GeomId::Def::SMRD::Bar::kYokeLSB);
1630 }
1631 
1632 int ND::GeomId::SMRD::GetBarLayer(TGeometryId i) {
1633  TSettableGeometryId id(i);
1634  if (ND::GeomId::Def::kSMRD
1635  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
1636  ND::GeomId::Def::kDetectorIdLSB)) return -1;
1637  if (ND::GeomId::Def::SMRD::kScintillator
1638  != id.GetField(ND::GeomId::Def::SMRD::kSeqIdMSB,
1639  ND::GeomId::Def::SMRD::kSeqIdLSB)) return -1;
1640 
1641  return id.GetField(ND::GeomId::Def::SMRD::Bar::kLayerMSB,
1642  ND::GeomId::Def::SMRD::Bar::kLayerLSB);
1643 }
1644 
1645 int ND::GeomId::SMRD::GetBarSlot(TGeometryId i) {
1646  TSettableGeometryId id(i);
1647  if (ND::GeomId::Def::kSMRD
1648  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
1649  ND::GeomId::Def::kDetectorIdLSB)) return -1;
1650  if (ND::GeomId::Def::SMRD::kScintillator
1651  != id.GetField(ND::GeomId::Def::SMRD::kSeqIdMSB,
1652  ND::GeomId::Def::SMRD::kSeqIdLSB)) return -1;
1653 
1654  return id.GetField(ND::GeomId::Def::SMRD::Bar::kSlotMSB,
1655  ND::GeomId::Def::SMRD::Bar::kSlotLSB);
1656 }
1657 
1658 int ND::GeomId::SMRD::GetBarNumber(TGeometryId i) {
1659  TSettableGeometryId id(i);
1660  if (ND::GeomId::Def::kSMRD
1661  != id.GetField(ND::GeomId::Def::kDetectorIdMSB,
1662  ND::GeomId::Def::kDetectorIdLSB)) return -1;
1663  if (ND::GeomId::Def::SMRD::kScintillator
1664  != id.GetField(ND::GeomId::Def::SMRD::kSeqIdMSB,
1665  ND::GeomId::Def::SMRD::kSeqIdLSB)) return -1;
1666 
1667  return id.GetField(ND::GeomId::Def::SMRD::Bar::kBarMSB,
1668  ND::GeomId::Def::SMRD::Bar::kBarLSB);
1669 }
1670 
1671 ////////////////////////////////////////////////////////////
1672 // INGRID Geometry Identifiers.
1673 ////////////////////////////////////////////////////////////
1674 
1675 ND::TGeometryId ND::GeomId::INGRID::Scintillator(
1676  int objType, int obj, int trk, int proj, int scinti){
1677  TSettableGeometryId id;
1678  id.SetField(ND::GeomId::Def::kINGRID,
1679  ND::GeomId::Def::kDetectorIdMSB,
1680  ND::GeomId::Def::kDetectorIdLSB);
1681  id.SetField(objType,
1682  ND::GeomId::Def::INGRID::kIngridObjIDMSB,
1683  ND::GeomId::Def::INGRID::kIngridObjIDLSB);
1684  id.SetField(obj,
1685  ND::GeomId::Def::INGRID::kIngridModNumMSB,
1686  ND::GeomId::Def::INGRID::kIngridModNumLSB);
1687  id.SetField(trk,
1688  ND::GeomId::Def::INGRID::kIngridTrkNumMSB,
1689  ND::GeomId::Def::INGRID::kIngridTrkNumLSB);
1690  id.SetField(proj,
1691  ND::GeomId::Def::INGRID::kIngridProjMSB,
1692  ND::GeomId::Def::INGRID::kIngridProjLSB);
1693  id.SetField(scinti,
1694  ND::GeomId::Def::INGRID::kIngridSciNumMSB,
1695  ND::GeomId::Def::INGRID::kIngridSciNumLSB);
1696  return id;
1697 }
1698 
1699 ND::TGeometryId ND::GeomId::INGRID::ModScintillator(
1700  int obj, int trk, int proj, int scinti){
1701  return ND::GeomId::INGRID::Scintillator(
1702  ND::GeomId::Def::INGRID::kIngridModule, obj, trk, proj, scinti);
1703 }
1704 
1705 ND::TGeometryId ND::GeomId::INGRID::VertVetoScintillator(int obj, int scinti){
1706  return ND::GeomId::INGRID::Scintillator(
1707  ND::GeomId::Def::INGRID::kIngridVeto, obj, 0,
1708  ND::GeomId::Def::INGRID::kVertical, scinti);
1709 }
1710 
1711 ND::TGeometryId ND::GeomId::INGRID::HorzVetoScintillator(int obj, int scinti){
1712  return ND::GeomId::INGRID::Scintillator(
1713  ND::GeomId::Def::INGRID::kIngridVeto, obj, 0,
1714  ND::GeomId::Def::INGRID::kHorizontal, scinti);
1715 }
ND::GeomId::FGD::IsFGD1
bool IsFGD1(TGeometryId id)
Definition: ND280GeomId.cxx:607
ND::GeomId::INGRID::VertVetoScintillator
TGeometryId VertVetoScintillator(int veto, int scinti)
The geometry of a vertical veto scintillator.
Definition: ND280GeomId.cxx:1705
ND::GeomId::P0D::Bar
TGeometryId Bar(int p0dule, int layer, int bar)
Definition: ND280GeomId.cxx:262
ND::GeomId::TPC::Module
TGeometryId Module(int mod)
Definition: ND280GeomId.cxx:334
ND::GeomId::ECal::Container
TGeometryId Container(int ecal, int clam, int module)
Definition: ND280GeomId.cxx:855
ND::GeomId::PECal::GetRadiatorClam
int GetRadiatorClam(TGeometryId id)
Definition: ND280GeomId.cxx:1445
ND::GeomId::PECal::Layer
TGeometryId Layer(int clam, int module, int layer)
The layer within the ecal.
Definition: ND280GeomId.cxx:1400
ND::GeomId::P0D::ECalRadiator
TGeometryId ECalRadiator(int radiator)
Definition: ND280GeomId.cxx:198
ND::GeomId::P0D::GetTargetRadiator
int GetTargetRadiator(TGeometryId id)
Definition: ND280GeomId.cxx:247
ND::GeomId::DSECal::GetLayer
int GetLayer(TGeometryId id)
Definition: ND280GeomId.cxx:1266
ND::GeomId::TECal::GetBarClam
int GetBarClam(TGeometryId id)
Definition: ND280GeomId.cxx:1375
ND::GeomId::TPC::MicroMega
TGeometryId MicroMega(int tpc, int half, int mm)
Definition: ND280GeomId.cxx:386
ND::GeomId::Def::kDetectorIdLSB
const GeomIdBit kDetectorIdLSB
The least significant bit of the detector identifier.
Definition: ND280GeomIdDef.hxx:50
ND::GeomId::TECal::Module
TGeometryId Module(int clam, int module)
The geometry identifier for the tracker ecal modules.
Definition: ND280GeomId.cxx:1305
ND::GeomId::Def::INGRID::kIngridSciNumMSB
const GeomIdBit kIngridSciNumMSB
Definition: ND280GeomIdDef.hxx:759
ND::GeomId::PECal::Radiator
TGeometryId Radiator(int clam, int module, int radiator)
The geometry identifier for the radiators.
Definition: ND280GeomId.cxx:1405
ND::GeomId::P0D::GetSuperP0Dule
int GetSuperP0Dule(TGeometryId id)
Definition: ND280GeomId.cxx:118
ND::GeomId::Def::INGRID::kIngridObjIDMSB
const GeomIdBit kIngridObjIDMSB
Definition: ND280GeomIdDef.hxx:723
ND::GeomId::PECal::Bar
TGeometryId Bar(int clam, int module, int layer, int bar)
The geometry identifiers for the scintillators.
Definition: ND280GeomId.cxx:1410
ND::GeomId::P0D::GetP0Dule
int GetP0Dule(TGeometryId id)
Definition: ND280GeomId.cxx:151
ND::GeomId::SMRD::GetModuleClam
int GetModuleClam(TGeometryId id)
Definition: ND280GeomId.cxx:1512
ND::GeomId::FGD::GetTarget
int GetTarget(TGeometryId id)
Definition: ND280GeomId.cxx:710
ND::GeomId::SMRD::Bar
TGeometryId Bar(int clam, int yoke, int layer, int slot, int bar)
The geometry identifer for a SMRD bar.
Definition: ND280GeomId.cxx:1579
ND::GeomId::DSECal::Bar
TGeometryId Bar(int layer, int bar)
The geometry identifiers for DSECal scintillators.
Definition: ND280GeomId.cxx:1280
ND::GeomId::P0D::GetTarget
int GetTarget(TGeometryId id)
Definition: ND280GeomId.cxx:183
ND::GeomId::Def::FGD::kSeqIdMSB
const GeomIdBit kSeqIdMSB
Definition: ND280GeomIdDef.hxx:301
ND::GeomId::FGD::Target
TGeometryId Target(int i)
The geometry identifier for the targets.
Definition: ND280GeomId.cxx:693
ND::GeomId::SMRD::Module
TGeometryId Module(int clam, int yoke, int layer, int slot)
The geometry identifier for a SMRD module.
Definition: ND280GeomId.cxx:1485
ND::GeomId::P0D::GetBarLayer
int GetBarLayer(TGeometryId id)
Definition: ND280GeomId.cxx:299
ND::GeomId::TECal::Radiator
TGeometryId Radiator(int clam, int module, int radiator)
The geometry identifier for the radiators.
Definition: ND280GeomId.cxx:1315
ND::GeomId::TECal::GetRadiatorClam
int GetRadiatorClam(TGeometryId id)
Definition: ND280GeomId.cxx:1355
ND::GeomId::ECal::Module
TGeometryId Module(int ecal, int clam, int module)
Definition: ND280GeomId.cxx:803
ND::GeomId::FGD::Layer
TGeometryId Layer(int fgd, int module, int layer)
Definition: ND280GeomId.cxx:725
ND::GeomId::TPC::GetMicroMegaTPC
int GetMicroMegaTPC(TGeometryId id)
Definition: ND280GeomId.cxx:481
ND::GeomId::Def::FGD::Global::kFieldMSB
const GeomIdBit kFieldMSB
The sequence value for this type of global volume.
Definition: ND280GeomIdDef.hxx:334
ND::GeomId::FGD::GetLayerFGD
int GetLayerFGD(TGeometryId id)
Definition: ND280GeomId.cxx:754
ND::GeomId::TPC::IsTPC1
bool IsTPC1(TGeometryId id)
Check to see if the id is for a TPCn.
Definition: ND280GeomId.cxx:374
ND::GeomId::FGD::FGD
TGeometryId FGD(int mod)
Definition: ND280GeomId.cxx:568
ND::GeomId::Def::INGRID::kIngridTrkNumMSB
const GeomIdBit kIngridTrkNumMSB
Definition: ND280GeomIdDef.hxx:737
ND::GeomId::Def::TPC::Pad::kHalfMSB
const GeomIdBit kHalfMSB
The half of the drift volume.
Definition: ND280GeomIdDef.hxx:257
ND::GeomId::ECal::Radiator
TGeometryId Radiator(int ecal, int clam, int module, int radiator)
The geometry identifier for the radiators.
Definition: ND280GeomId.cxx:1042
ND::GeomId::FGD::FGD1
TGeometryId FGD1()
The geometry identifier for FGD1.
Definition: ND280GeomId.cxx:589
ND::GeomId::Empty
TGeometryId Empty()
Return the geometry id indication None.
Definition: ND280GeomId.cxx:63
ND::GeomId::P0D::GetBarNumber
int GetBarNumber(TGeometryId id)
Definition: ND280GeomId.cxx:311
ND::GeomId::ECal::Bar
TGeometryId Bar(int ecal, int clam, int module, int layer, int bar)
The geometry identifiers for the scintillators.
Definition: ND280GeomId.cxx:1144
ND::GeomId::PECal::GetLayerClam
int GetLayerClam(TGeometryId id)
Definition: ND280GeomId.cxx:1430
ND::GeomId::Def::kDetectorIdMSB
const GeomIdBit kDetectorIdMSB
Definition: ND280GeomIdDef.hxx:47
ND::GeomId::INGRID::Scintillator
TGeometryId Scintillator(int objType, int obj, int trk, int proj, int scinti)
Definition: ND280GeomId.cxx:1675
ND::GeomId::SMRD::IsSMRD
bool IsSMRD(TGeometryId id)
Test if this id is for the SMRD.
Definition: ND280GeomId.cxx:1478
ND::GeomId::TECal::GetModuleClam
int GetModuleClam(TGeometryId id)
Access the information for an ECal Module.
Definition: ND280GeomId.cxx:1330
ND::GeomId::Def::P0D::Global::kP0D
A special case for one of a kind global volumes.
Definition: ND280GeomIdDef.hxx:137
ND::GeomId::Def::P0D::Global::kTargetRadiator
One of the brass target radiators.
Definition: ND280GeomIdDef.hxx:147
ND::GeomId::Def::P0D::Global::kTarget
One of the water target modules.
Definition: ND280GeomIdDef.hxx:143
ND::GeomId::Def::ECal::kModuleMSB
const GeomIdBit kModuleMSB
Definition: ND280GeomIdDef.hxx:473
ND::GeomId::DSECal::IsDSECal
bool IsDSECal(TGeometryId id)
Check if the id corresponds to the DSECal.
Definition: ND280GeomId.cxx:1251
ND::GeomId::Def::P0D::Global::kFieldMSB
const GeomIdBit kFieldMSB
The sequence value for this type of global volume.
Definition: ND280GeomIdDef.hxx:151
ND::GeomId::PECal::IsPECal
bool IsPECal(TGeometryId id)
Check if the id corresponds to the PECal.
Definition: ND280GeomId.cxx:1388
ND::GeomId::TPC::GetPadTPC
int GetPadTPC(TGeometryId id)
Definition: ND280GeomId.cxx:485
ND::GeomId::DSECal::GetRadiator
int GetRadiator(TGeometryId id)
Definition: ND280GeomId.cxx:1275
ND::GeomId::ECal::GetModuleECal
int GetModuleECal(TGeometryId id)
Access the information for an ECal Module.
Definition: ND280GeomId.cxx:903
si_syst::field
const TVector3 field
Constants.
Definition: SecondaryInteractionSystematic.hxx:33
ND::GeomId::DSECal::Layer
TGeometryId Layer(int layer)
The layer within the ecal.
Definition: ND280GeomId.cxx:1262
ND::GeomId::TPC::IsTPC
bool IsTPC(TGeometryId id)
Check to see if the id is for a TPC.
Definition: ND280GeomId.cxx:327
ND::GeomId::P0D::Target
TGeometryId Target(int target)
Definition: ND280GeomId.cxx:166
ND::GeomId::TPC::Pad
TGeometryId Pad(int tpc, int half, int mm, int pad)
Definition: ND280GeomId.cxx:439
ND::GeomId::PECal::Module
TGeometryId Module(int clam, int module)
The geometry identifier for the tracker ecal modules.
Definition: ND280GeomId.cxx:1395
ND::GeomId::P0D::TargetRadiator
TGeometryId TargetRadiator(int radiator)
Definition: ND280GeomId.cxx:230
ND::GeomId::FGD::GetBarFGD
int GetBarFGD(TGeometryId id)
Definition: ND280GeomId.cxx:645
ND::GeomId::TECal::GetLayerClam
int GetLayerClam(TGeometryId id)
Definition: ND280GeomId.cxx:1340
ND::GeomId::DSECal::GetBarLayer
int GetBarLayer(TGeometryId id)
Definition: ND280GeomId.cxx:1285
ND::GeomId::INGRID::ModScintillator
TGeometryId ModScintillator(int mod, int trk, int proj, int scinti)
The geometry of a module scintillator.
Definition: ND280GeomId.cxx:1699
ND::GeomId::TPC::GetModule
int GetModule(TGeometryId id)
Definition: ND280GeomId.cxx:355
ND::GeomId::PECal::GetBarClam
int GetBarClam(TGeometryId id)
Definition: ND280GeomId.cxx:1465
ND::GeomId::Def::INGRID::kIngridProjMSB
const GeomIdBit kIngridProjMSB
Definition: ND280GeomIdDef.hxx:752
ND::GeomId::DSECal::Detector
TGeometryId Detector()
The geometry identifier for the downstream ecal detector.
Definition: ND280GeomId.cxx:1258
ND::GeomId::PECal::GetModuleClam
int GetModuleClam(TGeometryId id)
Access the information for an ECal Module.
Definition: ND280GeomId.cxx:1420
ND::GeomId::INGRID::HorzVetoScintillator
TGeometryId HorzVetoScintillator(int veto, int scinti)
The geometry of an horizontal veto scintillator.
Definition: ND280GeomId.cxx:1711
ND::GeomId::P0D::SuperP0Dule
TGeometryId SuperP0Dule(int sP0Dule)
Definition: ND280GeomId.cxx:101
ND::GeomId::SMRD::GetBarClam
int GetBarClam(TGeometryId id)
Definition: ND280GeomId.cxx:1606
ND::GeomId::FGD::GetFGD
int GetFGD(TGeometryId id)
Get the FGD number. This returns -1 if the id is invalid.
Definition: ND280GeomId.cxx:592
ND::GeomId::TPC::TPC2
TGeometryId TPC2()
The geometry identifier for TPC2.
Definition: ND280GeomId.cxx:371
ND::GeomId::P0D::Detector
TGeometryId Detector()
Define a geometry identifier to the entire P0D detector.
Definition: ND280GeomId.cxx:84
ND::GeomId::Def::INGRID::kIngridModNumMSB
const GeomIdBit kIngridModNumMSB
Definition: ND280GeomIdDef.hxx:730
ND::GeomId::Def::ECal::Global::kLayer
One of the ECal scintillator layers.
Definition: ND280GeomIdDef.hxx:502
ND::GeomId::Def::SMRD::Global::kModule
The sequence for all of the modules in the SMRD.
Definition: ND280GeomIdDef.hxx:593
ND::GeomId::TPC::IsPad
bool IsPad(TGeometryId id)
Check if a TGeometryId is a PAD.
Definition: ND280GeomId.cxx:543
ND::GeomId::ECal::Layer
TGeometryId Layer(int ecal, int clam, int module, int layer)
The layer within the ecal.
Definition: ND280GeomId.cxx:940
ND::GeomId::P0D::GetECalRadiator
int GetECalRadiator(TGeometryId id)
Definition: ND280GeomId.cxx:215
ND::GeomId::FGD::Bar
TGeometryId Bar(int fgd, int module, int layer, int bar)
Definition: ND280GeomId.cxx:615
ND::GeomId::ECal::GetRadiatorECal
int GetRadiatorECal(TGeometryId id)
Definition: ND280GeomId.cxx:1093
ND::GeomId::TECal::IsTECal
bool IsTECal(TGeometryId id)
Check if the id corresponds to the TECal.
Definition: ND280GeomId.cxx:1298
ND::GeomId::Def::ECal::kSeqIdMSB
const GeomIdBit kSeqIdMSB
Definition: ND280GeomIdDef.hxx:449
ND::GeomId::TPC::TPC3
TGeometryId TPC3()
The geometry identifier for TPC3.
Definition: ND280GeomId.cxx:372
ND::GeomId::Def::ECal::Global::kFieldMSB
const GeomIdBit kFieldMSB
The sequence value for this type of global volume.
Definition: ND280GeomIdDef.hxx:513
ND::GeomId::TPC::IsMicroMega
bool IsMicroMega(TGeometryId id)
Check if a TGeometryId is a MicroMega.
Definition: ND280GeomId.cxx:424
ND::GeomId::Def::TPC::kSeqIdMSB
const GeomIdBit kSeqIdMSB
Definition: ND280GeomIdDef.hxx:215
ND::GeomId::P0D::GetBarP0Dule
int GetBarP0Dule(TGeometryId id)
Definition: ND280GeomId.cxx:287
ND::GeomId::Def::SMRD::kSeqIdMSB
const GeomIdBit kSeqIdMSB
Definition: ND280GeomIdDef.hxx:569
ND::GeomId::GetSubSystem
int GetSubSystem(TGeometryId id)
Get the sub-system type.
Definition: ND280GeomId.cxx:67
ND::GeomId::ECal::GetBarECal
int GetBarECal(TGeometryId id)
Definition: ND280GeomId.cxx:1195
ND::GeomId::DSECal::Radiator
TGeometryId Radiator(int radiator)
The geometry identifier for the radiators.
Definition: ND280GeomId.cxx:1271
ND::GeomId::ECal::GetLayerECal
int GetLayerECal(TGeometryId id)
Definition: ND280GeomId.cxx:991
ND::GeomId::Def::ECal::Global::kECal
A special case for one of a kind global volumes.
Definition: ND280GeomIdDef.hxx:500
ND::GeomId::Def::TPC::Pad::kTPCMSB
const GeomIdBit kTPCMSB
The TPC withing the tracker. 0) TPC1 1) TPC2 2) TPC3.
Definition: ND280GeomIdDef.hxx:251
ND::GeomId::FGD::IsFGD
bool IsFGD(TGeometryId id)
Check if this is an FGD id.
Definition: ND280GeomId.cxx:561
ND::GeomId::Def::SMRD::Global::Module::kSlotMSB
const GeomIdBit kSlotMSB
The slot within a layer.
Definition: ND280GeomIdDef.hxx:637
ND::GeomId::P0D::P0Dule
TGeometryId P0Dule(int p0dule)
Definition: ND280GeomId.cxx:134
ND::GeomId::FGD::FGD2
TGeometryId FGD2()
The geometry identifier for FGD2. FGD2 has the water targets.
Definition: ND280GeomId.cxx:590
ND::GeomId::TECal::Layer
TGeometryId Layer(int clam, int module, int layer)
The layer within the ecal.
Definition: ND280GeomId.cxx:1310
ND::GeomId::P0D::IsP0D
bool IsP0D(TGeometryId id)
Check if the id is for the P0D.
Definition: ND280GeomId.cxx:77
ND::GeomId::Def::P0D::kSeqIdMSB
const GeomIdBit kSeqIdMSB
Definition: ND280GeomIdDef.hxx:112
ND::GeomId::TPC::TPC1
TGeometryId TPC1()
The geometry identifier for TPC1.
Definition: ND280GeomId.cxx:370
ND::GeomId::TECal::Bar
TGeometryId Bar(int clam, int module, int layer, int bar)
The geometry identifiers for the scintillators.
Definition: ND280GeomId.cxx:1320
ND::GeomId::Def::SMRD::Bar::kSlotMSB
const GeomIdBit kSlotMSB
The slot within a layer.
Definition: ND280GeomIdDef.hxx:679
Generated by   doxygen 1.8.13