first steps to set the covariance matrix from the errors calculated in conformal...
authorrichterm <richterm@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 9 Jul 2008 18:20:46 +0000 (18:20 +0000)
committerrichterm <richterm@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 9 Jul 2008 18:20:46 +0000 (18:20 +0000)
HLT/TPCLib/AliHLTTPCConfMapFit.cxx
HLT/TPCLib/AliHLTTPCConfMapTrack.cxx
HLT/TPCLib/AliHLTTPCTrack.cxx
HLT/TPCLib/AliHLTTPCTrack.h

index 3efbdaa..be07e55 100644 (file)
@@ -414,6 +414,8 @@ Int_t AliHLTTPCConfMapFit::FitCircle()
   fTrack->SetCenterX(acent);
   fTrack->SetCenterY(bcent);
 
+  //set error for pT and Y. psi, Z and Tgl are set.
+
   //
 //    Get errors from fast fit
 //
index 7e413b7..58c457c 100644 (file)
@@ -232,6 +232,7 @@ void AliHLTTPCConfMapTrack::Fill(AliHLTTPCVertex *vertex,Double_t max_Dca)
       //fPt = (Double_t)(AliHLTTPCTransform::GetBFieldValue() * fRadius);
       Double_t pt = (Double_t)(AliHLTTPCTransform::GetBFieldValue() * GetRadius());
       SetPt(pt);
+      //Set Error here?????
       
       if(GetPt() > max_Dca) //go for fit of helix in real space
        {
index 60edf9b..3382986 100644 (file)
@@ -57,6 +57,7 @@ AliHLTTPCTrack::AliHLTTPCTrack()
   fPterr(0),
   fPsierr(0),
   fZ0err(0),
+  fY0err(0),
   fTanlerr(0),
 
   fPhi0(0),
@@ -657,11 +658,11 @@ int AliHLTTPCTrack::Convert2AliKalmanTrack()
   
   //covariance matrix
   Double_t cov[15]={
-    0.,
-    0.,  0.,
-    0.,  0.,  0.,
-    0.,  0.,  0.,  0.,
-    0.,  0.,  0.,  0.,  0.
+    GetY0err(),                         //Error in Y (Y and X are the same)
+    0.,  GetZ0err(),                    //Error in Z
+    0.,  0.,  GetPsierr(),              //Error for Psi
+    0.,  0.,  0.,  GetTglerr(),         //Error for Tgl
+    0.,  0.,  0.,  0.,  GetPterr()      //Error for Pt
   };
 
   Int_t nCluster = GetNHits();
index c9f2388..3839e0a 100644 (file)
@@ -97,7 +97,9 @@ class AliHLTTPCTrack : public AliTPCtrack {
   Double_t GetPterr() const {return fPterr;}
   Double_t GetPsierr() const {return fPsierr;}
   Double_t GetTglerr() const {return fTanlerr;}
-  
+  Double_t GetZ0err() const {return fZ0err;}
+  Double_t GetY0err() const {return fY0err;}
+
   Double_t GetKappa() const {return fKappa;}
   Double_t GetRadius() const {return fRadius;}
   Double_t GetCenterX() const {return fCenterX;}
@@ -141,6 +143,7 @@ class AliHLTTPCTrack : public AliTPCtrack {
   void SetPterr(Double_t f) {fPterr = f;}
   void SetPsierr(Double_t f) {fPsierr = f;}
   void SetZ0err(Double_t f) {fZ0err = f;}
+  void SetY0err(Double_t f) {fY0err = f;}  
   void SetTglerr(Double_t f) {fTanlerr = f;}
   void SetKappa(Double_t f) {fKappa = f;}
   void SetNHits(Int_t f) {fNHits = f;}
@@ -187,6 +190,7 @@ class AliHLTTPCTrack : public AliTPCtrack {
   Double_t fPterr;   //error in pt
   Double_t fPsierr;  //error in psi
   Double_t fZ0err;   //error in first point
+  Double_t fY0err;   //error in first point
   Double_t fTanlerr; //error in tanl
 
   Double_t fPhi0; //azimuthal angle of the first point