Double_t &fC43=fC[13];
Double_t &fC44=fC[14];
- Double_t p2=(1.+ fP3*fP3)/(fP4*fP4);
+ Double_t p=GetP();
+ Double_t p2=p*p;
Double_t beta2=p2/(p2 + mass*mass);
d*=TMath::Sqrt((1.+ fP3*fP3)/(1.- fP2*fP2));
const Double_t kB2C=0.299792458e-3;
const Double_t kAlmost0Field=1.e-13;
-const Double_t kVeryBigConvConst=1/kB2C/kAlmost0Field;
-const Double_t kMostProbableMomentum=0.35;
+const Double_t kMostProbablePt=0.35;
class AliESDVertex;
fC[1] = 0.; fC[2]*= s2;
fC[3] = 0.; fC[4] = 0.; fC[5]*= s2;
fC[6] = 0.; fC[7] = 0.; fC[8] = 0.; fC[9]*= s2;
- fC[10]= 0.; fC[11]= 0.; fC[12]= 0.; fC[13]= 0.; fC[14]*=10.;
+ fC[10]= 0.; fC[11]= 0.; fC[12]= 0.; fC[13]= 0.; fC[14]*=s2;
}
const Double_t *GetParameter() const {return fP;}
#include <TClass.h>
#include <TMath.h>
+#include "AliMagF.h"
#include "AliTracker.h"
#include "AliCluster.h"
#include "AliKalmanTrack.h"
Bool_t AliTracker::fgUniformField=kTRUE;
-Double_t AliTracker::fgBz=0.;
+Double_t AliTracker::fgBz=kAlmost0Field;
const AliMagF *AliTracker::fgkFieldMap=0;
ClassImp(AliTracker)
//Double_t bz=-b[2];
Double_t bz=-map->SolenoidField();
- fgBz=TMath::Sign(1e-13,bz) + bz;
+ fgBz=TMath::Sign(kAlmost0Field,bz) + bz;
}
c->Use();
}
}
+
+Double_t AliTracker::GetBz(Float_t *r) {
+ //------------------------------------------------------------------
+ // Returns Bz (kG) at the point "r" .
+ //------------------------------------------------------------------
+ Float_t b[3]; fgkFieldMap->Field(r,b);
+ Double_t bz=-Double_t(b[2]);
+ return (TMath::Sign(kAlmost0Field,bz) + bz);
+}
// Origin: Iouri Belikov, CERN, Jouri.Belikov@cern.ch
//-------------------------------------------------------------------------
#include <TObject.h>
-#include "AliMagF.h"
+class AliMagF;
class AliCluster;
class TTree;
class AliKalmanTrack;
static void SetFieldMap(const AliMagF* map, Bool_t uni);
static const AliMagF *GetFieldMap() {return fgkFieldMap;}
- static Double_t GetBz(Float_t *r) {
- Float_t b[3]; fgkFieldMap->Field(r,b);
- Double_t bz=-Double_t(b[2]);
- return (TMath::Sign(1e-13,bz) + bz);
- }
+ static Double_t GetBz(Float_t *r);
static Double_t GetBz(Double_t *r) {
Float_t rr[]={r[0],r[1],r[2]};
return GetBz(rr);