Usage of the error and shape parameterization
authormarian <marian@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sat, 8 Nov 2008 10:32:14 +0000 (10:32 +0000)
committermarian <marian@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sat, 8 Nov 2008 10:32:14 +0000 (10:32 +0000)
(Marian)

TPC/AliTPCTracklet.cxx
TPC/AliTPCcalibCalib.cxx
TPC/AliTPCcalibCalib.h

index e8b77d5..9eddd5a 100755 (executable)
@@ -326,6 +326,9 @@ void AliTPCTracklet::FitKalman(const AliTPCseed *seed,Int_t sector) {
     if (!c) continue;
     Double_t r[3]={c->GetX(),c->GetY(),c->GetZ()};
     Double_t cov[3]={0.01,0.,0.01}; //TODO: correct error parametrisation
+    AliTPCseed::GetError(c, track,cov[0],cov[2]);
+    cov[0]*=cov[0];
+    cov[2]*=cov[2];
     if (!track->PropagateTo(r[0])) {
       isOK=kFALSE;
       break;
@@ -342,7 +345,10 @@ void AliTPCTracklet::FitKalman(const AliTPCseed *seed,Int_t sector) {
     AliTPCclusterMI *c=track->GetClusterPointer(i);
     if (!c) continue;
     Double_t r[3]={c->GetX(),c->GetY(),c->GetZ()};
-    Double_t cov[3]={0.01,0.,0.01}; //TODO: correct error parametrisation
+    Double_t cov[3]={0.01,0.,0.01}; 
+    AliTPCseed::GetError(c, track,cov[0],cov[2]);
+    cov[0]*=cov[0];
+    cov[2]*=cov[2];
     if (!track->PropagateTo(r[0])) {
       isOK=kFALSE;
       break;
@@ -359,7 +365,10 @@ void AliTPCTracklet::FitKalman(const AliTPCseed *seed,Int_t sector) {
     AliTPCclusterMI *c=track->GetClusterPointer(i);
     if (!c) continue;
     Double_t r[3]={c->GetX(),c->GetY(),c->GetZ()};
-    Double_t cov[3]={0.01,0.,0.01}; //TODO: correct error parametrisation
+    Double_t cov[3]={0.01,0.,0.01}; 
+    AliTPCseed::GetError(c, track,cov[0],cov[2]);
+    cov[0]*=cov[0];
+    cov[2]*=cov[2];
     if (!track->PropagateTo(r[0])) {
       isOK=kFALSE;
       break;
index a0b89b9..25034e2 100644 (file)
@@ -118,7 +118,7 @@ void     AliTPCcalibCalib::Process(AliESDEvent *event){
   
   if (GetDebugLevel()>20) printf("Hallo world: Im here\n");
   Int_t ntracks=event->GetNumberOfTracks();   
-  AliTPCcalibDB::Instance()->SetExBField(fMagF);
+  //AliTPCcalibDB::Instance()->SetExBField(fMagF);
 
   //
   //
@@ -227,7 +227,11 @@ Bool_t  AliTPCcalibCalib::RefitTrack(AliESDtrack * track, AliTPCseed *seed){
     if (TMath::Abs(dalpha)>0.01)
       trackOut.Rotate(TMath::DegToRad()*(sector%18*20.+10.));
     Double_t r[3]={cl->GetX(),cl->GetY(),cl->GetZ()};
-    Double_t cov[3]={0.01,0.,0.01}; //TODO: correct error parametrisation
+
+    Double_t cov[3]={0.01,0.,0.01}; //TODO: correct error parametrisation    
+    AliTPCseed::GetError(cl, &trackOut,cov[0],cov[2]);
+    cov[0]*=cov[0];
+    cov[2]*=cov[2];
     trackOut.GetXYZ(xyz);
     Double_t bz = AliTracker::GetBz(xyz);
     if (trackOut.PropagateTo(r[0],bz)) nclOut++;
@@ -248,7 +252,10 @@ Bool_t  AliTPCcalibCalib::RefitTrack(AliESDtrack * track, AliTPCseed *seed){
       trackIn.Rotate(TMath::DegToRad()*(sector%18*20.+10.));
     Double_t r[3]={cl->GetX(),cl->GetY(),cl->GetZ()};
     Double_t cov[3]={0.01,0.,0.01}; //TODO: correct error parametrisation
-    trackOut.GetXYZ(xyz);
+    AliTPCseed::GetError(cl, &trackIn,cov[0],cov[2]);
+    cov[0]*=cov[0];
+    cov[2]*=cov[2];
+    trackIn.GetXYZ(xyz);
     Double_t bz = AliTracker::GetBz(xyz);
 
     if (trackIn.PropagateTo(r[0],bz)) nclIn++;
@@ -314,3 +321,4 @@ Bool_t AliTPCcalibCalib::RejectCluster(AliTPCclusterMI* cl, AliExternalTrackPara
 }
 
 
+
index bb9b211..012e774 100644 (file)
@@ -31,7 +31,6 @@ public:
   Bool_t  RejectCluster(AliTPCclusterMI* cl, AliExternalTrackParam * param);
   void     Process(AliESDtrack *track, Int_t runNo=-1){AliTPCcalibBase::Process(track,runNo);};
   void     Process(AliTPCseed *track){return AliTPCcalibBase::Process(track);}
-
 protected: 
 private:
   ClassDef(AliTPCcalibCalib,1)