Some of the coding violations corrected
[u/mrichter/AliRoot.git] / TPC / AliTPCExBBShape.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 // AliTPCExBBShape class                                                  //
18 ////////////////////////////////////////////////////////////////////////////
19
20 #include <AliMagF.h>
21 #include "TGeoGlobalMagField.h"
22 #include "AliTPCcalibDB.h"
23 #include "AliTPCParam.h"
24 #include "AliLog.h"
25
26 #include "AliTPCExBBShape.h"
27
28 AliTPCExBBShape::AliTPCExBBShape()
29   : AliTPCCorrection("exb_bshape","ExB B-shape"),
30     fC1(0.),fC2(0.),
31     fBField(0)
32 {
33   //
34   // default constructor
35   //
36 }
37
38 AliTPCExBBShape::~AliTPCExBBShape() {
39   //
40   // virtual destructor
41   //
42 }
43
44 void AliTPCExBBShape::Init() {
45   //
46   // Initialization funtion
47   //
48   
49   AliMagF* magF= (AliMagF*)TGeoGlobalMagField::Instance()->GetField();
50   if (!magF) AliError("Magneticd field - not initialized");
51   Double_t bzField = magF->SolenoidField()/10.; //field in T
52   SetBField(magF);
53   AliTPCParam *param= AliTPCcalibDB::Instance()->GetParameters();
54   if (!param) AliError("Parameters - not initialized");
55   Double_t vdrift = param->GetDriftV()/1000000.; // [cm/us]   // From dataBase: to be updated: per second (ideally)
56   Double_t ezField = 400; // [V/cm]   // to be updated: never (hopefully)
57   Double_t wt = -10.0 * (bzField*10) * vdrift / ezField ; 
58   // Correction Terms for effective omegaTau; obtained by a laser calibration run
59   SetOmegaTauT1T2(wt,fT1,fT2);
60
61
62 }
63
64 void AliTPCExBBShape::Update(const TTimeStamp &/*timeStamp*/) {
65   //
66   // Update function 
67   //
68   AliMagF* magF= (AliMagF*)TGeoGlobalMagField::Instance()->GetField();
69   if (!magF) AliError("Magneticd field - not initialized");
70   Double_t bzField = magF->SolenoidField()/10.; //field in T
71   SetBField(magF);
72   AliTPCParam *param= AliTPCcalibDB::Instance()->GetParameters();
73   if (!param) AliError("Parameters - not initialized");
74   Double_t vdrift = param->GetDriftV()/1000000.; // [cm/us]   // From dataBase: to be updated: per second (ideally)
75   Double_t ezField = 400; // [V/cm]   // to be updated: never (hopefully)
76   Double_t wt = -10.0 * (bzField*10) * vdrift / ezField ; 
77   // Correction Terms for effective omegaTau; obtained by a laser calibration run
78   SetOmegaTauT1T2(wt,fT1,fT2);
79
80
81 }
82
83
84
85 void AliTPCExBBShape::GetCorrection(const Float_t x[],const Short_t roc,Float_t dx[]) {
86   //
87   // Calculates the space point corrections of the B field inperfections (B field shape) 
88   //
89
90   if (!fBField) {
91     for (Int_t j=0;j<3;++j) dx[j]=0.;
92     return;
93   }
94
95   const Double_t xStart[3]={ x[0], x[1], x[2] };
96   const Double_t xEnd[3]={ x[0],  x[1],  roc%36<18?fgkTPCZ0:-fgkTPCZ0 };
97
98   Double_t intBStart[3];
99   Double_t intBEnd[3];
100
101   fBField->GetTPCRatInt(xStart,intBStart);
102   fBField->GetTPCRatInt(xEnd,  intBEnd  );
103
104   const Float_t intBxOverBz=(intBEnd[0]-intBStart[0]);
105   const Float_t intByOverBz=(intBEnd[1]-intBStart[1]);
106   
107   dx[0]=fC2*intBxOverBz-fC1*intByOverBz;
108   dx[1]=fC1*intBxOverBz+fC2*intByOverBz;
109   dx[2]=0.;
110
111
112 }
113
114 void AliTPCExBBShape::GetBxAndByOverBz(const Float_t x[],const Short_t roc,Float_t BxByOverBz[]) {
115   //
116   // This function is purely for calibration purposes
117   // Returns the via AliMagF obtaind B field integrals  
118   // 
119
120   if (!fBField) {
121     for (Int_t j=0;j<3;++j) BxByOverBz[j]=0.;
122     return;
123   }
124
125   const Double_t xStart[3]={ x[0], x[1], x[2] };
126   const Double_t xEnd[3]={ x[0],  x[1],  roc%36<18?fgkTPCZ0:-fgkTPCZ0 };
127
128   Double_t intBStart[3];
129   Double_t intBEnd[3];
130
131   fBField->GetTPCRatInt(xStart,intBStart);
132   fBField->GetTPCRatInt(xEnd,  intBEnd  );
133
134   const Float_t intBxOverBz=(intBEnd[0]-intBStart[0]);
135   const Float_t intByOverBz=(intBEnd[1]-intBStart[1]);
136   
137   BxByOverBz[0]=intBxOverBz;
138   BxByOverBz[1]=intByOverBz;
139
140 }
141
142 void AliTPCExBBShape::Print(Option_t* option) const {
143   //
144   // Print function to check the settings (e.g. voltage offsets)
145   // option=="a" prints details of the B field settings and the 
146   // C0 and C1 coefficents (for calibration purposes)
147   //
148   TString opt = option; opt.ToLower();
149   printf("%s\t%s\n - B field settings:\n",GetTitle(),GetName());
150   fBField->Print(option);
151   //  printf(" - B field: X-Twist: %1.5lf rad, Y-Twist: %1.5lf rad \n",fBField->Print(option));
152   if (opt.Contains("a")) { // Print all details  
153     printf(" - T1: %1.4f, T2: %1.4f \n",fT1,fT2);
154     printf(" - C1: %1.4f, C2: %1.4f \n",fC1,fC2);
155   }    
156 }
157
158 Double_t AliTPCExBBShape::GetBFieldXYZ(Double_t gx, Double_t gy, Double_t gz, Int_t axisType){
159   //
160   // return B field at given x,y,z
161   // 
162   AliMagF* field = (AliMagF*)TGeoGlobalMagField::Instance()->GetField();
163   if (!field) return 0;
164   Double_t xyz[3]={gx,gy,gz};
165   Double_t bxyz[3]={0};
166   field->Field(xyz,bxyz);
167   //
168   Double_t r=TMath::Sqrt(gx*gx+gy*gy);
169   //  Double_t b=TMath::Sqrt(bxyz[0]*bxyz[0]+bxyz[1]*bxyz[1]);
170   if (axisType==0) {
171     return (xyz[0]*bxyz[1]-xyz[1]*bxyz[0])/(bxyz[2]*r);
172   }
173   if (axisType==1){
174     return (xyz[0]*bxyz[0]+xyz[1]*bxyz[1])/(bxyz[2]*r);
175   }
176   if (axisType==2) return bxyz[2];
177   if (axisType==3) return bxyz[0];
178   if (axisType==4) return bxyz[1];
179   if (axisType==5) return bxyz[2];
180   return bxyz[2];
181 }