Bug fix in PID object definition
authorcblume <cblume@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 9 Apr 2008 12:11:27 +0000 (12:11 +0000)
committercblume <cblume@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 9 Apr 2008 12:11:27 +0000 (12:11 +0000)
TRD/AliTRDcalibDB.cxx
TRD/AliTRDrecoParam.cxx
TRD/AliTRDseedV1.cxx

index e7a45c1..5ace68f 100644 (file)
@@ -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;
 
 }
 
index 92802dc..87a9410 100644 (file)
@@ -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)
index bb31e5b..bae8ec3 100644 (file)
@@ -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]));