Fixes for coding violations
[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 /* $Id$ */
17
18 ///////////////////////////////////////////////////////////////////////
19 //  Paramter class for AliDetector                                   //
20 //                                                                   //
21 //  Origin:  Marian Ivanov, Uni. of Bratislava, ivanov@fmph.uniba.sk // 
22 //                                                                   //  
23 ///////////////////////////////////////////////////////////////////////
24
25 #include <TMath.h>
26 #include <TObject.h>
27
28 #include "AliDetectorParam.h"
29
30
31 AliDetectorParam::AliDetectorParam()
32                    :TNamed(),
33                     fBField(0.),
34                     fNPrimLoss(0.),
35                     fNTotalLoss(0.)
36 {
37   //
38   //  default constructor
39   //
40 }
41
42 Float_t * AliDetectorParam::GetAnglesAccMomentum(Float_t *x, Int_t * /*index*/, Float_t *momentum, Float_t *angle)
43 {
44   //
45   //calculates deflection angle of particle with longitudinal
46   //longitudinal  momentum[0] and transversal momentum momentum[1]
47   //at position (x,y,z) = (x[0],x[1],x[2]) 
48   //angle[0] - deep angle
49   //angle[1] - magnetic deflection angle 
50   if (momentum==0) {
51     Float_t rtotal =TMath::Sqrt(x[0]*x[0]+x[1]*x[1]);
52     if (rtotal==0) angle[0]=0;
53     else    
54       angle[0] = TMath::ATan(x[2]/rtotal);
55     angle[1]=0;
56     return angle;
57   }
58   Float_t mtotal =TMath::Sqrt(momentum[0]*momentum[0]+momentum[1]*momentum[1]);
59   if (mtotal==0) {
60     angle[0]= 0;
61     angle[1]=0;
62     return angle;
63   }
64   angle[0]= TMath::ATan(momentum[2]/mtotal);
65   Float_t radius1 = TMath::Sqrt(x[0]*x[0]+x[1]*x[1]); //axial symetry in z
66   Float_t radius2 = 1000*mtotal/(3*fBField);
67   if (radius1<radius2)
68     angle[1]= TMath::ASin(radius1/radius2);
69   else 
70     angle[1]=0;
71   return angle;
72
73
74
75
76
77
78 ClassImp(AliDetectorParam)