]>
Commit | Line | Data |
---|---|---|
1 | /************************************************************************** | |
2 | * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
3 | * * | |
4 | * Author: The ALICE Off-line Project. * | |
5 | * Contributors are mentioned in the code where appropriate. * | |
6 | * * | |
7 | * Permission to use, copy, modify and distribute this software and its * | |
8 | * documentation strictly for non-commercial purposes is hereby granted * | |
9 | * without fee, provided that the above copyright notice appears in all * | |
10 | * copies and that both the copyright notice and this permission notice * | |
11 | * appear in the supporting documentation. The authors make no claims * | |
12 | * about the suitability of this software for any purpose. It is * | |
13 | * provided "as is" without express or implied warranty. * | |
14 | **************************************************************************/ | |
15 | ||
16 | /* | |
17 | $Log$ | |
18 | Revision 1.4 2000/10/16 01:16:53 cblume | |
19 | Changed timebin 0 to be the one closest to the readout | |
20 | ||
21 | Revision 1.3 2000/10/15 23:40:01 cblume | |
22 | Remove AliTRDconst | |
23 | ||
24 | Revision 1.2 2000/10/06 16:49:46 cblume | |
25 | Made Getters const | |
26 | ||
27 | Revision 1.1.2.2 2000/10/04 16:34:58 cblume | |
28 | Replace include files by forward declarations | |
29 | ||
30 | Revision 1.1.2.1 2000/09/22 14:47:52 cblume | |
31 | Add the tracking code | |
32 | ||
33 | */ | |
34 | ||
35 | #include <TObject.h> | |
36 | ||
37 | #include "AliRun.h" | |
38 | ||
39 | #include "AliTRD.h" | |
40 | #include "AliTRDgeometry.h" | |
41 | #include "AliTRDcluster.h" | |
42 | #include "AliTRDtimeBin.h" | |
43 | #include "AliTRDtrackingSector.h" | |
44 | ||
45 | ClassImp(AliTRDtrackingSector) | |
46 | ||
47 | //_______________________________________________________ | |
48 | ||
49 | AliTRDtrackingSector::~AliTRDtrackingSector() | |
50 | { | |
51 | // | |
52 | // Destructor | |
53 | // | |
54 | ||
55 | delete[] fTimeBin; | |
56 | ||
57 | } | |
58 | ||
59 | //_______________________________________________________ | |
60 | AliTRDtimeBin &AliTRDtrackingSector::operator[](Int_t i) | |
61 | { | |
62 | // | |
63 | // Index operator | |
64 | // | |
65 | ||
66 | return *(fTimeBin+i); | |
67 | ||
68 | } | |
69 | ||
70 | //_______________________________________________________ | |
71 | ||
72 | void AliTRDtrackingSector::SetUp() | |
73 | { | |
74 | AliTRD *TRD = (AliTRD*) gAlice->GetDetector("TRD"); | |
75 | fGeom = TRD->GetGeometry(); | |
76 | ||
77 | fTimeBinSize = fGeom->GetTimeBinSize(); | |
78 | ||
79 | fN = AliTRDgeometry::Nplan() * (Int_t(AliTRDgeometry::DrThick() | |
80 | /fTimeBinSize) + 1); | |
81 | ||
82 | fTimeBin = new AliTRDtimeBin[fN]; | |
83 | ||
84 | } | |
85 | ||
86 | //______________________________________________________ | |
87 | ||
88 | Double_t AliTRDtrackingSector::GetX(Int_t l) const | |
89 | { | |
90 | if( (l<0) || (l>fN-1)) { | |
91 | fprintf(stderr,"AliTRDtrackingSector::GetX: TimeBin index is out of range !\n"); | |
92 | return -99999.; | |
93 | } | |
94 | else { | |
95 | ||
96 | Int_t tb_per_plane = fN/AliTRDgeometry::Nplan(); | |
97 | Int_t plane = l/tb_per_plane; | |
98 | Int_t time_slice = l%(Int_t(AliTRDgeometry::DrThick() | |
99 | /fTimeBinSize) + 1); | |
100 | ||
101 | Float_t t0 = fGeom->GetTime0(plane); | |
102 | Double_t x = t0 - time_slice * fTimeBinSize; | |
103 | ||
104 | return x; | |
105 | } | |
106 | } | |
107 | ||
108 | //______________________________________________________ | |
109 | ||
110 | Double_t AliTRDtrackingSector::GetMaxY(Int_t l) const | |
111 | { | |
112 | ||
113 | if((l<(fN-1)) && (l>-1)) { | |
114 | Int_t tb_per_plane = fN/AliTRDgeometry::Nplan(); | |
115 | Int_t plane = l/tb_per_plane; | |
116 | return fGeom->GetChamberWidth(plane); | |
117 | } | |
118 | else { | |
119 | fprintf(stderr, | |
120 | "AliTRDtrackingSector::GetMaxY: TimeBin index is out of range !\n"); | |
121 | if(l<0) return fGeom->GetChamberWidth(0); | |
122 | else return fGeom->GetChamberWidth(AliTRDgeometry::Nplan()-1); | |
123 | } | |
124 | } | |
125 | ||
126 | //______________________________________________________ | |
127 | ||
128 | Int_t AliTRDtrackingSector::GetTimeBinNumber(Double_t x) const | |
129 | { | |
130 | Float_t r_out = fGeom->GetTime0(AliTRDgeometry::Nplan()-1); | |
131 | Float_t r_in = fGeom->GetTime0(0) - AliTRDgeometry::DrThick(); | |
132 | ||
133 | if(x >= r_out) return fN-1; | |
134 | if(x <= r_in) return 0; | |
135 | ||
136 | Float_t gap = fGeom->GetTime0(1) - fGeom->GetTime0(0); | |
137 | ||
138 | Int_t plane = Int_t((x - r_in + fTimeBinSize/2)/gap); | |
139 | ||
140 | Int_t local_tb = Int_t((fGeom->GetTime0(plane)-x)/fTimeBinSize + 0.5); | |
141 | ||
142 | ||
143 | Int_t tb_per_plane = fN/AliTRDgeometry::Nplan(); | |
144 | ||
145 | Int_t time_bin = plane * (Int_t(AliTRDgeometry::DrThick()/fTimeBinSize) + 1) | |
146 | + (tb_per_plane - 1 - local_tb); | |
147 | ||
148 | return time_bin; | |
149 | } | |
150 | ||
151 | //______________________________________________________ | |
152 | ||
153 | Int_t AliTRDtrackingSector::GetTimeBin(Int_t det, Int_t local_tb) const | |
154 | { | |
155 | Int_t plane = fGeom->GetPlane(det); | |
156 | ||
157 | Int_t tb_per_plane = fN/AliTRDgeometry::Nplan(); | |
158 | ||
159 | Int_t time_bin = plane * (Int_t(AliTRDgeometry::DrThick()/fTimeBinSize) + 1) | |
160 | + (tb_per_plane - 1 - local_tb); | |
161 | ||
162 | return time_bin; | |
163 | } | |
164 |