-Double_t AliTPCCorrection::GetCorrSector(Double_t sector, Double_t localX, Double_t kZ, Int_t axisType, Int_t corrType){
+Double_t AliTPCCorrection::GetCorrSector(Double_t sector, Double_t r, Double_t kZ, Int_t axisType, Int_t corrType){
//
// calculate the correction at given position - check the geffCorr
//
AliTPCCorrection *corr = (AliTPCCorrection*)fgVisualCorrection->At(corrType);
if (!corr) return 0;
Double_t phi=sector*TMath::Pi()/9.;
- Double_t gx = localX*TMath::Cos(phi);
- Double_t gy = localX*TMath::Sin(phi);
- Double_t gz = localX*kZ;
+ Double_t gx = r*TMath::Cos(phi);
+ Double_t gy = r*TMath::Sin(phi);
+ Double_t gz = r*kZ;
Int_t nsector=(gz>0) ? 0:18;
//
//
void FastSimDistortedVertex(Double_t orgVertex[3], Int_t nTracks, AliESDVertex &aV, AliESDVertex &avOrg, AliESDVertex &cV, AliESDVertex &cvOrg, TTreeSRedirector * const pcstream, Double_t etaCuts);
static void AddVisualCorrection(AliTPCCorrection* corr, Int_t position);
- static Double_t GetCorrSector(Double_t sector, Double_t localX, Double_t kZ, Int_t axisType, Int_t corrType=0);
+ static Double_t GetCorrSector(Double_t sector, Double_t r, Double_t kZ, Int_t axisType, Int_t corrType=0);
static Double_t GetCorrXYZ(Double_t gx, Double_t gy, Double_t gz, Int_t axisType, Int_t corrType=0);
protected:
axisName[1]="sector"; axisTitle[1]="Sector Number";
binsTrack[1]=180; xminTrack[1]=0; xmaxTrack[1]=18;
//
- axisName[2]="localX"; axisTitle[2]="x (cm)";
+ axisName[2]="R"; axisTitle[2]="r (cm)";
binsTrack[2]=53; xminTrack[2]=85.; xmaxTrack[2]=245.;
//
//
trackOut.GetXYZ(xyz);
resVector[1]= 9.*TMath::ATan2(xyz[1],xyz[0])/TMath::Pi();
if (resVector[1]<0) resVector[1]+=18;
- resVector[2]= cl->GetX();
+ resVector[2]= TMath::Sqrt(cl->GetX()*cl->GetX()+cl->GetY()*cl->GetY());
resVector[3]= cl->GetZ()/cl->GetX();
//
resVector[0]= cl->GetY()-trackOut.GetY();
trackIn.GetXYZ(xyz);
resVector[1]= 9.*TMath::ATan2(xyz[1],xyz[0])/TMath::Pi();
if (resVector[1]<0) resVector[1]+=18;
- resVector[2]= cl->GetX();
+ resVector[2]= TMath::Sqrt(cl->GetX()*cl->GetX()+cl->GetY()*cl->GetY());
resVector[3]= cl->GetZ()/cl->GetX();
//
resVector[0]= cl->GetY()-trackIn.GetY();
if (fCurrentTrack->IsOn(0x4)&&TMath::Abs(vImpact[0])<1&&TMath::Abs(vImpact[1])<1){
Double_t resVector[5];
resVector[1]= 9.*gphi/TMath::Pi();
- resVector[2]= c->GetX();
+ resVector[2]= TMath::Sqrt(c->GetX()*c->GetX()+c->GetY()*c->GetY());
resVector[3]= c->GetZ()/c->GetX();
//
//