#include <AliITSVertexer3DTapan.h>
#include <AliITSRecPoint.h>
-#include <AliITSgeom.h>
+#include <AliITSgeomTGeo.h>
#include <AliESDVertex.h>
ClassImp(AliITSVertexer3DTapan)
// local cluster X"-axis
//
- // Double_t rot[9]; fITSgeom->GetRotMatrix(i,rot);
- Int_t lay,lad,det; fITSgeom->GetModuleId(i,lay,lad,det);
+ Int_t lay,lad,det; AliITSgeomTGeo::GetModuleId(i,lay,lad,det);
if (lay>2) break; //load the SPD clusters only
- /*
- Float_t tx,ty,tz; fITSgeom->GetTrans(lay,lad,det,tx,ty,tz);
-
- Double_t alpha=TMath::ATan2(rot[1],rot[0])+TMath::Pi();
- Double_t phi=TMath::Pi()/2+alpha;
-
- if (lay==1) phi+=TMath::Pi();
- Double_t cp=TMath::Cos(phi), sp=TMath::Sin(phi);
- Double_t r=tx*cp+ty*sp;
- */
-
Int_t ncl=clusters->GetEntriesFast();
- Float_t hPhi;
+ Float_t hPhi=0.;
while (ncl--) {
AliITSRecPoint *c=(AliITSRecPoint*)clusters->UncheckedAt(ncl);
Float_t pos[3];
//
LoadClusters(cTree);
- Double_t pos[3], postemp[3], sigpos[3];
+ Double_t pos[3], postemp[3];
+ Double_t sigpos[3]={0.,0.,0.};
Int_t ncontr, ncontrtemp;
Float_t cuts[3];
Int_t vtxstatus=0;
}
AliInfo(Form("Final step: %d %f %f %f st=%d",ncontr,pos[0],pos[1],pos[2],vtxstatus));
- return new AliESDVertex(pos,sigpos,(Double_t)vtxstatus,ncontr,"AliITSVertexer3DTapan");
+ Double_t covma[6]={0.,0.,0.,0.,0.,0.};
+ covma[0]=sigpos[0];
+ covma[2]=sigpos[1];
+ covma[5]=sigpos[2];
+ return new AliESDVertex(pos,covma,(Double_t)vtxstatus,ncontr,"AliITSVertexer3DTapan");
}