Bugfix, seems more like trial and error.

[u/mrichter/AliRoot.git] / JETAN / AliJetParticlesReaderHLT.cxx

diff --git a/JETAN/AliJetParticlesReaderHLT.cxx b/JETAN/AliJetParticlesReaderHLT.cxx
index 967b43c528e6fce65ac7f900fa7241b0d3fdd2f8..28cbfccc8704e853605350312022bedb4f7287c1 100644 (file)
--- a/JETAN/AliJetParticlesReaderHLT.cxx
+++ b/JETAN/AliJetParticlesReaderHLT.cxx
@@ -110,12 +110,13 @@ Int_t AliJetParticlesReaderHLT::ReadESD(AliESD* esd)
    }
 
   fEventParticles->SetVertex(vertexpos[0],vertexpos[1],vertexpos[2]);
+  //cout << vertexpos[0] << " " << vertexpos[1] << " " << vertexpos[2] << endl;
 
   AliL3Track l3;
   AliL3Vertex v;
   v.SetX(vertexpos[0]);
-  v.SetY(vertexpos[0]);
-  v.SetZ(vertexpos[0]);
+  v.SetY(vertexpos[1]);
+  v.SetZ(vertexpos[2]);
   Double_t xc=0.,yc=0.,zc=0.;
   for (Int_t i = 0;i<ntr; i++) {
     AliESDHLTtrack *kesdtrack;
@@ -145,6 +146,7 @@ Int_t AliJetParticlesReaderHLT::ReadESD(AliESD* esd)
 
     if(fTrackerType){
       //if(!kesdtrack->ComesFromMainVertex()) continue;
+      //cout << kesdtrack->GetPx() << " " << kesdtrack->GetPy() << " " << kesdtrack->GetPz() << endl;
       l3.SetFirstPoint(kesdtrack->GetFirstPointX(),kesdtrack->GetFirstPointY(),kesdtrack->GetFirstPointZ());
       l3.SetLastPoint(kesdtrack->GetLastPointX(),kesdtrack->GetLastPointY(),kesdtrack->GetLastPointZ());
       l3.SetCharge(kesdtrack->GetCharge());
@@ -155,10 +157,12 @@ Int_t AliJetParticlesReaderHLT::ReadESD(AliESD* esd)
       l3.GetClosestPoint(&v,xc,yc,zc);
       if(TMath::Abs(zc)>10.) continue;
       l3.SetFirstPoint(vertexpos[0],vertexpos[1],vertexpos[2]);
+      //l3.CalculateHelix();
       l3.UpdateToFirstPoint();
       px=l3.GetPx();
       py=l3.GetPy();
       pz=l3.GetPz();
+      //cout << px << " " << py << " " << pz << endl;
     }
 
     const Float_t kpt=TMath::Sqrt(px*px+py*py);