Put back the two bellows in front of the absorber.
[u/mrichter/AliRoot.git] / TPC / AliDetectorParam.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.1.4.2  2000/04/10 11:36:13  kowal2
19
20 New Detector parameters handling class
21
22 */
23
24 ///////////////////////////////////////////////////////////////////////
25 //  Paramter class for AliDetector                                   //
26 //                                                                   //
27 //  Origin:  Marian Ivanov, Uni. of Bratislava, ivanov@fmph.uniba.sk // 
28 //                                                                   //  
29 ///////////////////////////////////////////////////////////////////////
30
31 #include "TObject.h"
32 #include "AliDetectorParam.h"
33
34
35
36
37 Float_t * AliDetectorParam::GetAnglesAccMomentum(Float_t *x, Int_t * index, Float_t *momentum, Float_t *angle)
38 {
39   //
40   //calculates deflection angle of particle with longitudinal
41   //longitudinal  momentum[0] and transversal momentum momentum[1]
42   //at position (x,y,z) = (x[0],x[1],x[2]) 
43   //angle[0] - deep angle
44   //angle[1] - magnetic deflection angle 
45   if (momentum==0) {
46     Float_t rtotal =TMath::Sqrt(x[0]*x[0]+x[1]*x[1]);
47     if (rtotal==0) angle[0]=0;
48     else    
49       angle[0] = TMath::ATan(x[2]/rtotal);
50     angle[1]=0;
51     return angle;
52   }
53   Float_t mtotal =TMath::Sqrt(momentum[0]*momentum[0]+momentum[1]*momentum[1]);
54   if (mtotal==0) {
55     angle[0]= 0;
56     angle[1]=0;
57     return angle;
58   }
59   angle[0]= TMath::ATan(momentum[2]/mtotal);
60   Float_t radius1 = TMath::Sqrt(x[0]*x[0]+x[1]*x[1]); //axial symetry in z
61   Float_t radius2 = 1000*mtotal/(3*fBField);
62   if (radius1<radius2)
63     angle[1]= TMath::ASin(radius1/radius2);
64   else 
65     angle[1]=0;
66   return angle;
67
68
69
70
71
72
73 ClassImp(AliDetectorParam)