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 **************************************************************************/
17 ///////////////////////////////////////////////////////////////////////////////
19 // This class is a container class for the laser beam positions of //
20 // the TPC laser system (ALICE-INT-2002-022 v.1 ). //
21 // The system consits of 6 Laser Rods, each contains four mirror bundels //
22 // on each side of the TPC. Each bundle has seven mirrors. //
24 // The TPC side 0 corresponds to the "Shaft Side (A)", side 1 corresponds //
25 // to the "Muon Side (C)".
26 // The laser rods are counted counter clockwise starting with the one //
27 // with the smalles phi angle. //
28 // The mirror bundles in one rod for each side are counted starting from //
29 // form the greatest z value to the smalles. //
30 // The beams in each bundel are counted counter clockwise. Having all beams //
31 // pointing downward, the first beam will be the leftmost. //
33 ///////////////////////////////////////////////////////////////////////////////
36 Create File with design positions
38 .L AliTPCLaserTrack.cxx+
39 AliTPCLaserTracks ltracks
40 ltracks.WriteTreeDesignData()
41 TFile f("LaserTracksDesing.root");
42 TTree * tree =(TTree*)f.Get("LaserTracks");
46 #include <Riostream.h>
48 #include <TPolyLine3D.h>
52 #include <TEventList.h>
60 #include "AliTPCLaserTracks.h"
62 ////////////////////////////////////////////////////////////////////////
63 // Class AliTPCLaserTracks
64 ////////////////////////////////////////////////////////////////////////
66 ClassImp(AliTPCLaserTracks)
68 AliTPCLaserTracks::AliTPCLaserTracks():TObject(),
88 // AliTPCLaserTracks default constructor
92 //_______________________________________________________________________
93 AliTPCLaserTracks::AliTPCLaserTracks(Int_t npoints):TObject(),
113 // AliTPCLaserTracks constructor to initialise array of points
118 //______________________________________________________________________
119 AliTPCLaserTracks::AliTPCLaserTracks(const AliTPCLaserTracks ¶m):TObject(),
142 //______________________________________________________________________
143 AliTPCLaserTracks & AliTPCLaserTracks::operator=(const AliTPCLaserTracks & param)
146 // assignment operator - dummy
152 //_______________________________________________________________________
153 Double_t AliTPCLaserTracks::FindBeamLength(TVector3 vF, TVector3 vP)
156 // Find distance between mirror and the intersection between
157 // inner and outer field cage respectively
159 // use the pq formula to find the intersection point between
160 // the beam and the inner and outer fieldcage cylinders
161 // The formulae solved are
164 // xPoint = vF.X() + n*vP.X() = r*cos(phi)
165 // yPoint = vF.Y() + n*vP.Y() = r*sin(phi)
166 // zPoint = vF.Y() + n*vP.Y()
168 // get n from the first two equations, calculate x,y,zPoint
170 // vF is the mirror position and vP the beam direction
171 Double_t r[2]; //smalles r[0] and largest r[1] pad radius
172 r[0] = 80.0; //smalles pad radius
173 r[1] = 260.0; //largest pad radius
175 Double_t fxpxfypy = vF.X()*vP.X()+vF.Y()*vP.Y();
176 Double_t px2py2 = vP.X()*vP.X()+vP.Y()*vP.Y();
178 for (Int_t i=0; i<2; i++){
179 Double_t rad = (fxpxfypy/px2py2)*(fxpxfypy/px2py2)-
180 (vF.X()*vF.X()+vF.Y()*vF.Y()-r[i]*r[i])/px2py2;
183 Double_t n1 = -(fxpxfypy)/px2py2+TMath::Sqrt(rad);
184 TVector3 vI1(vF.X()+n1*vP.X(),
188 Double_t n2 = -(fxpxfypy)/px2py2-TMath::Sqrt(rad);
189 TVector3 vI2(vF.X()+n2*vP.X(),
194 //if we cross two boarders on our way return the closer boarder
196 if ( (vF-vI1).Mag() <= (vF-vI2).Mag() )
197 return (vF-vI1).Mag();
199 return (vF-vI2).Mag();
202 return (vF-vI1).Mag();
205 return (vF-vI2).Mag();
210 return 3.4e38; //not found!!
213 //_______________________________________________________________________
214 void AliTPCLaserTracks::WriteTreeDesignData()
217 // Write a tree with the design data of the laser track positions
220 TFile *f = TFile::Open("LaserTracksDesing.root","recreate");
222 AliTPCLaserTracks *ltp = new AliTPCLaserTracks(2);
224 TTree *t = new TTree("LaserTracks","LaserTracks");
225 t->Branch("LaserTracks","AliTPCLaserTracks",<p);
228 Double_t phiBeam[7] = {-31.8,-16.,-9.2,2.5,9.2,16.,31.8};
229 Double_t phiRod[6] = {40.,100.,160.,220.,280.,340.}; //-20. for the muon side
230 // Double_t zBundleCorse[2][4] = {{10,79,163,241},{13,85,169,247}};
231 Double_t zBundleCorse[2][4] = {{241.,163.,79.,10.},{247.,169.,85.,13.}};
232 Double_t zBeamFine[7] = {1.45,1.3,1.15,1.3,1.15,1.3,1.45};
235 //loop over TPC side -- 0:A (Shaft) side -- 1:C (Muon) side
236 for (Int_t side=0; side<2; side++){
237 //loop over laser rods -- counterclockwise
238 for (Int_t rod=0; rod<6; rod++){
242 Double_t phiRod2 = phiRod[rod]-side*20;
243 vRod.SetMagPhi(254.25, //center of the rod at 254.25cm
245 (phiRod2) //-20 deg on C-Side
248 //loop over bundle -- counted from large z to small z
249 for (Int_t bundle=0; bundle<4; bundle++){
250 //center of the bundle; not yet rotated
251 Int_t bundleS = bundle;
252 // if ( side == 1 ) bundleS = 4-bundle;
253 if ( side == 1 ) bundleS = 3-bundle;
256 vBundle.SetMagPhi(.65,
258 //? TMath::Power(-1,side)*
259 (phiRod2+180.-90.*bundleS)
264 for (Int_t beam=0; beam<7; beam++){
269 vBeam.SetMagPhi(1.,TMath::DegToRad()*(phiRod2+i*60));
273 TVector2 xyMirror = vRod+vBundle+vBeam;
274 Double_t zMirror = (zBundleCorse[rod%2][bundleS]+zBeamFine[beam])*TMath::Power(-1,side);
275 TVector3 v3Mirror(xyMirror.X(),xyMirror.Y(),zMirror);
277 Double_t phi = TMath::DegToRad()*(phiRod2+180+phiBeam[beam]*TMath::Power(-1,side));
278 Double_t theta = TMath::DegToRad()*90;
281 v3Beam.SetMagThetaPhi(1,theta,phi); //Direction Vector
282 v3Beam.SetMagThetaPhi(FindBeamLength(v3Mirror,v3Beam), theta, phi);
289 ltp->SetBundle(bundleS);
292 ltp->SetX(xyMirror.X());
293 ltp->SetY(xyMirror.Y());
297 ltp->SetTheta(theta);
299 // ltp->InitPoints(2);
300 ltp->SetPoint(0,xyMirror.X(),xyMirror.Y(),zMirror);
301 ltp->SetPoint(1,xyMirror.X()+v3Beam.X(),xyMirror.Y()+v3Beam.Y(),zMirror+v3Beam.Z());
306 << " Bundle: " << bundleS
314 // cout << "Filled!!!" << endl;
327 //_______________________________________________________________________
328 TPolyLine3D* AliTPCLaserTracks::GetLine()
331 // Make a polilyne object oout of the points
333 if ( fNpoints ) return new TPolyLine3D(fNpoints,fXarr,fYarr,fZarr);
338 //_______________________________________________________________________
339 TObjArray* AliTPCLaserTracks::GetLines(Char_t* file, Char_t *cuts)
342 // Read the input files with the laser track
343 // Make a array of polylines for visualization
345 TObjArray *array = new TObjArray;
347 TFile *f = TFile::Open(file,"read");
348 TTree *tree = (TTree*)f->Get("LaserTracks");
350 AliTPCLaserTracks *ltp = new AliTPCLaserTracks();
351 // TEventList *evList = new TEventList("evList");
352 TEventList evList("evList");
354 tree->SetBranchAddress("LaserTracks",<p);
356 tree->Draw(">>evList",cuts,"goff");
358 // cout << "N: " << evList.GetN() << endl;
360 for (Int_t ev = 0; ev < evList.GetN(); ev++){
361 // cout << ev << endl;
362 tree->GetEntry( evList.GetEntry(ev) );
363 array->Add(ltp->GetLine());
365 // cout << "delte f"<< endl;
368 // cout << "evlist" << endl;
373 //_______________________________________________________________________
374 void AliTPCLaserTracks::InitPoints()
381 fXarr = new Double_t[fMaxSize];
382 fYarr = new Double_t[fMaxSize];
383 fZarr = new Double_t[fMaxSize];
386 //_______________________________________________________________________
387 Int_t AliTPCLaserTracks::SetPoint(Int_t point, Double_t x, Double_t y, Double_t z)
390 // Set point to point arrays
393 if ( !fXarr || !fYarr || !fZarr ) return 1;
394 if ( point > fNpoints-1 ) return 1;
402 //_______________________________________________________________________
403 Int_t AliTPCLaserTracks::FindMirror(Char_t *file, Double_t x, Double_t y, Double_t z, Double_t phi)
406 // return the id of the mirror in 'file' that maches best the given parameters
409 TFile *f = TFile::Open(file,"read");
410 TTree *tree = (TTree*)f->Get("LaserTracks");
412 AliTPCLaserTracks *ltp = new AliTPCLaserTracks();
413 TEventList evList("evList");
415 tree->SetBranchAddress("LaserTracks",<p);
418 s = "abs(fX-"; s+= x;
419 s+= ")<3&&abs(fY-"; s+=y;
420 s+= ")<3&&abs(fZ-"; s+= z;
422 // s+= "&&((abs(fPhi-";s+= phi;
423 // s+= ")<.06)||(abs(fPhi-";s+= (phi-TMath::DegToRad()*180);
425 cout << s.Data() << endl;
427 tree->Draw(">>evList",s.Data());
429 Double_t dphiMin=TMath::Pi();
432 cout << "nev: " << evList.GetN() << endl;
433 for (Int_t i=0; i<evList.GetN(); i++){
434 tree->GetEntry( evList.GetEntry(i) );
435 Double_t dphi = TMath::Abs(ltp->GetPhi() - phi);
436 if ( dphi > TMath::DegToRad()*180 ) dphi-=TMath::DegToRad()*180;
437 if ( dphi<dphiMin ) {
446 //_______________________________________________________________________
447 AliTPCLaserTracks::~AliTPCLaserTracks()
450 // standard destructor
453 // if ( fP ) delete fP;
454 if ( fXarr ) delete [] fXarr;
455 if ( fYarr ) delete [] fYarr;
456 if ( fZarr ) delete [] fZarr;