Bug fix in PID object definition

diff --git a/TRD/AliTRDcalibDB.cxx b/TRD/AliTRDcalibDB.cxx
index e7a45c1136a46ffac9aab6c5eac2917b9f823724..5ace68fa96f5fe05d95d668493a2beaa02b880b0 100644 (file)
--- a/TRD/AliTRDcalibDB.cxx
+++ b/TRD/AliTRDcalibDB.cxx
@@ -846,12 +846,12 @@ const AliTRDCalPID *AliTRDcalibDB::GetPIDObject(const Int_t method)
 
   switch(method) {
     case 0: return dynamic_cast<const AliTRDCalPID *> 
-                              (GetCachedCDBObject(kIDPIDNN));
-    case 1: return dynamic_cast<const AliTRDCalPID *>
                               (GetCachedCDBObject(kIDPIDLQ));
+    case 1: return dynamic_cast<const AliTRDCalPID *>
+                              (GetCachedCDBObject(kIDPIDNN));
   }
 
-  return 0;
+  return 0x0;
 
 }
 

diff --git a/TRD/AliTRDrecoParam.cxx b/TRD/AliTRDrecoParam.cxx
index 92802dc54b75a3f32d5a2821e6a676455b79bd2a..87a9410dda20667d9ec0d625e6b9f98ff85ffe50 100644 (file)
--- a/TRD/AliTRDrecoParam.cxx
+++ b/TRD/AliTRDrecoParam.cxx
@@ -34,7 +34,7 @@ ClassImp(AliTRDrecoParam)
 AliTRDrecoParam::AliTRDrecoParam()
   :AliDetectorRecoParam()
   ,fkClusterSharing(0)
-  ,fkPIDMethod(1) // LQ PID
+  ,fkPIDMethod(0) // LQ PID
   ,fkMaxTheta(1.0)
   ,fkMaxPhi(2.0)
   ,fkRoad0y(6.0)

diff --git a/TRD/AliTRDseedV1.cxx b/TRD/AliTRDseedV1.cxx
index bb31e5be92b5167eaca22a30cb4e815014dc3be4..bae8ec3ff6722c3f2ebd4142915365a8753af2ec 100644 (file)
--- a/TRD/AliTRDseedV1.cxx
+++ b/TRD/AliTRDseedV1.cxx
@@ -221,13 +221,15 @@ void AliTRDseedV1::CookdEdx(Int_t nslices)
        } // End of loop over clusters
 
        // calculate mean charge per slice
-       for(int is=0; is<nslices; is++) if(nclusters[is]) fdEdx[is] /= nclusters[is];
+       for(int is=0; is<nslices; is++){ 
+    if(nclusters[is]) fdEdx[is] /= nclusters[is];
+  }
 }
 
 //____________________________________________________________________
 Float_t AliTRDseedV1::GetdQdl(Int_t ic) const
 {
-       return fClusters[ic] ? TMath::Abs(fClusters[ic]->GetQ()) /fdX / TMath::Sqrt(1. + fYfit[1]*fYfit[1] + fZfit[1]*fZfit[1]) : 0.;
+       return fClusters[ic] ? TMath::Abs(fClusters[ic]->GetQ()) /fdX / TMath::Sqrt(1. + fYfit[1]*fYfit[1] + fZref[1]*fZref[1]) : 0.;
 }
 
 //____________________________________________________________________
@@ -256,7 +258,8 @@ Double_t* AliTRDseedV1::GetProbability()
     AliError("No access to AliTRDCalPID object");
     return 0x0;
   }
-       
+       //AliInfo(Form("Method[%d] : %s", AliTRDReconstructor::RecoParam()->GetPIDMethod(), pd->IsA()->GetName()));
+
        // calculate tracklet length TO DO
   Float_t length = (AliTRDgeometry::AmThick() + AliTRDgeometry::DrThick());
   /// TMath::Sqrt((1.0 - fSnp[iPlane]*fSnp[iPlane]) / (1.0 + fTgl[iPlane]*fTgl[iPlane]));