During simulation: fill STU region w/ non null time sums
[u/mrichter/AliRoot.git] / TPC / AliTPCExBTwist.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 // AliTPCExBTwist class                                                   //
18 ////////////////////////////////////////////////////////////////////////////
19
20
21 #include "AliMagF.h"
22 #include "TGeoGlobalMagField.h"
23 #include "AliTPCcalibDB.h"
24 #include "AliTPCParam.h"
25 #include "AliLog.h"
26
27 #include "AliTPCExBTwist.h"
28
29 AliTPCExBTwist::AliTPCExBTwist()
30   : AliTPCCorrection("exb_twist","ExB twist"),
31     fC1(0.),fC2(0.),
32     fXTwist(0.),fYTwist(0.)
33 {
34   //
35   // default constructor
36   //
37 }
38
39 AliTPCExBTwist::~AliTPCExBTwist() {
40   //
41   // default destructor
42   //
43 }
44
45
46
47 void AliTPCExBTwist::Init() {
48   //
49   // Initialization funtion
50   //
51   
52   AliMagF* magF= (AliMagF*)TGeoGlobalMagField::Instance()->GetField();
53   if (!magF) AliError("Magneticd field - not initialized");
54   Double_t bzField = magF->SolenoidField()/10.; //field in T
55   AliTPCParam *param= AliTPCcalibDB::Instance()->GetParameters();
56   if (!param) AliError("Parameters - not initialized");
57   Double_t vdrift = param->GetDriftV()/1000000.; // [cm/us]   // From dataBase: to be updated: per second (ideally)
58   Double_t ezField = 400; // [V/cm]   // to be updated: never (hopefully)
59   Double_t wt = -10.0 * (bzField*10) * vdrift / ezField ; 
60   // Correction Terms for effective omegaTau; obtained by a laser calibration run
61   SetOmegaTauT1T2(wt,fT1,fT2);
62
63
64 }
65
66 void AliTPCExBTwist::Update(const TTimeStamp &/*timeStamp*/) {
67   //
68   // Update function 
69   //
70   AliMagF* magF= (AliMagF*)TGeoGlobalMagField::Instance()->GetField();
71   if (!magF) AliError("Magneticd field - not initialized");
72   Double_t bzField = magF->SolenoidField()/10.; //field in T
73   AliTPCParam *param= AliTPCcalibDB::Instance()->GetParameters();
74   if (!param) AliError("Parameters - not initialized");
75   Double_t vdrift = param->GetDriftV()/1000000.; // [cm/us]   // From dataBase: to be updated: per second (ideally)
76   Double_t ezField = 400; // [V/cm]   // to be updated: never (hopefully)
77   Double_t wt = -10.0 * (bzField*10) * vdrift / ezField ; 
78   // Correction Terms for effective omegaTau; obtained by a laser calibration run
79   SetOmegaTauT1T2(wt,fT1,fT2);
80
81
82 }
83
84
85
86 void AliTPCExBTwist::GetCorrection(const Float_t x[],const Short_t roc,Float_t dx[]) {
87   //
88   // Calculates the correction of a mismatch between the E and B field axis
89   // 
90   
91   const Float_t zstart=x[2];
92   const Float_t zend  =(roc%36<18?fgkTPCZ0:-fgkTPCZ0);
93   const Float_t zdrift=zstart-zend;
94   
95   dx[0]=(fC2*fXTwist-fC1*fYTwist)*zdrift;
96   dx[1]=(fC1*fXTwist+fC2*fYTwist)*zdrift;
97   dx[2]=0.;
98 }
99
100 void AliTPCExBTwist::Print(const Option_t* option) const {
101   //
102   // Print function to check the settings (e.g. the twist in the X direction)
103   // option=="a" prints the C0 and C1 coefficents for calibration purposes
104   //
105
106   TString opt = option; opt.ToLower();
107   printf("%s\n",GetTitle());
108   
109   printf(" - Twist settings: X-Twist: %1.5f rad, Y-Twist: %1.5f rad \n",fXTwist,fYTwist);
110   if (opt.Contains("a")) { // Print all details
111     printf(" - T1: %1.4f, T2: %1.4f \n",fT1,fT2);
112     printf(" - C1: %1.4f, C2: %1.4f \n",fC1,fC2);
113   }    
114  
115  
116 }