1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
24 #include "AliITShit.h"
25 #include "AliITSmodule.h"
26 #include "AliITSgeomTGeo.h"
28 ClassImp(AliITSmodule)
30 //_______________________________________________________________________
32 // Impementation of class AliITSmodule
34 // created by: A. Bouchm, W. Peryt, S. Radomski, P. Skowronski
35 // R. Barbers, B. Batyunia, B. S. Nilsen
36 // ver 1.0 CERN 16.09.1999
37 //_______________________________________________________________________
38 //________________________________________________________________________
40 // Constructors and deconstructor
41 //________________________________________________________________________
43 AliITSmodule::AliITSmodule():
52 //_________________________________________________________________________
53 AliITSmodule::AliITSmodule(Int_t index):
61 fHitsM = new TObjArray();
62 fTrackIndex = new TArrayI(16);
63 fHitIndex = new TArrayI(16);
64 fITS = (AliITS*)(gAlice->GetDetector("ITS"));
66 //__________________________________________________________________________
67 AliITSmodule::~AliITSmodule() {
68 // The destructor for AliITSmodule. Before destoring AliITSmodule
69 // we must first destroy all of it's members.
78 fITS = 0; // We don't delete this pointer since it is just a copy.
80 //____________________________________________________________________________
81 AliITSmodule::AliITSmodule(const AliITSmodule &source):TObject(source),
83 fIndex(source.fIndex),
84 fHitsM(source.fHitsM),
85 fTrackIndex(source.fTrackIndex),
86 fHitIndex(source.fHitIndex){
87 ////////////////////////////////////////////////////////////////////////
89 ////////////////////////////////////////////////////////////////////////
92 //_____________________________________________________________________________
93 AliITSmodule& AliITSmodule::operator=(const AliITSmodule &source){
94 ////////////////////////////////////////////////////////////////////////
95 // Assignment operator
96 ////////////////////////////////////////////////////////////////////////
97 this->~AliITSmodule();
98 new(this) AliITSmodule(source);
101 //_________________________________________________________________________
104 //__________________________________________________________________________
105 Int_t AliITSmodule::AddHit(AliITShit* hit,Int_t t,Int_t h) {
108 //printf("AddHit: beginning hit %p t h %d %d\n",hit,t,h);
109 fHitsM->AddLast(new AliITShit(*hit));
110 Int_t fNhitsM = fHitsM->GetEntriesFast();
111 if(fNhitsM-1>=fTrackIndex->GetSize()){ // need to expand the TArrayI
112 fTrackIndex->Set(fNhitsM+64);
114 if(fNhitsM-1>=fHitIndex->GetSize()){ // need to expand the TArrayI
115 fHitIndex->Set(fNhitsM+64);
117 (*fTrackIndex)[fNhitsM-1] = t;
118 (*fHitIndex)[fNhitsM-1] = h;
121 //___________________________________________________________________________
122 Double_t AliITSmodule::PathLength(const AliITShit *itsHit1,
123 const AliITShit *itsHit2) {
129 itsHit1->GetPositionG(x1g,y1g,z1g);
130 itsHit2->GetPositionG(x2g,y2g,z2g);
132 s = TMath::Sqrt( ((Double_t)(x2g-x1g)*(Double_t)(x2g-x1g)) +
133 ((Double_t)(y2g-y1g)*(Double_t)(y2g-y1g)) +
134 ((Double_t)(z2g-z1g)*(Double_t)(z2g-z1g)) );
137 //___________________________________________________________________________
138 void AliITSmodule::PathLength(Float_t x,Float_t y,Float_t z,
139 Int_t status,Int_t &nseg,
140 Float_t &x1,Float_t &y1,Float_t &z1,
141 Float_t &dx1,Float_t &dy1,Float_t &dz1,
144 static Float_t x0,y0,z0;
146 if ((status&0x0002)!=0){ // entering
160 if ((status&0x0004)!=0) flag = 0; //exiting
161 if ((status&0x0001)!=0) flag = 2; // inside
162 else flag = 2; //inside ?
168 //___________________________________________________________________________
169 Bool_t AliITSmodule::LineSegmentL(Int_t hitindex,Double_t &a,Double_t &b,
170 Double_t &c,Double_t &d,
171 Double_t &e,Double_t &f,Double_t &de){
176 if(hitindex>= fHitsM->GetEntriesFast()) return kFALSE;
178 h1 = (AliITShit *) (fHitsM->At(hitindex));
179 if(h1->StatusEntering()){ // if track entering volume, get index for next
182 } // end if StatusEntering()
184 de = h1->GetIonization();
185 h1->GetPositionL0(a,c,e,t);
186 h1->GetPositionL(b,d,f);
192 //___________________________________________________________________________
193 Bool_t AliITSmodule::LineSegmentG(Int_t hitindex,Double_t &a,Double_t &b,
194 Double_t &c,Double_t &d,
195 Double_t &e,Double_t &f,Double_t &de){
200 if(hitindex>= fHitsM->GetEntriesFast()) return kFALSE;
202 h1 = (AliITShit *) (fHitsM->At(hitindex));
203 if(h1->StatusEntering()){ // if track entering volume, get index for next
206 } // end if StatusEntering()
208 de = h1->GetIonization();
209 h1->GetPositionG0(a,c,e,t);
210 h1->GetPositionG(b,d,f);
216 //___________________________________________________________________________
217 Bool_t AliITSmodule::LineSegmentL(Int_t hitindex,Double_t &a,Double_t &b,
218 Double_t &c,Double_t &d,
219 Double_t &e,Double_t &f,
220 Double_t &de,Int_t &track){
225 if(hitindex>= fHitsM->GetEntriesFast()) return kFALSE;
227 h1 = (AliITShit *) (fHitsM->At(hitindex));
228 if(h1->StatusEntering()){ // if track entering volume, get index for next
230 track = h1->GetTrack();
232 } // end if StatusEntering()
234 de = h1->GetIonization();
235 h1->GetPositionL0(a,c,e,t);
236 h1->GetPositionL(b,d,f);
240 track = h1->GetTrack();
243 //___________________________________________________________________________
244 Bool_t AliITSmodule::LineSegmentG(Int_t hitindex,Double_t &a,Double_t &b,
245 Double_t &c,Double_t &d,
246 Double_t &e,Double_t &f,
247 Double_t &de,Int_t &track){
252 if(hitindex>= fHitsM->GetEntriesFast()) return kFALSE;
254 h1 = (AliITShit *) (fHitsM->At(hitindex));
255 if(h1->StatusEntering()){ // if track entering volume, get index for next
257 track = h1->GetTrack();
259 } // end if StatusEntering()
261 de = h1->GetIonization();
262 h1->GetPositionG0(a,c,e,t);
263 h1->GetPositionG(b,d,f);
267 track = h1->GetTrack();
270 //______________________________________________________________________
271 Bool_t AliITSmodule::MedianHitG(AliITShit *h1,AliITShit *h2,
272 Float_t &x,Float_t &y,Float_t &z){
273 // Computes the mean hit location for a set of hits that make up a track
274 // passing through a volume. Returns kFALSE untill the the track leaves
277 Float_t x1l=0.,y1l=0.,z1l=0.;
278 Float_t x2l=0.,y2l=0.,z2l=0.;
279 Float_t xMl,yMl=0,zMl;
282 h1->GetPositionG(x1l,y1l,z1l);
283 h2->GetPositionG(x2l,y2l,z2l);
286 xMl = (-y1l / (y2l-y1l))*(x2l-x1l) + x1l;
287 zMl = (-y1l / (y2l-y1l))*(z2l-z1l) + z1l;
296 AliITSgeomTGeo::LocalToGlobal(h1->GetModule(),l,g);
302 //___________________________________________________________________________
303 void AliITSmodule::MedianHitG(Int_t index,
304 Float_t hitx1,Float_t hity1,Float_t hitz1,
305 Float_t hitx2,Float_t hity2,Float_t hitz2,
306 Float_t &xMg, Float_t &yMg, Float_t &zMg){
310 Float_t xMl,yMl=0,zMl;
316 AliITSgeomTGeo::GlobalToLocal(index,g,l);
324 AliITSgeomTGeo::GlobalToLocal(index,g,l);
330 xMl = (-y1l / (y2l-y1l))*(x2l-x1l) + x1l;
331 zMl = (-y1l / (y2l-y1l))*(z2l-z1l) + z1l;
340 AliITSgeomTGeo::LocalToGlobal(index,l,g);
345 //___________________________________________________________________________
346 Bool_t AliITSmodule::MedianHitL( AliITShit *itsHit1,
348 Float_t &xMl, Float_t &yMl, Float_t &zMl) const{
353 itsHit1->GetPositionL(x1l,y1l,z1l);
354 itsHit2->GetPositionL(x2l,y2l,z2l);
358 xMl = (-y1l / (y2l-y1l))*(x2l-x1l) + x1l;
359 zMl = (-y1l / (y2l-y1l))*(z2l-z1l) + z1l;
366 //___________________________________________________________________________
367 void AliITSmodule::MedianHit(Int_t index,
368 Float_t xg,Float_t yg,Float_t zg,
370 Float_t &xMg,Float_t &yMg,Float_t &zMg,
373 static Float_t x1,y1,z1;
375 if ((status&0x0002)!=0){ // entering
380 } else if ((status&0x0004)!=0){ // exiting
381 MedianHitG(index,x1,y1,z1,xg,yg,zg,xMg,yMg,zMg);
386 //___________________________________________________________________________
387 void AliITSmodule::GetID(Int_t &lay,Int_t &lad,Int_t &det){
389 AliITSgeomTGeo::GetModuleId(fIndex,lay,lad,det);