Changed timebin 0 to be the one closest to the readout
[u/mrichter/AliRoot.git] / TRD / AliTRDtrackingSector.cxx
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.3  2000/10/15 23:40:01  cblume
19 Remove AliTRDconst
20
21 Revision 1.2  2000/10/06 16:49:46  cblume
22 Made Getters const
23
24 Revision 1.1.2.2  2000/10/04 16:34:58  cblume
25 Replace include files by forward declarations
26
27 Revision 1.1.2.1  2000/09/22 14:47:52  cblume
28 Add the tracking code
29
30 */                        
31                                 
32 #include <TObject.h>
33
34 #include "AliRun.h"
35
36 #include "AliTRD.h" 
37 #include "AliTRDgeometry.h" 
38 #include "AliTRDcluster.h" 
39 #include "AliTRDtimeBin.h" 
40 #include "AliTRDtrackingSector.h" 
41
42 ClassImp(AliTRDtrackingSector) 
43
44 //_______________________________________________________
45
46 AliTRDtrackingSector::~AliTRDtrackingSector()
47 {
48   //
49   // Destructor
50   //
51
52   delete[] fTimeBin;
53
54 }
55
56 //_______________________________________________________
57 AliTRDtimeBin &AliTRDtrackingSector::operator[](Int_t i)
58 {
59   //
60   // Index operator 
61   //
62
63   return *(fTimeBin+i);
64
65 }
66
67 //_______________________________________________________
68
69 void AliTRDtrackingSector::SetUp()
70
71
72   AliTRD *TRD = (AliTRD*) gAlice->GetDetector("TRD");
73   fGeom = TRD->GetGeometry();
74
75   fTimeBinSize = fGeom->GetTimeBinSize();
76
77   fN = AliTRDgeometry::Nplan() * (Int_t(AliTRDgeometry::DrThick()
78                                        /fTimeBinSize) + 1);
79
80   fTimeBin = new AliTRDtimeBin[fN]; 
81
82 }
83
84 //______________________________________________________
85
86 Double_t AliTRDtrackingSector::GetX(Int_t l) const
87 {
88   if( (l<0) || (l>fN-1)) {
89     fprintf(stderr,"AliTRDtrackingSector::GetX: TimeBin index is out of range !\n");
90     return -99999.;
91   }
92   else { 
93     
94     Int_t tb_per_plane = fN/AliTRDgeometry::Nplan();
95     Int_t plane = l/tb_per_plane;
96     Int_t time_slice = l%(Int_t(AliTRDgeometry::DrThick()
97                                /fTimeBinSize) + 1);
98
99     Float_t t0 = fGeom->GetTime0(plane);
100     Double_t x = t0 + time_slice * fTimeBinSize;
101
102     //    cerr<<"plane, tb, x = "<<plane<<","<<time_slice<<","<<x<<endl;
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   //              + AliTRDgeometry::DrThick(); 
132   // Changed to new time0 (CBL)
133   Float_t r_out = fGeom->GetTime0(AliTRDgeometry::Nplan()-1); 
134   Float_t r_in = fGeom->GetTime0(0);
135   //  cerr<<"GetTimeBinNumber: r_in,r_out = "<<r_in<<","<<r_out<<endl;
136
137   if(x >= r_out) return fN-1;
138   if(x <= r_in) return 0;
139
140   Float_t gap = fGeom->GetTime0(1) - fGeom->GetTime0(0);
141   //  cerr<<"GetTimeBinNumber: gap = "<<gap<<endl;
142
143   Int_t plane = Int_t((x - r_in + fTimeBinSize/2)/gap);
144   //  cerr<<"GetTimeBinNumber: plane="<<plane<<endl;
145
146   Int_t local_tb = Int_t((x-fGeom->GetTime0(plane))/fTimeBinSize + 0.5);
147   //  cerr<<"GetTimeBinNumber: local_tb="<<local_tb<<endl;
148
149   Int_t time_bin = plane * (Int_t(AliTRDgeometry::DrThick()
150                                  /fTimeBinSize) + 1) + local_tb;
151   //  cerr<<"GetTimeBinNumber: time_bin = "<<time_bin<<endl;
152   
153
154   return time_bin;
155 }
156
157 //______________________________________________________
158
159 Int_t AliTRDtrackingSector::GetTimeBin(Int_t det, Int_t local_tb) const 
160 {
161   Int_t plane = fGeom->GetPlane(det);
162
163   Int_t time_bin = plane * (Int_t(AliTRDgeometry::DrThick()
164                                  /fTimeBinSize) + 1) + local_tb;  
165   return time_bin;
166 }
167