TPC/TPCbase/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 Bool_t AliTPCExBTwist::AddCorrectionCompact(AliTPCCorrection* corr, Double_t weight){
  46   //
  47   // Add correction  and make them compact
  48   // Assumptions:
  49   //  - origin of distortion/correction are additive
  50   //  - only correction ot the same type supported ()
  51   if (corr==NULL) {
  52     AliError("Zerro pointer - correction");
  53     return kFALSE;
  54   }
  55   AliTPCExBTwist * corrC = dynamic_cast< AliTPCExBTwist*>(corr);
  56   if (corrC == NULL) return kFALSE;
  57   fXTwist+=weight*corrC->fXTwist;        // Twist of E to B field in X-Z [rad]
  58   fYTwist+=weight*corrC->fYTwist;        // Twist of E to B field in Y-Z [rad]
  59   return kTRUE;
  60 }
  61
  62
  63
  64 void AliTPCExBTwist::Init() {
  65   //
  66   // Initialization funtion
  67   //
  68
  69   AliMagF* magF= (AliMagF*)TGeoGlobalMagField::Instance()->GetField();
  70   if (!magF) AliError("Magneticd field - not initialized");
  71   Double_t bzField = magF->SolenoidField()/10.; //field in T
  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 void AliTPCExBTwist::Update(const TTimeStamp &/*timeStamp*/) {
  84   //
  85   // Update function
  86   //
  87   AliMagF* magF= (AliMagF*)TGeoGlobalMagField::Instance()->GetField();
  88   if (!magF) AliError("Magneticd field - not initialized");
  89   Double_t bzField = magF->SolenoidField()/10.; //field in T
  90   AliTPCParam *param= AliTPCcalibDB::Instance()->GetParameters();
  91   if (!param) AliError("Parameters - not initialized");
  92   Double_t vdrift = param->GetDriftV()/1000000.; // [cm/us]   // From dataBase: to be updated: per second (ideally)
  93   Double_t ezField = 400; // [V/cm]   // to be updated: never (hopefully)
  94   Double_t wt = -10.0 * (bzField*10) * vdrift / ezField ;
  95   // Correction Terms for effective omegaTau; obtained by a laser calibration run
  96   SetOmegaTauT1T2(wt,fT1,fT2);
  97
  98
  99 }
 100
 101
 102
 103 void AliTPCExBTwist::GetCorrection(const Float_t x[],const Short_t roc,Float_t dx[]) {
 104   //
 105   // Calculates the correction of a mismatch between the E and B field axis
 106   //
 107
 108   const Float_t zstart=x[2];
 109   const Float_t zend  =(roc%36<18?fgkTPCZ0:-fgkTPCZ0);
 110   const Float_t zdrift=zstart-zend;
 111
 112   dx[0]=(fC2*fXTwist-fC1*fYTwist)*zdrift;
 113   dx[1]=(fC1*fXTwist+fC2*fYTwist)*zdrift;
 114   dx[2]=0.;
 115 }
 116
 117 void AliTPCExBTwist::Print(const Option_t* option) const {
 118   //
 119   // Print function to check the settings (e.g. the twist in the X direction)
 120   // option=="a" prints the C0 and C1 coefficents for calibration purposes
 121   //
 122
 123   TString opt = option; opt.ToLower();
 124   printf("%s\n",GetTitle());
 125
 126   printf(" - Twist settings: X-Twist: %1.5f rad, Y-Twist: %1.5f rad \n",fXTwist,fYTwist);
 127   if (opt.Contains("a")) { // Print all details
 128     printf(" - T1: %1.4f, T2: %1.4f \n",fT1,fT2);
 129     printf(" - C1: %1.4f, C2: %1.4f \n",fC1,fC2);
 130   }
 131
 132
 133 }