X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=TOF%2FAliTOFtrackerMI.cxx;h=8a9359000ab8c85e461c36f24a84bd7a211927f3;hb=cf9f6481fd65a9c59660c943ed4f44aaad41d3b5;hp=a2b56493ebc977e56e06701bd87856bb1ef45fdb;hpb=de60fa8a10952e57c00273dc54ea095d26062666;p=u%2Fmrichter%2FAliRoot.git diff --git a/TOF/AliTOFtrackerMI.cxx b/TOF/AliTOFtrackerMI.cxx index a2b56493ebc..8a9359000ab 100644 --- a/TOF/AliTOFtrackerMI.cxx +++ b/TOF/AliTOFtrackerMI.cxx @@ -70,6 +70,8 @@ AliTOFtrackerMI::AliTOFtrackerMI(): { //AliTOFtrackerMI main Ctor + for (Int_t ii=0; iiDelete(); + } //_____________________________________________________________________________ void AliTOFtrackerMI::GetPidSettings(AliESDpid *esdPID) { @@ -271,12 +277,12 @@ void AliTOFtrackerMI::MatchTracks( Bool_t /*mLastStep*/) const { void AliTOFtrackerMI::MatchTracksMI(Bool_t mLastStep){ //Match ESD tracks to clusters in TOF - const Float_t kTofOffset = 26; // time offset + const Float_t kTofOffset = 0; // time offset const Float_t kMinQuality = -6.; // minimal quality const Float_t kMaxQualityD = 1.; // max delta quality if cluster used const Float_t kForbiddenR = 0.1; // minimal PID according TPC - static const Double_t kMasses[]={ + static const Double_t kMasses[AliPID::kSPECIES+1]={ 0.000511, 0.105658, 0.139570, 0.493677, 0.938272, 1.875613 }; @@ -297,7 +303,7 @@ void AliTOFtrackerMI::MatchTracksMI(Bool_t mLastStep){ Float_t mintimedist[kNclusterMax]; Float_t likelihood[kNclusterMax]; Float_t length[kNclusterMax]; - Double_t tpcpid[5]; + Double_t tpcpid[AliPID::kSPECIES+1]; // overrun_static - coverity warning dist3D[0][0]=1; for (Int_t i=0; iGetTPCpid(tpcpid); Double_t sumpid=0; - for (Int_t ipid=0;ipid<5;ipid++){ + for (Int_t ipid=0;ipid0) tpcpid[ipid]/=sumpid; else{ tpcpid[ipid]=0.2; @@ -433,7 +439,8 @@ void AliTOFtrackerMI::MatchTracksMI(Bool_t mLastStep){ Double_t tof2=AliTOFGeometry::TdcBinWidth()*cluster->GetTDC()+kTofOffset; // in ps // Float_t tgamma = TMath::Sqrt(cluster->GetR()*cluster->GetR()+cluster->GetZ()*cluster->GetZ())/0.03; //time for "primary" gamma //if (trackTOFin->GetPt()<0.7 && TMath::Abs(tgamma-tof2)<350) continue; // gamma conversion candidate - TEMPORARY - for(Int_t j=0;j<=5;j++){ + for(Int_t j=0;jkForbiddenR){ @@ -792,8 +799,8 @@ void AliTOFtrackerMI::GetLikelihood(Float_t dy, Float_t dz, const Double_t *cov, Float_t normwidth, normd, p0,p1; - Float_t sigmay = TMath::Max(TMath::Sqrt(cov[0]+kMeanSigmaY*kMeanSigmaY),kMaxSigmaY); - Float_t sigmaz = TMath::Max(TMath::Sqrt(cov[2]+kMeanSigmaZ*kMeanSigmaZ),kMaxSigmaZ); + Float_t sigmay = TMath::Max(TMath::Sqrt(TMath::Abs(cov[0])+kMeanSigmaY*kMeanSigmaY),kMaxSigmaY); + Float_t sigmaz = TMath::Max(TMath::Sqrt(TMath::Abs(cov[2])+kMeanSigmaZ*kMeanSigmaZ),kMaxSigmaZ); py=0; pz=0;