#include "TMath.h"
#include "TTreeStream.h"
-#include "AliFieldMap.h"
#include "AliMagF.h"
#include "AliTPCExBExact.h"
ClassImp(AliTPCExBExact)
const Double_t AliTPCExBExact::fgkEM=1.602176487e-19/9.10938215e-31;
-const Double_t AliTPCExBExact::fgkDriftField=40.e3;
+const Double_t AliTPCExBExact::fgkDriftField=-40.e3;
AliTPCExBExact::AliTPCExBExact()
: fDriftVelocity(0),
- fkMap(0),fkField(0),fkN(0),
+ //fkMap(0),
+ fkField(0),fkN(0),
fkNX(0),fkNY(0),fkNZ(0),
fkXMin(-250.),fkXMax(250.),fkYMin(-250.),fkYMax(250.),
fkZMin(-250.),fkZMax(250.),
Double_t driftVelocity,
Int_t nx,Int_t ny,Int_t nz,Int_t n)
: fDriftVelocity(driftVelocity),
- fkMap(0),fkField(bField),fkN(n),
+ //fkMap(0),
+ fkField(bField),fkN(n),
fkNX(nx),fkNY(ny),fkNZ(nz),
fkXMin(-250.),fkXMax(250.),fkYMin(-250.),fkYMax(250.),
fkZMin(-250.),fkZMax(250.),
CreateLookupTable();
}
+/*
AliTPCExBExact::AliTPCExBExact(const AliFieldMap *bFieldMap,
Double_t driftVelocity,Int_t n)
: fDriftVelocity(driftVelocity),
CreateLookupTable();
}
+*/
AliTPCExBExact::~AliTPCExBExact() {
//
// corrected[2]=position[2];
}
+/*
void AliTPCExBExact::TestThisBeautifulObject(const AliFieldMap *bFieldMap,
const char* fileName) {
//
fkField=0;
TestThisBeautifulObjectGeneric(fileName);
}
+*/
void AliTPCExBExact::TestThisBeautifulObject(const AliMagF *bField,
const char* fileName) {
// Have a look at the common part "TestThisBeautifulObjectGeneric".
//
fkField=bField;
- fkMap=0;
+ //fkMap=0;
TestThisBeautifulObjectGeneric(fileName);
}
//
// Helper function returning the B field in SI units (T).
//
- Float_t xm[3];
+ Double_t xm[3];
// the beautiful m to cm (and the ugly "const_cast") and Double_t
// to Float_t read the NRs introduction!:
for (int i=0;i<3;++i) xm[i]=x[i]*100.;
- Float_t bf[3];
- if (fkMap!=0)
- fkMap->Field(xm,bf);
- else
- fkField->Field(xm,bf);
+ Double_t bf[3];
+ //if (fkMap!=0)
+ // fkMap->Field(xm,bf);
+ //else
+ ((AliMagF*)fkField)->Field(xm,bf);
for (int i=0;i<3;++i) b[i]=bf[i]/10.;
}