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("cAliTPCExBTwist","cAliTPCExBTwist",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