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 **************************************************************************/
16 /* $Id: AliITSUChip.cxx 53509 2011-12-10 18:55:52Z masera $ */
24 #include "AliITSUChip.h"
25 #include "AliITSUGeomTGeo.h"
26 using namespace TMath;
30 //_______________________________________________________________________
32 // Impementation of class AliITSUChip
34 //_______________________________________________________________________
35 //________________________________________________________________________
37 // Constructors and deconstructor
38 //________________________________________________________________________
40 AliITSUChip::AliITSUChip()
47 //_________________________________________________________________________
48 AliITSUChip::AliITSUChip(Int_t index, AliITSUGeomTGeo* tg)
49 :fHitsM(new TObjArray())
56 //__________________________________________________________________________
57 AliITSUChip::~AliITSUChip()
59 // The destructor for AliITSUChip. Before destoring AliITSUChip
60 // we must first destroy all of it's members.
61 if (fHitsM) fHitsM->Clear();
65 //____________________________________________________________________________
66 AliITSUChip::AliITSUChip(const AliITSUChip &source)
68 ,fHitsM(source.fHitsM)
69 ,fGeomTG(source.fGeomTG)
74 //_____________________________________________________________________________
75 AliITSUChip& AliITSUChip::operator=(const AliITSUChip &source)
77 // Assignment operator
80 new(this) AliITSUChip(source);
85 //___________________________________________________________________________
86 Double_t AliITSUChip::PathLength(const AliITSUHit *itsHit1,const AliITSUHit *itsHit2)
93 itsHit1->GetPositionG(x1g,y1g,z1g);
94 itsHit2->GetPositionG(x2g,y2g,z2g);
96 s = Sqrt( ((Double_t)(x2g-x1g)*(Double_t)(x2g-x1g)) +
97 ((Double_t)(y2g-y1g)*(Double_t)(y2g-y1g)) +
98 ((Double_t)(z2g-z1g)*(Double_t)(z2g-z1g)) );
102 //___________________________________________________________________________
103 void AliITSUChip::PathLength(Float_t x,Float_t y,Float_t z,Int_t status,Int_t &nseg,Float_t &x1,Float_t &y1,Float_t &z1,
104 Float_t &dx1,Float_t &dy1,Float_t &dz1,Int_t &flag) const
107 static Float_t x0,y0,z0;
109 if ((status&0x0002)!=0){ // entering
123 if ((status&0x0004)!=0) flag = 0; //exiting
124 if ((status&0x0001)!=0) flag = 2; // inside
125 else flag = 2; //inside ?
132 //___________________________________________________________________________
133 Bool_t AliITSUChip::LineSegmentL(Int_t hitindex,Double_t &a,Double_t &b,Double_t &c,Double_t &d,
134 Double_t &e,Double_t &f,Double_t &de)
140 if(hitindex>= fHitsM->GetEntriesFast()) return kFALSE;
142 h1 = (AliITSUHit *) (fHitsM->At(hitindex));
143 if(h1->StatusEntering()) return kFALSE; // if track entering volume, get index for next step
144 de = h1->GetIonization();
145 h1->GetPositionL0(a,c,e,t);
146 h1->GetPositionL(b,d,f);
153 //___________________________________________________________________________
154 Bool_t AliITSUChip::LineSegmentG(Int_t hitindex,Double_t &a,Double_t &b,Double_t &c,Double_t &d,
155 Double_t &e,Double_t &f,Double_t &de)
161 if(hitindex>= fHitsM->GetEntriesFast()) return kFALSE;
162 h1 = (AliITSUHit *) (fHitsM->At(hitindex));
163 if(h1->StatusEntering()) return kFALSE; // if track entering volume, get index for next step
164 de = h1->GetIonization();
165 h1->GetPositionG0(a,c,e,t);
166 h1->GetPositionG(b,d,f);
173 //___________________________________________________________________________
174 Bool_t AliITSUChip::LineSegmentL(Int_t hitindex,Double_t &a,Double_t &b,Double_t &c,Double_t &d,
175 Double_t &e,Double_t &f,Double_t &de,Double_t &tof, Int_t &track)
179 if(hitindex>= fHitsM->GetEntriesFast()) return kFALSE;
181 h1 = (AliITSUHit *) (fHitsM->At(hitindex));
182 if(h1->StatusEntering()){ // if track entering volume, get index for next
184 track = h1->GetTrack();
186 } // end if StatusEntering()
188 de = h1->GetIonization();
189 h1->GetPositionL0(a,c,e,tof);
190 h1->GetPositionL(b,d,f);
194 track = h1->GetTrack();
198 //___________________________________________________________________________
199 Bool_t AliITSUChip::LineSegmentG(Int_t hitindex,Double_t &a,Double_t &b,Double_t &c,Double_t &d,
200 Double_t &e,Double_t &f,Double_t &de,Double_t &tof, Int_t &track)
204 if(hitindex>= fHitsM->GetEntriesFast()) return kFALSE;
206 h1 = (AliITSUHit *) (fHitsM->At(hitindex));
207 if(h1->StatusEntering()){ // if track entering volume, get index for next
209 track = h1->GetTrack();
211 } // end if StatusEntering()
212 de = h1->GetIonization();
213 h1->GetPositionG0(a,c,e,tof);
214 h1->GetPositionG(b,d,f);
218 track = h1->GetTrack();
222 //______________________________________________________________________
223 Bool_t AliITSUChip::MedianHitG(AliITSUHit *h1,AliITSUHit *h2,Float_t &x,Float_t &y,Float_t &z)
225 // Computes the mean hit location for a set of hits that make up a track
226 // passing through a volume. Returns kFALSE untill the the track leaves
229 Float_t x1l=0.,y1l=0.,z1l=0.;
230 Float_t x2l=0.,y2l=0.,z2l=0.;
231 Float_t xMl,yMl=0,zMl;
234 h1->GetPositionG(x1l,y1l,z1l);
235 h2->GetPositionG(x2l,y2l,z2l);
238 xMl = (-y1l / (y2l-y1l))*(x2l-x1l) + x1l;
239 zMl = (-y1l / (y2l-y1l))*(z2l-z1l) + z1l;
248 fGeomTG->LocalToGlobal(h1->GetChip(),l,g);
255 //___________________________________________________________________________
256 void AliITSUChip::MedianHitG(Int_t index,Float_t hitx1,Float_t hity1,Float_t hitz1,Float_t hitx2,Float_t hity2,Float_t hitz2,
257 Float_t &xMg, Float_t &yMg, Float_t &zMg)
262 Float_t xMl,yMl=0,zMl;
268 fGeomTG->GlobalToLocal(index,g,l);
276 fGeomTG->GlobalToLocal(index,g,l);
282 xMl = (-y1l / (y2l-y1l))*(x2l-x1l) + x1l;
283 zMl = (-y1l / (y2l-y1l))*(z2l-z1l) + z1l;
292 fGeomTG->LocalToGlobal(index,l,g);
298 //___________________________________________________________________________
299 Bool_t AliITSUChip::MedianHitL( AliITSUHit *itsHit1,AliITSUHit *itsHit2, Float_t &xMl, Float_t &yMl, Float_t &zMl) const
305 itsHit1->GetPositionL(x1l,y1l,z1l);
306 itsHit2->GetPositionL(x2l,y2l,z2l);
310 xMl = (-y1l / (y2l-y1l))*(x2l-x1l) + x1l;
311 zMl = (-y1l / (y2l-y1l))*(z2l-z1l) + z1l;
319 //___________________________________________________________________________
320 void AliITSUChip::MedianHit(Int_t index,Float_t xg,Float_t yg,Float_t zg,Int_t status,
321 Float_t &xMg,Float_t &yMg,Float_t &zMg,Int_t &flag)
324 static Float_t x1,y1,z1;
326 if ((status&0x0002)!=0){ // entering
331 } else if ((status&0x0004)!=0){ // exiting
332 MedianHitG(index,x1,y1,z1,xg,yg,zg,xMg,yMg,zMg);