Skip to content
Snippets Groups Projects
Commit 52f8eeb3 authored by Vardan Tadevosyan's avatar Vardan Tadevosyan Committed by Stephen A. Wood
Browse files

Make counting of calorimeter tracks consistent with Engine.

Increment fNtracks++ each time a cluster matches a spectrometer
track. Only best matches were counted before. Move fNtracks++ from
THcShower::CoarceProcess to THcShower:MatchCluster.
parent 4060f7f9
No related branches found
No related tags found
No related merge requests found
...@@ -740,7 +740,7 @@ Int_t THcShower::CoarseProcess( TClonesArray& tracks) ...@@ -740,7 +740,7 @@ Int_t THcShower::CoarseProcess( TClonesArray& tracks)
Int_t mclust = MatchCluster(theTrack, ClusterList, Xtr, Ytr); Int_t mclust = MatchCluster(theTrack, ClusterList, Xtr, Ytr);
if (mclust >= 0) fNtracks++; // number of shower tracks // if (mclust >= 0) fNtracks++; // number of shower tracks (This is not consistent with engine!)
// Do this for the 1-st track only for now. // Do this for the 1-st track only for now.
// //
...@@ -907,7 +907,7 @@ Int_t THcShower::MatchCluster(THaTrack* Track, ...@@ -907,7 +907,7 @@ Int_t THcShower::MatchCluster(THaTrack* Track,
Double_t dx = TMath::Abs( (*cluster).clX() - XTrFront ); Double_t dx = TMath::Abs( (*cluster).clX() - XTrFront );
if (dx <= (0.5*BlockThick[0] + fSlop)) { if (dx <= (0.5*BlockThick[0] + fSlop)) {
fNtracks++; // number of shower tracks (Consistent with engine)
if (dx <= deltaX) { if (dx <= deltaX) {
mclust = i; mclust = i;
deltaX = dx; deltaX = dx;
......
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