Fixes for coding violations
[u/mrichter/AliRoot.git] / TPC / AliTPCExBTwist.h
1 #ifndef ALI_TPC_EX_B_TWIST_H
2 #define ALI_TPC_EX_B_TWIST_H
3
4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5  * See cxx source for full Copyright notice                               */
6
7 // _________________________________________________________________
8 //
9 // Begin_Html
10 //   <h2> AliTPCExBTwist class    </h2>       
11 //   The class calculates the space point distortions due to a mismatch     
12 //   of the E and B field axis (original code from STAR)                    
13 //   <p>
14 //   It assumes a homogeneous magnetic and electric field. The calculations performed 
15 //   within this class can be written in a compact form as well (compare notes of Jim Thomas)
16 //   <p>
17 //   The inputs to this class are the twist angles in the x and y plane ($\theta_x,\theta_y$) 
18 //   in cartesian coordinates. They can be set via the member functions SetXTwist and SetYTwist.
19 // End_Html
20 //
21 // Begin_Macro(source) 
22 //   {
23 //   gROOT->SetStyle("Plain"); gStyle->SetPalette(1);
24 //   TCanvas *c2 = new TCanvas("c2","c2",500,300); 
25 //   AliTPCExBTwist twist;
26 //   twist.SetXTwist(0.001);  // x angle in [rad]
27 //   twist.SetXTwist(0.0005); // y angle in [rad]
28 //   twist.SetOmegaTauT1T2(0.32,1,1);
29 //   twist.CreateHistoDRPhiinXY(1.)->Draw("surf2"); // A side
30 //   return c2;
31 //   } 
32 // End_Macro
33 //
34 // Begin_Html
35 //   <p>
36 // Date: 27/04/2010  <br>
37 // Authors: Jim Thomas, Magnus Mager, Stefan Rossegger 
38 // End_Html 
39 // _________________________________________________________________
40
41 #include "AliTPCCorrection.h"
42
43 class AliTPCExBTwist : public AliTPCCorrection {
44 public:
45   AliTPCExBTwist();
46   virtual ~AliTPCExBTwist();
47
48   // initialization and update functions
49   virtual void Init();
50   virtual void Update(const TTimeStamp &timeStamp);
51
52
53   // common setters and getters for ExB
54   virtual void SetOmegaTauT1T2(Float_t omegaTau,Float_t t1,Float_t t2) {
55     fT1=t1; fT2=t2; 
56     const Float_t wt1=t1*omegaTau;     fC1=wt1/(1.+wt1*wt1);
57     const Float_t wt2=t2*omegaTau;     fC2=wt2*wt2/(1.+wt2*wt2);
58   };
59   void SetC1C2(Float_t c1,Float_t c2) {fC1=c1;fC2=c2;} // CAUTION: USE WITH CARE
60   Float_t GetC1() const {return fC1;}
61   Float_t GetC2() const {return fC2;}
62
63   // setters and getters for twist
64   void SetXTwist(Float_t xTwist) {fXTwist=xTwist;}
65   void SetYTwist(Float_t yTwist) {fYTwist=yTwist;}
66   Float_t GetXTwist() const {return fXTwist;}
67   Float_t GetYTwist() const {return fYTwist;}
68
69   virtual void Print(const Option_t* option="") const;
70
71 protected:
72   virtual void GetCorrection(const Float_t x[],const Short_t roc,Float_t dx[]);
73
74 private:
75   Float_t fC1; // coefficient C1                 (compare Jim Thomas's notes for definitions)
76   Float_t fC2; // coefficient C2                 (compare Jim Thomas's notes for definitions)
77
78   Float_t fXTwist;               // Twist of E to B field in X-Z [rad]
79   Float_t fYTwist;               // Twist of E to B field in Y-Z [rad]
80
81   ClassDef(AliTPCExBTwist,1);
82 };
83
84 #endif