1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
16 ////////////////////////////////////////////////////////////////////////////
17 // AliTPCExBTwist class //
18 ////////////////////////////////////////////////////////////////////////////
22 #include "TGeoGlobalMagField.h"
23 #include "AliTPCcalibDB.h"
24 #include "AliTPCParam.h"
27 #include "AliTPCExBTwist.h"
29 AliTPCExBTwist::AliTPCExBTwist()
30 : AliTPCCorrection("exb_twist","ExB twist"),
32 fXTwist(0.),fYTwist(0.)
35 // default constructor
39 AliTPCExBTwist::~AliTPCExBTwist() {
45 Bool_t AliTPCExBTwist::AddCorrectionCompact(AliTPCCorrection* corr, Double_t weight){
47 // Add correction and make them compact
49 // - origin of distortion/correction are additive
50 // - only correction ot the same type supported ()
52 AliError("Zerro pointer - correction");
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]
64 void AliTPCExBTwist::Init() {
66 // Initialization funtion
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);
83 void AliTPCExBTwist::Update(const TTimeStamp &/*timeStamp*/) {
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);
103 void AliTPCExBTwist::GetCorrection(const Float_t x[],const Short_t roc,Float_t dx[]) {
105 // Calculates the correction of a mismatch between the E and B field axis
108 const Float_t zstart=x[2];
109 const Float_t zend =(roc%36<18?fgkTPCZ0:-fgkTPCZ0);
110 const Float_t zdrift=zstart-zend;
112 dx[0]=(fC2*fXTwist-fC1*fYTwist)*zdrift;
113 dx[1]=(fC1*fXTwist+fC2*fYTwist)*zdrift;
117 void AliTPCExBTwist::Print(const Option_t* option) const {
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
123 TString opt = option; opt.ToLower();
124 printf("%s\n",GetTitle());
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);