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
195 if ( (n1>0) && (n2>0) ) {
196 if ( (vF-vI1).Mag() <= (vF-vI2).Mag() )
197 return (vF-vI1).Mag();
199 return (vF-vI2).Mag();
203 return (vF-vI1).Mag();
206 return (vF-vI2).Mag();
211 return 3.4e38; //not found!!
214 //_______________________________________________________________________
215 void AliTPCLaserTracks::WriteTreeDesignData()
218 // Write a tree with the design data of the laser track positions
221 TFile *f = TFile::Open("LaserTracksDesing.root","recreate");
223 AliTPCLaserTracks *ltp = new AliTPCLaserTracks(2);
225 TTree *t = new TTree("LaserTracks","LaserTracks");
226 t->Branch("LaserTracks","AliTPCLaserTracks",<p);
229 Double_t phiBeam[7] = {-31.8,-16.,-9.2,2.5,9.2,16.,31.8};
230 Double_t phiRod[6] = {40.,100.,160.,220.,280.,340.}; //-20. for the muon side
231 // Double_t zBundleCorse[2][4] = {{10,79,163,241},{13,85,169,247}};
232 Double_t zBundleCorse[2][4] = {{241.,163.,79.,10.},{247.,169.,85.,13.}};
233 Double_t zBeamFine[7] = {1.45,1.3,1.15,1.3,1.15,1.3,1.45};
236 //loop over TPC side -- 0:A (Shaft) side -- 1:C (Muon) side
237 for (Int_t side=0; side<2; side++){
238 //loop over laser rods -- counterclockwise
239 for (Int_t rod=0; rod<6; rod++){
243 Double_t phiRod2 = phiRod[rod]-side*20;
244 vRod.SetMagPhi(254.25, //center of the rod at 254.25cm
246 (phiRod2) //-20 deg on C-Side
249 //loop over bundle -- counted from large z to small z
250 for (Int_t bundle=0; bundle<4; bundle++){
251 //center of the bundle; not yet rotated
252 Int_t bundleS = bundle;
253 // if ( side == 1 ) bundleS = 4-bundle;
254 if ( side == 1 ) bundleS = 3-bundle;
257 vBundle.SetMagPhi(.65,
259 //? TMath::Power(-1,side)*
260 (phiRod2+180.-90.*bundleS)
265 for (Int_t beam=0; beam<7; beam++){
270 vBeam.SetMagPhi(1.,TMath::DegToRad()*(phiRod2+i*60));
274 TVector2 xyMirror = vRod+vBundle+vBeam;
275 Double_t zMirror = (zBundleCorse[rod%2][bundleS]+zBeamFine[beam])*TMath::Power(-1,side);
276 TVector3 v3Mirror(xyMirror.X(),xyMirror.Y(),zMirror);
278 Double_t phi = TMath::DegToRad()*(phiRod2+180+phiBeam[beam]*TMath::Power(-1,side));
279 Double_t theta = TMath::DegToRad()*90;
282 v3Beam.SetMagThetaPhi(1,theta,phi); //Direction Vector
283 v3Beam.SetMagThetaPhi(FindBeamLength(v3Mirror,v3Beam), theta, phi);
290 ltp->SetBundle(bundleS);
293 ltp->SetX(xyMirror.X());
294 ltp->SetY(xyMirror.Y());
298 ltp->SetTheta(theta);
300 // ltp->InitPoints(2);
301 ltp->SetPoint(0,xyMirror.X(),xyMirror.Y(),zMirror);
302 ltp->SetPoint(1,xyMirror.X()+v3Beam.X(),xyMirror.Y()+v3Beam.Y(),zMirror+v3Beam.Z());
307 << " Bundle: " << bundleS
315 // cout << "Filled!!!" << endl;
328 //_______________________________________________________________________
329 TPolyLine3D* AliTPCLaserTracks::GetLine()
332 // Make a polilyne object oout of the points
334 if ( fNpoints ) return new TPolyLine3D(fNpoints,fXarr,fYarr,fZarr);
339 //_______________________________________________________________________
340 TObjArray* AliTPCLaserTracks::GetLines(const Char_t* file, const Char_t *cuts)
343 // Read the input files with the laser track
344 // Make a array of polylines for visualization
346 TObjArray *array = new TObjArray;
348 TFile *f = TFile::Open(file,"read");
349 TTree *tree = (TTree*)f->Get("LaserTracks");
351 AliTPCLaserTracks *ltp = new AliTPCLaserTracks();
352 // TEventList *evList = new TEventList("evList");
353 TEventList evList("evList");
355 tree->SetBranchAddress("LaserTracks",<p);
357 tree->Draw(">>evList",cuts,"goff");
359 // cout << "N: " << evList.GetN() << endl;
361 for (Int_t ev = 0; ev < evList.GetN(); ev++){
362 // cout << ev << endl;
363 tree->GetEntry( evList.GetEntry(ev) );
364 array->Add(ltp->GetLine());
366 // cout << "delte f"<< endl;
369 // cout << "evlist" << endl;
374 //_______________________________________________________________________
375 void AliTPCLaserTracks::InitPoints()
382 fXarr = new Double_t[fMaxSize];
383 fYarr = new Double_t[fMaxSize];
384 fZarr = new Double_t[fMaxSize];
387 //_______________________________________________________________________
388 Int_t AliTPCLaserTracks::SetPoint(Int_t point, Double_t x, Double_t y, Double_t z)
391 // Set point to point arrays
394 if ( !fXarr || !fYarr || !fZarr ) return 1;
395 if ( point > fNpoints-1 ) return 1;
403 //_______________________________________________________________________
404 Int_t AliTPCLaserTracks::FindMirror(Char_t *file, Double_t x, Double_t y, Double_t z, Double_t phi)
407 // return the id of the mirror in 'file' that maches best the given parameters
410 TFile *f = TFile::Open(file,"read");
411 TTree *tree = (TTree*)f->Get("LaserTracks");
413 AliTPCLaserTracks *ltp = new AliTPCLaserTracks();
414 TEventList evList("evList");
416 tree->SetBranchAddress("LaserTracks",<p);
419 s = "abs(fX-"; s+= x;
420 s+= ")<3&&abs(fY-"; s+=y;
421 s+= ")<3&&abs(fZ-"; s+= z;
423 // s+= "&&((abs(fPhi-";s+= phi;
424 // s+= ")<.06)||(abs(fPhi-";s+= (phi-TMath::DegToRad()*180);
426 cout << s.Data() << endl;
428 tree->Draw(">>evList",s.Data());
430 Double_t dphiMin=TMath::Pi();
433 cout << "nev: " << evList.GetN() << endl;
434 for (Int_t i=0; i<evList.GetN(); i++){
435 tree->GetEntry( evList.GetEntry(i) );
436 Double_t dphi = TMath::Abs(ltp->GetPhi() - phi);
437 if ( dphi > TMath::DegToRad()*180 ) dphi-=TMath::DegToRad()*180;
438 if ( dphi<dphiMin ) {
447 //_______________________________________________________________________
448 AliTPCLaserTracks::~AliTPCLaserTracks()
451 // standard destructor
454 // if ( fP ) delete fP;
455 if ( fXarr ) delete [] fXarr;
456 if ( fYarr ) delete [] fYarr;
457 if ( fZarr ) delete [] fZarr;