]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TPC/AliTPCtrackerParam.cxx
New task to check TPC-ITS track prolongation eff with cosmics
[u/mrichter/AliRoot.git] / TPC / AliTPCtrackerParam.cxx
index 022fa1eeb326878dc0351c705e857879acbfe85e..60983c3e8b8d95ebadda013e92694f8c1b65a1ff 100644 (file)
@@ -73,6 +73,7 @@
 #include <TLine.h>
 #include <TMatrixD.h>
 #include <TParticle.h>
+#include <TDatabasePDG.h>
 #include <TStyle.h>
 #include <TSystem.h>
 #include <TFile.h>
@@ -276,9 +277,7 @@ Int_t AliTPCtrackerParam::Init() {
     
   } else printf("\n ! Creating ALL TRUE tracks at TPC inner radius !\n\n");
 
-
-  // Check if value in the galice file is equal to selected one (fBz)
-  AliMagF *fiel = (AliMagF*)gAlice->Field();
+  AliMagF *fiel = (AliMagF*)TGeoGlobalMagField::Instance()->GetField();
   Double_t fieval=TMath::Abs((Double_t)fiel->SolenoidField()/10.);
   printf("Magnetic field is %6.2f Tesla\n",fieval);
   if(fBz!=fieval) {
@@ -288,9 +287,6 @@ Int_t AliTPCtrackerParam::Init() {
     return 1;
   }
 
-  // Set the conversion constant between curvature and Pt
-  AliTracker::SetFieldMap(fiel,kTRUE);
-
   return 0;
 }
 //-----------------------------------------------------------------------------
@@ -1113,7 +1109,7 @@ void AliTPCtrackerParam::CompareTPCtracks(
   TFile *galiceFile = TFile::Open(galiceName);
 
   // get the AliRun object
-  AliRun *gAlice = (AliRun*)galiceFile->Get("gAlice");
+  AliRun *lAlice = (AliRun*)galiceFile->Get("gAlice");
 
 
   // create the tree for comparison results
@@ -1164,7 +1160,7 @@ void AliTPCtrackerParam::CompareTPCtracks(
     sprintf(tname,"TreeT_TPC_%d",evt);
     
     // particles from TreeK
-    const Int_t knparticles = gAlice->GetEvent(evt);
+    const Int_t knparticles = lAlice->GetEvent(evt);
 
     Int_t *kalLab = new Int_t[knparticles];
     for(Int_t i=0; i<knparticles; i++) kalLab[i] = -1; 
@@ -1223,7 +1219,7 @@ void AliTPCtrackerParam::CompareTPCtracks(
       geatree->GetEvent(j);
       
       label = geatrack->GetLabel();
-      part = (TParticle*)gAlice->GetMCApp()->Particle(label);
+      part = (TParticle*)lAlice->GetMCApp()->Particle(label);
       
       // use only injected tracks with fixed values of pT
       ptgener = part->Pt();
@@ -1373,7 +1369,7 @@ void AliTPCtrackerParam::CompareTPCtracks(
   WriteEffs(tpceffrootName);
 
   // delete AliRun object
-  delete gAlice; gAlice=0;
+  delete lAlice; lAlice=0;
   
   // close all input files
   kalFile->Close();