Skip to content
Snippets Groups Projects
Commit 169209c4 authored by fjlinz's avatar fjlinz Committed by Administrator
Browse files

Remove broken hits from eloss calculation in TRD

MR 1634
(cherry picked from commit 01616065)
parent 22fd47e9
No related branches found
No related tags found
No related merge requests found
...@@ -92,7 +92,6 @@ void CbmTrdTracksConverter::ProcessData(CbmEvent* event) ...@@ -92,7 +92,6 @@ void CbmTrdTracksConverter::ProcessData(CbmEvent* event)
trd_track->GetParamLast()->Momentum(mom_last); trd_track->GetParamLast()->Momentum(mom_last);
track.SetMomentum3(mom); track.SetMomentum3(mom);
track.SetField(int(trd_track->GetNofHits()), i_n_hits_);
track.SetField(float(trd_track->GetPidLikeEL()), i_pid_like_); track.SetField(float(trd_track->GetPidLikeEL()), i_pid_like_);
track.SetField(float(trd_track->GetPidLikePI()), i_pid_like_ + 1); track.SetField(float(trd_track->GetPidLikePI()), i_pid_like_ + 1);
...@@ -109,14 +108,20 @@ void CbmTrdTracksConverter::ProcessData(CbmEvent* event) ...@@ -109,14 +108,20 @@ void CbmTrdTracksConverter::ProcessData(CbmEvent* event)
track.SetField(0.f, i_e_loss_i_ + i); track.SetField(0.f, i_e_loss_i_ + i);
} }
int trd_hits = trd_track->GetNofHits();
for (Int_t ihit = 0; ihit < trd_track->GetNofHits(); ihit++) { for (Int_t ihit = 0; ihit < trd_track->GetNofHits(); ihit++) {
Int_t idx = trd_track->GetHitIndex(ihit); Int_t idx = trd_track->GetHitIndex(ihit);
auto* hit = (CbmTrdHit*) cbm_trd_hits_->At(idx); auto* hit = (CbmTrdHit*) cbm_trd_hits_->At(idx);
if (hit) { if (hit) {
// std::cout << hit->GetELoss()*1e6 << " " << hit->GetPlaneId() << std::endl; if (hit->GetELoss() > 0) {
track.SetField(float(hit->GetELoss() * 1e6), i_e_loss_i_ + hit->GetPlaneId()); track.SetField(float(hit->GetELoss() * 1e6), i_e_loss_i_ + hit->GetPlaneId());
} else {
trd_hits--;
}
} }
} }
track.SetField(trd_hits, i_n_hits_);
if (rec_tracks_map.empty()) { continue; } if (rec_tracks_map.empty()) { continue; }
const Int_t stsTrackIndex = global_track->GetStsTrackIndex(); const Int_t stsTrackIndex = global_track->GetStsTrackIndex();
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment