#include <TLine.h>
#include <TMatrixD.h>
#include <TParticle.h>
+#include <TDatabasePDG.h>
#include <TStyle.h>
#include <TSystem.h>
#include <TFile.h>
} 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) {
return 1;
}
- // Set the conversion constant between curvature and Pt
- AliTracker::SetFieldMap(fiel,kTRUE);
-
return 0;
}
//-----------------------------------------------------------------------------
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
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;
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();
WriteEffs(tpceffrootName);
// delete AliRun object
- delete gAlice; gAlice=0;
+ delete lAlice; lAlice=0;
// close all input files
kalFile->Close();