From 889fd9eb55c60fc308c2d457efbf75d9232cc4a6 Mon Sep 17 00:00:00 2001
From: Luisa Faber <l_fabe01@uni-muenster.de>
Date: Fri, 5 Apr 2024 14:48:09 +0200
Subject: [PATCH] Fix likelihood calculation for tracks with hits with
incomplete clusters
---
reco/detectors/trd/pid/CbmTrdSetTracksPidLike.cxx | 5 ++---
1 file changed, 2 insertions(+), 3 deletions(-)
diff --git a/reco/detectors/trd/pid/CbmTrdSetTracksPidLike.cxx b/reco/detectors/trd/pid/CbmTrdSetTracksPidLike.cxx
index 6100edd615..647103d7cb 100644
--- a/reco/detectors/trd/pid/CbmTrdSetTracksPidLike.cxx
+++ b/reco/detectors/trd/pid/CbmTrdSetTracksPidLike.cxx
@@ -13,12 +13,10 @@
#include "CbmTrdGas.h"
#include "CbmTrdHit.h"
#include "CbmTrdTrack.h"
-
#include "FairParamList.h"
#include "FairRootManager.h"
#include "FairRunAna.h"
#include "FairRuntimeDb.h"
-
#include "TClonesArray.h"
#include "TH1.h"
#include "TH2.h"
@@ -27,6 +25,7 @@
#include "TObjArray.h"
#include "TROOT.h"
#include "TString.h"
+
#include <TFile.h>
#include <iostream>
@@ -344,7 +343,7 @@ void CbmTrdSetTracksPidLike::Exec(Option_t*)
for (Int_t iTRD = 0; iTRD < pTrack->GetNofHits(); iTRD++) {
Int_t index = pTrack->GetHitIndex(iTRD);
CbmTrdHit* trdHit = (CbmTrdHit*) fTrdHitArray->At(index);
-
+ if (trdHit->GetELoss() < 0.) continue;
dEdx = trdHit->GetELoss() * 1.e+6; //GeV->keV
dEsum += trdHit->GetELoss() * 1.e+6; //GeV->keV
--
GitLab